A systematic screen for co-option of transposable elements across the fungal kingdom.

How novel protein functions are acquired is a central question in molecular biology. Key paths to novelty include gene duplications, recombination or horizontal acquisition. Transposable elements (TEs) are increasingly recognized as a major source of novel domain-encoding sequences. However, the impact of TE coding sequences on the evolution of the proteome remains understudied. Here, we analyzed 1237 genomes spanning the phylogenetic breadth of the fungal kingdom. We scanned proteomes for evidence of co-occurrence of TE-derived domains along with other conventional protein functional domains. We detected more than 13,000 predicted proteins containing potentially TE-derived domain, of which 825 were identified in more than five genomes, indicating that many host-TE fusions may have persisted over long evolutionary time scales. We used the phylogenetic context to identify the origin and retention of individual TE-derived domains. The most common TE-derived domains are helicases derived from Academ, Kolobok or Helitron. We found putative TE co-options at a higher rate in genomes of the Saccharomycotina, providing an unexpected source of protein novelty in these generally TE depleted genomes. We investigated in detail a candidate host-TE fusion with a heterochromatic transcriptional silencing function that may play a role in TE and gene regulation in ascomycetes. The affected gene underwent multiple full or partial losses within the phylum. Overall, our work establishes a kingdom-wide view of putative host-TE fusions and facilitates systematic investigations of candidate fusion proteins.

Recent reactivation of a pathogenicity-associated transposable element is associated with major chromosomal rearrangements in a fungal wheat pathogen.

Transposable elements (TEs) are key drivers of genomic variation contributing to recent adaptation in most species. Yet, the evolutionary origins and insertion dynamics within species remain poorly understood. We recapitulate the spread of the pathogenicity-associated Styx element across five species that last diverged ∼11 000 years ago. We show that the element likely originated in the Zymoseptoria fungal pathogen genus and underwent multiple independent reactivation events. Using a global 900-genome panel of the wheat pathogen Zymoseptoria tritici, we assess Styx copy number variation and identify renewed transposition activity in Oceania and South America. We show that the element can mobilize to create additional Styx copies in a four-generation pedigree. Importantly, we find that new copies of the element are not affected by genomic defenses suggesting minimal control against the element. Styx copies are preferentially located in recombination breakpoints and likely triggered multiple types of large chromosomal rearrangements. Taken together, we establish the origin, diversification and reactivation of a highly active TE with likely major consequences for chromosomal integrity and the expression of disease.

A population-level invasion by transposable elements triggers genome expansion in a fungal pathogen.

Genome evolution is driven by the activity of transposable elements (TEs). The spread of TEs can have deleterious effects including the destabilization of genome integrity and expansions. However, the precise triggers of genome expansions remain poorly understood because genome size evolution is typically investigated only among deeply divergent lineages. Here, we use a large population genomics dataset of 284 individuals from populations across the globe of Zymoseptoria tritici, a major fungal wheat pathogen. We built a robust map of genome-wide TE insertions and deletions to track a total of 2456 polymorphic loci within the species. We show that purifying selection substantially depressed TE frequencies in most populations, but some rare TEs have recently risen in frequency and likely confer benefits. We found that specific TE families have undergone a substantial genome-wide expansion from the pathogen's center of origin to more recently founded populations. The most dramatic increase in TE insertions occurred between a pair of North American populations collected in the same field at an interval of 25 years. We find that both genome-wide counts of TE insertions and genome size have increased with colonization bottlenecks. Hence, the demographic history likely played a major role in shaping genome evolution within the species. We show that both the activation of specific TEs and relaxed purifying selection underpin this incipient expansion of the genome. Our study establishes a model to recapitulate TE-driven genome evolution over deeper evolutionary timescales.

Machine-learning predicts genomic determinants of meiosis-driven structural variation in a eukaryotic pathogen.

Species harbor extensive structural variation underpinning recent adaptive evolution. However, the causality between genomic features and the induction of new rearrangements is poorly established. Here, we analyze a global set of telomere-to-telomere genome assemblies of a fungal pathogen of wheat to establish a nucleotide-level map of structural variation. We show that the recent emergence of pesticide resistance has been disproportionally driven by rearrangements. We use machine learning to train a model on structural variation events based on 30 chromosomal sequence features. We show that base composition and gene density are the major determinants of structural variation. Retrotransposons explain most inversion, indel and duplication events. We apply our model to Arabidopsis thaliana and show that our approach extends to more complex genomes. Finally, we analyze complete genomes of haploid offspring in a four-generation pedigree. Meiotic crossover locations are enriched for new rearrangements consistent with crossovers being mutational hotspots. The model trained on species-wide structural variation accurately predicts the position of >74% of newly generated variants along the pedigree. The predictive power highlights causality between specific sequence features and the induction of chromosomal rearrangements. Our work demonstrates that training sequence-derived models can accurately identify regions of intrinsic DNA instability in eukaryotic genomes.

Rapid sequence evolution driven by transposable elements at a virulence locus in a fungal wheat pathogen.

BACKGROUND: Plant pathogens cause substantial crop losses in agriculture production and threaten food security. Plants evolved the ability to recognize virulence factors and pathogens have repeatedly escaped recognition due rapid evolutionary change at pathogen virulence loci (i.e. effector genes). The presence of transposable elements (TEs) in close physical proximity of effector genes can have important consequences for gene regulation and sequence evolution. Species-wide investigations of effector gene loci remain rare hindering our ability to predict pathogen evolvability. RESULTS: Here, we performed genome-wide association studies (GWAS) on a highly polymorphic mapping population of 120 isolates of Zymoseptoria tritici, the most damaging pathogen of wheat in Europe. We identified a major locus underlying significant variation in reproductive success of the pathogen and damage caused on the wheat cultivar Claro. The most strongly associated locus is intergenic and flanked by genes encoding a predicted effector and a serine-type endopeptidase. The center of the locus contained a highly dynamic region consisting of multiple families of TEs. Based on a large global collection of assembled genomes, we show that the virulence locus has undergone substantial recent sequence evolution. Large insertion and deletion events generated length variation between the flanking genes by a factor of seven (5-35 kb). The locus showed also strong signatures of genomic defenses against TEs (i.e. RIP) contributing to the rapid diversification of the locus. CONCLUSIONS: In conjunction, our work highlights the power of combining GWAS and population-scale genome analyses to investigate major effect loci in pathogens.

Emergence and diversification of a highly invasive chestnut pathogen lineage across southeastern Europe.

Invasive microbial species constitute a major threat to biodiversity, agricultural production and human health. Invasions are often dominated by one or a small number of genotypes, yet the underlying factors driving invasions are poorly understood. The chestnut blight fungus Cryphonectria parasitica first decimated the North American chestnut, and a more recent outbreak threatens European chestnut stands. To unravel the chestnut blight invasion of southeastern Europe, we sequenced 230 genomes of predominantly European strains. Genotypes outside of the invasion zone showed high levels of diversity with evidence for frequent and ongoing recombination. The invasive lineage emerged from the highly diverse European genotype pool rather than a secondary introduction from Asia or North America. The expansion across southeastern Europe was mostly clonal and is dominated by a single mating type, suggesting a fitness advantage of asexual reproduction. Our findings show how an intermediary, highly diverse bridgehead population gave rise to an invasive, largely clonally expanding pathogen.

The rise and fall of genes: origins and functions of plant pathogen pangenomes.

Plant pathogens can rapidly overcome resistance of their hosts by mutating key pathogenicity genes encoding for effectors. Pathogen adaptation is fuelled by extensive genetic variability in populations and different strains may not share the same set of genes. Recently, such an intra-specific variation in gene content became formalized as pangenomes distinguishing core genes (i.e. shared) and accessory genes (i.e. lineage or strain-specific). Across pathogens species, key effectors tend to be part of the rapidly evolving accessory genome. Here, we show how the construction and analysis of pathogen pangenomes provide deep insights into the dynamic host adaptation process. We also discuss how pangenomes should ideally be built and how geography, niche and lifestyle likely determine pangenome sizes.

A 19-isolate reference-quality global pangenome for the fungal wheat pathogen Zymoseptoria tritici.

BACKGROUND: The gene content of a species largely governs its ecological interactions and adaptive potential. A species is therefore defined by both core genes shared between all individuals and accessory genes segregating presence-absence variation. There is growing evidence that eukaryotes, similar to bacteria, show intra-specific variability in gene content. However, it remains largely unknown how functionally relevant such a pangenome structure is for eukaryotes and what mechanisms underlie the emergence of highly polymorphic genome structures. RESULTS: Here, we establish a reference-quality pangenome of a fungal pathogen of wheat based on 19 complete genomes from isolates sampled across six continents. Zymoseptoria tritici causes substantial worldwide losses to wheat production due to rapidly evolved tolerance to fungicides and evasion of host resistance. We performed transcriptome-assisted annotations of each genome to construct a global pangenome. Major chromosomal rearrangements are segregating within the species and underlie extensive gene presence-absence variation. Conserved orthogroups account for only ~ 60% of the species pangenome. Investigating gene functions, we find that the accessory genome is enriched for pathogenesis-related functions and encodes genes involved in metabolite production, host tissue degradation and manipulation of the immune system. De novo transposon annotation of the 19 complete genomes shows that the highly diverse chromosomal structure is tightly associated with transposable element content. Furthermore, transposable element expansions likely underlie recent genome expansions within the species. CONCLUSIONS: Taken together, our work establishes a highly complex eukaryotic pangenome providing an unprecedented toolbox to study how pangenome structure impacts crop-pathogen interactions.

Stress-Driven Transposable Element De-repression Dynamics and Virulence Evolution in a Fungal Pathogen.

Transposable elements (TEs) are drivers of genome evolution and affect the expression landscape of the host genome. Stress is a major factor inducing TE activity; however, the regulatory mechanisms underlying de-repression are poorly understood. Plant pathogens are excellent models to dissect the impact of stress on TEs. The process of plant infection induces stress for the pathogen, and virulence factors (i.e., effectors) located in TE-rich regions become expressed. To dissect TE de-repression dynamics and contributions to virulence, we analyzed the TE expression landscape of four strains of the major wheat pathogen Zymoseptoria tritici. We experimentally exposed strains to nutrient starvation and host infection stress. Contrary to expectations, we show that the two distinct conditions induce the expression of different sets of TEs. In particular, the most highly expressed TEs, including miniature inverted-repeat transposable element and long terminal repeat-Gypsy element, show highly distinct de-repression across stress conditions. Both the genomic context of TEs and the genetic background stress (i.e., different strains harboring the same TEs) were major predictors of de-repression under stress. Gene expression profiles under stress varied significantly depending on the proximity to the closest TEs and genomic defenses against TEs were largely ineffective to prevent de-repression. Next, we analyzed the locus encoding the Avr3D1 effector. We show that the insertion and subsequent silencing of TEs in close proximity likely contributed to reduced expression and virulence on a specific wheat cultivar. The complexity of TE responsiveness to stress across genetic backgrounds and genomic locations demonstrates substantial intraspecific genetic variation to control TEs with consequences for virulence.

Expression polymorphism at the ARPC4 locus links the actin cytoskeleton with quantitative disease resistance to Sclerotinia sclerotiorum in Arabidopsis thaliana.

Quantitative disease resistance (QDR) is a form of plant immunity widespread in nature, and the only one active against broad host range fungal pathogens. The genetic determinants of QDR are complex and largely unknown, and are thought to rely partly on genes controlling plant morphology and development. We used genome-wide association mapping in Arabidopsis thaliana to identify ARPC4 as associated with QDR against the necrotrophic fungal pathogen Sclerotinia sclerotiorum. Mutants impaired in ARPC4 showed enhanced susceptibility to S. sclerotiorum, defects in the structure of the actin filaments and in their responsiveness to S. sclerotiorum. Disruption of ARPC4 also alters callose deposition and the expression of defense-related genes upon S. sclerotiorum infection. Analysis of ARPC4 diversity in A. thaliana identified one haplotype (ARPC4R ) showing a c. 1 kbp insertion in ARPC4 regulatory region and associated with higher level of QDR. Accessions from the ARPC4R haplotype showed enhanced ARPC4 expression upon S. sclerotiorum challenge, indicating that polymorphisms in ARPC4 regulatory region are associated with enhanced QDR. This work identifies a novel actor of plant QDR against a fungal pathogen and provides a prime example of genetic mechanisms leading to the recruitment of cell morphology processes in plant immunity.

Parallel evolution of the POQR prolyl oligo peptidase gene conferring plant quantitative disease resistance.

Plant pathogens with a broad host range are able to infect plant lineages that diverged over 100 million years ago. They exert similar and recurring constraints on the evolution of unrelated plant populations. Plants generally respond with quantitative disease resistance (QDR), a form of immunity relying on complex genetic determinants. In most cases, the molecular determinants of QDR and how they evolve is unknown. Here we identify in Arabidopsis thaliana a gene mediating QDR against Sclerotinia sclerotiorum, agent of the white mold disease, and provide evidence of its convergent evolution in multiple plant species. Using genome wide association mapping in A. thaliana, we associated the gene encoding the POQR prolyl-oligopeptidase with QDR against S. sclerotiorum. Loss of this gene compromised QDR against S. sclerotiorum but not against a bacterial pathogen. Natural diversity analysis associated POQR sequence with QDR. Remarkably, the same amino acid changes occurred after independent duplications of POQR in ancestors of multiple plant species, including A. thaliana and tomato. Genome-scale expression analyses revealed that parallel divergence in gene expression upon S. sclerotiorum infection is a frequent pattern in genes, such as POQR, that duplicated both in A. thaliana and tomato. Our study identifies a previously uncharacterized gene mediating QDR against S. sclerotiorum. It shows that some QDR determinants are conserved in distantly related plants and have emerged through the repeated use of similar genetic polymorphisms at different evolutionary time scales.

Codon optimization underpins generalist parasitism in fungi.

The range of hosts that parasites can infect is a key determinant of the emergence and spread of disease. Yet, the impact of host range variation on the evolution of parasite genomes remains unknown. Here, we show that codon optimization underlies genome adaptation in broad host range parasites. We found that the longer proteins encoded by broad host range fungi likely increase natural selection on codon optimization in these species. Accordingly, codon optimization correlates with host range across the fungal kingdom. At the species level, biased patterns of synonymous substitutions underpin increased codon optimization in a generalist but not a specialist fungal pathogen. Virulence genes were consistently enriched in highly codon-optimized genes of generalist but not specialist species. We conclude that codon optimization is related to the capacity of parasites to colonize multiple hosts. Our results link genome evolution and translational regulation to the long-term persistence of generalist parasitism.

Emerging Trends in Molecular Interactions between Plants and the Broad Host Range Fungal Pathogens Botrytis cinerea and Sclerotinia sclerotiorum.

Fungal plant pathogens are major threats to food security worldwide. Sclerotinia sclerotiorum and Botrytis cinerea are closely related Ascomycete plant pathogens causing mold diseases on hundreds of plant species. There is no genetic source of complete plant resistance to these broad host range pathogens known to date. Instead, natural plant populations show a continuum of resistance levels controlled by multiple genes, a phenotype designated as quantitative disease resistance. Little is known about the molecular mechanisms controlling the interaction between plants and S. sclerotiorum and B. cinerea but significant advances were made on this topic in the last years. This minireview highlights a selection of nine themes that emerged in recent research reports on the molecular bases of plant-S. sclerotiorum and plant-B. cinerea interactions. On the fungal side, this includes progress on understanding the role of oxalic acid, on the study of fungal small secreted proteins. Next, we discuss the exchanges of small RNA between organisms and the control of cell death in plant and fungi during pathogenic interactions. Finally on the plant side, we highlight defense priming by mechanical signals, the characterization of plant Receptor-like proteins and the hormone abscisic acid in the response to B. cinerea and S. sclerotiorum, the role of plant general transcription machinery and plant small bioactive peptides. These represent nine trends we selected as remarkable in our understanding of fungal molecules causing disease and plant mechanisms associated with disease resistance to two devastating broad host range fungi.

Common protein sequence signatures associate with Sclerotinia borealis lifestyle and secretion in fungal pathogens of the Sclerotiniaceae.

Fungal plant pathogens produce secreted proteins adapted to function outside fungal cells to facilitate colonization of their hosts. In many cases such as for fungi from the Sclerotiniaceae family the repertoire and function of secreted proteins remains elusive. In the Sclerotiniaceae, whereas Sclerotinia sclerotiorum and Botrytis cinerea are cosmopolitan broad host-range plant pathogens, Sclerotinia borealis has a psychrophilic lifestyle with a low optimal growth temperature, a narrow host range and geographic distribution. To spread successfully, S. borealis must synthesize proteins adapted to function in its specific environment. The search for signatures of adaptation to S. borealis lifestyle may therefore help revealing proteins critical for colonization of the environment by Sclerotiniaceae fungi. Here, we analyzed amino acids usage and intrinsic protein disorder in alignments of groups of orthologous proteins from the three Sclerotiniaceae species. We found that enrichment in Thr, depletion in Glu and Lys, and low disorder frequency in hot loops are significantly associated with S. borealis proteins. We designed an index to report bias in these properties and found that high index proteins were enriched among secreted proteins in the three Sclerotiniaceae fungi. High index proteins were also enriched in function associated with plant colonization in S. borealis, and in in planta-induced genes in S. sclerotiorum. We highlight a novel putative antifreeze protein and a novel putative lytic polysaccharide monooxygenase identified through our pipeline as candidate proteins involved in colonization of the environment. Our findings suggest that similar protein signatures associate with S. borealis lifestyle and with secretion in the Sclerotiniaceae. These signatures may be useful for identifying proteins of interest as targets for the management of plant diseases.