MPMI PhytoFrontiers Phytobiomes all journals
Spotlight: Disease Control and Integrated ManagementFree Access icon

Characterization of High Fludioxonil Resistance in Botrytis cinerea Isolates from Calibrachoa Flowers

    Affiliations
    Authors and Affiliations
    • Madeline Dowling1
    • Jhulia Gelain2
    • Louise Larissa May De Mio2
    • Guido Schnabel1
    1. 1Department of Plant and Environmental Sciences, Clemson University, Clemson, SC 29634, U.S.A.
    2. 2Department of Plant Sciences, Universidade Federal do Paraná, Curitiba, Brazil

    Published Online:https://doi.org/10.1094/PHYTO-07-20-0268-R

    Abstract

    The fungicide fludioxonil is one of the most effective single-site fungicides available for managing flower blight caused by Botrytis cinerea on fruit and ornamental crops. Although low and moderate levels of resistance to fludioxonil have been reported in the pathogen across the United States and Europe, high resistance has been reported only from greenhouses in China. In this study, two B. cinerea isolates with high resistance (half maximal effective concentration >100 µg/ml) to fludioxonil were detected on ornamental calibrachoa flowers grown in a greenhouse. These isolates exhibited stable resistance for >20 generations, produced symptoms on calibrachoa flowers sprayed with label rates of fludioxonil, and displayed in vitro fitness penalties with decreased mycelial growth (P < 0.0001) and sporulation (P < 0.0001) compared with sensitive isolates. Highly resistant isolates were identified as MDR1h, containing the ΔL/V497 deletion in mrr1. However, resistance levels and in vitro fitness parameter characteristics were not consistent with this phenotype. One isolate contained the mutation L267V between HAMP domains 1 and 2 of the Bos-1 gene, and both isolates exhibited high osmotic sensitivity and reduced glycerol accumulation in the presence of fludioxonil, indicating that high resistance of these isolates may be associated with the high-osmolarity glycerol mitogen-activated protein kinase pathway.

    The phenylpyrrole fungicide fludioxonil is an important mainstay for managing diseases caused by Botrytis cinerea on both ornamental and food crops before and after harvest (Förster et al. 2007; Rupp et al. 2017). Fludioxonil currently combines fungicidal efficacy with limited resistance development (Fernández-Ortuño et al. 2014; Kim et al. 2016; Petit et al. 2011; Rupp et al. 2017). This combination is particularly important considering the widespread resistance of B. cinerea to most major chemical classes, with some isolates having simultaneous resistance to seven chemical classes (Cosseboom et al. 2019; Fernández-Ortuño et al. 2014, 2015).

    Fludioxonil’s mode of action has not been fully delineated (Kilani and Fillinger 2016; Li et al. 2014). However, exposure to fludioxonil is associated with the activation of the high-osmolarity glycerol (HOG1) mitogen activated protein kinase (MAPK) pathway (Furukawa et al. 2012; Zhang et al. 2002). This pathway includes a two-component regulatory system composed of the sensor kinase, OS-1, which detects osmotic stress conditions, and a response regulator that receives signals from OS-1 and adjusts gene expression to regulate the cell’s response to osmotic stress (Furukawa et al. 2012; Viaud et al. 2006). Fludioxonil is believed to mimic an osmotic stress signal by binding to OS-1, triggering inappropriate activation, which results in overaccumulation of glycerol, ion fluxes, and abnormally high turgor pressure (Lew 2010; Zhang et al. 2002). Fungi exposed to fludioxonil often exhibit swollen hyphae with ruptured apexes, indicating that this increased turgor pressure ultimately leads to cell lysis and death (Leroux 1996; Zhang et al. 2002).

    Few documented instances of fungal isolates highly resistant (HR) to fludioxonil in the field have occurred to date (Kilani and Fillinger 2016). Most of these instances occurred in Alternaria spp., where high resistance to fludioxonil (half maximal effective concentration [EC50 > 10] µg/ml) has been reported in Europe and the United States (Avenot and Michailides 2015; Malandrakis et al. 2015). In one case, dual fludioxonil- and iprodione-resistant isolates made up 14% of a population (Malandrakis et al. 2015). For B. cinerea, fludioxonil HR isolates (EC50 > 5 µg/ml) have never been reported in the field but have been documented from greenhouses in China on strawberry, cucumber, and tomato (Ren et al. 2016; Sang et al. 2018; Zhou et al. 2020). In one analysis in China, HR isolates made up <3% of the Botrytis population from greenhouses included in the study (Ren et al. 2016). Mutations were detected in HR B. cinerea isolates in the HAMP domains of the OS-1 gene, called Bos-1 in B. cinerea, and some of these isolates were cross-resistant to dicarboximide fungicides (Ren et al. 2016). Generally, both laboratory and greenhouse HR isolates exhibited fitness penalties affecting sporulation and osmotic sensitivity (Ren et al. 2016; Sang et al. 2018; Zhou et al. 2020). To the best of our knowledge, high resistance to fludioxonil has not been reported outside China for B. cinerea on any crop.

    Although B. cinerea isolates highly resistant to fludioxonil are rarely found in the field, low and moderately resistant isolates of B. cinerea are often detected in ornamental flower greenhouses, strawberry fields, and vineyards in the United States and Europe (Cosseboom et al. 2019; Fernández-Ortuño et al. 2014; Kretschmer et al. 2009; Leroch et al. 2011, 2012; Li et al. 2014; Rupp et al. 2017). This resistance is not specific to fludioxonil but is drug efflux based multidrug resistance (MDR) associated with overexpression of the ABC transporter AtrB. This pump is regulated by the transcription factor mrr1, and mutations in the mrr1 gene result in two main phenotypes: MDR1, conferring low resistance, and MDR1h, conferring moderate resistance to fludioxonil (Hahn 2014) (Table 1).

    TABLE 1. Comparison of traits for isolates resistant or sensitive to fludioxonil from this study and previous literature reportsz

    While performing a Botrytis resistance monitoring program in 2019, we detected two isolates with abnormally high resistance to fludioxonil. The goal of this study was to characterize these isolates by investigating the level of resistance to fludioxonil and the dicarboximide fungicide iprodione, the stability of resistance, the field relevance of this resistance, the potential molecular sources of resistance, and the fitness of resistant isolates.

    MATERIALS AND METHODS

    Fungal isolates and culture conditions.

    B. cinerea spores were collected on sample swabs from infected calibrachoa leaf tissue in a commercial greenhouse in North Carolina. The plants had developed gray mold disease despite routine fungicide applications. Spores were transferred to media amended with discriminatory dosages of fungicides as described previously (Fernández-Ortuño et al. 2014). Based on their abnormally rapid mycelial growth on a 0.5 µg/ml discriminatory dose of fludioxonil, samples 5.8.19 and 5.10.19 were stored. For isolation, spores were transferred to malt extract agar (MEA) plates amended with 100 µg/ml fludioxonil to distinguish resistant spores. Spore cultures were incubated in the dark at 22°C for approximately 24 h. Individual germinating spores were then transferred to new MEA amended with 100 µg/ml fludioxonil to ensure collection of only highly resistant isolates (Table 2).

    TABLE 2. Name, origin, host, fludioxonil resistance phenotype, amino acid alterations, and source of isolates used in this study

    For comparison, three previously characterized Botrytis isolates from strawberry (Fernández-Ortuño et al. 2015) were revived (HP9, NC4, and FLOR5). These isolates were sensitive (S) to fludioxonil (no growth at 0.5 µg/ml) and iprodione (10 µg/ml) discriminatory dosages (Table 2).

    Characterization of sensitivity to fludioxonil and iprodione.

    To characterize fludioxonil and iprodione resistance, spore germination and mycelial growth were tested for isolates 5.8.19 and 5.10.19 (HR). For spore germination, 15- to 20-day-old cultures sporulating on 90-mm Petri dishes containing approximately 15 ml of MEA were washed with deionized sterile water, and the spore concentration adjusted to 3 × 105 spores/ml. Then 50 µl of this spore suspension was spread onto fungicide amended plates with a Drigalski spatula, and cultures were incubated 24 h at 22°C in the dark. For mycelial growth, 5-mm plugs were transferred from the margins of 3- to 5-day-old cultures to fungicide amended media and colony diameters measured after 2 days of incubation at 22°C in the dark. Both mycelial growth and spore germination tests were performed on MEA and potato dextrose agar (PDA) either unamended or amended with 100 µg/ml fludioxonil (Scholar; Syngenta, Greensboro, NC) and 100 µg/ml iprodione (Rovral; BayerCrop Science, Research Triangle Park, NC). Three replicate plates were used for each fungicide concentration and medium type, and the experiment was conducted twice.

    The isolates’ stability of resistance to fludioxonil and iprodione was determined by transferring them 22 times on unamended MEA. The 5-mm plugs were transferred from colony margins every 2 to 3 days. Resistance of isolates in storage was also analyzed for isolates incubated in the dark at 22°C for 3 months and transferred once each month. For both tests, after the final transfer, isolates were transferred to plates amended with 100 µg/ml of fludioxonil and 100 µg/ml of iprodione to detect any loss of resistance.

    DNA extraction, amplification, and sequencing Bos-1 and mrr1 genes.

    Because the Botrytis cinerea OS-1 (Bos-1) and the multidrug resistance regulator transcription factor (mrr1) genes have both been associated with fludioxonil resistance, their sequences were analyzed for potential resistance mutations in HR isolates. DNA was extracted according to the protocol by Chi et al. (2009).

    The complete Bos-1 gene was amplified with the primer pairs BF1/BR1, BF2/BR2, BF3/BR3, BF4/BR4, and BF5/BR5, described by Ma et al. (2007). Amplification was performed with the following parameters: initial denaturation at 95°C for 3 min, followed by 40 cycles at 94°C for 1 min, 60°C for 1 min, and 72°C for 1 min and 30 s, with final elongation at 72°C for 5 min.

    The majority of the mrr1 gene was amplified with primer pairs MRR1_atg/TF1-2-new and TF1-3/TF1-4 (Leroch et al. 2012). The amplification protocol consisted of an initial denaturation at 94°C for 5 min, followed by 40 cycles of 94°C for 30 s, 50°C for 1 min, and 72°C for 2 min, with a final elongation at 72°C for 10 min.

    PCR products were purified and sequenced by Arizona State University Core Laboratories. Resulting DNA sequences were aligned and translated, and then translations were aligned in SEQman Pro, SEQbuilder Pro (DNASTAR, Inc., Madison, WI), and MEGA7 (Kumar et al. 2016), respectively.

    Mycelium osmotic sensitivity determination.

    To determine whether fludioxonil resistance in HR isolates resulted in higher sensitivity to osmotic stress, S and HR isolates were grown on PDA for 3 to 5 days. Then 5 mm plugs were transferred from colony margins to PDA amended with no fungicide, 2% NaCl, and 4% NaCl, with three replicate plates for each concentration. After 2 days, relative growth was determined. The experiment was conducted twice.

    Mycelium intracellular glycerol accumulation assay.

    Because fludioxonil resistance has been associated with inhibition of the HOG-MAPK pathway and limitation of glycerol production, intracellular glycerol levels were determined in the presence and absence of fludioxonil for S and HR isolates. Mycelial hyphae were transferred on the tip of a toothpick from 1- to 2-week-old cultures to 50-ml Falcon tubes containing 20 ml potato dextrose broth. Three replicate tubes were used for each isolate and treatment. These liquid cultures were incubated at 22°C for 72 h and shaken at 100 rpm. At the end of this incubation period, fludioxonil was added to treatment tubes to a final concentration of 1 µg/ml, and tubes were incubated for an additional 4 h at 22°C and shaken at 100 rpm. Mycelium was then washed with sterile water, vacuum filtered, and weighed. Mycelium was then placed in 1 ml centrifuge tubes, and sterile water was added at a ratio of 0.1 g of mycelium per milliliter of water. After tissue was pulverized with an electric grinder, the mixture was centrifuged for 10 min at 12,000 rpm. Glycerol content in the supernatant was determined with the EnzyChrom Glycerol Assay Kit (BioAssay Systems, San Francisco, CA) according to manufacturer specifications.

    Pathogenicity and fungicide sensitivity detached flower assays.

    To determine pathogenicity of isolates and field relevance of fludioxonil resistance, detached flower assays were performed based on a protocol previously developed for petunia flowers (Bennett 2019). Twenty-three 5.7-liter pots with commercially grown calibrachoa plants, cultivar ‘MiniFamous Neo Royal Blue,’ were kept in a climate-controlled greenhouse in Clemson, South Carolina, with heating and cooling set points at 22 and 27°C, respectively, and retractable curtains providing 55% shade when external solar radiation exceeded 800 W/m2. No fungicides or insecticides were applied to plants or flowers for a month before experiments began. To ensure use of flowers at a similar growth stage and to decrease the likelihood of flower senescence during experimental trials, all open flowers were removed from plants before each trial. One to three days later, newly opened flowers were collected and placed in 7.6-cm floral tubes (1.27 cm diameter opening) filled to the brim with water. Treatment flowers were sprayed with fludioxonil (Scholar; Syngenta) at 1.5 mg a.i./ml, iprodione (Rovral; Bayer CropScience) at 1.5 mg a.i./ml, and chlorothalonil (Daconil; Syngenta) at 1.7 mg a.i./ml. These dosages represent the field rates recommended for gray mold control according to the labels of products registered for ornamental use: Medallion (fludioxonil; Syngenta), Chipco (Iprodione; Bayer), and Daconil (chlorothalonil; Syngenta). Sterile water was applied to untreated control flowers. Sprayed flowers in tubes were left overnight in a sterile hood, then placed in plastic boxes measuring 16.7 × 46.2 × 30.48 cm (L × W × D; Hefty, Reynolds Consumer Products, Lake Forest, IL) filled with sterilized sand to a depth of 15 cm. A total of 600 ml of water was added to the sand in each box just before the experiment to produce a humid chamber with approximately 100% relative humidity. Three floral tubes containing flowers were inserted into sand to a depth of approximately 3 cm, and three sprays (approximately 0.4 ml) of 2 to 3 × 105 spores/ml suspension were applied to each flower with a fine spray hand mister. Boxes were then sealed, randomized, and incubated in a growth chamber with temperature set to 22°C with 12 h light and dark cycles. Flower infection incidence was rated based on the presence or absence of disease symptoms, results were read after 2 days, and the entire experiment was performed twice.

    Fitness assessment.

    Fitness of HR isolates was assessed through mycelial growth and sporulation in vitro. In vitro mycelial growth was analyzed for HR and S isolates on both PDA and MEA media. Two perpendicular colony diameter measurements were taken for each plate after 2 days of growth. Three replicate plates were used for each isolate and medium type, and the experiment was conducted twice. One experimental replicate of S isolate FLOR5 had abnormally low growth and was therefore replaced with a replicate performed at a later date.

    For in vitro sporulation trials, colonies were grown at 22°C in 12 h light/dark cycles for 30 days, then spores were counted. To determine spore counts, all mycelium was scraped from the surface of plates, washed, and vortexed with 2 ml of water. Then 10 µl of spore suspension was removed and spores counted twice with a hemacytometer. Spore counts were performed for all isolates except isolate HP9, which did not sporulate in trials.

    To investigate the ability of isolates to sporulate on calibrachoa flowers, untreated flowers were inoculated with 0.4 ml of 2 to 3 × 105 spores/ml suspension in humid chambers. After 3 days of incubation, flowers were removed from humid chambers and placed in Petri dishes at 22°C in a light/dark incubator to decrease humidity and allow spore production. After 15 days, samples were examined for spore production under a microscope.

    Data analysis.

    Statistical models were developed to determine the relationships between the responses of mycelial growth, osmotic sensitivity, glycerol concentration, flower infection rating, and fitness parameters with the factors of experimental replicate, medium type, isolate, resistance, and their combinations. All models were full factorial general linear regression models. The data from two independent experiments for each experimental response were combined, and significant differences determined via analysis of variance. For mean separations, least squares mean differences Student’s t tests were used. All statistical analyses were performed with α = 0.05 and run in JMP Pro version 14 (SAS Institute Inc., Cary, NC).

    RESULTS

    Characterization of sensitivity to fludioxonil and iprodione.

    EC50 values for mycelial growth and germination.

    Both isolates exhibited high resistance to fludioxonil and iprodione in vitro, with spore germination and mycelial growth EC50 values >100 µg/ml. On both MEA and PDA media, germination was not significantly inhibited by fludioxonil or iprodione (P = 0.8647). No significant effects or interactions were detected in statistical analysis of spore germination.

    Mycelial growth of isolate 5.8.19 was not significantly inhibited by fludioxonil on either MEA (P = 0.1435) or PDA (P = 0.4615). Mycelial growth of isolate 5.10.19 was not significantly inhibited by fludioxonil on PDA (P = 0.3031). However, it was inhibited on MEA (P < 0.0001), with an average relative growth of 73%. Isolate 5.8.19 mycelial growth was inhibited by iprodione on PDA (P < 0.0001), with a relative growth of 91%. However, on MEA (P = 0.0066), isolate 5.8.19 grew more on iprodione than on the control. Isolate 5.10.19 mycelial growth was inhibited by iprodione on MEA (P < 0.0001), with a relative growth of 78% but not on PDA (P = 0.1738). Significant statistical effects of isolate (P < 0.0001), medium (P = 0.0023), and treatment (P = 0.0491) factors on the colony diameter response and significant interactions between isolate by medium by treatment (P < 0.0001); isolate by treatment (P < 0.0001); and medium by treatment (P = 0.0226) were observed.

    Stability of resistance.

    High resistance to both fludioxonil and iprodione was stable after 22 transfers on unamended MEA. On the last transfer, both HR isolates grew between 90 and 110% relative to the control on media amended with 100 µg/ml fludioxonil and 100 µg/ml iprodione, retaining their EC50 >100 µg/ml. EC50 values >100 µg/ml were observed in isolates transferred on unamended MEA for 3 months and for isolates stored in filter paper cultures in the −20°C freezer for 3 months.

    Bos-1 and mrr1 gene sequences.

    Sequence data, protein sequence.

    Mutations in the Botrytis two-component osmosensing histidine-kinase gene (Bos-1) have been shown to confer high resistance to both fludioxonil and iprodione in previous studies (Ren et al. 2016; Sang et al. 2018). The Bos-1 gene was sequenced for both HR isolates to detect potential fludioxonil resistance-causing mutations. The regions sequenced were 4,492 bp long and included the entire coding sequence of the gene (accession nos. MT497335 and MT497336). The sequences for both isolates contained previously reported iprodione resistance mutations F127S (ATC → AAC) and I365N (TTT → TCT). However, isolate 5.8.19 additionally contained the mutation L267V between HAMP domains 1 and 2 (Table 2).

    A variety of mutations in the multidrug resistance regulator 1 gene (mrr1) can result in low to moderate fludioxonil resistance but not cross-resistance to iprodione (Table 1). To determine whether these mutations were present, the majority of the mrr1 coding region was sequenced for HR isolates 5.8.19 and 5.10.19. The resulting nucleotide sequences were 2,400 bp long and were identical for both isolates (accession nos. MT497337 and MT497338). These sequences were 100% identical to previously characterized mrr1 sequences of a representative group S isolate from Germany, accession no. JX266765 (Leroch et al. 2012), and contained the 3-bp deletion ΔL497, associated with the MDR1h phenotype (Table 1).

    Osmotic sensitivity.

    Because resistance associated with the MAPK pathway generally confers increased osmotic sensitivity, the sensitivity of HR isolates to salt stress was investigated on 2 and 4% NaCl media. Growth of both HR isolates was significantly inhibited compared with sensitive isolates by 2% NaCl on both MEA (P < 0.0001) and PDA (P < 0.0001) and completely inhibited by 4% NaCl on both MEA and PDA (Fig. 1). Resistant isolates were consistently inhibited on both media, although some differences between media were observed for sensitive isolates. (Fig. 1). Significant effects of experimental replicate (P = 0.0328), isolate (P < 0.0001), medium (P < 0.0001), and treatment (P < 0.0001) on the relative growth response were detected. Additionally, significant interactions were detected between experimental replicate and isolate (P < 0.0001); experimental replicate and medium (P = 0.0011); isolate and medium (P < 0.0001); experimental replicate, isolate, and medium (P = 0.0027); and isolate and treatment (P < 0.0001). The significant effect of experimental replicate (P = 0.0328) was detected because sensitive isolate FLOR5 grew inconsistently between replications on MEA. However, its growth was consistently significantly greater than that of resistant isolates, so replications were still combined.

    Fig. 1.

    Fig. 1. Osmotic sensitivity of isolates highly resistant and sensitive to fludioxonil on 2 and 4% NaCl on A, potato dextrose agar and B, malt extract agar media. Statistical comparisons were performed separately for 2 and 4% salt treatments, as indicated by upper and lowercase letters.

    Download as PowerPoint

    Intracellular glycerol accumulation.

    Because fludioxonil resistance associated with the MAPK pathway typically prevents accumulation of intracellular glycerol in the presence of fungicide, glycerol content of HR and S isolates was analyzed in the presence and absence of fludioxonil. In the absence of fludioxonil, all isolates had low average glycerol concentrations, ranging between 0.05 and 0.13 mM. Neither resistant isolate significantly increased production of glycerol in the presence of fludioxonil (P = 0.3658). However, all three sensitive isolates significantly increased glycerol production in the presence of fludioxonil (P < 0.0001) (Fig. 2). Significant effects of isolate (P < 0.0001), treatment (P < 0.0001), and experimental replicate (P = 0.0003) factors were observed on the response of glycerol concentration. Additionally, significant interactions between isolate and treatment (P < 0.0001); isolate and experimental replicate (P = 0.0004); treatment and experimental replicate (P = 0.0108); and isolate, treatment, and experimental replicate (P = 0.0171) were detected. The significant effect of experimental replicate (P = 0.00027) was detected because of variation in glycerol accumulation between replicates of the S isolate HP9 treated with fludioxonil. However, fludioxonil-treated replications were always significantly greater than the untreated replications for this isolate, so replications were still combined for analyses.

    Fig. 2.

    Fig. 2. Glycerol concentration in mycelium of isolates treated and untreated with fludioxonil, with statistical comparisons across all isolates and treatments. Letters designating statistical differences compare glycerol concentrations within but not between isolates.

    Download as PowerPoint

    Pathogenicity and fungicide resistance ex vivo on detached flowers.

    To assess practical relevance of high fludioxonil resistance, detached flower assays were performed, and isolates were exposed to recommended label rates of fludioxonil and iprodione. No fungal growth or symptoms were observed on any uninoculated flowers. On flowers not treated with fungicide, the percentage of infected flowers ranged from 67 to 100%, with an average of 94.5%. No isolate was consistently less pathogenic than other isolates. After treatment with label rates of fludioxonil and iprodione, sensitive isolate FLOR5 did not form any disease symptoms on flowers. Highly resistant isolates 5.8.19 and 5.10.19 caused significantly more or the same amount of disease as on flowers not treated with fungicide (P = 0.00238). Chlorothalonil generally reduced disease incidence significantly compared with the water-treated control or the resistant isolates treated with fungicide, though not as dramatically as fludioxonil and iprodione. Exceptions were a single repetition of isolate 5.8.19 treated with chlorothalonil that did not significantly reduce incidence and individual replicates of isolates 5.8.19 and FLOR5 that chlorothalonil inhibited to the same extent as fludioxonil and iprodione (P = 0.1798). (Fig. 3, Supplementary Fig. S1). Significant effects of the factors treatment (P < 0.0001) and isolate (P < 0.0001) on presence or absence of infection symptoms were detected, as well as interactions between experimental replicate and isolate (P = 0.0147); treatment and isolate (P < 0.0001); and experimental replicate, treatment, and isolate (P = 0.0024).

    Fig. 3.

    Fig. 3. Disease incidence on uninoculated flowers and flowers inoculated with highly resistant isolates 5.8.19, 5.10.19, and sensitive isolate FLOR5. Letters designating statistical differences compare flower disease incidence within but not between treatments.

    Download as PowerPoint

    Fitness of isolates resistant and sensitive to fludioxonil.

    To determine whether high resistance to fludioxonil and iprodione resulted in fitness penalties, growth rate and spore production were analyzed for S and HR isolates. In both experimental replicates, HR isolates grew significantly more slowly than sensitive isolates on both MEA (P < 0.0001) and PDA (P < 0.0001) after 2 days of growth. Isolate 5.8.19 had significantly larger colony diameter than 5.10.19 on MEA (P < 0.0001), but the difference between the two isolates was insignificant on PDA (P = 0.0534), although isolate 5.8.19’s average diameter was numerically larger for both media (Table 3). Significant effects of the factors of experimental replicate (P < 0.0001), isolate (P < 0.0001), and medium (P < 0.0001) on colony diameter were detected, as were significant interactions between experimental replicate and isolate (P < 0.0001); isolate and medium (P < 0.0001); and experimental replicate, isolate, and medium (P = 0.370). The significant effect of experimental replicate (P < 0.0001) was detected through variations between S isolate diameters between replications. However, on both media types, the sensitive isolates consistently grew significantly more than resistant isolates, so replications were combined.

    TABLE 3. Fitness components of isolates sensitive and highly resistant to fludioxonily

    HR and S isolates were able to produce spores both in vitro and ex vivo on detached flowers. In vitro, sensitive isolates NC4 and FLOR5 produced significantly more spores than HR isolates 5.8.19 and 5.10.19 on both media types (P < 0.0001). No significant difference was detected in sporulation between isolates 5.8.19 and 5.10.19 (P = 0.8563) or between isolates NC4 and FLOR5 (P = 0.1543) on either medium. Significant effects of the factor of isolate (P < 0.0001) on spore count were detected, as were interactions between assay and medium (P = 0.0003) and assay, isolate, and medium (P = 0.0072).

    DISCUSSION

    Although MDR-based low and moderate resistance to fludioxonil have occurred in fields for years, detection of HR B. cinerea isolates is unusual and potentially concerning. Fludioxonil has had an unprecedented history of limited resistance development and is often the primary single-site inhibitor used to control emergency outbreaks in commercial production of strawberries and other crops (Kilani and Fillinger 2016). However, this article is the fourth report of high resistance to fludioxonil in greenhouse production worldwide in 5 years (Ren et al. 2016; Sang et al. 2018; Zhou et al. 2020) (Table 1).

    High resistance to fludioxonil has been consistently linked to increased osmotic sensitivity, lack of glycerol accumulation in the presence of fludioxonil, and (in most cases) mutations in the HAMP domains of the Bos-1 gene associated with iprodione or fludioxonil resistance (Fillinger et al. 2012; Ren et al. 2016; Sang et al. 2018; Zhou et al. 2020). Our isolates followed similar trends, lacking fludioxonil-induced glycerol accumulation and exhibiting increased osmotic sensitivity. Both HR isolates exhibited iprodione resistance ex vivo and in vitro and contained the documented dicarboximide resistance mutations F127S and I365N (Fernández-Ortuño et al. 2015; Grabke et al. 2014; Ma et al. 2007). Isolate 5.8.19 also contained the mutation L267V between HAMP domains 1 and 2, near the mutations G262S and G265D, associated with fludioxonil resistance by Ren et al. (2016) in greenhouse isolates. A knockout and molecular modeling study revealed that the nearby G278D mutation in the Bos-1 gene of a lab mutant conferred high fludioxonil and iprodione resistance and osmotic sensitivity (Fillinger et al. 2012) (Table 1).

    Isolate 5.10.19 did not contain mutations other than F127S and I365N in the Bos-1 gene. HR to fludioxonil and osmotic sensitivity in the absence of Bos-1 mutations has been reported before for a laboratory mutant but never for a greenhouse or field isolate (Zhou et al. 2020). Other genes such as Sak1 (Segmüller et al. 2007), BcOS4 (Yang et al. 2012), bos5 (Yan et al. 2010), and BRRG-1 (Yan et al. 2011), associated with the HOG1 MAPK pathway in B. cinerea, have been reported, but knockouts of these genes resulted in unchanged or increased sensitivity to fludioxonil (Segmüller et al. 2007; Yan et al. 2011, 2010; Yang et al. 2012), and isolate 5.10.19 exhibits HR to fludioxonil. Previous literature has suggested that Bos-1 regulates the gene Sak-1 but that a separate pathway is responsible for dicarboximide and phenylpyrrole resistance (Liu et al. 2008). It is possible that malfunction of a gene in this separate pathway may be responsible for resistance in isolate 5.10.19 (Liu et al. 2008).

    Although HR isolates 5.8.19 and 5.10.19 exhibited the 3-bp deletion ΔL/V497 in the gene mrr1, conferring the MDR1h phenotype, it is unlikely that the MDR1h phenotype is responsible for the fludioxonil HR observed. The highest fludioxonil EC50 values for B. cinerea isolates with the MDR1h phenotype are typically about 3.1, and HR (EC50 > 5 µg/ml) to fludioxonil caused by atrB efflux has never been reported for B. cinerea (Fernández-Ortuño et al. 2015; Fillinger et al. 2012; Hu et al. 2019) (Table 1). Also, our HR isolates exhibited changes in glycerol accumulation and osmotic sensitivity indicative of HOG1 MAPK interruption. Osmotic sensitivity has been reported for MDR1h isolates, but inhibition on 4% NaCl PDA media was <50% compared with the complete inhibition observed for our HR isolates (Fernández-Ortuño et al. 2015). Additionally, mycelial growth (Supplementary Fig. S2) and sporulation of our HR isolates were significantly less than those of sensitive isolates, and the literature reports that MDR1h isolates do not significantly differ from sensitive isolates in these traits when grown at 22°C (Fernández-Ortuño et al. 2015). It is also unlikely that mycelial growth and osmotic fitness penalties are caused by resistance to other chemical classes of fungicides. Chen et al. (2016) previously showed that mycelial growth of B. cinerea isolates with resistance to five or six chemical classes at the same time was not significantly different from that of fully sensitive controls. Additionally, isolates with resistance to five or six chemical classes were only slightly more sensitive to osmotic pressure than fully sensitive isolates at 22°C (Chen et al. 2016).

    Although recent reports highlight the presence of highly resistant fludioxonil isolates in China, and this study shows that HR isolates are also present in the United States, their importance to commercial production is uncertain (Ren et al. 2016; Sang et al. 2018; Zhou et al. 2020). HR isolates consistently exhibit lower growth rates, sporulation, and greater osmotic sensitivity than sensitive isolates. Also, HR isolates have never been reported outside greenhouse conditions, where high humidity and plant concentration provide the ideal atmosphere for resistance development and spread. Even within greenhouses, outbreaks of HR isolates appear to be well contained, making up small percentages of the population and not resulting in complete loss of fludioxonil efficacy in the locations where resistance has been reported (Ren et al. 2016). Therefore, fitness penalties may prevent these isolates from causing imminent widespread crop loss through lack of efficacy. However, these findings indicate the need to use fludioxonil with caution and use resistance management techniques.

    ACKNOWLEDGMENTS

    We thank Melissa Muñoz and William Bridges for their assistance with statistical analysis, and P. K. Bryson and Sydney Lykins (Clemson University) for technical support.

    The author(s) declare no conflict of interest.

    LITERATURE CITED

    The author(s) declare no conflict of interest.

    Funding: Technical contribution number 6899 of the Clemson University Experiment Station. This material is based on work supported by NIFA/USDA, under project number SC-1700530.