Utilizing Gene Tree Variation to Identify Candidate Effector Genes in Zymoseptoria tritici

Fungal plant pathogen "mutagenomics" reveals tagged and untagged mutations in Zymoseptoria tritici and identifies SSK2 as key morphogenesis and stress-responsive virulence factor

"Mutagenomics" is the combination of random mutagenesis, phenotypic screening, and whole-genome re-sequencing to uncover all tagged and untagged mutations linked with phenotypic changes in an organism. In this study, we performed a mutagenomics screen on the wheat pathogenic fungus Zymoseptoria tritici for altered morphogenetic switching and stress sensitivity phenotypes using Agrobacterium-mediated "random" T-DNA mutagenesis (ATMT). Biological screening identified four mutants which were strongly reduced in virulence on wheat. Whole genome re-sequencing defined the positions of the T-DNA insertion events and revealed several unlinked mutations potentially affecting gene functions. Remarkably, two independent reduced virulence mutant strains, with similarly altered stress sensitivities and aberrant hyphal growth phenotypes, were found to have a distinct loss of function mutations in the ZtSSK2 MAPKKK gene. One mutant strain had a direct T-DNA insertion affecting the predicted protein's N-terminus, while the other possessed an unlinked frameshift mutation towards the C-terminus. We used genetic complementation to restore both strains' wild-type (WT) function (virulence, morphogenesis, and stress response). We demonstrated that ZtSSK2 has a non-redundant function with ZtSTE11 in virulence through the biochemical activation of the stress-activated HOG1 MAPK pathway. Moreover, we present data suggesting that SSK2 has a unique role in activating this pathway in response to specific stresses. Finally, dual RNAseq-based transcriptome profiling of WT and SSK2 mutant strains revealed many HOG1-dependent transcriptional changes in the fungus during early infection and suggested that the host response does not discriminate between WT and mutant strains during this early phase. Together these data define new genes implicated in the virulence of the pathogen and emphasise the importance of a whole genome sequencing step in mutagenomic discovery pipelines.

The identification of a transposon affecting the asexual reproduction of the wheat pathogen Zymoseptoria tritici

Zymoseptoria tritici, the causal agent of Septoria tritici blotch, is a fungal wheat pathogen that causes significant global yield losses. Within Z. tritici populations, quantitative differences in virulence among different isolates are commonly observed; however, the genetic components that underpin these differences remain elusive. In this study, intraspecific comparative transcriptomic analysis was used to identify candidate genes that contribute to differences in virulence on the wheat cultivar WW2449. This led to the identification of a multicopy gene that was not expressed in the high-virulence isolate when compared to the medium- and low-virulence isolates. Further investigation suggested this gene resides in a 7.9-kb transposon. Subsequent long-read sequencing of the isolates used in the transcriptomic analysis confirmed that this gene did reside in an active Class II transposon, which is composed of four genes named REP9-1 to -4. Silencing and overexpression of REP9-1 in two distinct genetic backgrounds demonstrated that its expression alone reduces the number of pycnidia produced by Z. tritici during infection. The REP9-1 gene identified within a Class II transposon is the first discovery of a gene in a transposable element that influences the virulence of Z. tritici. This discovery adds further complexity to genetic loci that contribute to quantitative virulence in this important pathogen.

Recent loss of the Dim2 DNA methyltransferase decreases mutation rate in repeats and changes evolutionary trajectory in a fungal pathogen

DNA methylation is found throughout all domains of life, yet the extent and function of DNA methylation differ among eukaryotes. Strains of the plant pathogenic fungus Zymoseptoria tritici appeared to lack cytosine DNA methylation (5mC) because gene amplification followed by Repeat-Induced Point mutation (RIP) resulted in the inactivation of the dim2 DNA methyltransferase gene. 5mC is, however, present in closely related sister species. We demonstrate that inactivation of dim2 occurred recently as some Z. tritici isolates carry a functional dim2 gene. Moreover, we show that dim2 inactivation occurred by a different path than previously hypothesized. We mapped the genome-wide distribution of 5mC in strains with or without functional dim2 alleles. Presence of functional dim2 correlates with high levels of 5mC in transposable elements (TEs), suggesting a role in genome defense. We identified low levels of 5mC in strains carrying non-functional dim2 alleles, suggesting that 5mC is maintained over time, presumably by an active Dnmt5 DNA methyltransferase. Integration of a functional dim2 allele in strains with mutated dim2 restored normal 5mC levels, demonstrating de novo cytosine methylation activity of Dim2. To assess the importance of 5mC for genome evolution, we performed an evolution experiment, comparing genomes of strains with high levels of 5mC to genomes of strains lacking functional dim2. We found that presence of a functional dim2 allele alters nucleotide composition by promoting C to T transitions (C→T) specifically at CpA (CA) sites during mitosis, likely contributing to TE inactivation. Our results show that 5mC density at TEs is a polymorphic trait in Z. tritici populations that can impact genome evolution.

Interspecific Gene Exchange Introduces High Genetic Variability in Crop Pathogen

Genome analyses have revealed a profound role of hybridization and introgression in the evolution of many eukaryote lineages, including fungi. The impact of recurrent introgression on fungal evolution however remains elusive. Here, we analyzed signatures of introgression along the genome of the fungal wheat pathogen Zymoseptoria tritici. We applied a comparative population genomics approach, including genome data from five Zymoseptoria species, to characterize the distribution and composition of introgressed regions representing segments with an exceptional haplotype pattern. These regions are found throughout the genome, comprising 5% of the total genome and overlapping with > 1,000 predicted genes. We performed window-based phylogenetic analyses along the genome to distinguish regions which have a monophyletic or nonmonophyletic origin with Z. tritici sequences. A majority of nonmonophyletic windows overlap with the highly variable regions suggesting that these originate from introgression. We verified that incongruent gene genealogies do not result from incomplete lineage sorting by comparing the observed and expected length distribution of haplotype blocks resulting from incomplete lineage sorting. Although protein-coding genes are not enriched in these regions, we identify 18 that encode putative virulence determinants. Moreover, we find an enrichment of transposable elements in these regions implying that hybridization may contribute to the horizontal spread of transposable elements. We detected a similar pattern in the closely related species Zymoseptoria ardabiliae, suggesting that hybridization is widespread among these closely related grass pathogens. Overall, our results demonstrate a significant impact of recurrent hybridization on overall genome evolution of this important wheat pathogen.

Rapid Parallel Evolution of Azole Fungicide Resistance in Australian Populations of the Wheat Pathogen Zymoseptoria tritici

Zymoseptoria tritici is a globally distributed fungal pathogen which causes Septoria tritici blotch on wheat. Management of the disease is attempted through the deployment of resistant wheat cultivars and the application of fungicides. However, fungicide resistance is commonly observed in Z. tritici populations, and continuous monitoring is required to detect breakdowns in fungicide efficacy. We recently reported azole-resistant isolates in Australia; however, it remained unknown whether resistance was brought into the continent through gene flow or whether resistance emerged independently. To address this question, we screened 43 isolates across five Australian locations for azole sensitivity and performed whole-genome sequencing on 58 isolates from seven locations to determine the genetic basis of resistance. Population genomic analyses showed extremely strong differentiation between the Australian population recovered after azoles began to be used and both Australian populations recovered before azoles began to be used and populations on different continents. The apparent absence of recent gene flow between Australia and other continents suggests that azole fungicide resistance has evolved de novo and subsequently spread within Tasmania. Despite the isolates being distinct at the whole-genome level, we observed combinations of nonsynonymous substitutions at the CYP51 locus identical to those observed elsewhere in the world. We observed nine previously reported nonsynonymous mutations as well as isolates that carried a combination of the previously reported L50S, S188N, A379G, I381V, Y459DEL, G460DEL, and N513K substitutions. Assays for the 50% effective concentration against a subset of isolates exposed to the tebuconazole and epoxiconazole fungicides showed high levels of azole resistance. The rapid, parallel evolution of a complex CYP51 haplotype that matches a dominant European haplotype demonstrates the enormous potential for de novo resistance emergence in pathogenic fungi.IMPORTANCE Fungicides are essential to control diseases in agriculture because many crops are highly susceptible to pathogens. However, many pathogens rapidly evolve resistance to fungicides. A large body of studies have described specific mutations conferring resistance and have often made inferences about the origins of resistance based on sequencing data from the target gene alone. Here, we show the de novo acquisition of resistance to the ubiquitously used azole fungicides in genetically isolated populations of the wheat pathogen Zymoseptoria tritici in Tasmania, Australia. We confirm evidence for parallel evolution through genome-scale analyses of representative worldwide populations. The emergence of complex resistance haplotypes following a well-documented recent introduction of azoles into Australian farming practices demonstrates how rapidly chemical resistance evolves in agricultural ecosystems.

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.

Extraordinary Genome Instability and Widespread Chromosome Rearrangements During Vegetative Growth

The haploid genome of the pathogenic fungus Zymoseptoria tritici is contained on "core" and "accessory" chromosomes. While 13 core chromosomes are found in all strains, as many as eight accessory chromosomes show presence/absence variation and rearrangements among field isolates. The factors influencing these presence/absence polymorphisms are so far unknown. We investigated chromosome stability using experimental evolution, karyotyping, and genome sequencing. We report extremely high and variable rates of accessory chromosome loss during mitotic propagation in vitro and in planta Spontaneous chromosome loss was observed in 2 to >50% of cells during 4 weeks of incubation. Similar rates of chromosome loss in the closely related Zymoseptoria ardabiliae suggest that this extreme chromosome dynamic is a conserved phenomenon in the genus. Elevating the incubation temperature greatly increases instability of accessory and even core chromosomes, causing severe rearrangements involving telomere fusion and chromosome breakage. Chromosome losses do not affect the fitness of Zymoseptoria tritici in vitro, but some lead to increased virulence, suggesting an adaptive role of this extraordinary chromosome instability.

Origin and Spread of Spliceosomal Introns: Insights from the Fungal Clade Zymoseptoria

Spliceosomal introns are a key feature of eukaryote genome architecture and have been proposed to originate from selfish group II introns from an endosymbiotic bacterium, that is, the ancestor of mitochondria. However, the mechanisms underlying the wide spread of spliceosomal introns across eukaryotic genomes have been obscure. In this study, we characterize the dynamic evolution of spliceosomal introns in the fungal genus Zymoseptoria at different evolutionary scales, that is, within a genome, among conspecific strains within species, and between different species. Within the genome, spliceosomal introns can proliferate in unrelated genes and intergenic regions. Among conspecific strains, spliceosomal introns undergo rapid turnover (gains and losses) and frequent sequence exchange between geographically distinct strains. Furthermore, spliceosomal introns could undergo introgression between distinct species, which can further promote intron invasion and proliferation. The dynamic invasion and proliferation processes of spliceosomal introns resemble the life cycles of mobile selfish (group I/II) introns, and these intron movements, at least in part, account for the dramatic processes of intron gain and intron loss during eukaryotic evolution.

Distinct Trajectories of Massive Recent Gene Gains and Losses in Populations of a Microbial Eukaryotic Pathogen

Differences in gene content are a significant source of variability within species and have an impact on phenotypic traits. However, little is known about the mechanisms responsible for the most recent gene gains and losses. We screened the genomes of 123 worldwide isolates of the major pathogen of wheat Zymoseptoria tritici for robust evidence of gene copy number variation. Based on orthology relationships in three closely related fungi, we identified 599 gene gains and 1,024 gene losses that have not yet reached fixation within the focal species. Our analyses of gene gains and losses segregating in populations showed that gene copy number variation arose preferentially in subtelomeres and in proximity to transposable elements. Recently lost genes were enriched in virulence factors and secondary metabolite gene clusters. In contrast, recently gained genes encoded mostly secreted protein lacking a conserved domain. We analyzed the frequency spectrum at loci segregating a gene presence–absence polymorphism in four worldwide populations. Recent gene losses showed a significant excess in low-frequency variants compared with genome-wide single nucleotide polymorphism, which is indicative of strong negative selection against gene losses. Recent gene gains were either under weak negative selection or neutral. We found evidence for strong divergent selection among populations at individual loci segregating a gene presence–absence polymorphism. Hence, gene gains and losses likely contributed to local adaptation. Our study shows that microbial eukaryotes harbor extensive copy number variation within populations and that functional differences among recently gained and lost genes led to distinct evolutionary trajectories.

Pangenome analyses of the wheat pathogen Zymoseptoria tritici reveal the structural basis of a highly plastic eukaryotic genome

Background

Structural variation contributes substantially to polymorphism within species. Chromosomal rearrangements that impact genes can lead to functional variation among individuals and influence the expression of phenotypic traits. Genomes of fungal pathogens show substantial chromosomal polymorphism that can drive virulence evolution on host plants. Assessing the adaptive significance of structural variation is challenging, because most studies rely on inferences based on a single reference genome sequence.

Results

We constructed and analyzed the pangenome of Zymoseptoria tritici, a major pathogen of wheat that evolved host specialization by chromosomal rearrangements and gene deletions. We used single-molecule real-time sequencing and high-density genetic maps to assemble multiple genomes. We annotated the gene space based on transcriptomics data that covered the infection life cycle of each strain. Based on a total of five telomere-to-telomere genomes, we constructed a pangenome for the species and identified a core set of 9149 genes. However, an additional 6600 genes were exclusive to a subset of the isolates. The substantial accessory genome encoded on average fewer expressed genes but a larger fraction of the candidate effector genes that may interact with the host during infection. We expanded our analyses of the pangenome to a worldwide collection of 123 isolates of the same species. We confirmed that accessory genes were indeed more likely to show deletion polymorphisms and loss-of-function mutations compared to core genes.

Conclusions

The pangenome construction of a highly polymorphic eukaryotic pathogen showed that a single reference genome significantly underestimates the gene space of a species. The substantial accessory genome provides a cradle for adaptive evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-017-0457-4) contains supplementary material, which is available to authorized users.

A fungal wheat pathogen evolved host specialization by extensive chromosomal rearrangements

Fungal pathogens can rapidly evolve virulence towards resistant crops in agricultural ecosystems. Gains in virulence are often mediated by the mutation or deletion of a gene encoding a protein recognized by the plant immune system. However, the loci and the mechanisms of genome evolution enabling rapid virulence evolution are poorly understood. We performed genome-wide association mapping on a global collection of 106 strains of Zymoseptoria tritici, the most damaging pathogen of wheat in Europe, to identify polymorphisms linked to virulence on two wheat varieties. We found 25 distinct genomic loci associated with reproductive success of the pathogen. However, no locus was shared between the host genotypes, suggesting host specialization. The main locus associated with virulence encoded a highly expressed, small secreted protein. Population genomic analyses showed that the gain in virulence was explained by a segregating gene deletion polymorphism. The deletion was likely adaptive by preventing detection of the encoded protein. Comparative genomics of closely related species showed that the locus emerged de novo since speciation. A large cluster of transposable elements in direct proximity to the locus generated extensive rearrangements leading to multiple independent gene losses. Our study demonstrates that rapid turnover in the chromosomal structure of a pathogen can drive host specialization.

Rapid emergence of pathogens in agro-ecosystems: global threats to agricultural sustainability and food security

Agricultural ecosystems are composed of genetically depauperate populations of crop plants grown at a high density and over large spatial scales, with the regional composition of crop species changing little from year to year. These environments are highly conducive for the emergence and dissemination of pathogens. The uniform host populations facilitate the specialization of pathogens to particular crop cultivars and allow the build-up of large population sizes. Population genetic and genomic studies have shed light on the evolutionary mechanisms underlying speciation processes, adaptive evolution and long-distance dispersal of highly damaging pathogens in agro-ecosystems. These studies document the speed with which pathogens evolve to overcome crop resistance genes and pesticides. They also show that crop pathogens can be disseminated very quickly across and among continents through human activities. In this review, we discuss how the peculiar architecture of agro-ecosystems facilitates pathogen emergence, evolution and dispersal. We present four example pathosystems that illustrate both pathogen specialization and pathogen speciation, including different time frames for emergence and different mechanisms underlying the emergence process. Lastly, we argue for a re-design of agro-ecosystems that embraces the concept of dynamic diversity to improve their resilience to pathogens. This article is part of the themed issue 'Tackling emerging fungal threats to animal health, food security and ecosystem resilience'.

The Evolution of Orphan Regions in Genomes of a Fungal Pathogen of Wheat

Fungal plant pathogens rapidly evolve virulence on resistant hosts through mutations in genes encoding proteins that modulate the host immune responses. The mutational spectrum likely includes chromosomal rearrangements responsible for gains or losses of entire genes. However, the mechanisms creating adaptive structural variation in fungal pathogen populations are poorly understood. We used complete genome assemblies to quantify structural variants segregating in the highly polymorphic fungal wheat pathogen Zymoseptoria tritici The genetic basis of virulence in Z. tritici is complex, and populations harbor significant genetic variation for virulence; hence, we aimed to identify whether structural variation led to functional differences. We combined single-molecule real-time sequencing, genetic maps, and transcriptomics data to generate a fully assembled and annotated genome of the highly virulent field isolate 3D7. Comparative genomics analyses against the complete reference genome IPO323 identified large chromosomal inversions and the complete gain or loss of transposable-element clusters, explaining the extensive chromosomal-length polymorphisms found in this species. Both the 3D7 and IPO323 genomes harbored long tracts of sequences exclusive to one of the two genomes. These orphan regions contained 296 genes unique to the 3D7 genome and not previously known for this species. These orphan genes tended to be organized in clusters and showed evidence of mutational decay. Moreover, the orphan genes were enriched in genes encoding putative effectors and included a gene that is one of the most upregulated putative effector genes during wheat infection. Our study showed that this pathogen species harbored extensive chromosomal structure polymorphism that may drive the evolution of virulence.

IMPORTANCE

Pathogen outbreak populations often harbor previously unknown genes conferring virulence. Hence, a key puzzle of rapid pathogen evolution is the origin of such evolutionary novelty in genomes. Chromosomal rearrangements and structural variation in pathogen populations likely play a key role. However, identifying such polymorphism is challenging, as most genome-sequencing approaches only yield information about point mutations. We combined long-read technology and genetic maps to assemble the complete genome of a strain of a highly polymorphic fungal pathogen of wheat. Comparisons against the reference genome of the species showed substantial variation in the chromosome structure and revealed large regions unique to each assembled genome. These regions were enriched in genes encoding likely effector proteins, which are important components of pathogenicity. Our study showed that pathogen populations harbor extensive polymorphism at the chromosome level and that this polymorphism can be a source of adaptive genetic variation in pathogen evolution.