Preventing multi-resistance: New insights for managing fungal adaptation.

Sustainable crop protection is vital for food security, yet it is under threat due to the adaptation of a diverse and evolving pathogen population. Resistance can be managed by maximising the diversity of selection pressure through dose variation and the spatial and temporal combination of active ingredients. This study explores the interplay between operational drivers for maximising the sustainability of management strategies in relation to the resistance status of fungal populations. We applied an experimental evolution approach to three artificial populations of Zymoseptoria tritici, an economically significant wheat pathogen, each differing in initial resistance status. Our findings reveal that diversified selection pressure curtails the selection of resistance in naïve populations and those with low frequencies of single resistance. Increasing the number of modes of action most effectively delays resistance development, surpassing the increase in the number of fungicides, fungicide choice based on resistance risk, and temporal variation in fungicide exposure. However, this approach favours generalism in the evolved populations. The prior presence of multiple resistant isolates and their subsequent selection in populations override the effects of diversity in management strategies, thereby invalidating any universal ranking. Therefore, the initial resistance composition must be specifically considered in sustainable resistance management to address real-world field situations.

A thousand-genome panel retraces the global spread and adaptation of a major fungal crop pathogen.

Human activity impacts the evolutionary trajectories of many species worldwide. Global trade of agricultural goods contributes to the dispersal of pathogens reshaping their genetic makeup and providing opportunities for virulence gains. Understanding how pathogens surmount control strategies and cope with new climates is crucial to predicting the future impact of crop pathogens. Here, we address this by assembling a global thousand-genome panel of Zymoseptoria tritici, a major fungal pathogen of wheat reported in all production areas worldwide. We identify the global invasion routes and ongoing genetic exchange of the pathogen among wheat-growing regions. We find that the global expansion was accompanied by increased activity of transposable elements and weakened genomic defenses. Finally, we find significant standing variation for adaptation to new climates encountered during the global spread. Our work shows how large population genomic panels enable deep insights into the evolutionary trajectory of a major crop pathogen.

Fungicide Resistance Evolution and Detection in Plant Pathogens: Plasmopara viticola as a Case Study.

The use of single-site fungicides to control plant pathogens in the agroecosystem can be associated with an increased selection of resistance. The evolution of resistance represents one of the biggest challenges in disease control. In vineyards, frequent applications of fungicides are carried out every season for multiple years. The agronomic risk of developing fungicide resistance is, therefore, high. Plasmopara viticola, the causal agent of grapevine downy mildew, is a high risk pathogen associated with the development of fungicide resistance. P. viticola has developed resistance to most of the fungicide classes used and constitutes one of the most important threats for grapevine production. The goals of this review are to describe fungicide resistance evolution in P. viticola populations and how to conduct proper monitoring activities. Different methods have been developed for phenotyping and genotyping P. viticola for fungicide resistance and the different phases of resistance evolution and life cycles of the pathogen are discussed, to provide a full monitoring toolkit to limit the spread of resistance. A detailed revision of the available tools will help in shaping and harmonizing the monitoring activities between countries and organizations.

A dispensable paralog of succinate dehydrogenase subunit C mediates standing resistance towards a subclass of SDHI fungicides in Zymoseptoria tritici.

Succinate dehydrogenase inhibitor (SDHI) fungicides are widely used for the control of a broad range of fungal diseases. This has been the most rapidly expanding fungicide group in terms of new molecules discovered and introduced for agricultural use over the past fifteen years. A particular pattern of differential sensitivity (resistance) to the stretched heterocycle amide SDHIs (SHA-SDHIs), a subclass of chemically-related SDHIs, was observed in naïve Zymoseptoria tritici populations not previously exposed to these chemicals. Subclass-specific resistance was confirmed at the enzyme level but did not correlate with the genotypes of the succinate dehydrogenase (SDH) encoding genes. Mapping and characterization of the molecular mechanisms responsible for standing SHA-SDHI resistance in natural field isolates identified a gene paralog of SDHC, termed ZtSDHC3, which encodes for an alternative C subunit of succinate dehydrogenase, named alt-SDHC. Using reverse genetics, we showed that alt-SDHC associates with the three other SDH subunits, leading to a fully functional enzyme and that a unique Qp-site residue within the alt-SDHC protein confers SHA-SDHI resistance. Enzymatic assays, computational modelling and docking simulations for the two SQR enzymes (altC-SQR, WT_SQR) enabled us to describe enzyme-inhibitor interactions at an atomistic level and to propose rational explanations for differential potency and resistance across SHA-SDHIs. European Z. tritici populations displayed a presence (20-30%) / absence polymorphism of ZtSDHC3, as well as differences in ZtSDHC3 expression levels and splicing efficiency. These polymorphisms have a strong impact on SHA-SDHI resistance phenotypes. Characterization of the ZtSDHC3 promoter in European Z. tritici populations suggests that transposon insertions are associated with the strongest resistance phenotypes. These results establish that a dispensable paralogous gene determines SHA-SDHIs fungicide resistance in natural populations of Z. tritici. This study paves the way to an increased awareness of the role of fungicidal target paralogs in resistance to fungicides and demonstrates the paramount importance of population genomics in fungicide discovery.

Comparative Transcriptome Analyses in Zymoseptoria tritici Reveal Significant Differences in Gene Expression Among Strains During Plant Infection.

Zymoseptoria tritici is an ascomycete fungus that causes Septoria tritici blotch, a globally distributed foliar disease on wheat. Z. tritici populations are highly polymorphic and exhibit significant quantitative variation for virulence. Despite its importance, the genes responsible for quantitative virulence in this pathogen remain largely unknown. We investigated the expression profiles of four Z. tritici strains differing in virulence in an experiment conducted under uniform environmental conditions. Transcriptomes were compared at four different infection stages to characterize the regulation of gene families thought to be involved in virulence and to identify new virulence factors. The major components of the fungal infection transcriptome showed consistent expression profiles across strains. However, strain-specific regulation was observed for many genes, including some encoding putative virulence factors. We postulate that strain-specific regulation of virulence factors can determine the outcome of Z. tritici infections. We show that differences in gene expression may be major determinants of virulence variation among Z. tritici strains, adding to the already known contributions to virulence variation based on differences in gene sequence and gene presence/absence polymorphisms.

Comparative transcriptomic analyses of Zymoseptoria tritici strains show complex lifestyle transitions and intraspecific variability in transcription profiles.

Zymoseptoria tritici causes Septoria tritici blotch (STB) on wheat. The disease interaction is characterized by clearly defined temporal phases of infection, ultimately resulting in the death of host tissue. Zymoseptoria tritici is a highly polymorphic species with significant intraspecific variation in virulence profiles. We generated a deep transcriptomic sequencing dataset spanning the entire time course of an infection using a previously uncharacterized, highly virulent Z. tritici strain isolated from a Swiss wheat field. We found that seven clusters of gene transcription profiles explained the progression of the infection. The earliest highly up-regulated genes included chloroperoxidases, which may help the fungus cope with plant defences. The onset of necrotrophy was characterized by a concerted up-regulation of proteases, plant cell wall-degrading enzymes and lipases. Functions related to nutrition and growth characterized late necrotrophy and the transition to saprotrophic growth on dead plant tissue. We found that the peak up-regulation of genes essential for mating coincided with the necrotrophic phase. We performed an intraspecies comparative transcriptomics analysis using a comparable time course infection experiment of the genome reference isolate IPO323. Major components of the fungal infection transcriptome were conserved between the two strains. However, individual small, secreted proteins, proteases and cell wall-degrading enzymes showed strongly differentiated transcriptional profiles between isolates. Our analyses illustrate that successful STB infections involve complex transcriptomic remodelling to up-regulate distinct gene functions. Heterogeneity in transcriptomes among isolates may explain some of the considerable variation in virulence and host specialization found within the species.

Microbial mediation of complex subterranean mineral structures.

Helictites--an enigmatic type of mineral structure occurring in some caves--differ from classical speleothems as they develop with orientations that defy gravity. While theories for helictite formation have been forwarded, their genesis remains equivocal. Here, we show that a remarkable suite of helictites occurring in Asperge Cave (France) are formed by biologically-mediated processes, rather than abiotic processes as had hitherto been proposed. Morphological and petro-physical properties are inconsistent with mineral precipitation under purely physico-chemical control. Instead, microanalysis and molecular-biological investigation reveals the presence of a prokaryotic biofilm intimately associated with the mineral structures. We propose that microbially-influenced mineralization proceeds within a gliding biofilm which serves as a nucleation site for CaCO3, and where chemotaxis influences the trajectory of mineral growth, determining the macroscopic morphology of the speleothems. The influence of biofilms may explain the occurrence of similar speleothems in other caves worldwide, and sheds light on novel biomineralization processes.

Zymoseptoria tritici: A major threat to wheat production, integrated approaches to control.

Syngenta is one of the major agrochemical companies with enormous breadth of technologies in Crop Protection, Seeds and Seed Care. Through an exceptionally broad product range and research investment, we are not only able to provide the grower with integrated offers now but also truly innovative and transformative technologies in the future. In this commentary Syngenta scientists give their views on the key wheat pathogen Zymoseptoria tritici from its business importance in Europe, the way we screen new Z. tritici fungicides, the way we monitor the evolution of fungicide resistance and breed for Z. tritici resistance. These four points are continuously revisited and adapted during the development of new fungicides, and academic collaborations are critically important to stay at the fore front of developments in cell biology, physiology and genetic research.

Hitchhiking selection is driving intron gain in a pathogenic fungus.

The variability of intron density among eukaryotes is puzzling and still debated. Most previous studies have been limited because of the near absence of intron presence-absence polymorphism (IPAP) within species or because comparisons could be made only between distantly related species. We conducted population genetic analyses on eight loci showing IPAP to investigate the effect of natural selection on intron dynamics in a global collection of the panmictic fungal plant pathogen Zymoseptoria tritici and its very close relatives. Five of these loci likely represent recent intron gains because their absence is fixed among the closest relatives of Z. tritici, and three likely represent recent intron losses because their presence is fixed among the close relatives. We analyzed signatures of selection by comparing allele frequencies, nucleotide diversities, and rates of recombination and found compelling evidence that at least two out of the five intron-gain loci, a SWIM zinc-finger gene and a sugar transporter, are under directional selection favoring alleles that gained the intron. Our results suggest that the intron-present alleles of these loci are sweeping to fixation, providing a genetic hitchhiking mechanism to explain rapid intron gain in Z. tritici. The overall findings are consistent with the hypothesis that intron gains are more likely to be driven by selection while intron losses are more likely to be due to neutral processes such as genetic drift.

Comparative analysis of mitochondrial genomes from closely related Rhynchosporium species reveals extensive intron invasion.

We sequenced and annotated the complete mitochondrial (mt) genomes of four closely related Rhynchosporium species that diverged ∼14,000-35,000years ago. During this time frame, three of the mt genomes expanded significantly due to an invasion of introns into three genes (cox1, cox2, and nad5). The enlarged mt genomes contained ∼40% introns compared to 8.1% in uninvaded relatives. Many intron gains were accompanied by co-conversion of flanking exonic regions. The comparative analysis revealed a highly variable set of non-intronic, free-standing ORFs of unknown function (uORFs). This is consistent with a rapidly evolving accessory compartment in the mt genome of these closely related species. Only one free-standing uORF was shared among all mt genomes analyzed. This uORF had a mutation rate similar to the core mt protein-encoding genes, suggesting conservation of function among the species. The nucleotide composition of the core protein-encoding genes significantly differed from those of introns and uORFs. The mt mutation rate was 77 times higher than the nuclear mutation rate, indicating that the phylogeny inferred from mt genes may better resolve the phylogenetic relationships among closely related Rhynchosporium species than phylogenies inferred from nuclear genes.

Coevolution and life cycle specialization of plant cell wall degrading enzymes in a hemibiotrophic pathogen.

Zymoseptoria tritici is an important fungal pathogen on wheat that originated in the Fertile Crescent. Its closely related sister species Z. pseudotritici and Z. ardabiliae infect wild grasses in the same region. This recently emerged host-pathogen system provides a rare opportunity to investigate the evolutionary processes shaping the genome of an emerging pathogen. Here, we investigate genetic signatures in plant cell wall degrading enzymes (PCWDEs) that are likely affected by or driving coevolution in plant-pathogen systems. We hypothesize four main evolutionary scenarios and combine comparative genomics, transcriptomics, and selection analyses to assign the majority of PCWDEs in Z. tritici to one of these scenarios. We found widespread differential transcription among different members of the same gene family, challenging the idea of functional redundancy and suggesting instead that specialized enzymatic activity occurs during different stages of the pathogen life cycle. We also find that natural selection has significantly affected at least 19 of the 48 identified PCWDEs. The majority of genes showed signatures of purifying selection, typical for the scenario of conserved substrate optimization. However, six genes showed diversifying selection that could be attributed to either host adaptation or host evasion. This study provides a powerful framework to better understand the roles played by different members of multigene families and to determine which genes are the most appropriate targets for wet laboratory experimentation, for example, to elucidate enzymatic function during relevant phases of a pathogen's life cycle.

Evidence for extensive recent intron transposition in closely related fungi.

Though spliceosomal introns are a major structural component of most eukaryotic genes and intron density varies by more than three orders of magnitude among eukaryotes [1-3], the origins of introns are poorly understood, and only a few cases of unambiguous intron gain are known [4-8]. We utilized population genomic comparisons of three closely related fungi to identify crucial transitory phases of intron gain and loss. We found 74 intron positions showing intraspecific presence-absence polymorphisms (PAPs) for the entire intron. Population genetic analyses identified intron PAPs at different stages of fixation and showed that intron gain or loss was very recent. We found direct support for extensive intron transposition among unrelated genes. A substantial proportion of highly similar introns in the genome either were recently gained or showed a transient phase of intron PAP. We also identified an intron transfer among paralogous genes that created a new intron. Intron loss was due mainly to homologous recombination involving reverse-transcribed mRNA. The large number of intron positions in transient phases of either intron gain or loss shows that intron evolution is much faster than previously thought and provides an excellent model to study molecular mechanisms of intron gain.

Finished genome of the fungal wheat pathogen Mycosphaerella graminicola reveals dispensome structure, chromosome plasticity, and stealth pathogenesis.

The plant-pathogenic fungus Mycosphaerella graminicola (asexual stage: Septoria tritici) causes septoria tritici blotch, a disease that greatly reduces the yield and quality of wheat. This disease is economically important in most wheat-growing areas worldwide and threatens global food production. Control of the disease has been hampered by a limited understanding of the genetic and biochemical bases of pathogenicity, including mechanisms of infection and of resistance in the host. Unlike most other plant pathogens, M. graminicola has a long latent period during which it evades host defenses. Although this type of stealth pathogenicity occurs commonly in Mycosphaerella and other Dothideomycetes, the largest class of plant-pathogenic fungi, its genetic basis is not known. To address this problem, the genome of M. graminicola was sequenced completely. The finished genome contains 21 chromosomes, eight of which could be lost with no visible effect on the fungus and thus are dispensable. This eight-chromosome dispensome is dynamic in field and progeny isolates, is different from the core genome in gene and repeat content, and appears to have originated by ancient horizontal transfer from an unknown donor. Synteny plots of the M. graminicola chromosomes versus those of the only other sequenced Dothideomycete, Stagonospora nodorum, revealed conservation of gene content but not order or orientation, suggesting a high rate of intra-chromosomal rearrangement in one or both species. This observed "mesosynteny" is very different from synteny seen between other organisms. A surprising feature of the M. graminicola genome compared to other sequenced plant pathogens was that it contained very few genes for enzymes that break down plant cell walls, which was more similar to endophytes than to pathogens. The stealth pathogenesis of M. graminicola probably involves degradation of proteins rather than carbohydrates to evade host defenses during the biotrophic stage of infection and may have evolved from endophytic ancestors.

Evolutionary history of the mitochondrial genome in Mycosphaerella populations infecting bread wheat, durum wheat and wild grasses.

Plant pathogens emerge in agro-ecosystems following different evolutionary mechanisms over different time scales. Previous analyses based on sequence variation at six nuclear loci indicated that Mycosphaerella graminicola diverged from an ancestral population adapted to wild grasses during the process of wheat domestication approximately 10,500 years ago. We tested this hypothesis by conducting coalescence analyses based on four mitochondrial loci using 143 isolates that included four closely related pathogen species originating from four continents. Pathogen isolates from bread and durum wheat were included to evaluate the emergence of specificity towards these hosts in M. graminicola. Although mitochondrial and nuclear genomes differed greatly in degree of genetic variability, their coalescence was remarkably congruent, supporting the proposed origin of M. graminicola through host tracking. The coalescence analysis was unable to trace M. graminicola host specificity through recent evolutionary time, indicating that the specificity towards durum or bread wheat emerged following the domestication of the pathogen on wheat.

Intraspecific comparison and annotation of two complete mitochondrial genome sequences from the plant pathogenic fungus Mycosphaerella graminicola.

The mitochondrial genomes of two isolates of the wheat pathogen Mycosphaerella graminicola were sequenced completely and compared to identify polymorphic regions. This organism is of interest because it is phylogenetically distant from other fungi with sequenced mitochondrial genomes and it has shown discordant patterns of nuclear and mitochondrial diversity. The mitochondrial genome of M. graminicola is a circular molecule of approximately 43,960bp containing the typical genes coding for 14 proteins related to oxidative phosphorylation, one RNA polymerase, two rRNA genes and a set of 27 tRNAs. The mitochondrial DNA of M. graminicola lacks the gene encoding the putative ribosomal protein (rps5-like), commonly found in fungal mitochondrial genomes. Most of the tRNA genes were clustered with a gene order conserved with many other ascomycetes. A sample of 35 additional strains representing the known global mt diversity was partially sequenced to measure overall mitochondrial variability within the species. Little variation was found, confirming previous RFLP-based findings of low mitochondrial diversity. The mitochondrial sequence of M. graminicola is the first reported from the family Mycosphaerellaceae or the order Capnodiales. The sequence also provides a tool to better understand the development of fungicide resistance and the conflicting pattern of high nuclear and low mitochondrial diversity in global populations of this fungus.

Dothideomycete plant interactions illuminated by genome sequencing and EST analysis of the wheat pathogen Stagonospora nodorum.

Stagonospora nodorum is a major necrotrophic fungal pathogen of wheat (Triticum aestivum) and a member of the Dothideomycetes, a large fungal taxon that includes many important plant pathogens affecting all major crop plant families. Here, we report the acquisition and initial analysis of a draft genome sequence for this fungus. The assembly comprises 37,164,227 bp of nuclear DNA contained in 107 scaffolds. The circular mitochondrial genome comprises 49,761 bp encoding 46 genes, including four that are intron encoded. The nuclear genome assembly contains 26 classes of repetitive DNA, comprising 4.5% of the genome. Some of the repeats show evidence of repeat-induced point mutations consistent with a frequent sexual cycle. ESTs and gene prediction models support a minimum of 10,762 nuclear genes. Extensive orthology was found between the polyketide synthase family in S. nodorum and Cochliobolus heterostrophus, suggesting an ancient origin and conserved functions for these genes. A striking feature of the gene catalog was the large number of genes predicted to encode secreted proteins; the majority has no meaningful similarity to any other known genes. It is likely that genes for host-specific toxins, in addition to ToxA, will be found among this group. ESTs obtained from axenic mycelium grown on oleate (chosen to mimic early infection) and late-stage lesions sporulating on wheat leaves were obtained. Statistical analysis shows that transcripts encoding proteins involved in protein synthesis and in the production of extracellular proteases, cellulases, and xylanases predominate in the infection library. This suggests that the fungus is dependant on the degradation of wheat macromolecular constituents to provide the carbon skeletons and energy for the synthesis of proteins and other components destined for the developing pycnidiospores.

Significant difference in pathogenicity between MAT1-1 and MAT1-2 isolates in the wheat pathogen Mycosphaerella graminicola.

Five Mycosphaerella graminicola populations from four geographic regions (Australia, Israel, Switzerland, and the USA) were assayed for neutral RFLP markers and mating type idiomorphs. On average, 25-30 genetically distinct isolates were selected from each population and their pathogenicity was measured on two wheat cultivars in a common garden experiment conducted in a greenhouse. A significant difference in pathogenicity was found between MAT1-1 and MAT1-2 isolates. On average, MAT1-1 isolates had 14-22% greater pathogenicity than MAT1-2 isolates. The pattern of higher pathogenicity in MAT1-1 isolates was consistent across four geographical populations and on two wheat cultivars. A uniform and continuous variation in pathogenicity was found among isolates within each mating type, but no genetic differentiation in selectively neutral RFLP loci was found between mating types, consistent with the hypothesis that differences in pathogenicity were not due to the effects of specific pathogenicity genes or non-random genetic backgrounds.