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.

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.

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.

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.

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 (h2 = 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.

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 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.

Natural selection drives population divergence for local adaptation in a wheat pathogen.

Evolution favors the emergence of locally-adapted optimum phenotypes that are likely to differ across a wide array of environmental conditions. The emergence of favorable adaptive characteristics is accelerated in agricultural pathogens due to the unique properties of agro-ecosystems. We performed a QST - FST comparison using 164 strains of Parastagonospora nodorum sampled from eight global field populations to disentangle the predominant evolutionary forces driving population divergence in a wheat pathogen. We used digital image analysis to obtain quantitative measurements of growth rate and melanization at different temperatures and under different fungicide concentrations in a common garden experiment. FST measures were based on complete genome sequences obtained for all 164 isolates. Our analyses indicated that all measured traits were under selection. Growth rates at 18 °C and 24 °C were under stabilizing selection (QST < FST), while diversifying selection (QST > FST) was the predominant evolutionary force affecting growth under fungicide and high temperature stress. Stabilizing selection (QST < FST) was the predominant force affecting melanization across the different environments. Melanin production increased at 30 °C but was negatively correlated with higher growth rates, consistent with a trade-off under heat stress. Our results demonstrate that global populations of P. nodorum possess significant evolutionary potential to adapt to changing local conditions, including warmer temperatures and applications of fungicides.

A tradeoff between tolerance and resistance to a major fungal pathogen in elite wheat cultivars.

Tolerance and resistance represent two strategies that hosts evolved to protect themselves from pathogens. Tolerance alleviates the reduction in host fitness due to infection without reducing a pathogen's growth, whereas resistance reduces pathogen growth. We investigated the tolerance of wheat to the major fungal pathogen Zymoseptoria tritici in 335 elite wheat cultivars. We used a novel digital phenotyping approach that included 11 152 infected leaves and counted 2069 048 pathogen fruiting bodies. We discovered a new component of tolerance that is based on the relationship between the green area remaining on a leaf and the number of pathogen fruiting bodies. We found a negative correlation between tolerance and resistance among intolerant cultivars, presenting the first compelling evidence for a tradeoff between tolerance and resistance to plant pathogens. Surprisingly, the tradeoff arises due to limits in the host resources available to the pathogen and not due to metabolic constraints, contrary to what ecological theory suggests. The mechanism underlying this tradeoff may be relevant for many plant diseases in which the amount of host resources available to the pathogen can limit the pathogen population. Our analysis indicates that European wheat breeders may have selected for tolerance instead of resistance to an important pathogen.

The emergence of the multi-species NIP1 effector in Rhynchosporium was accompanied by high rates of gene duplications and losses.

Plant pathogens secrete effector proteins to manipulate the host and facilitate infection. Cognate hosts trigger strong defence responses upon detection of these effectors. Consequently, pathogens and hosts undergo rapid coevolutionary arms races driven by adaptive evolution of effectors and receptors. Because of their high rate of turnover, most effectors are thought to be species-specific and the evolutionary trajectories are poorly understood. Here, we investigate the necrosis-inducing protein 1 (NIP1) effector in the multihost pathogen genus Rhynchosporium. We retraced the evolutionary history of the NIP1 locus using whole-genome assemblies of 146 strains covering four closely related species. NIP1 orthologues were present in all species but the locus consistently segregated presence-absence polymorphisms suggesting long-term balancing selection. We also identified previously unknown paralogues of NIP1 that were shared among multiple species and showed substantial copy-number variation within R. commune. The NIP1A paralogue was under significant positive selection suggesting that NIP1A is the dominant effector variant coevolving with host immune receptors. Consistent with this prediction, we found that copy number variation at NIP1A had a stronger effect on virulence than NIP1B. Our analyses unravelled the origins and diversification mechanisms of a pathogen effector family shedding light on how pathogens gain adaptive genetic variation.

Transposable element insertions shape gene regulation and melanin production in a fungal pathogen of wheat.

BACKGROUND: Fungal plant pathogens pose major threats to crop yield and sustainable food production if they are highly adapted to their host and the local environment. Variation in gene expression contributes to phenotypic diversity within fungal species and affects adaptation. However, very few cases of adaptive regulatory changes have been reported in fungi and the underlying mechanisms remain largely unexplored. Fungal pathogen genomes are highly plastic and harbor numerous insertions of transposable elements, which can potentially contribute to gene expression regulation. In this work, we elucidated how transposable elements contribute to variation in melanin accumulation, a quantitative trait in fungi that affects survival under stressful conditions. RESULTS: We demonstrated that differential transcriptional regulation of the gene encoding the transcription factor Zmr1, which controls expression of the genes in the melanin biosynthetic gene cluster, is responsible for variation in melanin accumulation in the fungal plant pathogen Zymoseptoria tritici. We show that differences in melanin levels between two strains of Z. tritici are due to two levels of transcriptional regulation: (1) variation in the promoter sequence of Zmr1 and (2) an insertion of transposable elements upstream of the Zmr1 promoter. Remarkably, independent insertions of transposable elements upstream of Zmr1 occurred in 9% of Z. tritici strains from around the world and negatively regulated Zmr1 expression, contributing to variation in melanin accumulation. CONCLUSIONS: Our studies identified two levels of transcriptional control that regulate the synthesis of melanin. We propose that these regulatory mechanisms evolved to balance the fitness costs associated with melanin production against its positive contribution to survival in stressful environments.

Widespread signatures of selection for secreted peptidases in a fungal plant pathogen.

BACKGROUND: Fungal plant pathogens secrete a large arsenal of hydrolytic enzymes during the course of infection, including peptidases. Secreted peptidases have been extensively studied for their role as effectors. In this study, we combined transcriptomics, comparative genomics and evolutionary analyses to investigate all 39 secreted peptidases in the fungal wheat pathogen Zymoseptoria tritici and its close relatives Z. pseudotritici and Z. ardabiliae. RESULTS: RNA-seq data revealed that a majority of the secreted peptidases displayed differential transcription during the course of Z. tritici infection, indicative of specialization for different stages in the life cycle. Evolutionary analyses detected widespread evidence of adaptive evolution acting on at least 28 of the peptidases. A few peptidases displayed lineage-specific rates of molecular evolution, suggesting altered selection pressure in Z. tritici following host specialization on domesticated wheat. The peptidases belonging to MEROPS families A1 and G1 emerged as a particularly interesting group that may play key roles in host-pathogen co-evolution, host adaptation and pathogenicity. Sister genes in the A1 and G1 families showed accelerated substitution rates after gene duplications. CONCLUSIONS: These results suggest widespread evolution of secreted peptidases leading to novel gene functions, consistent with predicted models of "escape from adaptive conflict" and "neo-functionalization". Our analyses identified candidate genes worthy of functional analyses that may encode effector functions, for example by suppressing plant defenses during the biotrophic phase of infection.

Comparative Transcriptomics Reveals How Wheat Responds to Infection by Zymoseptoria tritici.

The fungus Zymoseptoria tritici causes septoria tritici blotch (STB) on wheat, an important disease globally and the most damaging wheat disease in Europe. Despite the global significance of STB, the molecular basis of wheat defense against Z. tritici is poorly understood. Here, we use a comparative transcriptomic study to investigate how wheat responds to infection by four distinct strains of Z. tritici. We examined the response of wheat across the entire infection cycle, identifying both shared responses to the four strains and strain-specific responses. We found that the early asymptomatic phase is characterized by strong upregulation of genes encoding receptor-like kinases and pathogenesis-related proteins, indicating the onset of a defense response. We also identified genes that were differentially expressed among the four fungal strains, including genes related to defense. Genes involved in senescence were induced during both the asymptomatic phase and at late stages of infection, suggesting manipulation of senescence processes by both the plant and the pathogen. Our findings illustrate the need, when identifying important genes affecting disease resistance in plants, to include multiple pathogen strains.

Evidence for local adaptation and pleiotropic effects associated with melanization in a plant pathogenic fungus.

We combined a common garden experimental design with digital image analysis to determine how melanization responds to temperature and fungicide stress in 126 strains of Rhynchosporium commune sampled from nine global field populations. We found that different temperatures and fungicide stress significantly affected the degree of melanization. The nine field populations showed similar patterns in response to the different temperatures. Significant correlations were found between the degree of melanization and the local environment, including mean annual temperature, latitude, and relative humidity, suggesting that melanization is a locally adaptive trait. We also found that melanization is positively correlated with both virulence and fungicide resistance. These findings suggest that melanization has pleiotropic effects in Rhynchosporium commune.

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

Cultivar-strain specificity in the wheat-Zymoseptoria tritici pathosystem determines the infection outcome and is controlled by resistance genes on the host side, many of which have been identified. On the pathogen side, however, the molecular determinants of specificity remain largely unknown. We used genetic mapping, targeted gene disruption and allele swapping to characterise the recognition of the new avirulence factor Avr3D1. We then combined population genetic and comparative genomic analyses to characterise the evolutionary trajectory of Avr3D1. Avr3D1 is specifically recognised by wheat cultivars harbouring the Stb7 resistance gene, triggering a strong defence response without preventing pathogen infection and reproduction. Avr3D1 resides in a cluster of putative effector genes located in a genome region populated by independent transposable element insertions. The gene was present in all 132 investigated strains and is highly polymorphic, with 30 different protein variants identified. We demonstrated that specific amino acid substitutions in Avr3D1 led to evasion of recognition. These results demonstrate that quantitative resistance and gene-for-gene interactions are not mutually exclusive. Localising avirulence genes in highly plastic genomic regions probably facilitates accelerated evolution that enables escape from recognition by resistance proteins.

Genome-wide evidence for divergent selection between populations of a major agricultural pathogen.

The genetic and environmental homogeneity in agricultural ecosystems is thought to impose strong and uniform selection pressures. However, the impact of this selection on plant pathogen genomes remains largely unknown. We aimed to identify the proportion of the genome and the specific gene functions under positive selection in populations of the fungal wheat pathogen Zymoseptoria tritici. First, we performed genome scans in four field populations that were sampled from different continents and on distinct wheat cultivars to test which genomic regions are under recent selection. Based on extended haplotype homozygosity and composite likelihood ratio tests, we identified 384 and 81 selective sweeps affecting 4% and 0.5% of the 35 Mb core genome, respectively. We found differences both in the number and the position of selective sweeps across the genome between populations. Using a XtX-based outlier detection approach, we identified 51 extremely divergent genomic regions between the allopatric populations, suggesting that divergent selection led to locally adapted pathogen populations. We performed an outlier detection analysis between two sympatric populations infecting two different wheat cultivars to identify evidence for host-driven selection. Selective sweep regions harboured genes that are likely to play a role in successfully establishing host infections. We also identified secondary metabolite gene clusters and an enrichment in genes encoding transporter and protein localization functions. The latter gene functions mediate responses to environmental stress, including interactions with the host. The distinct gene functions under selection indicate that both local host genotypes and abiotic factors contributed to local adaptation.

Ranking Quantitative Resistance to Septoria tritici Blotch in Elite Wheat Cultivars Using Automated Image Analysis.

Quantitative resistance is likely to be more durable than major gene resistance for controlling Septoria tritici blotch (STB) on wheat. Earlier studies hypothesized that resistance affecting the degree of host damage, as measured by the percentage of leaf area covered by STB lesions, is distinct from resistance that affects pathogen reproduction, as measured by the density of pycnidia produced within lesions. We tested this hypothesis using a collection of 335 elite European winter wheat cultivars that was naturally infected by a diverse population of Zymoseptoria tritici in a replicated field experiment. We used automated image analysis of 21,420 scanned wheat leaves to obtain quantitative measures of conditional STB intensity that were precise, objective, and reproducible. These measures allowed us to explicitly separate resistance affecting host damage from resistance affecting pathogen reproduction, enabling us to confirm that these resistance traits are largely independent. The cultivar rankings based on host damage were different from the rankings based on pathogen reproduction, indicating that the two forms of resistance should be considered separately in breeding programs aiming to increase STB resistance. We hypothesize that these different forms of resistance are under separate genetic control, enabling them to be recombined to form new cultivars that are highly resistant to STB. We found a significant correlation between rankings based on automated image analysis and rankings based on traditional visual scoring, suggesting that image analysis can complement conventional measurements of STB resistance, based largely on host damage, while enabling a much more precise measure of pathogen reproduction. We showed that measures of pathogen reproduction early in the growing season were the best predictors of host damage late in the growing season, illustrating the importance of breeding for resistance that reduces pathogen reproduction in order to minimize yield losses caused by STB. These data can already be used by breeding programs to choose wheat cultivars that are broadly resistant to naturally diverse Z. tritici populations according to the different classes of resistance.

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

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

Reversing resistance: different routes and common themes across pathogens.

Resistance spreads rapidly in pathogen or pest populations exposed to biocides, such as fungicides and antibiotics, and in many cases new biocides are in short supply. How can resistance be reversed in order to prolong the effectiveness of available treatments? Some key parameters affecting reversion of resistance are well known, such as the fitness cost of resistance. However, the population biological processes that actually cause resistance to persist or decline remain poorly characterized, and consequently our ability to manage reversion of resistance is limited. Where do susceptible genotypes that replace resistant lineages come from? What is the epidemiological scale of reversion? What information do we need to predict the mechanisms or likelihood of reversion? Here, we define some of the population biological processes that can drive reversion, using examples from a wide range of taxa and biocides. These processes differ primarily in the origin of revertant genotypes, but also in their sensitivity to factors such as coselection and compensatory evolution that can alter the rate of reversion, and the likelihood that resistance will re-emerge upon re-exposure to biocides. We therefore argue that discriminating among different types of reversion allows for better prediction of where resistance is most likely to persist.