key: cord-0004996-4rbrzv5o
authors: Choudhary, Manohar Lal; Anand, Siddharth P.; Sonawane, Nupoor S.; Chadha, Mandeep S.
title: Development of real-time RT-PCR for detection of human metapneumovirus and genetic analysis of circulating strains (2009-2011) in Pune, India
date: 2013-08-09
journal: Arch Virol
DOI: 10.1007/s00705-013-1812-6
sha: 1222247f98c199285e36df6eb36cd819454b7b37
doc_id: 4996
cord_uid: 4rbrzv5o

Human metapneumovirus (HMPV) is an important respiratory virus implicated in respiratory infections. The purpose of this study was to develop a one-step real-time RT-PCR assay that can detect all four lineages of HMPV and to identify the HMPV lineages circulating in Pune, India. Conserved regions of the nucleoprotein gene were used to design real-time primers and a probe. A total of 224 clinical samples that were positive for different respiratory viruses (including 51 samples that were positive for HMPV) were tested using the real time RT-PCR assay, and the specificity of the assay was observed to be 100 %. Using in vitro-synthesized RNA, the sensitivity of the assay was ascertained to be 100 copies of the target gene per reaction. Phylogenetic analysis of the nucleoprotein (N) and attachment glycoprotein (G) genes confirmed that this assay detected all lineages of HMPV. A2, B1 and B2 strains were observed during the study period. Our assay is highly sensitive and specific for all known lineages of HMPV, making it a valuable tool for rapid detection of the virus. A2 and B2 were the predominant subtypes circulating in Pune, Western India.

Human metapneumovirus (HMPV), first isolated from children with acute lower-respiratory-tract infections (ALRTI) in the Netherlands in 2001, is an enveloped, non-segmented RNA virus that belongs to the family Paramyxoviridae and the genus Metapneumovirus [1] . Subsequently, it has been reported globally [2] [3] [4] [5] [6] [7] [8] . It frequently causes both upper and lower acute respiratory tract infections (ARTIs) in people of all age groups [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] and in immunocompromised patients [20, 21] . Several groups have described reverse transcription polymerase chain reaction (RT-PCR)-based assays to detect HMPV targeting the nucleoprotein, fusion or polymerase genes [22] [23] [24] [25] [26] [27] [28] [29] .

MPV has been classified into two broad groups: group A, with subgroups A1 and A2, and group B, with subgroups B1 and B2 [30] [31] [32] . More recently, the A2 group has been split into two clades (A2a and A2b) [33] . However, there is limited information on the prevalence and genetic diversity of human metapneumovirus (HMPV) strains circulating in India [34] [35] [36] .

It is important to establish diagnostic HMPV assays which detect all subtypes of the virus. This study describes the development of real time RT-PCR assay for the detection of all four HMPV lineages (A1, A2, B1, B2) in respiratory specimens and genotyping of circulating strains from July 2009 to August 2011 in Pune, western India.

Clinical samples (nasal/throat swab) referred to the National Institute of Virology (NIV), Pune, for diagnosis of H1N1pdm09 from July 2009 to August 2011 were retrospectively used in this study with the approval of NIV's Human Ethics Committee [110(01)/EC-I/1056] as per Indian Council of Medical Research (ICMR) guidelines. This included outpatient and admitted cases in Pune, Maharashtra, India. A total of 224 clinical specimens were tested to measure the sensitivity and specificity of the assay.

RNA was extracted using a MagMax-96 viral RNA isolation kit according to the manufacturer's instructions. Fifty microlitres of clinical specimen was used as starting material, and the RNA was finally eluted in 50 ll of elution buffer.

Primer and probe design Primers and probe were designed using the HMPV sequences available in the NCBI database. Representative sequences of nucleoprotein genes of prototype sequences of A1, A2, B1 and B2 subtypes [31, 36] were aligned, and primers were designed from conserved regions of the consensus sequence to ensure that all lineages would be amplified and detected using this assay. Primers were checked for cross-reactivity, self-dimerization, hairpin formation and secondary structure formation using NCBI BLAST, OligoCalc software (http:// www.basic.northwestern.edu/biotools/oligocalc.html) and OligoAnalyzer. Primer and probe sequences are given in Table 1 . Sequencing primers were designed for the nucleoprotein (N) and attachment glycoprotein (G) genes as shown in Table 1 .

The entire nucleoprotein gene was amplified using T7HMPVNFP1 and HMPVNRP3 primers ( Table 1 ). The forward primer was tagged with a universal T7 promoter sequence at its 5' end. PCR amplification was performed on a GeneAmp PCR System 9700 using an Invitrogen Superscript III one step RT-PCR kit. The master mix for RT-PCR consisted of 25 ll 2x buffer, 2 ll enzyme mix, 20 lM forward primer, 20 lM reverse primer, and 5 ll of RNA template to make a 50-ll reaction. 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, 15 s at 52°C, 60 s at 68°C -and final extension at 68°C for 5 min. In vitro-synthesized RNA was made using a T7 riboprobe kit (Promega, Leiden, The Netherlands) as per manufacturer's instructions. In vitro-synthesized RNA was quantified using a Nanophotometer (IMPLEN, Germany), and the RNA copy number was calculated. The in vitro-transcribed RNA was then serially diluted tenfold and used to determine the analytical sensitivity of the real-time RT-PCR.

Real-time RT-PCR assay was performed on an Applied Biosystems ABI 7500 machine using an Invitrogen Superscript III one-step qRT-PCR kit (Invitrogen, California, USA). Different parameters of the assay, such as annealing temperature and time and primer and probe concentrations, were optimized. A standard curve was plotted using tenfold serially diluted in vitro-synthesized RNA, testing each concentration in duplicate. A typical 25-ll PCR reaction consisted of 10 lM forward and 10 lM reverse primer, 2.5 lM of TaqMan probe, 12.5 ll 29 buffer, 1 ll Superscript III enzyme, and 5 ll RNA template. Real-time RT-PCR thermal cycling conditions were as follows: 50°C for 30 min for reverse transcription, initial denaturation at 94°C for 5 min and 45 cycles of 15 s at 94°C and 45 s at 55°C.

Sequencing primers were designed for the nucleoprotein (N) and attachment glycoprotein (G) genes. The glycoprotein sequence is highly variable, and hence, subgroupspecific primers were designed as shown in Table 1 . The nucleoprotein gene was amplified as described earlier, and the PCR product was purified using a QIAquick PCR Purification Kit according to the manufacturer's instructions. Sublineages were determined using BLAST analysis of N gene sequences. After determining the subgroup, the glycoprotein gene was amplified using a subgroup-specific forward primer (HMPVGAFP1 for A1 and A2) (HMPVGBFP1 for B1 and B2) and a reverse primer HMPVLRP (Table 1) . PCR reaction setup and conditions for the G gene were the same as those used for the N genes.

Internal primers ( 

Studies for assessing the detection limit of the assay were performed using in vitro-transcribed RNA from template of lineage A and B viruses. The standard curve was plotted using tenfold serial dilutions of the in vitro-synthesized RNA in duplicate from 5 9 10 6 copies per reaction to 0.5 copies per reaction to check the PCR efficiency of the assay, and it was found to be approximately 95 % (Fig. 1) . The slope of the standard curve was found to be -3.434 (acceptable value = -3.3 to -3.8). The analytical sensitivity of the assay was 100 copies of in vitro-synthesized RNA per reaction for both lineages. The sensitivity and specificity were determined by testing HMPV-positive samples (n=51), known positive samples of respiratory viruses other than HMPV (n=128), and samples that were negative for any of the respiratory viruses (n=45). Clinical specimens that were confirmed previously as positive by multiplex RT-PCR [37] were used in this study. A total of 224 samples were tested, as shown in Table 2 . No cross-reactivity was observed with other respiratory viruses, and all samples that were positive for HMPV by multiplex RT-PCR were positive in the realtime RT-PCR assay. One known RSV-positive sample also tested positive for HMPV in the assay. Both RSV and HMPV were found to be present in the sample, as shown by sequencing, confirming a dual infection.

Phylogenetic analysis of the 39 nucleoprotein gene sequences showed homology with HMPV reference sequences, as shown in Fig. 2 For the nucleoprotein gene, 96 % sequence homology was observed between the samples and reference sequence (CAN97-83) at the nucleotide level, and 97-99 % at the amino acid level, for the A2 genotype. Also, the sample sequences showed 100 % amino acid sequence identity to Indian strain IND/06-12/313OP. For the B1 lineage, sequence identities between the sample sequences and reference strains (NL/1/99) were between 92 and 98 % at the nucleotide level and 93-99 % at the amino acid level. The amino acid sequence identity between the sample sequences and Indian strain IND/06-9/294OP was between 91 and 100 %. For the B2 lineage, sequence identity between the sample sequences and reference strains (CAN98-75) was 98 % at the nucleotide level and between 96 and 100 % at the amino acid level. As expected, the nucleoprotein sequences were highly conserved, both at the nucleotide and amino acid level, showing 95-99 % identity between the sample sequences and reference sequences.

Twenty-seven randomly selected HMPV-positive samples were sequenced for the glycoprotein gene. Of the 27 sequences, 11 sequences aligned with the A2 subtype, 3 with the B1 subtype, and 13 with the B2 subtype (Fig. 3) . Of the 27 sequences, 13 were from 2009, 13 were from 2011, and one was from 2010. From the two phylogenetic trees it was observed that the B1 subtype was observed only in 2011. However in 2011, the A2 and B2 subtypes were also found in co-circulation with the B1 subtype, as shown in Fig. 2 .

Glycoprotein gene sequence alignments of the A2 subtype showed 88 % sequence identity between sample and reference sequences (CAN97-83, NL/17/00, KOL/1038/07) at the nucleotide level and 82 % with CAN97-83 and KOL/ 1038/07, and 80 % with Ind/06-15 at the amino acid level. The B1 subtype sequences for the G gene showed 92 % sequence identity with NL1/99 at the nucleotide level and 89 % at the amino acid level, but only 87 % at the nucleotide level and 82.6 % at the amino acid level with a strain from eastern India (KOL/925/07). The sequence homologies are similar to those reported by Agrawal et al. [34] . The B2 subtype G gene sequences showed 93 % identity to the reference sequence HMPV73-1998 at the nucleotide level and 86 % to CAN98-75 at the amino acid level. Again, homology to an Indian strain (Ind/07-27) was lower, showing 77 % sequence identity at the amino acid level. The deduced G proteins of different HMPV strains had different lengths. Changes in the stop codon were observed among strains of different lineages. Those from subgroup A2a were 219 aa in length with the stop codon UAG (nt 658), A2b had 217 aa and utilized stop codon UAA at nt position 652, whereas those from subgroup B1 were 231 aa in length and terminated at nt position 694 using the stop codon UGA. Strains from the B2 subgroup had two different stop codons, resulting in a G protein with 208 aa using a UAG stop codon at nt position 625 and one with 238 aa using a UAA stop codon at nt position 715.

The attachment glycoprotein of the B2 subtype has 41 potential O-glycosylation sites on average. The region spanning amino acids 77 to 208 (extracellular domain) was checked for potential glycosylation sites. G-values of C 0.5 were considered to indicate potential O-glycosylation sites. The extracellular domain (amino acids 66-208) of the A2 subtype showed an average of 43 O-glycosylation sites. Also, A2 subtype sequences exhibited a higher ratio of serine to threonine than those belonging to the B2 subtypes.

There are four distinct genetic lineages of HMPV, A1, A2, B1 and B2, which exhibit substantial diversity. Kuypers et al. [26] developed an assay in which detection of all Fig. 3 Phylogenetic analysis of the attachment glycoprotein (G) gene. A phylogenetic tree was constructed with MEGA 5.1 software using the maximum-likelihood method. Bootstrap probabilities greater than 85 % are shown at the branch nodes. Reference sequences for each genotype were obtained from GenBank b known subtypes of HMPV was achieved by using two primer/probe sets for the fusion gene in a one-step realtime RT-PCR reaction. Each set was designed to amplify isolates belonging to subtype A or B but to be used together in one reaction mix. Pabbaraju et al. [27] developed a twostep real-time RT-PCR using a fusion gene with a lineagespecific primer-probe and a universal probe for detection of all HMPV lineages. Maertzdorf et al. [24] designed two sets of primers and probes within the nucleoprotein gene to detect all known genetic lineages of HMPV in a one-step real-time RT-PCR. The original NL-N assay developed by Maertzdorf et al. [24] performed poorly with subgroup B viruses, and it was further modified by Klemenc et al. [29] to improve its sensitivity for detecting all four lineages of HMPV. We have designed a primer-probe set using conserved regions of the nucleoprotein gene to detect all known genetic lineages of HMPV in a one-step real-time RT-PCR. Our assay detected 100 copies of in vitro-synthesized RNA per reaction for both lineages. A panel of 51 clinical samples that had previously been found to be positive for HMPV by multiplex PCR [37] was confirmed to be positive by our real-time RT-PCR assay. This assay was also concordant with an additional panel of samples that were negative for HMPV (Table 2 ). These results suggest that our one-step real-time RT-PCR assay is sensitive for detecting all genetic lineages and specific, as it did not exhibit cross-reactivity with other respiratory viruses.

Phylogenetic analysis of the nucleoprotein and glycoprotein genes in the present study showed that both group A and groups B HMPV were co-circulating in Pune, Western India, during the study period. Phylogenetic analysis reveals that in both 2009 and 2011, A2, B1 and B2 strains were seen in circulation, with predominance of the A2 and B2 subtypes. During the same period, the A2 lineage was predominantly detected in Southwest China [8] , Okinawa, Japan [38] and New York, USA [39] . A2 and B2 strains were found to be co-circulating in our study population in 2009 and 2011, but the B1 subtype was seen only in 2011 in our study. Of all the samples sequenced, only one sequence belonging to A1 subtype was seen. Agrawal et al. [34] reported a higher prevalence of subgroup A2 (77 %) than subgroup B1 ( [40] . In Central and South America, subtypes A2 and B2 were the predominant strains circulating in 2008-2009 [41] . In Italy, B1 was the predominant subtype, followed by A2 and B2 in 2008-2009, and in 2009-2010, lineage A1 was the predominantly circulating lineage, followed by A2 [42] . The HMPV lineage A2 has been reported to cause HMPV epidemics globally [43] . This is the first study from India where the complete N gene from all circulating subtypes has been analysed. Sequence analysis of the G gene revealed homology to reference strains similar to that described previously [34] . The predicted G proteins of different HMPV strains had different lengths, which may be attributed to amino acid substitutions and insertions, deletions, and/or changes in the stop codon. Five different lengths of G proteins were observed in this study, with 217 or 219 amino acids in subgroup A2, 208 or 238 amino acids in the B2 subgroup, and 231 in B1. This was the first report showing an HMPV strain with 208 amino acids circulating in Pune, India. We need to further monitor HMPV cases from Pune and other parts of India to gain insights into the circulating patterns of this particular strain of HMPV. This will help us to better understand the mechanism of emergence of new strains of the virus.

In conclusion, we have developed a one-step real-time RT-PCR assay for the detection of HMPV from all four genetic lineages. Our assay is sensitive and highly specific and can be used for future implementation in a diagnostic setting. Phylogenetic analysis shows that the A2 and B2 lineages are predominantly circulating in Pune.

A newly discovered human pneumovirus isolated from young children with respiratory tract disease

Children with respiratory disease associated with metapneumovirus in Hong Kong

Human metapneumovirus infection in Japanese children

First detection of human metapneumovirus in children with acute respiratory infection in India: a preliminary report

Human metapneumovirus genetic variability

Genetic diversity of human Real-time RT-PCR for human metapneumovirus 223

metapneumovirus over 4 consecutive years in Australia

Genetic variability and circulation pattern of human metapneumovirus isolated in Italy over five epidemic seasons

Detection and genetic diversity of human metapneumovirus in hospitalized children with acute respiratory infections in Southwest China

Virological features and clinical manifestations associated with human metapneumovirus: a new paramyxovirus responsible for acute respiratory-tract infections in all age groups

Ruuskanen O (2002) Metapneumovirus and acute wheezing in children

Human metapneumovirus infections in young and elderly adults

Human metapneumovirus in severe respiratory syncytial virus bronchiolitis

Human metapneumovirus among children hospitalized for acute respiratory illness

Seasonality and clinical features of human metapneumovirus infection in children in Northern Alberta

Epidemiology of human metapneumovirus

The role of human metapneumovirus in upper respiratory tract infections in children: a 20-year experience

Study of human metapneumovirus-associated lower respiratory tract infections in Egyptian adults

An outbreak of severe respiratory tract infection due to human metapneumovirus in a long-term care facility for the elderly in Oregon

Rates of hospitalizations for respiratory syncytial virus, human metapneumovirus, and influenza virus in older adults

Severe metapneumovirus infections among immunocompetent and immunocompromised patients admitted to hospital with respiratory infection

Humane metapneumovirus (HMPV) associated pulmonary infections in immunocompromised adults-initial CT findings, disease course and comparison to respiratory-syncytialvirus (RSV) induced pulmonary infections

Comparative evaluation of realtime PCR assays for detection of the human metapneumovirus

Molecular assays for detection of human metapneumovirus

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

Detection of human metapneumovirus RNA sequences in nasopharyngeal aspirates of young French children with acute bronchiolitis by real-time reverse transcriptase PCR and phylogenetic analysis

Detection and quantification of human metapneumovirus in pediatric specimens by real-time RT-PCR

Diagnosis and epidemiological studies of human metapneumovirus using real-time PCR

Detection and characterisation of human metapneumovirus from children with acute respiratory symptoms in northwest England

Real-time reverse transcriptase PCR assay for improved detection of human metapneumovirus

Analysis of the genomic sequence of a human metapneumovirus

Antigenic and genetic variability of human metapneumoviruses

Global genetic diversity of human metapneumovirus fusion gene

Novel human metapneumovirus sublineage

Genetic variability of attachment (G) and Fusion (F) protein genes of human metapneumovirus strains circulating during

Detection and genetic diversity of human metapneumovirus in hospitalized children with acute respiratory infections in India

Genomic analysis of four human metap-neumovirus prototypes

Development of a multiplex one step RT-PCR that detects eighteen respiratory viruses in clinical specimens and comparison with real time RT-PCR

Molecular epidemiology of human metapneumovirus from

Phylogenetic analysis of human metapneumovirus from New York State patients during February through

Genotype variability and clinical features of human metapneumovirus isolated from Korean children

Human metapneumovirus strains circulating in Latin America

Human metapneumovirus-associated hospital admissions over five consecutive epidemic seasons: evidence for alternating circulation of different genotypes

Evolutionary Dynamics Analysis of Human Metapneumovirus Subtype A2: Genetic Evidence for Its Dominant Epidemic

Real-time RT-PCR for human metapneumovirus 225

Acknowledgments The authors would like to thank the Director for support, and the Indian Council of Medical Research (ICMR) for funding the study. We would like to acknowledge Mr. Santosh Kumar Jadhav from the bioinformatics group and the technical team of the influenza group for their help.