A fungal avirulence factor encoded in a highly plastic genomic region triggers partial resistance to septoria tritici blotch

The Zymoseptoria tritici Avirulence Factor AvrStb6 Accumulates in Hyphae Close to Stomata and Triggers a Wheat Defense Response Hindering Fungal Penetration

Zymoseptoria tritici, the causal agent of Septoria tritici blotch, is one of Europe's most damaging wheat pathogens, causing significant economic losses. Genetic resistance is a common strategy to control the disease, Stb6 being a resistance gene used for more than 100 years in Europe. This study investigates the molecular mechanisms underlying Stb6-mediated resistance. Utilizing confocal microscopy imaging, we determined that Z. tritici epiphytic hyphae mainly accumulate the corresponding avirulence factor AvrStb6 in close proximity to stomata. Consequently, the progression of AvrStb6-expressing avirulent strains is hampered during penetration. The fungal growth inhibition co-occurs with a transcriptional reprogramming in wheat characterized by an induction of immune responses, genes involved in stomatal regulation, and cell wall-related genes. Overall, we shed light on the gene-for-gene resistance mechanisms in the wheat-Z. tritici pathosystem at the cytological and transcriptomic level, and our results highlight that stomatal penetration is a critical process for pathogenicity and resistance. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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.

Quantitative pathogenicity and host adaptation in a fungal plant pathogen revealed by whole-genome sequencing

Knowledge of genetic determinism and evolutionary dynamics mediating host-pathogen interactions is essential to manage fungal plant diseases. Studies on the genetic architecture of fungal pathogenicity often focus on large-effect effector genes triggering strong, qualitative resistance. It is not clear how this translates to predominately quantitative interactions. Here, we use the Zymoseptoria tritici-wheat model to elucidate the genetic architecture of quantitative pathogenicity and mechanisms mediating host adaptation. With a multi-host genome-wide association study, we identify 19 high-confidence candidate genes associated with quantitative pathogenicity. Analysis of genetic diversity reveals that sequence polymorphism is the main evolutionary process mediating differences in quantitative pathogenicity, a process that is likely facilitated by genetic recombination and transposable element dynamics. Finally, we use functional approaches to confirm the role of an effector-like gene and a methyltransferase in phenotypic variation. This study highlights the complex genetic architecture of quantitative pathogenicity, extensive diversifying selection and plausible mechanisms facilitating pathogen adaptation.

Multiplexed effector screening for recognition by endogenous resistance genes using positive defense reporters in wheat protoplasts

Plant resistance (R) and pathogen avirulence (Avr) gene interactions play a vital role in pathogen resistance. Efficient molecular screening tools for crops lack far behind their model organism counterparts, yet they are essential to rapidly identify agriculturally important molecular interactions that trigger host resistance. Here, we have developed a novel wheat protoplast assay that enables efficient screening of Avr/R interactions at scale. Our assay allows access to the extensive gene pool of phenotypically described R genes because it does not require the overexpression of cloned R genes. It is suitable for multiplexed Avr screening, with interactions tested in pools of up to 50 Avr candidates. We identified Avr/R-induced defense genes to create a promoter-luciferase reporter. Then, we combined this with a dual-color ratiometric reporter system that normalizes read-outs accounting for experimental variability and Avr/R-induced cell death. Moreover, we introduced a self-replicative plasmid reducing the amount of plasmid used in the assay. Our assay increases the throughput of Avr candidate screening, accelerating the study of cellular defense signaling and resistance gene identification in wheat. We anticipate that our assay will significantly accelerate Avr identification for many wheat pathogens, leading to improved genome-guided pathogen surveillance and breeding of disease-resistant crops.

GWAS reveals a rapidly evolving candidate avirulence effector in the Cercospora leaf spot pathogen

The major resistance gene BvCR4 recently bred into sugar beet hybrids provides a high level of resistance to Cercospora leaf spot caused by the fungal pathogen Cercospora beticola. The occurrence of pathogen strains that overcome BvCR4 was studied using field trials in Switzerland conducted under natural disease pressure. Virulence of a subset of these strains was evaluated in a field trial conducted under elevated artificial disease pressure. We created a new C. beticola reference genome and mapped whole genome sequences of 256 isolates collected in Switzerland and Germany. These were combined with virulence phenotypes to conduct three separate genome-wide association studies (GWAS) to identify candidate avirulence genes. We identified a locus associated with avirulence containing a putative avirulence effector gene named AvrCR4. All virulent isolates either lacked AvrCR4 or had nonsynonymous mutations within the gene. AvrCR4 was present in all 74 isolates from non-BvCR4 hybrids, whereas 33 of 89 isolates from BvCR4 hybrids carried a deletion. We also mapped genomic data from 190 publicly available US isolates to our new reference genome. The AvrCR4 deletion was found in only one of 95 unique isolates from non-BvCR4 hybrids in the United States. AvrCR4 presents a unique example of an avirulence effector in which virulent alleles have only recently emerged. Most likely these were selected out of standing genetic variation after deployment of BvCR4. Identification of AvrCR4 will enable real-time screening of C. beticola populations for the emergence and spread of virulent isolates.

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.

Combined reference-free and multi-reference based GWAS uncover cryptic variation underlying rapid adaptation in a fungal plant pathogen

Microbial pathogens often harbor substantial functional diversity driven by structural genetic variation. Rapid adaptation from such standing variation threatens global food security and human health. Genome-wide association studies (GWAS) provide a powerful approach to identify genetic variants underlying recent pathogen adaptation. However, the reliance on single reference genomes and single nucleotide polymorphisms (SNPs) obscures the true extent of adaptive genetic variation. Here, we show quantitatively how a combination of multiple reference genomes and reference-free approaches captures substantially more relevant genetic variation compared to single reference mapping. We performed reference-genome based association mapping across 19 reference-quality genomes covering the diversity of the species. We contrasted the results with a reference-free (i.e., k-mer) approach using raw whole-genome sequencing data in a panel of 145 strains collected across the global distribution range of the fungal wheat pathogen Zymoseptoria tritici. We mapped the genetic architecture of 49 life history traits including virulence, reproduction and growth in multiple stressful environments. The inclusion of additional reference genome SNP datasets provides a nearly linear increase in additional loci mapped through GWAS. Variants detected through the k-mer approach explained a higher proportion of phenotypic variation than a reference genome-based approach and revealed functionally confirmed loci that classic GWAS approaches failed to map. The power of GWAS in microbial pathogens can be significantly enhanced by comprehensively capturing structural genetic variation. Our approach is generalizable to a large number of species and will uncover novel mechanisms driving rapid adaptation of pathogens.

Virulent strains of Zymoseptoria tritici suppress the host immune response and facilitate the success of avirulent strains in mixed infections

Plants interact with a plethora of pathogenic microorganisms in nature. Pathogen-plant interaction experiments focus mainly on single-strain infections, typically ignoring the complexity of multi-strain infections even though mixed infections are common and critical for the infection outcome. The wheat pathogen Zymoseptoria tritici forms highly diverse fungal populations in which several pathogen strains often colonize the same leaf. Despite the importance of mixed infections, the mechanisms governing interactions between a mixture of pathogen strains within a plant host remain largely unexplored. Here we demonstrate that avirulent pathogen strains benefit from being in mixed infections with virulent strains. We show that virulent strains suppress the wheat immune response, allowing avirulent strains to colonize the apoplast and to reproduce. Our experiments indicate that virulent strains in mixed infections can suppress the plant immune system, probably facilitating the persistence of avirulent pathogen strains in fields planted with resistant host plants.

Sas3-mediated histone acetylation regulates effector gene activation in a fungal plant pathogen

IMPORTANCE

Pathogen infections require the production of effectors that enable host colonization. Effectors have diverse functions and are only expressed at certain stages of the infection cycle. Thus, effector genes are tightly regulated by several mechanisms, including chromatin remodeling. Here, we investigate the role of histone acetylation in effector gene activation in the fungal wheat pathogen Zymoseptoria tritici. We demonstrate that lysine acetyltransferases (KATs) are essential for the spatiotemporal regulation of effector genes. We show that the KAT Sas3 is involved in leaf symptom development and pycnidia formation. Importantly, our results indicate that Sas3 controls histone acetylation of effector loci and is a regulator of effector gene activation during stomatal penetration. Overall, our work demonstrates the key role of histone acetylation in regulating gene expression associated with plant infection.

Epiphytic proliferation of Zymoseptoria tritici isolates on resistant wheat leaves

The wheat pathogen Zymoseptoria tritici is capable of a long period of pre-invasive epiphytic growth. Studies have shown that virulent isolates vary in the extent, duration and growth form of this epiphytic growth, and the fungus has been observed to undergo behaviours such as asexual reproduction by budding and vegetative fusion of hyphae on the leaf surface. This epiphytic colonisation has been investigated very little during interactions in which an isolate of Z. tritici is unable to colonise the apoplast, as occurs during avirulence. However, avirulent isolates have been seen to undergo sexual crosses in the absense of leaf penetration, and it is widely accepted that the main point of distinction between virulent and avirulent isolates occurs at the point of attempted leaf penetration or attempted apoplastic growth, which fails in the avirulent case. In this work, we describe extensive epiphytic growth in three isolates which are unable or have very limited ability to invade the leaf, and show that growth form is as variable as for fully virulent isolates. We demonstrate that during certain interactions, Z. tritici isolates rarely invade the leaf and form pycnidia, but induce necrosis. These isolates are able to achieve higher epiphytic biomass than fully virulent isolates during asymptomatic growth, and may undergo very extensive asexual reproduction on the leaf surface. These findings have implications for open questions such as whether and how Z. tritici obtains nutrients on the leaf surface and the nature of its interaction with wheat defences.

A transcription factor and a phosphatase regulate temperature-dependent morphogenesis in the fungal plant pathogen Zymoseptoria tritici

Naturally fluctuating temperatures provide a constant environmental stress that requires adaptation. Some fungal pathogens respond to heat stress by producing new morphotypes that maximize their overall fitness. The fungal wheat pathogen Zymoseptoria tritici responds to heat stress by switching from its yeast-like blastospore form to hyphae or chlamydospores. The regulatory mechanisms underlying this switch are unknown. Here, we demonstrate that a differential heat stress response is ubiquitous in Z. tritici populations around the world. We used QTL mapping to identify a single locus associated with the temperature-dependent morphogenesis and we found two genes, the transcription factor ZtMsr1 and the protein phosphatase ZtYvh1, regulating this mechanism. We find that ZtMsr1 regulates repression of hyphal growth and induces chlamydospore formation whereas ZtYvh1 is required for hyphal growth. We next showed that chlamydospore formation is a response to the intracellular osmotic stress generated by the heat stress. This intracellular stress stimulates the cell wall integrity (CWI) and high-osmolarity glycerol (HOG) MAPK pathways resulting in hyphal growth. If cell wall integrity is compromised, however, ZtMsr1 represses the hyphal development program and may induce the chlamydospore-inducing genes as a stress-response survival strategy. Taken together, these results suggest a novel mechanism through which morphological transitions are orchestrated in Z. tritici - a mechanism that may also be present in other pleomorphic fungi.

A secreted protease-like protein in Zymoseptoria tritici is responsible for avirulence on Stb9 resistance gene in wheat

Zymoseptoria tritici is the fungal pathogen responsible for Septoria tritici blotch on wheat. Disease outcome in this pathosystem is partly determined by isolate-specific resistance, where wheat resistance genes recognize specific fungal factors triggering an immune response. Despite the large number of known wheat resistance genes, fungal molecular determinants involved in such cultivar-specific resistance remain largely unknown. We identified the avirulence factor AvrStb9 using association mapping and functional validation approaches. Pathotyping AvrStb9 transgenic strains on Stb9 cultivars, near isogenic lines and wheat mapping populations, showed that AvrStb9 interacts with Stb9 resistance gene, triggering an immune response. AvrStb9 encodes an unusually large avirulence gene with a predicted secretion signal and a protease domain. It belongs to a S41 protease family conserved across different filamentous fungi in the Ascomycota class and may constitute a core effector. AvrStb9 is also conserved among a global Z. tritici population and carries multiple amino acid substitutions caused by strong positive diversifying selection. These results demonstrate the contribution of an 'atypical' conserved effector protein to fungal avirulence and the role of sequence diversification in the escape of host recognition, adding to our understanding of host-pathogen interactions and the evolutionary processes underlying pathogen adaptation.

Quantitative and qualitative plant-pathogen interactions call upon similar pathogenicity genes with a spectrum of effects

Septoria leaf blotch is a foliar wheat disease controlled by a combination of plant genetic resistances and fungicides use. R-gene-based qualitative resistance durability is limited due to gene-for-gene interactions with fungal avirulence (Avr) genes. Quantitative resistance is considered more durable but the mechanisms involved are not well documented. We hypothesize that genes involved in quantitative and qualitative plant-pathogen interactions are similar. A bi-parental population of Zymoseptoria tritici was inoculated on wheat cultivar 'Renan' and a linkage analysis performed to map QTL. Three pathogenicity QTL, Qzt-I05-1, Qzt-I05-6 and Qzt-I07-13, were mapped on chromosomes 1, 6 and 13 in Z. tritici, and a candidate pathogenicity gene on chromosome 6 was selected based on its effector-like characteristics. The candidate gene was cloned by Agrobacterium tumefaciens-mediated transformation, and a pathology test assessed the effect of the mutant strains on 'Renan'. This gene was demonstrated to be involved in quantitative pathogenicity. By cloning a newly annotated quantitative-effect gene in Z. tritici that is effector-like, we demonstrated that genes underlying pathogenicity QTL can be similar to Avr genes. This opens up the previously probed possibility that 'gene-for-gene' underlies not only qualitative but also quantitative plant-pathogen interactions in this pathosystem.

Natural variation in Avr3D1 from Zymoseptoria sp. contributes to quantitative gene-for-gene resistance and to host specificity

Successful host colonization by plant pathogens requires the circumvention of host defense responses, frequently through sequence modifications in secreted pathogen proteins known as avirulence factors (Avrs). Although Avr sequences are often polymorphic, the contribution of these polymorphisms to virulence diversity in natural pathogen populations remains largely unexplored. We used molecular genetic tools to determine how natural sequence polymorphisms of the avirulence factor Avr3D1 in the wheat pathogen Zymoseptoria tritici contributed to adaptive changes in virulence. We showed that there is a continuous distribution in the magnitude of resistance triggered by different Avr3D1 isoforms and demonstrated that natural variation in an Avr gene can lead to a quantitative resistance phenotype. We further showed that homologues of Avr3D1 in two nonpathogenic sister species of Z. tritici are recognized by some wheat cultivars, suggesting that Avr-R gene-for-gene interactions can contribute to nonhost resistance. We suggest that the mechanisms underlying host range, qualitative resistance, and quantitative resistance are not exclusive.

Genetics and breeding for resistance against four leaf spot diseases in wheat (Triticum aestivum L.)

In wheat, major yield losses are caused by a variety of diseases including rusts, spike diseases, leaf spot and root diseases. The genetics of resistance against all these diseases have been studied in great detail and utilized for breeding resistant cultivars. The resistance against leaf spot diseases caused by each individual necrotroph/hemi-biotroph involves a complex system involving resistance (R) genes, sensitivity (S) genes, small secreted protein (SSP) genes and quantitative resistance loci (QRLs). This review deals with resistance for the following four-leaf spot diseases: (i) Septoria nodorum blotch (SNB) caused by Parastagonospora nodorum; (ii) Tan spot (TS) caused by Pyrenophora tritici-repentis; (iii) Spot blotch (SB) caused by Bipolaris sorokiniana and (iv) Septoria tritici blotch (STB) caused by Zymoseptoria tritici.

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.

Recent transposable element bursts are associated with the proximity to genes in a fungal plant pathogen

The activity of transposable elements (TEs) contributes significantly to pathogen genome evolution. TEs often destabilize genome integrity but may also confer adaptive variation in pathogenicity or resistance traits. De-repression of epigenetically silenced TEs often initiates bursts of transposition activity that may be counteracted by purifying selection and genome defenses. However, how these forces interact to determine the expansion routes of TEs within a pathogen species remains largely unknown. Here, we analyzed a set of 19 telomere-to-telomere genomes of the fungal wheat pathogen Zymoseptoria tritici. Phylogenetic reconstruction and ancestral state estimates of individual TE families revealed that TEs have undergone distinct activation and repression periods resulting in highly uneven copy numbers between genomes of the same species. Most TEs are clustered in gene poor niches, indicating strong purifying selection against insertions near coding sequences, or as a consequence of insertion site preferences. TE families with high copy numbers have low sequence divergence and strong signatures of defense mechanisms (i.e., RIP). In contrast, small non-autonomous TEs (i.e., MITEs) are less impacted by defense mechanisms and are often located in close proximity to genes. Individual TE families have experienced multiple distinct burst events that generated many nearly identical copies. We found that a Copia element burst was initiated from recent copies inserted substantially closer to genes compared to older copies. Overall, TE bursts tended to initiate from copies in GC-rich niches that escaped inactivation by genomic defenses. Our work shows how specific genomic environments features provide triggers for TE proliferation in pathogen genomes.

A highly multiplexed assay to monitor pathogenicity, fungicide resistance and gene flow in the fungal wheat pathogen Zymoseptoria tritici

Crop pathogens pose severe risks to global food production due to the rapid rise of resistance to pesticides and host resistance breakdowns. Predicting future risks requires monitoring tools to identify changes in the genetic composition of pathogen populations. Here we report the design of a microfluidics-based amplicon sequencing assay to multiplex 798 loci targeting virulence and fungicide resistance genes, and randomly selected genome-wide markers for the fungal pathogen Zymoseptoria tritici. The fungus causes one of the most devastating diseases on wheat showing rapid adaptation to fungicides and host resistance. We optimized the primer design by integrating polymorphism data from 632 genomes of the same species. To test the performance of the assay, we genotyped 192 samples in two replicates. Analysis of the short-read sequence data generated by the assay showed a fairly stable success rate across samples to amplify a large number of loci. The performance was consistent between samples originating from pure genomic DNA as well as material extracted directly from infected wheat leaves. In samples with mixed genotypes, we found that the assay recovers variations in allele frequencies. We also explored the potential of the amplicon assay to recover transposable element insertion polymorphism relevant for fungicide resistance. As a proof-of-concept, we show that the assay recovers the pathogen population structure across French wheat fields. Genomic monitoring of crop pathogens contributes to more sustainable crop protection and yields.

Magnaporthe oryzae pathotype Triticum (MoT) can act as a heterologous expression system for fungal effectors with high transcript abundance in wheat

Plant pathogens deliver effector proteins to reprogramme a host plants circuitry, supporting their own growth and development, whilst thwarting defence responses. A subset of these effectors are termed avirulence factors (Avr) and can be recognised by corresponding host resistance (R) proteins, creating a strong evolutionary pressure on pathogen Avr effectors that favours their modification/deletion to evade the immune response. Hence, identifying Avr effectors and tracking their allele frequencies in a population is critical for understanding the loss of host recognition. However, the current systems available to confirm Avr effector function, particularly for obligate biotrophic fungi, remain limited and challenging. Here, we explored the utility of the genetically tractable wheat blast pathogen Magnaporthe oryzae pathotype Triticum (MoT) as a suitable heterologous expression system in wheat. Using the recently confirmed wheat stem rust pathogen (Puccina graminis f. sp. tritici) avirulence effector AvrSr50 as a proof-of-concept, we found that delivery of AvrSr50 via MoT could elicit a visible Sr50-dependant cell death phenotype. However, activation of Sr50-mediated cell death correlated with a high transgene copy number and transcript abundance in MoT transformants. This illustrates that MoT can act as an effective heterologous delivery system for fungal effectors from distantly related fungal species, but only when enough transgene copies and/or transcript abundance is achieved.

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.

Automated Quantitative Measurement of Yellow Halos Suggests Activity of Necrotrophic Effectors in Septoria tritici Blotch

Many necrotrophic plant pathogens utilize host-selective toxins or necrotrophic effectors during the infection process. We hypothesized that the chlorotic yellow halos frequently observed around necrotic lesions caused by the wheat pathogen Zymoseptoria tritici could result from the activity of necrotrophic effectors interacting with the products of toxin sensitivity genes. As an initial step toward testing this hypothesis, we developed an automated image analysis (AIA) workflow that could quantify the degree of yellow halo formation occurring in wheat leaves naturally infected by a highly diverse pathogen population under field conditions. This AIA based on statistical learning was applied to more than 10,000 naturally infected leaves collected from 335 wheat cultivars grown in a replicated field experiment. We estimated a high heritability (h² = 0.71) for the degree of yellow halo formation, suggesting that this quantitative trait has a significant genetic component. Using genome-wide association mapping, we identified six chromosome segments significantly associated with the yellow halo phenotype. Most of these segments contained candidate genes associated with targets of necrotrophic effectors in other necrotrophic pathogens. Our findings conform with the hypothesis that toxin sensitivity genes could account for a significant fraction of the observed variation in quantitative resistance to Septoria tritici blotch. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Heterochromatin in the fungal plant pathogen, Zymoseptoria tritici: Control of transposable elements, genome plasticity and virulence

Heterochromatin is a repressive chromatin state that plays key roles in the functional organisation of eukaryotic genomes. In fungal plant pathogens, effector genes that are required for host colonization tend to be associated with heterochromatic regions of the genome that are enriched with transposable elements. It has been proposed that the heterochromatin environment silences effector genes in the absence of host and dynamic chromatin remodelling facilitates their expression during infection. Here we discuss this model in the context of the key wheat pathogen, Zymoseptoria tritici. We cover progress in understanding the deposition and recognition of heterochromatic histone post translational modifications in Z. tritici and the role that heterochromatin plays in control of genome plasticity and virulence.

Quantitative Genetic Analysis of Interactions in the Pepper-Phytophthora capsici Pathosystem

The development of pepper cultivars with durable resistance to the oomycete Phytophthora capsici has been challenging due to differential interactions between the species that allow certain pathogen isolates to cause disease on otherwise resistant host genotypes. Currently, little is known about the pathogen genes involved in these interactions. To investigate the genetic basis of P. capsici virulence on individual pepper genotypes, we inoculated sixteen pepper accessions, representing commercial varieties, sources of resistance, and host differentials, with 117 isolates of P. capsici, for a total of 1,864 host-pathogen combinations. Analysis of disease outcomes revealed a significant effect of inter-species genotype-by-genotype interactions, although these interactions were quantitative rather than qualitative in scale. Isolates were classified into five pathogen subpopulations, as determined by their genotypes at over 60,000 single-nucleotide polymorphisms (SNPs). While absolute virulence levels on certain pepper accessions significantly differed between subpopulations, a multivariate phenotype reflecting relative virulence levels on certain pepper genotypes compared with others showed the strongest association with pathogen subpopulation. A genome-wide association study (GWAS) identified four pathogen loci significantly associated with virulence, two of which colocalized with putative RXLR effector genes and another with a polygalacturonase gene cluster. All four loci appeared to represent broad-spectrum virulence genes, as significant SNPs demonstrated consistent effects regardless of the host genotype tested. Host genotype-specific virulence variants in P. capsici may be difficult to map via GWAS with all but excessively large sample sizes, perhaps controlled by genes of small effect or by multiple allelic variants that have arisen independently. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Zymoseptoria tritici white-collar complex integrates light, temperature and plant cues to initiate dimorphism and pathogenesis

Transitioning from spores to hyphae is pivotal to host invasion by the plant pathogenic fungus Zymoseptoria tritici. This dimorphic switch can be initiated by high temperature in vitro (~27 °C); however, such a condition may induce cellular heat stress, questioning its relevance to field infections. Here, we study the regulation of the dimorphic switch by temperature and other factors. Climate data from wheat-growing areas indicate that the pathogen sporadically experiences high temperatures such as 27 °C during summer months. However, using a fluorescent dimorphic switch reporter (FDR1) in four wild-type strains, we show that dimorphic switching already initiates at 15-18 °C, and is enhanced by wheat leaf surface compounds. Transcriptomics reveals 1261 genes that are up- or down-regulated in hyphae of all strains. These pan-strain core dimorphism genes (PCDGs) encode known effectors, dimorphism and transcription factors, and light-responsive proteins (velvet factors, opsins, putative blue light receptors). An FDR1-based genetic screen reveals a crucial role for the white-collar complex (WCC) in dimorphism and virulence, mediated by control of PCDG expression. Thus, WCC integrates light with biotic and abiotic cues to orchestrate Z. tritici infection.

Asexual reproductive potential trumps virulence as a predictor of competitive ability in mixed infections

Natural infections frequently involve several co-infecting pathogen strains. These mixed infections can affect the extent of the infection, the transmission success of the pathogen and the eventual epidemic outcome. To date, few studies have investigated how mixed infections affect transmission between hosts. Zymoseptoria tritici is a highly diverse wheat pathogen in which multiple strains often coexist in the same lesion. Here we demonstrate that the most competitive strains often exclude their competitors during serial passages of mixed infections. The outcome of the competition depended on both the host genotype and the genotypes of the competing pathogen strains. Differences in virulence among the strains were not associated with competitive advantages during transmission, while differences in reproductive potential had a strong effect on strain competitive ability. Overall, our findings suggest that host specialization is determined mainly by the ability to successfully transmit offspring to new hosts during mixed infections.

Soil composition and plant genotype determine benzoxazinoid-mediated plant-soil feedbacks in cereals

Plant-soil feedbacks refer to effects on plants that are mediated by soil modifications caused by the previous plant generation. Maize conditions the surrounding soil by secretion of root exudates including benzoxazinoids (BXs), a class of bioactive secondary metabolites. Previous work found that a BX-conditioned soil microbiota enhances insect resistance while reducing biomass in the next generation of maize plants. Whether these BX-mediated and microbially driven feedbacks are conserved across different soils and response species is unknown. We found the BX-feedbacks on maize growth and insect resistance conserved between two arable soils, but absent in a more fertile grassland soil, suggesting a soil-type dependence of BX feedbacks. We demonstrated that wheat also responded to BX-feedbacks. While the negative growth response to BX-conditioning was conserved in both cereals, insect resistance showed opposite patterns, with an increase in maize and a decrease in wheat. Wheat pathogen resistance was not affected. Finally and consistent with maize, we found the BX-feedbacks to be cultivar-specific. Taken together, BX-feedbacks affected cereal growth and resistance in a soil and genotype-dependent manner. Cultivar-specificity of BX-feedbacks is a key finding, as it hides the potential to optimize crops that avoid negative plant-soil feedbacks in rotations.

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.

Tolerance to oxidative stress is associated with both oxidative stress response and inherent growth in a fungal wheat pathogen

Reactive oxygen species are toxic byproducts of aerobic respiration that are also important in mediating a diversity of cellular functions. Reactive oxygen species form an important component of plant defenses to inhibit microbial pathogens during pathogen-plant interactions. Tolerance to oxidative stress is likely to make a significant contribution to the viability and pathogenicity of plant pathogens, but the complex network of oxidative stress responses hinders identification of the genes contributing to this trait. Here, we employed a forward genetic approach to investigate the genetic architecture of oxidative stress tolerance in the fungal wheat pathogen Zymoseptoria tritici. We used quantitative trait locus (QTL) mapping of growth and melanization under axenic conditions in two cross-populations to identify genomic regions associated with tolerance to oxidative stress. We found that QTLs associated with growth under oxidative stress as well as inherent growth can affect oxidative stress tolerance, and we identified two uncharacterized genes in a major QTL associated with this trait. Our data suggest that melanization does not affect tolerance to oxidative stress, which differs from what was found for animal pathogens. This study provides a whole-genome perspective on the genetic basis of oxidative stress tolerance in a plant pathogen.

Characterization of Effector-Target Interactions in Necrotrophic Pathosystems Reveals Trends and Variation in Host Manipulation

Great strides have been made in defining the details of the plant defense response involving biotrophic fungal and bacterial pathogens. The groundwork for the current model was laid by H.H. Flor and others who defined the gene-for-gene hypothesis, which is now known to involve effector-triggered immunity (ETI). PAMP-triggered immunity (PTI) is also a highly effective response to most pathogens because of the recognition of common pathogen molecules by pattern recognition receptors. In this article, we consider the three pathogens that make up the foliar disease complex of wheat, Zymoseptoria tritici, Pyrenophora tritici-repentis, and Parastagonospora nodorum, to review the means by which necrotrophic pathogens circumvent, or outright hijack, the ETI and PTI pathways to cause disease.

A gene-for-gene interaction involving a 'late' effector contributes to quantitative resistance to the stem canker disease in Brassica napus

The control of stem canker disease of Brassica napus (rapeseed), caused by the fungus Leptosphaeria maculans is based largely on plant genetic resistance: single-gene specific resistance (Rlm genes) or quantitative, polygenic, adult-stage resistance. Our working hypothesis was that quantitative resistance partly obeys the gene-for-gene model, with resistance genes 'recognizing' fungal effectors expressed during late systemic colonization. Five LmSTEE (stem-expressed effector) genes were selected and placed under the control of the AvrLm4-7 promoter, an effector gene highly expressed at the cotyledon stage of infection, for miniaturized cotyledon inoculation test screening of a gene pool of 204 rapeseed genotypes. We identified a rapeseed genotype, 'Yudal', expressing hypersensitive response to LmSTEE98. The LmSTEE98-RlmSTEE98 interaction was further validated by inactivation of the LmSTEE98 gene with a CRISPR-Cas9 approach. Isolates with mutated versions of LmSTEE98 induced more severe stem symptoms than the wild-type isolate in 'Yudal'. This single-gene resistance was mapped in a 0.6 cM interval of the 'Darmor_bzh' × 'Yudal' genetic map. One typical gene-for-gene interaction contributes partly to quantitative resistance when L. maculans colonizes the stems of rapeseed. With numerous other effectors specific to stem colonization, our study provides a new route for resistance gene discovery, elucidation of quantitative resistance mechanisms and selection for durable resistance.

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.

Mapping for Adult-Plant Resistance Against Septoria Tritici Blotch in a Common Wheat Line Murga

Septoria tritici blotch (STB) is a major foliar disease globally that is notorious for quickly developing fungicide resistance, making host resistance an indispensable component in mitigating STB. The International Maize and Wheat Improvement Center (CIMMYT) wheat line Murga is well known for its high, durable, and broad-spectrum resistance against STB infection. This study aimed to investigate the resistance mechanism of Murga to facilitate its utilization in breeding. A recombinant inbred line population was derived from a cross between Murga and STB-susceptible line Huirivis#1, comprising 297 progenies. The population was evaluated for adult-plant STB resistance in Toluca, Mexico (from 2017 to 2019), and in La Estanzuela, Uruguay (from 2016 to 2018). Genotyping was performed with the DArTseq platform. Quantitative trait locus (QTL) mapping indicated a major and stable QTL on chromosome 3DL, explaining a phenotypic variation for STB of 41.2 to 62.5% in Mexico and 27.5 to 40.3% in Uruguay. This QTL was regarded as Stb16 based on the comparison of its physical position, the possible origin from synthetic wheat, and its broad-spectrum resistance. Additional QTLs with minor effects were identified on chromosomes 2B, 2D, 3A, 3B, and 5B. The QTL on 5BS was significant in four of the six environments and must be new. Murga was the resistant donor for all QTLs except for those on 2B and 3A. Being an elite breeding line, Stb16 carrier Murga could be used as a promising STB resistance donor. Rational employment of Stb16 could contribute to STB management yet avoid the rapid emergence of Stb16-virulent isolates.

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.

QTL mapping in Fusarium graminearum identified an allele of FgVe1 involved in reduced aggressiveness

Fusarium graminearum is one of the most frequent causal agents of the Fusarium Head Blight, a cereal disease spread throughout the world, reducing grain production and quality. F. graminearum isolates are genetically and phenotypically highly diverse. Notably, remarkable variations of aggressiveness between isolates have been observed, which could reflect an adaptive potential of this pathogen. In this study, we aimed to characterize the genetic basis of aggressiveness variation observed in an F1 population (n = 94), for which genome sequences of both parental strains are available. Aggressiveness was assessed by a panel of in planta and in vitro proxies during two phenotyping trials including, among others, disease severity and mycotoxin accumulation in wheat spike. One major and single QTL was mapped for all the traits measured, on chromosome I, that explained up to 90% of the variance for disease severity. The confidence interval at the QTL spanned 1.2 Mb and contained 428 genes on the reference genome. Of these, four candidates were selected based on the postulate that a non-synonymous mutation affecting protein function may be responsible for phenotypic differences. Finally, a new mutation was identified and functionally validated in the gene FgVe1, coding for a velvet protein known to be involved in pathogenicity and secondary metabolism production in several fungi.

Mapping the adaptive landscape of a major agricultural pathogen reveals evolutionary constraints across heterogeneous environments

The adaptive potential of pathogens in novel or heterogeneous environments underpins the risk of disease epidemics. Antagonistic pleiotropy or differential resource allocation among life-history traits can constrain pathogen adaptation. However, we lack understanding of how the genetic architecture of individual traits can generate trade-offs. Here, we report a large-scale study based on 145 global strains of the fungal wheat pathogen Zymoseptoria tritici from four continents. We measured 50 life-history traits, including virulence and reproduction on 12 different wheat hosts and growth responses to several abiotic stressors. To elucidate the genetic basis of adaptation, we used genome-wide association mapping coupled with genetic correlation analyses. We show that most traits are governed by polygenic architectures and are highly heritable suggesting that adaptation proceeds mainly through allele frequency shifts at many loci. We identified negative genetic correlations among traits related to host colonization and survival in stressful environments. Such genetic constraints indicate that pleiotropic effects could limit the pathogen's ability to cause host damage. In contrast, adaptation to abiotic stress factors was likely facilitated by synergistic pleiotropy. Our study illustrates how comprehensive mapping of life-history trait architectures across diverse environments allows to predict evolutionary trajectories of pathogens confronted with environmental perturbations.

Torres D, Usami T, Schouten HJ, Bai Y, Seidl MF, Thomma BPHJ. Comparative genomics reveals the in planta-secreted Verticillium dahliae Av2 effector protein recognized in tomato plants that carry the V2 resistance locus

Plant pathogens secrete effector molecules during host invasion to promote colonization. However, some of these effectors become recognized by host receptors to mount a defence response and establish immunity. Recently, a novel resistance was identified in wild tomato, mediated by the single dominant V2 locus, to control strains of the soil-borne vascular wilt fungus Verticillium dahliae that belong to race 2. With comparative genomics of race 2 strains and resistance-breaking race 3 strains, we identified the avirulence effector that activates V2 resistance, termed Av2. We identified 277 kb of race 2-specific sequence comprising only two genes encoding predicted secreted proteins that are expressed during tomato colonization. Subsequent functional analysis based on genetic complementation into race 3 isolates and targeted deletion from the race 1 isolate JR2 and race 2 isolate TO22 confirmed that one of the two candidates encodes the avirulence effector Av2 that is recognized in V2 tomato plants. Two Av2 allelic variants were identified that encode Av2 variants that differ by a single acid. Thus far, a role in virulence could not be demonstrated for either of the two variants.

Mixed infections alter transmission potential in a fungal plant pathogen

Infections by more than one strain of a pathogen predominate under natural conditions. Mixed infections can have significant, though often unpredictable, consequences for overall virulence, pathogen transmission and evolution. However, effects of mixed infection on disease development in plants often remain unclear and the critical factors that determine the outcome of mixed infections remain unknown. The fungus Zymoseptoria tritici forms genetically diverse infections in wheat fields. Here, for a range of pathogen traits, we experimentally decompose the infection process to determine how the outcomes and consequences of mixed infections are mechanistically realized. Different strains of Z. tritici grow in close proximity and compete in the wheat apoplast, resulting in reductions in growth of individual strains and in pathogen reproduction. We observed different outcomes of competition at different stages of the infection. Overall, more virulent strains had higher competitive ability during host colonization, and less virulent strains had higher transmission potential. We showed that within‐host competition can have a major effect on infection dynamics and pathogen population structure in a pathogen and host genotype‐specific manner. Consequently, mixed infections likely have a major effect on the development of septoria tritici blotch epidemics and the evolution of virulence in Z. tritici.

A small secreted protein from Zymoseptoria tritici interacts with a wheat E3 ubiquitin ligase to promote disease

Septoria tritici blotch (STB), caused by the ascomycete fungus Zymoseptoria tritici, is a major threat to wheat production worldwide. The Z. tritici genome encodes many small secreted proteins (ZtSSPs) that are likely to play a key role in the successful colonization of host tissues. However, few of these ZtSSPs have been functionally characterized for their role during infection. In this study, we identified and characterized a small, conserved cysteine-rich secreted effector from Z. tritici which has homologues in other plant pathogens in the Dothideomycetes. ZtSSP2 was expressed throughout Z. tritici infection in wheat, with the highest levels observed early during infection. A yeast two-hybrid assay revealed an interaction between ZtSSP2 and wheat E3 ubiquitin ligase (TaE3UBQ) in yeast, and this was further confirmed in planta using bimolecular fluorescence complementation and co-immunoprecipitation. Down-regulation of this wheat E3 ligase using virus-induced gene silencing increased the susceptibility of wheat to STB. Together, these results suggest that TaE3UBQ is likely to play a role in plant immunity to defend against Z. tritici.

A small secreted protein from Zymoseptoria tritici interacts with a wheat E3 ubiquitin ligase to promote disease

Septoria tritici blotch (STB), caused by the ascomycete fungus Zymoseptoria tritici, is a major threat to wheat production worldwide. The Z. tritici genome encodes many small secreted proteins (ZtSSPs) that are likely to play a key role in the successful colonization of host tissues. However, few of these ZtSSPs have been functionally characterized for their role during infection. In this study, we identified and characterized a small, conserved cysteine-rich secreted effector from Z. tritici which has homologues in other plant pathogens in the Dothideomycetes. ZtSSP2 was expressed throughout Z. tritici infection in wheat, with the highest levels observed early during infection. A yeast two-hybrid assay revealed an interaction between ZtSSP2 and wheat E3 ubiquitin ligase (TaE3UBQ) in yeast, and this was further confirmed in planta using bimolecular fluorescence complementation and co-immunoprecipitation. Down-regulation of this wheat E3 ligase using virus-induced gene silencing increased the susceptibility of wheat to STB. Together, these results suggest that TaE3UBQ is likely to play a role in plant immunity to defend against Z. tritici.

Maintenance of variation in virulence and reproduction in populations of an agricultural plant pathogen

Genetic diversity within pathogen populations is critically important for predicting pathogen evolution, disease outcomes and prevalence. However, we lack a good understanding of the processes maintaining genetic variation and constraints on pathogen life-history evolution. Here, we analysed interactions between 12 wheat host genotypes and 145 strains of Zymoseptoria tritici from five global populations to investigate the evolution and maintenance of variation in pathogen virulence and reproduction. We found a strong positive correlation between virulence (amount of leaf necrosis) and reproduction (pycnidia density within lesions), with substantial variation in both traits maintained within populations. On average, highly virulent isolates exhibited higher reproduction, which might increase transmission potential in agricultural fields planted to homogeneous hosts at a high density. We further showed that pathogen strains with a narrow host range (i.e. specialists) for reproduction were on average less virulent, and those with a broader host range (i.e. generalists) were on average less fecund on a given specific host. These costs associated with adaptation to different host genotypes might constrain the emergence of generalists by disrupting the directional evolution of virulence and fecundity. We conclude that selection favouring pathogen strains that are virulent across diverse hosts, coupled with selection that maximizes fecundity on specific hosts, may explain the maintenance of these pathogenicity traits within and among populations.

A wheat cysteine-rich receptor-like kinase confers broad-spectrum resistance against Septoria tritici blotch

The poverty of disease resistance gene reservoirs limits the breeding of crops for durable resistance against evolutionary dynamic pathogens. Zymoseptoria tritici which causes Septoria tritici blotch (STB), represents one of the most genetically diverse and devastating wheat pathogens worldwide. No fully virulent Z. tritici isolates against synthetic wheats carrying the major resistant gene Stb16q have been identified. Here, we use comparative genomics, mutagenesis and complementation to identify Stb16q, which confers broad-spectrum resistance against Z. tritici. The Stb16q gene encodes a plasma membrane cysteine-rich receptor-like kinase that was recently introduced into cultivated wheat and which considerably slows penetration and intercellular growth of the pathogen.

The complex genomic basis of rapid convergent adaptation to pesticides across continents in a fungal plant pathogen

Convergent evolution leads to identical phenotypic traits in different species or populations. Convergence can be driven by standing variation allowing selection to favour identical alleles in parallel or the same mutations can arise independently. However, the molecular basis of such convergent adaptation remains often poorly resolved. Pesticide resistance in agricultural ecosystems is a hallmark of convergence in phenotypic traits. Here, we analyse the major fungal pathogen Zymoseptoria tritici causing serious losses on wheat and with fungicide resistance emergence across several continents. We sampled three population pairs each from a different continent spanning periods early and late in the application of fungicides. To identify causal loci for resistance, we combined knowledge from molecular genetics work and performed genome-wide association studies (GWAS) on a global set of isolates. We discovered yet unknown factors in azole resistance including a gene encoding membrane associated functions. We found strong support for the "hotspot" model of resistance evolution with convergent changes in a small set of loci but additional loci showed more population-specific allele frequency changes. Genome-wide scans of selection showed that half of all known resistance loci were overlapping a selective sweep region. Hence, the application of fungicides was one of the major selective agents acting on the pathogen over the past decades. Furthermore, loci identified through GWAS showed the highest overlap with selective sweep regions underlining the importance to map phenotypic trait variation in evolving populations. Our population genomic analyses highlighted that both de novo mutations and gene flow contributed to convergent pesticide adaptation.

Chromatin Dynamics Contribute to the Spatiotemporal Expression Pattern of Virulence Genes in a Fungal Plant Pathogen

Fungal plant pathogens possess a large repertoire of genes encoding putative effectors, which are crucial for infection. Many of these genes are expressed at low levels in the absence of the host but are strongly induced at specific stages of the infection. The mechanisms underlying this transcriptional reprogramming remain largely unknown. We investigated the role of the genomic environment and associated chromatin modifications of effector genes in controlling their expression pattern in the fungal wheat pathogen Zymoseptoria tritici. Depending on their genomic location, effector genes are epigenetically repressed in the absence of the host and during the initial stages of infection. Derepression of effector genes occurs mainly during and after penetration of plant leaves and is associated with changes in histone modifications. Our work demonstrates the role of chromatin in shaping the expression of virulence components and, thereby, the interaction between fungal pathogens and their plant hosts.

Dynamic changes in transcription profiles are key for the success of pathogens in colonizing their hosts. In many pathogens, genes associated with virulence, such as effector genes, are located in regions of the genome that are rich in transposable elements and heterochromatin. The contribution of chromatin modifications to gene expression in pathogens remains largely unknown. Using a combination of a reporter gene-based approach and chromatin immunoprecipitation, we show that the heterochromatic environment of effector genes in the fungal plant pathogen Zymoseptoria tritici is a key regulator of their specific spatiotemporal expression patterns. Enrichment in trimethylated lysine 27 of histone H3 dictates the repression of effector genes in the absence of the host. Chromatin decondensation during host colonization, featuring a reduction in this repressive modification, indicates a major role for epigenetics in effector gene induction. Our results illustrate that chromatin modifications triggered during host colonization determine the specific expression profile of effector genes at the cellular level and, hence, provide new insights into the regulation of virulence in fungal plant pathogens.

Extracellular vesicles from the apoplastic fungal wheat pathogen Zymoseptoria tritici

BACKGROUND

The fungal pathogen Zymoseptoria tritici is a significant constraint to wheat production in temperate cropping regions around the world. Despite its agronomic impacts, the mechanisms allowing the pathogen to asymptomatically invade and grow in the apoplast of wheat leaves before causing extensive host cell death remain elusive. Given recent evidence of extracellular vesicles (EVs)-secreted, membrane-bound nanoparticles containing molecular cargo-being implicated in extracellular communication between plants and fungal pathogen, we have initiated an in vitro investigation of EVs from this apoplastic fungal wheat pathogen. We aimed to isolate EVs from Z. tritici broth cultures and examine their protein composition in relation to the soluble protein in the culture filtrate and to existing fungal EV proteomes.

RESULTS

Zymoseptoria tritici EVs were isolated from broth culture filtrates using differential ultracentrifugation (DUC) and examined with transmission electron microscopy (TEM) and nanoparticle tracking analysis (NTA). Z. tritici EVs were observed as a heterogeneous population of particles, with most between 50 and 250 nm. These particles were found in abundance in the culture filtrates of viable Z. tritici cultures, but not heat-killed cultures incubated for an equivalent time and of comparable biomass. Bottom-up proteomic analysis using LC-MS/MS, followed by stringent filtering revealed 240 Z. tritici EV proteins. These proteins were distinct from soluble proteins identified in Z. tritici culture filtrates, but were similar to proteins identified in EVs from other fungi, based on sequence similarity analyses. Notably, a putative marker protein recently identified in Candida albicans EVs was also consistently detected in Z. tritici EVs.

CONCLUSION

We have shown EVs can be isolated from the devastating fungal wheat pathogen Z. tritici and are similar to protein composition to previously characterised fungal EVs. EVs from human pathogenic fungi are implicated in virulence, but the role of EVs in the interaction of phytopathogenic fungi and their hosts is unknown. These in vitro analyses provide a basis for expanding investigations of Z. tritici EVs in planta, to examine their involvement in the infection process of this apoplastic wheat pathogen and more broadly, advance understanding of noncanonical secretion in filamentous plant pathogens.

Evolution and Plasticity of the Transcriptome Under Temperature Fluctuations in the Fungal Plant Pathogen Zymoseptoria tritici

Most species live in a variable environment in nature. Yet understanding the evolutionary processes underlying molecular adaptation to fluctuations remains a challenge. In this study we investigate the transcriptome of the fungal wheat pathogen Zymoseptoria tritici after experimental evolution under stable or fluctuating temperature, by comparing ancestral and evolved populations simultaneously. We found that temperature regimes could have a large and pervasive effect on the transcriptome evolution, with as much as 38% of the genes being differentially expressed between selection regimes. Although evolved lineages showed different changes of gene expression based on ancestral genotypes, we identified a set of genes responding specifically to fluctuation. We found that transcriptome evolution in fluctuating conditions was repeatable between parallel lineages initiated from the same genotype for about 60% of the differentially expressed genes. Further, we detected several hotspots of significantly differentially expressed genes in the genome, in regions known to be enriched in repetitive elements, including accessory chromosomes. Our findings also evidenced gene expression evolution toward a gain of robustness (loss of phenotypic plasticity) associated with the fluctuating regime, suggesting robustness is adaptive in changing environment. This work provides valuable insight into the role of transcriptional rewiring for rapid adaptation to abiotic changes in filamentous plant pathogens.

The role of vegetative cell fusions in the development and asexual reproduction of the wheat fungal pathogen Zymoseptoria tritici

Background

The ability of fungal cells to undergo cell-to-cell communication and anastomosis, the process of vegetative hyphal fusion, allows them to maximize their overall fitness. Previous studies in a number of fungal species have identified the requirement of several signaling pathways for anastomosis, including the so far best characterized soft (So) gene, and the MAPK pathway components MAK-1 and MAK-2 of Neurospora crassa. Despite the observations of hyphal fusions’ involvement in pathogenicity and host adhesion, the connection between cell fusion and fungal lifestyles is still unclear. Here, we address the role of anastomosis in fungal development and asexual reproduction in Zymoseptoria tritici, the most important fungal pathogen of wheat in Europe.

Results

We show that Z. tritici undergoes self-fusion between distinct cellular structures, and its mechanism is dependent on the initial cell density. Contrary to other fungi, cell fusion in Z. tritici only resulted in cytoplasmic mixing but not in multinucleated cell formation. The deletion of the So orthologous ZtSof1 disrupted cell-to-cell communication affecting both hyphal and germling fusion. We show that Z. tritici mutants for MAPK-encoding ZtSlt2 (orthologous to MAK-1) and ZtFus3 (orthologous to MAK-2) genes also failed to undergo anastomosis, demonstrating the functional conservation of this signaling mechanism across species. Additionally, the ΔZtSof1 mutant was severely impaired in melanization, suggesting that the So gene function is related to melanization. Finally, we demonstrated that anastomosis is dispensable for pathogenicity, but essential for the pycnidium development, and its absence abolishes the asexual reproduction of Z. tritici.

Conclusions

We demonstrate the role for ZtSof1, ZtSlt2, and ZtFus3 in cell fusions of Z. tritici. Cell fusions are essential for different aspects of the Z. tritici biology, and the ZtSof1 gene is a potential target to control septoria tritici blotch (STB) disease.

Genome comparisons suggest an association between Ceratocystis host adaptations and effector clusters in unique transposable element families

Ceratocystis fimbriata is a host specific fungal pathogen of sweet potato (Ipomoea batatas). The closely related species, C. manginecans, is an important pathogen of trees (e.g. Acacia mangium and Mangifera indica) but has never been isolated from tuber crops. The genetic factors that determine the host range and host specificity of these species have not been determined. The aim of this study was to compare the genomes of C. fimbriata and C. manginecans in order to identify species-specific genetic differences that could be associated with host specificity. This included whole-genome alignments as well as comparisons of gene content and transposable elements (TEs). The genomes of the two species were found to be very similar, sharing similar catalogues of CAZymes, peptidases and lipases. However, the genomes of the two species also varied, harbouring species-specific genes (e.g. small secreted effectors, nutrient processing proteins and stress response proteins). A portion of the TEs identified (17%) had a unique distribution in each species. Transposable elements appeared to have played a prominent role in the divergence of the two species because they were strongly associated with chromosomal translocations and inversions as well as with unique genomic regions containing species-specific genes. Two large effector clusters, with unique TEs in each species, were identified. These effectors displayed non-synonymous mutations and deletions, conserved within a species, and could serve as mutational hot-spots for the development of host specificity in the two species.

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.