Population-level transposable element expression dynamics influence trait evolution in a fungal crop pathogen.

The rapid adaptive evolution of microbes is driven by strong selection pressure acting on genetic variation. How adaptive genetic variation is generated within species and how such variation influences phenotypic trait expression is often not well understood though. We focused on the recent activity of transposable elements (TEs) using deep population genomics and transcriptomics analyses of a fungal plant pathogen with a highly active content of TEs in the genome. Zymoseptoria tritici causes one of the most damaging diseases on wheat, with recent adaptation to the host and environment being facilitated by TE-associated mutations. We obtained genomic and RNA-sequencing data from 146 isolates collected from a single wheat field. We established a genome-wide map of TE insertion polymorphisms in the population by analyzing recent TE insertions among individuals. We quantified the locus-specific transcription of individual TE copies and found considerable population variation at individual TE loci in the population. About 20% of all TE copies show transcription in the genome suggesting that genomic defenses such as repressive epigenetic marks and repeat-induced polymorphisms are at least partially ineffective at preventing the proliferation of TEs in the genome. A quarter of recent TE insertions are associated with expression variation of neighboring genes providing broad potential to influence trait expression. We indeed found that TE insertions are likely responsible for variation in virulence on the host and potentially diverse components of secondary metabolite production. Our large-scale transcriptomics study emphasizes how TE-derived polymorphisms segregate even in individual microbial populations and can broadly underpin trait variation in pathogens.IMPORTANCEPathogens can rapidly adapt to new hosts, antimicrobials, or changes in the environment. Adaptation arises often from mutations in the genome; however, how such variation is generated remains poorly understood. We investigated the most dynamic regions of the genome of Zymoseptoria tritici, a major fungal pathogen of wheat. We focused on the transcription of transposable elements. A large proportion of the transposable elements not only show signatures of potential activity but are also variable within a single population of the pathogen. We find that this variation in activity is likely influencing many important traits of the pathogen. Hence, our work provides insights into how a microbial species can adapt over the shortest time periods based on the activity of transposable elements.

Genome-wide expression QTL mapping reveals the highly dynamic regulatory landscape of a major wheat pathogen.

BACKGROUND: In agricultural ecosystems, outbreaks of diseases are frequent and pose a significant threat to food security. A successful pathogen undergoes a complex and well-timed sequence of regulatory changes to avoid detection by the host immune system; hence, well-tuned gene regulation is essential for survival. However, the extent to which the regulatory polymorphisms in a pathogen population provide an adaptive advantage is poorly understood. RESULTS: We used Zymoseptoria tritici, one of the most important pathogens of wheat, to generate a genome-wide map of regulatory polymorphism governing gene expression. We investigated genome-wide transcription levels of 146 strains grown under nutrient starvation and performed expression quantitative trait loci (eQTL) mapping. We identified cis-eQTLs for 65.3% of all genes and the majority of all eQTL loci are within 2kb upstream and downstream of the transcription start site (TSS). We also show that polymorphism in different gene elements contributes disproportionally to gene expression variation. Investigating regulatory polymorphism in gene categories, we found an enrichment of regulatory variants for genes predicted to be important for fungal pathogenesis but with comparatively low effect size, suggesting a separate layer of gene regulation involving epigenetics. We also show that previously reported trait-associated SNPs in pathogen populations are frequently cis-regulatory variants of neighboring genes with implications for the trait architecture. CONCLUSIONS: Overall, our study provides extensive evidence that single populations segregate large-scale regulatory variation and are likely to fuel rapid adaptation to resistant hosts and environmental change.

Genome-wide association mapping reveals genes underlying population-level metabolome diversity in a fungal crop pathogen.

BACKGROUND: Fungi produce a wide range of specialized metabolites (SMs) involved in biotic interactions. Pathways for the production of SMs are often encoded in clusters of tightly arranged genes identified as biosynthetic gene clusters. Such gene clusters can undergo horizontal gene transfers between species and rapid evolutionary change within species. The acquisition, rearrangement, and deletion of gene clusters can generate significant metabolome diversity. However, the genetic basis underlying variation in SM production remains poorly understood. RESULTS: Here, we analyzed the metabolite production of a large population of the fungal pathogen of wheat, Zymoseptoria tritici. The pathogen causes major yield losses and shows variation in gene clusters. We performed untargeted ultra-high performance liquid chromatography-high resolution mass spectrometry to profile the metabolite diversity among 102 isolates of the same species. We found substantial variation in the abundance of the detected metabolites among isolates. Integrating whole-genome sequencing data, we performed metabolite genome-wide association mapping to identify loci underlying variation in metabolite production (i.e., metabolite-GWAS). We found that significantly associated SNPs reside mostly in coding and gene regulatory regions. Associated genes encode mainly transport and catalytic activities. The metabolite-GWAS identified also a polymorphism in the 3'UTR region of a virulence gene related to metabolite production and showing expression variation. CONCLUSIONS: Taken together, our study provides a significant resource to unravel polymorphism underlying metabolome diversity within a species. Integrating metabolome screens should be feasible for a range of different plant pathogens and help prioritize molecular studies.

Histone H3K27 Methylation Perturbs Transcriptional Robustness and Underpins Dispensability of Highly Conserved Genes in Fungi.

Epigenetic modifications are key regulators of gene expression and underpin genome integrity. Yet, how epigenetic changes affect the evolution and transcriptional robustness of genes remains largely unknown. Here, we show how the repressive histone mark H3K27me3 underpins the trajectory of highly conserved genes in fungi. We first performed transcriptomic profiling on closely related species of the plant pathogen Fusarium graminearum species complex. We determined transcriptional responsiveness of genes across environmental conditions to determine expression robustness. To infer evolutionary conservation, we used a framework of 23 species across the Fusarium genus including three species covered with histone methylation data. Gene expression variation is negatively correlated with gene conservation confirming that highly conserved genes show higher expression robustness. In contrast, genes marked by H3K27me3 do not show such associations. Furthermore, highly conserved genes marked by H3K27me3 encode smaller proteins, exhibit weaker codon usage bias, higher levels of hydrophobicity, show lower intrinsically disordered regions, and are enriched for functions related to regulation and membrane transport. The evolutionary age of conserved genes with H3K27me3 histone marks falls typically within the origins of the Fusarium genus. We show that highly conserved genes marked by H3K27me3 are more likely to be dispensable for survival during host infection. Lastly, we show that conserved genes exposed to repressive H3K27me3 marks across distantly related Fusarium fungi are associated with transcriptional perturbation at the microevolutionary scale. In conclusion, we show how repressive histone marks are entangled in the evolutionary fate of highly conserved genes across evolutionary timescales.

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.

Optimization of parameters to improve transformation efficiency of elephant foot yam (Amorphophallus paeoniifolius (Dennst.) Nicolson.

Elephant foot yam (Amorphophallus paeoniifolius (Dennst.) Nicolson), is an important edible tropical tuber crop, belonging to the family Araceae. Corms produced by this plant is very big and they are rich in starch, protein, mineral, vitamins, and dietary fiber but has acridity problem. This crop is susceptible to virus and phytoplasma diseases which affects crop growth and corm yield. Even though this crop has high commercial value, the problems like susceptibility to viral diseases, acridity problems, and lack of genetic diversity made hindrance in their exploitation. These issues can be resolved only by improving the characters through genetic transformation. To achieve genetic transformation in this important crop, a study was conducted to optimize various parameters for efficient Agrobacterium-mediated genetic transformation using embryogenic calli with vectors having gus reporter gene. Calli were developed using petiole and leaves of in vitro plantlets of elephant foot yam cultivar Gajendra and experiments were conducted to evaluate the sensitivity of calli to different doses of antibiotics viz. geneticin, hygromycin, ticarcillin. It was observed that complete death and discoloration of the calli were obtained with 25 mgl-1 geneticin and 10 mgl-1 hygromycin. The lowest lethal concentration of ticarcillin against Agrobacterium growth was found to be 500 mgl-1 which did not affect calli growth. Optimized parameters for efficient transformation in elephant foot yam include 100 µM acetosyringone concentration with 2 days of co-cultivation at temperature 22 °C using LBA4404 strain. The putative transformants were characterized for the integration of the gus gene using PCR and nucleic acid spot hybridization. The optimized protocol is simple and reproducible and may be adapted for other cultivars also. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s13205-021-02824-6.

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.

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.