Multilocus resistance evolution to azole fungicides in fungal plant pathogen populations

Target and non-target site mechanisms of fungicide resistance and their implications for the management of crop pathogens

Fungicides are indispensable for high-quality crops, but the rapid emergence and evolution of fungicide resistance have become the most important issues in modern agriculture. Hence, the sustainability and profitability of agricultural production have been challenged due to the limited number of fungicide chemical classes. Resistance to site-specific fungicides has principally been linked to target and non-target site mechanisms. These mechanisms change the structure or expression level, affecting fungicide efficacy and resulting in different and varying resistance levels. This review provides background information about fungicide resistance mechanisms and their implications for developing anti-resistance strategies in plant pathogens. Here, our purpose was to review changes at the target and non-target sites of quinone outside inhibitor (QoI) fungicides, methyl-benzimidazole carbamate (MBC) fungicides, demethylation inhibitor (DMI) fungicides, and succinate dehydrogenase inhibitor (SDHI) fungicides and to evaluate if they may also be associated with a fitness cost on crop pathogen populations. The current knowledge suggests that understanding fungicide resistance mechanisms can facilitate resistance monitoring and assist in developing anti-resistance strategies and new fungicide molecules to help solve this issue. © 2023 Society of Chemical Industry.

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.

A comprehensive review on potential applications of metallic nanoparticles as antifungal therapies to combat human fungal diseases

Human pathogenic fungi are responsible for causing a range of infection types including mucosal, skin, and invasive infections. Life-threatening and invasive fungal infections (FIs) are responsible for mortality and morbidity, especially for individuals with compromised immune function. The number of currently available therapeutic agents against invasive FIs is limited compared to that against bacterial infections. In addition, the increased mortality and morbidity caused by FIs are linked to the limited number of available antifungal agents, antifungal resistance, and the increased toxicity of these agents. Currently available antifungal agents have several drawbacks in efficiency, efficacy, toxicity, activity spectrum, and selectivity. It has already been demonstrated with numerous metallic nanoparticles (MNPs) that these nanoparticles can serve as an effective and alternative solution as fungicidal agents. MNPs have great potential owing to their intrinsic antifungal properties and potential to deliver antifungal drugs. For instance, gold nanoparticles (AuNPs) have the capacity to disturb mitochondrial calcium homeostasis induced AuNP-mediated cell death in Candida albicans. In addition, both copper nanoparticles and copper oxide nanoparticles exerted significant suppressive properties against pathogenic fungi. Silver nanoparticles showed strong antifungal properties against numerous pathogenic fungi, such as Stachybotrys chartarum, Mortierella alpina, Chaetomium globosum, A. fumigatus, Cladosporium cladosporioides, Penicillium brevicompactum, Trichophyton rubrum, C. tropicalis, and C. albicans. Iron oxide nanoparticles showed potent antifungal activities against A. niger and P. chrysogenum. It has also been reported that zinc oxide nanoparticles can significantly inhibit fungal growth. These NPs have already exerted potent antifungal properties against a number of pathogenic fungal species including Candida, Aspergillus, Fusarium, and many others. Several strategies are currently used for the research and development of antifungal NPs including chemical modification of NPs and combination with the available drugs. This review has comprehensively presented the current and innovative antifungal approach using MNPs. Moreover, different types of MNPs, their physicochemical characteristics, and production techniques have been summarized in this review.

Emerging Polymer-Based Nanosystem Strategies in the Delivery of Antifungal Drugs

Nanosystems-based antifungal agents have emerged as an effective strategy to address issues related to drug resistance, drug release, and toxicity. Among the diverse materials employed for antifungal drug delivery, polymers, including polysaccharides, proteins, and polyesters, have gained significant attention due to their versatility. Considering the complex nature of fungal infections and their varying sites, it is crucial for researchers to carefully select appropriate polymers based on specific scenarios when designing antifungal agent delivery nanosystems. This review provides an overview of the various types of nanoparticles used in antifungal drug delivery systems, with a particular emphasis on the types of polymers used. The review focuses on the application of drug delivery systems and the release behavior of these systems. Furthermore, the review summarizes the critical physical properties and relevant information utilized in antifungal polymer nanomedicine delivery systems and briefly discusses the application prospects of these systems.

Synergism: biocontrol agents and biostimulants in reducing abiotic and biotic stresses in crop

In today's fast-shifting climate change scenario, crops are exposed to environmental pressures, abiotic and biotic stress. Hence, these will affect the production of agricultural products and give rise to a worldwide economic crisis. The increase in world population has exacerbated the situation with increasing food demand. The use of chemical agents is no longer recommended due to adverse effects towards the environment and health. Biocontrol agents (BCAs) and biostimulants, are feasible options for dealing with yield losses induced by plant stresses, which are becoming more intense due to climate change. BCAs and biostimulants have been recommended due to their dual action in reducing both stresses simultaneously. Although protection against biotic stresses falls outside the generally accepted definition of biostimulant, some microbial and non-microbial biostimulants possess the biocontrol function, which helps reduce biotic pressure on crops. The application of synergisms using BCAs and biostimulants to control crop stresses is rarely explored. Currently, a combined application using both agents offer a great alternative to increase the yield and growth of crops while managing stresses. This article provides an overview of crop stresses and plant stress responses, a general knowledge on synergism, mathematical modelling used for synergy evaluation and type of in vitro and in vivo synergy testing, as well as the application of synergism using BCAs and biostimulants in reducing crop stresses. This review will facilitate an understanding of the combined effect of both agents on improving crop yield and growth and reducing stress while also providing an eco-friendly alternative to agroecosystems.

Heritability and gene functions associated with sclerotia formation of Rhizoctonia solani AG-7 using whole genome sequencing and genome-wide association study

Sclerotia are specialized fungal structures formed by pigmented and aggregated hyphae, which can survive under unfavourable environmental conditions and serve as the primary inocula for several phytopathogenic fungi including Rhizoctonia solani. Among 154 R. solani anastomosis group 7 (AG-7) isolates collected in fields, the sclerotia-forming capability regarding sclerotia number and sclerotia size varied in the fungal population, but the genetic makeup of these phenotypes remained unclear. As limited studies have focused on the genomics of R. solani AG-7 and the population genetics of sclerotia formation, this study completed the whole genome sequencing and gene prediction of R. solani AG-7 using the Oxford NanoPore and Illumina RNA sequencing. Meanwhile, a high-throughput image-based method was established to quantify the sclerotia-forming capability, and the phenotypic correlation between sclerotia number and sclerotia size was low. A genome-wide association study identified three and five significant SNPs associated with sclerotia number and size in distinct genomic regions, respectively. Of these significant SNPs, two and four showed significant differences in the phenotypic mean separation for sclerotia number and sclerotia size, respectively. Gene ontology enrichment analysis focusing on the linkage disequilibrium blocks of significant SNPs identified more categories related to oxidative stress for sclerotia number, and more categories related to cell development, signalling and metabolism for sclerotia size. These results indicated that different genetic mechanisms may underlie these two phenotypes. Moreover, the heritability of sclerotia number and sclerotia size were estimated for the first time to be 0.92 and 0.31, respectively. This study provides new insights into the heritability and gene functions related to the development of sclerotia number and sclerotia size, which could provide additional knowledge to reduce fungal residues in fields and achieve sustainable disease management.

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.

The population genetics of adaptation through copy‐number variation in a fungal plant pathogen

Microbial pathogens can adapt rapidly to changing environments such as the application of pesticides or host resistance. Copy number variations (CNVs) are a major source of adaptive genetic variation for recent adaptation. Here, we analyse how a major fungal pathogen of barley, Rhynchosporium commune, has adapted to the host environment and fungicide applications. We screen the genomes of 125 isolates sampled across a worldwide set of populations and identify a total of 7,879 gene duplications and 116 gene deletions. Most gene duplications result from segmental chromosomal duplications. Although CNVs are generally under negative selection, we find that genes affected by CNVs are enriched in functions related to host exploitation (i.e., effectors and cell-wall-degrading enzymes). We perform genome-wide association studies (GWAS) and identify a large segmental duplication of CYP51A that has contributed to the emergence of azole resistance and a duplication encompassing an effector gene affecting virulence. We show that the adaptive CNVs were probably created by recently active transposable element families. Moreover, we find that specific transposable element families are important drivers of recent gene CNV. Finally, we use a genome-wide single nucleotide polymorphism data set to replicate the GWAS and contrast it with the CNV-focused analysis. Together, our findings show how extensive segmental duplications create the raw material for recent adaptation in global populations of a fungal pathogen.

Striking Back against Fungal Infections: The Utilization of Nanosystems for Antifungal Strategies

Fungal infections have become a major health concern, given that invasive infections by Candida, Cryptococcus, and Aspergillus species have led to millions of mortalities. Conventional antifungal drugs including polyenes, echinocandins, azoles, allylamins, and antimetabolites have been used for decades, but their limitations include off-target toxicity, drug-resistance, poor water solubility, low bioavailability, and weak tissue penetration, which cannot be ignored. These drawbacks have led to the emergence of novel antifungal therapies. In this review, we discuss the nanosystems that are currently utilized for drug delivery and the application of antifungal therapies.

Genome-Wide Association and Selective Sweep Studies Reveal the Complex Genetic Architecture of DMI Fungicide Resistance in Cercospora beticola

The rapid and widespread evolution of fungicide resistance remains a challenge for crop disease management. The demethylation inhibitor (DMI) class of fungicides is a widely used chemistry for managing disease, but there has been a gradual decline in efficacy in many crop pathosystems. Reliance on DMI fungicides has increased resistance in populations of the plant pathogenic fungus Cercospora beticola worldwide. To better understand the genetic and evolutionary basis for DMI resistance in C. beticola, a genome-wide association study (GWAS) and selective sweep analysis were conducted for the first time in this species. We performed whole-genome resequencing of 190 C. beticola isolates infecting sugar beet (Beta vulgaris ssp. vulgaris). All isolates were phenotyped for sensitivity to the DMI tetraconazole. Intragenic markers on chromosomes 1, 4, and 9 were significantly associated with DMI fungicide resistance, including a polyketide synthase gene and the gene encoding the DMI target CbCYP51. Haplotype analysis of CbCYP51 identified a synonymous mutation (E170) and nonsynonymous mutations (L144F, I387M, and Y464S) associated with DMI resistance. Genome-wide scans of selection showed that several of the GWAS mutations for fungicide resistance resided in regions that have recently undergone a selective sweep. Using radial plate growth on selected media as a fitness proxy, we did not find a trade-off associated with DMI fungicide resistance. Taken together, we show that population genomic data from a crop pathogen can allow the identification of mutations conferring fungicide resistance and inform about their origins in the pathogen population.

Assessment and modeling using machine learning of resistance to scald (Rhynchosporium commune) in two specific barley genetic resources subsets

Barley production worldwide is limited by several abiotic and biotic stresses and breeding of highly productive and adapted varieties is key to overcome these challenges. Leaf scald, caused by Rhynchosporium commune is a major disease of barley that requires the identification of novel sources of resistance. In this study two subsets of genebank accessions were used: one extracted from the Reference set developed within the Generation Challenge Program (GCP) with 191 accessions, and the other with 101 accessions selected using the filtering approach of the Focused Identification of Germplasm Strategy (FIGS). These subsets were evaluated for resistance to scald at the seedling stage under controlled conditions using two Moroccan isolates, and at the adult plant stage in Ethiopia and Morocco. The results showed that both GCP and FIGS subsets were able to identify sources of resistance to leaf scald at both plant growth stages. In addition, the test of independence and goodness of fit showed that FIGS filtering approach was able to capture higher percentages of resistant accessions compared to GCP subset at the seedling stage against two Moroccan scald isolates, and at the adult plant stage against four field populations of Morocco and Ethiopia, with the exception of Holetta nursery 2017. Furthermore, four machine learning models were tuned on training sets to predict scald reactions on the test sets based on diverse metrics (accuracy, specificity, and Kappa). All models efficiently identified resistant accessions with specificities higher than 0.88 but showed different performances between isolates at the seedling and to field populations at the adult plant stage. The findings of our study will help in fine-tuning FIGS approach using machine learning for the selection of best-bet subsets for resistance to scald disease from the large number of genebank accessions.

Machine-learning predicts genomic determinants of meiosis-driven structural variation in a eukaryotic pathogen

Species harbor extensive structural variation underpinning recent adaptive evolution. However, the causality between genomic features and the induction of new rearrangements is poorly established. Here, we analyze a global set of telomere-to-telomere genome assemblies of a fungal pathogen of wheat to establish a nucleotide-level map of structural variation. We show that the recent emergence of pesticide resistance has been disproportionally driven by rearrangements. We use machine learning to train a model on structural variation events based on 30 chromosomal sequence features. We show that base composition and gene density are the major determinants of structural variation. Retrotransposons explain most inversion, indel and duplication events. We apply our model to Arabidopsis thaliana and show that our approach extends to more complex genomes. Finally, we analyze complete genomes of haploid offspring in a four-generation pedigree. Meiotic crossover locations are enriched for new rearrangements consistent with crossovers being mutational hotspots. The model trained on species-wide structural variation accurately predicts the position of >74% of newly generated variants along the pedigree. The predictive power highlights causality between specific sequence features and the induction of chromosomal rearrangements. Our work demonstrates that training sequence-derived models can accurately identify regions of intrinsic DNA instability in eukaryotic genomes.

Genome-Wide Association Analysis for Triazole Resistance in Aspergillus fumigatus

Aspergillus fumigatus is a ubiquitous fungus and the main agent of aspergillosis, a common fungal infection in the immunocompromised population. Triazoles such as itraconazole and voriconazole are the common first-line drugs for treating aspergillosis. However, triazole resistance in A. fumigatus has been reported in an increasing number of countries. While most studies of triazole resistance have focused on mutations in the triazole target gene cyp51A, >70% of triazole-resistant strains in certain populations showed no mutations in cyp51A. To identify potential non-cyp51A mutations associated with triazole resistance in A. fumigatus, we analyzed the whole genome sequences and triazole susceptibilities of 195 strains from 12 countries. These strains belonged to three distinct clades. Our genome-wide association study (GWAS) identified a total of six missense mutations significantly associated with itraconazole resistance and 18 missense mutations with voriconazole resistance. In addition, to investigate itraconazole and pan-azole resistance, Fisher's exact tests revealed 26 additional missense variants tightly linked to the top 20 SNPs obtained by GWAS, of which two were consistently associated with triazole resistance. The large number of novel mutations related to triazole resistance should help further investigations into their molecular mechanisms, their clinical importance, and the development of a comprehensive molecular diagnosis toolbox for triazole resistance in A. fumigatus.

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.

Increasing temperature elevates the variation and spatial differentiation of pesticide tolerance in a plant pathogen

Climate change and pesticide resistance are two of the most imminent challenges human society is facing today. Knowledge of how the evolution of pesticide resistance may be affected by climate change such as increasing air temperature on the planet is important for agricultural production and ecological sustainability in the future but is lack in scientific literatures reported from empirical research. Here, we used the azoxystrobin-Phytophthora infestans interaction in agricultural systems to investigate the contributions of environmental temperature to the evolution of pesticide resistance and infer the impacts of global warming on pesticide efficacy and future agricultural production and ecological sustainability. We achieved this by comparing azoxystrobin sensitivity of 180 P. infestans isolates sampled from nine geographic locations in China under five temperature schemes ranging from 13 to 25°C. We found that local air temperature contributed greatly to the difference of azoxystrobin tolerance among geographic populations of the pathogen. Both among-population and within-population variations in azoxystrobin tolerance increased as experimental temperatures increased. We also found that isolates with higher azoxystrobin tolerance adapted to a broader thermal niche. These results suggest that global warming may enhance the risk of developing pesticide resistance in plant pathogens and highlight the increased challenges of administering pesticides for effective management of plant diseases to support agricultural production and ecological sustainability under future thermal conditions.

Photodynamic inactivation of Lasiodiplodia theobromae: lighting the way towards an environmentally friendly phytosanitary treatment

The fungus Lasiodiplodia theobromae is one of the main causal agents of trunk canker and dieback of grapevine. The objective of this work was to evaluate the efficiency of photodynamic inactivation (PDI) of L. theobromae with synthetic and natural photosensitizers and irradiation with either sunlight or artificial photosynthetically active radiation. Although the growth of the mycelium could not be completely prevented with natural sunlight irradiation, phenothiazine dyes (methylene blue, MB; toluidine blue O, TBO), riboflavin and a cationic porphyrin (Tetra-Py⁺-Me) caused complete inhibition under continuous irradiation with artificial light. Free radicals were the main cytotoxic agents in the PDI with MB, indicating the predominance of the type I mechanism. PDI with MB or Tetra-Py⁺-Me may represent a promising approach for the sanitation of vine material in greenhouse nurseries, in order to reduce the risk of infection upon grafting.

Non-Target Site Mechanisms of Fungicide Resistance in Crop Pathogens: A Review

The rapid emergence of resistance in plant pathogens to the limited number of chemical classes of fungicides challenges sustainability and profitability of crop production worldwide. Understanding mechanisms underlying fungicide resistance facilitates monitoring of resistant populations at large-scale, and can guide and accelerate the development of novel fungicides. A majority of modern fungicides act to disrupt a biochemical function via binding a specific target protein in the pathway. While target-site based mechanisms such as alternation and overexpression of target genes have been commonly found to confer resistance across many fungal species, it is not uncommon to encounter resistant phenotypes without altered or overexpressed target sites. However, such non-target site mechanisms are relatively understudied, due in part to the complexity of the fungal genome network. This type of resistance can oftentimes be transient and noninheritable, further hindering research efforts. In this review, we focused on crop pathogens and summarized reported mechanisms of resistance that are otherwise related to target-sites, including increased activity of efflux pumps, metabolic circumvention, detoxification, standing genetic variations, regulation of stress response pathways, and single nucleotide polymorphisms (SNPs) or mutations. In addition, novel mechanisms of drug resistance recently characterized in human pathogens are reviewed in the context of nontarget-directed resistance.

The complex genomic basis of rapid convergent adaptation to pesticides across continents in a fungal plant pathogen

Convergent evolution leads to identical phenotypic traits in different species or populations. Convergence can be driven by standing variation allowing selection to favour identical alleles in parallel or the same mutations can arise independently. However, the molecular basis of such convergent adaptation remains often poorly resolved. Pesticide resistance in agricultural ecosystems is a hallmark of convergence in phenotypic traits. Here, we analyse the major fungal pathogen Zymoseptoria tritici causing serious losses on wheat and with fungicide resistance emergence across several continents. We sampled three population pairs each from a different continent spanning periods early and late in the application of fungicides. To identify causal loci for resistance, we combined knowledge from molecular genetics work and performed genome-wide association studies (GWAS) on a global set of isolates. We discovered yet unknown factors in azole resistance including a gene encoding membrane associated functions. We found strong support for the "hotspot" model of resistance evolution with convergent changes in a small set of loci but additional loci showed more population-specific allele frequency changes. Genome-wide scans of selection showed that half of all known resistance loci were overlapping a selective sweep region. Hence, the application of fungicides was one of the major selective agents acting on the pathogen over the past decades. Furthermore, loci identified through GWAS showed the highest overlap with selective sweep regions underlining the importance to map phenotypic trait variation in evolving populations. Our population genomic analyses highlighted that both de novo mutations and gene flow contributed to convergent pesticide adaptation.

Fungicide SYP-14288 Inducing Multidrug Resistance in Rhizoctonia solani

Rhizoctonia solani is a widely distributed soilborne plant pathogen, and can cause significant economic losses to crop production. In chemical controls, SYP-14288 is highly effective against plant pathogens, including R. solani. To examine the sensitivity to SYP-14288, 112 R. solani isolates were collected from infected rice plants. An established baseline sensitivity showed that values of effective concentration for 50% growth inhibition (EC₅₀) ranged from 0.0003 to 0.0138 μg/ml, with an average of 0.0055 ± 0.0030 μg/ml. The frequency distribution of the EC₅₀ was unimodal and the range of variation factor (the ratio of maximal over minimal EC₅₀) was 46.03, indicating that all wild-type strains were sensitive to SYP-14288. To examine the risk of fungicide resistance, 20 SYP-14288-resistant mutants were generated on agar plates amended with SYP-14288. Eighteen mutants remained resistant after 10 transfers, and their fitness was significantly different from the parental strain. All of the mutants grew more slowly but showed high virulence to rice plants, though lower than the parental strain. A cross-resistance assay demonstrated that there was a positive correlation between SYP-14288 and fungicides having or not having the same mode of action with SYP-14288, including fluazinam, fentin chloride, fludioxonil, difenoconazole, cyazofamid, chlorothalonil, and 2,4-dinitrophen. This result showed a multidrug resistance induced by SYP-14288, which could be a concern in increasing the spectrum of resistance in R. solani to commonly used fungicides.

The Genetic Architecture of Emerging Fungicide Resistance in Populations of a Global Wheat Pathogen

Containing fungal diseases often depends on the application of fungicidal compounds. Fungicides can rapidly lose effectiveness due to the rise of resistant individuals in populations. However, the lack of knowledge about resistance mutations beyond known target genes challenges investigations into pathways to resistance. We used whole-genome sequencing data and association mapping to reveal the multilocus genetic architecture of fungicide resistance in a global panel of 159 isolates of Parastagonospora nodorum, an important fungal pathogen of wheat. We found significant differences in azole resistance among global field populations. The populations evolved distinctive combinations of resistance alleles which can interact when co-occurring in the same genetic background. We identified 34 significantly associated single nucleotide polymorphisms located in close proximity to genes associated with fungicide resistance in other fungi, including a major facilitator superfamily transporter. Using fungal colony growth rates and melanin production at different temperatures as fitness proxies, we found no evidence that resistance was constrained by genetic trade-offs. Our study demonstrates how genome-wide association studies of a global collection of pathogen strains can recapitulate the emergence of fungicide resistance. The distinct complement of resistance mutations found among populations illustrates how the evolutionary trajectory of fungicide adaptation can be complex and challenging to predict.

Recent insights into barley and Rhynchosporium commune interactions

Rhynchosporium commune is the causal pathogen of scald in barley (Hordeum vulgare), a foliar disease that can reduce yield by up to 40% in susceptible cultivars. R. commune is found worldwide in all temperate growing regions and is regarded as one of the most economically important barley pathogens. It is a polycyclic pathogen with the ability to rapidly evolve new virulent strains in response to resistance genes deployed in commercial cultivars. Hence, introgression and pyramiding of different loci for resistance (qualitative or quantitative) through marker-assisted selection is an effective way to improve scald resistance in barley. This review summarizes all 148 resistance quantitative trait loci reported at the date of submission of this review and projects them onto the barley physical map, where it is clear many loci co-locate on chromosomes 3H and 7H. We have summarized the major named resistance loci and reiterated the renaming of Rrs15 (CI8288) to Rrs17. This review provides a comprehensive resource for future discovery and breeding efforts of qualitative and quantitative scald resistance loci.

Temperature-Mediated Plasticity Regulates the Adaptation of Phytophthora infestans to Azoxystrobin Fungicide

Fungicide is one of the main approaches used in agriculture to manage plant diseases for food production, but their effectiveness can be reduced due to the evolution of plant pathogens. Understanding the genetics and evolutionary processes responsible for the development of fungicide resistance is a key to food production and social sustainability. In this study, we used a common garden experiment to examine the source of genetic variation, natural selection, and temperature contributing to the development of azoxystrobin resistance in Phytophthora infestans and infer sustainable ways of plant disease management in future. We found that plasticity contributed to ~40% of phenotypic variation in azoxystrobin sensitivity while heritability accounted for 16%. Further analysis indicated that overall population differentiation in azoxystrobin sensitivity (QST) was significantly greater than the overall population differentiation in simple sequence repeat (SSR) marker (FST), and the P. infestans isolates demonstrated higher level of azoxystrobin sensitivity at the higher experimental temperature. These results suggest that changes in target gene expression, enzymatic activity, or metabolic rate of P. infestans play a more important role in the adaptation of the pathogen to azoxystrobin resistance than that of mutations in target genes. The development of azoxystrobin resistance in P. infestans is likely driven by diversifying selection for local adaptation, and elevated temperature associated with global warming in the future may increase the effectiveness of using azoxystrobin to manage P. infestans. The sustainable approaches for increasing disease control effectiveness and minimizing the erosion of the fungicide efficacy are proposed.

Out in the Cold: Identification of Genomic Regions Associated With Cold Tolerance in the Biocontrol Fungus Clonostachys rosea Through Genome-Wide Association Mapping

There is an increasing importance for using biocontrol agents in combating plant diseases sustainably and in the long term. As large scale genomic sequencing becomes economically viable, the impact of single nucleotide polymorphisms (SNPs) on biocontrol-associated phenotypes can be easily studied across entire genomes of fungal populations. Here, we improved a previously reported genome assembly of the biocontrol fungus Clonostachys rosea strain IK726 using the PacBio sequencing platform, which resulted in a total genome size of 70.7 Mbp and 21,246 predicted genes. We further performed whole-genome re-sequencing of 52 additional C. rosea strains isolated globally using Illumina sequencing technology, in order to perform genome-wide association studies in conditions relevant for biocontrol activity. One such condition is the ability to grow at lower temperatures commonly encountered in cryic or frigid soils in temperate regions, as these will be prevalent for protecting growing crops in temperate climates. Growth rates at 10°C on potato dextrose agar of the 53 sequenced strains of C. rosea were measured and ranged between 0.066 and 0.413 mm/day. Performing a genome wide association study, a total of 1,478 SNP markers were significantly associated with the trait and located in 227 scaffolds, within or close to (< 1000 bp distance) 265 different genes. The predicted gene products included several chaperone proteins, membrane transporters, lipases, and proteins involved in chitin metabolism with possible roles in cold tolerance. The data reported in this study provides a foundation for future investigations into the genetic basis for cold tolerance in fungi, with important implications for biocontrol.

Low evolutionary risk of iprovalicarb resistance in Phytophthora infestans

Iprovalicarb is a carboxylic acid amide (CAA) fungicide, highly effective in the control of potato late blight, causing by Phytophthora infestans. Due to cross-resistance with other CAA fungicides and moderate resistance risk of P. capsici to iprovalicarb, the evolutionary risk of P. infestans resistance to this fungicide and the contribution of inherited genes and environmental effect was evaluated using a common garden experiment. The results showed that the ratio of heritability and plasticity of iprovalicarb in the seven populations equaled 1.0, indicating both inherited genes and environmental factors were essential to iprovalicarb sensitivity in P. infestans. The pairwise population differentiation determined by SSR loci (F_ST) between populations ranged from 0.007 to 0.133 and the overall F_ST was significantly higher than population differentiation in RGR (Q_ST), suggesting constraining selection acting on iprovalicarb sensitivity. We also found a new indicator of growth rate inhibition (GRI) for fungicide sensitivity, which was negatively correlated to growth rate in the absence of iprovalicarb, indicating a trade-off between iprovalicarb resistance and pathogen's growth. The constraining selection plus a trade-off between GRI and growth rate revealed low risk of P. infestans evolving resistance to iprovalicarb.

Genetic Diversity and Population Structure of Pythium irregulare from Soybean and Corn Production Fields in Ohio

High levels of genetic diversity have been described within the Pythium irregulare complex from several host plants; however, little is known about the population structure in fields used for grain production. Therefore, the objective of this study was to evaluate the genetic diversity and population structure of 53 isolates baited from 28 soybean and corn production fields from 25 counties in Ohio. Genetic diversity was characterized based on sequence analysis of the internal transcribed spacer (ITS1-5.8S-ITS2) region and with 21 simple sequence repeat (SSR) markers. In addition, aggressiveness on soybean, optimum growth temperature, and sensitivity to metalaxyl fungicide were determined. ITS sequence analysis indicated that four isolates clustered with P. cryptoirregulare, whereas the remaining isolates clustered with P. irregulare that was subdivided into two groups (1 and 2). Cluster analysis of SSR data revealed a similar subdivision, which was also supported by structure analysis. The isolates from group 2 grew at a slower rate, but both groups of P. irregulare and P. cryptoirregulare recovered in this study had the same optimum growth at 27°C. Variability of aggressiveness and sensitivity toward metalaxyl fungicide was also observed among isolates within each group. The results from this study will help in the selection of isolates to be used in screening for resistance, assessment of fungicide efficacy, and disease management recommendations.

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.

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.

Genome-Wide Detection of Genes Under Positive Selection in Worldwide Populations of the Barley Scald Pathogen

Coevolution between hosts and pathogens generates strong selection pressures to maintain resistance and infectivity, respectively. Genomes of plant pathogens often encode major effect loci for the ability to successfully infect specific host genotypes. Hence, spatial heterogeneity in host genotypes coupled with abiotic factors could lead to locally adapted pathogen populations. However, the genetic basis of local adaptation is poorly understood. Rhynchosporium commune, the pathogen causing barley scald disease, interacts at least partially in a gene-for-gene manner with its host. We analyzed global field populations of 125 R. commune isolates to identify candidate genes for local adaptation. Whole genome sequencing data showed that the pathogen is subdivided into three genetic clusters associated with distinct geographic and climatic regions. Using haplotype-based selection scans applied independently to each genetic cluster, we found strong evidence for selective sweeps throughout the genome. Comparisons of loci under selection among clusters revealed little overlap, suggesting that ecological differences associated with each cluster led to variable selection regimes. The strongest signals of selection were found predominantly in the two clusters composed of isolates from Central Europe and Ethiopia. The strongest selective sweep regions encoded protein functions related to biotic and abiotic stress responses. Selective sweep regions were enriched in genes encoding functions in cellular localization, protein transport activity, and DNA damage responses. In contrast to the prevailing view that a small number of gene-for-gene interactions govern plant pathogen evolution, our analyses suggest that the evolutionary trajectory is largely determined by spatially heterogeneous biotic and abiotic selection pressures.

Slow and temperature-mediated pathogen adaptation to a nonspecific fungicide in agricultural ecosystem

The spread of antimicrobial resistance and global change in air temperature represent two major phenomena that are exerting a disastrous impact on natural and social issues but investigation of the interaction between these phenomena in an evolutionary context is limited. In this study, a statistical genetic approach was used to investigate the evolution of antimicrobial resistance in agricultural ecosystem and its association with local air temperature, precipitation, and UV radiation. We found no resistance to mancozeb, a nonspecific fungicide widely used in agriculture for more than half a century, in 215 Alternaria alternata isolates sampled from geographic locations along a climatic gradient and cropping system representing diverse ecotypes in China, consistent with low resistance risk in many nonspecific fungicides. Genetic variance accounts for ~35% of phenotypic variation, while genotype-environment interaction is negligible, suggesting that heritability plays a more important role in the evolution of resistance to mancozeb in plant pathogens than phenotypic plasticity. We also found that tolerance to mancozeb in agricultural ecosystem is under constraining selection and significantly associated with local air temperature, possibly resulting from a pleiotropic effect of resistance with thermal and other ecological adaptations. The implication of these results for fungicide and other antimicrobial management in the context of global warming is discussed.

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.

A Cationic Polymer That Shows High Antifungal Activity against Diverse Human Pathogens

Invasive fungal diseases are generally difficult to treat and often fatal. The therapeutic agents available to treat fungi are limited, and there is a critical need for new agents to combat these deadly infections. Antifungal compound development has been hindered by the challenge of creating agents that are highly active against fungal pathogens but not toxic to the host. Host defense peptides (HDPs) are produced by eukaryotes as a component of the innate immune response to pathogens and have served as inspiration for the development of many new antibacterial compounds. HDP mimics, however, have largely failed to exhibit potent and selective antifungal activity. Here, we present an HDP-like nylon-3 copolymer that is effective against diverse fungi while displaying only mild to moderate toxicity toward mammalian cells. This polymer is active on its own and in synergy with existing antifungal drugs against multiple species of Candida and Cryptococcus, reaching levels of efficacy comparable to those of the clinical agents amphotericin B and fluconazole in some cases. In addition, the polymer acts synergistically with azoles against different species of Aspergillus, including some azole-resistant strains. These findings indicate that nylon-3 polymers are a promising lead for development of new antifungal therapeutic strategies.

High-density genetic mapping identifies the genetic basis of a natural colony morphology mutant in the root rot pathogen Armillaria ostoyae

Filamentous fungi exhibit a broad spectrum of heritable growth patterns and morphological variations reflecting the adaptation of the different species to distinct ecological niches. But also within species, isolates show considerable variation in growth rates and other morphological characteristics. The genetic basis of this intraspecific variation in mycelial growth and morphology is currently poorly understood. By chance, a growth mutant in the root rot pathogen Armillaria ostoyae was discovered. The mutant phenotype was characterized by extremely compact and slow growth, as well as shorter aerial hyphae and hyphal compartments in comparison to the wildtype phenotype. Genetic analysis revealed that the abnormal phenotype is caused by a recessive mutation, which segregates asa single locus in sexual crosses. In order to identify the genetic basis of the mutant phenotype, we performed a quantitative trait locus (QTL) analysis. A mapping population of 198 haploid progeny was genotyped at 11,700 genome-wide single nucleotide polymorphisms (SNPs) making use of double digest restriction site associated DNA sequencing (ddRADseq). In accordance with the genetic analysis, a single significant QTL was identified for the abnormal growth phenotype. The QTL confidence interval spans a narrow, gene dense region of 87kb in the A. ostoyae genome which contains 37 genes. Overall, our study reports the first high-density genetic map for an Armillaria species and shows its successful application in forward genetics by resolving the genetic basis of a mutant phenotype with a severe defect in hyphal growth.