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

Multi-computational screening identifies homovanillic acid as a potential SAP5 inhibitor against Candida albicans biofilms

This work aims to find inhibitors of SAP5, a virulence factor in Candida albicans polymicrobial biofilms. The methodology included docking simulations, MMGBSA calculations, and molecular dynamics simulations. Of the 107 phenolic acids retrieved from PubChem, 20 passed ADMET screening. The research finds homovanillic acid to be a possible SAP5 inhibitor, with a binding energy of -19.92 kcal/mol as shown by molecular docking and MMGBSA analysis. The compound showed favorable ADMET properties, indicating low toxicity and high drug-likeness. Molecular dynamics simulations over 100 nanoseconds confirmed stable protein-ligand interactions. These findings suggest homovanillic acid's potential in treating AMR-associated biofilms and establish a foundation for experimental validation. The study demonstrates how computational methods can accelerate the discovery of novel antifungal medicines targeting polymicrobial infections.

Genome-wide Evidence of Host Specialization in Wild and Farmland Populations of the Fungal Leaf Spot Pathogen, Cercospora beticola

One of the most recent crop species to be domesticated is sugar beet (Beta vulgaris L. ssp. vulgaris Doell.), which was bred for high sucrose content within the last few centuries in Europe. Crop domestication can also lead to the evolution of novel pathogens, which may spread across large geographical distances with their crop host. In this study, we addressed the recent evolution of the fungal pathogen causing the disease Cercospora leaf spot, Cercospora beticola. This pathogen has become increasingly important in sugar beet and table beet production worldwide. We used genome sequences of 326 C. beticola isolates collected from 4 continents from 4 closely related Beta subspecies (3 domesticated and 1 wild). We applied population genomic analyses to identify signatures of population differentiation and host specialization in C. beticola populations derived from the cultivated and wild hosts. We found evidence that C. beticola populations in agro-ecosystems likely originate from sea beet-infecting isolates. Intriguingly, host jumps from wild to cultivated beet occurred in at least 2 independent events as evidenced by our population data of C. beticola from wild beet collected in the Mediterranean and the UK. We explore the occurrence of genetic variants associated with fungicide resistance and virulence and show that standing genetic variation in C. beticola populations from both wild and domesticated plants may serve as a reservoir of functionally important alleles. Overall, our results highlight the ability of C. beticola to invade the agro-ecosystem and establish new populations, demonstrating the rapid adaptation potential of the species.

Genomic Surveillance and Molecular Evolution of Fungicide Resistance in European Populations of Wheat Powdery Mildew

Fungicides are used in agriculture to manage fungal infections and maintain crop yield and quality. In Europe, their application on cereals increased drastically starting from the mid 1970s, contributing to a significant improvement in yields. However, extensive usage has led to the rapid evolution of resistant pathogen populations within just a few years of fungicide deployment. Here we focus on wheat powdery mildew, a disease caused by the ascomycete fungus Blumeria graminis forma specialis tritici (Bgt). Previous research on Bgt documented the emergence of resistance to different fungicides and identified various resistance mechanisms. Yet, the frequency, distribution and evolutionary dynamics of fungicide resistance in Bgt populations remain largely unexplored. In this study, we leveraged extensive sampling and whole-genome sequencing of Bgt populations in Europe and the Mediterranean to investigate the population genetics and molecular epidemiology of fungicide resistance towards five major fungicide classes. We analysed gene sequences and copy number variation of eight known fungicide target genes in 415 Bgt isolates sampled between 1980 and 2023. We observed that mutations conferring resistance to various fungicides increased in frequency over time and had distinct geographic distributions, probably due to diverse deployment of fungicides across different regions. For demethylation inhibitor fungicides, we identified multiple independent events of resistance emergence with distinct mutational profiles, and we tracked their rapid spread in the last decades. Overall, we revealed the evolutionary and epidemiological dynamics of fungicide resistance mutations in European Bgt populations. These results underscore the potential of genomic surveillance and population genetics to enhance our understanding of fungicide resistance.

Molecular Monitoring of Fungicide Resistance in Crop Phytopathogens

The fight against crop pathogens relies mainly on host genetics and chemistry; however, both areas are compromised by the evolution of resistance in the pathogen population. Fungicide resistance is an ongoing challenge to global food security, as it threatens these important crop protection chemistries. One critical component of resistance management is an effective detection and monitoring program, which needs to be agile, scalable, sensitive, accurate, and cost effective. A rapidly evolving suite of molecular tools are being developed for the detection of fungicide resistance mutations in phytopathogen populations, including high-throughput PCR-based quantitative assays and cutting-edge third-generation DNA sequencing. A single "silver bullet" detection technology that will satisfy all study objectives does not exist; thus, every tool has a niche in an integrated detection and monitoring program. This review presents an overview of the rapidly changing landscape of fungicide resistance detection, illustrates how molecular techniques are being exploited to combat fungicide resistance in cereal crop phytopathogens, and highlights challenges and future research directions to aid in the design of effective monitoring systems that aim to apply fungicides strategically and minimize the cost of resistance. [Formula: see text] Copyright © 2025 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Repeat-induced point mutations driving Parastagonospora nodorum genomic diversity are balanced by selection against non-synonymous mutations

Parastagonospora nodorum is necrotrophic fungal pathogen of wheat with significant genomic resources. Population-level pangenome data for 173 isolates, of which 156 were from Western Australia (WA) and 17 were international, were examined for overall genomic diversity and effector gene content. A heterothallic core population occurred across all regions of WA, with asexually-reproducing clonal clusters in dryer northern regions. High potential for SNP diversity in the form of repeat-induced point mutation (RIP)-like transitions, was observed across the genome, suggesting widespread 'RIP-leakage' from transposon-rich repetitive sequences into non-repetitive regions. The strong potential for RIP-like mutations was balanced by negative selection against non-synonymous SNPs, that was observed within protein-coding regions. Protein isoform profiles of known effector loci (SnToxA, SnTox1, SnTox3, SnTox267, and SnTox5) indicated low-levels of non-synonymous and high-levels of silent RIP-like mutations. Effector predictions identified 186 candidate secreted predicted effector proteins (CSEPs), 69 of which had functional annotations and included confirmed effectors. Pangenome-based effector isoform profiles across WA were distinct from global isolates and were conserved relative to population structure, and may enable new approaches for monitoring crop disease pathotypes.

CbCyp51-Mediated Demethylation Inhibitor Resistance Is Modulated by Codon Bias

Cercospora leaf spot, caused by the fungus Cercospora beticola, is the most destructive foliar disease of sugarbeet worldwide. Resistance to the sterol demethylation inhibitor (DMI) fungicide tetraconazole has been previously correlated with synonymous and nonsynonymous mutations in CbCyp51. Here, we extend these analyses to the DMI fungicides prothioconazole, difenoconazole, and mefentrifluconazole in addition to tetraconazole to confirm whether the synonymous and nonsynonymous mutations at amino acid positions 144 and 170 are associated with resistance to these fungicides. Nearly half of the 593 isolates of C. beticola collected in the Red River Valley of North Dakota and Minnesota in 2021 were resistant to all four DMIs. Another 20% were resistant to tetraconazole and prothioconazole but sensitive to difenoconazole and mefentrifluconazole. A total of 13% of isolates were sensitive to all DMIs tested. We found five CbCyp51 haplotypes and associated them with phenotypes to the four DMIs. The most predominant haplotype (E170_A/L144F_C) correlated with resistance to all four DMIs with up to 97.6% accuracy. The second most common haplotype (E170_A/L144) consisted of isolates associated with resistance phenotypes to tetraconazole and prothioconazole while also exhibiting sensitive phenotypes to difenoconazole and mefentrifluconazole with up to 98.4% accuracy. Quantitative PCR did not identify differences in CbCyp51 expression between haplotypes. This study offers an understanding of the importance of codon usage in fungicide resistance and provides crop management acuity for fungicide application decision-making.

Deciphering fungicide resistance in Phytophthora: mechanisms, prevalence, and sustainable management approaches

The genus Phytophthora contains more than 100 plant pathogenic species that parasitize a wide range of plants, including economically important fruits, vegetables, cereals, and forest trees, causing significant losses. Global agriculture is seriously threatened by fungicide resistance in Phytophthora species, which makes it imperative to fully comprehend the mechanisms, frequency, and non-chemical management techniques related to resistance mutations. The mechanisms behind fungicide resistance, such as target-site mutations, efflux pump overexpression, overexpression of target genes and metabolic detoxification routes for fungicides routinely used against Phytophthora species, are thoroughly examined in this review. Additionally, it assesses the frequency of resistance mutations in various Phytophthora species and geographical areas, emphasizing the rise of strains that are resistant to multiple drugs. The effectiveness of non-chemical management techniques, including biological control, host resistance, integrated pest management plans, and cultural practices, in reducing fungicide resistance is also thoroughly evaluated. The study provides important insights for future research and the development of sustainable disease management strategies to counter fungicide resistance in Phytophthora species by synthesizing current information and identifying knowledge gaps.

Genetic diversity, mating type and pathogenicity of two Phytophthora species infecting black pepper in India

Phytophthora capsici and P. tropicalis are the two species of Phytophthora associated with foot rot disease of black pepper in India. High genetic diversity amongst the Phytophthora species contributes to its wide host range and variability in the virulence pattern. In the present study, genetic diversity of Phytophthora species infecting black pepper was analysed using RAMS (Random Amplified Microsatellites) and REP (Repetitive Extragenic Palindromic)-PCR fingerprinting. Forty-eight isolates, 24 each of P. capsici and P. tropicalis collected from major black pepper growing states, such as Karnataka, Kerala, Tamil Nadu and Goa, were used in the study. The analyses revealed a total of 160 loci of which 150 (93.75%) were polymorphic. UPGMA cluster and PCoA analysis based on combined RAMS and REP-PCR data clearly grouped the P. capsici and P. tropicalis isolates into two clusters which were further divided into four sub-clusters viz., I & II (P. capsici) and III & IV (P. tropicalis). The study clearly indicated that all the isolates were genetically unique and the entire population was heterogeneous. REP-PCR primers showed more polymorphic loci than RAMS primers. Further, sixteen isolates were selected for morphological and infectivity analyses under in vitro conditions. The isolates exhibited varied colony morphology, sporangial shapes and belonged to A1 mating type. Under in vitro conditions, all the sixteen black pepper Phytophthora isolates could infect nutmeg, tomato, chilli, pumpkin, and cucumber and few of the isolates could infect cardamom. None of the isolates could infect coconut, areca nut and vanilla.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03843-1.

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.

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.

Important wheat diseases in the US and their management in the 21st century

Wheat is a crop of historical significance, as it marks the turning point of human civilization 10,000 years ago with its domestication. Due to the rapid increase in population, wheat production needs to be increased by 50% by 2050 and this growth will be mainly based on yield increases, as there is strong competition for scarce productive arable land from other sectors. This increasing demand can be further achieved using sustainable approaches including integrated disease pest management, adaption to warmer climates, less use of water resources and increased frequency of abiotic stress tolerances. Out of 200 diseases of wheat, 50 cause economic losses and are widely distributed. Each year, about 20% of wheat is lost due to diseases. Some major wheat diseases are rusts, smut, tan spot, spot blotch, fusarium head blight, common root rot, septoria blotch, powdery mildew, blast, and several viral, nematode, and bacterial diseases. These diseases badly impact the yield and cause mortality of the plants. This review focuses on important diseases of the wheat present in the United States, with comprehensive information of causal organism, economic damage, symptoms and host range, favorable conditions, and disease management strategies. Furthermore, major genetic and breeding efforts to control and manage these diseases are discussed. A detailed description of all the QTLs, genes reported and cloned for these diseases are provided in this review. This study will be of utmost importance to wheat breeding programs throughout the world to breed for resistance under changing environmental conditions.

Genome-Wide Identification of Variants Associated with Antifungal Drug Resistance

Genomic approaches are widely applied to study the genetic basis of antifungal drug resistance in clinical isolates and experimental studies. Whole-genome sequencing of clinical isolates can comprehensively identify mutations associated with drug resistance and their frequency across fungal populations. In addition, genome comparison of serially collected isolates, such as from patient samples or in vitro drug selection experiments, will identify a small number of changes that can be evaluated for association with drug resistance. Here, we provide a detailed protocol for the computational analysis of genome sequences to identify variants associated with drug resistance.

New insights from short and long reads sequencing to explore cytochrome b variants in Plasmopara viticola populations collected from vineyards and related to resistance to complex III inhibitors

Downy mildew is caused by Plasmopara viticola, an obligate oomycete plant pathogen, a devasting disease of grapevine. To protect plants from the disease, complex III inhibitors are among the fungicides widely used. They specifically target the mitochondrial cytochrome b (cytb) of the pathogen to block cellular respiration mechanisms. In the French vineyard, P. viticola has developed resistance against a first group of these fungicides, the Quinone outside Inhibitors (QoI), with a single amino acid substitution G143A in its cytb mitochondrial sequence. The use of QoI was limited and another type of fungicide, the Quinone inside Inhibitors, targeting the same gene and highly effective against oomycetes, was used instead. Recently however, less sensitive P. viticola populations were detected after treatments with some inhibitors, in particular ametoctradin and cyazofamid. By isolating single-sporangia P. viticola strains resistant to these fungicides, we characterized new variants in the cytb sequences associated with cyazofamid resistance: a point mutation (L201S) and more strikingly, two insertions (E203-DE-V204, E203-VE-V204). In parallel with the classical tools, pyrosequencing and qPCR, we then benchmarked short and long-reads NGS technologies (Ion Torrent, Illumina, Oxford Nanopore Technologies) to sequence the complete cytb with a view to detecting and assessing the proportion of resistant variants of P. viticola at the scale of a field population. Eighteen populations collected from French vineyard fields in 2020 were analysed: 12 showed a variable proportion of G143A, 11 of E203-DE-V204 and 7 populations of the S34L variant that confers resistance to ametoctradin. Interestingly, the long reads were able to identify variants, including SNPs, with confidence and to detect a small proportion of P. viticola with multiple variants along the same cytb sequence. Overall, NGS appears to be a promising method for assessing fungicide resistance of pathogens linked to cytb modifications at the field population level. This approach could rapidly become a robust decision support tool for resistance management in the future.

Global genomic analyses of wheat powdery mildew reveal association of pathogen spread with historical human migration and trade

The fungus Blumeria graminis f. sp. tritici causes wheat powdery mildew disease. Here, we study its spread and evolution by analyzing a global sample of 172 mildew genomes. Our analyses show that B.g. tritici emerged in the Fertile Crescent during wheat domestication. After it spread throughout Eurasia, colonization brought it to America, where it hybridized with unknown grass mildew species. Recent trade brought USA strains to Japan, and European strains to China. In both places, they hybridized with local ancestral strains. Thus, although mildew spreads by wind regionally, our results indicate that humans drove its global spread throughout history and that mildew rapidly evolved through hybridization.

A Secondary Metabolite Secreted by Penicillium citrinum Is Able to Enhance Parastagonospora nodorum Sensitivity to Tebuconazole and Azoxystrobin

Parastagonospora nodorum causes glume and leaf blotch of wheat, a harmful disease resulting in serious losses in grain yield. In many countries including Russia, fungicidal formulations based on triazoles and on triazoles combined with strobilurins are used to control this fungus. However, their prolonged application may promote the selection of fungicide-resistant strains of P. nodorum leading to significant attenuation or even loss of fungicidal effect. Chemosensitization of plant pathogenic fungi with natural compounds represents a promising strategy for mitigating fungicide resistance and other negative impacts of fungicides. In this work, we applied a chemosensitization approach towards P. nodorum strains non-resistant or resistant to tebuconazole or azoxystrobin using 6-demethylmevinolin (6-DMM), a metabolite of Penicillium citrinum. The resistant strains were obtained by the mutagenesis and subsequent culturing on agar media incorporated with increasing doses of Folicur® EC 250 (i.e., tebuconazole) or Quadris® SC 250 (i.e., azoxystrobin). Test strains m8-4 and kd-18, most resistant to tebuconazole and azoxystrobin, respectively, were selected for sensitization experiments. These experiments demonstrated that combining 6-DMM with Folicur® enhanced fungicidal effectiveness in vitro and in vivo in addition to attenuating the resistance of P. nodorum to tebuconazole in vitro. 6-DMM was also found to augment Quadris® efficacy towards kd-18 when applied on detached wheat leaves inoculated with this strain. Experiments on P. nodorum sensitization under greenhouse conditions included preventive (applying test compounds simultaneously with inoculation) or post-inoculation spraying of wheat seedlings with 6-DMM together with Folicur® at dose rates (DR) amounting to 10% and 20% of DR recommended for field application (RDR). Combined treatments were run in parallel with using the same DR of the fungicide and sensitizer, alone. A synergistic effect was observed in both preventive and post-inoculation treatments, when the sensitizer was co-applied with the fungicide at 10% of the RDR. In this case, disease reduction significantly exceeded the protective effect of Folicur® at 10% or 20% of the RDR, alone, and also a calculated additive effect. Collectively, our findings suggest that 6-DMM is promising as a putative component for formulations with triazole and strobilurin fungicides. Such new formulations would improve fungicide efficacy and, potentially, lower rates of fungicides needed for plant pathogen control.

Biology and molecular interactions of Parastagonospora nodorum blotch of wheat

Parastagonospora nodorum is one of the important necrotrophic pathogens of wheat which causes severe economical loss to crop yield. So far, a number of effectors of Parastagonospora nodorum origin and their target interacting genes on the host plant have been characterized. Since targeting effector-sensitive gene carefully can be helpful in breeding for resistance. Therefore, constant efforts are required to further characterize the effectors, their interacting genes, and underlying biochemical pathways. Furthermore, to develop effective counter-strategies against emerging diseases, continuous efforts are required to determine the qualitative resistance that demands to screen of diverse genotypes for host resistance. Stagonospora nodorum blotch also refers to as Stagonospora glume blotch and leaf is caused by Parastagonospora nodorum. The pathogen deploys necrotrophic effectors for the establishment and development on wheat plants. The necrotrophic effectors and their interaction with host receptors lead to the establishment of infection on leaves and extensive lesions formation which either results in host cell death or suppression/activation of host defence mechanisms. The wheat Stagonospora nodorum interaction involves a set of nine host gene-necrotrophic effector interactions. Out of these, Snn1-SnTox1, Tsn1-SnToxA and Snn-SnTox3 are one of the most studied interaction, due to its role in the suppression of reactive oxygen species production, regulating the cytokinin content through ethylene-dependent wayduring initial infection stage. Further, although the molecular basis is not fully unveiled, these effectors regulate the redox state and influence the ethylene biosynthesis in infected wheat plants. Here, we have discussed the biology of the wheat pathogen Parastagonospora nodorum, role of its necrotrophic effectors and their interacting sensitivity genes on the redox state, how they hijack the resistance mechanisms, hormonal regulated immunity and other signalling pathways in susceptible wheat plants. The information generated from effectors and their corresponding sensitivity genes and other biological processes could be utilized effectively for disease management strategies.

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.

Population genomics of transposable element activation in the highly repressive genome of an agricultural pathogen

The activity of transposable elements (TEs) can be an important driver of genetic diversity with TE-mediated mutations having a wide range of fitness consequences. To avoid deleterious effects of TE activity, some fungi have evolved highly sophisticated genomic defences to reduce TE proliferation across the genome. Repeat-induced point mutation (RIP) is a fungal-specific TE defence mechanism efficiently targeting duplicated sequences. The rapid accumulation of RIPs is expected to deactivate TEs over the course of a few generations. The evolutionary dynamics of TEs at the population level in a species with highly repressive genome defences is poorly understood. Here, we analyse 366 whole-genome sequences of Parastagonospora nodorum, a fungal pathogen of wheat with efficient RIP. A global population genomics analysis revealed high levels of genetic diversity and signs of frequent sexual recombination. Contrary to expectations for a species with RIP, we identified recent TE activity in multiple populations. The TE composition and copy numbers showed little divergence among global populations regardless of the demographic history. Miniature inverted-repeat transposable elements (MITEs) and terminal repeat retrotransposons in miniature (TRIMs) were largely underlying recent intra-species TE expansions. We inferred RIP footprints in individual TE families and found that recently active, high-copy TEs have possibly evaded genomic defences. We find no evidence that recent positive selection acted on TE-mediated mutations rather that purifying selection maintained new TE insertions at low insertion frequencies in populations. Our findings highlight the complex evolutionary equilibria established by the joint action of TE activity, selection and genomic repression.

Genome-scale phylogenies reveal relationships among Parastagonospora species infecting domesticated and wild grasses

Several plant pathogenic Parastagonospora species have been identified infecting wheat and other cereals over the past 50 years. As new lineages were discovered, naming conventions grew unwieldy and the relationships with previously recognized species remained unclear. We used genome sequencing to clarify relationships among these species and provided new names for most of these species. Six of the nine described Parastagonospora species were recovered from wheat, with five of these species coming from Iran. Genome sequences revealed that three strains thought to be hybrids between P. nodorum and P. pseudonodorum were not actually hybrids, but rather represented rare gene introgressions between those species. Our data are consistent with the hypothesis that P. nodorum originated as a pathogen of wild grasses in the Fertile Crescent, then emerged as a wheat pathogen via host-tracking during the domestication of wheat in the same region. The discovery of a diverse array of Parastagonospora species infecting wheat in Iran suggests that new wheat pathogens could emerge from this region in the future. Citation: Croll D, Crous PW, Pereira D, et al. 2021. Genome-scale phylogenies reveal relationships among Parastagonospora species infecting domesticated and wild grasses. Persoonia 46: 116-128. https://doi.org/10.3767/persoonia.2021.46.04.

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.