The genomic determinants of adaptive evolution in a fungal pathogen

Unravelling the strength, frequency, and distribution of selective variants along the genome as well as the underlying factors shaping this distribution are fundamental goals of evolutionary biology. Antagonistic host-pathogen coevolution is thought to be a major driver of genome evolution between interacting species. While rapid evolution of pathogens has been documented in several model organisms, the genetic mechanisms of their adaptation are still poorly understood and debated, particularly the role of sexual reproduction. Here, we apply a population genomic approach to infer genome-wide patterns of selection among 13 isolates of Zymoseptoria tritici, a fungal pathogen characterized by extremely high genetic diversity, gene density, and recombination rates. We report that the genome of Z. tritici undergoes a high rate of adaptive substitutions, with 44% of nonsynonymous substitutions being adaptive on average. This fraction reaches 68% in so-called effector genes encoding determinants of pathogenicity, and the distribution of fitness effects differs in this class of genes as they undergo adaptive mutations with stronger positive fitness effects, but also more slightly deleterious mutations. Besides the globally high rate of adaptive substitutions, we report a negative relationship between pN/pS and the fine-scale recombination rate and a strong positive correlation between the rate of adaptive nonsynonymous substitutions (ωa) and recombination rate. This result suggests a pervasive role of both background selection and Hill-Robertson interference even in a species with an exceptionally high recombination rate (60 cM/Mb on average). While transposable elements (TEs) have been suggested to contribute to adaptation by creating compartments of fast-evolving genomic regions, we do not find a significant effect of TEs on the rate of adaptive mutations. Overall our study suggests that sexual recombination is a significant driver of genome evolution, even in rapidly evolving organisms subject to recurrent mutations with large positive effects.

Highly flexible infection programs in a specialized wheat pathogen

Many filamentous plant pathogens exhibit high levels of genomic variability, yet the impact of this variation on host-pathogen interactions is largely unknown. We have addressed host specialization in the wheat pathogen Zymoseptoria tritici. Our study builds on comparative analyses of infection and gene expression phenotypes of three isolates and reveals the extent to which genomic variation translates into phenotypic variation. The isolates exhibit genetic and genomic variation but are similarly virulent. By combining confocal microscopy, disease monitoring, staining of ROS, and comparative transcriptome analyses, we conducted a detailed comparison of the infection processes of these isolates in a susceptible wheat cultivar. We characterized four core infection stages: establishment, biotrophic growth, lifestyle transition, and necrotrophic growth and asexual reproduction that are shared by the three isolates. However, we demonstrate differentiated temporal and spatial infection development and significant differences in the expression profiles of the three isolates during the infection stages. More than 20% of the genes were differentially expressed and these genes were located significantly closer to transposable elements, suggesting an impact of epigenetic regulation. Further, differentially expressed genes were enriched in effector candidates suggesting that isolate-specific strategies for manipulating host defenses are present in Z. tritici. We demonstrate that individuals of a host-specialized pathogen have highly differentiated infection programs characterized by flexible infection development and functional redundancy. This illustrates how high genetic diversity in pathogen populations results in highly differentiated infection phenotypes, which fact needs to be acknowledged to understand host-pathogen interactions and pathogen evolution.

Histone modifications rather than the novel regional centromeres of Zymoseptoria tritici distinguish core and accessory chromosomes

BACKGROUND

Supernumerary chromosomes have been found in many organisms. In fungi, these "accessory" or "dispensable" chromosomes are present at different frequencies in populations and are usually characterized by higher repetitive DNA content and lower gene density when compared to the core chromosomes. In the reference strain of the wheat pathogen, Zymoseptoria tritici, eight discrete accessory chromosomes have been found. So far, no functional role has been assigned to these chromosomes; however, they have existed as separate entities in the karyotypes of Zymoseptoria species over evolutionary time. In this study, we addressed what-if anything-distinguishes the chromatin of accessory chromosomes from core chromosomes. We used chromatin immunoprecipitation combined with high-throughput sequencing ("ChIP-seq") of DNA associated with the centromere-specific histone H3, CENP-A (CenH3), to identify centromeric DNA, and ChIP-seq with antibodies against dimethylated H3K4, trimethylated H3K9 and trimethylated H3K27 to determine the relative distribution and proportion of euchromatin, obligate and facultative heterochromatin, respectively.

RESULTS

Centromeres of the eight accessory chromosomes have the same sequence composition and structure as centromeres of the 13 core chromosomes and they are of similar length. Unlike those of most other fungi, Z. tritici centromeres are not composed entirely of repetitive DNA; some centromeres contain only unique DNA sequences, and bona fide expressed genes are located in regions enriched with CenH3. By fluorescence microscopy, we showed that centromeres of Z. tritici do not cluster into a single chromocenter during interphase. We found dramatically higher enrichment of H3K9me3 and H3K27me3 on the accessory chromosomes, consistent with the twofold higher proportion of repetitive DNA and poorly transcribed genes. In contrast, no single histone modification tested here correlated with the distribution of centromeric nucleosomes.

CONCLUSIONS

All centromeres are similar in length and composed of a mixture of unique and repeat DNA, and most contain actively transcribed genes. Centromeres, subtelomeric regions or telomere repeat length cannot account for the differences in transfer fidelity between core and accessory chromosomes, but accessory chromosomes are greatly enriched in nucleosomes with H3K27 trimethylation. Genes on accessory chromosomes appear to be silenced by trimethylation of H3K9 and H3K27.

Genomes and transcriptomes of partners in plant-fungal-interactions between canola (Brassica napus) and two Leptosphaeria species

Leptosphaeria maculans ‘brassicae’ is a damaging fungal pathogen of canola (Brassica napus), causing lesions on cotyledons and leaves, and cankers on the lower stem. A related species, L. biglobosa ‘canadensis’, colonises cotyledons but causes few stem cankers. We describe the complement of genes encoding carbohydrate-active enzymes (CAZys) and peptidases of these fungi, as well as of four related plant pathogens. We also report dual-organism RNA-seq transcriptomes of these two Leptosphaeria species and B. napus during disease. During the first seven days of infection L. biglobosa ‘canadensis’, a necrotroph, expressed more cell wall degrading genes than L. maculans ‘brassicae’, a hemi-biotroph. L. maculans ‘brassicae’ expressed many genes in the Carbohydrate Binding Module class of CAZy, particularly CBM50 genes, with potential roles in the evasion of basal innate immunity in the host plant. At this time, three avirulence genes were amongst the top 20 most highly upregulated L. maculans ‘brassicae’ genes in planta. The two fungi had a similar number of peptidase genes, and trypsin was transcribed at high levels by both fungi early in infection. L. biglobosa ‘canadensis’ infection activated the jasmonic acid and salicylic acid defence pathways in B. napus, consistent with defence against necrotrophs. L. maculans ‘brassicae’ triggered a high level of expression of isochorismate synthase 1, a reporter for salicylic acid signalling. L. biglobosa ‘canadensis’ infection triggered coordinated shutdown of photosynthesis genes, and a concomitant increase in transcription of cell wall remodelling genes of the host plant. Expression of particular classes of CAZy genes and the triggering of host defence and particular metabolic pathways are consistent with the necrotrophic lifestyle of L. biglobosa ‘canadensis’, and the hemibiotrophic life style of L. maculans ‘brassicae’.

[Transposable elements reshaping genomes and favouring the evolutionary and adaptive potential of fungal phytopathogens]

Phytopathogenic fungi are a major threat for global food security and show an extreme plasticity in pathogenicity behaviours. They often have a high adaptive potential allowing them to rapidly counteract the control methods used by men in agrosystems. In this paper, we evaluate the link between genome plasticity and adaptive potential using genomics and comparative genomics approaches. Our model is the evolutionary series Leptosphaeria maculans-Leptosphaeria biglobosa, encompassing five distinct entities, whose conspecificity or heterospecificity status is unclear, and which all are pathogens of cruciferous plants. They however differ by their host range and pathogenicity. Compared to other species of the species complex, the species best adapted to oilseed rape, L. maculans "brassicae", causing important losses in the crop, has a genome that was submitted to a recent and massive burst of transposition by a few families of transposable elements (TEs). Whether the genome invasion contributed to speciation is still unclear to-date but there is a coincidence between this burst of TEs and divergence between two species. This TE burst contributed to diversification of effector proteins and thus to generation of novel pathogenic specificities. In addition, the location of effector genes within genome regions enriched in TEs has direct consequences on adaptation to plant resistance and favours a multiplicity of mutation events allowing "breakdown" of resistance. These data are substantiated by other examples in the literature showing that fungi tend to have a "two-speed" genome, in which a plastic compartment enriched in TE host genes is involved in pathogenicity and adaptation to host.

The dispensable chromosome of Leptosphaeria maculans shelters an effector gene conferring avirulence towards Brassica rapa

Phytopathogenic fungi frequently contain dispensable chromosomes, some of which contribute to host range or pathogenicity. In Leptosphaeria maculans, the stem canker agent of oilseed rape (Brassica napus), the minichromosome was previously suggested to be dispensable, without evidence for any role in pathogenicity. Using genetic and genomic approaches, we investigated the inheritance and molecular determinant of an L. maculans-Brassica rapa incompatible interaction. Single gene control of the resistance was found, while all markers located on the L. maculans minichromosome, absent in the virulent parental isolate, co-segregated with the avirulent phenotype. Only one candidate avirulence gene was identified on the minichromosome, validated by complementation experiments and termed AvrLm11. The minichromosome was frequently lost following meiosis, but the frequency of isolates lacking it remained stable in field populations sampled at a 10-yr time interval, despite a yearly sexual stage in the L. maculans life cycle. This work led to the cloning of a new 'lost in the middle of nowhere' avirulence gene of L. maculans, interacting with a B. rapa resistance gene termed Rlm11 and introgressed into B. napus. It demonstrated the dispensability of the L. maculans minichromosome and suggested that its loss generates a fitness deficit.

Effector diversification within compartments of the Leptosphaeria maculans genome affected by Repeat-Induced Point mutations

Fungi are of primary ecological, biotechnological and economic importance. Many fundamental biological processes that are shared by animals and fungi are studied in fungi due to their experimental tractability. Many fungi are pathogens or mutualists and are model systems to analyse effector genes and their mechanisms of diversification. In this study, we report the genome sequence of the phytopathogenic ascomycete Leptosphaeria maculans and characterize its repertoire of protein effectors. The L. maculans genome has an unusual bipartite structure with alternating distinct guanine and cytosine-equilibrated and adenine and thymine (AT)-rich blocks of homogenous nucleotide composition. The AT-rich blocks comprise one-third of the genome and contain effector genes and families of transposable elements, both of which are affected by repeat-induced point mutation, a fungal-specific genome defence mechanism. This genomic environment for effectors promotes rapid sequence diversification and underpins the evolutionary potential of the fungus to adapt rapidly to novel host-derived constraints.

FONZIE: An optimized pipeline for minisatellite marker discovery and primer design from large sequence data sets

BACKGROUND

Micro-and minisatellites are among the most powerful genetic markers known to date. They have been used as tools for a large number of applications ranging from gene mapping to phylogenetic studies and isolate typing. However, identifying micro-and minisatellite markers on large sequence data sets is often a laborious process.

RESULTS

FONZIE was designed to successively 1) perform a search for markers via the external software Tandem Repeat Finder, 2) exclude user-defined specific genomic regions, 3) screen for the size and the percent matches of each relevant marker found by Tandem Repeat Finder, 4) evaluate marker specificity (i.e., occurrence of the marker as a single copy in the genome) using BLAST2.0, 5) design minisatellite primer pairs via the external software Primer3, and 6) check the specificity of each final PCR product by BLAST. A final file returns to users all the results required to amplify markers. A biological validation of the approach was performed using the whole genome sequence of the phytopathogenic fungus Leptosphaeria maculans, showing that more than 90% of the minisatellite primer pairs generated by the pipeline amplified a PCR product, 44.8% of which showed agarose-gel resolvable polymorphism between isolates. Segregation analyses confirmed that the polymorphic minisatellites corresponded to single-locus markers.

CONCLUSION

FONZIE is a stand-alone and user-friendly application developed to minimize tedious manual operations, reduce errors, and speed up the search for efficient minisatellite and microsatellite markers departing from whole-genome sequence data. This pipeline facilitates the integration of data and provides a set of specific primer sequences for PCR amplification of single-locus markers. FONZIE is freely downloadable at: http://www.versailles-grignon.inra.fr/bioger/equipes/leptosphaeria_maculans/outils_d_analyses/fonzie.