key: cord-292831-oihcay6w
authors: Choudhary, Manohar L.; Anand, Siddharth P.; Heydari, Mostafa; Rane, Grishma; Potdar, Varsha A.; Chadha, Mandeep S.; Mishra, Akhilesh C.
title: Development of a multiplex one step RT-PCR that detects eighteen respiratory viruses in clinical specimens and comparison with real time RT-PCR
date: 2013-01-08
journal: J Virol Methods
DOI: 10.1016/j.jviromet.2012.12.017
sha: 
doc_id: 292831
cord_uid: oihcay6w

Rapid and accurate diagnosis of viral respiratory infections is crucial for patient management. Multiplex reverse transcriptase polymerase chain reaction (mRT-PCR) is used increasingly to diagnose respiratory infections and has shown to be more sensitive than viral culture and antigen detection. Objective of the present study was to develop a one-step mRT-PCR that could detect 18 respiratory viruses in three sets. The method was compared with real time RT-PCR (rRT-PCR) for its sensitivity and specificity. Clinical specimens from 843 pediatric patients with respiratory symptoms were used in the study. 503 (59.7%) samples were detected positive by mRT-PCR. Of these 462 (54.8%) exhibited presence of a single pathogen and 41 (4.9%) had multiple pathogens. rRT-PCR detected 439 (52.1%) positive samples, where 419 (49.7%) exhibited one virus and 20 (2.4%) showed co-infections. Concordance between mRT-PCR and rRT-PCR was 91.9% and kappa correlation 0.837. Sensitivity and specificity of mRT-PCR were 99.5% and 83.7% while that of rRT-PCR was 86.9% and 99.4% respectively. Rhinovirus (17.2%) was the most frequently detected virus followed by respiratory syncytial virus B (15.4%), H1N1pdm09 (8.54%), parainfluenza virus-3 (5.8%) and metapneumovirus (5.2%). In conclusion, mRT-PCR is a rapid, cost effective, specific and highly sensitive method for detection of respiratory viruses.

Respiratory viral disease is a serious public health problem globally and has social and economic impact. Respiratory viruses cause disease that ranges from mild respiratory illness to fatal pneumonia and contribute significantly to morbidity and mortality worldwide (Williams et al., 2002) . Clinical presentations of respiratory viral infections are very similar to each other and etiological diagnosis based on clinical parameters is difficult. Rapid diagnosis is essential for prompt patient management and public health action.

Diagnosis of viral respiratory infections has been based on the use of conventional methods such as isolation by cell culture and antigen detection. Although these methods are effective and often complementary, they have certain limitations. Cell culture is considered to be the "gold standard" for virus detection, but it is too laborious and time consuming. Antigen detection is insufficiently sensitive and/or specific.

Sensitivity and specificity of respiratory virus detection have improved considerably with the advent of nucleic acid amplification tests (NATs; Bellau-Pujol et al., 2005; Fox, 2007; Freymuth et al., 2006; Liolios et al., 2001; Vabret et al., 2000; Weinberg et al., 2004) . To enable a rapid response to a potential outbreak, it is desirable to have fast, accurate and comprehensive diagnostic methods that are capable of detecting simultaneously and subtyping viruses. However, it is very expensive to use virus specific PCR assays. Several groups have developed multiplex RT-PCRs to identify respiratory viruses in clinical samples (Bellau-Pujol et al., 2005; Bharaj et al., 2009; Brittain-Long et al., 2010; Coiras et al., 2004; Gunson et al., 2005; Jansen et al., 2011; Kim et al., 2009; Lee et al., 2007; Mahony et al., 2007; Pabbaraju et al., 2008; Wang et al., 2009) . The preference of one test over the other depends on its specificity, sensitivity and turnaround time as well as cost in resource limited settings.

The purpose of this study was to develop a one step mRT-PCR that could detect respiratory viruses including influenza A viruses, H1N1pdm09, seasonal H1N1, H3N2, influenza B viruses, respiratory syncytial virus (RSV) A and B, human metapneumovirus (HMPV), parainfluenza viruses (PIV) 1, 2, 3, 4, rhinovirus, enterovirus, corona viruses OC43, 229E, NL63 and HKU1 in three sets in human clinical samples and to compare it with rRT-PCR.

Primers used in this assay were designed using dual priming oligonucleotide (DPO) technology which allows specific detection of viruses without any non specific amplification (Chun et al., 2007) . Conserved regions of target genes were chosen to design the forward and reverse primers. For mRT-PCR, different viral genes with their PCR product sizes are given in Table 1 . Primer sequences can be provided on request. The above mRT-PCR primers were obtained from Sigma (Bangalore, India). Multiplex RT PCRs were carried out in 3 tubes (sets) (set 1: influenza A, subtype seasonal H1, H3 and pandemic H1 (2009), influenza B; set 2: RSV-A, RSV-B, HMPV, PIV-1, -2, -3; and set 3: PIV-4, corona virus OC43/HKU1, corona virus 229E/NL63, rhinovirus/enterovirus). The mRT-PCR did not differentiate between corona virus OC43 and HKU1; 229E and NL63; and between rhinovirus and enterovirus. Monoplex one step real time RT PCRs (rRT-PCR) for influenza A, subtype seasonal H1, H3, pandemic H1 (2009), influenza B and internal control RnaseP were carried out as per protocol provided by Centers for Disease Control and Prevention, USA (CDC) (CDC protocol, 2009 ). For RSV-A and B, PIV-1, -2, -3, corona virus OC43, 229E, NL63 and rhinovirus, the rRT-PCR assay was performed as described (Gunson et al., 2005) , for HMPV (Maertzdorf et al., 2004) and for PIV-4 as described on the following website www3.appliedbiosystems.com/cms/groups/mcb/cms 088565.pdf. For rRT-PCR, all influenza, RSV-A, rhinovirus, PIV-1, and coronavirus NL63 TaqMan probes were labelled with FAM (fluorescein amidite); RSV-B, PIV-2, -3, -4, coronavirus OC43 and 229E were labelled with VIC. Real time PCR primers and TaqMan probes were obtained from Applied Biosystems, USA.

Tissue culture grown viruses used as positive controls for influenza B, seasonal influenza A H1N1, H3N2, H1N1pdm09, RSV-A and B, HMPV, PIV-1, -2, -3 and enterovirus. Clinical samples positive for PIV-4, rhinovirus, human corona virus OC43/HKU1 and 229E/NL63 were confirmed by sequencing and used to check primers cross reactivity. Nasal, nasopharyngeal or throat swab specimens of ILI (influenza like illness)/SARI (severe acute respiratory illness) patients collected from July 2009 to August 2011 and referred to National Institute of Virology by different hospitals in Pune for diagnosis of H1N1pdm09 were used in this study. The study was approved by the institutional ethics committee. A total of 843 clinical specimens of pediatric patients were tested by conventional mRT-PCR and monoplex rRT-PCR. One hundred and fourteen influenza positive samples (27 seasonal H3N2, 27 seasonal H1N1, 29 H1N1pdm09 and 31 influenza B) were used for validation of the multiplex assay set 1.

RNA was extracted using the MagMax automated RNA extractor as per manufacturer's instructions. Fifty microlitres of clinical sample was used for RNA extraction and finally eluted in 50 l of elution buffer.

Amplification for conventional mRT-PCR was carried out on GeneAmp PCR System 9700 using Invitrogen Superscript III one step RT-PCR kit. Master mix for mRT-PCR comprised of 25 l 2× buffer, 2 l enzyme mix, 20 M of forward and reverse primer each and 19 l of RNA template to make a 50 l reaction. Negative and positive controls were run with each experiment. Thermal cycling conditions were 50 • C for 30 min for reverse transcription, initial denaturation at 94 • C for 10 min, 45 cycles of three steps −15 s at 94 • C, 60 s at 55 • C, 45 s at 68 • C, and final extension at 68 • C for 5 min. 10 l of PCR products were mixed with gel loading dye and run on a 2% agarose gel containing ethidium bromide and visualized on a UV transilluminator.

Real time RT-PCR assay was performed on Applied Biosystems' ABI 7500 machine using Invitrogen Superscript III one step qRT-PCR kit. A typical 25 l PCR reaction comprised of 10 M of each forward and reverse primer, 5 M of TaqMan probe, 12.5 l 2× buffer, 0.5 l Superscript III enzyme and 5 l RNA template. Thermal cycling conditions for rRT-PCR were 50 • C for 30 min for reverse transcription, initial denaturation at 94 • C for 5 min, 45 cycles of 15 s at 94 • C, 45 s at 55 • C.

The specificity of the multiplex PCR assay was evaluated by cross reaction tests with known viral isolates and different panels of sequence confirmed known clinical respiratory samples as reference material. The assay set 1 was also tested on influenza A H5N1, H11N1, H9N2, and newcastle disease virus procured from National Institute of Virology's virus repository and WHO (World Health Organization) influenza panel for QA/QC.

To analyze sensitivity, known positive samples were subjected to monoplex PCR. The forward primer was tagged with T7 promoter and by gel electrophoresis in a 1% agarose gel, bands of expected size were observed. These bands were extracted from agarose gel using Qiagen gel extraction kit (Qiagen, Germany). Purified PCR product was quantified using a Nanoplex nanophotometer (Implant, Germany). In vitro transcription was carried out using RiboMax T7 in vitro transcription system from Promega as per manufacturer's instructions (Promega, USA). Further, RNA was purified using Qiagen's RNA extraction kit and estimated using the Nanoplex nanophotometer. RNA copy number was calculated and 10 fold serial dilutions of the in vitro RNA were performed. Copy number referred to the number of copies of the target gene used to ascertain the limit of detection of the assay. The concentrations tested for each target in set 1 were 10 6 to 10 copies and in sets 2 and 3 were 10 7 copies to 1 copy of template RNA per reaction.

The limit of detection (LOD) was defined as the lowest concentration at which each target could be detected consistently. PCR products were detected by gel electrophoresis. Positive controls were prepared by pooling in vitro RNA of each set.

Initial PCR was set up using monoplex primers for the relevant target gene with the same mastermix composition and thermal cycling conditions as described for the mRT-PCR. PCR product was purified using Qiagen's PCR purification kit according to the manufacturer's instructions. The purified product was used to set up cycle sequencing reaction with each sequencing primer as follows: mastermix: 5× sequencing buffer -2 l, cycle sequencing reaction mix RR100 -2 l, primer -1 l, DDW -3 l and template -2 l. Thermal cycling conditions: 95 • C for 5 min; 25 cycles of 94 • C −15 s, 50 • C for 10 s and 60 • C for 4 min. The above cycle sequencing product was purified using Qiagen's DyeEx cycle sequencing purification kit as per manufacturer's instructions. Applied Biosystems 3730xl 96 capillary sequencer was used for sequencing and results were analyzed using sequencing analysis 5.2 software. The sequences were confirmed using NCBI BLAST.

Statistical analysis was carried out using PASW Statistics 18 software. The agreement between two diagnostic tests was calculated by means of concordance and kappa statistics. The sensitivity and specificity of one test with reference to the other and vice versa were also computed.

Specificity of the mRT-PCR assay was evaluated by cross reaction tests against known respiratory virus isolates/positive samples and a WHO QA/QC panel for influenza viruses showed no cross reactivity amongst different viruses. A product of the expected size was obtained for each viral target in the presence of all the primers of the respective multiplex PCR and the specific products could be distinguished on agarose gel electrophoresis based on their size. For confirmation of the observed bands, sequencing of the amplified PCR products was carried out. For validation of set 1 of the multiplex assay, 114 respiratory clinical specimens previously positive for influenza viruses by rRT-PCR were used and results were 100% concordant. Co-infection was observed in one sample where both seasonal influenza A/H1N1 and A/H3N2 detected, this was confirmed by sequencing of the hemagglutinin (HA) and neuraminidase (NA) genes. Set 1 of the mRT-PCR detected specifically only H1N1pdm09, seasonal H1N1, H3N2 and all other influenza-A viruses were detected as a universal influenza-A product of 450 bp. In the WHO Influenza panel, 3 samples exhibited bands corresponding to influenza-A only. These samples were confirmed as influenza A H5N1 by sequencing the matrix gene.

Analytical sensitivity of the conventional multiplex assay was determined by testing 10 fold serial dilutions of the quantified RNA transcripts of each target. Set 1 detected 100 RNA copies per reaction for each target. Set 2 detected RSV-A at 10 copies per reaction; HMPV, PIV-1, -2 and -3 at 100 copies per reaction and RSV-B at 1000 RNA copies per reaction. Set 3 detected enterovirus at 1 copy per reaction; PIV4, rhinovirus and corona OC43/HKU1, 229E/NL63 at 1000 RNA copies per reaction in the mRT-PCR. The experiments were repeated three times using all three sets to obtain consistent limits of detection for each virus in the mRT-PCR Table 1 ).

Supplementary data associated with this article can be found, in the online version, at http://dx.doi.org/10.1016/ j.jviromet.2012.12.017.

Overall, concordance between both the methods was 91.9% and kappa correlation was 0.837. Sensitivity and specificity of mRT-PCR were 99.5% and 83.7% respectively considering rRT-PCR as the "gold standard". Conversely sensitivity and specificity of rRT-PCR were 86.9% and 99.4% respectively considering mRT-PCR as the "gold standard" (Table 2) . For seasonal H3N2, H1N1pdm09, influenza B, RSV-A and B and PIV-2, 100% concordance were observed using both the methods. The sensitivity and specificity for rhinovirus using mRT-PCR were 96.6 and 92.2% and rRT-PCR were 59.3% and 99.6% respectively. Thirty seven out of 59 samples positive for rhinovirus/enterovirus in mRT-PCR but negative by rRT-PCR were selected randomly and sequenced for the 5 -UTR region. Sequencing result confirmed 34 samples as rhinovirus (91.9%) and 3 as enterovirus (8.1%).

For HMPV, sensitivity and specificity were 100% and 97.6% for mRT-PCR while the sensitivity and specificity were 54.6% and 100% respectively for rRT-PCR. The 20 samples positive in the mRT-PCR, but negative by rRT-PCR was confirmed to be positive for HMPV by sequencing. Sensitivity and specificity by mRT-PCR were 100% and 99.4% for PIV-3 and those for rRT-PCR were 89.8% and 100% respectively. Among PIV-3 positive clinical samples, approximately 15% samples showed amplification bands corresponding to PIV-1 and PIV-2 which were tested further using monoplex primers. They were negative for PIV-1 and PIV-2 indicating that the bands were nonspecific. For PIV-4, the sensitivity and specificity were 50% and 100% for rRT-PCR respectively. Four out of the 5 samples negative by rRT-PCR were sequenced and found to be PIV-4. Randomly selected 8 corona virus OC43/HKU1 positive samples were sequenced and confirmed as corona virus OC43. Rhinovirus [145/843 (17.2%) ] was the respiratory virus detected most frequently followed by RSV-B in 130 samples (15.4%), H1N1pdm09 in 72/843 (8.5%), PIV-3 in 49 samples (5.8%) and HMPV in 44 samples (5.2%). Corona virus 229E was detected in 4 samples by rRT-PCR. As shown in supplementary Table 1 , rhinovirus was detected most frequently in co-infection followed by RSV. All coinfections were confirmed by performing monoplex PCRs with their respective primers.

Major respiratory viral pathogens recorded globally are RSV, influenza virus, metapneumovirus, parainfluenza virus, rhinovirus, enterovirus, adenovirus and bocavirus. A number of studies have attempted for development and evaluation of multiplex PCRs for detection of various respiratory viruses (Auburn et al., 2011; Bellau-Pujol et al., 2005; Boonsuk et al., 2008; Brittain-Long et al., 2010; Coiras et al., 2004; Jansen et al., 2011; Kim et al., 2009; Mahony et al., 2007; Suwannakarn et al., 2008; Wang et al., 2009) . This is the first study using single step mRT-PCR for the detection of 18 respiratory viruses from clinical specimens.

For the conventional mRT-PCR, conserved regions of genes were selected (Table 1) . Each set in the assay produced distinctively sized PCR products which could be visualized and easily differentiated by agarose gel electrophoresis. The only limitation of this multiplex PCR is lack of an internal control to check quality of the samples.

The mRT-PCR is sensitive as determined by testing 10 fold serial dilutions of quantified RNA transcripts of each target. Multiplex PCR set 1 detected 100 RNA copies/reaction for each target which is comparable with multiplex PCR developed by Kim et al. (2009) . An influenza typing kit available commercially, the Seeplex A/B typing kit has a sensitivity and specificity of 83.7% and 100% in detecting H1N1pdm09 when compared to the CDC real time RT-PCR (Choi et al., 2010) while sensitivity and specificity of our mRT-PCR was 100%, showing excellent concordance with CDC's real time RT-PCR.

Multiplex RT-PCR detected 503 (59.7%) samples positive for one or more of the respiratory viruses listed above and real time PCR detected 439 (52.1%) samples as positive (supplementary Table 1 ). This low detection rate by rRT-PCR was mainly due to the low sensitivity of detecting rhinovirus and metapneumovirus. Previously, detection rates between 43% and 63% have been reported (Brittain-Long et al., 2010; Broor et al., 2007; Coiras et al., 2004; Kim et al., 2009; Mahony et al., 2007; Maitreyi et al., 2000; Yeolekar et al., 2008) . Co-infections have been described previously in approximately 5-20% of infected patients in different studies (Brittain-Long et al., 2010; Coiras et al., 2004 , Kim et al., 2009 Mahony et al., 2007; Wang et al., 2010) .

One of the biggest advantages of mRT-PCR is its ability to detect co-infections which are often undetected in viral culture with immunofluorescence detection. Several cases of co-infection were detected using both the PCR methods. Multiplex RT-PCR detected 41 (4.9%) samples with co-infection including 3 samples in which 3 viruses were detected whereas rRT-PCR detected 20 (2.4%) samples with co-infection. One case of co-infection with influenza A/H1N1 and A/H3N2 viruses was also detected. Co-infection with different influenza viruses occurs naturally and plays an important role in evolution, epidemiology and pathogenicity of the virus.

Concordance between both the methods was 100% for seasonal H1N1, H3N2, H1N1pdm09, influenza B, PIV-2, RSV-A and B. Lower detection rates by rRT-PCR were due to the lower sensitivity of this system for rhinovirus and HMPV. This may be due to sequence variation in primer probe binding regions in Indian subtypes. Frequency of PIV-1, PIV-2, PIV-4 and coronavirus was low; a larger number of samples are required to better assess their clinical sensitivity and specificity.

The multiplex RT-PCR scores over the monoplex rRT-PCR in a resource limited settings. The cost of mRT-PCR for the 18 viruses discussed here is approximately US$ 27 whereas for the real time PCRs it is US$ 45. Moreover, the real time PCR machines are very expensive when compared to a conventional thermal cycler.

The one step multiplex RT-PCR assay developed in the present study provides a simple, rapid and cost effective method for the detection of 18 different respiratory viruses from human clinical specimens. The mRT-PCR is as sensitive and as specific as rRT-PCR based assay. This mRT-PCR assay can be used for respiratory virus surveillance as well as diagnosis.

Detection of nine respiratory RNA viruses using three multiplex RT PCR assays incorporating a novel RNA internal control transcript

Development of three multiplex RT PCR assays for the detection of 12 respiratory RNA viruses

Respiratory viral infections detected by multiplex PCR among pediatric patients with lower respiratory tract infections seen at an urban hospital in Delhi from

Detection of influenza virus types A and B and type A subtypes (H1, H3, and H5) by multiplex polymerase chain reaction

Prospective evaluation of a novel multiplex real-time PCR assay for detection of fifteen respiratory pathogens-duration of symptoms significantly affects detection rate

A prospective three-year cohort study of the epidemiology and virology of acute respiratory infections of children in rural India. PLoS ONE 6, e491. CDC protocol of real-time RTPCR for influenza A (H1N1)

Comparative analysis of the multiple test methods for the detection of pandemic influenza A/H1N1 2009 virus

Dual priming oligonucleotide system for the multiplex detection of respiratory viruses and SNP genotyping of CYP2C19 gene

Simultaneous detection of fourteen respiratory viruses in clinical specimens by two multiplex reverse transcription nested-PCR assays

Nucleic acid amplification tests for detection of respiratory viruses

Comparison of multiplex PCR assays and conventional techniques for the diagnostic of respiratory virus infections in children admitted to hospital with an acute respiratory illness

Real time RT PCR detection of 12 respiratory viral infections in four triplex reactions

Development and evaluation of a four-tube real time multiplex PCR assay covering fourteen respiratory viruses, and comparison to its corresponding single target counterparts

Rapid detection and identification of 12 respiratory viruses using a dual priming oligonucleotide system-based multiplex PCR assay

High-throughput, sensitive, and accurate multiplex PCR-microsphere flow cytometry system for large-scale comprehensive detection of respiratory viruses

Comparison of a multiplex reverse transcription-PCR-enzyme hybridization assay with conventional viral culture and immunofluorescence techniques for the detection of seven viral respiratory pathogens

Real-time reverse transcriptase PCR assay for detection of human metapneumoviruses from all known genetic lineages

Development of a respiratory virus panel test for detection of twenty human respiratory viruses by use of multiplex PCR and a fluid microbead-based assay

Rapid detection of respiratory viruses by centrifugation enhanced cultures from children with acute lower respiratory tract infections

Comparison of the Luminex xTAG respiratory viral panel with in-house nucleic acid amplification tests for diagnosis of respiratory virus infections

Typing (A/B) and subtyping (H1/H3/H5) of influenza-A viruses by multiplex real-time RT PCR assays

Comparison of three non-nested RT PCR for the detection of influenza-A viruses

Molecular monitoring of causative viruses in child acute respiratory infection in endemo-epidemic situations in Shanghai

Simultaneous detection of respiratory viruses in children with acute respiratory infection using two different multiplex reverse transcription-PCR assays

Superiority of reverse-transcription polymerase chain reaction to conventional viral culture in the diagnosis of acute respiratory tract infections in children

Estimates of worldwide distribution of child deaths from acute respiratory infections

Respiratory viruses in acute respiratory tract infections in Western India

The authors would like to thank Indian Council of Medical Research (ICMR) for funding the study. Technical team of influenza group and Atul Walimbe for statistical analysis is acknowledged.