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

Large-scale study validates that regional fungicide applications are major determinants of resistance evolution in the wheat pathogen Zymoseptoria tritici in France

In modern cropping systems, the near-universal use of plant protection products selects for resistance in pest populations. The emergence and evolution of this adaptive trait threaten treatment efficacy. We identified determinants of fungicide resistance evolution and quantified their effects at a large spatiotemporal scale. We focused on Zymoseptoria tritici, which causes leaf blotch in wheat. Phenotypes of qualitative or quantitative resistance to various fungicides were monitored annually, from 2004 to 2017, at about 70 sites throughout 22 regions of France (territorial units of 25 000 km² on average). We modelled changes in resistance frequency with regional anti-Septoria fungicide use, yield losses due to the disease and the regional area under organic wheat. The major driver of resistance dynamics was fungicide use at the regional scale. We estimated its effect on the increase in resistance and relative apparent fitness of each resistance phenotype. The predictions of the model replicated the spatiotemporal patterns of resistance observed in field populations (R² from 0.56 to 0.82). The evolution of fungicide resistance is mainly determined at the regional scale. This study therefore showed that collective management at the regional scale could effectively complete local actions.

Diversity and Virulence of Alternaria spp. Causing Potato Early Blight and Brown Spot in Wisconsin

Early blight, caused by Alternaria solani, along with brown spot, caused by A. alternata, have the potential to reduce quality and yield in potato production globally. Prior to this study, the incidence, disease impact, and fungicide resistance attributes of A. alternata in Wisconsin were poorly understood. Potato pathogens were isolated from foliar lesions at three commercial locations in Wisconsin in 2012 and 2017 and were initially morphologically identified as A. solani (n = 33) and A. alternata (n = 40). Identifications were further corroborated with the phylogenetic analysis of the internal transcribed spacer (ITS), translation elongation factor 1 (TEF1), gapdh, Alt a 1, and OPA10-2. A multigene phylogeny of ITS, TEF1, gapdh, and Alt a 1 showed five genotypes of A. alternata and one single genotype of A. solani. We demonstrated that the A. alternata isolates were virulent on potato cultivars Russet Burbank (P < 0.013) and Atlantic (P < 0.0073), though they caused less disease than A. solani (P < 0.0001 and P < 0.0001, respectively). A. alternata caused little disease on the breeding line 24-24-12 (P = 0.9929), and A. solani caused fewer disease symptoms on 24-24-12 than on Russet Burbank (P < 0.0001) or Atlantic (P < 0.0001). Breeding line 24-24-12 may be a promising source of potential resistance for the two diseases. There was no significant difference in virulence of different A. alternata genotypes, and no significant difference in virulence or genotype clustering among isolates from the three locations. Isolates of A. alternata that induced chlorosis caused larger lesion areas than isolates that did not in Russet Burbank (P < 0.0001), Atlantic (P < 0.0001), and 24-24-12 (P = 0.0365). There was no significant difference in virulence between quinone outside inhibitor (QoI)-sensitive and QoI-resistant isolates of A. alternata. This study enhanced our understanding of potato early blight and brown spot in Wisconsin, and suggested that A. alternata in addition to A. solani should be carefully monitored and possibly uniquely managed in order to achieve overall disease control.

Fitness Penalties in the Evolution of Fungicide Resistance

The evolution of resistance poses an ongoing threat to crop protection. Fungicide resistance provides a selective advantage under fungicide selection, but resistance-conferring mutations may also result in fitness penalties, resulting in an evolutionary trade-off. These penalties may result from the functional constraints of an evolving target site or from the resource allocation costs of overexpression or active transport. The extent to which such fitness penalties are present has important implications for resistance management strategies, determining whether resistance persists or declines between treatments, and for resistance risk assessments for new modes of action. Experimental results have proven variable, depending on factors such as temperature, nutrient status, osmotic or oxidative stress, and pathogen life-cycle stage. Functional genetics tools allow pathogen genetic background to be controlled, but this in turn raises the question of epistatic interactions. Combining fitness penalties under various conditions into a field-realistic scenario poses an important future challenge.

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.

Variable competitive effects of fungicide resistance in field experiments with a plant pathogenic fungus

Classic evolutionary theory suggests that mutations associated with antimicrobial and pesticide resistance result in a fitness cost in the absence of the selective antimicrobial agent or pesticide. There is experimental evidence to support fitness costs associated with resistance to anti-microbial compounds and pesticides across many biological disciplines, including human pathology, entomology, plant sciences, and plant pathology. However, researchers have also found examples of neutral and increased fitness associated with resistance, where the effect of a given resistance mutation depends on environmental and biological factors. We used Zymoseptoria tritici, a model evolutionary plant pathogenic fungus, to compare the competitive ability of fungicide-resistant isolates to fungicide-sensitive isolates. We conducted four large-scale inoculated winter wheat experiments at Oregon State University agriculture experiment stations. We found a significant change in the frequency of fungicide resistance over time in all four experiments. The direction and magnitude of these changes, however, differed by experimental location, year of experiment, and inoculum resistance treatment (fungicide-resistant, resistant/sensitive mixture, and fungicide-sensitive). These results suggest that the competitive ability of resistant isolates relative to sensitive isolates varied depending upon environmental conditions, including the initial frequency of resistant individuals in the population.

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.

Multilocus resistance evolution to azole fungicides in fungal plant pathogen populations

Evolution of fungicide resistance is a major threat to food production in agricultural ecosystems. Fungal pathogens rapidly evolved resistance to all classes of fungicides applied to the field. Resistance to the commonly used azole fungicides is thought to be driven mainly by mutations in a gene (CYP51) encoding a protein of the ergosterol biosynthesis pathway. However, some fungi gained azole resistance independently of CYP51 mutations and the mechanisms leading to CYP51-independent resistance are poorly understood. We used whole-genome sequencing and genome-wide association studies (GWAS) to perform an unbiased screen of azole resistance loci in Rhynchosporium commune, the causal agent of the barley scald disease. We assayed cyproconazole resistance in 120 isolates collected from nine populations worldwide. We found that mutations in highly conserved genes encoding the vacuolar cation channel YVC1, a transcription activator, and a saccharopine dehydrogenase made significant contributions to fungicide resistance. These three genes were not previously known to confer resistance in plant pathogens. However, YVC1 is involved in a conserved stress response pathway known to respond to azoles in human pathogenic fungi. We also performed GWAS to identify genetic polymorphism linked to fungal growth rates. We found that loci conferring increased fungicide resistance were negatively impacting growth rates, suggesting that fungicide resistance evolution imposed costs. Analyses of population structure showed that resistance mutations were likely introduced into local populations through gene flow. Multilocus resistance evolution to fungicides shows how pathogen populations can evolve a complex genetic architecture for an important phenotypic trait within a short time span.