Evidence of Selection for Fungicide Resistance in Zymoseptoria tritici Populations on Wheat in Western Oregon

Detection of Maternal and Cytoplasmic Effects on Resistance to Zymoseptoria tritici in Durum Wheat

Septoria tritici blotch (STB) is a major disease problem of wheat worldwide. To optimize the introgression of resistance genes in elite genotypes throughout traditional or molecular breeding programs, a full understanding of the quantitative inheritance of resistance to Zymoseptoria tritici, plant height (PH), and thousand kernel weight (TKW) is needed. In this study, maternal and cytoplasmic effects of resistance to STB were investigated using P1 (susceptible, high-yielding line) and P2 (resistant, low-yielding line) durum wheat lines and their F₁, RF₁, F₂, RF₂, BC₁, RBC₁, BC₂, and RBC₂ progeny, assessed for resistance to STB during three growing seasons. Duncan mean's analysis revealed significant differences between generation means for STB, PH, and TKW. The two parents had an extreme pattern. The F₁ and RF₁ segregated close to their respective parents, suggesting the presence of cytoplasmic and maternal genetic effects for Z. tritici resistance, PH, and TKW. Separate generation mean's analysis confirmed the results of the Duncan test. A three-parameter model was found to be not adequate for all traits in all three growing years; while a digenic epistatic model with cytoplasmic or/and maternal effect was adequate for all cases. Narrow-sense heritability was in the range of 50–60%, 30–69%, and 28–31% for STB, PH, and TKW, respectively. For STB, high heritability and the presence of fixable epistatic effect is encouraging and could lead to creating varieties with the right female parent to exploit cytoplasmic and maternal effects in order to improve resistance to Z. tritici in durum wheat.

Adaptation of Biotrophic Leaf Pathogens to Fertilization-Mediated Changes in Plant Traits: A Comparison of the Optimization Principle to Invasion Fitness

One of the conclusions of evolutionary ecology applied to agroecosystem management is that sustainable disease management strategies must be adaptive to overcome the immense adaptive potential of crop pathogens. In this context, knowledge of how pathogens adapt to changes in cultural practices is necessary. In this article we address the issue of the evolutionary response of biotrophic crop pathogens to changes in fertilization practices. For this purpose, we compare predictions of latent period evolution based on three empirical fitness measures (seasonal spore production, within-season exponential growth rate, and area under disease progress curves [AUDPCs]) with predictions based on the concept of invasion fitness from adaptive dynamics. We use pairwise invisibility plots to identify the evolutionarily stable strategies (ESSs) of the pathogen latent period. We find that the ESS latent period is in between the latent periods that maximize the seasonal spore production and the within-season exponential growth rate of the pathogen. The latent periods that maximize the AUDPC are similar to those of the ESS latent periods. The AUDPC may therefore be a critical variable to determine the issue of between-strain competition and shape pathogen evolution.

Differential dynamics of microbial community networks help identify microorganisms interacting with residue-borne pathogens: the case of Zymoseptoria tritici in wheat

BACKGROUND

Wheat residues are a crucial determinant of the epidemiology of Septoria tritici blotch, as they support the sexual reproduction of the causal agent Zymoseptoria tritici. We aimed to characterize the effect of infection with this fungal pathogen on the microbial communities present on wheat residues and to identify microorganisms interacting with it. We used metabarcoding to characterize the microbiome associated with wheat residues placed outdoors, with and without preliminary Z. tritici inoculation, comparing the first set of residues in contact with the soil and a second set without contact with the soil, on four sampling dates in two consecutive years.

RESULTS

The diversity of the tested conditions, leading to the establishment of different microbial communities according to the origins of the constitutive taxa (plant only, or plant and soil), highlighted the effect of Z. tritici on the wheat residue microbiome. Several microorganisms were affected by Z. tritici infection, even after the disappearance of the pathogen. Linear discriminant analyses and ecological network analyses were combined to describe the communities affected by the infection. The number of fungi and bacteria promoted or inhibited by inoculation with Z. tritici decreased over time and was smaller for residues in contact with the soil. The interactions between the pathogen and other microorganisms appeared to be mostly indirect, despite the strong position of the pathogen as a keystone taxon in networks. Direct interactions with other members of the communities mostly involved fungi, including other wheat pathogens. Our results provide essential information about the alterations to the microbial community in wheat residues induced by the mere presence of a fungal pathogen, and vice versa. Species already described as beneficial or biocontrol agents were found to be affected by pathogen inoculation.

CONCLUSIONS

The strategy developed here can be viewed as a proof-of-concept focusing on crop residues as a particularly rich ecological compartment, with a high diversity of fungal and bacterial taxa originating from both the plant and soil compartments, and for Z. tritici-wheat as a model pathosystem. By revealing putative antagonistic interactions, this study paves the way for improving the biological control of residue-borne diseases.

Effect of Long-Term Fungicide Applications on Virulence and Diversity of Colletotrichum spp. Associated to Olive Anthracnose

In this study, the presence and variability of Colletotrichum spp. was evaluated by comparing fungal isolates obtained from olive trees under long-time phytosanitary treatments with trees without any phytosanitary treatments (treated and untreated, respectively). Olive fruits of trees of the highly susceptible ‘Galega vulgar’ cultivar growing in the Alentejo region were used as samples. From the 210 olive trees sampled (half from treated and half from untreated orchards), 125 (59.5%) presented Colletotrichum spp., with a significant lower number of infected trees in treated (39) when compared to untreated orchards (86). The alignment and analysis of beta-tubulin (tub2), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), actin (ACT), chitin synthase (CHS-1) and histone H3 (HIS-3) gene sequences allowed the identification of all 125 isolates as belonging to the C. acutatum complex. The vast majority of the isolates (124) were identified as C. nymphaeae and one isolate, from an untreated tree, was identified as C. godetiae. Isolates were divided into five different groups: Group A: 39 isolates from treated trees matched in 100% with C. nymphaeae sequences from the database; Group B: 76 isolates from untreated trees matched in 100% with C. nymphaeae sequences from the database; Group C: one isolate from untreated trees presenting a single nucleotidic difference in the HIS-3 sequence; Group D: eight isolates from untreated trees presenting differences in two nucleotides in the tub2 sequences that changed the protein structure, together with differences in two specific nucleotides of the GAPDH sequences; Group E: one isolate, from untreated olive trees, matched 100% with C. godetiae sequences from the database in all genes. Considering the similarities of the sampled areas, our results show that the long-time application of fungicides may have caused a reduction in the number of olive trees infected with Colletotrichum spp. but an increase in the number of fruits positive to Colletotrichum spp. within each tree, which may suggest different degrees of virulence of Colletotrichum isolates from trees growing different management regimes. It is imperative that the fungicides described as causing resistance are applied at appropriate times and intervals, since their efficiency decreases when applied incorrectly and new and more virulent species may arise.

Spatiotemporal dynamics of fungicide resistance in the wheat pathogen Zymoseptoria tritici in France

BACKGROUND

Management of pesticide resistance is a major issue in modern agricultural systems, particularly in the context of the broader challenge of reducing pesticide use. However, such management must be adapted to resistance dynamics, which remains challenging to predict due to its dependence on many biological traits of pests, interactions with the environment and pesticide use. We retrospectively studied the evolution of reported resistances to four modes of action (benzimidazoles, quinone outside inhibitors, sterol demethylation inhibitors and succinate dehydrogenase inhibitors), in French populations of the wheat pathogen Zymoseptoria tritici.

RESULTS

We used statistical models to analyse the Performance trial network data set (2004-2017; ∼ 70 locations in France yearly). They highlighted contrasting behaviours between phenotypes, for example: (i) stable spatial distributions and colonization front structures over time, and (ii) different frequency growth rates at the national scale and between regions.

CONCLUSION

We provide a quantitative description of the spatiotemporal patterns of resistance evolution for fungicides with several modes of action. Moreover, we highlight some unexpected resistance dynamics in France, with major differences between the north and south. This complex pattern of resistance evolution in French populations is consistent with previous descriptions of dynamics at the European scale. These results should make it easier to anticipate evolution locally and to improve the management of resistance. © 2019 Society of Chemical Industry.

A whole genome scan of SNP data suggests a lack of abundant hard selective sweeps in the genome of the broad host range plant pathogenic fungus Sclerotinia sclerotiorum

The pathogenic fungus Sclerotinia sclerotiorum infects over 600 species of plant. It is present in numerous environments throughout the world and causes significant damage to many agricultural crops. Fragmentation and lack of gene flow between populations may lead to population sub-structure. Within discrete recombining populations, positive selection may lead to a ‘selective sweep’. This is characterised by an increase in frequency of a favourable allele leading to reduction in genotypic diversity in a localised genomic region due to the phenomenon of genetic hitchhiking. We aimed to assess whether isolates of S. sclerotiorum from around the world formed genotypic clusters associated with geographical origin and to determine whether signatures of population-specific positive selection could be detected. To do this, we sequenced the genomes of 25 isolates of S. sclerotiorum collected from four different continents–Australia, Africa (north and south), Europe and North America (Canada and the northen United States) and conducted SNP based analyses of population structure and selective sweeps. Among the 25 isolates, there was evidence for two major population clusters. One of these consisted of 11 isolates from Canada, the USA and France (population 1), and the other consisted of nine isolates from Australia and one from Morocco (population 2). The rest of the isolates were genotypic outliers. We found that there was evidence of outcrossing in these two populations based on linkage disequilibrium decay. However, only a single candidate selective sweep was observed, and it was present in population 2. This sweep was close to a Major Facilitator Superfamily transporter gene, and we speculate that this gene may have a role in nutrient uptake from the host. The low abundance of selective sweeps in the S. sclerotiorum genome contrasts the numerous examples in the genomes of other fungal pathogens. This may be a result of its slow rate of evolution and low effective recombination rate due to self-fertilisation and vegetative reproduction.

Precision Phenotyping Reveals Novel Loci for Quantitative Resistance to Septoria Tritici Blotch

Accurate, high-throughput phenotyping for quantitative traits is a limiting factor for progress in plant breeding. We developed an automated image analysis to measure quantitative resistance to septoria tritici blotch (STB), a globally important wheat disease, enabling identification of small chromosome intervals containing plausible candidate genes for STB resistance. 335 winter wheat cultivars were included in a replicated field experiment that experienced natural epidemic development by a highly diverse but fungicide-resistant pathogen population. More than 5.4 million automatically generated phenotypes were associated with 13,648 SNP markers to perform the GWAS. We identified 26 chromosome intervals explaining 1.9-10.6% of the variance associated with four independent resistance traits. Sixteen of the intervals overlapped with known STB resistance intervals, suggesting that our phenotyping approach can identify simultaneously (i.e., in a single experiment) many previously defined STB resistance intervals. Seventeen of the intervals were less than 5 Mbp in size and encoded only 173 genes, including many genes associated with disease resistance. Five intervals contained four or fewer genes, providing high priority targets for functional validation. Ten chromosome intervals were not previously associated with STB resistance, perhaps representing resistance to pathogen strains that had not been tested in earlier experiments. The SNP markers associated with these chromosome intervals can be used to recombine different forms of quantitative STB resistance that are likely to be more durable than pyramids of major resistance genes. Our experiment illustrates how high-throughput automated phenotyping can accelerate breeding for quantitative disease resistance.

Phosphopantetheinyl transferase (Ppt)-mediated biosynthesis of lysine, but not siderophores or DHN melanin, is required for virulence of Zymoseptoria tritici on wheat

Zymoseptoria tritici is the causal agent of Septoria tritici blotch (STB) disease of wheat. Z. tritici is an apoplastic fungal pathogen, which does not penetrate plant cells at any stage of infection, and has a long initial period of symptomless leaf colonisation. During this phase it is unclear to what extent the fungus can access host plant nutrients or communicate with plant cells. Several important primary and secondary metabolite pathways in fungi are regulated by the post-translational activator phosphopantetheinyl transferase (Ppt) which provides an essential co-factor for lysine biosynthesis and the activities of non-ribosomal peptide synthases (NRPS) and polyketide synthases (PKS). To investigate the relative importance of lysine biosynthesis, NRPS-based siderophore production and PKS-based DHN melanin biosynthesis, we generated deletion mutants of ZtPpt. The ∆ZtPpt strains were auxotrophic for lysine and iron, non-melanised and non-pathogenic on wheat. Deletion of the three target genes likely affected by ZtPpt loss of function (Aar- lysine; Nrps1-siderophore and Pks1- melanin), highlighted that lysine auxotrophy was the main contributing factor for loss of virulence, with no reduction caused by loss of siderophore production or melanisation. This reveals Ppt, and the lysine biosynthesis pathway, as potential targets for fungicides effective against Z. tritici.

Genome-wide associations for multiple pest resistances in a Northwestern United States elite spring wheat panel

Northern areas of the western United States are one of the most productive wheat growing regions in the United States. Increasing productivity through breeding is hindered by several biotic stresses which slow and constrain targeted yield improvement. In order to understand genetic variation for stripe rust (Puccinia striiformis f. sp. tritici), Septoria tritici blotch (Mycosphaerella graminicola), and Hessian fly (Mayetiola destructor) in regional germplasm, a panel of 408 elite spring wheat lines was characterized and genotyped with an Illumina 9K wheat single nucleotide polymorphism (SNP) chip to enable genome-wide association study (GWAS) analyses. Significant marker-trait associations were identified for stripe rust (38 loci), Septoria tritici blotch (8) and Hessian fly (9) resistance. Many of the QTL corresponded with previously reported gene locations or QTL, but we also discovered new resistance loci for each trait. We validated one of the stripe rust resistance loci detected by GWAS in a bi-parental mapping population, which confirmed the detection of Yr15 in the panel. This study elucidated well-defined chromosome regions for multiple pest resistances in elite Northwest germplasm. Newly identified resistance loci, along with SNPs more tightly linked to previously reported genes or QTL will help future breeding and marker assisted selection efforts.

Temporal Dynamics and Spatial Variation of Azoxystrobin and Propiconazole Resistance in Zymoseptoria tritici: A Hierarchical Survey of Commercial Winter Wheat Fields in the Willamette Valley, Oregon

Fungicide resistance can cause disease control failure in agricultural systems, and is particularly concerning with Zymoseptoria tritici, the causal agent of Septoria tritici blotch of wheat. In North America, the first quinone outside inhibitor resistance in Z. tritici was discovered in the Willamette Valley of Oregon in 2012, which prompted this hierarchical survey of commercial winter wheat fields to monitor azoxystrobin- and propiconazole-resistant Z. tritici. Surveys were conducted in June 2014, January 2015, May 2015, and January 2016. The survey was organized in a hierarchical scheme: regions within the Willamette Valley, fields within the region, transects within the field, and samples within the transect. Overall, frequency of azoxystrobin-resistant isolates increased from 63 to 93% from June 2014 to January 2016. Resistance to azoxystrobin increased over time even within fields receiving no strobilurin applications. Propiconazole sensitivity varied over the course of the study but, overall, did not significantly change. Sensitivity to both fungicides showed no regional aggregation within the Willamette Valley. Greater than 80% of spatial variation in fungicide sensitivity was at the smallest hierarchical scale (within the transect) of the survey for both fungicides, and the resistance phenotypes were randomly distributed within sampled fields. Results suggest a need for a better understanding of the dynamics of fungicide resistance at the landscape level.

How Knowledge of Pathogen Population Biology Informs Management of Septoria Tritici Blotch

Zymoseptoria tritici (previously Mycosphaerella graminicola) causes Septoria tritici blotch (STB) on wheat. The population biology of Z. tritici has been exceptionally well characterized as a result of intensive studies conducted over nearly 30 years. These studies provided important insights into the biology, epidemiology and evolutionary history of Z. tritici that will prove useful for management of STB. The well-documented, rapid adaptation of Z. tritici populations to fungicide applications and deployment of wheat cultivars carrying both major gene and quantitative resistance reflects the high evolutionary potential predicted by the large effective population size, high degree of gene flow and high levels of recombination found in field populations of Z. tritici globally. QST studies that assessed the global diversity for several important quantitative traits confirmed the adaptive potential of field populations and laid the groundwork for quantitative trait loci (QTL) mapping studies. QTL mapping elucidated the genetic architecture of each trait and led to identification of candidate genes affecting fungicide resistance, thermal adaptation, virulence, and host specialization. The insights that emerged through these analyses of Z. tritici population biology can now be used to generate actionable disease management strategies aimed at sustainably reducing losses due to STB. The high evolutionary potential found in field populations of Z. tritici requires deployment of a corresponding dynamically diverse set of control measures that integrate cultural, chemical, biological and resistance breeding strategies. In this review, we describe and prioritize STB control strategies based on current knowledge of Z. tritici population biology and propose a future research agenda oriented toward long-term STB management.

Reduced Virulence of Azoxystrobin-Resistant Zymoseptoria tritici Populations in Greenhouse Assays

The development of resistance to multiple fungicide classes is currently limiting disease management options for many pathogens, while the discovery of new fungicide classes may become less frequent. In light of this, more research is needed to quantify virulence trade-offs of fungicide resistance in order to more fully understand the implications of fungicide resistance on pathogen fitness. The purpose of this study was to measure the virulence of azoxystrobin-resistant and -sensitive Zymoseptoria tritici populations collected from North and South Willamette Valley, Oregon, in 2012 and 2015. Inoculum mixtures of known fungicide-resistant phenotypes were used to simulate natural field conditions, where multiple genotypes exist and interact in close proximity. Six greenhouse inoculations were conducted over 2 years, and virulence of the isolate mixtures was evaluated in planta. We considered virulence to be "the degree of pathology caused by the organism" and visually estimated the percent area of leaf necrosis as a measure of virulence. In greenhouse conditions, a consistent association of reduced virulence with azoxystrobin-resistant Z. tritici isolate mixtures was observed. North Willamette Valley and South Willamette Valley populations did not differ in virulence.

An Improved Method for Measuring Quantitative Resistance to the Wheat Pathogen Zymoseptoria tritici Using High-Throughput Automated Image Analysis

Zymoseptoria tritici causes Septoria tritici blotch (STB) on wheat. An improved method of quantifying STB symptoms was developed based on automated analysis of diseased leaf images made using a flatbed scanner. Naturally infected leaves (n = 949) sampled from fungicide-treated field plots comprising 39 wheat cultivars grown in Switzerland and 9 recombinant inbred lines (RIL) grown in Oregon were included in these analyses. Measures of quantitative resistance were percent leaf area covered by lesions, pycnidia size and gray value, and pycnidia density per leaf and lesion. These measures were obtained automatically with a batch-processing macro utilizing the image-processing software ImageJ. All phenotypes in both locations showed a continuous distribution, as expected for a quantitative trait. The trait distributions at both sites were largely overlapping even though the field and host environments were quite different. Cultivars and RILs could be assigned to two or more statistically different groups for each measured phenotype. Traditional visual assessments of field resistance were highly correlated with quantitative resistance measures based on image analysis for the Oregon RILs. These results show that automated image analysis provides a promising tool for assessing quantitative resistance to Z. tritici under field conditions.