Next Article in Journal
Fungal Endophytes: A Potential Source of Antibacterial Compounds
Previous Article in Journal
Genomic Analysis of Stropharia rugosoannulata Reveals Its Nutritional Strategy and Application Potential in Bioremediation
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JoF
Volume 8
Issue 2
10.3390/jof8020163
jof-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Detection and Molecular Phylogenetic-Morphometric Characterization of Rhizoctonia tuliparum, Causal Agent of Gray Bulb Rot of Tulips and Bulbous Iris

by
Katie Coats
1,
Annie DeBauw
1,
Dilip K. Lakshman
2,*,
Daniel P. Roberts
2,
Adnan Ismaiel
2 and
Gary Chastagner
1,*
1
Department of Plant Pathology, Washington State University Research & Extension Center, 2606 West Pioneer, Puyallup, WA 98371, USA
2
Sustainable Agricultural Systems Laboratory, USDA-ARS, 10300 Baltimore Avenue, Beltsville, MD 20705, USA
*
Authors to whom correspondence should be addressed.
J. Fungi 2022, 8(2), 163; https://doi.org/10.3390/jof8020163
Submission received: 18 January 2022 / Revised: 31 January 2022 / Accepted: 3 February 2022 / Published: 8 February 2022
(This article belongs to the Section Fungal Genomics, Genetics and Molecular Biology)
Browse Figures
Review Reports Versions Notes

Abstract

:
Gray bulb rot of tulips and bulbous iris is caused by the soil-borne fungal pathogen, Rhizoctonia tuliparum (Rtul). Sclerotia present in infected bulbs, as well as overwintering sclerotia in soil and field debris, are the primary sources of infection. A method for accurate and sensitive detection of Rtul from soil and infected bulbs, and estimation of inoculum threshold levels, is needed for the management of disease caused by this pathogen. We designed a unique set of primers targeting the ITS2 region of the Rtul genome and developed a highly sensitive quantitative PCR (qPCR)-based method for Rtul identification using these primers, where the threshold of detection was approximately 1 fg Rtul DNA. The assay was more sensitive with sclerotia collected from the field (natural) than with those grown in the lab, and more sensitive with natural-light than natural-dark sclerotia. Also, the detection method was more sensitive when sclerotia were extracted from soil than from bulb tissue. The qPCR method was highly specific, as no PCR amplification was detected when genomic DNA from 62 non-Rtul Rhizoctonia isolates from a wide range of anastomosis groups were tested. To understand the evolutionary relationships and genomic diversity of Rtul, we performed phylogenetics of the ITS1-5.8S-ITS2 region and ITS2-molecular morphometric characterization (MMC) of Rtul isolates. The three Rtul isolates whose ITS sequences were available in GenBank formed a distinct phylogenetic clade with Ceratobasidium anceps as the nearest relative. Furthermore, MMC analysis revealed genetic divergence among these three Rtul isolates.

1. Introduction

Gray bulb rot of tulips and bulbous iris is caused by the soil-borne fungal pathogen Rhizoctonia tuliparum (Rtul) [Whetzel and Arthur; Synonym: Sclerotium tuliparum (Kleb)] [1,2,3,4,5,6]. The disease was initially described by Wakker [7] in the Netherlands, and the pathogen was described by Whetzel and Arthur [8]. Gray bulb rot is primarily a problem in field or landscape-grown tulips and irises, but the pathogen has also been reported in daffodil, crocus, gladiolus, and lily [9]. The fungus forms felty masses between the scales, and infected bulbs rot from the top down. Both the bulbs and roots turn grayish in color and progressively dry and shrivel at the site of infection. Infected shoots either fail to emerge or die shortly after infection during late winter. If a plant does grow, it is usually slow to develop and withers and dies before flowering [9]. Diseased plants tend to occur in patches in field and landscape plantings, but this disease can completely destroy the crop after several rotations if unmanaged.
Other than volunteer plants, overwintering sclerotia of Rtul residing in soil and host debris are important sources of inoculum. Sclerotia range in size from 1–10 mm in diameter and are formed on infected plant material. They turn from white to dark brown or black as they ripen and can be found on bulbs and in surrounding soil. According to Coley-Smith et al. [10], about 10% of the sclerotia can remain viable in soil for at least ten years. Sclerotia usually remain attached to soil and plant debris when the bulbs are dug, allowing them to be easily spread from field to field by cultivating equipment. The pathogen can also be spread by the movement of infected planting stock.
Currently, there are no reliable and rapid assay methods to assess inoculum levels of Rtul in the soil and on infected bulbs [11,12,13]. The ability to quantify inoculum levels of Rtul in the soil and on planting stock, coupled with information on inoculum threshold levels necessary for disease development, may allow for the development of predictive diagnostic tests to identify high-risk fields and planting stocks where inoculum of this pathogen is above threshold values [14,15,16,17]. Access to this type of information would assist growers in making management decisions relating to the application of chemical or other control measures and potentially reduce the use of fungicides in the production of tulips and iris [13]. Therefore, we developed a sensitive and specific quantitative PCR (qPCR) assay [18] to detect Rtul and tested its effectiveness with bulb and soil samples.
The nuclear status and the teleomorph of the anamorphic fungus Rtul have not been unambiguously resolved; isolates from the UK [19] and Netherlands [20] having been reported to be binucleate, potentially placing Rtul as Ceratobasidium sp. Therefore, we determined the phylogenetic status of Rtul within the Ceratobasidiaceae using the ribosomal internal transcribed spacer (ITS1-5.8S-ITS2). To gain information on the evolutionary relationships and genomic diversity of Rtul, we also characterized the ITS2-primary and secondary structure-based molecular morphometrics (MMC) of isolates currently available in GenBank [21,22].

2. Materials and Methods

2.1. Fungal isolates

The fungal isolates used in this study are listed in Table 1. Isolate I-399 of Rtul was obtained in 1998 from infested soil from a commercial bulbous iris field near Mount Vernon, WA. The other Rtul isolates were collected in 2008 and 2012 from symptomatic bulbous irises grown near Puyallup, WA. Agar plugs of Rtul were retrieved from long-term storage at 4 °C in water, transferred to potato dextrose agar (PDA), and cultures of Rtul maintained on PDA at room temperature prior to DNA extraction. The sixty-two non-Rtul isolates used for specificity testing in qPCR assays were stored on filter paper and held at room temperature prior to DNA extraction.

2.2. DNA Extraction from Fungal Isolates, Soil, and Plant Tissue

DNA was extracted from homogenized fungal isolates using a Gentra Puregene Tissue Kit (QIAGEN, Valencia, CA, USA), following the manufacturer’s instructions. Sclerotia and/or mycelia were homogenized (FastPrep-24 Homogenizer, MP Biomedicals, Santa Ana, CA, USA) at maximum speed for one min in 2 mL microcentrifuge tubes containing six 3.5-mm-diameter glass beads (BioSpec Products, Bartlesville, OK, USA).
Soil samples were collected from the top 10 cm of a Puyallup fine sandy loam field site at the Washington State University Research and Extension Center in Puyallup, air dried, passed through a 2-mm sieve, and stored at room temperature. Soil was autoclaved at 125 °C for 40 min prior to DNA extraction. DNA was extracted using the Wizard Magnetic DNA Purification System for Food (Cat. No. PRFF3750; Promega, Madison, WI, USA) with Lysis Buffer A replaced by Soil Extraction Buffer (SEB; 120 mM sodium phosphate buffer, pH 8, 2% CTAB, 1.5 M NaCl, 2% antifoam B emulsion). Soil samples were mixed to homogeneity and triplicate 4-g samples weighed and added to 50 mL tubes containing 16 mL SEB and 5 mL 3.5-mm-diameter glass beads. Samples were then homogenized at maximum speed for one min and centrifuged at 2000× g for three min. Supernatant (500 µL) was transferred to 2 mL microcentrifuge tubes containing 250 µL Buffer B (Cat. No. Z3191, Promega) using 200 mg of starting material (manufacturer’s instructions were followed starting from Step 3). DNA samples were stored at 4 °C until qPCR analysis.
Tulip bulbs (Tulipa gesneriana var. Advent) from a Washington farm were harvested in 2013. Bulbs were surface sterilized with a 0.825% hypochlorite solution, rinsed twice with water, and homogenized to a textured paste consistency. Tulip tissue was kept frozen at –20 °C prior to DNA extraction. DNA was extracted with the Wizard Magnetic DNA Purification System for Food with SEB modification, as described above.

2.3. Primer Design and qPCR Assay for R. tuliparum

Sequence of the ribosomal DNA ITS region of the Rtul isolates (GenBank accessions KX767078, KX767079, KX767080) were compared among themselves and to a variety of Rhizoctonia sequences available in GenBank. Several sequence spans were identified that were conserved within the Rtul isolates but variable in other Rhizoctonia species. Three sets of primers and TaqMan® MGB probes potentially specific to Rtul were designed using Primer Express Software Version 3.0 (Applied Biosystems, Foster City, CA, USA) and tested. The primer/probe set selected for further experimentation was designated as “Rtul” and amplified a 100-bp fragment of ITS (Table 2; Figure S1). The Rtul hybridization probe was labeled at the 5′ end with tetrachlorofluorescein (TET™) and at the 3′ end with a minor groove binder nonfluorescent quencher (MGB-NFQ) following the manufacturer’s instructions (Applied Biosystems, Bedford, MA, USA). qPCR was carried out at the default cycling parameters of 50 °C for 2 min, 95 °C for 10 min, 40 cycles at 95 °C for 15 s, and 60 °C for 1 min on a 7500 Real Time PCR System (Applied Biosystems) using Sequence Detection Software Version 1.3.1. An exogenous amplification control (EAC) of salmon testes DNA spiked into each reaction to measure potential amplification inhibition was detected by Sketa primers and a Sketa TaqMan® probe [23] labeled at the 5′ end with fluorescein (6-FAM) dye and at the 3′ end with a tetramethylrhodamine (TAMRA™) quencher (Thermo Fisher Scientific, Waltham, MA, USA). The qPCR mixture (25 μL) included 12.5 µL TaqMan® Universal Master Mix II no AmpErase UNG; 100 nM each Rtul primers, Rtul probe, EAC Sketa primers and EAC Sketa probe; and 0.04 ng µL−1 salmon testes DNA (Cat. No. D7656, Sigma Chemical Co., St. Louis, MO, USA). Standard curves were generated based on linear regression of Ct values produced by a dilution series of I-399 Rtul DNA containing 7.7 × 10−2 to 7.7 × 10−6 ng DNA. Rtul DNA quantity in unknown samples was calculated by applying their cycle threshold (Ct) value as the independent variable in the resulting linear regression equation. Sterile water in place of the DNA template was run as the negative control. To determine the lowest threshold of detection with the Rtul primer/probe set, the concentration of DNA sample extracted from a pure culture of Rtul (DGF2-VV) was estimated with a fluorometer (Thermo Fisher Scientific) and used to make a 1:10 serial dilution series from 1.2 ng μL−1–0.12 fg μL−1 of DNA in the Rtul qPCR assay (Figure 1). The amplification efficiency (E) was calculated using the formula: E = 10(−1/slope) – 1 (Ruijter et al. [24]). DNA samples, including standards and controls, were assayed in triplicate (n = 3).

2.4. Detection of Lab-Grown and Field-Collected (Natural) Sclerotia of R. tuliparum

Lab-grown sclerotia were prepared by transferring small cuttings of I-399 Rtul culture, containing mycelia and sclerotia, onto 10-cm-diameter PDA plates and incubating at 23 °C in the dark. Sclerotia were removed from cultures on PDA plates with a sterile pipet tip after 10 weeks and stored at room temperature. Natural sclerotia were obtained by manually sieving air-dried field soil (Puyallup fine sandy loam) which was collected 9 months after iris bulbs were planted in a field site that had been infested with Rtul inoculum on rice grains. Sclerotia were then stored in a vial at ambient temperature. Natural sclerotia ranged from brown to black in coloration and were categorized as light-colored or dark-colored (Figure S2). Due to the variability of sclerotia size and weight, 189 lab-grown sclerotia were weighed to determine the average weight of one sclerotium to be 1.78 mg. The amount of sclerotial material used in each experiment was a weight equivalent based on this number. DNA extracted from six replicates of eight (14.25 mg) lab-grown, natural-light, and natural-dark sclerotia were analyzed in triplicate by Rtul qPCR undiluted and diluted with sterile water 1:10, 1:100, and 1:1000. Experiments were repeated and combined prior to analysis as there were no significant differences in results.

2.5. Detection of R. tuliparum in Soil and Plant Samples

For detection of Rtul in soil, the weight equivalent of one, two, four, eight, 16 and 32 lab-grown sclerotia were added to 4 g soil (Puyallup fine sandy loam field soil) samples in triplicate, prior to DNA extraction. DNA samples were diluted 1:1000 and analyzed with the Rtul/Sketa qPCR method described above in triplicate and the entire test repeated and results combined prior to analysis.
For the detection of Rtul in tulip tissue, the weight equivalent of eight lab-grown sclerotia was added to the weight equivalent of 1/8, 1/4, and 1/2 tulip bulb in triplicate, and DNA was extracted as described above. The tunic was removed from thirty bulbs prior to determining mean bulb weight to be 20.78 g. DNA samples were analyzed by Rtul/Sketa qPCR and compared to the equivalent treatment of sclerotia in the absence of bulb tissue.

2.6. Phylogenetic Analyses

The Rtul ITS sequences of our three GenBank submissions (GenBank accessions KX767078, KX767079, KX767080) were found to be identical through alignment and henceforth will be considered as the same haplotype, accession KX767078 being the largest of the sequences. Accession KX767078 was subjected to a BLAST search using the National Center for Bioinformatics Institute (NCBI) website. Sequences with high homology to KX767078, along with ITS sequences from known species of Ceratobasidum, Rhizoctonia solani (Thanatephorus cucumeris), and Waitea circinata were obtained in FASTA format, aligned using Clustal X [25], and manually adjusted using Mesquite editing software [26]. The final alignment file contained sequence from 34 isolates. Phylogenetic trees were constructed using Bayesian inference (BI) and maximum parsimony (MP) analyses. The Bayesian tree was obtained using MrBayes, Version 3.2.7 [27], using the general time reversible (GTR + I + G, nst = 6) model with gamma distribution and invariable sites with default priors. Four chains and one million generations were run with a tree selected every 1000 generations. The initial 25% of trees were discarded during the burn-in phase. The remaining trees were used to produce a 50% majority rule consensus tree. This tree was edited and posterior probabilities above 0.95 placed above the branches. The MP tree was obtained with PAUP Version b10 software [28] using a heuristic search with a starting tree obtained via 1000 random stepwise addition sequences, and tree-bisection-reconnection as the branch-swapping algorithm, with MULTrees in effect. Bootstrap values were calculated with 1000 replicates. For both the Bayesian and MP trees, W. circinata var. agrostis was used as an outgroup based on the distance of the isolate from the rest of the isolates in the tree.

2.7. Morphometric Analyses of the ITS2 Region

For ITS2 structure analysis, the ITS1-5.8S-ITS2 sequences were retrieved for C. anceps isolates CBS 152.32 (GenBank accession MH855251) and CBS 152.32 (AJ427402), and three Rtul isolates: isolate I-399 (KX767078), the European isolate 29,792 (EU191041), and isolate from New York CBS206.25 (MH854847) [29]. The corresponding ITS2 sequences were identified utilizing the annotation utility of the ITS2 database [30]. Secondary structure was predicted with the prediction utility of the ITS2 database [30] and manually analyzed for morphometric features like compensatory base change (CBC), hemi-compensatory base change (hCBC), indels of base pairing, pseudoknot formation, and conserved loop sequences. Percent GC was estimated utilizing tools from PAUP Version b10 software [28].

3. Results

3.1. R. tulipia qPCR Detection Assay Sensitivity and Specificity

Analysis of the serial dilution standard curve showed a linear relationship (R2 = 0.994) between the log of the DNA concentration of Rtul and the cycle threshold (Ct) value over the seven most concentrated DNAs (i.e., 7.7 × 10−2 to 7.7 × 10−6 ng) (Figure 1). The least concentrated DNA sample used in the standard curve, at 0.12 fg μL−1, was not detected and most likely indicates that the threshold of detection for this assay was approximately 1 fg Rtul DNA. For determination of assay specificity, 62 non-Rtul and six Rtul isolates of Rhizoctonia (Table 1) were analyzed to ensure that amplification occurs for all Rtul isolates and does not occur for any non-Rtul isolate. Tested isolates were from R. solani anastomosis groups (AGs) 2-1, 3, 4, 5, 8, and 9, and Ceratobasidium isolates from AGs A, E, and I, as well as W. circinata and several unidentified Rhizoctonia species. No amplification was observed with the 62 non-Rtul DNA samples, and Ct values between 13.93 and 33.51 were observed for various Rtul DNA samples (Table 1).

3.2. Detection of Lab-Grown and Natural Sclerotia

Analysis of dilutions of DNA extracted from three types of sclerotia (lab-grown, natural-dark, and natural-light) revealed that the Rtul detection assay was more sensitive, in that it more efficiently detected Rtul DNA, when the DNA was diluted in the range of 1:10 to 1:1000 before qPCR (Figure 2). The Sketa internal amplification control confirmed that amplification was suppressed at higher DNA concentration levels. Three replicates each of DNA extracted from the weight equivalent of eight lab-grown, natural-dark, and natural-light sclerotia were diluted 1:1000 and analyzed by the Rtul qPCR assay to determine which type of sclerotia yielded more Rtul DNA (Figure 3). The natural-light sclerotia yielded an average of two and seven times more DNA than the natural-dark and the lab grown sclerotia, respectively.

3.3. Detection of R. tuliparum in Soil

Since sclerotia overwinter in soil, the utility of the qPCR method was evaluated by determining the amount of Rtul DNA from various numbers of sclerotia extracted from fixed quantities of soil. Weight equivalents of one, two, four, eight, and 16 lab-grown sclerotia did not yield significantly different Rtul qPCR results when extracted alone or in 4 g soil (Figure 4). In contrast, the weight equivalent of 32 sclerotia was detected more efficiently in soil; about 62% more Rtul sclerotial DNA was detected in soil compared to pure sclerotia not mixed with anything. The detection of lab-grown sclerotia became unreliable when the number of sclerotia per gram of soil was less than two.

3.4. Detection of R. tuliparum in Bulbs

To determine the detection limit of the Rtul assay using infected bulb tissue, an increasing amount of tulip bulb tissue was extracted with eight lab-grown Rtul sclerotia. An inverse sensitivity of detection of Rtul was observed when sclerotial DNA was extracted in the presence of increasing amounts of bulb tissue. At the maximum amount tested, one-half a tulip bulb, detection of DNA from eight lab-grown sclerotia decreased by an average of 84% (calculated as 100 – [100 × 0.45 bulb fraction/2.8 ng recovered DNA]), indicating possible carryover of PCR inhibitors in sample preparations in the presence of bulb tissues (Figure 5).

3.5. Phylogenetic Analyses

The Bayesian and MP trees were essentially identical in topology; therefore, we only present the Bayesian tree. The bootstrap support values of 75% or higher from the MP analyses are indicated on the Bayesian tree with the symbol (*) under the branches (Figure 6). The Bayesian tree showed that all three geographically separated isolates of Rtul [KX767078 (Washington, DC, USA), EU191041 (Europe) and MH854847 (Ithaca, NY, USA)] formed a monophyletic group with high posterior probability and bootstrap values. The clade had a strong sister relationship with the clade of two C. anceps isolates, with posterior probability of 0.99 and bootstrap value of >70%. No further relationship can be inferred from the tree due to lack of phylogenetic support with other clades.

3.6. Morphometric Analyses of the ITS2 Region

The ITS2 sequence length was 237 bases for C. anceps, and 238 bases for each of the three Rtul accessions. The GC content of the ITS2 region was 42.6% for C. anceps, and 46.2%, 45.8%, and 45.8% for the isolates of Rtul with GenBank Accessions KX767078, EU191041, and MH854847, respectively (Table S1).
In comparison with the secondary structure of the ITS2 region of C. anceps isolate CBS 152.32 (MH855251), the ITS2 secondary structure of Rtul isolate I-399 (KX767078) had a CBC and a hCBC in helix II, lost three base pairs in helix III, and lost one base pair in helix IV (Figure 7; Table S2). The presence of the CBC, hCBC, and a loss/gain of base pairs amongst the two species is commensurate with two distinct biological species [31]. This point is supported by the fact that they formed two distinct phylogenetic clades with strong bootstrap value at the primary ITS2 sequence level (Figure 6). However, phylogenetically, the two species are closer to one another than the many other species of Ceratobasidium. The C. anceps and the three isolates of Rtul have maintained a conserved AAGT in the terminal loop of helix III. The conserved triplets of AAT in helix I, CAC in helix II, and TAG in helix IV were observed in respective ITS2 terminal loops of C. anceps and the three Rtul isolates. (Figure 7; Table S3).
In comparison with the ITS2 secondary structure of the Rtul Washington isolate I-399 (KX767078), the corresponding structures of European isolate 29,792 (EU191041) only have a point mutation in the loop region (transversion, C → A) of helix III, resulting in the gain of a base pair (A-T) (Figure 8; Table S4). However, in comparison with the ITS2 secondary structure of isolate I-399, the ITS2 secondary structure of isolate CBS206.25 (MH854847) showed structural changes in helices II, III, and IV. Thus, a G→C transversion in helix II resulted in the loss of a G-C base pairing in helix II. In helix III, (i) a transversion in the loop region (C→A) resulted in the gain of a base pair (A-T), (ii) an A→C transversion resulted in the loss of a base pairing (A-T) needed for a single base pair pseudoknot-like structure, and (iii) a C→G transversion resulted in the net loss of a C-G base pair in helix III. Therefore, a net loss of one base pair was observed in helix III of isolate CBS206.25. A missing pseudoknot in helix III of the CBS206.25 isolate of Rtul may be a significant genetic shift. Moreover, a G→A transition in helix IV of the latter isolate resulted in the loss of a C-G base pair.
We also noted that the single base pair pseudoknot-like structure was formed in helix III of Ceratobasidium AGs A, Ba, Bb, B(o), C, G, H, I, L, Q, uninucleate Rhizoctonia 1(UNR1) isolates (Table S5), C. anceps, C. fluccosa, R. alpine, R. butini, C. pseudocornigerum, and C. noxium isolates which phylogenetically belonged in the extended clade containing the Rtul isolates (Figure 6). Interestingly, a pseudoknot-like structure did not form in the ITS2 helix III from AG groups E, Fa, Fb, K, P, O, R, S, U, and UNR2 accessions AF200515 and AF200517 (Table S5).

4. Discussion

We developed a robust TaqMan® real-time qPCR assay to specifically detect and quantify Rtul DNA based on a unique sequence of the internal transcribed spacer region. Specificity to Rtul with the exclusion of other Rhizoctonia genera, species, AGs, and subgroups was established by testing the assay on nontarget isolates, of which none amplified. The Rtul PCR primers and the probe exactly matched the Rtul ITS2 sequence. Moreover, the NCBI BLAST search (dated 14 January 2022) of the Rtul probe showed 100% sequence homology with ITS sequence region of C. pseudocornigerum (MH861653), C. angustisporum (NR_154601), and several species of mycorrhizal Ceratobasidium, mostly from orchid hosts (Table S6). The Rtul probe also matched 13 out of 17 bases (data not shown) with its close phylogenic relative, C. anceps (Figure 6), which is a pathogen of persimmon (Diospyros kaki) and tea [32]. Because of host differences, the presence of above pathogens in proximity to Rtul pathogen would be remote. Therefore, the qPCR assay should be considered specific for Rtul when tested with soils from iris fields. All Rtul samples tested were amplified, and a dilution series of Rtul DNA revealed a limit of detection of 1 fg DNA. The Rtul qPCR assay presents a signification advance in molecular diagnostics of the pathogen and is appropriate for identification of isolates or sclerotia, as well as detection of the pathogen in soil or plant tissue.
Most Rhizoctonia pathogenicity studies use inoculum grown on culture plates, which produce sclerotia which are slightly different in appearance from those naturally occurring in soil. In our observation, lab-grown sclerotia are brown and corky in texture, while natural sclerotia have a darker and harder outer rind. In addition, natural sclerotia range from light brown to black in coloration. Because the intended application of the qPCR assay is to detect natural sclerotia, but validation experiments were performed with lab-grown sclerotia, the assay’s ability to detect Rtul DNA in both naturally found and lab-grown sclerotia was compared. The sensitivity of Rtul detection from light brown sclerotia was two-fold and seven-fold higher then corresponding natural black sclerotia and lab-grown sclerotia, respectively; thus, confirming the suitability of the qPCR assay to detect Rtul in the field.
The presence of soil did not affect the efficiency of DNA extraction or detection of sclerotia using qPCR, as observed by Guo et al. [33] for the detection and quantification of R. cerealis in soil. In fact, sclerotial DNA extraction in the presence of soil seemed to remove or neutralize PCR inhibitors. Fungi sequester many secondary metabolites in sclerotia that act as a chemical defense system against competitors, fungivorous predators and parasites [34], although the effect of such chemicals on PCR has not been investigated [35].
The ability to detect Rtul DNA decreased with increased tulip bulb tissue, such that the equivalent of almost ten sclerotia would need to be present if one whole tulip bulb was extracted. This limited sensitivity has the potential to limit the effectiveness of the Rtul assay in quantifying inoculum levels in planting stock. It is likely the decreased sensitivity of the Rtul assay associated with increased bulb tissue may be due to the presence of qPCR inhibitors in the dormant tulip bulb tissues [36]. A variety of protocols have been developed to remove inhibitors from nucleic acids before PCR [37]. Additional research is needed to determine if removal of potential inhibitors from tulip bulb tissues to increase the sensitivity of the whole bulb assay, and thereby maximize the likelihood of detecting bulb infections, would be beneficial. Moreover, Tellenbach et al. [38] reported increased sensitivity of qPCR detection from sclerotial DNA when probes were made from a single-copy gene as opposed to a multicopy gene like ITS. This enhancement effect was attributed to dilution of PCR inhibitors after the first amplification step. With the genome information of Rtul potentially available, future enhancement of pathogen detection from sclerotial DNA may be possible with probes designed from single copy genes [39].
We also conducted ITS (ITS1-5.8S-ITS2)-based phylogenetic characterization and ITS2-based molecular morphometric characterization of three geographically distinct Rtul isolates and demonstrated divergent relationships among them as well as their taxonomic position with fungi belonging to Ceratobasidiaceae. The Bayesian inference of phylogeny revealed that the three isolates of Rtul investigated have originated as a distinct clade, but closer to another distinct clade consisting of a binucleate Rhizoctonia, C. anceps. Besides, the Rtul isolates were more distant from other Ceratobasidium spp., including Ceratobasidium AGs A to U. This hypothesis is well-supported from the ITS2 molecular morphometric analyses. ITS2 of C. anceps has lower GC content and is one nucleotide shorter than those of the Rtul isolates (Table S1). Also, secondary structures of the two species have at least one CBC and one hCBC in helix II, indel of three base pairs in helix III, and indel of one base pair in helix IV (Figure 7).
Recent studies have demonstrated that ITS2 sequence evolves very quickly amongst eukaryotes, making it an excellent marker that is widely used in low level phylogenetic analyses and DNA barcoding of organisms, including fungi [40,41,42,43]. By the same token, any hypervariable feature may also be problematic for investigating deeper phylogeny. This could be overcome by its conserved secondary structure. The conservation of ITS2 secondary structure is a critical component of rRNA processing [41]. Thus, any compensatory base change (CBC) and hemi-compensatory base change (hCBC) in the two most conserved helices, namely helix II and helix III, have been correlated with delineation of species [44]. Speciation in eukaryotes, particularly amongst plants and fungi, has been correlated with even one CBC in the relatively conserved pairing positions of the ITS2 transcript secondary structure [45]. The probability to discriminate between two species with a single CBC between the members of an ITS2 sequence-structure pair has a confidence of 93% [31,45]. On the other hand, if there is no CBC, then individuals may belong to the same species with 76% confidence [45]; Wolf et al. [31]. Molecular morphometric analyses utilize the ITS2 primary sequence as well as secondary structural constraints, homologous locations, indels and substitution to distinguish closely related fungal and other species [21,22]. As a proof of concept for Rhizoctonia species, Ahvenniemi et al. [46] showed that the potato (PT) and tobacco (TB) subgroups of R. solani AG 3 differ in the ITS2 secondary structure helix I by two CBCs. Moreover, the CBCs in potato strains of AG 3 distinguish them from other AGs of R. solani, confirming the hypothesis that the latter is functionally a composite species. Interestingly, our investigations on the ITS2 molecular morphometrics of the three geographically distant Rtul isolates demonstrated that major genetic shifts have taken place resulting in significant alteration of respective ITS2 secondary structures. Thus, although the ITS (ITS-5.8S-ITS2) primary sequences of all the three isolates of Rtul formed a well-supported clade in the phylogenetic tree (Figure 6), molecular morphometric analyses clearly differentiated them as three potentially distinct sexual species within R. tuliparum, sensu lato, per Müller et al. [45] and Wolf et al. [31] (Figure 8; Table S4), and call for further investigation. Moreover, from disease epidemiology and management viewpoints, our findings call for accurate identification of pathogen genotypes as well as determination of sensitivities of genetically diverse isolates of Rhizoctonia pathogens to commonly used fungicides [13,43,47,48,49,50].
The significance of the formation of single base pair pseudoknot-like structure in the stem III of ITS2 in Rtul and phylogenetically related species (Figure 6; Table S5) is not understood at present but may serve as an additional criterion to differentiate groups or strains of Ceratobasidium spp. Sharon et al. [51] phylogenetically clustered the Ceratobasidium isolates AG-A to AG-U in six groups. In this regard, it is interesting to note that the structure was present in group 4 (AG-D), group 5 [AG-B (including its three subgroups) and AG-Q], and group 6 (AG-C, AG-H, and AG-I) and absent in group 3 (AG-E, AG-F, AG-P AG-U, AG-R, and AG-S). Mixed results were also noted because in group 1 (e.g., AG-A has the pseudoknot, but not in AG-K), and in group 2 (e.g., AG-G and AG-L has the pseudoknot but absent in AG-O).
Taken together, both the development of a sensitive qPCR detection method and determination of the genetic diversity of the Rtul pathogens as documented in this investigation are significant for disease management. While cultural practices such as deep plowing to bury sclerotia, shallow planting, and flooding help control this disease, treatments such as soil fumigation, fungicide bulb dips or in-furrow or broadcast applications of fungicides to the soil at the time of planting are also often used as methods of disease control [3,4,9,12,52,53]. In this context, additional research is needed to determine soil-borne inoculum thresholds through qPCR, as developed in this investigation, for disease development. Given the spatial distribution of the disease, additional research is also needed to define a sampling protocol that could be used along with this qPCR diagnostic assay to accurately predict the risk of disease development. Furthermore, the presence of genetic diversity of Rtul pathogens as determined in this report also calls for assaying differential sensitivity of Rtul isolates to various fungicides and finding the most efficient fungicide for each Rtul genotype. A sensitive detection method, along with the use of the most effective fungicide, would ultimately give rise to environmentally and economically sound management of gray bulb rot in tulips and bulbus irises [54,55,56].

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/jof8020163/s1. Table S1: The ITS2 sequence length and base composition of C. anceps and three R. tuliparum (Rtul) isolates, Table S2: Base-pairing in the ITS2 RNA secondary structures of Rtul Washington isolate (KX767078) in comparison to that of C. anceps (MH855251), Table S3: The conserved terminal loop sequences in all isolates of C. anceps (MH855251 and AJ427402), and Rtul (KX767078, EU191041, and MH854847), Table S4: Base-pairing in the ITS2 RNA secondary structures and the conserved terminal loop sequences of the European (EU191041) and the New York (MH854847) isolates of Rtul with respect to the Washington isolate (KX767078), Table S5: Presence or absence of pseudoknot-like structure in helix III of Ceratobasidium binucleate and uninucleate isolates within the greater clade containing Rtul in Figure 6. Table S6: GenBank Accessions of the ITS (ITS1-5.8S-ITS2) regions showing complete sequence match with the Rtul Q-PCR Probe. Figure S1: Selected ITS regions from alignment file of Rhizoctonia pathogens for designing primers and probes for the Rtul Q-PCR detection, Figure S2: Morphologies of “lab-grown” sclerotia and a mixture of “natural light” and “natural dark” sclerotia.

Author Contributions

K.C. conceptualization, data generation, analysis, interpretation and writing initial draft; A.D. data generation and analysis; D.K.L. data generation, analysis, interpretation, and writing; A.I. data generation, analysis, interpretation, and writing; D.P.R. data curation and editing; G.C. conceptualization, resources and funding acquisition, writing, editing. All authors have read and agreed to the published version of the manuscript.

Funding

A portion of the research was funded by Cooperative Agreements between USDA-ARS and Washington State University (58-8042-6-091 and 58-8042-1-059) and funding from the Northwest Agricultural Research Foundation and Wally Staatz Turfgrass/Ornamental Endowment. Funding was also provided by USDA-ARS in-house project 8042 21220 181 000 D.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All relevant data generated or analyzed during this study are included in this article.

Acknowledgments

In-kind support from Washington Bulb Company and Washington State Bulb Grower Association is appreciated.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Whetzel, H.H.; Arthur, J.M. Gray Bulb-Rot of Tulips Caused By Rhizoctonia tuliparum (Klebh.) N. Comb; Cornell University: Ithaca, NY, USA, 1925. [Google Scholar]
  2. Buddin, W. The grey bulb rot of tulips and its control. J. Minist. Agric. Fish. 1937, 44, 1154–1159. [Google Scholar]
  3. Moore, W.C. Grey Bulb Rot of tulip. J. Hortic. Soc. 1950, 75, 113–117. [Google Scholar]
  4. Moore, W.C. Diseases of Bulbs, 2nd ed.; Dickens, J.S.W., Ed.; Her Majesty’s Stationary Office: London, UK, 1979; 205p. [Google Scholar]
  5. Boerema, G.H.; Hamers, M.E.C. Check-list for scientific names of common parasitic fungi. Series 3a: Fungi on bulbs: Liliaceae. Nether. J. Plant Pathol. 1988, 94, 1–29. [Google Scholar] [CrossRef]
  6. Gladders, P.; Coley-Smith, J.R. Rhizoctonia tuliparum: A winter-active pathogen. Trans. Brit. Mycol. Soc. 1978, 71, 129–139. [Google Scholar] [CrossRef]
  7. Wakker, J.H. Onderzoek der Ziekten van Hyacinthen en Andere Bol- En Knolgewassen; Algemeene Vereeniging voor Bloembollencultuur: Haarlem, The Netherlands, 1884; 26p. [Google Scholar]
  8. Whetzel, H.H.; Arthur, J.M. The grey bulb-rot of tulips. Phytopathology 1924, 14, 30–31. [Google Scholar]
  9. McGovern, R.J.; Elmer, W.H. Diseases of Tulip. In Handbook of Florists’ Crops Diseases. Handbook of Plant Disease Management; McGovern, R., Elmer, W., Eds.; Springer: Cham, Switzerland, 2017. [Google Scholar] [CrossRef]
  10. Coley-Smith, J.R.; Humphreys-Jones, D.R.; Gladders, P. Long-term survival of sclerotia of Rhizoctonia tuliparum. Plant Pathol. 1979, 28, 128–130. [Google Scholar] [CrossRef]
  11. Blok, W.; Coenen, T.C.M.; Termorshiuzen, A.J.; Lamers, J.G. The potential of biological soil disinfestation to manage fusarium foot and root rot in asparagus. Acta Hortic. 2008, 776, 135–144. [Google Scholar] [CrossRef] [Green Version]
  12. Chastagner, G.A.; DeBauw, A. Alternatives to PCNB for controlling Gray Bulb Rot on tulips. Acta Hort. (ISHS) 2011, 886, 311–317. [Google Scholar] [CrossRef]
  13. Lakshman, D.K.; Cloyd, R.A.; Chastagner, G.A. Integrated management of diseases and pests on ornamental geophytes: Challenges and progress. Acta Hortic. 2019, 1237, 13–32. [Google Scholar] [CrossRef]
  14. Leclerc, M.; Doré, T.; Gilligan, C.A.; Lucas, P.; Filipe, J.A. Host growth can cause invasive spread of crops by soilborne pathogens. PLoS ONE 2013, 8, e63003. [Google Scholar] [CrossRef] [Green Version]
  15. Termorshuizen, A.J.; Jeger, M.J. Assessing inoculum of soilborne plant pathogens: Theory and practice. Acta Hortic. 2014, 1044, 75–80. [Google Scholar] [CrossRef]
  16. Wei, F.; Fan, R.; Dong, H.T.; Shang, W.-J.; Xu, X.M.; Zhu, H.Q.; Yang, J.R.; Hu, X.P. Threshold microsclerotial inoculum for cotton Verticillium wilt determined through wet-sieving and real-time quantitative PCR. Phytopathology 2015, 105, 220–229. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Panth, M.; Hassler, S.C.; Baysal-Gurel, F. Methods for management of soilborne diseases in crop production. Agriculture 2020, 10, 16. [Google Scholar] [CrossRef] [Green Version]
  18. Heid, C.A.; Stevens, J.; Livak, K.J.; Williams, P.M. Real time quantitative PCR. Genome Res. 1996, 10, 986–994. [Google Scholar] [CrossRef] [Green Version]
  19. Burton, R.J.; Coley-Smith, J.R. Antibiotics in sclerotia and mycelium of Rhizoctonia species. Trans. Br. Mycol. Soc. 1985, 85, 447–453. [Google Scholar] [CrossRef]
  20. Dijst, G.; Schneider, J.H.M. Flower bulb diseases incited by Rhizoctonia species. In Rhizoctonia Species: Taxonomy, Molecular Biology, Ecology, Pathology and Disease Control; Sneh, B., Jabaji-Hare, S., Neate, S.M., Dijst, G., Eds.; Kluwer: Dordrecht, The Netherlands, 1996; pp. 279–288. [Google Scholar]
  21. Sundaresan, N.; Jagan, E.G.; Kathamuthu, G.; Pandi, M. Internal transcribed spacer 2 (ITS2) molecular morphometric analysis-based species delimitation of foliar endophytic fungi from Aglaia elaeagnoidea, Flacourtia inermis and Premna serratifolia. PLoS ONE 2019, 14, e0215024. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Billoud, B.; Guerrucci, M.A.; Masselot, M.; Deutsch, J.S. Cirrepede phylogeny using a novel approach: Molecular morphometrics. Mol. Biol. Evol. 2020, 17, 1435–1445. [Google Scholar] [CrossRef] [Green Version]
  23. Haugland, R.A.; Siefring, S.C.; Wymer, L.J.; Brenner, K.P.; Dufour, A.P. Comparison of Enterococcus measurements in freshwater at two recreational beaches by quantitative Polymerase Chain Reaction and membrane filter culture analysis. Water Res. 2005, 39, 559–568. [Google Scholar] [CrossRef]
  24. Ruijter, J.M.; Ramakers, C.; Hoogaars, W.M.; Karlen, Y.; Bakker, O.; van den Hoff, M.J.; Moorman, A.F. Amplification efficiency: Linking baseline and bias in the analysis of quantitative PCR data. Nucleic Acids Res. 2009, 37, e45. [Google Scholar] [CrossRef] [Green Version]
  25. Thompson, J.D.; Gibson, T.J.; Plewniak, F.; Jeanmougin, F.; Higgins, D.G. The Clustal X Windows interface; flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucl. Acids Res. 1997, 24, 4876–4882. [Google Scholar] [CrossRef] [Green Version]
  26. Maddison, W.P.; Maddison, D.R. Mesquite: A Modular System for Evolutionary Analysis. 2019, Version 3.61. Available online: http://www.mesquiteproject.org (accessed on 10 November 2021).
  27. Huelsenbeck, J.P.; Ronquist, F. MrBayes: Bayesian inference of phylogeny trees. Bioinformatics 2001, 17, 754–755. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Swofford, D.L. PAUP*: Phylogenetic Analysis Using Parsimony (*and Other Methods); Version 4.0b10; Sinauer Associates: Sunderland, MA, USA, 2002. [Google Scholar] [CrossRef]
  29. Vu, D.; Groenewald, M.; de Vries, M.; Gehrmann, T.; Stielow, B.; Eberhardt, U.; Al-Hatmi, A.; Groenewald, J.Z.; Cardinal, G.; Houbraken, J.; et al. Large-scale generation and analysis of filamentous fungal DNA barcodes boosts coverage for kingdom fungi and reveals thresholds for fungal species and higher taxon delimitation. Stud. Mycol. 2019, 92, 135–154. [Google Scholar] [CrossRef]
  30. Ankenbrand, M.J.; Keller, A.; Wolf, M.; Schultz, J.; Förster, F. ITS2 database V: Twice as much. Mol. Biol. Evol. 2015, 32, 3030–3032. [Google Scholar] [CrossRef] [PubMed]
  31. Wolf, M.; Chen, S.; Song, J.; Ankenbrand, M.; Müller, T. Compensatory base changes in ITS2 secondary structures correlate with the biological species concept despite intragenomic variability in ITS2 sequences—A proof of concept. PLoS ONE 2013, 8, e66726. [Google Scholar] [CrossRef] [PubMed]
  32. Ceresini, P.C.; Costa-Souza, E.; Zala, M.; Furtado, E.L.; Souza, N.L. Evidence that the Ceratobasidium-like white-thread blight and black rot fungal pathogens from persimmon and tea crops in the Brazilian Atlantic Forest agroecosystem are two distinct phylospecies. Genet. Mol. Biol. 2012, 35, 480–497. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Guo, Y.; Sun, L.W.H.; Wang, N.; Yu, H.; Chen, H. Detection and quantification of Rhizoctonia cerealis in soil using real-time PCR. J. Gen. Plant Pathol. 2012, 78, 247–254. [Google Scholar] [CrossRef]
  34. Calvo, A.M.; Cary, J.W. Association of fungal secondary metabolism and sclerotial biology. Front. Microbiol. 2015, 6, 62. [Google Scholar] [CrossRef] [Green Version]
  35. Shu, C.; Sun, S.; Chen, J.; Chen, J.; Zhou, E. Comparison of different methods for total RNA extraction from sclerotia of Rhizoctonia solani. Elec. J. Biotechnol. 2014, 17, 50–54. [Google Scholar] [CrossRef] [Green Version]
  36. Wei, T.; Lebas, B.S.M.; Shiller, J.B.; Quinn, B.D.; Clover, G.R.G. Detection of five viruses infecting dormant bulbs by TaqMan-based real-time RT-PCR. Australas. Plant Pathol. 2012, 41, 93–98. [Google Scholar] [CrossRef]
  37. Schrader, C.; Schielke, A.; Ellerbroek, L.; Johne, R. PCR inhibitors—Occurrence, properties, and removal. J. Appl. Microbiol. 2012, 113, 1014–1026. [Google Scholar] [CrossRef]
  38. Tellenbach, C.; Grünig, C.R.; Sieber, T.N. Suitability of quantitative real-time PCR to estimate biomass of fungal root endophytes. Appl. Environ. Microbiol. 2010, 76, 5764–5772. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Feau, N.; Decourcelle, T.; Husson, C.; Desprez-Loustau, M.-L.; Dutech, C. Finding single copy genes out of sequenced genomes for multilocus phylogenetics in non-model fungi. PLoS ONE 2011, 6, e18803. [Google Scholar] [CrossRef] [PubMed]
  40. Yang, R.H.; Su, J.H.; Shang, J.J.; Wu, Y.Y.; Li, Y.; Bao, D.P.; Yao, Y.J. Evaluation of the ribosomal DNA internal transcribed spacer (ITS), specifically ITS1 and ITS2, for the analysis of fungal diversity by deep sequencing. PLoS ONE 2018, 13, e0206428. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Mbareche, H.; Veillette, M.; Bilodeau, G.; Duchaine, C. Comparison of the performance of ITS1 and ITS2 as barcodes in amplicon-based sequencing of bioaerosols. PeerJ 2020, 8, e8523. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Zhang, W.; Tian, W.; Gao, Z.; Wang, G.; Zhao, H. Phylogenetic utility of rRNA ITS2 sequence-structure under functional constraint. Int. J. Mol. Sci. 2020, 21, 6395–6407. [Google Scholar] [CrossRef] [PubMed]
  43. Singh, V.; Lakshman, D.K.; Roberts, D.P.; Ismaiel, A.; Hooda, K.S.; Gogoi, R. Morphopathological and molecular morphometric characterization of Waitea circinata var. prodigus causing a novel sheath spot disease of maize in India. Plant Dis. 2021. [Google Scholar] [CrossRef]
  44. Coleman, A.W. Is there a molecular key to the level of “biological species” in eukaryotes? A DNA guide. Mol. Phylogenet. Evol. 2009, 50, 197–203. [Google Scholar] [CrossRef]
  45. Müller, T.; Philippi, N.; Dandekar, T.; Schultz, J.; Wolf, M. Distinguishing species. RNA 2007, 13, 1469–1472. [Google Scholar] [CrossRef] [Green Version]
  46. Ahvenniemi, P.; Wolf, M.; Lehtonen, M.J.; Wilson, P.; German-Kinnari, M.; Valkonen, J.P.T. Evolutionary diversification indicated by compensatory base changes in ITS2 secondary structures in a complex fungal species, Rhizoctonia solani. J. Mol. Evol. 2009, 69, 150–163. [Google Scholar] [CrossRef]
  47. Kammerer, S.J.; Burpee, L.L.; Harmon, P.F. Identification of a new Waitea circinata variety causing basal leaf blight of seashore paspalum. Plant Dis. 2011, 95, 515–522. [Google Scholar] [CrossRef] [Green Version]
  48. Amaradasa, B.S.; Lakshman, D.K.; Horvath, B.J.; Amundsen, K.L. Development of SCAR markers and UP-PCR cross-hybridization method for specific detection of four major subgroups of Rhizoctonia from infected turfgrasses. Mycologia 2014, 106, 163–172. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Amaradasa, B.S.; Lakshman, D.K.; McCall, D.; Horvath, B.J. In vitro fungicide sensitivity of Rhizoctonia and Waitea isolates collected from turfgrasses. J. Environ. Hort. 2014, 32, 126–132. [Google Scholar] [CrossRef]
  50. Amaradasa, B.S.; Lakshman, D.K.; Amundsen, K. AFLP fingerprinting for identification of infra-species groups of Rhizoctonia solani and Waitea circinata. J. Plant Protect. Microbiol. 2015, 6, 3. [Google Scholar] [CrossRef] [Green Version]
  51. Sharon, M.S.; Kuninaga, S.; Hyakumachi, M.; Naito, S.; Sneh, B. Classification of Rhizoctonia spp., using rDNA-ITS sequence analysis supports the genetic basis of the classical anastomosis grouping. Mycoscience 2008, 49, 93–114. [Google Scholar] [CrossRef]
  52. van Klaveren, C. Ziegnten in bloembollen. Maandbl. Landbouw Voorl. 1949, 6, 147–153. [Google Scholar]
  53. Muller, P.J.; Vink, P.; Van Zaagen, A. Flooding caused loss in viability and pathogenicity of sclerotia of Rhizoctonia tuliparum. Eur. J. Plant Pathol. 1988, 94, 45–47. [Google Scholar] [CrossRef]
  54. Chastagner, G.; Hanks, G.; Daughtrey, M.; Yedidia, I.; Miller, T.; Pappu, H. Sustainable Production and Integrated Management: Environmental Issues. In Ornamental Geophytes: From Basic Science to Sustainable Horticultural Production; Kamenetsky, R., Okubo, H., Eds.; CRC Press: Boca Raton, FL, USA, 2012; pp. 363–420. [Google Scholar] [CrossRef]
  55. Conijn, C. Tulip Diseases. Roodbont Publishers B.V.: Zutphen, The Netherlands, 2016; p. 149. [Google Scholar]
  56. Chastagner, G.A.; Garfinkel, A. Diseases Affecting Ornamental Geophytes and Their Control. In Achieving Sustainable Cultivation of Ornamental Plants; Reid, M., Ed.; Burleigh Dodds Science Publishing: Cambridge, UK, 2020; pp. 367–414. [Google Scholar] [CrossRef]
Jof 08 00163 g001 550
Figure 1. Determining the lowest threshold of detection with the Rtul primer/probe set used in this study. 1:10 serial dilution series (1.2 ng μL−1–0.12 fg μL−1 of DNA) of Rhizoctonia tuliparum (DGF2-VV) DNA analyzed in triplicate with qPCR. The log DNA concentration versus cycle threshold was plotted and a regression line drawn for the data.
Figure 1. Determining the lowest threshold of detection with the Rtul primer/probe set used in this study. 1:10 serial dilution series (1.2 ng μL−1–0.12 fg μL−1 of DNA) of Rhizoctonia tuliparum (DGF2-VV) DNA analyzed in triplicate with qPCR. The log DNA concentration versus cycle threshold was plotted and a regression line drawn for the data.
Jof 08 00163 g001
Jof 08 00163 g002 550
Figure 2. Effect of dilution on the sensitivity of detection of DNA from Rhizoctonia tuliparum sclerotia. DNAs extracted from lab-grown sclerotia in triplicate were diluted and analyzed by the Rtul qPCR assay, and the detected DNA for each dilution plotted. Error bars represent one standard deviation from the mean. The assays were more precise when the DNA was diluted in the range of 1:10 to 1:1000 before qPCR.
Figure 2. Effect of dilution on the sensitivity of detection of DNA from Rhizoctonia tuliparum sclerotia. DNAs extracted from lab-grown sclerotia in triplicate were diluted and analyzed by the Rtul qPCR assay, and the detected DNA for each dilution plotted. Error bars represent one standard deviation from the mean. The assays were more precise when the DNA was diluted in the range of 1:10 to 1:1000 before qPCR.
Jof 08 00163 g002
Jof 08 00163 g003 550
Figure 3. Sensitivity of detection of Rhizoctonia tuliparum DNA from lab-grown, natural-light, and natural-dark sclerotia. DNAs were extracted from 14.25 mg sclerotia (equivalent of eight average size sclerotia) of three different types in triplicate diluted 1/1000 and analyzed by qPCR. The detected DNAs for each type were plotted in a bar chart with error bars representing one standard deviation from the mean.
Figure 3. Sensitivity of detection of Rhizoctonia tuliparum DNA from lab-grown, natural-light, and natural-dark sclerotia. DNAs were extracted from 14.25 mg sclerotia (equivalent of eight average size sclerotia) of three different types in triplicate diluted 1/1000 and analyzed by qPCR. The detected DNAs for each type were plotted in a bar chart with error bars representing one standard deviation from the mean.
Jof 08 00163 g003
Jof 08 00163 g004 550
Figure 4. The Rtul qPCR assay was used to evaluate the effect of soil on the detection of Rhizoctonia tuliparum sclerotial DNA. The average amount of R. tuliparum DNA detected when the weight-equivalent of one to 32 sclerotia was added to four grams of soil was fairly consistent, with the exception of the highest amount of sclerotia that yielded more DNA in the presence of soil.
Figure 4. The Rtul qPCR assay was used to evaluate the effect of soil on the detection of Rhizoctonia tuliparum sclerotial DNA. The average amount of R. tuliparum DNA detected when the weight-equivalent of one to 32 sclerotia was added to four grams of soil was fairly consistent, with the exception of the highest amount of sclerotia that yielded more DNA in the presence of soil.
Jof 08 00163 g004
Jof 08 00163 g005 550
Figure 5. Analysis of the effect of tulip bulb tissue on sensitivity of qPCR detection of Rhizoctonia tuliparum DNA. 14.25 mg sclerotia (equivalent of 8 average size lab-grown sclerotia) were added to various fractions of tulip bulb tissue to determine the limit of detection of sclerotial DNA in infected bulbs by Rtul qPCR. The amount of R. tuliparum DNA detected was compared to the amount of DNA in sclerotia in the absence of bulb tissue. The presence of bulb tissue reduced the ability to detect R. tuliparum.
Figure 5. Analysis of the effect of tulip bulb tissue on sensitivity of qPCR detection of Rhizoctonia tuliparum DNA. 14.25 mg sclerotia (equivalent of 8 average size lab-grown sclerotia) were added to various fractions of tulip bulb tissue to determine the limit of detection of sclerotial DNA in infected bulbs by Rtul qPCR. The amount of R. tuliparum DNA detected was compared to the amount of DNA in sclerotia in the absence of bulb tissue. The presence of bulb tissue reduced the ability to detect R. tuliparum.
Jof 08 00163 g005
Jof 08 00163 g006 550
Figure 6. Bayesian tree obtained based on internal transcribed spacer sequences illustrating phylogeny of Rhizoctonia tuliparum within closely related species. One million generations with 25% of the initial trees discarded in the burn-in phase were used for analysis. Posterior probabilities >0.95 are shown above the branches. The tree was rooted with two isolates of Waitea circinata. The symbol (*) under branches indicates bootstrap support ≥ 70% obtained from parsimony analysis.
Figure 6. Bayesian tree obtained based on internal transcribed spacer sequences illustrating phylogeny of Rhizoctonia tuliparum within closely related species. One million generations with 25% of the initial trees discarded in the burn-in phase were used for analysis. Posterior probabilities >0.95 are shown above the branches. The tree was rooted with two isolates of Waitea circinata. The symbol (*) under branches indicates bootstrap support ≥ 70% obtained from parsimony analysis.
Jof 08 00163 g006
Jof 08 00163 g007 550
Figure 7. The ITS2 secondary structure of Ceratobasidium anceps, isolate CBS 152.32 (GenBank Accession MH855251). as predicted by the ITS2 DB ‘Predict’ utility. The nucleotides are highlighted in green (C), orange (G), violet (T), and pink (A). The structural changes in ITS2 of Rhizoctonia tuliparum isolate I-399 (GenBank Accession KX767078) are indicated by compensatory base change (CBC) in the blue box, the hemi-compensatory base change (hCBC) in the orange box, and the loss of base pairs in light green boxes. The conserved terminal loop sequences are marked with curved lines and indicated in black oval boxes in all four helices.
Figure 7. The ITS2 secondary structure of Ceratobasidium anceps, isolate CBS 152.32 (GenBank Accession MH855251). as predicted by the ITS2 DB ‘Predict’ utility. The nucleotides are highlighted in green (C), orange (G), violet (T), and pink (A). The structural changes in ITS2 of Rhizoctonia tuliparum isolate I-399 (GenBank Accession KX767078) are indicated by compensatory base change (CBC) in the blue box, the hemi-compensatory base change (hCBC) in the orange box, and the loss of base pairs in light green boxes. The conserved terminal loop sequences are marked with curved lines and indicated in black oval boxes in all four helices.
Jof 08 00163 g007
Jof 08 00163 g008 550
Figure 8. The ITS2 secondary structure of Rhizoctonia tuliparum isolate I-399 from Washington State (GenBank Accession KX767078), as predicted by the ITS2 DB ‘Predict’ utility. The nucleotides are highlighted in green (C), orange (G), violet (T), and pink (A). The structural changes in ITS2 of the R. tuliparum isolate 29,792 from Europe [EU19104], and the isolate CBS206.25 from New York [MH854847] are indicated as gain of a base pair due to a hemi-compensatory base change (hCBC) in the blue box and loss of base pairs due to hCBCs in the light green boxes. A hCBC resulting in loss of a pseudoknot in the isolate CBS206.25 [MH854847] is indicated in the dark green box.
Figure 8. The ITS2 secondary structure of Rhizoctonia tuliparum isolate I-399 from Washington State (GenBank Accession KX767078), as predicted by the ITS2 DB ‘Predict’ utility. The nucleotides are highlighted in green (C), orange (G), violet (T), and pink (A). The structural changes in ITS2 of the R. tuliparum isolate 29,792 from Europe [EU19104], and the isolate CBS206.25 from New York [MH854847] are indicated as gain of a base pair due to a hemi-compensatory base change (hCBC) in the blue box and loss of base pairs due to hCBCs in the light green boxes. A hCBC resulting in loss of a pseudoknot in the isolate CBS206.25 [MH854847] is indicated in the dark green box.
Jof 08 00163 g008
Table
Table 1. Fungal isolates used to test species specificity of the TaqMan qPCR assay for Rhizoctonia tuliparum, and their resulting Rtul cycle threshold (Ct) value.
Table 1. Fungal isolates used to test species specificity of the TaqMan qPCR assay for Rhizoctonia tuliparum, and their resulting Rtul cycle threshold (Ct) value.
IsolateSpeciesHost/SourceOriginCt Value
PC76Rhizoctonia solani AG-2-1TulipWashingtonN/A a
Rh060811 cR. solani AG 2-1PeaWashingtonN/A
Rh070913 cR. solani AG 2-1OnionWashingtonN/A
Rh070937 cR. solani AG 2-1PeaWashingtonN/A
Rh070933 cR. solani AG 3PotatoWashingtonN/A
Rh070934 cR. solani AG 3PotatoWashingtonN/A
Rh070935 cR. solani AG 3OnionWashingtonN/A
Rh070942 cR. solani AG 3PeaWashingtonN/A
Rh070912 cR. solani AG 3OnionWashingtonN/A
Rh010901 cR. solani AG 4PeaWashingtonN/A
Rh070908 cR. solani AG 4PeaWashingtonN/A
Rh070909 cR. solani AG 4PeaWashingtonN/A
Rh070910 cR. solani AG 4PeaWashingtonN/A
Rh070915 cR. solani AG 4OnionWashingtonN/A
Rh070929 cR. solani AG 4OnionWashingtonN/A
Rh070939 cR. solani AG 4OnionWashingtonN/A
Rh070940 cR. solani AG 4PotatoWashingtonN/A
Rh070930 cR. solani AG 5PotatoWashingtonN/A
Rh070931 cR. solani AG 5PeaWashingtonN/A
Rh070932 cR. solani AG 5PotatoWashingtonN/A
Rh080714 cR. solani AG 8WheatWashingtonN/A
Rh010911 cR. solani AG 8OnionWashingtonN/A
Rh070922 cR. solani AG 8Soil dWashingtonN/A
Rh070927 cR. solani AG 8OnionWashingtonN/A
Rh070943 cR. solani AG 8PeaWashingtonN/A
Rh100557 cR. solani AG 8WheatWashingtonN/A
Rh070921 cR. solani AG 9PeaWashingtonN/A
Rh070938 cR. solani AG 9PeaWashingtonN/A
Rh090801 cCeratobasidium AG-AOnionWashingtonN/A
Rh010913 cCeratobasidium AG-AOnionWashingtonN/A
Rh070923 cCeratobasidium AG-ESoil dWashingtonN/A
Rh110712 cCeratobasidium AG-IWheatWashingtonN/A
Rh070716 cCeratobasidium AG-IWheatWashingtonN/A
Rh070914 cCeratobasidium AG-IPeaWashingtonN/A
Rh010905 cCeratobasidium sp.OnionWashingtonN/A
Rh010909 cWaitea circinata (R. oryzae grp. 1)OnionWashingtonN/A
Rh070924 cW. circinata (var. circinata)Soil dWashingtonN/A
Rh070925 cW. circinata (var. circinata)Soil dWashingtonN/A
Rh070936 cW. circinata (var. circinata)OnionWashingtonN/A
Rh070911 cCeratobasidium sp.PeaWashingtonN/A
Rh070926 cRhizoctonia spp.Soil dWashingtonN/A
Rh070928 cRhizoctonia spp.OnionWashingtonN/A
VSP 05-01 A cR. solani AG 4, HG IIOnionColumbia BasinN/A
VSP 05-01 B cR. solani AG 4, HG IIOnionColumbia BasinN/A
VSP 05-37 A c Swiss ChardWashingtonN/A
VSP 05-37 B c Swiss ChardWashingtonN/A
VSP 06-26 B cR. solani AG 4, HG IIOnionColumbia BasinN/A
VSP 08-13 A c OnionColumbia BasinN/A
VSP 08-14 B c OnionColumbia BasinN/A
VSP 08-19 A c OnionColumbia BasinN/A
VSP 08-19 G c OnionColumbia BasinN/A
VSP 08-33 B, plt 4 c SkullcapColumbia BasinN/A
06SPFWLA cR. solani AG 3SpinachWashingtonN/A
VSP 08-33 A, plt 6 cR. solaniSkullcapColumbia BasinN/A
VSP 10-20 Tamara 1 cR. solaniOnionColumbia BasinN/A
VSP 10-20 Tamara 2 cR. solaniOnionColumbia BasinN/A
VSP 10-20 Field 2-2 cR. solaniOnionColumbia BasinN/A
VSP 10-20 Field 2-3 cR. solaniOnionColumbia BasinN/A
LP 3 cR. solaniPeaOregonN/A
LP 10 cR. solaniPeaOregonN/A
LP 12 cR. solaniPeaOregonN/A
LP 13 cR. solaniPeaOregonN/A
MBL12140R. tuliparumIrisWashington19.60
MBL12141R. tuliparumIrisWashington16.10
MBL12145R. tuliparumIrisWashington27.25
MBL12146R. tuliparumIrisWashington17.84
MBL11082 (DGF2_VV)R. tuliparumIrisWashington33.51
Rh_TulipR. tuliparumTulipWashington13.93
I-399 bR. tuliparumIrisWashington
I-399R. tuliparum 1:10IrisWashington18.79
I-399R. tuliparum 1:100IrisWashington21.32
I-399R. tuliparum 1:1000IrisWashington24.90
No template control e---------N/A
a N/A = no amplification; b I-399 DNA used for standard curve; c Cultures provided by Lindsey du Toit, Washington State University, Mount Vernon Northwestern Washington Research and Extension Center; d Isolated from wheat cultivated soil; e No DNA template control.
Table
Table 2. Primers and probes used for the TaqMan® quantitative PCR diagnostic assay for detection and quantification of Rhizoctonia tuliparum and the associated exogenous amplification control a.
Table 2. Primers and probes used for the TaqMan® quantitative PCR diagnostic assay for detection and quantification of Rhizoctonia tuliparum and the associated exogenous amplification control a.
Primer/Probe NameSequence (5′-3′)Target RegionStudy
R. tuliparum primers
Rtul_forCGAGGTCGACTTTTTGTTTTGGITSaThis study
Rtul_revCCGAGTGGAACCGAGTTCACITS This study
R. tuliparum probe
R. tuliparum probe
Rtul_probe[TET] TTTGCGGATTCACGTCC [MGB-NFQ]ITS This study
Sketa Exogenous Amplification Control (EAC) primers
SketaF2GGTTTCCGCAGCTGGGITS [23]
SketaR3CCGAGCCGTCCTGGTCTAITS [23]
Sketa probe
SketaP2[6-FAM] AGTCGCAGGCGGCCACCGT [TAMRA]ITS[23]
a ITS, Internal Transcribed Spacer region of ribosomal DNA; TET, tetrachlorofluorescein; MGB-NFQ, minor groove binder nonfluorescent quencher; 6-FAM, fluorescein (6-FAM) dye; TAMRA, tetramethylrhodamine quencher.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Coats, K.; DeBauw, A.; Lakshman, D.K.; Roberts, D.P.; Ismaiel, A.; Chastagner, G. Detection and Molecular Phylogenetic-Morphometric Characterization of Rhizoctonia tuliparum, Causal Agent of Gray Bulb Rot of Tulips and Bulbous Iris. J. Fungi 2022, 8, 163. https://doi.org/10.3390/jof8020163

AMA Style

Coats K, DeBauw A, Lakshman DK, Roberts DP, Ismaiel A, Chastagner G. Detection and Molecular Phylogenetic-Morphometric Characterization of Rhizoctonia tuliparum, Causal Agent of Gray Bulb Rot of Tulips and Bulbous Iris. Journal of Fungi. 2022; 8(2):163. https://doi.org/10.3390/jof8020163

Chicago/Turabian Style

Coats, Katie, Annie DeBauw, Dilip K. Lakshman, Daniel P. Roberts, Adnan Ismaiel, and Gary Chastagner. 2022. "Detection and Molecular Phylogenetic-Morphometric Characterization of Rhizoctonia tuliparum, Causal Agent of Gray Bulb Rot of Tulips and Bulbous Iris" Journal of Fungi 8, no. 2: 163. https://doi.org/10.3390/jof8020163

APA Style

Coats, K., DeBauw, A., Lakshman, D. K., Roberts, D. P., Ismaiel, A., & Chastagner, G. (2022). Detection and Molecular Phylogenetic-Morphometric Characterization of Rhizoctonia tuliparum, Causal Agent of Gray Bulb Rot of Tulips and Bulbous Iris. Journal of Fungi, 8(2), 163. https://doi.org/10.3390/jof8020163

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop