Development and Evaluation of ELISA and qRT-PCR for Identification of Squash vein yellowing virus in Cucurbits
- Craig G. Webster
- William W. Turechek , United States Department of Agriculture–Agricultural Research Service (USDA-ARS), U.S. Horticultural Research Laboratory, Fort Pierce, FL 34945
- Weimin Li , Citrus Research and Education Center, University of Florida, Lake Alfred 33850; and Biotechnology Research Institute, Chinese Academy of Agricultural Sciences, Beijing 100081
- Chandrasekar S. Kousik , USDA-ARS, U.S. Vegetable Laboratory, Charleston, SC 29414
- Scott Adkins , USDA-ARS, U.S. Horticultural Research Laboratory, Fort Pierce, FL
Abstract
Squash vein yellowing virus (SqVYV) causes viral watermelon vine decline. To facilitate detection of SqVYV, enzyme linked-immunosorbent assay (ELISA) and quantitative reverse-transcription polymerase chain reaction (qRT-PCR) diagnostic methods were developed. Both methods were capable of detecting SqVYV in a wide range of cucurbit hosts. ELISA was able to detect virus in infected host tissue diluted to at least 1:2,560, which was sufficient for detection in symptomatic squash and watermelon plants. The qRT-PCR method was capable of reliably detecting as few as 3.4 copies of a cloned fragment of SqVYV genomic RNA with an average cycle threshold (Ct) value of 36.4. The sensitivities and specificities for each detection method were estimated by latent class analysis for a set of inoculated squash and watermelon plants at two sampling scales. The scales were hierarchical, with individual plants representing the upper scale and samples from the plant representing the lower scale. The number of samples per plant varied from 1 to 8, and a plant was diagnosed positive if any of its samples tested positive. For all analyses, a cutoff Ct of 35 was chosen for qRT-PCR, which is approximately 2.5 cycles lower than the lowest Ct value achieved for mock-inoculated plants (presumed to be a false positive). qRT-PCR showed high sensitivities (≥0.99) at both sampling scales for squash and watermelon, whereas the sensitivities for ELISA ranged from 0.58 to 0.76. The specificities for both tests were very similar (≥0.94), with ELISA sometimes outperforming qRT-PCR. These diagnostic methods provide additional tools for the identification of SqVYV and management of SqVYV-induced watermelon vine decline.
Squash vein yellowing virus (SqVYV) is a whitefly-transmitted member of the genus Ipomovirus (Adkins et al. 2007; Li et al. 2008; Webb et al. 2012). Watermelon (Citrullus lanatus var. lanatus) is the most economically important host of SqVYV, in which it causes viral watermelon vine decline. SqVYV was first detected in watermelon in 2003 in Florida. The initial symptoms of infection in watermelon are a leaf wilt and petiole necrosis, which progress rapidly to a complete collapse of infected plants. Another characteristic symptom in watermelon fruit is the development of internal discoloration and necrosis in the rind and degradation of the flesh (Adkins et al. 2013). All other reported hosts of SqVYV are cucurbits, which typically display only mild vein yellowing symptoms on leaves and no fruit symptoms. These symptoms may be transient in appearance in some cucurbits (e.g., cucumber). Wilt and necrosis leading to plant death also occur in a limited number of other vining cucurbits (e.g., Cucurbita maxima) upon infection with SqVYV (Webster et al. 2013).
First found in Florida (Adkins et al. 2007), the known geographic range of SqVYV has increased with findings in Indiana, South Carolina, Georgia, California, Puerto Rico, Guatemala, and Israel (Acevedo et al. 2013; Batuman et al. 2015; Egel and Adkins 2007; Jeyaprakash et al. 2015; Reingold et al. 2016; Webster and Adkins 2012). Despite its fairly large geographical distribution, the genetic diversity of the virus is relatively low, with most Florida isolates and all Indiana, Georgia, and Puerto Rico isolates belonging to one predominant genotype (Webster and Adkins 2012). However, a second genotype has occasionally been detected in Florida (Webster and Adkins 2012), and only this genotype is reported from California, Guatemala, and Israel (Batuman et al. 2015; Jeyaprakash et al. 2015; Reingold et al. 2016). The sudden collapse and wilt of infected watermelon vines is consistently observed upon infection with both SqVYV genotypes. Collectively, this suggests that both SqVYV genotypes were recently introduced to the western hemisphere (Gilbertson et al. 2015; Webster and Adkins 2012).
For hosts such as watermelon and squash, experienced individuals can often correctly diagnose SqVYV infection based on characteristic symptoms such as the rapid progression of foliar wilt or necrosis and discoloration of the rind of watermelon, or vein yellowing in squash. Diagnosis based on rind necrosis is problematic because rind necrosis does not occur in all hosts and it requires destructive sampling of fruit (Webster et al. 2013). Visual diagnosis is further complicated by the frequent occurrence and coinfections with aphid-transmitted Papaya ringspot virus (PRSV) or whitefly-transmitted Cucurbit leaf crumple virus (CuLCrV) and Cucurbit yellow stunting disorder virus (CYSDV) that may be confused with or alter the symptoms of SqVYV infection (Turechek et al. 2010). Asymptomatic SqVYV infections are common in other cucurbit hosts and require laboratory methods to diagnose; for example, smell melon (Adkins et al. 2009), balsam apple (Adkins et al. 2008), or melon (Webster et al. 2013)].
Conventional reverse-transcription polymerase chain reaction (RT-PCR) (Adkins et al. 2008) and tissue-blot nucleic acid hybridization assays (TB) (Turechek et al. 2010) have both been developed for SqVYV detection. Although these methods are useful, there are some limitations. For instance, RT-PCR is currently the most accurate test for SqVYV (i.e., it is both sensitive and specific) but is impractical for large-scale testing of many samples due to its costs (Turechek et al. 2013). Conversely, TB can be used to test large numbers of samples but it is less sensitive and specific than RT-PCR. For some cucurbits, including bottle gourd (Lagenaria siceraria), TB is prone to producing false-positive results (Ling et al. 2013; Turechek et al. 2013; Webster et al. 2013). Both RT-PCR and TB are also capable of detecting both SqVYV genotypes.
In this study, we developed an enzyme linked-immunosorbent assay (ELISA) as a cost-effective, high-throughput system for screening large numbers of samples, and a real-time or quantitative RT-PCR (qRT-PCR) as a highly sensitive method for more targeted surveys for diagnosis of SqVYV. Performances of both assays were quantified by latent class analysis (LCA) using data that represented an array of sample types to generate a wide range of possible test outcomes. Finally, we discuss the potential applicability of these methods to research applications and to management of viral watermelon vine decline.
Materials and Methods
Virus isolates.
Initial development of diagnostic tests compared mock-inoculated and SqVYV-infected tissue from a variety of mature, fruit-bearing cucurbits collected in a previous study (Webster et al. 2013) to determine the likelihood of these tests to cross react with plant proteins (ELISA) or nucleic acids (qRT-PCR) of the various cucurbit species (Table 1). Florida isolates of commonly encountered cucurbit viruses Cucumber mosaic virus (CMV), CuLCrV, CYSDV, PRSV type W (PRSV-W), Tobacco streak virus (TSV), Watermelon leaf mottle virus (WLMV), Watermelon mosaic virus (WMV), Zucchini tigré mosaic virus (ZTMV), and Zucchini yellow mosaic virus (ZYMV) were kindly provided in infected cucurbit tissue by Susan Webb (Department of Entomology and Nematology, University of Florida) and used to determine the likelihood of these viruses to cross react with these SqVYV diagnostic reagents. In addition, 58 field samples of cucurbits were collected from southwest and west-central Florida during this study and used to qualitatively compare conventional and qRT-PCR.
As discussed in greater detail below, a series of experiments was used to further characterize qRT-PCR and ELISA test performances. A single isolate of the predominant Florida SqVYV genotype collected from watermelon in 2007 in Hendry County in southwest Florida (Webster and Adkins 2012) was used for these experiments. It was maintained in squash (C. pepo) ‘Prelude II’ (Seminis Seeds, Oxnard, CA) and watermelon ‘Crimson Sweet’ (Bunton Seed Co., Louisville, KY) by mechanical inoculation using 20 mM sodium phosphate buffer (pH 7.0) containing 0.1% (wt/vol) sodium sulfite and 1% (wt/vol) Celite. Plants were maintained in a temperature-controlled greenhouse under natural lighting with a daytime high temperature of 30°C. Plants were drenched with the insecticide imidacloprid (Admire Pro; Bayer Crop Science, Research Triangle Park, NC) at 2 weeks postseeding and subsequently as needed to prevent whitefly infestation.
Development of ELISA for SqVYV detection.
Prior attempts to develop useful antiserum from partially purified SqVYV virions and to use bacteria for expression of the SqVYV coat protein (CP) for use as an antigen were not successful, thus requiring exploration of alternative antigen sources. The SqVYV CP sequence was analyzed for unique amino acid sequences for synthesis of peptides to be used in antibody production. Two peptides were selected from within the SqVYV CP and synthesized (Biosynthesis, Lewisville, TX) for use as antigens. The first peptide (NH2-AFDLAAKEAEIQKLRDEWDANKPSVTLSPF-OH) corresponds to SqVYV CP amino acids 1 to 30 (GenBank accession number EU259611) and was unique to SqVYV isolates. The second peptide (NH2-EGEEDSDGEDEAKDKGKNVTTPDKKQILKG-OH) corresponds to SqVYV CP amino acids 67 to 96 and is also present in Cucumber vein yellowing virus (CVYV), the most closely related ipomovirus. Both peptides were highly conserved in SqVYV isolates, with the first peptide being identical and the second having at least 80% sequence identity to all known SqVYV isolates.
Individual male New Zealand white rabbits were injected with one or both of the synthetic peptides in an equal volume of Freund’s Complete Adjuvant (Cocalico Biologicals, Reamstown, PA) using recommended protocols to produce three different polyclonal antisera. Booster injections were given to the rabbits with an equal volume of Freund’s Incomplete Adjuvant subcutaneously or intramuscularly at 1 and 2 weeks. Test bleeds were made at 5 and 8 weeks, with an exsanguination bleed at 12 weeks.
All three antisera were tested in multiple formats. Only the antiserum produced from the mixture of the two synthetic peptides was useful, and only in an indirect plate-trapped antigen (PTA)-ELISA format. Briefly, tissue was homogenized in coating buffer (15 mM Na2CO3, 35 mM NaHCO3, and 0.5 mM polyvinyl pyrrolidone) at a ratio of 1:10 (wt/vol), with 100 μl of homogenate added per well to a polystyrene plate and incubated overnight at 4°C. The following day, primary antibody was diluted 1:250 (vol/vol) and 100 μl was added and incubated at 37°C for 3 h, followed by incubation with alkaline phosphatase conjugated goat-antirabbit antiserum (GAR-AP; Bio-Rad, Hercules, CA) diluted 1:3000 (vol/vol). Both primary and secondary antibodies were diluted in 1× phosphate-buffered saline-Tween (1× PBST; 140 mM NaCl, 1 mM KH2PO4, 8 mM Na2HPO4, 2.5 mM KCl, and 0.5% [vol/vol] Tween 20 [pH 7.4]) containing 0.5 mM polyvinylpyrrolidone. Plates were washed three times with 1× PBST between all steps. Plates were developed (100 μl of 1 mg/ml p-nitrophenyl phosphate [Sigma-Aldrich, St. Louis] dissolved in 0.01 M diethanolamine buffer [pH 9.6]) and read using an MRX II plate reader (Dynex, Chantilly, VA), with readings twice that of mock-inoculated plants considered positive.
Development of qRT-PCR for SqVYV detection.
A Taqman assay (Applied Biosystems, Foster City, CA) for real-time detection and quantification of the nuclear inclusion b (NIb) gene of SqVYV was developed. Primers and probes were created to amplify a 133-bp region (nucleotides 7,095 to 7,227 of accession number EU259611) within the NIb gene in the 3′ portion of the SqVYV genome. Primers (SqVYV NIb-F, 5′-GTGAAGTCCAGAAGCAAG-3′; SqVYV NIb-R, 5′-ATGGAAGGTCCATACCTC-3′) and probe (SqVYV NIb-P, 5′-FAM-ATCAAGGGCAAGCGAGCTTCTTTC-BHQ-3′) were designed using an online development tool (PrimerQuest; Integrated DNA Technologies, Coralville, IA; www.idtdna.com) using default parameters and were synthesized by Biosearch Technologies (Novarto, CA).
Total RNA was extracted from 100 mg of plant tissue using Trizol (Invitrogen, Carlsbad, CA). First-strand cDNA was synthesized using a degenerate potyvirus 3′ terminus primer (Gibbs and Mackenzie 1997), as previously described (Adkins et al. 2008). Reactions (20 μl) were prepared in 96-well plates using 1 μl of first-strand cDNA, Express 2× Fast master mix (Invitrogen), and primers and probe at final concentrations of 300 and 100 μM, respectively. The passive reference dye carboxy-X-rhodamine was added to a final concentration of 500 nM. All samples, standards, and controls were tested in duplicate or triplicate on an ABI Prism 7500 Fast thermocycler (Applied Biosystems), Eppendorf Mastercycler ep realplex4 (Eppendorf, New York), or Bio-Rad IQ5 Thermocycler (Bio-Rad) using a program of 95.0°C for 20 s followed by 40 cycles at 95°C for 3 s and 60°C for 30 s, with data collection made in the 60°C for 30 s step.
Minimum detectable concentration and efficiency of qRT-PCR.
Segments of the SqVYV genome containing the region amplified by the qRT-PCR and conventional RT-PCR methods (6,368 to 7,963 and 8,493 to 9,563 nucleotides [nt], respectively, of accession number EU259611) were amplified using primers as previously described (Adkins et al. 2008; Li et al. 2008). Both amplicons were excised and purified (Qiagen, Valencia CA) from a 1% agarose gel, and eluted in 30 μl of elution buffer and quantified on a NanoDrop spectrophotometer (NanoDrop Technologies, Wilmington, DE). The copy number of the PCR products was estimated at 3.4 × 1011 copies/μl (NIb) and 8.3 × 1010 copies/μl (CP) using the lengths 1,596 nt for NIb and 1,071 nt for CP, Avogadro’s number (6.02 × 1023 copies/mole), and average nucleotide weight (660 g/mol) in the calculations. Ten-fold dilutions of each product were made from 10−4 to 10−10 (estimated at 3.4 × 107 to 340 copies/μl for NIb and 8.3 × 106 to 8 copies/μl for CP) and used to construct a standard curve to estimate the efficiency of the assay. Standards containing from 3.4 × 103 to 3.4 × 108 copies/reaction were run on each plate of samples and were used to estimate the viral RNA concentration from cycle threshold (Ct) values for unknown samples.
Construction of datasets and statistical analysis.
Appropriate analysis of a diagnostic tool requires data containing the likely range of possible outcomes collected over a diverse set of conditions. Data should consist of a population of samples that have been exposed to the pathogen, typically referred to as “cases”, and another population that has not, known as the “controls”. Using these as guidelines, four sets of data were used to evaluate the ELISA and qRT-PCR tests for SqVYV. All plants discussed below were mechanically inoculated on a single, mature leaf, as described above. In addition, samples were also tested for SqVYV by TB and conventional RT-PCR, and through simple visual assessment for comparative and analytical purposes, also described below (Adkins et al. 2009; Turechek et al. 2010).
The first dataset was constructed to account for differences that might result from dissimilarities in the spatial distribution and accumulation of virus within plants (Turechek et al. 2010). Eight samples were collected along the length of each plant: the petiole and leaf tissue (squash) or vine segment (watermelon) from the inoculated leaf; the leaf immediately above the inoculated leaf; the youngest, fully expanded leaf; the crown tissue; and the growing tip. Each sample was first visually rated for symptoms of SqVYV infection. Samples were then processed and tested by TB and ELISA, and RNA was extracted for subsequent testing by RT-PCR within a few days of collection. The remaining cDNA was stored frozen until tested by qRT-PCR once this assay was developed, approximately 3 years later. In total, 64 samples (8 each from eight plants) were tested.
The second dataset was collected from 110 plants (55 watermelon and 55 squash) and was constructed to account for differences that might result from dissimilarities in the temporal distribution and accumulation of virus within plants from inoculation to 14 days postinoculation (dpi). However, 15 of the 55 squash plants were used in a preliminary study. At 0 and 1 dpi, the petiole of the inoculated leaf (squash) or the crown (watermelon) were rated for symptoms of SqVYV and then tested for SqVYV by TB, ELISA, RT-PCR, and qRT-PCR as above. At 4, 5, 6, 7, 8, and 14 dpi, two samples per plant—petioles of the inoculated and youngest leaves (squash) or crown and growing tip (watermelon)—were collected, rated, and tested for SqVYV as above.
The third dataset consisted of 26 PRSV-inoculated watermelon plants and 7 PRSV-inoculated squash plants. PRSV commonly infects both watermelon and squash in Florida, and testing of this group of plants was done to measure the effect of PRSV infection, if any, on test performance. Plant samples were first visually examined for symptoms that can be confused as SqVYV and then tested immediately for SqVYV by TB, ELISA, RT-PCR, and qRT-PCR (i.e., cDNA was not stored frozen). For each watermelon plant, the crown, a segment from a noninoculated vine, and the growing tip of the plant were tested. On each squash plant, a petiole segment from the mock-inoculated and the newest and second noninoculated leaves were tested. In total, 99 samples from 33 plants were tested.
The final dataset consisted of 18 mock-inoculated squash plants. Plant samples were rated for symptoms that can be confused as SqVYV and then tested immediately for SqVYV by TB, ELISA, RT-PCR, and qRT-PCR. On each squash plant, a petiole segment from the mock-inoculated and the newest and second noninoculated leaves were tested. In total, 54 samples from 18 plants were tested.
LCA.
Sensitivities, specificities, and positive predictive values for the five SqVYV diagnostic tests were estimated at two sampling scales using LCA: the individual sample scale and the entire plant scale (Turechek et al. 2013). The sensitivity of a test is defined as the conditional probability of a positive test result given that the test subject is infected, the specificity is the conditional probability of a negative test result given that the subject is not infected, and the positive predictive value is the probability of infection for a given set of test results.
In a latent class model, it is assumed that the responses of the manifest variables (i.e., test results) are the result of an unobserved (i.e., latent) structure. In our case, the latent structure is the disease status of the samples or plants. LCA further requires that the manifest variables be categorical variables. Thus, qPCR results were converted to binary values by labeling samples with Ct ≤ 35 as positive for SqVYV; otherwise, the sample was considered negative. With d = 5 diagnostic tests, each giving two possible results, there are 2d = 32 possible response patterns (Table 2).
In the diagnostic setting, the latent class methodology takes advantage of the cross-classified test results and uses a maximum-likelihood approach to assign individuals into one of two mutually exclusive latent classes (infected and noninfected) to estimate the individual test sensitivities and specificities. At the sample scale, test characteristics were determined directly from the test results of the individual samples. At the plant scale, a plant was labeled positive for SqVYV if any of the sample-scale samples originating from that plant tested positive; otherwise, the plant was considered free of infection.
A fuller account of LCA can be found in Turechek et al. (2013). We used PROC LCA in SAS (version 9.3; SAS Institute, Cary, NC) to fit the latent class models by treating disease status as a two-class latent variable (Lanza et al. 2007, 2011). PROC LCA is an add-on procedure for SAS that can be obtained free of charge through Pennsylvania State University’s Methodology Center’s website (https://methodology.psu.edu/ra/lca).
Results
Development of ELISA for SqVYV detection.
Antiserum background and specificity.
Antiserum produced from the mixture of two synthetic peptides proved highly reliable for the detection of SqVYV in a range of host species when tested in an indirect PTA-ELISA format, with no significant cross reactivity observed with mock-inoculated controls (Table 1). The test did not cross react with other common cucurbit viruses, including aphid-transmitted CMV, PRSV-W, WMV, ZTMV, and ZYMV; or whitefly-transmitted CuLCrV and CYSDV (data not shown). For regulatory reasons, we were not able to directly test ELISA on CVYV or the SqVYV genotype typical in California (Batuman et al. 2015), Guatemala (Jeyaprakash et al. 2015), or Israel (Reingold et al. 2016); however, based on alignments of the deduced CP amino acid sequences (>80% identity and minimum of 10 consecutive identical amino acids), our antibodies were predicted to detect this SqVYV genotype. In addition, we provided our SqVYV antiserum to researchers in California who directly demonstrated detection of their SqVYV isolate (Batuman et al. 2015).
Minimum detectable concentration.
When infected plant material was homogenized, diluted, and tested by ELISA, the virus was consistently detectable in watermelon crown and squash leaf samples at dilutions of 1 part in 2,560 (absorbance at 405 nm [A405] of 0.839 and 0.393, respectively; Fig. 1). The virus was often detectable in watermelon crown samples at the next highest dilution of 1 part in 10,240 (A405 of 0.410), and greater absorbance values were observed for infected watermelon crown compared with squash leaf samples (Fig. 1). Absorbance values from mock-inoculated control plant tissue (A405 of 0.122) were similar to those of buffer-only controls (A405 of 0.102), indicating a low level of background noise in the test.
Development of qRT-PCR for SqVYV detection.
PCR efficiency, primer cross reactivity, and performance on field samples.
The qRT-PCR assay was determined to be an efficient assay with a high linear correlation between the logarithm of virus concentration and Ct value, as determined from the standard curve (R2 = 0.997; Fig. 2). Equivalent results were obtained with all three thermocyclers tested (data not shown) but the ABI 7500 machine was used for all experiments presented here. The efficiency of the reaction, as estimated from the slope of the standard curve, was 95.3%. Additionally, an estimated 3.4 copies of the target region per reaction could be detected, which gave, on average, a Ct of 36.4. This detection limit was approximately 2,440 times lower than the existing conventional RT-PCR method (Adkins et al. 2008), which required an estimated 8,300 copies of the SqVYV-CP gene (data not shown).
SqVYV qRT-PCR was tested against the same common cucurbit viruses used in ELISA above and showed no amplification after 40 cycles, except for CYSDV, which gave a Ct of 38.7. However, several mock-inoculated cucurbit species also showed amplification at Ct >37.5 (Table 1); therefore, the cross reaction seen with CYSDV may be partly associated with the cucurbit background rather than the virus. Low-level amplification from noninfected plant samples is common and is the reason why cutoff Ct values are typically set at 35. This is approximately 2.5 Ct lower than the lowest reading in Table 1.
When tested on a collection of 58 cucurbit samples collected from southwest and west-central Florida, 20 samples tested positive by conventional RT-PCR; the remaining 38 tested negative. For qRT-PCR, the range of Ct values observed among the positive RT-PCR samples was 15.5 to 28.0, whereas the range observed among the negative RT-PCR samples was 32.5 to 40.0. Nearly all of the SqVYV-negative samples were coinfected with CuLCrV, CYSDV, or PRSV-W.
Estimation of the sensitivity and specificity of ELISA and qRT-PCR by LCA.
Complete agreement among all five test results—either all positive or all negative—was found in approximately 57% of the sample-scale samples for both squash ([19 all positive + 96 all negative]/198 total) and watermelon ([10 all positive + 109 all negative]/207 total), and in 69% ([24 all positive + 30 all negative]/78 total) and 61% ([8 all positive + 44 all negative]/85 total) of the plant-scale samples for squash and watermelon, respectively (Table 2). There was a greater number of response patterns to consider at the sample scale than the plant scale, which is likely due to the greater number of samples at that scale.
At both the sample and plant scales, the sensitivities and specificities of qRT-PCR and RT-PCR were greater than 0.95 and often greater than 0.97 (Table 3), as expected. For the other three tests, the presence of typical virus symptoms was the least sensitive diagnostic, except on squash at the plant scale, where both ELISA and TB had lower sensitivity. This can be attributed to the highly characteristic vein yellowing that develops with SqVYV-infected squash plants, along with the skill of the raters and, perhaps, knowing that plants were inoculated. ELISA tended to be more sensitive than TB, except on watermelon at the plant scale. The specificities were all greater than 0.91 for the remaining tests, except for TB in watermelon (Table 3). The positive predictive values were greater than 0.95 when two or more tests gave positive results, as long as one of those tests was qRT-PCR (Table 4). The only exception was with TB, where the positive predictive values were 0.79 and 0.61 for squash and watermelon, respectively, at the sample scale.
Discussion
ELISA and qRT-PCR methods for SqVYV detection were developed and validated by determining minimum detectable concentrations, the likelihood to cross react with other common cucurbit viruses or cucurbit host tissue, and through testing of a set of field samples shown previously to contain SqVYV. Furthermore, test performances were quantified in a series of experiments designed to provide the data necessary to allow estimation of test sensitivity and specificity through LCA. qRT-PCR was shown to have superior sensitivity to ELISA, TB, RT-PCR, and symptoms for identification of infected material in a wide range of hosts (Table 3). However, its specificity (i.e., the ability to correctly identify virus-free samples) was similar to RT-PCR and ELISA and on par with a number of the other tests when Ct = 35.0 was selected as the cutoff.
When available, ELISA is often a preferred diagnostic method because it is simple to perform, provides rapid test results, can be used to process large numbers of samples, and is relatively inexpensive. The ELISA developed was less sensitive than both of the RT-PCR methods but its specificity was 0.94 or greater in the applications tested here, meaning that false positives occurred infrequently. Thus, a positive reaction by ELISA is strong indication of SqVYV infection. The ELISA method has distinct advantages over the existing TB method for SqVYV detection. First and foremost, the higher sensitivity and specificity of ELISA over TB in combination with the robust and simple interpretation of results, such as automated spectrophotometric plate reading, make the determination of infected samples with ELISA less subjective than TB. The TB method can be prone to difficulties in interpreting faint signals and nonspecific artifacts of blot development. For example, in previous studies, TB routinely produced a nonspecific signal from bottle gourd germplasm on blots (Ling et al. 2013; Turechek et al. 2013; Webster et al. 2013). In comparison, no nonspecific signal from bottle gourd was noted in ELISA testing in the current study (Table 1). Bottle gourd and other cucurbit species show resistance to SqVYV-induced vine decline and are of interest to breeding programs to introduce SqVYV resistance to watermelon germplasm (Kousik et al. 2009; Webster et al. 2013). Thus, it is essential for SqVYV diagnostic methods to function across a wide selection of cucurbit species.
The qRT-PCR method developed targeted the SqVYV NIb gene, making it independent of our other CP gene-based techniques. With a 40-cycle run, no detection of other cucurbit viruses commonly found in Florida was observed with the qRT-PCR tests. However, nonspecific reactions were occasionally seen with mock-inoculated cucurbits, manifested as high Ct values (Table 1). In plant pathology, qPCR is generally used as a qualitative method (presence or absence), with the cutoff Ct being chosen to reduce the likelihood of encountering such false positives as well as to provide a simple means to distinguish positive from negative samples (Chandelier et al. 2010). In our case, the cutoff Ct was set to 35.0, which is approximately 2.5 cycles lower than the lowest Ct value achieved for the mock-inoculated plants shown in Table 1. Even with this low cutoff, the sensitivity of the test was still 0.99 at both the sample and plant scales for squash and watermelon. The cutoff Ct of 35.0 was sufficient to allow SqVYV detection prior to symptom development in inoculated plants, when viral titer is presumably low, and allowed detection of approximately 3.4 copies of the viral RNA per reaction.
Comparing the performance of qRT-PCR to the RT-PCR method published for SqVYV (Adkins et al. 2008) on a collection of 58 field samples, there were six instances where qRT-PCR identified viral infection (with a cutoff Ct of 35) but RT-PCR did not. For these six tests, the Ct values ranged between 32.5 and 34.7. At a minimum, this represents a 10-fold difference in sensitivity based on the standard curve. Unfortunately, there were no field samples among the 58 with Ct values between 28.0 (the maximum value among RT-PCR-positive samples) and 32.5 (the minimum value among RT-PCR-negative samples) to better compare the two sensitivities. Direct comparison of the two procedures is further complicated by the fact that the procedures target different regions of the SqVYV genome, as discussed above.
Dichotomizing continuous test results is a practice that is not necessarily recommended, largely because of the loss of information, the increased likelihood of misclassification, and the loss of statistical power that occurs with artificial grouping (Dawson and Weiss 2012; Wang and Turechek 2016). Nonetheless, the practice is rampant for many reasons. Perhaps the most obvious is the desire for a yes or no answer, rather than a probability, and the fact that some common analyses (e.g., LCA) do not easily accommodate mixed data types (i.e., nominal, ordinal, binary, and so on). Numerous methods exist for cutoff selection and the topic is too large for discussion here (Caraguel et al. 2011; Chandelier et al. 2010). However, it should be noted that the reason that we did not endeavor to identify an “optimal” cutoff is that selection is largely based on the type of misclassifications one wishes to minimize. Those wishing to minimize false-positive classifications should establish lower values for the cutoff Ct, those wishing to minimize false-negative classifications should establish higher values for the cutoff Ct, and those desiring a balance between the two can identify a cutoff based on one of several indices or graphical methods (Gönen 2007).
In the statistical analysis, it should be clear that the choice of the cutoff Ct will affect the estimates of qRT-PCR sensitivities and specificities; however, that choice also affects the estimates of the other tests because of how LCA uses cross-classified data (Wang and Turechek, 2016). For example, when a cutoff Ct of 36.5 was chosen (rather than 35), the specificity of qRT-PCR at the sample scale for squash decreased from 0.989 to 0.964. This, in turn, altered the sensitivities and specificities of the other tests, typically by less than 0.01, and also lowered the positive predictive values for the five test results (data not shown).
The data were analyzed at two sampling scales that were hierarchical by design. At the plant scale, a plant was rated as positive if any of the plant’s parts (= samples) tested positive. If the assumption is made that plant inoculation is 100% effective, then disease status at the plant scale was presumably known. However, for tests other than symptoms, confirmation of disease status depends on whether the plant sample being tested contains the pathogen. Because vascular-invading pathogens can be nonuniformly distributed in the plant, multiple samples from a plant are often analyzed to determine a plant’s disease status. The largest differences in test performance at the two scales were found between the two tests that required subjective interpretation. For squash, the sensitivity of symptom expression at the sample scale was 41% lower than at the plant scale (Table 3). In general, it is simpler to identify the disease status of a plant when visualizing the entire plant rather than just a sample from that plant. Similarly, the sensitivity of TB was 35% lower at the sample scale than the plant scale for watermelon. This is likely a reflection of a combination of both the lower sensitivity of TB and the uneven distribution of virus in the plant. The specificity of most tests at both scales on both species was greater than 0.94 and often greater than 0.98, with the exception of TB. Again, the propensity for false signal and subjective interpretation are the likely culprits in the reduction in test specificity.
In conclusion, two additional diagnostic methods were developed for identification of SqVYV. The advantages of ELISA are its user friendliness, quick turnaround of test results, and scalability to process large numbers of samples at a reasonably low cost. On the other hand, qRT-PCR can be multiplexed with tests for other common cucurbit-infecting viruses of Florida (e.g., CYSDV), as has been done for hosts commonly infected by several viruses (Mortimer-Jones et al. 2009). Alternatively, the assay could be easily modified for detection of SqVYV within the whitefly vector, as was done for CVYV (Gil-Salas et al. 2007), which would be useful for monitoring viruliferous vector populations as a part of efforts to predict virus outbreaks. Both tests were shown to be accurate and robust and, thus, provide additional tools for the detection and management of SqVYV-induced watermelon vine decline.
The choice of which test to use will depend on the purpose (routine diagnostics or experimental), number of samples (few or many), and type of results desired (qualitative or quantitative). For routine SqVYV diagnostic purposes, we use conventional RT-PCR for small numbers of samples and ELISA for large numbers of samples. Conventional RT-PCR and ELISA both detect all known genotypes of SqVYV, as does TB. In contrast, qRT-PCR targets the less-conserved SqVYV NIb gene and is most useful for the predominant Florida SqVYV genotype. Accordingly, qRT-PCR offers a means to distinguish SqVYV genotypes without sequencing. A major use for qRT-PCR is in controlled experiments with the predominant Florida SqVYV genotype.
Acknowledgments
We thank S. Webb for providing cucurbit viruses for use in specificity testing; C. Vanderspool, B. Burns, and N. Miller for excellent technical assistance; and M. Hilf, K. Sims, and C. Latza for helpful advice and useful discussions on qRT-PCR. Financial support was provided, in part, by the United States Department of Agriculture Cooperative State Research, Education, and Extension Service Specialty Crop Research Initiative grant 2008-04890 and the National Watermelon Association.
Literature Cited
- 2013. First report of Squash vein yellowing virus affecting watermelon and bitter gourd in Puerto Rico. Plant Dis. 97:1516. https://doi.org/10.1094/PDIS-03-13-0322-PDN LinkWeb of ScienceGoogle Scholar
- 2013. Physiological effects of Squash vein yellowing virus infection on watermelon. Plant Dis. 97:1137-1148. https://doi.org/10.1094/PDIS-01-13-0075-RE LinkWeb of ScienceGoogle Scholar
- 2007. Identification and characterization of a novel whitefly-transmitted member of the family Potyviridae isolated from cucurbits in Florida. Phytopathology 97:145-154. https://doi.org/10.1094/PHYTO-97-2-0145 LinkWeb of ScienceGoogle Scholar
- 2008. Squash vein yellowing virus detection using nested polymerase chain reaction demonstrates that the cucurbit weed Momordica charantia is a reservoir host. Plant Dis. 92:1119-1123. https://doi.org/10.1094/PDIS-92-7-1119 LinkWeb of ScienceGoogle Scholar
- 2009. Detection of three whitefly-transmitted viruses infecting the cucurbit weed Cucumis melo var. dudaim in Florida. Plant Health Prog. doi:10.1094/PHP-2009-1118-01-BR. LinkGoogle Scholar
- 2015. First report of an ipomovirus infecting cucurbits in the Imperial Valley of California. Plant Dis. 99:1042. https://doi.org/10.1094/PDIS-12-14-1248-PDN LinkWeb of ScienceGoogle Scholar
- 2011. Selection of a cutoff value for real-time polymerase chain reaction results to fit a diagnostic purpose: Analytical and epidemiological approaches. J. Vet. Diagn. Invest. 23:2-15. https://doi.org/10.1177/104063871102300102 CrossrefWeb of ScienceGoogle Scholar
- 2010. Determination of cycle cut off in real-time PCR for the detection of regulated plant pathogens. Bull. OEPP/EPPO Bull. 40:52-58. https://doi.org/10.1111/j.1365-2338.2009.02354.x CrossrefGoogle Scholar
- 2012. Dichotomizing continuous variables in statistical analysis: A practice to avoid. Med. Decis. Making 32:225-226. https://doi.org/10.1177/0272989X12437605 CrossrefWeb of ScienceGoogle Scholar
- 2007. Squash vein yellowing virus identified in watermelon (Citrullus lanatus) in Indiana. Plant Dis. 91:1056. https://doi.org/10.1094/PDIS-91-8-1056B LinkWeb of ScienceGoogle Scholar
- 1997. A primer pair for amplifying part of the genome of all potyvirids by RT-PCR. J. Virol. Methods 63:9-16. https://doi.org/10.1016/S0166-0934(96)02103-9 CrossrefWeb of ScienceGoogle Scholar
- 2015. Role of the insect supervectors Bemisia tabaci and Frankliniella occidentalis in the emergence and global spread of plant viruses. Annu. Rev. Virol. 2:67-93. https://doi.org/10.1146/annurev-virology-031413-085410 CrossrefWeb of ScienceGoogle Scholar
- 2007. Development of a real-time RT-PCR assay for the detection of Cucumber vein yellowing virus (CVYV) and Cucurbit yellow stunting disorder virus (CYSDV) in the whitefly vector Bemisia tabaci. J. Virol. Methods 146:45-51. https://doi.org/10.1016/j.jviromet.2007.05.032 CrossrefWeb of ScienceGoogle Scholar
- 2007. Analyzing Receiver Operating Characteristic Curves with SAS. SAS Institute Inc., Cary, NC. Google Scholar
- 2015. First report of Squash vein yellowing virus in watermelon in Guatemala. Plant Health Prog. 16:113-114. LinkGoogle Scholar
- 2009. Sources of resistance in U.S. plant introductions to watermelon vine decline caused by Squash vein yellowing virus. HortScience 44:1-7. CrossrefWeb of ScienceGoogle Scholar
- 2007. PROC LCA: A SAS procedure for latent class analysis. Struct. Equation Model. 14:671-694. https://doi.org/10.1080/10705510701575602 CrossrefWeb of ScienceGoogle Scholar
- 2011. Proc LCA & Proc LTA Users’ Guide (Version 1.2.6). The Methodology Center, Pennsylvania State University, University Park. https://methodology.psu.edu Google Scholar
- 2008. Presence of P1b and absence of HC-Pro in Squash vein yellowing virus suggests a general feature of the genus Ipomovirus in the family Potyviridae. Virus Res. 135:213-219. https://doi.org/10.1016/j.virusres.2008.03.015 CrossrefWeb of ScienceGoogle Scholar
- 2013. Development and field evaluation of multiple virus-resistant bottle gourd (Lagenaria siceraria). Plant Dis. 97:1057-1062. https://doi.org/10.1094/PDIS-07-12-0639-RE LinkWeb of ScienceGoogle Scholar
- 2009. A single tube, quantitative real-time PCR assays that detects four potato viruses simultaneously. J. Virol. Methods 161:289-296. https://doi.org/10.1016/j.jviromet.2009.06.027 CrossrefWeb of ScienceGoogle Scholar
- 2016. Watermelon fruit rot disease in Israel is caused by a distinct Squash vein yellowing virus (SqVYV) strain. Plant Dis. 100:1176-1183. https://doi.org/10.1094/PDIS-09-15-1040-RE LinkWeb of ScienceGoogle Scholar
- 2010. Distribution of four viruses in single and mixed infections within infected watermelon plants in Florida. Phytopathology 100:1194-1203. https://doi.org/10.1094/PHYTO-01-10-0018 LinkWeb of ScienceGoogle Scholar
- 2013. The use of latent class analysis to estimate the sensitivities and specificities of diagnostic tests for Squash vein yellowing virus in cucurbit species when there is no gold standard. Phytopathology 103:1243-1251. https://doi.org/10.1094/PHYTO-03-13-0071-R LinkWeb of ScienceGoogle Scholar
- 2016. A loop-mediated isothermal amplification assay and sample preparation procedure for sensitive detection of Xanthomonas fragariae in strawberry. PLoS One 11:e0147122. https://doi.org/10.1371/journal.pone.0147122 CrossrefWeb of ScienceGoogle Scholar
- 2012. Semipersistent whitefly transmission of Squash vein yellowing virus, causal agent of viral watermelon vine decline. Plant Dis. 96:839-844. https://doi.org/10.1094/PDIS-09-11-0761 LinkWeb of ScienceGoogle Scholar
- 2012. Low genetic diversity of Squash vein yellowing virus in wild and cultivated cucurbits in the U.S. suggests a recent introduction. Virus Res. 163:520-527. https://doi.org/10.1016/j.virusres.2011.11.017 CrossrefWeb of ScienceGoogle Scholar
- 2013. Squash vein yellowing virus infection of vining cucurbits and the vine decline response. Plant Dis. 97:1149-1157. https://doi.org/10.1094/PDIS-01-13-0076-RE LinkWeb of ScienceGoogle Scholar
Disclaimer: the use of trade, firm or corporation names in this publication is for the information and convenience of the reader. Such use does not constitute an official endorsement or approval by the United States Department of Agriculture or the Agricultural Research Service of any product or service to the exclusion of others that may be suitable.