INTRODUCTION
Phymatotrichopsis root rot, also known as cotton root rot, Texas root rot, or Ozonium or Phymatotrichum root rot, is an important soilborne disease of over 2,000 dicotyledonous species, including cotton, alfalfa, vegetable crops, and fruit and nut trees (
1). The causal fungus,
Phymatotrichopsis omnivora, is a member of the family Rhizinaceae (Ascomycota: Pezizomycetes) and allied to the genera
Psilopezia and
Rhizina (
2). The disease occurs in most of Texas (excluding the panhandle), southern Oklahoma, and New Mexico, southern and western Arizona, and northern Mexico (
1,
2). Because of its extensive host range,
P. omnivora represents a significant threat to agricultural productivity throughout its geographic range. Phymatotrichopsis root rot is a considerable economic concern, causing up to $100 million in annual losses to the U.S. cotton crop alone (
2). Due to its high economic impact and broad host range,
P. omnivora has been included in lists of regulated organisms by the European and Mediterranean Plant Protection Organization (EPPO A1 list no. 21 [
http://www.eppo.org/QUARANTINE/listA1.htm]), the California Department of Food & Agriculture (no. 3261 in the
CDFA Plant Quarantine Manual, 1989 [
http://pi.cdfa.ca.gov/pqm/manual/pdf/309.pdf]), and the United Nations Security Council's monitoring program of Iraq (S/1995/208 [
http://www.fas.org/news/un/iraq/s/s1995-0208.htm]). Isolation of
P. omnivora from infected plants and infested soils is difficult and often not possible. Hence, culture-independent (molecular) approaches to pathogen identification are desirable. Accurate and sensitive detection and discrimination assays for early diagnosis would facilitate disease management. Methods developed for plant pathogen detection include enzyme-linked immunosorbent assay (ELISA) (
3), loop-mediated isothermal amplification (LAMP) (
4), oligonucleotide array (
5), endpoint PCR (
6), and real-time quantitative PCR (qPCR) (
7). PCR-based techniques are generally more sensitive than immunological methods and have high specificity and discriminatory capabilities. Real-time qPCR offers greater sensitivity and speed than endpoint PCR in the detection of target DNA (
7,
8). TaqMan and SYBR green qPCR are the two most popular qPCR formats and chemistries (
9). SYBR green qPCR detects all amplified double-stranded DNA, including nonspecific reaction products, while TaqMan qPCR detects only specific amplification products because its amplification-dependent cleavage of probes, incorporating reporter and quencher dyes, results in increased fluorescence (
10).
Currently, methods for on-site, rapid, reliable, and sensitive detection of
P. omnivora are not available. On-site accurate and sensitive detection of fungi is a challenge. Available PCR assays require PCR inhibitor-free DNA purification from fungus-infected plants and access to laboratory facilities equipped with PCR and qPCR machines. Simple on-site DNA purification and PCR protocols for accurate pathogen identification would facilitate diagnosis and regulation, as well as disease monitoring and management. Ideally, these protocols would use a portable, battery-operated real-time qPCR platform designed for on-site molecular testing, which allows plant pathogen detection by minimally trained operators in the absence of laboratory facilities and conditions, including electricity, centrifuges, liquid nitrogen, water baths, incubators, and hazardous chemicals. Portable instruments developed previously for on-site pathogen detection include the SmartCycler (Cepheid, Sunnyvale, CA), the R.A.P.I.D. system (Idaho Technologies, Salt Lake City, UT), the LightCycler (Roche Applied Science, Indianapolis, IN), and the Bio-Seeq instrument (Smiths Detection, Edgewood, MD). Recently, Tomlinson and coworkers (
10,
11) demonstrated the on-site detection of
Phytophthora ramorum using qPCR (Cepheid SmartCycler II) and loop-mediated isothermal amplification. Use of the SmartCycler has been successful for on-site detection of the bacterium
Xylella fastidiosa in grapevines (
12) and for diagnosis of foot-and-mouth disease (
13). These methods take about 2 h and detect a single target gene. The Razor Ex BioDetection system (Idaho Technology, Inc., Salt Lake City, UT) was designed originally for military use to identify biological-threat organisms on-site. The Razor Ex BioDetection system offers ready-to-use, freeze-dried reagent pouches and bar code-based PCR cycling program upload. These features make the Razor Ex BioDetection system user-friendly. Matero and coworkers (
14) developed 40-min on-site detection methods based on the Razor Ex for anthrax, brucellosis, tularemia, and plague.
Nucleic acid-based detection of fungi in plant material requires the purification of genomic DNA (
15), which can often be contaminated with phenolics and other PCR inhibitors. Therefore, on-site molecular detection requires both a battery-operated real-time qPCR thermocycler and a rapid and PCR inhibitor-free genomic DNA extraction method.
In this paper, we describe the development of TaqMan qPCR assays that may be performed rapidly, on-site, in about 30 min using the field-deployable Razor Ex BioDetection system protocols for reliable, sensitive, and accurate detection of P. omnivora using three target fungal genes. We also describe a modified magnetic-bead-based method for on-site DNA extraction (∼10 min) from fungus-infected plant roots. These tools provide enhanced investigative capability for applications in plant disease diagnostics and management, pathogen population monitoring, agricultural biosecurity, and microbial forensics.
MATERIALS AND METHODS
Plant materials and P. omnivora genomic DNA.
Asymptomatic (apparently healthy), symptomatic (wilted), and dead plants as well as soil samples were collected from affected alfalfa fields at the Samuel Roberts Noble Foundation in Ardmore, OK, and cotton/alfalfa fields from other areas of southwest Oklahoma (
Table 1). All healthy plants included in the exclusivity panel (
Table 2) were grown in a biosafety level two (BSL-2) greenhouse at the Noble Research Center, Oklahoma State University (OSU), Stillwater, except orange, rose, peach, pecan, grape, and Boston fern, which were collected directly as plant tissue from different sources (
Table 2). The genomic DNA of
P. omnivora isolates (
Table 3) used in the inclusivity panel were received from C. Garzon, Department of Entomology and Plant Pathology, Oklahoma State University. The members of the microbial exclusivity panel are listed in
Table 2.
DNA isolation from plant, microbes, and soil.
Genomic DNA from microbes and plants (
Tables 2 and
3) was purified using a DNeasy Plant Mini Kit (Qiagen, Valencia, CA) according to the manufacturer's instructions. Soil DNA was isolated using a MoBio soil extraction kit (MoBio, Carlsbad, CA). DNA concentrations were determined using a NanoDrop v.2000 spectrophotometer (Thermo Fisher Scientific Inc., Worcester, MA). Control
P. omnivora-free soil was collected from Stillwater, OK, an area in which Phymatotrichopsis root rot has not been reported.
DNA isolation from infected plant roots for Razor Ex BioDetection system application.
A Dynabeads DNA Direct Universal kit (Invitrogen, Carlsbad, CA) was used with modifications to the manufacturer's instructions. The target P. omnivora templates were prepared using 10 to 30 mg infected plant root tissues and 100 to 150 μl Tris-EDTA (TE) buffer (Promega, Madison, WI). Samples were macerated using an Eppendorf tube and pestle without liquid nitrogen. The 40 μl of macerated supernatant and 200 μl of Dynabeads were mixed and incubated for 5 min. Tubes containing the macerated supernatant and Dynabead mix were placed in a magnetic rack until the beads formed a tight pellet, following which the supernatant was carefully removed and discarded, followed by two rinses with wash buffer. Finally, 40 μl of kit suspension buffer was added to the washed beads to release and suspend isolated DNA. These modifications permitted DNA isolation in the field.
Primer and probe design.
Sequences of the three target genes, encoding rRNA internal transcribed spacers (ITS-ribosomal DNA [rDNA]), beta-tubulin, and the second-largest subunit of RNA polymerase II (RPB2), were used to design three primer and probe sets (
Table 4). Multiple sequences of the
P. omnivora genes were retrieved from the NCBI GenBank database (
http://www.ncbi.nlm.nih.gov/). Sequences from
P. omnivora ITS-rDNA GenBank accession numbers
FJ914885,
EF494043,
EF494042,
EF494041,
EF494040,
EF494039,
EF494038,
EF494037,
EF442000,
EF441999,
EF441998,
EF441995,
EF441994,
EF441991,
AY549457,
AY549456,
AY549455, and
AY549454 were aligned and used to design primers PO4 F (forward) and PO4 R (reverse) and probe PO4. Sequences from accession numbers
EF494066,
EF494065,
EF494064,
EF494063,
EF494062, and
EF494061 were used to design the primer/probe set Pobt1, and those of
EF494070,
EF494069,
EF494068, and
EF494067 were used to design the primer/probe set PORPB2-2. Sequences of each gene were aligned using Clustal X v.2 (
16) and examined for conserved regions. Percent identity matrices and nucleotide sequence alignments were constructed using GeneDoc (
17). Each primer and probe set, targeting a different gene, was designed using Primer3 (
18) from a consensus sequence within the target gene of
P. omnivora. Primer thermodynamics, internal structures, and self-dimer formation were examined
in silico with mFold (
19). The specificity was confirmed
in silico by screening the primer and probe sequences with BLASTn on the NCBI GenBank database (
20) (
Table 4). Primers and double-quencher probes 5′ 6-carboxyfluorescein (FAM)/ZEN/3′ Iowa Black FQ (5′ 6-FAM/ZEN/3′ IBFQ) were synthesized by IDT (Integrated DNA Technologies, Inc., Coralville, IA).
Positive controls.
For each primer set, positive controls carrying target gene segments of P. omnivora were generated. Endpoint PCR amplicons were eluted from the agarose gel using Quantum Prep Freeze 'N Squeeze Spin Columns (Bio-Rad, Hercules, CA) and cloned using a TOPO-TA cloning kit (Invitrogen). Plasmid DNA carrying the target sequence for each primer set was purified from overnight bacterial cultures using a QIAprep Spin Miniprep kit (Qiagen). The concentrations of total plasmid DNA were determined using a NanoDrop v.2000 spectrophotometer. Inserts from isolated plasmids were sequenced at the OSU Recombinant DNA/Protein Core Facility using M13F and M13R primers. Amplicon sequences were compared against sequences available on the GenBank nucleotide database using BLASTn.
Endpoint PCR amplification.
Preliminary PCR assays to assess the three primer sets were carried out in 20-μl reaction mixtures containing 10 μl GoTaq Green Master Mix (Promega), 1 μl of each forward and reverse primer (5 μM), 1 μl of DNA template, and 7 μl nuclease-free water (Ambion, Austin, TX) in an Eppendorf thermal cycler (Eppendorf, Hauppauge, NY). The cycling parameters were as follows: initial denaturation for 3 min at 94�C followed by 35 cycles of denaturation at 94�C for 20 s, annealing at 58�C for 20 s, and extension at 72�C for 20 s, followed by a final extension at 72�C for 3 min. Positive (plasmid DNA; carrying the target gene) and negative (nontemplate [water]) controls were included in each PCR amplification. A volume of 20 μl of amplified PCR product was electrophoresed in a 1.5% agarose gel in 1� Tris-acetate-EDTA (TAE) buffer, and amplicon sizes were estimated using 1Kb Plus ladders (Invitrogen).
qPCR amplification.
Amplifications with each primer set were carried out in 20-μl reaction mixtures containing 10 μl of Platinum Quantitative PCR SuperMix-UDG (Invitrogen), 0.8 μl (5 μM) of each forward and reverse primer, 0.8 μl (5 μM) probe, 0.12 μl bovine serum albumin (BSA) (50 mg/ml), 1 μl of template genomic DNA, and 6.48 μl of nuclease-free water (Ambion). Positive (plasmid DNA, carrying the target gene) and negative (nontemplate [water]) controls were included in each round of qPCR amplification, and each reaction was performed in three replicates. Cycling parameters included two initial holds, each for 2 min at 50�C and 95�C, followed by 40 cycles at 95�C for 15 s and 60�C for 60 s. The assays were performed in a Rotor-Gene 6000 thermocycler, and data analysis was done using the Rotor-Gene 6000 series software 1.7 (Built 87) (Corbett Research, Sydney, Australia) with a manual cycle threshold (CT) of 0.2.
qPCR sensitivity and spiked assays.
To determine the detection limits of all three primer and probe sets, four sensitivity assays were performed with each set in the Rotor-Gene 6000 thermocycler. Plasmid or genomic DNA was serially diluted in 10-fold increments and used at 1 fg to 10 ng per reaction mixture. Each mixed assay was done by adding 1 μl (per reaction mixture) of soil or cotton leaf extract (uninfected soil or a healthy cotton leaf was macerated in 1 ml of TE buffer and clarified by a 2-min centrifugation at 14,000 rpm; the supernatant was used to spike the DNA preparation) into serially diluted P. omnivora genomic DNA. Each reaction was performed in three replicates.
Razor Ex BioDetection system amplification.
Empty pouches with regular TaqMan reagent were used in place of freeze-dried reagent pouches. Amplification with each primer set was carried out in 150-μl reaction mixtures containing 75 μl of Platinum Quantitative PCR SuperMix-UDG, 6.0 μl (5 μM) of each forward (biotinylated) and reverse primer, 6.0 μl (5 μM) probe, 4 μl of template DNA, and 53 μl of nuclease-free water. Positive and negative controls were included in each round of Razor Ex BioDetection system amplification. Each reaction was performed in duplicate for infected plant samples due to the limited number (12) of well slots. Cycling parameters included one initial hold for 2 min at 50�C and a first cycle at 94�C for 4 min and 60�C for 15 s, followed by 54 cycles at 91�C for 3 s and 60�C for 15 s. The program was uploaded using a bar code (
Fig. 1). The assays were performed with a Razor Ex BioDetection system.
PCR sensitivity comparison with and without primer modification. (i) Effect of primer biotinylation.
The forward primer for each target gene was labeled with biotin. Two concentrations (100 pg and 10 pg) of plasmid DNA carrying the target gene sequences were used for each set of primers in TaqMan qPCR.
(ii) Effect of the addition of two different 5′ AT-rich flaps to primers.
Three qPCR assays were performed to compare the sensitivity obtained using primers designed (i) with a 5′ AT-rich flap or (ii) without flaps. Two different flaps (5′-AATAAATCATAA-3′ and 5′-AAAATTATTTT-3′) were used at the 5′ position of the primers (
21,
22).
DISCUSSION
We report the development of a reliable, accurate, and sensitive qPCR assay for detection of the plant-pathogenic fungus P. omnivora, using three pairs of primers and probes targeting three different genes, ITS-rDNA, beta-tubulin, and RPB2 genes. These primer and probe sets appear to be ready for use in the field-deployable Razor Ex BioDetection system. ZEN double-quenched probes having the internal ZEN quencher with the traditional 3′ end quencher shortened the distance between the reporter and quencher, leading to highly efficient quenching with significantly less background signal and greater precision (GEN, New Rochelle, NY).
High assay specificity, accuracy, and reliability are critical in biosecurity, quarantine, and microbial forensics applications. The use of a multigene format maximizes reliability, specificity, and broad-range detection within
P. omnivora variants by minimizing the risk of false positives and false negatives because each gene acts as an internal control for the two other genes, which enhances the assay's applicability in microbial forensics. No nonspecific cross-reactivity (false positives) was observed for the three primer and probe sets when tested against plant and microbial exclusivity panels (
Table 2). Positive PCR results with genomic DNA of 13 diverse
P. omnivora isolates collected from different locations confirmed broad-range detection of the target pathogen (
Table 3). All 14 symptomatic, asymptomatic, and dead plant roots were PCR positive when tested by all three primer and probe sets. Thirty-one samples of 52 symptomatic and asymptomatic plant and soil samples from alfalfa and cotton fields were qPCR positive for
P. omnivora using the PO4 primers and probe. A few of the positive samples included asymptomatic or infected/dead weed plants collected from the alfalfa field.
The developed assays are sensitive, detecting as little as 1 fg of plasmid DNA (positive control) and 10 fg (PO4) to 100 fg (Pobt1 and PoRPB2-2) of
P. omnivora genomic DNA. The higher sensitivity achieved with primer set PO4 may be due to the greater number of rDNA copies per genome (50 to 200 copies per fungal genome [
23]) relative to the single or low-copy-number protein-encoding beta-tubulin and RPB2 genes. The high copy number of rDNA allowed detection even from fragmented
P. omnivora DNA. The genome size of
P. omnivora is estimated to be ∼115 Mb (
24). Using a previously published formula (
25), 1 pg = 978 Mb (or 1 Mb = 1.022 � 10
−3 pg), the detection limit of a single genome of
P. omnivora should be around 117.5 fg. However, the PO4 primers and probe detect down to 10 fg. Using SYBR green qPCR, Arif and coworkers (
26) detected purified
P. omnivora genomic DNA at amounts as low as 1 fg. However, SYBR green qPCR can result in false positives due to the ability of SYBR green dye to bind to any double-stranded DNA, including primer dimers, limiting the application of this protocol for forensic applications and quarantine inspections.
Previous work had shown that the addition of an AT-rich sequence at the 5′ terminus of primer sequences could lead to increases in PCR yields, sensitivity, and fluorescence (
21,
22,
27). But in this study, the addition of 5′ AT-rich flaps made no difference in
CT values, indicating that the designed primers and probes already had optimal PCR thermodynamics.
An important capability for biosecurity, quarantine, and microbial forensics applications is the ability to use the procedure in the field. Tomlinson and coworkers (
4,
10) developed on-site, loop-mediated isothermal amplification for phytoplasmas and SmartCycler protocols for
Phytophthora ramorum. The SmartCycler has been used also for
Xylella fastidiosa detection (
12). Razor Ex BioDetection system-based methods are more rapid than those employing the SmartCycler (
10,
12,
14). Matero and coworkers (
14) found no difference in sensitivity or accuracy between the Razor Ex and Applied Biosystem 7300/7500 formats. Recently, Idaho Technology, Inc., announced the development of the Razor CRP BioThreat-X kit, which represents the first U.S. Department of Defense system capable of testing for 10 biothreat agents at one time in the field within 30 min. In our investigation, the Razor Ex BioDetection system protocol took about 30 min from sample preparation (DNA extraction, ∼10 min) to final detection (∼20 to 25 min). The rapid DNA isolation method, unlike many laboratory DNA extraction methods, does not require centrifugation, incubation, organic solvents, or the use of liquid nitrogen for sample homogenization. Commercially available Razor Ex BioDetection system pouches generally contain lyophilized PCR reagents to minimize contamination and eliminate the need for cold storage, but in our experiment we injected regular TaqMan qPCR components into the pouch using a disposable syringe. The primers for the Razor Ex BioDetection system were biotinylated, so as to capture the amplified product using streptavidin magnetic beads for further confirmation, if required. No assay inhibition was observed with biotinylated primers.
Multigene TaqMan qPCR and field-deployable Razor Ex BioDetection system-based detection assays for P. omnivora are rapid, reliable, sensitive, and efficient, and their use can speed phytosanitary diagnostics and pathogen detection. Early detection can help to prevent pathogen dissemination through plant materials during interstate or international commerce. The assays also have potential applications for farm management and monitoring of resistance in plant breeding programs. The modified Dynabeads-based rapid DNA isolation method, validated in this work for P. omnivora, can be applied to other plant pathogens, as well as for field and routine diagnostics.
ACKNOWLEDGMENTS
We thank George Opit and Trenna Blagden (Oklahoma State University) for reviewing the manuscript. We extend our appreciation to Carolyn Young and Eddie Funderburg (Noble Foundation) for providing P. omnivora-infected samples used for validation of the techniques developed in this study. Thanks to Rimvydas Vasaitis (Swedish University of Agricultural Sciences [SLU]) for kindly providing dried R. undulata apothecia for primer validation. We are grateful to Carla Garzon (Oklahoma State University) for providing the genomic DNA of P. omnivora, R. undulate, and P. aphanidermatum isolates and reviewing the manuscript. We thank Ouyang Ping (Oklahoma State University) for technical assistance. We also appreciate the help provided by Astri Wayadande, Mindy James, Shefali Dobhal, and Donna Ria Cassi (Oklahoma State University).
This work was supported by the Oklahoma Agricultural Experiment Station, Oklahoma State University (project number OKL 02773).
The mention of trade names or commercial products in this publication does not imply recommendation or endorsement by Oklahoma State University.