Previous Article | Next Article ![]()
Journal of Clinical Microbiology, April 2004, p. 1511-1518, Vol. 42, No. 4
0095-1137/04/$08.00+0 DOI: 10.1128/JCM.42.4.1511-1518.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Department of Microbiology and Immunology, The University of Oklahoma Health Sciences Center, Oklahoma City, Oklahoma 73104
Received 11 October 2003/ Returned for modification 26 November 2003/ Accepted 10 January 2004
|
|
|---|
|
|
|---|
Originally isolated from the blood of a febrile Ugandan woman in 1937, WNV is widely distributed throughout Africa, the Middle East, areas of Europe, and Asia (31). The summer of 1999 marked the first incidence of WNV on the North American continent, in the northeastern United States (23). During this outbreak, more than 60 people became clinically ill, and seven deaths occurred due to encephalitis (12). WNV has continued to cause seasonal epidemics and to spread westward. In 2002, the virus was detected in 43 of 50 states, and more than 3,800 human cases were reported; 225 of these were fatal (8). Currently, no specific drug or vaccine against WNV is approved for human use, so that palliative care, nonspecific antiviral treatment (alpha interferon and ribavirin), surveillance, and screening of potentially contaminated products (such as blood) are the only available measures against WNV.
At the present time, enzyme-linked immunosorbent assays and PCR assays are used to monitor the presence of WNV (17, 22, 28). A real-time quantitative reverse transcriptase PCR (RT-QPCR) assay for WNV has been reported to enable routine high-throughput screening for WNV (19, 22). This assay is based on TaqMan technology; in addition to two specific outer primers, a third, fluorescence-labeled oligonucleotide (TaqMan) is used for detection. This assay has been proven to be as sensitive as gel-based reverse transcriptase PCR (RT-PCR), with a detection limit of 0.1 PFU/ml.
We have developed an alternative real-time RT-QPCR assay with the intercalating dye SYBR green. This assay is less expensive than TaqMan or beacon-based real-time quantitative PCR but is still faster than gel-based, single, or nested RT-PCR and should be economical for large-scale routine testing of clinical samples and blood products. Importantly, we found that this new assay is insensitive to nucleotide variations within the amplified region and thus has a substantially lower false-negative rate than prior assays. Furthermore, through dissociation curve analysis, the SYBR green-based RT-QPCR assay allows for the identification of novel WNV strains.
|
|
|---|
RNA isolation and reverse transcription. RNA was isolated from the homogenate by using an Absolutely RNA Micro-Prep kit (Stratagene, La Jolla, Calif.) according to the manufacturer's instructions. Reverse transcription was performed as previously described (13). Briefly, 500 ng of RNA was reverse transcribed in a 20-µl reaction mixture with 100 U of SupercriptII reverse transcriptase (Invitrogen Inc., Carlsbad, Calif.), 2 mM deoxyribonucleoside triphosphates, 2.5 mM MgCl2, 1 U of RNasin (the last three from Applied Biosystems, Foster City, Calif.), and 0.5 µg of random hexanucleotide primers (Amersham Pharmacia Biotech, Piscataway, N.J.). The reaction mixture was sequentially incubated at 42°C for 45 min, 52°C for 30 min, and 70°C for 10 min. The reverse transcription reaction was stopped by heating to 95°C for 5 min. Next, 0.5 U of RNase H (Invitrogen) was added, and the reaction mixture was incubated at 37°C for an additional 30 min. Afterward, the cDNA pool was diluted 25-fold with diethyl pyrocarbonate-treated, distilled H2O and stored at 80°C.
Real-time quantitative PCR. The primers used for WNV in this study were previously described by Lanciotti et al. (22). They are 5'-TCA CGC ATC TCT CCA CCA AAG-3' and 5'-GGG TCA GCA CGT TTG TCA TTG-3' and are specific for the WNV env region. Positive control RNA (strain WN-NY99) was obtained from the Centers for Disease Control and Prevention Reference Laboratory (Ft. Collins, Colo.). Ninety-six variant 60-mer oligonucleotides were synthesized at The University of Oklahoma Health Sciences Center genomics facility. Real-time PCR was conducted with previously established procedures (11, 13). The final PCR mixture contained 2.5 µl each of forward and reverse primers (final concentration of each, 166 nM), 7.5 µl of 2x SYBR PCR mix (Applied Biosystems), and 5 µl of sample. The PCR was set up in a dedicated room with a CAS-1200 pipetting robot (Pheonix Research, Hayward, Calif.). The CAS-1200 robot uses filtered carbon-graphite pipette tips (Tecan Inc.) for liquid level sensing, allowing for a pipetting accuracy of 0.1 µl and eliminating carryover contamination. All surfaces were cleaned with 10% bleach and exposed to UV light overnight on a daily basis. Gowns, gloves, face masks, and equipment were required for all work. Real-time PCR was performed with an ABI Prism 5700 or ABI Prism 7700 machine (Applied Biosystems) and universal cycling conditions (2 min at 50°C, 10 min at 95°C, 40 cycles of 15 s at 95°C, and 1 min at 60°C). Cycle threshold (CT) values were determined by automated threshold analysis with ABI Prism version 1.0 software. The amplification efficiencies were determined by serial dilution and calculated as E = exp1/m, where E is the amplification efficiency and m is the slope of the dilution curve. Dissociation curves were recorded after each run, and the amplified products were visualized by 2% agarose gel electrophoresis. The entire procedure from sample receipt to data analysis was routinely completed in 4 h.
Statistical analysis. Calculations were performed with Excel version 10.1 (Microsoft Inc., Redwood, Wash.) and SPSS version 12.0 (SPSS, Chicago, Ill.). Analysis of dissociation profiles was performed with ABI Prism version 1.0 software. The distance between dissociation profiles was calculated as a weighted sum of squares of differences between pairs of normalized data points. When necessary, linear interpolation was used between neighboring data points (10). Initially, dissociation profiles from 96 identical samples were compared to establish a detection threshold and overall measurement error. In later comparisons, all distances above the threshold were reported to be different from the reference.
|
|
|---|
To compare the sensitivity of the SYBR green-based assay to that of the TaqMan assay, serial dilutions of strain WN-NY99 were tested in both assays with the WNV env probe and primer pairs previously described by Lanciotti et al. (22). Both assays performed with equal sensitivities and had similar linear ranges (Fig. 1). CT values for both assays were similar, with y-axis intercept points at the greatest dilution of 41.23 ± 3.26 (mean and standard error [SE]) and 40.72 ± 2.99 for TaqMan-based detection and SYBR green-based detection, respectively. The amplification efficiencies were identical, as evidenced by the identical slopes of the regression lines. The regression coefficients (R2) for the TaqMan assay and the SYBR green-based assay were 0.9688 and 0.9625, respectively (Fig. 1). The parallel testing of these two detection methods establishes that the SYBR green-based assay is as sensitive as the TaqMan assay with the same (outer) PCR primers.
![]() View larger version (15K): [in a new window] |
FIG. 1. Linearity and sensitivity comparisons of TaqMan and SYBR green-based real-time quantitative PCRs. Linear regression plots were generated for TaqMan and SYBR green-based real-time quantitative PCRs of WN-NY99 dilutions with the same WNV env primer set. The sample dilution on the horizontal axis is plotted against the CT on the vertical axis (n = 3). (A) Linear regression of TaqMan assay. (B) Linear regression of SYBR green-based assay.
|
![]() View larger version (25K): [in a new window] |
FIG. 2. Detection of WN virus in mosquito pools by SYBR green-based real-time RT-QPCR. (A) Plot of amplification of mosquito pools with the WNV env primer set. Cycle number on the horizontal axis is plotted against the relative fluorescence (Rn) on the vertical axis (log scale). Dotted lines indicate the two positive control reactions. Positive and negative mosquito pools are indicated by plus and minus signs, respectively. (B) Agarose (2%) gel electrophoresis of amplification products stained with ethidium bromide. Lane M, molecular weight markers. Lanes +, positive control (in duplicate). Positive reactions showed a single amplification product of approximately 60 bp. (C) Dissociation plot of amplification products for the positive control and mosquito pools. Temperature on the x axis is plotted against the first derivative of the measured fluorescence [d(F)] on the y axis. The asterisk indicates the peak position of the reference dissociation curve.
|
SYBR green-based real-time RT-QPCR can detect polymorphic targets, a property which lowers the overall false-negative rate. In screening natural WNV isolates, we encountered samples with a substantial viral load (as evidenced by low CT values) but divergent dissociation profiles. These samples likely represented naturally occurring WNV variants with sequence variations within the amplicon and thus may have escaped detection by the TaqMan assay, yielding false-negative results. These samples yielded products which were indistinguishable from those yielded by the WN-NY99 sample (used as a positive control) in gel electrophoresis, since the readout in conventional or nested PCR is a band on a 2% agarose gel that can differentiate only gross variations in size (>5 bp, depending on the gel system) and not nucleotide substitutions. In contrast, the SYBR green-based assay provides an amplicon dissociation profile as an added measure of specificity. With this method, even single nucleotide changes can yield an altered profile, which should allow for the rapid routine detection of sequence variations.
To compare directly the ability of the TaqMan assay and the SYBR green-based assay to detect polymorphisms, we synthesized 60-mer oligonucleotides containing every possible single point substitution in the central 28-bp amplicon region as well as selected double and triple mutations (Table 1). These oligonucleotides were each used as templates in both assays. As expected, all targets yielded a signal with SYBR green, but only 47 (53%) of 88 single-point-substitution amplicons were detected by the probe-based assay (Fig. 3A and B). This finding represents a false-negative rate of 47% for the TaqMan assay for single nucleotide polymorphisms across the probe-binding site. Of the eight double- and triple-substitution amplicons in our experiment, only five (63%) were detected by the TaqMan assay, even though a PCR product was amplified in all instances, as determined by agarose gel electrophoresis (Fig. 3C). This outcome indicates that WNV variants with multiple substitutions will escape detection by TaqMan real-time quantitative PCR with a high probability, as would insertion or deletion mutants.
|
View this table: [in a new window] |
TABLE 1. Sequences of 60-mer mutant oligonucleotides that served as targets for amplification with WN env primers
|
![]() View larger version (32K): [in a new window] |
FIG. 3. Comparison of the effects of nucleotide variations on the TaqMan and SYBR green-based assays. (A) Plot of amplification of 96 mutant oligonucleotides in the SYBR green-based assay. The cycle number on the horizontal axis is plotted against the relative measured fluorescence (Rn) on the vertical axis. (B) CT values obtained for individual mutant target oligonucleotides in either the SYBR green-based assay (horizontal axis) or the TaqMan assay (vertical axis). Samples that did not yield any signal were assigned a CT value of 45. (C) Agarose (2%) gel electrophoresis of the amplification products from the TaqMan assay stained with ethidium bromide. The asterisk indicates an empty lane. (D) Dissociation plot of 96 identical samples (WN-NY99) after amplification in the SYBR green-based assay with the WNV env primer set. (E) Dissociation plot of the 96 mutant oligonucleotides listed in Table 1 after amplification in the SYBR green based assay with the WNV env primer set. The derivative of the measured fluorescence [d(F)] on the vertical axis is plotted against the temperature at which it was sampled on the horizontal axis for panels D and E.
|
G mutations (which are also transversions), 75% of which could be detected by melting curve analysis. No T
A or C
G mutation could be detected, although we are currently optimizing the pattern recognition software (W. Vahrson and D. P. Dittmer, unpublished data). To assess the potential contributions of variations within the instrument to differences in Tms, 96 reactions with WN-NY99 as the template were performed (Fig. 3D). Little to no variation (mean Tm, 83.9°C; SD, ±0.1°C; n = 96) was observed, indicating that differences in Tms are due to differences in sequences and not to machine variations. To determine how close the real-time PCR target area in WNV is to those in other flaviviruses or any other sequence, we performed a blastn (1) search of GenBank (Table 2). Within the first 100 hits based on sequence identity, 78 entries belonged to WNV sequences. Of these, 76 differed by no more than two nucleotides (96%) from the real-time quantitative PCR target sequence. The remaining two WNV sequences showed identity for 25 nucleotides (42%) and 27 nucleotides (45%). The next closest sequences were for three isolates classified as Kunjin virus, with 46 identical nucleotides (76%). JE virus isolates exhibited sequence identity in 20 nucleotide positions (33%) on the right side of this region only and thus would not have been recognized by the left-hand PCR primer. All other flaviviruses exhibited no significant sequence similarity at the nucleotide level.
|
View this table: [in a new window] |
TABLE 2. Flavivirus sequences with significant similarity in the PCR target regiona
|
![]() View larger version (20K): [in a new window] |
FIG. 4. Titration of WNV from an avian sample in a Vero cell culture. (A) Plot of amplification of serially diluted supernatants from infected Vero cell cultures. The cycle number on the horizontal axis is plotted against the relative measured fluorescence (Rn) on the vertical axis. (B) Dissociation plot of amplification products. The temperature on the horizontal axis is plotted against the derivative of the measured fluorescence [d(F)] on the vertical axis.
|
G substitutions were detected, while no T
A substitutions could be detected. It recognized 50% of double nucleotide substitutions in the target sequence and 100% of all triple nucleotide substitutions. It is clear that dissociation profile analysis does not substitute for single-nucleotide-polymorphism analysis for the detection of specific mutants or defined genotypes (14, 16, 18, 30). However, in the context of routine screening, samples with an aberrant dissociation profile indicate variant WNV strains, which then may be selected for detailed characterization by sequence analysis. |
|
|---|
0.1 PFU/ml). We found that real-time quantitative PCR with SYBR green also has a comparable performance profile. Comparisons at the extreme lower end of the linear range between different amplification methods are difficult to interpret, since different primer pairs with associated different amplification efficiencies (27) and different reagents (13) have been used. Previous studies estimated the limit of detection for the probe-based TaqMan RT-QPCR assay for WNV to be
0.1 PFU/ml (2, 22, 26, 28), well below the detection limit for culture-based assays. Using serial dilutions on Vero cells (Fig. 4), we confirmed these observations for SYBR green-based real-time quantitative PCR. The possibility of false-negative test results poses a substantial problem in diagnostic TaqMan assays because mutations within the probe-binding site can prevent annealing of the probe and subsequent detection. While this problem has been documented for herpes simplex virus (3), it is a particular concern for WNV, severe acute respiratory syndrome virus, human immunodeficiency virus, hepatitis C virus, and other RNA viruses, since they exhibit higher sequence variability than DNA viruses and since it is not always possible to identify regions of the genome which are highly conserved. By synthesizing and testing every possible point mutation in the target region for the widely used WNV TaqMan assay, we encountered a potential false-negative rate of 47% for the TaqMan assay due to sequence variations in the probe-binding site. In comparison, the SYBR green-based assay produced a false-positive rate of 0% (Fig. 3). We did not explore the effects of mutations within the binding sites for the forward and reverse primers, since these would affect the TaqMan assay and the SYBR green-based assay (as well as gel-based PCR assays) in the same way. The high-stringency annealing temperature used in real-time quantitative PCR (60°C) for the two outer primers leaves open the possibility that viral variants with mutations in the forward or reverse primer-binding sites may escape detection. Although not tested here, it is likely that the frequency of such false-negatives can be reduced by the use of multiple primer sets for each virus.
Could a SYBR green-based real-time PCR assay be used to screen for novel WNV isolates? A number of studies investigated WNV sequence variability (2, 4, 20, 23, 29). For this flavivirus, the average sequence variability did not exceed 3% at the nucleotide level, and the majority of mutations did not lead to amino acid substitutions. For the 1999 outbreak in New York, Connecticut, and New Jersey, 99% sequence identity was reported over a 1,278-nucleotide region in 13 isolates from avian species, humans, and mosquitoes. Many of the same mutations were shared among isolates from mosquitoes, birds, and humans. Presumably, the need to enter and replicate in cells of avian as well as arthropod origins places severe constraints on the variability of this virus. In contrast, human retroviruses, such as human immunodeficiency virus type 1 (15), exhibit amino acid variability of up to 20% within the same subtype. We found that a preliminary measure of natural variability can be obtained by Tm or dissociation profile analysis (Fig. 3). SYBR green-based real-time quantitative PCR incorporates dissociation profile analysis for each sample at no extra cost. This feature constitutes a considerable improvement over gel-based or TaqMan-based PCR, neither of which can identify the presence of single nucleotide mutations or small (<5-bp) insertions or deletions. In the TaqMan assay, these alterations led to a complete loss of signal (false-negative) or produced signals that were indiscernible from the wild-type signal (Fig. 3); in the gel-based assay, the signals produced by these mutations were indiscernible from the wild-type signal. Since real-time quantitative PCR amplicons are small (<100 bp), even a single base-pair change can result in a distinguishable change in the Tm (7), although in this study only changes of >3 nucleotides were identified with 100% accuracy. Currently, we are trying to improve upon the bioinformatics tools for this analysis. While not as reliable as sequence analysis, dissociation profile analysis has the potential to be used for rapid initial screening for the identification of viral mutants.
This work was supported by NIH grants EB00983 and CA97951 to D.P.D. J.F.P. was supported by NIH training grant T32 AI07364 to the Department of Microbiology and Immunology, The University of Oklahoma Health Sciences Center.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»