Epidemiological trade-off between intra- and interannual scales in the evolution of aggressiveness in a local plant pathogen population

Modeling plant diseases under climate change: evolutionary perspectives

Infectious plant diseases are a major threat to global agricultural productivity, economic development, and ecological integrity. There is widespread concern that these social and natural disasters caused by infectious plant diseases may escalate with climate change and computer modeling offers a unique opportunity to address this concern. Here, we analyze the intrinsic problems associated with current modeling strategies and highlight the need to integrate evolutionary principles into polytrophic, eco-evolutionary frameworks to improve predictions. We particularly discuss how evolutionary shifts in functional trade-offs, relative adaptability between plants and pathogens, ecosystems, and climate preferences induced by climate change may feedback to future plant disease epidemics and how technological advances can facilitate the generation and integration of this relevant knowledge for better modeling predictions.

Multiple scenarios for sexual crosses in the fungal pathogen Zymoseptoria tritici on wheat residues: Potential consequences for virulence gene transmission

Little is known about the impact of host immunity on sexual reproduction in fungal pathogens. In particular, it is unclear whether crossing requires both sexual partners to infect living plant tissues. We addressed this issue in a three-year experiment investigating different scenarios of Zymoseptoria tritici crosses according to the virulence ('vir') or avirulence ('avr') of the parents against a qualitative resistance gene. Co-inoculations ('vir × vir', 'avr × vir', 'avr × avr') and single inoculations were performed on a wheat cultivar carrying the Stb16q resistance gene (Cellule) and a susceptible cultivar (Apache), in the greenhouse. We assessed the intensity of asexual reproduction by scoring disease severity, and the intensity of sexual reproduction by counting the ascospores discharged from wheat residues. As expected, disease severity was more intense on Cellule for 'vir × vir' co-inoculations than for 'avr × vir' co-inoculations, with no disease for 'avr × avr'. However, all types of co-inoculation yielded sexual offspring, whether or not the parental strains caused plant symptoms. Parenthood was confirmed by genotyping (SSR markers), and the occurrence of crosses between (co-)inoculated and exogenous strains (other strains from the experiment, or from far away) was determined. We showed that symptomatic asexual infection was not required for a strain to participate in sexual reproduction, and deduced from this result that avirulent strains could be maintained asymptomatically "on" or "in" leaf tissues of plants carrying the corresponding resistant gene for long enough to reproduce sexually. In two of the three years, the intensity of sexual reproduction did not differ between the three types of co-inoculation in Cellule, suggesting that crosses involving avirulent strains are not anecdotal. We discuss the possible mechanisms explaining the maintenance of avirulence in Z. tritici populations and the potential impact of particular resistance deployments such as cultivar mixtures for limiting resistance breakdown.

Genome-wide association mapping reveals genes underlying population-level metabolome diversity in a fungal crop pathogen

Background

Fungi produce a wide range of specialized metabolites (SMs) involved in biotic interactions. Pathways for the production of SMs are often encoded in clusters of tightly arranged genes identified as biosynthetic gene clusters. Such gene clusters can undergo horizontal gene transfers between species and rapid evolutionary change within species. The acquisition, rearrangement, and deletion of gene clusters can generate significant metabolome diversity. However, the genetic basis underlying variation in SM production remains poorly understood.

Results

Here, we analyzed the metabolite production of a large population of the fungal pathogen of wheat, Zymoseptoria tritici. The pathogen causes major yield losses and shows variation in gene clusters. We performed untargeted ultra-high performance liquid chromatography-high resolution mass spectrometry to profile the metabolite diversity among 102 isolates of the same species. We found substantial variation in the abundance of the detected metabolites among isolates. Integrating whole-genome sequencing data, we performed metabolite genome-wide association mapping to identify loci underlying variation in metabolite production (i.e., metabolite-GWAS). We found that significantly associated SNPs reside mostly in coding and gene regulatory regions. Associated genes encode mainly transport and catalytic activities. The metabolite-GWAS identified also a polymorphism in the 3′UTR region of a virulence gene related to metabolite production and showing expression variation.

Conclusions

Taken together, our study provides a significant resource to unravel polymorphism underlying metabolome diversity within a species. Integrating metabolome screens should be feasible for a range of different plant pathogens and help prioritize molecular studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01422-z.

How large and diverse are field populations of fungal plant pathogens? The case of Zymoseptoria tritici

Pathogen populations differ in the amount of genetic diversity they contain. Populations carrying higher genetic diversity are thought to have a greater evolutionary potential than populations carrying less diversity. We used published studies to estimate the range of values associated with two critical components of genetic diversity, the number of unique pathogen genotypes and the number of spores produced during an epidemic, for the septoria tritici blotch pathogen Zymoseptoria tritici. We found that wheat fields experiencing typical levels of infection are likely to carry between 3.1 and 14.0 million pathogen genotypes per hectare and produce at least 2.1–9.9 trillion pycnidiospores per hectare. Given the experimentally derived mutation rate of 3 × 10⁻¹⁰ substitutions per site per cell division, we estimate that between 27 and 126 million pathogen spores carrying adaptive mutations to counteract fungicides and resistant cultivars will be produced per hectare during a growing season. This suggests that most of the adaptive mutations that have been observed in Z. tritici populations can emerge through local selection from standing genetic variation that already exists within each field. The consequences of these findings for disease management strategies are discussed.

Patterns of thermal adaptation in a globally distributed plant pathogen: Local diversity and plasticity reveal two-tier dynamics

Plant pathogen populations inhabit patchy environments with contrasting, variable thermal conditions. We investigated the diversity of thermal responses in populations sampled over contrasting spatiotemporal scales, to improve our understanding of their dynamics of adaptation to local conditions. Samples of natural populations of the wheat pathogen Zymoseptoria tritici were collected from sites within the Euro-Mediterranean region subject to a broad range of climatic conditions. We tested for local adaptation, by accounting for the diversity of responses at the individual and population levels on the basis of key thermal performance curve parameters and "thermotype" (groups of individuals with similar thermal responses) composition. The characterization of phenotypic responses and genotypic structure revealed the following: (i) a high degree of individual plasticity and variation in sensitivity to temperature conditions across spatiotemporal scales and populations; and (ii) geographic variation in thermal response among populations, with major alterations due to seasonal patterns over the wheat-growing season. The seasonal shifts in functional composition suggest that populations are locally structured by selection, contributing to adaptation patterns. Further studies combining selection experiments and modeling are required to determine how functional group selection drives population dynamics and adaptive potential in response to thermal heterogeneity.

Population Structure of Pyrenophora teres f. teres Barley Pathogens from Different Continents

Net form net blotch disease, caused by Pyrenophora teres f. teres, results in significant yield losses to barley industries. Up-to-date knowledge of the genetic diversity and structure of pathogen populations is critical for elucidating the disease epidemiology and unraveling pathogen survival and dispersal mechanisms. Thus, this study investigated long-distance dispersal and adaptation by analyzing the genetic structure of 250 P. teres f. teres isolates collected from Australia, Canada, Hungary, and Republic of South Africa (RSA), and historical isolates from Canada, Denmark, Japan, and Sweden. The population genetic structure detected by discriminant analysis of principal components, with the use of 5,890 Diversity Arrays Technology markers, revealed the presence of four clusters. Two of these contained isolates from all regions, and all isolates from RSA were grouped in these two. Australia and Hungary showed three clusters each. One of the Australian clusters contained only Australian isolates. One of the Hungarian clusters contained only Hungarian isolates and one Danish isolate. STRUCTURE analysis indicated that some isolates from Australia and Hungary shared recent ancestry with RSA, Canada, and historical isolates and were thus admixed. Subdivisions of the neighbor joining network indicated that isolates from distinct countries were closely related, suggesting that multiple introduction events conferred genetic heterogeneity in these countries. Through a neighbor joining analysis and amplification with form-specific DNA markers, we detected two hybrid isolates, CBS 281.31 from Japan and H-919 from Hungary, collected in 1931 and 2018, respectively. These results provide a foundation for exploring improved management of disease incursions and pathogen control through strategic deployment of resistance.

Population-level deep sequencing reveals the interplay of clonal and sexual reproduction in the fungal wheat pathogen Zymoseptoria tritici

Pathogens cause significant challenges to global food security. On annual crops, pathogens must re-infect from environmental sources in every growing season. Fungal pathogens have evolved mixed reproductive strategies to cope with the distinct challenges of colonizing growing plants. However, how pathogen diversity evolves during growing seasons remains largely unknown. Here, we performed a deep hierarchical sampling in a single experimental wheat field infected by the major fungal pathogen Zymoseptoria tritici. We analysed whole genome sequences of 177 isolates collected from 12 distinct cultivars replicated in space at three time points of the growing season to maximize capture of genetic diversity. The field population was highly diverse with 37 SNPs per kilobase, a linkage disequilibrium decay within 200-700 bp and a high effective population size. Using experimental infections, we tested a subset of the collected isolates on the dominant cultivar planted in the field. However, we found no significant difference in virulence of isolates collected from the same cultivar compared to isolates collected on other cultivars. About 20 % of the isolate genotypes were grouped into 15 clonal groups. Pairs of clones were disproportionally found at short distances (<5 m), consistent with experimental estimates for per-generation dispersal distances performed in the same field. This confirms predominant leaf-to-leaf transmission during the growing season. Surprisingly, levels of clonality did not increase over time in the field although reproduction is thought to be exclusively asexual during the growing season. Our study shows that the pathogen establishes vast and stable gene pools in single fields. Monitoring short-term evolutionary changes in crop pathogens will inform more durable strategies to contain diseases.

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.

A fungal pathogen induces systemic susceptibility and systemic shifts in wheat metabolome and microbiome composition

Yield losses caused by fungal pathogens represent a major threat to global food production. One of the most devastating fungal wheat pathogens is Zymoseptoria tritici. Despite the importance of this fungus, the underlying mechanisms of plant–pathogen interactions are poorly understood. Here we present a conceptual framework based on coinfection assays, comparative metabolomics, and microbiome profiling to study the interaction of Z. tritici in susceptible and resistant wheat. We demonstrate that Z. tritici suppresses the production of immune-related metabolites in a susceptible cultivar. Remarkably, this fungus-induced immune suppression spreads within the leaf and even to other leaves, a phenomenon that we term “systemic induced susceptibility”. Using a comparative metabolomics approach, we identify defense-related biosynthetic pathways that are suppressed and induced in susceptible and resistant cultivars, respectively. We show that these fungus-induced changes correlate with changes in the wheat leaf microbiome. Our findings suggest that immune suppression by this hemibiotrophic pathogen impacts specialized plant metabolism, alters its associated microbial communities, and renders wheat vulnerable to further infections.

The fungal plant pathogen Zymoseptoria tritici is a major threat to wheat yield. Here Seybold et al. show that Z. tritici can suppress immune responses not only in infected tissue but also on other leaves, a phenomenon termed “systemic induced susceptibility” that is correlated with systemic changes in metabolite accumulation.

Crop Residues in Wheat-Oilseed Rape Rotation System: a Pivotal, Shifting Platform for Microbial Meetings

Crop residues are a crucial ecological niche with a major biological impact on agricultural ecosystems. In this study, we used a combined diachronic and synchronic field experiment based on wheat-oilseed rape rotations to test the hypothesis that plant is a structuring factor of microbial communities in crop residues, and that this effect decreases over time with their likely progressive degradation and colonisation by other microorganisms. We characterised an entire fungal and bacterial community associated with 150 wheat and oilseed rape residue samples at a plurennial scale by metabarcoding. The impact of plant species on the residue microbiota decreased over time and our data revealed turnover, with the replacement of oligotrophs, often plant-specific genera (such as pathogens) by copiotrophs, belonging to more generalist genera. Within a single cropping season, the plant-specific genera and species were gradually replaced by taxa that are likely to originate from the soil. These changes occurred more rapidly for bacteria than for fungi, known to degrade complex compounds. Overall, our findings suggest that crop residues constitute a key fully-fledged microbial ecosystem. Taking into account this ecosystem, that has been neglected for too long, is essential, not only to improve the quantitative management of residues, the presence of which can be detrimental to crop health, but also to identify groups of beneficial microorganisms. Our findings are of particular importance, because the wheat-oilseed rape rotation, in which no-till practices are frequent, is particularly widespread in the European arable cropping systems.

Sexual Reproduction in the Fungal Foliar Pathogen Zymoseptoria tritici Is Driven by Antagonistic Density Dependence Mechanisms

This study provides empirical evidence for antagonistic density dependence mechanisms driving sexual reproduction in the wheat fungal pathogen Zymoseptoria tritici. Biparental crosses with 12 increasing inoculum concentrations, in controlled conditions, showed that sexual reproduction in Z. tritici was impacted by an Allee effect due to mate limitation and a competition with asexual multiplication for resource allocation. The highest number of ascospores discharged was reached at intermediate inoculum concentrations (from 5 × 10⁴ conidia mL^-1 to 10⁶ conidia mL^-1). Consistent with these results for controlled co-inoculation, we found that the intensity of sexual reproduction varied with both cropping period and the vertical position of the host tissues in the field, with a maximum between 25 and 35 cm above the ground. An optimal lesion density (disease severity of 30 to 45%) maximizing offspring (ascospores) number was established, and its eco-evolutionary consequences are considered here. Two ecological mechanisms may be involved: competition for resources between the two modes of reproduction (decrease in the host resources available for sexual reproduction due to their prior use in asexual multiplication), and competitive disequilibrium between the two parental isolates, due to differential interaction dynamics with the host, for example, leading to an imbalance between mating types. A conceptual model based on these results suggests that sexual reproduction plays a key role in the evolution of pathogenicity traits, including virulence and aggressiveness. Ecological knowledge about the determinants of sexual reproduction in Z. tritici may, therefore, open up new perspectives for the management of other fungal foliar pathogens with dual modes of reproduction.