A sensitive real-time PCR for detection and subgrouping of human respiratory syncytial virus

Highlights ► A novel real-time RT-PCR using specific locked nucleic acid probes is described. ► The assay is quantitative and distinguishes RSV subgroup A & B. ► Compared with a commercial multiplex PCR using 264 respiratory samplesin Vietnam. ► Sensitivity was significantly higher with detection a rate of 32 vs. 24%.


a b s t r a c t
Improved diagnostic tools for rapid detection, quantitation, and subgrouping of human respiratory syncytial virus (RSV) are needed to aid the development and evaluation of novel intervention strategies. A quantitative real-time RT-PCR using specific locked nucleic acid (LNA) probes was developed to identify RSV and to distinguish RSV subgroups A and B (RSV LNA assay). RSV subgroup diversity and the relationship between viral load and disease severity in confirmed RSV infections were also explored. 264 archived respiratory specimens from pediatric patients were tested in parallel using the commercial multiplex Seeplex TM RV detection kit (Seegene) and the novel RSV LNA assay. The LNA assay demonstrated a significantly higher sensitivity than Seeplex, improving overall detection rates from 24% (64/264) to 32% (84/264). Detection limits of 9.0 × 10 1 and 6.0 × 10 2 copies/mL were observed for RSV A and B, respectively. RSV A was detected in 53/84 (63%) cases, and 31/84 (37%) were positive for RSV B. This novel method offers a rapid, quantitative, highly specific and sensitive approach to laboratory diagnosis of RSV.

Introduction
Human respiratory syncytial virus (RSV) is a Pneumovirus and a member of the subfamily of Pneumovirinae within the family of Paramyxoviridae. Based on differences in the antigenicity and nucleotide sequences analysis, RSV strains have been classified in two major antigenic subgroups, denoted RSV A and RSV B (Sullender, 2000). RSV is a major cause of acute respiratory tract infections in children worldwide (Nair et al., 2010;Weber et al., 1998), and is responsible for a large proportion of hospital admissions during infancy and early childhood (Mejias et al., 2005). No vaccine is available, and treatment is limited to supportive care (Empey et al., 2010). RSV is also a frequent cause of nosocomial infections, and as such may significantly impact the duration of stay in hospital and healthcare costs during childhood (Macartney et al., 2000). Data on epidemiology and molecular virological characteristics of RSV isolates circulating in Southeast Asia are limited.
Diagnosis of RSV infection is typically performed on respiratory samples using virus culture, (rapid) antigen detection tests, or molecular assays (Henrickson, 2004). Real-time reverse transcriptase PCR (RT-PCR) provides a rapid and sensitive tool for detection of RSV compared to conventional techniques and enables (semi-)quantitation of viral load (Borg et al., 2003;Dewhurst-Maridor et al., 2004;do Nascimento et al., 2008;Empey et al., 2010;Falsey et al., 2003Falsey et al., , 2002Henrickson, 2004;Houben et al., 2010;Hu et al., 2003;Kuypers et al., 2004;Mentel et al., 2003;Perkins et al., 2005;van de Pol et al., 2006van de Pol et al., , 2010van Elden et al., 2003). The locked nucleic acid (LNA) technology increases thermal stability of probe-target interaction and improves discrimination of different nucleic acid targets (Reynisson et al., 2006). The N gene is the most conserved RSV gene across all known genetic clades, and is also expressed most abundantly during viral replication (David and Knipe, 2001;Johnson and Collins, 1989). For these reasons, the N gene was chosen as the target of a novel real-time RT-PCR diagnostic tool, to provide quantitative viral load analysis and discrimination between RSV A and B within a single assay. The quantitative assay was compared to a commercial conventional multiplex PCR method (Seeplex TM RV detection kit, Seegene, Inc., Seoul, Korea) (Kim et al., 2009;Roh et al., 2008)  respiratory samples from a study (Do et al., 2011) on acute respiratory infection in children at the Hospital for Tropical Diseases, Ho Chi Minh City, Vietnam. Here, the development and validation of this assay is reported, as well as initial observations regarding subgroup diversity within Vietnamese RSV strains, and relationships between viral load and disease severity in Vietnamese patients.

Clinical specimens
To develop and evaluate this novel RSV diagnostic, 264 respiratory samples, collected and stored as part of a study on acute respiratory infection among 309 children at the Hospital for Tropical Diseases, Ho Chi Minh City, Vietnam were used (Do et al., 2011). Specimens were kept at 4 • C for a maximum of 24 h and then aliquoted and stored at −80 • C until further processing. The samples comprised 240 nasopharyngeal aspirates, 23 throat swabs, and 1 nasal swab collected in viral transport medium, and represented cases admitted to the pediatric intensive care unit (ICU) (61/264, 23%) as well as the pediatric respiratory ward (203/264, 77%). One hundred and forty-one (53%) samples were from patients less than 2 years old, 113 (42%) from 2 to 5 years old, 11 (4%) over 5 years old.

RNA extraction and cDNA synthesis
RNA was extracted from 100 L of specimen and eluted in 60 L on the Easy MAG 2.0 system (bioMérieux, Marcy l'Étoile, France). An internal RNA-virus control (equine arteritis virus [EAV]) was added to each sample prior to extraction at a standard concentration yielding Ct values of 30-35 as previously described (Scheltinga et al., 2005). RNA was reverse transcribed using Superscript III reverse transcriptase (Invitrogen, Carlsbad (CA), USA) and random hexamers (Roche, Mannheim, Germany). Each 20 L reaction mixture contained 5 L extracted RNA, 1X RT-buffer (Invitrogen), 0.5 mM of each dNTP (Roche), 2 ng random hexamer, 10 mM DTT (Invitrogen), 1 U of RNAse inhibitor and 2 U RT Superscript III. cDNA synthesis was performed using an Eppendorf Master thermocycler gradient system (PerkinElmer Corporation, Foster City (CA), USA) under the following conditions: 10 min at 25 • C, 60 min at 50 • C and 15 min at 75 • C.

Assay development
The assay was designed to detect all known RSV strains and to differentiate RSV subgroups A and B. Primers (MTH1 and MTH2B) were selected using Primer Express software v2.0 (Applied Biosystems Inc., Foster city (CA), USA) and a dataset of 45N gene sequences available from the NCBI database (Table 1) to amplify 84 bp for both RSV A and RSV B. RSV subgroups A and B differ by 5 nucleotides within the highly conserved region of the N gene; two LNA probes, RSVA and RSV B (Sigma Proligo, Singapore), were designed to include three LNA residues at these sites (+T, +A and +T for RSV A; +C, +G and +A for RSV B). The choice of reporters and compatible quenchers was based on a Bio-Rad guideline (Biorade, 2006). Primers and probes were analyzed for false priming sites, hybridizations, and secondary structures at http://lnatools.com/hybridization. Plasmids containing the amplicons from RSV A and RSV B clinical strains were constructed for use as positive controls and generation of standard curves for quantitation. Briefly, target cDNAs amplified via RT-PCR from RSV A and RSV B-positive clinical specimens as identified by Seeplex TM RV detection kit, were cloned into pCRII ® -TOPO ® vector and transformed into the TOPO-10 bacterial strain using the TOPA TA cloning kit (Invitrogen). Clones were screened by PCR using M13 primers, and confirmed by sequencing (CEQ DTCS Quick Start Kit, CEQ 8000 Capillary Sequencer, Beckman Coulter, Inc., Fullerton (CA), USA). Quantitation of the plasmid preparations (ng/l) was conducted using PicoGreen dsDNA Quantitation Reagent (Invitrogen) in duplicate, and the corresponding copy number of target cDNA in the plasmid was calculated (copies/mL). Serial 10-fold dilutions were made for analysis of analytical sensitivity and used as standard curves for viral load quantitation. Viral load was expressed as copies per mL of sample, A Ct-value of 45 was chosen as a cut-off value for standard positivity.

Assay protocol
The RSV LNA assay was performed on a DNA Engine Peltier Thermocycler and Chromo 4 Real-time PCR system detector (Bio-Rad, Hercules, USA), in a total volume of 25 l containing 5 ul template cDNA, 0.05 U Hotstart Taq (Invitrogen), 1× PCR Buffer, 5 mM MgCl 2 , 0.2 mM of each dNTP (Roche), 0,4 M of each primer (RSV and EAV), 0.2 M of each probe for RSV, and 0.12 M of EAV probe (Table 1). The PCR employed the following thermal cycler settings: 15 min at 95 • C, followed by 45 cycles of 30 s at 95 • C and 1 min at 45 • C. Negative controls and standards were included in every assay run. The number of RSV RNA transcripts (copies/mL) of clinical samples was determined based on the standard curve. A Ct-value of 45 was chosen as the cut-off value for sample positivity. All clinical samples required a Ct value of 30 to 35 for the EAV internal control to be considered valid.
Potential interference between primers/probes for target (RSV) and internal control (EAV) was assessed by running serial dilutions of two clinical samples with confirmed high viral loads of RSV A and RSV B (data not shown).

Assay performance assessment
The analytical limit of detection was determined using six 10-fold dilutions of control plasmids from 9.0 × 10 1 to 9.0 × 10 6 copies/mL for RSV A and 6.0 × 10 2 to 6.0 × 10 7 copies/mL for RSV B. The viral load in clinical samples was assessed by running them in parallel with quantitative standards. Since the concentration of target cDNA in plasmid controls was known, the corresponding concentration of cDNA in the clinical sample could be estimated. Samples that tested positive by RSV LNA but not by Seeplex were further confirmed by a seminested RT-PCR targeted at the second hypervariable region of the RSV G gene as described previously (Parveen et al., 2006). Specificity of the assay was assessed by presence of non-RSV viruses within the sample set (as determined by Seeplex). Specificity was further assessed by attempted amplification of culture supernatants of selected human seasonal influenza isolates [A/Sydney/5/97(H3N2), A/New Caledonia/20/99(H1N1), B/Yamanashi/166/98].

Statistical analysis
Statistical comparisons between Seeplex RT-PCR and the RSV-LNA assay were performed using McNemar's test. Assay variability was quantified by the standard deviation (SD) and coefficient of variation (CV). Ct-values (of five 10-fold dilutions of control

Assay performance
The analytical limit of detection was determined using six 10-fold serial dilutions of control plasmid in sextuplicate. 100% detection was observed at concentrations ranging from 9.0 × 10 2 to 9.0 × 10 6 copies/mL for RSV A and 6.0 × 10 3 to 6.0 × 10 7 copies/mL for RSV B. The detection rate was 33% (2/6) for 9.0 × 10 1 copies/mL of RSV A and 50% (3/6) for 6.0 × 10 2 copies/mL of RSVB. Significant negative linear relationships were observed between copies/mL and Ct values for the first five 10-fold dilutions from 9.0 × 10 2 to 9.0 × 10 6 copies/mL for RSV A (adjusted R 2 = 0.97) and 6.0 × 10 3 to 6.0 × 10 7 copies/mL for RSV B (adjusted R 2 = 0.99). To assess the analytical intra-and inter-assay reproducibility of the amplification step of the assay, Ct values were compared between replicates of five 10-fold serial dilutions tested in the same batch (Table 2) or on different days (Table 3). Coefficient of variation (CVs) for analytical intra-assay variability of RSV A and B ranged from 0.9 to 3.3 and 0.6 to 3.2, respectively. Similarly, CVs for analytical inter-assay variability of both RSV A and B ranged from 1 to 3. Specificity of the assay was confirmed by failure to amplify from 81 specimens known to contain other respiratory viruses (as detected by Seeplex assay). These included influenza viruses (n = 38), rhinovirus (n = 10), human coronavirus (n = 14), human metapneumovirus (n = 11), parainfluenza virus (n = 14), and adenovirus (n = 6). In addition, no amplification was observed for culture supernatants containing high titer influenza viruses. The additional 20 samples positive by the RSV LNA assay were confirmed using an RSV G gene RT-PCR (Parveen et al., 2006). None of the samples negative by RSV LNA assay were positive by Seeplex. Fluorescence signals of EAV internal control were detected for all tested samples at Ct values ranging from 30 to 35.
between higher age and lower log-VL was also found (correlation −0.19; p = 0.09) among RSV positive patients, and younger age was associated with more severe disease (correlation 0.27, p = 0.02). When the dataset was adjusted for age, differences between viral load and severity lost statistical significance (data not shown). No significant differences in RSV viral load were observed between patients with single infections versus co-infections, or between patients requiring oxygen versus patients who did not.

Discussion
RSV was the pathogen detected most frequently among hospitalized children with acute respiratory infection in a prospective study conducted between 2004 and 2008 in pediatric ICUs and respiratory wards of the Hospital for Tropical Diseases in Ho Chi Minh City, Viet Nam (Do et al., 2011). Rapid and sensitive detection of RSV is important for installing infection control measures and thus preventing nosocomial spread, which has been recognized as a major risk in pediatric wards (Macartney et al., 2000). Several RT-PCR based methods and real-time RT-PCR have been described for RSV detection (Borg et al., 2003;Dewhurst-Maridor et al., 2004;do Nascimento et al., 2008;Falsey et al., 2003;Gueudin et al., 2003;Hu et al., 2003;Kuypers et al., 2006Kuypers et al., , 2004Mentel et al., 2003;Perkins et al., 2005;van de Pol et al., 2006van de Pol et al., , 2010van Elden et al., 2003). An improved PCR diagnostic tool using LNA probes targeting the most conserved region of the N gene was developed, to detect and distinguish simultaneously subgroup A and B and to quantify viral RNA load. An internal control virus (EAV) was added to each specimen in order to monitor efficiency of RNA extraction, reverse transcription, and amplification. This assay exhibited an increased diagnostic yield, from 24% to 32%, compared to a commercial conventional multiplex PCR assay, the Seeplex TM RV detection kit, and produced a diagnostic result within 4 h of sample receipt.
Previous studies have reported higher sensitivities of PCR methods targeting the N gene, compared to methods targeting other genes (Perkins et al., 2005;Rohwedder et al., 1998). Primers and probes were originally designed based on available N gene sequences from NCBI GenBank, including the two reference strains A2 and B. Twenty-four more recent submissions of N gene sequences from Cambodia showed perfect matches to these primers and probes. This indicates that these primers and probes theoretically could detect all currently known circulating strains of RSV. The 100% detection rate for 9.0 × 10 2 copies/mL for RSV A and for 6.0 × 10 3 copies/mL for RSV B are in the same range as previously published methods that used dilution series of plasmid DNA and, when re-calculated and assuming a 1:1 RT efficiency, as methods that used RNA transcripts (ranging from 10 3 to 10 6 copies/mL) (Borg et al., 2003;Gueudin et al., 2003;Kuypers et al., 2004). Significant linear correlations of Ct values and copies/mL were observed at levels of viral load typically found in clinical specimens. In addition, the assay exhibited good analytical reproducibility as evidenced by intra-and inter-assay CVs ranging from 1% to 3%. CVs from other reported diagnostic approaches (based on copies/mL) ranged from 1% to 35.8% (Borg et al., 2003;Gueudin et al., 2003;Kuypers et al., 2004). As RSV pathogenesis is multifactorial and varied (Collins and Graham, 2007), the relationships between disease severity, viral load and RSV subgroups were studied. Associations between RSV viral load and disease severity have been reported previously (Buckingham et al., 2000;Fodha et al., 2007;Houben et al., 2010), and several studies have definitively shown a negative correlation between RSV viral load and age (Borg et al., 2003;Gueudin et al., 2003;Kuypers et al., 2004). The results of the current study confirm a negative correlation between RSV viral load and age, and significant difference in viral load between severe and non-severe patients (defined by admission status to the pediatric ICU versus the Respiratory ward). However, when the dataset was adjusted for age, differences between viral load and severity lost statistical significance.
Previous studies have also suggested possible differences in pathogenicity between RSV A and B (Campanini et al., 2007;Gerna et al., 2008;Kuypers et al., 2004;Perkins et al., 2005). Here we report no significant differences in viral load or severity between RSV subgroups. Variations in sample collection may also affect observed viral loads, however sample collection for this cohort was conducted by experienced and well-trained hospital personnel using standard operating procedures In conclusion, this novel real-time RSV LNA assay offers a rapid, specific, highly reproducible and more sensitive alternative (32%) to a commercial multiplex RT-PCR (24%) for RSV detection, and in addition provides quantitation and subgrouping. This information is useful for epidemiological surveillance and clinical decisionmaking.