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

Transposon-Associated Small RNAs Involved in Plant Defense in Poplar

Utilizing high-throughput Illumina sequencing, we examined how small RNA (sRNA) profiles vary in Chinese white poplar (Populus tomentosa) across two pivotal infection stages by the rust fungus Melampsora larici-populina: the biotrophic growth phase (T02; 48 h post infection) and the urediniospore development and dispersal phase (T03; 168 h), both essential for plant colonization and prolonged biotrophic engagement. Far exceeding random expectations, siRNA clusters predominantly arose from transposon regions, with pseudogenes also contributing significantly, and infection-stage-specific variations were notably tied to these transposon-derived siRNAs. As the infection advanced, clusters of 24 nt siRNAs in transposon and intergenic regions exhibited pronounced abundance shifts. An analysis of targets indicated that Populus sRNAs potentially regulate 95% of Melampsora larici-populina genes, with pathogen effector genes showing heightened targeting by sRNAs during the biotrophic and urediniospore phases compared to controls, pointing to selective sRNA-target interactions. In contrast to conserved miRNAs across plant species, Populus-specific miRNAs displayed a markedly greater tendency to target NB-LRR genes. These observations collectively highlight the innovative roles of sRNAs in plant defense, their evolutionary roots, and their dynamic interplay with pathogen coevolution.

A Genome-Wide Association Study Identifies Markers and Candidate Genes Affecting Tolerance to the Wheat Pathogen Zymoseptoria tritici

Plants defend themselves against pathogens using either resistance, measured as the host's ability to limit pathogen multiplication, or tolerance, measured as the host's ability to reduce the negative effects of infection. Tolerance is a promising trait for crop breeding, but its genetic basis has rarely been studied and remains poorly understood. Here, we reveal the genetic basis of leaf tolerance to the fungal pathogen Zymoseptoria tritici that causes the globally important septoria tritici blotch (STB) disease on wheat. Leaf tolerance to Z. tritici is a quantitative trait that was recently discovered in wheat by using automated image analyses that quantified the symptomatic leaf area and counted the number of pycnidia found on the same leaf. A genome-wide association study identified four chromosome intervals associated with tolerance and a separate chromosome interval associated with resistance. Within these intervals, we identified candidate genes, including wall-associated kinases similar to Stb6, the first cloned STB resistance gene. Our analysis revealed a strong negative genetic correlation between tolerance and resistance to STB, indicative of a trade-off. Such a trade-off between tolerance and resistance would hinder breeding simultaneously for both traits, but our findings suggest a way forward using marker-assisted breeding. We expect that the methods described here can be used to characterize tolerance to other fungal diseases that produce visible fruiting bodies, such as speckled leaf blotch on barley, potentially unveiling conserved tolerance mechanisms shared among plant species. [Formula: see text] Copyright © 2025 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

ZymoSoups: A High-Throughput Forward Genetics Method for Rapid Identification of Virulence Genes in Zymoseptoria tritici

Septoria tritici blotch is caused by the fungus Zymoseptoria tritici and poses a major threat to wheat productivity. There are over 20 mapped loci in wheat that confer strong (gene-for-gene) resistance against this pathogen, but the corresponding genes in Z. tritici that confer virulence against distinct R genes remain largely unknown. In this study, we developed a rapid forward genetics methodology to identify genes that enable Z. tritici to gain virulence on previously resistant wheat varieties. We used the known gene-for-gene interaction between Stb6 and AvrStb6 as a proof of concept that this method could quickly recover single candidate virulence genes. We subjected the avirulent Z. tritici strain IPO323, which carries the recognized AvrStb6 allele, to ultraviolet (UV) mutagenesis and generated a library of over 66,000 surviving spores. We screened these survivors on leaves of the resistant wheat variety Cadenza in mixtures (soups) ranging from 100 to 500 survivors per soup. We identified five soups with a gain-of-virulence (GoV) phenotype relative to the IPO323 parental strain and re-sequenced 18 individual isolates, including four control isolates and two isolates lacking virulence, when screened individually. Of the 12 confirmed GoV isolates, one had a single nucleotide polymorphism (SNP) in the AvrStb6 coding region. The other 11 GoV isolates exhibited large (approximately 70 kb) deletions at the end of chromosome 5, including the AvrStb6 locus. Our findings demonstrate the efficiency of this forward genetic approach in elucidating the genetic basis of qualitative resistance to Z. tritici and the potential to rapidly identify other, currently unknown, Avr genes in this pathogen. [Formula: see text] Copyright © 2025 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

The Molecular Dialogue Between Zymoseptoria tritici and Wheat

Zymoseptoria tritici is a highly damaging pathogen that causes high wheat yield losses in temperate climates. Z. tritici emerged during the domestication of wheat in the Fertile Crescent and has been extensively used as a model system for population genetic and genomic studies. New genetic tools and resources have provided a better understanding of the molecular components involved in the wheat-Z. tritici interaction, which is highlighted by the cloning of three wheat resistance genes and four Z. tritici avirulence genes. Despite the considerable progress made in the last few years, the mechanisms that mediate Z. tritici colonization remain largely unknown. In this review, we summarize the latest advances in understanding the molecular components mediating wheat-Z. tritici interactions, and we discuss future research lines to close current knowledge gaps. [Formula: see text] Copyright © 2025 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Diversification, loss, and virulence gains of the major effector AvrStb6 during continental spread of the wheat pathogen Zymoseptoria tritici

Interactions between plant pathogens and their hosts are highly dynamic and mainly driven by pathogen effectors and plant receptors. Host-pathogen co-evolution can cause rapid diversification or loss of pathogen genes encoding host-exposed proteins. The molecular mechanisms that underpin such sequence dynamics remains poorly investigated at the scale of entire pathogen species. Here, we focus on AvrStb6, a major effector of the global wheat pathogen Zymoseptoria tritici, evolving in response to the cognate receptor Stb6, a resistance widely deployed in wheat. We comprehensively captured effector gene evolution by analyzing a global thousand-genome panel using reference-free sequence analyses. We found that AvrStb6 has diversified into 59 protein isoforms with a strong association to the pathogen spreading to new continents. Across Europe, we found the strongest differentiation of the effector consistent with high rates of Stb6 deployment. The AvrStb6 locus showed also a remarkable diversification in transposable element content with specific expansion patterns across the globe. We detected AvrStb6 gene losses and evidence for transposable element-mediated disruptions. We used virulence datasets of genome-wide association mapping studies to predict virulence changes across the global panel. Genomic predictions suggested marked increases in virulence on Stb6 cultivars concomitant with the spread of the pathogen to Europe and the subsequent spread to further continents. Finally, we genotyped French bread wheat cultivars for Stb6 and monitored resistant cultivar deployment concomitant with AvrStb6 evolution. Taken together, our data provides a comprehensive view of how a rapidly diversifying effector locus can undergo large-scale sequence changes concomitant with gains in virulence on resistant cultivars. The analyses highlight also the need for large-scale pathogen sequencing panels to assess the durability of resistance genes and improve the sustainability of deployment strategies.

Comparative Analyses of Compatible and Incompatible Host-Pathogen Interactions Provide Insight into Divergent Host Specialization of Closely Related Pathogens

Host-pathogen co-evolutionary dynamics drive constant changes in plant pathogens to thrive in their plant host. Factors that determine host specificity are diverse and range from molecular and morphological strategies to metabolic and reproductive adaptations. We applied an experimental approach and conducted comparative microscopy, transcriptome analyses, and functional analyses of selected pathogen traits to identify determinants of host specificity in an important wheat pathogen. We included three closely related fungal pathogens, Zymoseptoria tritici, Z. pseudotritici, and Z. ardabiliae, that establish compatible and incompatible interactions with wheat. Although infections of the incompatible species induce plant defenses during invasion of stomatal openings, we found a conserved early-infection program among the three species whereby only 9.2% of the 8,885 orthologous genes are significantly differentially expressed during initial infection. The genes upregulated in Z. tritici likely reflect specialization to wheat, whereas upregulated genes in the incompatible interaction may reflect processes to counteract cellular stress associated with plant defenses. We selected nine candidate genes encoding putative effectors and host-specificity determinants in Z. tritici and deleted these to study their functional relevance. Despite the particular expression patterns of the nine genes, only two mutants were impaired in virulence. We further expressed the Z. tritici proteins in Nicotiana benthamiana to investigate protein function and assess cell death reaction. Hereby, we identify three effectors with cell-death-inducing properties. From the functional analyses, we conclude that the successful infection of Z. tritici in wheat relies on an extensive redundancy of virulence determinants. [Formula: see text] Copyright © 2025 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Genome-wide association study reveals major loci for resistance to septoria tritici blotch in a Tunisian durum wheat collection

Septoria tritici blotch (STB) is a devastating fungal disease affecting durum and bread wheat worldwide. Tunisian durum wheat landraces are reported to be valuable genetic resources for resistance to STB and should prominently be deployed in breeding programs to develop new varieties resistant to STB disease. In this study, a collection of 367 old durum and 6 modern wheat genotypes previously assessed using single Tunisian Zymoseptoria tritici isolate TUN06 during 2016 and 2017 and TM220 isolate during 2017 were phenotyped for resistance to a mixture of isolates (BULK) under field conditions. Significant correlations for disease traits using the three different inoculums were observed. Using 7638 SNP markers, fifty-one marker-trait associations (MTAs) for STB resistance were identified by genome-wide association study (GWAS) at Bonferroni correction threshold of -log10(P) > 5.184 with phenotypic variance explained (PVE) reaching up to 58%. A total of eleven QTL were identified using TUN06 isolate mean disease scoring (TUNMeanD and TUNMeanA) including threeQTL controlling resistance to both isolates TUN06 and TM220. A major QTL was identified on each of chromosomes 1B, 4B, 5A, and 7B, respectively. The QTL on 7B chromosome colocalized with Stb8 identified in bread wheat. Four QTL including the major QTL identified on chromosome 1B were considered as novel. SNP linked to the significant QTL have the potential to be used in marker-assisted selection for breeding for resistance to STB.

The genetic architecture of resistance to septoria tritici blotch in French wheat cultivars

BACKGROUND

Septoria tritici blotch (STB) is one of the most damaging wheat diseases worldwide, and the development of resistant cultivars is of paramount importance for sustainable crop management. However, the genetic basis of the resistance present in elite wheat cultivars remains largely unknown, which limits the implementation of this strategy. A collection of 285 wheat cultivars originating mostly from France was challenged with ten Zymoseptoria tritici isolates at the seedling stage. The collection was further evaluated in seven field trials across France using artificial inoculation.

RESULTS

Genome-wide association study resulted in the detection of 57 wheat QTL, among which 40 were detected at the seedling stage. Three quarters of these QTL were in genomic regions previously reported for to confer resistance to Z. tritici, but 10 QTL are novel and may be of special interest as new sources of resistance. Some QTL colocalise with major Stb resistance genes, suggesting their presence in the French elite winter wheat germplasm. Among them, the three QTL with the strongest effect colocalize with Stb6, Stb9 and Stb18. There was minimal overlap between the QTL detected at the seedling and adult plant stages, with only 1 out of 20 seedling QTL also being detected in field trials inoculated with the same isolate. This suggests that different resistance genes are involved at the seedling and adult plant stages.

CONCLUSION

This work reveals the highly complex genetic architecture of French wheat resistance to STB and provides relatively small QTL intervals, which will be valuable for identifying the underlying causative genes and for marker-assisted selection.

The Zymoseptoria tritici effector Zt-11 contributes to aggressiveness in wheat

Zymoseptoria tritici is an ascomycete fungus and the causal agent of Septoria tritici leaf blotch (STB) in wheat. Z. tritici secretes an array of effector proteins that are likely to facilitate host infection, colonisation and pycnidia production. In this study we demonstrate a role for Zt-11 as a Z. tritici effector during disease progression. Zt-11 is upregulated during the transition of the pathogen from the biotrophic to necrotrophic phase of wheat infection. Deletion of Zt-11 delayed disease development in wheat, reducing the number and size of pycnidia, as well as the number of macropycnidiospores produced by Z. tritici. This delayed disease development by the ΔZt-11 mutants was accompanied by a lower induction of PR genes in wheat, when compared to infection with wildtype Z. tritici. Overall, these data suggest that Zt-11 plays a role in Z. tritici aggressiveness and STB disease progression possibly via a salicylic acid associated pathway.

An array of Zymoseptoria tritici effectors suppress plant immune responses

Zymoseptoria tritici is the most economically significant fungal pathogen of wheat in Europe. However, despite the importance of this pathogen, the molecular interactions between pathogen and host during infection are not well understood. Herein, we describe the use of two libraries of cloned Z. tritici effectors that were screened to identify effector candidates with putative pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI)-suppressing activity. The effectors from each library were transiently expressed in Nicotiana benthamiana, and expressing leaves were treated with bacterial or fungal PAMPs to assess the effectors' ability to suppress reactive oxygen species (ROS) production. From these screens, numerous effectors were identified with PTI-suppressing activity. In addition, some effectors were able to suppress cell death responses induced by other Z. tritici secreted proteins. We used structural prediction tools to predict the putative structures of all of the Z. tritici effectors and used these predictions to examine whether there was enrichment of specific structural signatures among the PTI-suppressing effectors. From among the libraries, multiple members of the killer protein-like 4 (KP4) and killer protein-like 6 (KP6) effector families were identified as PTI suppressors. This observation is intriguing, as these protein families were previously associated with antimicrobial activity rather than virulence or host manipulation. This data provides mechanistic insight into immune suppression by Z. tritici during infection and suggests that, similar to biotrophic pathogens, this fungus relies on a battery of secreted effectors to suppress host immunity during early phases of colonization.

Stomatal penetration: the cornerstone of plant resistance to the fungal pathogen Zymoseptoria tritici

BACKGROUND

Septoria tritici blotch (STB), caused by the foliar fungus Zymoseptoria tritici, is one of the most damaging disease of wheat in Europe. Genetic resistance against this fungus relies on different types of resistance from non-host resistance (NHR) and host species specific resistance (HSSR) to host resistance mediated by quantitative trait loci (QTLs) or major resistance genes (Stb). Characterizing the diversity of theses resistances is of great importance for breeding wheat cultivars with efficient and durable resistance. While the functional mechanisms underlying these resistance types are not well understood, increasing piece of evidence suggest that fungus stomatal penetration and early establishment in the apoplast are both crucial for the outcome of some interactions between Z. tritici and plants. To validate and extend these previous observations, we conducted quantitative comparative phenotypical and cytological analyses of the infection process corresponding to 22 different interactions between plant species and Z. tritici isolates. These interactions included four major bread wheat Stb genes, four bread wheat accessions with contrasting quantitative resistance, two species resistant to Z. tritici isolates from bread wheat (HSSR) and four plant species resistant to all Z. tritici isolates (NHR).

RESULTS

Infiltration of Z. tritici spores into plant leaves allowed the partial bypass of all bread wheat resistances and durum wheat resistance, but not resistances from other plants species. Quantitative comparative cytological analysis showed that in the non-grass plant Nicotiana benthamiana, Z. tritici was stopped before stomatal penetration. By contrast, in all resistant grass plants, Z. tritici was stopped, at least partly, during stomatal penetration. The intensity of this early plant control process varied depending on resistance types, quantitative resistances being the least effective. These analyses also demonstrated that Stb-mediated resistances, HSSR and NHR, but not quantitative resistances, relied on the strong growth inhibition of the few Z. tritici penetrating hyphae at their entry point in the sub-stomatal cavity.

CONCLUSIONS

In addition to furnishing a robust quantitative cytological assessment system, our study uncovered three stopping patterns of Z. tritici by plant resistances. Stomatal resistance was found important for most resistances to Z. tritici, independently of its type (Stb, HSSR, NHR). These results provided a basis for the functional analysis of wheat resistance to Z. tritici and its improvement.

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.