Fungicide efflux and the MgMFS1 transporter contribute to the multidrug resistance phenotype inZymoseptoria triticifield isolates

The Synergy of Chitosan and Azoxystrobin Against Fusarium graminearum Is Modulated by Selected ABC Transporters

The development of innovative and effective strategies to combat fungal pathogens is critical to sustainable crop protection. Fungicides have been used for over two centuries, with traditional copper- and sulfur-based formulations still in use due to their broad-spectrum, multisite mode of action, which minimizes the risk of pathogen resistance. In contrast, modern systemic fungicides, though potent, often target a single site of action, leading to the accelerated emergence of resistant fungal strains. This study explores synergistic interactions between chitosan (CS) and selected fungicides, focusing on their antifungal activity against Fusarium graminearum. Among the fungicides tested, azoxystrobin (Amistar) exhibited the highest 44.88 synergy score when combined with CS (30 kDa, degree of deacetylation ≥ 90), resulting in significantly improved antifungal efficacy. Furthermore, the combination of CS and Amistar with double-stranded RNA (dsRNA) targeting selected ABC transporter genes further amplified antifungal activity by silencing genes critical for fungal tolerance to treatment. This dual synergy highlights the potential of RNA interference (RNAi) as both a functional tool to investigate fungal physiology and an effective antifungal strategy. These findings reveal a promising and environmentally friendly approach to mitigate resistance while improving fungal control. Furthermore, the remarkable synergy between azoxystrobin and CS presents a novel mechanism with significant potential for sustainable agricultural applications, which warrants further investigation to elucidate its molecular basis.

Sterol regulatory element-binding proteins mediate intrinsic fungicide tolerance and antagonism in the fungal biocontrol agent Clonostachys rosea IK726

Sterol regulatory element-binding proteins (SREBPs) are transcription factors governing various biological processes in fungi, including virulence and fungicide tolerance, by regulating ergosterol biosynthesis and homeostasis. While studied in model fungal species, their role in fungal species used for biocontrol remains elusive. This study delves into the biological and regulatory function of SREBPs in the fungal biocontrol agent (BCA) Clonostachys rosea IK726, with a specific focus on fungicide tolerance and antagonism. Clonostachys rosea genome contains two SREBP coding genes (sre1 and sre2) with distinct characteristics. Deletion of sre1 resulted in mutant strains with pleiotropic phenotypes, including reduced C. rosea growth on medium supplemented with prothioconazole and boscalid fungicides, hypoxia mimicking agent CoCl2 and cell wall stressor SDS, and altered antagonistic abilities against Botrytis cinerea and Rhizoctonia solani. However, Δsre2 strains showed no significant effect. Consistent with the gene deletion results, overexpression of sre1 in Saccharomyces cerevisiae enhanced tolerance to prothioconazole. The functional differentiation between SRE1 and SRE2 was elucidated by the yeast-two-hybridization assay, which showed an interaction between SREBP cleavage-activating protein (SCAP) and SRE1 but not between SRE2 and SCAP. Transcriptome analysis of the Δsre1 strain unveiled SRE1-mediated expression regulation of genes involved in lipid metabolism, respiration, and xenobiotic tolerance. Notably, genes coding for antimicrobial compounds chitinases and polyketide synthases were downregulated, aligning with the altered antagonism phenotype. This study uncovers the role of SREBPs in fungal BCAs, providing insights for C. rosea IK726 application into integrated pest management strategies.

Overexpression of the CcCYP51A and CcCYP51B genes confer Corynespora cassiicola resistance to prochloraz

Cucumber Corynespora leaf spot caused by Corynespora cassiicola is the primary disease responsible for reducing cucumber yield, and prochloraz is the main fungicide used to control C. cassiicola. This study investigated the sensitivity and resistance mechanism of C. cassiicola isolates to prochloraz, and found that C. cassiicola has developed resistance to prochloraz. The prochloraz EC50 values ranged from 0.02 to 2.33 μg/mL, with a mean of 0.436 ± 0.447 μg/mL. In total, 36 of 146 isolates exhibited prochloraz resistance. The resistant isolates had no fitness cost and could not be completely controlled by 50 μg/mL prochloraz on detached leaves. Prochloraz exhibited positive cross-resistance with propiconazole and tebuconazole but not with difenoconazole, carbendazim, trifloxystrobin and pydiflumetofen. The sensitive isolates had significantly lower ergosterol content than the resistant isolates after prochloraz treatment. Compared to sensitive isolates, prochloraz-resistant isolates had no CcCYP51 gene mutation, but the CcCYP51A and CcCYP51B gene expression levels were significantly higher under the treatment of prochloraz. The overexpression of CcCYP51A and CcCYP51B were associated with prochloraz resistance in C. cassiicola.

Mutations across Diverse Domains of CjXDR1 Lead to Multidrug Resistance in Clarireedia jacksonii

Recently, Clarireedia jacksonii has emerged as a significant pathogen threatening turfgrass, and its escalating resistance to multiple drugs often undermines field interventions. This study highlighted the critical role of the fungus-specific transcription factor CjXDR1 (formerly ShXDR1) in regulating multidrug resistance (MDR) in C. jacksonii. This was demonstrated through experiments involving CjXDR1-knockout and CjXDR1-complemented strains. Our sequence analysis revealed five mutations in CjXDR1: G445D, K453E, S607F, D676H, and V690A. All five gain-of-function (GOF) mutations were confirmed to directly contribute to MDR against three different classes of fungicides (propiconazole: demethylation inhibitor, boscalid: succinate dehydrogenase inhibitor, and iprodione: dicarboximide) using the genetic transformation system and in vitro fungicide-sensitivity assay. Comparative transcriptome analysis revealed that CjXDR1 and its GOF mutations led to the overexpression of downstream genes encoding a Phase I metabolizing enzyme (CYP68) and two Phase III transporters (CjPDR1 and CjAtrD) previously reported. Knockout mutants of CYP68, CjPDR1, CjAtrD, and double-knockout mutants of CjPDR1 and CjAtrD exhibited increased sensitivity to all three fungicides tested. Among these, the CYP68-knockout mutants displayed the highest sensitivity to propiconazole, while the CjPDR1 knockout mutant exhibited significantly increased sensitivity to all three fungicides. Double-knockout mutants of CjPDR1 and CjAtrD displayed greater sensitivity than the single knockouts. In conclusion, multiple GOF mutants in CjXDR1 contribute to MDR by upregulating the expression of CjPDR1, CjAtrD, and CYP68. This study enhances our understanding of the molecular mechanisms underlying MDR in plant pathogenic fungi, providing valuable insights into GOF mutation structures and advancing the development of antifungal drugs.

Deciphering the mechanisms involved in reduced sensitivity to azoles and fengycin lipopeptide in Venturia inaequalis

Apple scab, caused by the hemibiotrophic fungus Venturia inaequalis, is currently the most common and damaging disease in apple orchards. Two strains of V. inaequalis (S755 and Rs552) with different sensitivities to azole fungicides and the bacterial metabolite fengycin were compared to determine the mechanisms responsible for these differences. Antifungal activity tests showed that Rs552 had reduced sensitivity to tebuconazole and tetraconazole, as well as to fengycin alone or in a binary mixture with other lipopeptides (iturin A, pumilacidin, lichenysin). S755 was highly sensitive to fengycin, whose activity was close to that of tebuconazole. Unlike fengycin, lipopeptides from the iturin family (mycosubtilin, iturin A) had similar activity on both strains, while those from the surfactin family (lichenysin, pumilacidin) were not active, except in binary mixtures with fengycin. The activity of lipopeptides varies according to their family and structure. Analyses to determine the difference in sensitivity to azoles (which target the CYP51 enzyme involved in the ergosterol biosynthesis pathway) showed that the reduced sensitivity in Rs552 is linked to (i) a constitutive increased expression of the Cyp51A gene caused by insertions in the upstream region and (ii) greater efflux by membrane pumps with the involvement of ABC transporters. Microscopic observations revealed that fengycin, known to interact with plasma membranes, induced morphological and cytological changes in cells from both strains. Sterol and phospholipid analyses showed a higher level of ergosta-7,22-dien-3-ol and a lower level of PI(C16:0/C18:1) in Rs552 compared with S755. These differences could therefore influence the composition of the plasma membrane and explain the differential sensitivity of the strains to fengycin. However, the similar antifungal activities of mycosubtilin and iturin A in the two strains indirectly indicate that sterols are probably not involved in the fengycin resistance mechanism. This leads to the conclusion that different mechanisms are responsible for the difference in susceptibility to azoles or fengycin in the strains studied.

Combination of Bacillus tequilensis with difenoconazole to control pear black spot and the related synergistic mechanism

Background

Pear black spot (PBS) is caused by Alternaria alternata and causes severe damage worldwide. It is particularly important to screen for synergistic fungicide combinations to address issues associated with the low efficacy of biocontrol agents, high dosage requirements and poor sustained effectiveness of chemical fungicides.

Methods

In vitro and in vivo studies were performed to determine the efficacy of a treatment for this important disease. Additionally, transcriptomic and metabolomic analyses were performed to determine the main molecular and biochemical mechanisms involved in the interaction.

Results

Bacillus tequilensis 2_2a has a significant synergistic effect with difenoconazole, causing hyphal entanglement and spore lysis and inhibiting the formation of PBS lesions in vitro. In the field, the control effect of the combination was greater than 95%. The pathways associated with the synergistic effect on the mycelia of A. alternata were divided into two main types: one included glycolysis, oxidative phosphorylation, and MAPK signal transduction, while the other included glycolysis, the TCA cycle, coenzyme A biosynthesis, sterol synthesis, and fatty acid degradation. Both types of pathways jointly affect the cell cycle. The main functions of the key genes and metabolites that have been verified as being affected are glucose synthesis and oxidative respiration, as well as citric acid synthesis, acetyl-CoA synthesis, and sterol synthesis. Both functions involve intracellular pyridine nucleotide metabolism and adenine nucleotide transformation.

Conclusion

This study helps to reveal the synergistic mechanisms underlying the combined efficacy of biological and chemical agents, providing a scientific basis for field applications.

Transcriptome Analysis Reveals Potential Regulators of DMI Fungicide Resistance in the Citrus Postharvest Pathogen Penicillium digitatum

Green mold, caused by Penicillium digitatum, is the major cause of citrus postharvest decay. Currently, the application of sterol demethylation inhibitor (DMI) fungicide is one of the main control measures to prevent green mold. However, the fungicide-resistance problem in the pathogen P. digitatum is growing. The regulatory mechanism of DMI fungicide resistance in P. digitatum is poorly understood. Here, we first performed transcriptomic analysis of the P. digitatum strain Pdw03 treated with imazalil (IMZ) for 2 and 12 h. A total of 1338 genes were up-regulated and 1635 were down-regulated under IMZ treatment for 2 h compared to control while 1700 were up-regulated and 1661 down-regulated under IMZ treatment for 12 h. The expression of about half of the genes in the ergosterol biosynthesis pathway was affected during IMZ stress. Further analysis identified that 84 of 320 transcription factors (TFs) were differentially expressed at both conditions, making them potential regulators in DMI resistance. To confirm their roles, three differentially expressed TFs were selected to generate disruption mutants using the CRISPR/Cas9 technology. The results showed that two of them had no response to IMZ stress while ∆PdflbC was more sensitive compared with the wild type. However, disruption of PdflbC did not affect the ergosterol content. The defect in IMZ sensitivity of ∆PdflbC was restored by genetic complementation of the mutant with a functional copy of PdflbC. Taken together, our results offer a rich source of information to identify novel regulators in DMI resistance.

Preventing multi-resistance: New insights for managing fungal adaptation

Sustainable crop protection is vital for food security, yet it is under threat due to the adaptation of a diverse and evolving pathogen population. Resistance can be managed by maximising the diversity of selection pressure through dose variation and the spatial and temporal combination of active ingredients. This study explores the interplay between operational drivers for maximising the sustainability of management strategies in relation to the resistance status of fungal populations. We applied an experimental evolution approach to three artificial populations of Zymoseptoria tritici, an economically significant wheat pathogen, each differing in initial resistance status. Our findings reveal that diversified selection pressure curtails the selection of resistance in naïve populations and those with low frequencies of single resistance. Increasing the number of modes of action most effectively delays resistance development, surpassing the increase in the number of fungicides, fungicide choice based on resistance risk, and temporal variation in fungicide exposure. However, this approach favours generalism in the evolved populations. The prior presence of multiple resistant isolates and their subsequent selection in populations override the effects of diversity in management strategies, thereby invalidating any universal ranking. Therefore, the initial resistance composition must be specifically considered in sustainable resistance management to address real-world field situations.

A systematic screen for co-option of transposable elements across the fungal kingdom

How novel protein functions are acquired is a central question in molecular biology. Key paths to novelty include gene duplications, recombination or horizontal acquisition. Transposable elements (TEs) are increasingly recognized as a major source of novel domain-encoding sequences. However, the impact of TE coding sequences on the evolution of the proteome remains understudied. Here, we analyzed 1237 genomes spanning the phylogenetic breadth of the fungal kingdom. We scanned proteomes for evidence of co-occurrence of TE-derived domains along with other conventional protein functional domains. We detected more than 13,000 predicted proteins containing potentially TE-derived domain, of which 825 were identified in more than five genomes, indicating that many host-TE fusions may have persisted over long evolutionary time scales. We used the phylogenetic context to identify the origin and retention of individual TE-derived domains. The most common TE-derived domains are helicases derived from Academ, Kolobok or Helitron. We found putative TE co-options at a higher rate in genomes of the Saccharomycotina, providing an unexpected source of protein novelty in these generally TE depleted genomes. We investigated in detail a candidate host-TE fusion with a heterochromatic transcriptional silencing function that may play a role in TE and gene regulation in ascomycetes. The affected gene underwent multiple full or partial losses within the phylum. Overall, our work establishes a kingdom-wide view of putative host-TE fusions and facilitates systematic investigations of candidate fusion proteins.

Quantitative trait locus mapping of osmotic stress response in the fungal wheat pathogen Zymoseptoria tritici

Osmotic stress is a ubiquitous and potent stress for all living organisms, but few studies have investigated the genetic basis of salt tolerance in filamentous fungi. The main aim of this study was to identify regions of the genome associated with tolerance to potassium chloride (KCl) in the wheat pathogen Zymoseptoria tritici. A secondary aim was to identify candidate genes affecting salt tolerance within the most promising chromosomal regions. We achieved these aims with a quantitative trait locus (QTL) mapping study using offspring from 2 crosses grown in vitro in the presence or absence of osmotic stress imposed by 0.75 M KCl. We identified significant QTL for most of the traits in both crosses. Several QTLs overlapped with QTL identified in earlier studies for other traits, and some QTL explained trait variation in both the control and salt stress environments. A significant QTL on chromosome 3 explained variation in colony radius at 8-day postinoculation (dpi) in the KCl environment as well as colony radius KCl tolerance at 8 dpi. The QTL peak had a high logarithm of the odds ratio (LOD) and encompassed an interval containing only 36 genes. Six of these genes present promising candidates for functional analyses. A gene ontology (GO) enrichment analysis of QTL unique to the KCl environment found evidence for the enrichment of functions involved in osmotic stress responses.

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.

Dose-Dependent Genetic Resistance to Azole Fungicides Found in the Apple Scab Pathogen

The evolution of azole resistance in fungal pathogens presents a major challenge in both crop production and human health. Apple orchards across the world are faced with the emergence of azole fungicide resistance in the apple scab pathogen Venturia inaequalis. Target site point mutations observed in this fungus to date cannot fully explain the reduction in sensitivity to azole fungicides. Here, polygenic resistance to tebuconazole was studied across a population of V. inaequalis. Genotyping by sequencing allowed Quantitative Trait Loci (QTLs) mapping to identify the genetic components controlling this fungicide resistance. Dose-dependent genetic resistance was identified, with distinct genetic components contributing to fungicide resistance at different exposure levels. A QTL within linkage group seven explained 65% of the variation in the effective dose required to reduce growth by 50% (ED50). This locus was also involved in resistance at lower fungicide doses (ED10). A second QTL in linkage group one was associated with dose-dependent resistance, explaining 34% of variation at low fungicide doses (ED10), but did not contribute to resistance at higher doses (ED50 and ED90). Within QTL regions, non-synonymous mutations were observed in several ATP-Binding Cassette and Major Facilitator SuperFamily transporter genes. These findings provide insight into the mechanisms of fungicide resistance that have evolved in horticultural pathogens. Identification of resistance gene candidates supports the development of molecular diagnostics to inform management practices.

Multiple and multidrug resistance in Botrytis cinerea: molecular mechanisms of MLR/MDR strains in Greece and effects of co-existence of different resistance mechanisms on fungicide sensitivity

Botrytis cinerea is a high-risk pathogen for fungicide resistance development. Within the fungal populations, strains have developed multiple mutations in different target genes leading to multiple resistance (MLR) or mutations associated with overexpression of efflux transporters leading to multidrug resistance (MDR). These types of resistance are a major threat, and their successful management is a major challenge. The current study was initiated to a) determine frequencies of MLR/MDR strains in populations originating from several crops, b) identify the types of MDR that occur in Greece, and c) determine interactions between MLR and MDR at the level of sensitivity to botryticides. The frequencies of MLR/MDR phenotypes were determined in 515 isolates subjected to bioassays using discriminatory concentrations of thiophanate-methyl, iprodione, cyprodinil, fenhexamid, boscalid, fluopyram, fludioxonil, pyraclostrobin, and tolnaftate. Interestingly, 7.8% and 31.3% of isolates from strawberry and rootstock seedlings were resistant to every single fungicide class, while MDR phenotypes from strawberries, rootstocks, and tomatoes accounted for 26%, 87%, and 13.4%, respectively. The MLR and MDR isolates were further molecularly analyzed regarding genes erg27, sdhB, Bcpos5, and Mrr1, responsible for resistance to fenhexamid, boscalid and fluopyram, cyprodinil, and MDR, respectively. The different mutations' presence was determined along with a new mutation in Mrr1 leading to MDR. MDR isolates were characterized as MDR1 or MDR1h based on the presence of a 3-bp deletion in Mrr1. MDR1h was predominant in isolates from rootstocks and MDR1 from tomatoes and strawberries, whereas the most frequent target-site mutations were F412S (erg27), H272R (sdhB), and L412F (Bcpos5). To determine whether the accumulation of target-site mutations along with MDR mutations exhibits an additive effect concerning fungicide resistance, the sensitivity of isolates possessing the predominant target-site mutations was calculated in both the presence and the absence of MDR-associated mutations. EC50 in cyprodinil and boscalid increased to about twofold in the presence of MDR mutations, while there was no difference for fenhexamid. In conclusion, MLR/MDR frequencies are notably high in heavily treated crops in Greece, and the combination of MLR and MDR mutations leads to even higher fungicide resistance levels, highlighting the importance of resistance management.

Deep amplicon sequencing reveals extensive allelic diversity in the erg11/CYP51 promoter and allows multi-population DMI fungicide resistance monitoring in the canola pathogen Leptosphaeria maculans

Continued use of fungicides provides a strong selection pressure towards strains with mutations to render these chemicals less effective. Previous research has shown that resistance to the demethylation inhibitor (DMI) fungicides, which target ergosterol synthesis, in the canola pathogen Leptosphaeria maculans has emerged in Australia and Europe. The change in fungicide sensitivity of individual isolates was found to be due to DNA insertions into the promoter of the erg11/CYP51 DMI target gene. Whether or not these were the only types of mutations and how prevalent they were in Australian populations was explored in the current study. New isolates with reduced DMI sensitivity were obtained from screens on DMI-treated plants, revealing eight independent insertions in the erg11 promoter. A novel deep amplicon sequencing approach applied to populations of ascospores fired from stubble identified an additional undetected insertion allele and quantified the frequencies of all known insertions, suggesting that, at least in the samples processed, the combined frequency of resistant alleles is between 0.0376% and 32.6%. Combined insertion allele frequencies positively correlated with population-level measures of in planta resistance to four different DMI treatments. Additionally, there was no evidence for erg11 coding mutations playing a role in conferring resistance in Australian populations. This research provides a key method for assessing fungicide resistance frequency in stubble-borne populations of plant pathogens and a baseline from which additional surveillance can be conducted in L. maculans. Whether or not the observed resistance allele frequencies are associated with loss of effective disease control in the field remains to be established.

Genome-wide transcriptional analyses of Clarireedia jacksonii isolates associated with multi-drug resistance

Emerging multidrug resistance (MDR) in Clarireedia spp. is a huge challenge to the management of dollar spot (DS) disease on turfgrass. Insight into the molecular basis of resistance mechanisms may help identify key molecular targets for developing novel effective chemicals. Previously, a MDR isolate (LT586) of C. jacksonii with significantly reduced sensitivities to propiconazole, boscalid, and iprodione, and a fungicide-sensitive isolate (LT15) of the same species were isolated from creeping bentgrass (Agrostis stolonifera L.). The present study aimed to further explore the molecular mechanisms of resistance by using genome-wide transcriptional analyses of the two isolates. A total of 619 and 475 differentially expressed genes (DEGs) were significantly down and upregulated in the MDR isolate LT586, compared with the sensitive isolate LT15 without fungicide treatment. Three hundreds and six and 153 DEGs showed significantly lower and higher expression in the MDR isolate LT586 than those in the sensitive isolate LT15, which were commonly induced by the three fungicides. Most of the 153 upregulated DEGs were xenobiotic detoxification-related genes and genes with transcriptional functions. Fifty and 17 upregulated DEGs were also commonly observed in HRI11 (a MDR isolate of the C. jacksonii) compared with the HRS10 (a fungicide-sensitive isolate of same species) from a previous study without and with the treatment of propiconazole, respectively. The reliability of RNA-seq data was further verified by qRT-PCR method using a few select potentially MDR-related genes. Results of this study indicated that there were multiple uncharacterized genes, possibly responsible for MDR phenotypes in Clarireedia spp., which may have important implications in understanding the molecular mechanisms underlying MDR resistance.

The Dynamic Fungal Genome: Polyploidy, Aneuploidy and Copy Number Variation in Response to Stress

Fungal species have dynamic genomes and often exhibit genomic plasticity in response to stress. This genome plasticity often comes with phenotypic consequences that affect fitness and resistance to stress. Fungal pathogens exhibit genome plasticity in both clinical and agricultural settings and often during adaptation to antifungal drugs, posing significant challenges to human health. Therefore, it is important to understand the rates, mechanisms, and impact of large genomic changes. This review addresses the prevalence of polyploidy, aneuploidy, and copy number variation across diverse fungal species, with special attention to prominent fungal pathogens and model species. We also explore the relationship between environmental stress and rates of genomic changes and highlight the mechanisms underlying genotypic and phenotypic changes. A comprehensive understanding of these dynamic fungal genomes is needed to identify novel solutions for the increase in antifungal drug resistance.

Multidrug resistance of Rhizoctonia solani determined by enhanced efflux for fungicides

Plant pathogens can develop multidrug resistance (MDR) through metabolomic and efflux activities. Although MDR has been observed in the field, its mechanisms are yet to be further studied. MDR in Rhizoctonia solani induced by the uncoupler SYP-14288, which involved efflux transporters including ATP binding cassette (ABC) and major facilitator superfamily (MFS) have been reported in our previous study. To confirm this, corresponding genes of the wild-type R. solani X19 and its derived MDR mutant X19-7 were compared through transcriptomics, RNA-Seq data validation, and heterologous expression. Genes encoding six ABC transporters and seven MFS transporters were identified to be associated with MDR and mostly showed a constitutively higher expression in X19-7 than in X19 regardless of SYP-14288 treatment. Eight ABC transporter-encoding genes and eight MFS transporter-encoding genes were further characterized by transferring into Saccharomyces cerevisiae. The sensitivity of transformants containing either ABC transporter-encoding gene AG1IA_06082 and MFS transporter-encoding gene AG1IA_08645 was significantly decreased in responses to fungicides having various modes of action including SYP-14288, fluazinam, chlorothalonil, and difenoconazole, indicating that these two genes were related to MDR. The roles of two genes were further confirmed by successfully detecting their protein products and high accumulation of SYP-14288 in yeast transformants. Thus, ABC and MFS transporters contributed to the development of MDR in R. solani. The result helps to understand the cause and mechanisms that influence the efficient use of fungicide.

Characterization and expression of fungal defensin in Escherichia coli and its antifungal mechanism by RNA-seq analysis

Introduction

Invasive fungal infections (IFIs) are fatally threatening to critical patients. The fungal defensin as an antifungal protein can widely inhibit fungi.

Methods

In this study, eight antifungal genes from different filamentous fungi were optimized by synonymous codon bias and heterologously expressed in Escherichia coli.

Results and discussion

Only the antifungal protein (AFP) from Aspergillus giganteus was produced, whereas the AFP from its mutation of the chitin-binding domain could not be expressed, thereby suggesting the importance of the motif for protein folding. In addition, the recombinant AFP (rAFP, 100 μg/mL) pre-heated at 50°C for 1 h effectively inhibited Paecilomyces variotii CICC40716 of IFIs by 55%, and no cell cytotoxicity was observed in RAW264.7 cells. After being pre-heated at 50°C for 8 h, the fluorescence emission intensity of the rAFP decreased and shifted from 343 nm to 335 nm. Moreover, the helix and β-turn of the rAFP gradually decreased with the pre-heated treatment temperature of 50°C via circular dichroism spectroscopy. Propidium iodide staining revealed that the rAFP could cause damage to the cell membrane. Moreover, the corresponding differentially expressed genes (DEGs) for downregulation such as amino sugar and nucleotide sugar metabolism, as well as mitogen-activated protein kinase (MAPK) signaling pathway involved in the cell wall integrity were found via the RNA-seq of rAFP treatment. By contrast, the upregulated DEGs were enriched in response to the oxidative stress of Biological Process by the Gene Ontology (GO) database. The encoding proteins of laccase, multicopper oxidase, and nitroreductase that contributed to reactive oxygen species (ROS) scavenging could be recognized. These results suggested that the rAFP may affect the integrity of the cell wall and cell membrane, and promote the increase in ROS, thereby resulting in fungal death. Consequently, drug development could be based on the inhibitory effect of the rAFP on IFIs.

A secreted protease-like protein in Zymoseptoria tritici is responsible for avirulence on Stb9 resistance gene in wheat

Zymoseptoria tritici is the fungal pathogen responsible for Septoria tritici blotch on wheat. Disease outcome in this pathosystem is partly determined by isolate-specific resistance, where wheat resistance genes recognize specific fungal factors triggering an immune response. Despite the large number of known wheat resistance genes, fungal molecular determinants involved in such cultivar-specific resistance remain largely unknown. We identified the avirulence factor AvrStb9 using association mapping and functional validation approaches. Pathotyping AvrStb9 transgenic strains on Stb9 cultivars, near isogenic lines and wheat mapping populations, showed that AvrStb9 interacts with Stb9 resistance gene, triggering an immune response. AvrStb9 encodes an unusually large avirulence gene with a predicted secretion signal and a protease domain. It belongs to a S41 protease family conserved across different filamentous fungi in the Ascomycota class and may constitute a core effector. AvrStb9 is also conserved among a global Z. tritici population and carries multiple amino acid substitutions caused by strong positive diversifying selection. These results demonstrate the contribution of an 'atypical' conserved effector protein to fungal avirulence and the role of sequence diversification in the escape of host recognition, adding to our understanding of host-pathogen interactions and the evolutionary processes underlying pathogen adaptation.

Aggressiveness and Patulin Production in Penicillium expansum Multidrug Resistant Strains with Different Expression Levels of MFS and ABC Transporters, in the Presence or Absence of Fludioxonil

Penicillium expansum is the most common postharvest pathogen of apple fruit, causing blue mold disease. Due to the extensive use of fungicides, strains resistant to multiple chemical classes have been selected. A previous study by our group proposed that the overexpression of MFS (major facilitator superfamily) and ABC (ATP binding cassette) transporters constitute an alternative resistance mechanism in Multi Drug resistant (MDR) strains of this pathogen. This study was initiated to determine two main biological fitness parameters of MDR strains: aggressiveness against apple fruit and patulin production. In addition, the expression pattern of efflux transporters and hydroxylase-encoding genes that belong to the patulin biosynthesis pathway, in the presence or absence of fludioxonil and under in vitro and in vivo conditions were investigated. Results showed that the MDR strains produced higher concentrations of patulin but showed a lower pathogenicity compared to the wild-type isolates. Moreover, expression analysis of patC, patM and patH genes indicated that the higher expression levels do not correlate with the detected patulin concentration. The selection of MDR strains in P. expansum populations and the fact that they produce more patulin, constitutes a serious concern not only for successful disease control but also for human health. The above-mentioned data represent the first report of MDR in P. expansum associated with its patulin-production ability and the expression level of patulin biosynthesis pathway genes.

Antifungal alternation can be beneficial for durability but at the cost of generalist resistance

The evolution of resistance to pesticides is a major burden in agriculture. Resistance management involves maximizing selection pressure heterogeneity, particularly by combining active ingredients with different modes of action. We tested the hypothesis that alternation may delay the build-up of resistance not only by spreading selection pressure over longer periods, but also by decreasing the rate of evolution of resistance to alternated fungicides, by applying an experimental evolution approach to the economically important crop pathogen Zymoseptoria tritici. Our results show that alternation is either neutral or slows the overall resistance evolution rate, relative to continuous fungicide use, but results in higher levels of generalism in evolved lines. We demonstrate that the nature of the fungicides, and therefore their relative intrinsic risk of resistance may underly this trade-off, more so than the number of fungicides and the rhythm of alternation. This trade-off is also dynamic over the course of resistance evolution. These findings open up new possibilities for tailoring resistance management effectively while optimizing interplay between alternation components.

Recent transposable element bursts are associated with the proximity to genes in a fungal plant pathogen

The activity of transposable elements (TEs) contributes significantly to pathogen genome evolution. TEs often destabilize genome integrity but may also confer adaptive variation in pathogenicity or resistance traits. De-repression of epigenetically silenced TEs often initiates bursts of transposition activity that may be counteracted by purifying selection and genome defenses. However, how these forces interact to determine the expansion routes of TEs within a pathogen species remains largely unknown. Here, we analyzed a set of 19 telomere-to-telomere genomes of the fungal wheat pathogen Zymoseptoria tritici. Phylogenetic reconstruction and ancestral state estimates of individual TE families revealed that TEs have undergone distinct activation and repression periods resulting in highly uneven copy numbers between genomes of the same species. Most TEs are clustered in gene poor niches, indicating strong purifying selection against insertions near coding sequences, or as a consequence of insertion site preferences. TE families with high copy numbers have low sequence divergence and strong signatures of defense mechanisms (i.e., RIP). In contrast, small non-autonomous TEs (i.e., MITEs) are less impacted by defense mechanisms and are often located in close proximity to genes. Individual TE families have experienced multiple distinct burst events that generated many nearly identical copies. We found that a Copia element burst was initiated from recent copies inserted substantially closer to genes compared to older copies. Overall, TE bursts tended to initiate from copies in GC-rich niches that escaped inactivation by genomic defenses. Our work shows how specific genomic environments features provide triggers for TE proliferation in pathogen genomes.

Fungicide sensitivity levels in the Lithuanian Zymoseptoria tritici population in 2021

Zymoseptoria tritici causes the disease known as septoria leaf blotch in winter wheat and is a major factor in yield loss worldwide. Farmers are inclined to use fungicides to protect their crops; however, the efficacy of these measures is rapidly decreasing due to the natural mechanisms of mutation emergence in pathogen populations. Increasing fungicide resistance is being recorded worldwide, therefore, screening of the current situation in Lithuania is essential to determine the subsequent steps of crop protection strategies. In this study, in vitro fungicide sensitivity tests, mutation detection, and field experiments were carried out. The mean EC₅₀ values for prothioconazole-desthio and mefentrifluconazole were 0.14 and 0.28 mg/l, respectively. Increased frequency of the mutation S524T, linked to DMIs resistance, was observed. Results revealed that the dominant point mutation in the gene CYP51 was I381V, and the most frequent CYP51 haplotype was D13 (V136C, I381V, Y461H, S524T). The mutation G143A, linked to QoI resistance, was detected in ¾ of the population. Mutations conferring resistance to SDHIs were not detected in single pycnidium isolates. Two-year field experiments likewise showed no decline in field efficacy of SDHI fungicide in Lithuania. Moreover, the baseline sensitivity of the Lithuanian Z. tritici population to QiI fungicide fenpicoxamid was established. The findings of this study provide an update on the current status of fungicide resistance in the Lithuanian Z. tritici population.

Do mitochondria use efflux pumps to protect their ribosomes from antibiotics?

Fungal environments are rich in natural and engineered antimicrobials, and this, combined with the fact that fungal genomes are rich in coding sequences for transporters, suggests that fungi are an intriguing group in which to search for evidence of antimicrobial efflux pumps in mitochondria. Herein, the range of protective mechanisms used by fungi against antimicrobials is introduced, and it is hypothesized, based on the susceptibility of mitochondrial and bacterial ribosomes to the same antibiotics, that mitochondria might also contain pumps that efflux antibiotics from these organelles. Preliminary evidence of ethidium bromide efflux is presented and several candidate efflux pumps are identified in fungal mitochondrial proteomes.

Baseline sensitivity of European Zymoseptoria tritici populations to the complex III respiration inhibitor fenpicoxamid

BACKGROUND

Fenpicoxamid is a recently developed fungicide belonging to the quinone inside inhibitor (QiI) group. This is the first fungicide within this group to be active against the Zymoseptoria tritici, which causes Septoria tritici blotch on wheat. The occurrence of pre-existing resistance mechanisms was monitored, using sensitivity assays and Illumina sequencing, in Z. tritici populations sampled in multiple European countries before the introduction of fenpicoxamid.

RESULTS

Although differences in sensitivity to all three fungicides tested (fenpicoxamid, fentin chloride and pyraclostrobin) existed between the isolate collections, no alterations associated with QiI resistance were detected. Among the isolates, a range in sensitivity to fenpicoxamid was observed (ratio between most sensitive/least sensitive = 53.1), with differences between the most extreme isolates when tested in planta following limited fenpicoxamid treatment. Sensitivity assays using fentin chloride suggest some of the observed differences in fenpicoxamid sensitivity are associated with multi-drug resistance. Detailed monitoring of the wider European population using Illumina-based partial sequencing of the Z. tritici also only detected the presence of G143A, with differences in frequencies of this alteration observed across the region.

CONCLUSIONS

This study provides a baseline sensitivity for European Z. tritici populations to fenpicoxamid. Target-site resistance appears to be limited or non-existing in European Z. tritici populations prior to the introduction of fenpicoxamid. Non-target site resistance mechanisms exist, but their impact in the field is predicted to be limited. © 2022 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Rapidly Increasing Boscalid Resistance in Corynespora cassiicola in China

Corynespora leaf spot caused by Corynespora cassiicola is an important foliar disease in cucumber. Succinate dehydrogenase inhibitors are the main fungicides used to control this disease. With the application of succinate dehydrogenase inhibitors (SDHIs) in the field, boscalid-resistant isolates have been continuously detected in the field. Resistance monitoring programs were performed to investigate the frequency and genotypes of resistant isolates. In our resistance monitoring, the frequency of resistant isolates rapidly increased from 9.68 to 85.88% in 2005 to 2020. Nine genotypes conferring SDHI resistance were found in resistant isolates, with different levels of resistance to SDHIs: B-H278R, B-H278L, B-H278Y, B-I280V, C-N75S, C-S73P, D-D95E, D-H105R, and D-G109V. The first sdh mutation was detected in Hebei Province in China, conferring an amino acid substitution at codon 278 in the sdhB subunit from histidine to tyrosine (B-H278Y), and it was the dominant resistance genotype in 2014 to 2015. Subsequently, other genotypes were gradually detected in the field, and the dominant mutations varied across years and across regions. The newest genotype (B-H278L) conferring SDHI resistance was found in 2020. To the best of our knowledge, this is the first report of C. cassiicola in cucumber. To date, multiple resistance to SDHIs, quinone outside inhibitors, benzimidazole fungicides, and dicarboximide fungicides have been detected, accounting for 75.64% of SDHI-resistant isolates. Therefore, the above four fungicides must be strictly restricted, and further monitoring work in other provinces with more isolates should be performed in the future.

RNA-Seq Provides Insights into the Mechanisms Underlying Ilyonectria robusta Responding to Secondary Metabolites of Bacillus methylotrophicus NJ13

(1) Background: Ilyonectria robusta can cause ginseng to suffer from rusty root rot. Secondary metabolites (SMs) produced by Bacillus methylotrophicus NJ13 can inhibit the mycelial growth of I. robusta. However, the molecular mechanism of the inhibition and response remains unclear. (2) Methods: Through an in vitro trial, the effect of B. methylotrophicus NJ13’s SMs on the hyphae and conidia of I. robusta was determined. The change in the physiological function of I. robusta was evaluated in response to NJ13’s SMs by measuring the electrical conductivity, malondialdehyde (MDA) content, and glucose content. The molecular interaction mechanism of I. robusta’s response to NJ13’s SMs was analyzed by using transcriptome sequencing. (3) Results: NJ13’s SMs exhibited antifungal activity against I. robusta: namely, the hyphae swelled and branched abnormally, and their inclusions leaked out due to changes in the cell membrane permeability and the peroxidation level; the EC50 value was 1.21% (v/v). In transcripts at 4 dpi and 7 dpi, the number of differentially expressed genes (DEGs) (|log2(fold change)| > 1, p adj ≤ 0.05) was 1960 and 354, respectively. NJ13’s SMs affected the glucose metabolism pathway, and the sugar-transporter-related genes were downregulated, which are utilized by I. robusta for energy production. The cell wall structure of I. robusta was disrupted, and chitin-synthase-related genes were downregulated. (4) Conclusions: A new dataset of functional responses of the ginseng pathogenic fungus I. robusta was obtained. The results will benefit the development of targeted biological fungicides for I. robusta and the study of the molecular mechanisms of interaction between biocontrol bacteria and phytopathogenic fungi.

LC-MS/MS-Based Fungicide Accumulation Assay to Demonstrate Efflux Activity in the Wheat Pathogen Zymoseptoria tritici

Increased drug efflux compromises the efficacy of a large panel of treatments in the clinic against cancer or bacterial, fungal, and viral diseases, and in agriculture due to the emergence of multidrug-resistant pathogenic fungi. Until recently, to demonstrate increased drug efflux, the use of labeled drugs or fluorescent dyes was necessary. With the increasing sensitivity of detection devices, direct assessment of drug efflux has become realistic. Here, we describe a medium-throughput method to assess the intracellular drug concentration in the plant pathogenic fungus Zymoseptoria tritici cultivated in the presence of a sublethal fungicide concentration. As a model fungicide, we used the succinate-dehydrogenase inhibitor boscalid. The boscalid concentration was assessed in the different culture fractions using mass spectrometry linked to liquid chromatography (LC-MS/MS). The ratio between the intracellular and total boscalid amount was used as an inversed proxy for the efflux activity. Using isogenic mutant strains known for their differential efflux capacities, we validated the negative correlation between the intracellular boscalid concentration and efflux activity. In addition, intra-cellular fungicide accumulation explains the susceptibility of the tested strains to boscalid. This assay may be useful in lead development when a new molecule displays good inhibitory activity against its isolated target protein but fails to control the target organism.

Comparative transcriptome analysis reveals the resistance regulation mechanism and fungicidal activity of the fungicide phenamacril in Fusarium oxysporum

Fusarium oxysporum (Fo) is an important species complex of soil-borne pathogenic fungi that cause vascular wilt diseases of agricultural crops and some opportunistic diseases of humans. The fungicide phenamacril has been extensively reported to have antifungal activity against Fusarium graminearum and Fusarium fujikuroi. In this study, we found that the amino acid substitutions (V151A and S418T) in Type I myosin FoMyo5 cause natural low resistance to phenamacril in the plant pathogenic Fo isolates. Therefore, we compared the transcriptomes of two phenamacril-resistant Fo isolates FoII5, Fo1st and one phenamacril-sensitive isolate Fo3_a after 1 μg/mL phenamacril treatment. Among the 2728 differentially expressed genes (DEGs), 14 DEGs involved in oxidation–reduction processes and MFS transporters, were significantly up-regulated in phenamacril-resistant isolates. On the other hand, 14 DEGs involved in ATP-dependent RNA helicase and ribosomal biogenesis related proteins, showed significantly down-regulated expression in both phenamacril-resistant and -sensitive isolates. These results indicated that phenamacril not only seriously affected the cytoskeletal protein binding and ATPase activity of sensitive isolate, but also suppressed ribosome biogenesis in all the isolates. Hence, this study helps us better understand resistance regulation mechanism and fungicidal activity of phenamacril and provide reference for the development of new fungicides to control Fo.

Cytochrome P450 and Glutathione S-Transferase Confer Metabolic Resistance to SYP-14288 and Multi-Drug Resistance in Rhizoctonia solani

SYP-14288 is a fungicide as an uncoupler of oxidative phosphorylation, which is effective in controlling fungal pathogens like Rhizoctonia solani. To determine whether R. solani can develop SYP-14288 resistance and possibly multi-drug resistance (MDR), an SYP-14288-resistant mutant of R. solani X19-7 was generated from wild-type strain X19, and the mechanism of resistance was studied through metabolic and genetic assays. From metabolites of R. solani treated with SYP-14288, three compounds including M1, M2, and M3 were identified according to UPLC-MS/MS analysis, and M1 accumulated faster than M2 and M3 in X19-7. When X19-7 was treated by glutathione-S-transferase (GST) inhibitor diethyl maleate (DEM) and SYP-14288 together, or by DEM plus one of tested fungicides that have different modes of action, a synergistic activity of resistance occurred, implying that GSTs promoted metabolic resistance against SYP-14288 and therefore led to MDR. By comparing RNA sequences between X19-7 and X19, six cytochrome P450s (P450s) and two GST genes were selected as a target, which showed a higher expression in X19-7 than X19 both before and after the exposure to SYP-14288. Furthermore, heterologous expression of P450 and GST genes in yeast was conducted to confirm genes involved in metabolic resistance. In results, the P450 gene AG1IA_05136 and GST gene AG1IA_07383 were related to fungal resistance to multiple fungicides including SYP-14288, fluazinam, chlorothalonil, and difenoconazole. It was the first report that metabolic resistance of R. solani to uncouplers was associated with P450 and GST genes.

Molecular Mechanisms Associated with the Resistance of Rhizoctonia solani AG-4 Isolates to the Succinate Dehydrogenase Inhibitor Thifluzamide

Thifluzamide, a succinate dehydrogenase (SDH) inhibitor, possesses high activity against Rhizoctonia. In this study, 144 Rhizoctonia solani AG-4 (4HGI, 4HGII, and 4HGIII) isolates, the predominate pathogen associated with sugar beet seedling damping-off, were demonstrated to be sensitive to thifluzamide with a calculated mean median effective concentration of 0.0682 ± 0.0025 μg/ml. Thifluzamide-resistant isolates were generated using fungicide-amended media, resulting in four AG-4HGI isolates and eight AG-4HGII isolates with stable resistance and almost no loss in fitness. Evaluation of cross-resistance of the 12 thifluzamide-resistant isolates and their corresponding parental-sensitive isolates revealed a moderately positive correlation between thifluzamide resistance and the level of resistance to eight other fungicides from three groups, the exception being fludioxonil. An active efflux of fungicide through ATP-binding cassette and major facilitator superfamily transporters was found to be correlated to the resistance of R. solani AG-4HGII isolates to thifluzamide based on RNA-sequencing and quantitative reverse transcription-PCR analyses. Sequence analysis of sdhA, sdhB, sdhC, and sdhD revealed replacement of isoleucine by phenylalanine at position 61 in SDHC in 9 of the 12 generated thifluzamide-resistant isolates. No other mutations were found in any of the other genes. Collectively, the data indicate that the active efflux of fungicide and a point mutation in sdhC may contribute to the resistance of R. solani AG-4HGI and AG-4HGII isolates to thifluzamide in vitro. This is the first characterization of the potential molecular mechanism associated with the resistance of R. solani AG-4 isolates to thifluzamide and provides practical guidance for the use of this fungicide.

Quantitative multiplexed proteomics analysis reveals reshaping of the lysine 2-hydroxyisobutyrylome in Fusarium graminearum by tebuconazole

Backgrounds

Lysine 2-hydroxyisobutyrylation (Khib) is a newly discovered posttranslational modification (PTM) and has been identified in several prokaryotic and eukaryotic organisms. Fusarium graminearum, a major pathogen of Fusarium head blight (FHB) in cereal crops, can cause considerable yield loss and produce various mycotoxins that threaten human health. The application of chemical fungicides such as tebuconazole (TEC) remains the major method to control this pathogen. However, the distribution of Khib in F. graminearum and whether Khib is remodified in response to fungicide stress remain unknown.

Results

Here, we carried out a proteome-wide analysis of Khib in F. graminearum, identifying the reshaping of the lysine 2-hydroxyisobutyrylome by tebuconazole, using the most recently developed high-resolution LC–MS/MS technique in combination with high-specific affinity enrichment. Specifically, 3501 Khib sites on 1049 proteins were identified, and 1083 Khib sites on 556 modified proteins normalized to the total protein content were changed significantly after TEC treatment. Bioinformatics analysis showed that Khib proteins are involved in a wide range of biological processes and may be involved in virulence and deoxynivalenol (DON) production, as well as sterol biosynthesis, in F. graminearum.

Conclusions

Here, we provided a wealth of resources for further study of the roles of Khib in the fungicide resistance of F. graminearum. The results enhanced our understanding of this PTM in filamentous ascomycete fungi and provided insight into the remodification of Khib sites during azole fungicide challenge in F. graminearum.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08372-4.

The Transcription Factor FgAtrR Regulates Asexual and Sexual Development, Virulence, and DON Production and Contributes to Intrinsic Resistance to Azole Fungicides in Fusarium graminearum

Simple Summary

Fusarium graminearum is a devastating plant pathogen that can cause wheat head blight. Azole fungicides are commonly used chemicals for control of this disease. However, F. graminearum strains resistant to these fungicides have emerged. To better understand the azole resistance mechanism of F. graminearum, we identified and characterized the Zn(II)2-Cys6 transcription factor FgAtrR in F. graminearum. We found that FgAtrR played critical roles in vegetative growth, conidia production, perithecium formation, and virulence on wheat heads and corn silks. FgAtrR was also involved in the resistance to azole antifungals by regulating the expression of the drug target FgCYP51s and efflux pump transporters. These results broadened our understanding of the azole resistance mechanisms of F. graminearum.

Abstract

Fusarium graminearum is the predominant causal agent of cereal Fusarium head blight disease (FHB) worldwide. The application of chemical fungicides such as azole antifungals is still the primary method for FHB control. However, to date, our knowledge of transcriptional regulation in the azole resistance of F. graminearum is quite limited. In this study, we identified and functionally characterized a Zn(II)2-Cys6 transcription factor FgAtrR in F. graminearum. We constructed a FgAtrR deletion mutant and found that deletion of FgAtrR resulted in faster radial growth with serious pigmentation defects, significantly reduced conidial production, and an inability to form perithecia. The pathogenicity of the ΔFgAtrR mutant on wheat spikes and corn silks was severely impaired with reduced deoxynivalenol production, while the tolerance to prochloraz and propiconazole of the deletion mutant was also significantly decreased. RNA-seq indicated that many metabolic pathways were affected by the deletion of FgAtrR. Importantly, FgAtrR could regulate the expression of the FgCYP51A and ABC transporters, which are the main contributors to azole resistance. These results demonstrated that FgAtrR played essential roles in asexual and sexual development, DON production, and pathogenicity, and contributed to intrinsic resistance to azole fungicides in F. graminearum. This study will help us improve the understanding of the azole resistance mechanism in F. graminearum.

Molecular characterization and overexpression of the difenoconazole resistance gene CYP51 in Lasiodiplodia theobromae field isolates

Stem-end rot (SER) caused by Lasiodiplodia theobromae is an important disease of mango in China. Demethylation inhibitor (DMI) fungicides are widely used for disease control in mango orchards. The baseline sensitivity to difenoconazole of 138 L. theobromae isolates collected from mango in the field in 2019 was established by the mycelial growth rate method. The cross-resistance to six site-specific fungicides with different modes of action were investigated using 20 isolates randomly selected. The possible mechanism for L. theobromae resistance to difenoconazole was preliminarily determined through gene sequence alignment and quantitative real-time PCR analysis. The results showed that the EC₅₀ values of 138 L. theobromae isolates to difenoconazole ranged from 0.01 to 13.72 µg/mL. The frequency of difenoconazole sensitivity formed a normal distribution curve when the outliers were excluded. Difenoconazole showed positive cross-resistance only with the DMI tebuconazole but not with non-DMI fungicides carbendazim, pyraclostrobin, fludioxonil, bromothalonil, or iprodione. Some multifungicide-resistant isolates of L. theobromae were found. Two amino acid substitutions (E209k and G207A) were found in the CYP51 protein, but they were unlikely to be related to the resistance phenotype. There was no alteration in the promoter region of the CYP51 gene. However, difenoconazole significantly increased the expression of the CYP51 gene in the resistant isolates compared to the susceptible isolates. These results are vital to develop effective mango disease management strategies to avoid the development of further resistance.

Spatio-temporal distribution of DMI and SDHI fungicide resistance of Zymoseptoria tritici throughout Europe based on frequencies of key target-site alterations

BACKGROUND

Over the past decade, demethylation inhibitor (DMI) and succinate dehydrogenase inhibitor (SDHI) fungicides have been extensively used to control to septoria tritici blotch, caused by Zymoseptoria tritici on wheat. This has led to the development and selection of alterations in the target-site enzymes (CYP51 and SDH, respectively).

RESULTS

Taking advantage of newly and previously developed qPCR assays, the frequency of key alterations associated with DMI (CYP51-S524T) and SDHI (SDHC-T79N/I, C-N86S and C-H152R) resistance was assessed in Z. tritici-infected wheat leaf samples collected from commercial crops (n = 140) across 14 European countries prior to fungicide application in the spring of 2019. This revealed the presence of a West to East gradient in the frequencies of the most common key alterations conferring azole (S524T) and SDHI resistance (T79N and N86S), with the highest frequencies measured in Ireland and Great Britain. These observations were corroborated by sequencing (CYP51 and SDH subunits) and sensitivity phenotyping (prothioconazole-desthio and fluxapyroxad) of Z. tritici isolates collected from a selection of field samples. Additional sampling made at the end of the 2019 season confirmed the continued increase in frequency of the targeted alterations. Investigations on historical leaf DNA samples originating from different European countries revealed that the frequency of all key alterations (except C-T79I) has been gradually increasing over the past decade.

CONCLUSION

Whilst these alterations are quickly becoming dominant in Ireland and Great Britain, scope still exists to delay their selection throughout the wider European population, emphasizing the need for the implementation of fungicide antiresistance measures. © 2021 Society of Chemical Industry.

Are Efficient-Dose Mixtures a Solution to Reduce Fungicide Load and Delay Evolution of Resistance? An Experimental Evolutionary Approach

Pesticide resistance poses a critical threat to agriculture, human health and biodiversity. Mixtures of fungicides are recommended and widely used in resistance management strategies. However, the components of the efficiency of such mixtures remain unclear. We performed an experimental evolutionary study on the fungal pathogen Z. tritici to determine how mixtures managed resistance. We compared the effect of the continuous use of single active ingredients to that of mixtures, at the minimal dose providing full control of the disease, which we refer to as the "efficient" dose. We found that the performance of efficient-dose mixtures against an initially susceptible population depended strongly on the components of the mixture. Such mixtures were either as durable as the best mixture component used alone, or worse than all components used alone. Moreover, efficient dose mixture regimes probably select for generalist resistance profiles as a result of the combination of selection pressures exerted by the various components and their lower doses. Our results indicate that mixtures should not be considered a universal strategy. Experimental evaluations of specificities for the pathogens targeted, their interactions with fungicides and the interactions between fungicides are crucial for the design of sustainable resistance management strategies.

Genome-Wide Analysis of Major Facilitator Superfamily and Its Expression in Response of Poplar to Fusarium oxysporum

The major facilitator superfamily (MFS) is one of the largest known membrane transporter families. MFSs are involved in many essential functions, but studies on the MFS family in poplar have not yet been reported. Here, we identified 41 MFS genes from Populus trichocarpa (PtrMFSs). We built a phylogenetic tree, which clearly divided members of PtrMFS into six groups with specific gene structures and protein motifs/domains. The promoter regions contain various cis-acting elements involved in stress and hormone responsiveness. Genes derived from segmental duplication events are unevenly distributed in 17 poplar chromosomes. Collinearity analysis showed that PtrMFS genes are conserved and homologous to corresponding genes from four other species. Transcriptome data indicated that 40 poplar MFS genes were differentially expressed when treated with Fusarium oxysporum. Co-expression networks and gene function annotations of MFS genes showed that MFS genes tightly co-regulated and closely related in function of transmembrane transport. Taken together, we systematically analyzed structure and function of genes and proteins in the PtrMFS family. Evidence indicated that poplar MFS genes play key roles in plant development and response to a biological stressor.

Florylpicoxamid, a new picolinamide fungicide with broad spectrum activity

BACKGROUND

Following the introduction of fenpicoxamid, a natural product-based fungicide targeting the Qi site of mitochondrial cytochrome bc1 complex, a second generation fully synthetic picolinamide, florylpicoxamid, was discovered and its biological activity and attributes were characterized.

RESULTS

In vitro fungal growth inhibition assays and in planta glasshouse biological activity evaluations showed florylpicoxamid was active against 21 different plant pathogenic fungi within the phyla Ascomycota and Basidiomycota. Among the pathogens evaluated, florylpicoxamid was most potent against Zymoseptoria tritici, the causal organism of wheat leaf blotch, providing 80% growth inhibition in vitro at 0.0046 mg L^-1 and 80% disease control in planta at 0.03 mg L^-1 when applied as a preventative treatment. Florylpicoxamid was more efficacious than epoxiconazole, fluxapyroxad, and benzovindiflupyr versus a Z. tritici wild-type isolate when applied as curative and preventative treatments, with superior 10-day curative reachback activity. Analytical studies and in planta tests demonstrated that florylpicoxamid partitioned into plants quickly and showed good systemicity and translaminar activity on both monocot and dicot plants. No cross-resistance was observed between florylpicoxamid and strobilurin or azole fungicides. Florylpicoxamid exerts its preventative effect by preventing spore germination on the leaf surface and curative activity by arresting mycelial growth and pycnidia development in leaf tissue.

CONCLUSIONS

With strong broad spectrum fungicidal activity, florylpicoxamid delivers an innovative solution for growers to sustain high productivity and quality of many crops, and also provides a new option for developing effective strategies for fungicide resistance management. © 2021 Society of Chemical Industry.

Molecular Mechanisms Underlying Fungicide Resistance in Citrus Postharvest Green Mold

The necrotrophic fungus Penicillium digitatum (Pd) is responsible for the green mold disease that occurs during postharvest of citrus and causes enormous economic losses around the world. Fungicides remain the main method used to control postharvest green mold in citrus fruit storage despite numerous occurrences of resistance to them. Hence, it is necessary to find new and more effective strategies to control this type of disease. This involves delving into the molecular mechanisms underlying the appearance of resistance to fungicides during the plant–pathogen interaction. Although mechanisms involved in resistance to fungicides have been studied for many years, there have now been great advances in the molecular aspects that drive fungicide resistance, which facilitates the design of new means to control green mold. A wide review allows the mechanisms underlying fungicide resistance in Pd to be unveiled, taking into account not only the chemical nature of the compounds and their target of action but also the general mechanism that could contribute to resistance to others compounds to generate what we call multidrug resistance (MDR) phenotypes. In this context, fungal transporters seem to play a relevant role, and their mode of action may be controlled along with other processes of interest, such as oxidative stress and fungal pathogenicity. Thus, the mechanisms for acquisition of resistance to fungicides seem to be part of a complex framework involving aspects of response to stress and processes of fungal virulence.

A population-level invasion by transposable elements triggers genome expansion in a fungal pathogen

Genome evolution is driven by the activity of transposable elements (TEs). The spread of TEs can have deleterious effects including the destabilization of genome integrity and expansions. However, the precise triggers of genome expansions remain poorly understood because genome size evolution is typically investigated only among deeply divergent lineages. Here, we use a large population genomics dataset of 284 individuals from populations across the globe of Zymoseptoria tritici, a major fungal wheat pathogen. We built a robust map of genome-wide TE insertions and deletions to track a total of 2456 polymorphic loci within the species. We show that purifying selection substantially depressed TE frequencies in most populations, but some rare TEs have recently risen in frequency and likely confer benefits. We found that specific TE families have undergone a substantial genome-wide expansion from the pathogen's center of origin to more recently founded populations. The most dramatic increase in TE insertions occurred between a pair of North American populations collected in the same field at an interval of 25 years. We find that both genome-wide counts of TE insertions and genome size have increased with colonization bottlenecks. Hence, the demographic history likely played a major role in shaping genome evolution within the species. We show that both the activation of specific TEs and relaxed purifying selection underpin this incipient expansion of the genome. Our study establishes a model to recapitulate TE-driven genome evolution over deeper evolutionary timescales.

Directed evolution predicts cytochrome b G37V target site modification as probable adaptive mechanism towards the QiI fungicide fenpicoxamid in Zymoseptoria tritici

Acquired resistance is a threat to antifungal efficacy in medicine and agriculture. The diversity of possible resistance mechanisms and highly adaptive traits of pathogens make it difficult to predict evolutionary outcomes of treatments. We used directed evolution as an approach to assess the resistance risk to the new fungicide fenpicoxamid in the wheat pathogenic fungus Zymoseptoria tritici. Fenpicoxamid inhibits complex III of the respiratory chain at the ubiquinone reduction site (Qi site) of the mitochondrially encoded cytochrome b, a different site than the widely used strobilurins which inhibit the same complex at the ubiquinol oxidation site (Q_o site). We identified the G37V change within the cytochrome b Q_i site as the most likely resistance mechanism to be selected in Z. tritici. This change triggered high fenpicoxamid resistance and halved the enzymatic activity of cytochrome b, despite no significant penalty for in vitro growth. We identified negative cross-resistance between isolates harbouring G37V or G143A, a Q_o site change previously selected by strobilurins. Double mutants were less resistant to both QiIs and quinone outside inhibitors compared to single mutants. This work is a proof of concept that experimental evolution can be used to predict adaptation to fungicides and provides new perspectives for the management of QiIs.

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.

Grapevine Powdery Mildew: Fungicides for Its Management and Advances in Molecular Detection of Markers Associated with Resistance

Grapevine powdery mildew is a principal fungal disease of grapevine worldwide. Even though it usually does not cause plant death directly, heavy infections can lead to extensive yield losses, and even low levels of the disease can negatively affect the quality of the wine. Therefore, intensive spraying programs are commonly applied to control the disease, which often leads to the emergence and spread of powdery mildew strains resistant to different fungicides. In this review, we describe major fungicide classes used for grapevine powdery mildew management and the most common single nucleotide mutations in target genes known to confer resistance to different classes of fungicides. We searched the current literature to review the development of novel molecular methods for quick detection and monitoring of resistance to commonly used single-site fungicides against Erysiphe necator. We analyze and compare the developed methods. From our investigation it became evident that this research topic has been strongly neglected and we hope that effective molecular methods will be developed also for resistance monitoring in biotroph pathogens.

Development of a PacBio Long-Read Sequencing Assay for High Throughput Detection of Fungicide Resistance in Zymoseptoria tritici

Fungicide resistance has become a challenging problem in management of Septoria tritici blotch (STB), caused by Zymoseptoria tritici, the most destructive disease of winter wheat throughout western and northern Europe. To ensure the continued effectiveness of those fungicides currently used, it is essential to monitor the development and spread of such resistance in field populations of the pathogen. Since resistance to the key families of fungicides used for STB control (demethyalation inhibitors or azoles, succinate dehydrogenase inhibitors or SDHIs and Quinone outside Inhibitors or QoIs) is conferred through target-site mutations, the potential exists to monitor resistance through the molecular detection of alterations in the target site genes. As more efficient fungicides were developed and applied, the pathogen has continuously adapted through accumulating multiple target-site alterations. In order to accurately monitor these changes in field populations, it is therefore becoming increasingly important to completely sequence the targeted genes. Here we report the development of a PacBio assay that facilitates the multiplex amplification and long-read sequencing of the target gene(s) for the azole (CYP51), SDHI (Sdh B, C, and D), and QoI (cytochrome b) fungicides. The assay was developed and optimised using three Irish Z. tritici collections established in spring 2017, which capture the range of fungicide resistance present in modern European populations of Z. tritici. The sequences obtained through the PacBio assay were validated using traditional Sanger sequencing and in vitro sensitivity screenings. To further exploit the long-read and high throughput potential of PacBio sequencing, an additional nine housekeeping genes (act, BTUB, cal, cyp, EF1, GAPDH, hsp80-1, PKC, TFC1) were sequenced and used to provide comprehensive Z. tritici strain genotyping.

Changes in DMI, SDHI, and QoI Fungicide Sensitivity in the Estonian Zymoseptoria tritici Population between 2019 and 2020

Zymoseptoria tritici (Zt) populations adapt under the selection pressure of fungicides applied for disease control. The primary objective of this study was to assess fungicide sensitivity in the Estonian Zt population. A total of 282 Zt isolates from 2019 and 2020 were tested for sensitivity to azoles (DMIs; prothioconazole-desthio, epoxiconazole, mefentrifluconazole) and succinate dehydrogenase inhibitors (SDHIs; boscalid, fluxapyroxad). The efficacy of the tested fungicides varied considerably between the Estonian counties, but the Zt population is mainly sensitive to DMIs. Additionally, the frequencies of CYP51 gene alterations varied; D134G, V136C, A379G, and S524T had increased, but V136A and I381V showed a moderate decrease in 2020 in comparison to 2019. Sensitivity to SDHIs was stable, but boscalid was less effective than fluxapyroxad. SdhC gene mutations C-T33N, C-T34N, and C-N86S were common, but not linked with SDHI fungicide sensitivity assay results. Otherwise, mutation B-N225I in the SdhB subunit occurred in isolates with reduced sensitivity to SDHIs. Sensitivity to strobilurins was evaluated by the mutation G143A in the CytB gene, which was present in nearly half of the population. The data presented confirm the ongoing evolution of fungicide sensitivity in the Zt population in Estonia and highlight the importance of knowledge-based decisions for optimizing anti-resistance strategies in the field.

Insights into the multitrophic interactions between the biocontrol agent Bacillus subtilis MBI 600, the pathogen Botrytis cinerea and their plant host

Botrytis cinerea is a plant pathogen causing the gray mold disease in a plethora of host plants. The control of the disease is based mostly on chemical pesticides, which are responsible for environmental pollution, while they also pose risks for human health. Furthermore, B. cinerea resistant isolates have been identified against many fungicide groups, making the control of this disease challenging. The application of biocontrol agents can be a possible solution, but requires deep understanding of the molecular mechanisms in order to be effective. In this study, we investigated the multitrophic interactions between the biocontrol agent Bacillus subtilis MBI 600, a new commercialized biopesticide, the pathogen B. cinerea and their plant host. Our analysis showed that this biocontrol agent reduced B. cinerea mycelial growth in vitro, and was able to suppress the disease incidence on cucumber plants. Moreover, treatment with B. subtilis led to induction of genes involved in plant immunity. RNA-seq analysis of B. cinerea transcriptome upon exposure to bacterial secretome, showed that genes coding for MFS and ABC transporters were highly induced. Deletion of the Bcmfs1 MFS transporter gene, using a CRISP/Cas9 editing method, affected its virulence and the tolerance of B. cinerea to bacterial secondary metabolites. These findings suggest that specific detoxification transporters are involved in these interactions, with crucial role in different aspects of B. cinerea physiology.

Contrasting levels of genetic predictability in the evolution of resistance to major classes of fungicides

The evolution of resistance has been seen across all major classes of xenobiotics, including antimicrobial drugs and agricultural pesticides. This repeated emergence of resistance is a case of phenotypic parallel evolution, but often the parallelism extends to the molecular level too, with multiple species gaining the same mutation in response to the same chemical treatment. We review the degree of repeatability in target-site resistance mutations affecting different classes of site-specific agricultural fungicides used in crop protection, comparing the extent to which resistance in different pathogen species has evolved via the same or different mutations. For all major fungicide target sites, substantial levels of molecular parallel evolution can be seen, with at least one mutation recurring in over 50% of species. Target-site mutations appear to be most repeatable in cytochrome b, target site of quinone-outside inhibitor fungicides, and least predictable for CYP51, target site of the azoles. Intermediate levels of repeatability are seen for the MBC target site β-tubulin, and the SDHI target site succinate dehydrogenase. Repeatability may be lower where there are selective trade-offs between resistance and pleiotropic fitness penalties, or differing levels of cross-resistance across members of a fungicide class; or where single mutations confer only partial resistance, and epistatic interactions between multiple mutations result in a rugged fitness landscape. This affects the predictive power of in vitro mutation studies, and has practical implications for resistance monitoring strategies and diagnostic methods.

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.

Difenoconazole Resistance Shift in Botrytis cinerea From Tomato in China Associated With Inducible Expression of CYP51

Gray mold caused by Botrytis cinerea is one of the most important diseases in tomato. It can be controlled effectively by demethylation inhibitor (DMI) fungicides, but their resistance status after long-term use in the field is unclear. The baseline sensitivity to difenoconazole of 142 B. cinerea isolates from China with no history of DMI usage was characterized, with a mean effective concentration for 50% mycelial growth inhibition (EC₅₀) of 0.97 ± 0.50 μg/ml. EC₅₀ values for difenoconazole sensitivity of another 248 isolates collected in 2011 and 2016 ranged from 0.04 to 11.99 μg/ml, and the frequency of difenoconazole sensitivity formed a nonnormal distribution curve. Detached fruit studies revealed that isolates with EC₅₀ values of approximately 6.00 μg/ml were not controlled effectively. The mean EC₅₀ of the resistant isolates changed from 6.74 to 8.65 μg/ml between 2011 and 2016. Positive cross-resistance was only observed between difenoconazole and two DMIs. One dual resistant isolate and one triple resistant isolate were found among the difenoconazole-resistant isolates collected in 2016, associated with point mutations in corresponding target proteins of the fungicides azoxystrobin and fludioxonil. This indicated that B. cinerea not only showed higher difenoconazole resistance levels but gradually changed from single to multiple fungicide resistance over time. No amino acid variation was found in the CYP51 protein. In the absence of difenoconazole, the relative expression of CYP51 was not significantly different in sensitive and resistant isolates. Induced expression of CYP51 is an important determinant of DMI resistance in B. cinerea from tomato. However, nucleotide variants found in the upstream region had no association with the fungicide resistance phenotype. These results will be helpful for the management of B. cinerea in the field.