Surveillance in eastern India (2007-2009) revealed reassortment event involving ns and PB1-F2 gene segments among co-circulating influenza a subtypes
© Sarkar et al; licensee BioMed Central Ltd. 2012
Received: 1 August 2011
Accepted: 5 January 2012
Published: 5 January 2012
Influenza A virus encodes for eleven proteins, of which HA, NA, NS1 and PB1-F2 have been implicated in viral pathogenicity and virulence. Thus, in addition to the HA and NA gene segments, monitoring diversity of NS1 and PB1-F2 is also important.
55 out of 166 circulating influenza A strains (31 H1N1 and 24 H3N2) were randomly picked during 2007-2009 and NS and PB1-F2 genes were sequenced. Phylogenetic analysis was carried out with reference to the prototype strains, concurrent vaccine strains and other reference strains isolated world wide.
Comparative analysis of both nucleotide and deduced amino acid sequences, revealed presence of NS gene with A/PR/8/34(H1N1)-like mutations (H4N, Q21R, A22V, K44R, N53D, C59R, V60A, F103S and M106I) in both RNA-binding and effector domain of NS1 protein, and G63E, the HPAI-H5N1-like mutation in NEP/NS2 of five A/H1N1 strains of 2007 and 2009. NS1 of other A/H1N1 strains clustered with concurrent A/H1N1 vaccine strains. Of 31 A/H1N1 strains, five had PB1-F2 similar to the H3N2 strains; six had non-functional PB1-F2 protein (11 amino acids) similar to the 2009 pandemic H1N1 strains and rest 20 strains had 57 amino acids PB1-F2 protein, similar to concurrent A/H1N1 vaccine strain. Interestingly, three A/H1N1 strains with H3N2-like PB1-F2 protein carried primitive PR8-like NS gene. Full gene sequencing of PB1 gene confirmed presence of H3N2-like PB1 gene in these A/H1N1 strains.
Overall the study highlights reassortment event involving gene segments other than HA and NA in the co-circulating A/H1N1 and A/H3N2 strains and their importance in complexity of influenza virus genetics. In contrast, NS and PB1-F2 genes of all A/H3N2 eastern India strains were highly conserved and homologous to the concurrent A/H3N2 vaccine strains suggesting that these gene segments of H3N2 viruses are evolutionarily more stable compared to H1N1 viruses.
KeywordsNS PB1-F2 A/H1N1 A/H3N2 Reassortment
Influenza A virus (IAV) is a cytolytic virus that is responsible for significant morbidity and mortality worldwide per year. The genome of IAV consists of eight single-stranded, negative- sense viral RNA segments encoding the subunits of the transcriptase complex (PB1, PB2, PA), nucleoprotein (NP), the matrix protein (M1), two non-structural proteins (NS1 and NS2/NEP), three integral membrane proteins (hemagglutinin (HA), neuraminidase (NA) and proton channel (M2)) and the eleventh gene product PB1-F2 which is encoded by an alternative ORF of segment 2 . Due to the segmented RNA genome, multiple subtypes, large number of hosts, IAVs cause yearly seasonal epidemics and have caused four pandemics in the last 100 years. Thus, there is an intense interest in understanding genomic diversity of virus encoded genes implicated in pathogenicity of diseases.
One such virulence factor is NS1, which is a multifunctional protein of IAV having role in suppression of host immune and apoptotic responses [2, 3]. The major role of NS1 is to antagonize the antiviral response of the host by preventing the activation of NF-κB and induction of alpha/beta interferon (IFN-α/β) . It is additionally involved in (i) inhibiting the pre-mRNA 3'- end processing by binding to two 3'- end processing factors, namely cleavage and polyadenylation specificity factor and poly(A)- binding protein II [5–7]; (ii) blocking the post-transcriptional processing and nuclear export of cellular mRNA ; (iii) stimulating the translation of matrix (M1) proteins [8, 9]; (iv) inhibiting the activation of a protein kinase that phosphorylates the eIF-2 translation initiation factor by binding to double stranded (ds) RNA [10, 11], (v) induction of the phosphatidylinositol-3-kinase (PI3K/Akt) signaling pathway in order to support viral replication . Additionally, a 15 kDA nuclear export protein (NEP, formally called NS2) translated from spliced mRNA of NS gene, mediates the export of viral ribonucleoproteins from the nucleus to the cytoplasm through nuclear export signals and is involved in independent interaction with human chromosome region maintenance protein Crm1 [13, 14], as well as in viral assembly through its interaction with the M1 protein . The second virulent factor PB1-F2 is encoded in the +1 reading frame of the PB1 gene and is translated from an AUG codon downstream of the PB1 start site, probably through a leaky ribosomal scanning . It has been shown to contribute to virulence both directly and indirectly, through modulation of responses to bacteria [17, 18].The exact mechanism(s) through which virulence is increased due to PB1-F2 expression is still not clear. Though based on overexpression studies, PB1-F2 has been shown to cause cell death in some cell types [1, 19], induce inflammation by recruitment of inflammatory cells in mice  and to bind to PB1 resulting in increased activity of the influenza virus polymerase in vitro .
Since NS1 and PB1-F2 proteins have important role in viral pathogenicity, the aim of this study was a comprehensive evaluation of the IAV gene sequences encoding NS1 and PB1-F2 (segment 8 and segment 2) to understand evolution and genetic diversity of PB1-F2 and NS1 as well as NEP/NS2 in A/H1N1 and A/H3N2 strains circulating in eastern India during 2007-2009.
Sequence analysis of the NS gene
Mutations observed in the NS genes of PR8-like A/H1N1 eastern India strains.
Corresponding amino acid mutations
Name of the strains
An alternative method of analysis of the sequence data involves comparison of the silent mutations in the gene sequence since these are not subjected to selective pressure and thus are predicted to be a more reliable marker for evolutionary analysis. Except A/PR/8/34-like eastern India strains, twenty seven silent base changes occurred during the evolution of the NS gene from A/PR/8/34 to A/2007 and/or A/2009 viruses. In essence, comparison of the silent mutations in NS gene sequence of all strains revealed similar evolutionary pattern as compared to one obtained when total nucleotide changes are used.
Sequence analysis of PB1-F2 gene
The complete nucleotide sequence of the NS gene and partial sequence of PB1 gene segment encoding full-length PB1-F2 of representative influenza A (H1N1/H3N2) positive samples collected from the out-patient departments (OPDs) of local hospitals were compared with the concurrent influenza A (H1N1/H3N2) strains, circulating worldwide. Cumulative point mutations and reassortment events due to segmented RNA genome contribute to continuous genetic and antigenic variation in circulating influenza viruses resulting in seasonal epidemics. Unlike HA and NA surface glycoproteins, mutations in the NS genes appeared to be sequential, suggesting that reassortment has probably not contributed significantly to the evolution of the NS gene of these human viruses which is in agreement with previous studies [28, 29]. In addition, due to relatively conserved nature of NS gene, reassortment events may have precluded detection. In 1978 recombinant H1N1 viruses with P1, P2, P3 and NP genes from H3N2 still carried HA, NA, M, and NS gene from parent H1N1 subtype . However, significant differences in the NS genes of the influenza A H1N1 and H3N2 subtypes during this study were identified, which allowed detection of an NS gene reassortment .
With respect to PB1-F2 protein coding region, five out of thirty-one H1N1 strains (2007-2009) with functional PB1-F2 were evolutionarily close to co-circulating A/H3N2 strains, whereas, corresponding NS gene showed H1N1 origin (Figure 1 and 5). Presumably these five 2007 H1N1 strains arose by reassortment between co-circulating H1N1 and H3N2 viruses in the region. For confirmation HA and NA genes (partial) and M1 (full length) were sequenced, which on analysis confirmed nucleotide identity with A/H1N1 strains. The PB1 gene segment of these strains, however, clustered with A/H3N2 strains suggesting that although these viruses were of H1N1 origin, they probably had derived PB1 segment from an H3N2 virus. The significance of selectively lateral transmission of PB1-F2 gene among co-circulating strains is not clear, but since PB1-F2 protein is associated with pathogenesis, it may confer improved infectivity or replication efficiency. As reported earlier by our group , six A/H1N1 strains had truncated 11 aa PB1-F2 similar to 2009 pH1N1 viruses. Rest twenty A/H1N1 strains with 57 aa PB1-F2 peptide were similar to the concurrent A/H1N1 vaccine strains (Figure 6).
Surprisingly, NS nucleotide sequences of five A/H1N1 strains was highly homologous (>97%) with the 1934 prototype strain [A/Puerto Rico/8/34(H1N1)] (Figure 2). In addition, these five strains contained G63E substitution in NEP, similar to the highly pathogenic avian influenza H5N1 viruses, which may confer higher pathogenicity . To verify possible cross contamination, BLAST search of HA, NA, M1 and NS1 gene sequences showed only NS1 having sole identity with A/PR/8/34(H1N1). Thus, a chance of cross contamination with laboratory PR8 strain was ruled out. Though the frequency of vaccination in India is very low but since the WHO approved vaccines with PR8 backbone are used, possibility of reassortment with the vaccine strain can not be ruled out. Moreover, 3 out of five PR8-like NS1 carrying A/H1N1 2007 strains, had H3N2-like PB1-F2 gene, whereas, 2/5 had non-functional PB1-F2 similar to pandemic A/H1N1 strains of 2009 (Figure 1 and 5). Thus the circulation of prototype NS gene carrying A/H1N1 strains in 2007 and 2009, with PB1-F2 gene from diverse origin underlines the complexity of influenza virus genetics and evolution. In contrast to A/H1N1 strains, all A/H3N2 (n = 24) strains analyzed in this study revealed highly conserved NS and PB1-F2 gene, with >98.5% homology to concurrent A/H3N2 strains circulating worldwide.
Thus, it can be hypothesized that NS and PB1-F2 gene segments of H3N2 viruses are evolutionarily more stable. This is in contrast to the analysis of HA and NA genes in the region, where comparative amino acid mutation rates were observed in both H1N1 and H3N2 strains . Reassortment events not involving the surface glycoproteins HA and NA largely remain undetected due to specific use of HA and NA specific antisera or sequencing primers for identification of circulating strains in most countries. This study highlights the existence of A/H1N1 and A/H3N2 viruses with viral virulence marker genes PB1 and NS from diverse origin co-circulating in the same geographical location. Therefore, analysis of gene segments other than HA and NA genes, is important to understand evolution of strains with variable pathogenic potential.
The study was approved by the Institutional Ethics Committee, National Institute of Cholera & Enteric Diseases (Indian Council of Medical Research). The concerned members of this committee are: Hon'ble Justice Pinaki Ghosh, Kolkata High Court (Chairman), Dr. D. Mahalanabis (Member), Prof. A. C. Ghosh (Member), Prof. Biswapati Mukherjee (Member), Prof. Mrinmoy Ghosh (Member), Prof. Mrinal Kanti Chatterjee (Member), Mr. M. A. Wohab (Member), Mr. Amitrajit Ukil (Member), Mrs. Debolina Sarkar (Member), Prof. Subir Kumar Dutta (Member) and Dr. Phalguni Dutta (Member Secretary). Written informed consent was obtained from all participants involved in this study.
Nasal and throat swabs were collected in Viral Transport Medium (VTM) from patients with influenza-like illness reporting in outpatient's ward of two referral hospitals; Dr. B.C. Roy Memorial Hospital for children (BCRMHC) and R.G. Kar Medical College and Hospital (RGKMCH) during 2007- 2009, as described previously . Of 166 influenza A positive samples, 55 samples were picked randomly (31 H1N1 and 24 H3N2) for sequencing of PB1-F2 and NS genes. For confirmation, total PB1 gene segment was sequenced in 19 strains chosen randomly. All sequences were submitted to Genbank and compared for nucleotide and amino acid homology.
Viral RNA Extraction
Extraction of viral RNA from the clinical samples was carried out using commercially available QiaAmp Viral RNA Mini Kit (Qiagen, GmbH, Hilden, Germany) according to the manufacturer's instruction.
Amplification of virus genes
To study the genetic diversity, full length NS, PB1, M1 and PB1-F2 encoding gene segments were amplified by RT-PCR using RevertAid™ First Strand cDNA Synthesis kit and DreamTaq™ DNA Polymerase (Fermentas Life Sciences, Burlington, Canada) as per kit protocol with specific primers for the mentioned gene segments [33, 34]. PCR products were purified by column purification using QIAquick PCR Purification Kit (Qiagen, GmbH, Hilden, Germany).
Sequencing and phylogenetic analysis
Genes were sequenced following dideoxynucleotides chain termination method of Sanger et al. (1977), in ABI Prism automated 3100 DNA sequencer (Applied Biosystem, Foster City, USA) using Big-Dye Terminator Chemistry. Sequences were compared with published cognate sequences of corresponding genes. DDBJ (DNA Data Bank of Japan) Clustal W system (version 1.83) was used for multiple sequence alignment of different nucleotide sequences of eastern India strains with other reference strains. Neighbour-Joining (N-J) trees were generated using pair-wise gap deletion, Maximum Composite Likelihood as distance measure and 1000 boot-strap replicates (generated with MEGA4) with boot strap values ≥70%.
Identification of Conserved Regions
To study the conservation among the NS1 and NS2 sequences of eastern India strains, amino acid (aa) sequences of eastern India strains were deduced by the DNA sequence translation tool EMBOSS-Transeq (EBI Group). Conserved regions were identified and mapped onto the protein structures using the web-based ConSurf server (http://consurf.tau.ac.il/) [35, 36] providing the multiple sequence alignment as input. The degree of conservation was subdivided into nine grades, with grade 1 being the least and grade 9 being the most conserved.
Accession Numbers for nucleotide sequences
Accession numbers of all newly reported sequences obtained for full length non-structural (NS) gene, full length polymerase basic protein 1 (PB1) gene, PB1-F2 gene and matrix (M1) gene in this study are: [GenBank: HQ853495 to HQ853549; HM460507, HM460510, HQ154076, HQ154078, HQ154084, HQ853550, HQ853552, HQ853555, HQ853556, HQ853557, HQ853559, HQ853567, HQ853571, HQ853573, HQ853574, HQ853579, HQ853586, HQ853589, HQ853591; HQ853551, HQ853553, HQ853554, HQ853558, HQ853560 to HQ853566, HQ853568 to HQ853570, HQ853572, HQ853575 to HQ853578, HQ853580 to HQ853585, HQ853587, HQ853588, HQ853590, HQ853592 to HQ853595, HM460508, HM460509, HM460511, HQ154075; JF903566 to JF903570].
Authors acknowledge Swati Ghosh and Malay Dey-Sarkar for technical support in sample processing and virus isolation in lab.
- Chen W, Calvo PA, Malide D, Gibbs J, Schubert U, Bacik I, Basta S, O'Neill R, Schickli J, Palase P, Henklein P, Yewdell JW: A novel influenza A virus mitochondrial protein that induces cell death. Nat Med 2001, 7:1306–1312.PubMedView Article
- Hale BG, Randall RE, Ortin J, Jackson D: The multifunctional NS1 protein of influenza A viruses. J Gen Virol 2008, 89:2359–2376.PubMedView Article
- Lin D, Lan J, Zhang Z: Structure and function of the NS1 protein of influenza A virus. Acta Biochim Biophys Sin (Shanghai) 2007, 39:155–162.View Article
- Wang X, Li M, Zheng HY, Muster T, Palese P, Beg AA, Garcia-Sastre A: Influenza A virus NS1 protein prevents activation of NF-κB and induction of alpha/beta interferon. J Virol 2000, 74:11566–11573.PubMedView Article
- Chen Z, Li Y, Krug RM: Influenza A virus NS1 protein targets poly (A)- binding protein II of the cellular 3'-end processing machinery. EMBO J 1999, 18:2273–2283.PubMedView Article
- Fortes P, Beleso A, Ortin J: Influenza virus NS1 protein inhibits pre-messenger RNA splicing and blocks messenger RNA nucleocytoplasmic transport. EMBO J 1994, 13:704–712.PubMed
- Numeroff ME, Barabino SML, Li YZ, Keller W, Krug RM: Influenza virus NS1 protein interacts with the cellular 30 kDa subunit of CPSF and inhibits 3'-end formation of cellular pre-mRNAs. Mol Cell 1998, 1:991–1000.View Article
- Enami K, Sato TA, Nakada S, Enami M: Influenza virus NS1 protein stimulates translation of the M1 protein. J Virol 1994, 68:1432–1437.PubMed
- Marión RM, Aragón T, Beleso A, Nieto A, Ortín J: The N-terminal half of the influenza virus NS1 protein is sufficient for nuclear retention of mRNA and enhancement of viral mRNA translation. Nucleic Acids Res 1997, 25:4271–4277.PubMedView Article
- Aragón T, de La luna S, Novoa I, Carrasco L, Ortín J, Nieto A: Eukaryotic translation initiation factor 4GI is a cellular target for NS1 protein, a translational activator of influenza virus. Mol Cell Biol 2000, 20:6259–6268.PubMedView Article
- Lu Y, Wambach M, Katze MG, Krug RM: Binding of the influenza virus NS1 protein to double-stranded RNA inhibits the activation of the protein kinase that phosphorylates the eIF-2 translation initiation factor. Virology 1995, 214:222–228.PubMedView Article
- Ehrhardt C, Ludwig S: A new player in a deadly game: influenza viruses and the PI3K/Akt signaling pathway. Cell Microbiol 2009, 11:863–871.PubMedView Article
- Newmann G, Hughs MT, Kawaoka Y: Influenza A virus NS2 protein mediates vRNP nuclear export through NES-independent interaction with hCRM1. EMBO J 2000, 19:6751–6758.View Article
- O'Neill RE, Talon J, Palese P: The influenza virus NEP (NS2 protein) mediates the nuclear export of viral ribonucleoproteins. EMBO J 1998, 17:288–296.PubMedView Article
- Schmitt AP, Lamb RA: Influenza virus assembly and budding at the viral budozone. Adv Virus Res 2005, 64:383–446.PubMedView Article
- McAuley JL, Zhang K, McCullers JA: The effects of influenza A virus PB1-F2 protein on polymerase activity are strain specific and do not impact pathogenesis. J Virol 2009, 84:558–564.View Article
- Conenello GM, Zamarin D, Perrone LA, Tumpey T, Palese P: A single mutation in the PB1-F2 of H5N1 (HK/97) and 1918 influenza A viruses contributes to increase virulence. PLoS Pathog 2007, 3:1414–1421.PubMedView Article
- McAuley JL, Hornung F, Boyd KL, Smith AM, McKeon R, Bennink J, Yewdell JW, McCullers JA: Expression of the 1918 influenza A virus PB1-F2 enhances the pathogenesis of viral and secondary bacterial pneumonia. Cell Host Microbe 2007, 2:240–249.PubMedView Article
- Gibbs JS, Malide D, Hornung F, Bennink JR, Yewdell JW: The influenza virus PB1 F2 protein targets the inner mitochondrial membrane via a predicted basic amphipathic helix that disrupts mitochondrial function. J Virol 2003, 77:7214–7224.PubMedView Article
- Mazur I, Anhlan D, Mitzner D, Wixler L, Schubert U, Ludwig S: The proapoptotic influenza A virus protein PB1-F2 regulates viral polymerase activity by interaction with the PB1 protein. Cell Microbiol 2008, 10:1140–1152.PubMedView Article
- Ludwig S, Schultz U, Mandler J, Fitch WM, Scholtissek C: Phylogenetic relationship of the nonstructural (NS) genes of influenza a virus. Virology 1991, 183:566–577.PubMedView Article
- Suarez DL, Perdue ML: Multiple alignment comparison of the non-structural genes of influenza A viruses. Virus Res 1998, 54:59–69.PubMedView Article
- Zohari S, Gyarmati P, Ejdersund A, Berglof U, Thoren P, Czifra G, Belak S, Waldenstrom J, Olsen B, Berg M: Phylogenetic analysis of the non-structural (NS) gene of influenza A viruses isolated from mallards in northern Europe in 2005. Virol J 2008, 5:147.PubMedView Article
- Akarsu H, Burnmeister WP, Petosa I, Petit CW, Muller R, Baudin F: Crystal structure of the M1 protein-binding domain of the influenza A virus nuclear export protein (NEP/NS2). EMBO J 2003, 22:4646–4655.PubMedView Article
- Sarkar M, Agrawal AS, Sharma Dey R, Chattopadhyay S, Mullick R, De P, Chakrabarti S, Chawla-Sarkar M: Molecular characterization and comparative analysis of pandemic H1N1/2009 strains with co-circulating seasonal H1N1/2009 strains from eastern India. Arch Virol 2011, 156:207–217.PubMedView Article
- Horisberger MA: The large P proteins of influenza a viruses are composed of one acidic and two basic polypeptides. Virology 1980, 107:302–305.PubMedView Article
- Ulmanen I, Broni BA, Krug RM: Role of two of the influenza virus core P proteins in recognizing cap 1 structures (m7GpppNm) on RNAs and in initiating viral RNA transcription. Proc Natl Acad Sci USA 1981, 78:7355–7359.PubMedView Article
- Scholtissek C, Rohde W, von Honingen V, Rott R: On the origin of the human influenza virus subtypes H2N2 and H3N2. Virology 1978, 87:13–20.PubMedView Article
- Scholtissek C, von Honingen-Huene V: Genetic relatedness of the gene which codes for the non-structural (NS) protein of different influenza A strains. Virology 1980, 102:13–20.PubMedView Article
- Young J, Palese P: Evolution of human influenza a viruses in nature: recombination contributes to genetic variation of H1N1 strains. Proc Natl Acad Sci USA 1979, 76:6547–6551.PubMedView Article
- Zohari S, Gyarmati P, Thorén P, Czifra G, Belak S, Bröjer C, Belák S, Berg M: Genetic characterization of the NS gene indicates co-circulation of two sub-lineages of highly pathogenic avian influenza virus of H5N1 subtype in northern Europe in 2006. Virus Genes 2008, 36:117–125.PubMedView Article
- Agrawal AS, Sarkar M, Ghosh S, Roy T, Chakrabarti S, Lal R, Mishra AC, Chadha MS, Chawla-Sarkar M: Genetic characterization of circulating seasonal Influenza A viruses (2005–2009) revealed introduction of oseltamivir resistant H1N1 strains during 2009 in eastern India. Infect Genet Evol 2010, 10:1188–1198.PubMedView Article
- Chan C-H, Lin K-L, Chan Y, Wang Y-L, Chi Y-T, Tu H-L, Shieh H-K, Liu W-T: Amplification of the entire genome of influenza A virus H1N1 and H3N2, subtypes by reverse-transcription polymerase chain reaction. J Virol Met 2006, 136:38–43.View Article
- Chen G-W, Yang C-C, Tsao K-C, Huang C-G, Lee L-A, Yang W-Z, Huang Y-L, Lin TY, Shih S-R: Influenza A virus PB1-F2 gene in recent Taiwanese isolates. Emerg Infect Dis 2004, 10:630–636.PubMed
- Glaser F, Pupko T, Paz I, Bell RE, Bechor-Shental D, Martz E, Ben-tal N: Consurf: identification of functional regions in proteins by surface-mapping of phylogenetic information. Bioinformatics 2003, 19:163–164.PubMedView Article
- Landau M, Mayrose I, Rosenberg Y, Glaser F, Martz E, Pupko T, Ben-Tal N: ConSurf 2005: the projection of evolutionary conservation scores of residues on protein structures. Nucleic Acids Res 2005, 33:W299-W302.PubMedView Article
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.