Genetic mutations in influenza H3N2 viruses from a 2012 epidemic in Southern China
© Zhong et al.; licensee BioMed Central Ltd. 2013
Received: 2 August 2013
Accepted: 15 November 2013
Published: 26 November 2013
An influenza H3N2 epidemic occurred throughout Southern China in 2012.
We analyzed the hemagglutinin (HA) and neuraminidase (NA) genes of influenza H3N2 strains isolated between 2011–2012 from Guangdong. Mutation sites, evolutionary selection, antigenic sites, and N-glycosylation within these strains were analyzed.
The 2011–2012 Guangdong strains contained the HA-A214S, HA-V239I, HA-N328S, NA-L81P, and NA-D93G mutations, similar to those seen in the A/ Perth/16/2009 influenza strain. The HA-NSS061–063 and NNS160–162 glycosylation sites were prevalent among the 2011–2012 Guangdong strains but the NA-NRS402–404 site was deleted. Antigenically, there was a four-fold difference between A/Perth/16/2009 -like strains and the 2011–2012 Guangdong strains.
Antigenic drift of the H3N2 subtype contributed to the occurrence of the Southern China influenza epidemic of 2012.
Influenza virus undergoes rapid evolution by both antigenic shift and antigenic drift and this presents a significant challenge for vaccine design to best match with viruses likely to circulate in the coming influenza season [1, 2]. Since their emergence in 1968, influenza H3N2 viruses tend to be highly prevalent most years. The two dominant proteins in this particular virus strains are the H3 hemagglutinin (HA) and N2 neuraminidase (NA) surface glycoproteins [1, 3]. Antigenic drift, also known as the result of positive selection, generally involves both of these proteins, resulting in evasion of the host immune response . Many amino acid variations occur at the antibody binding sites, creating diversity in the HA and NA proteins and allowing them to evade host antibodies. However, these variations do not significantly alter the stereochemical or functional properties of HA and NA.
In Canada, antigenic drift was observed in 2008 when the A/Brisbane/10/2007 strain (the vaccine strain used in 2008–2009 and 2009–2010) mutated into the A/Perth/16/2009 strain (vaccine strain used in 2010–2011) . Antigenic and molecular characterization of H3N2 viruses over those three seasons revealed that the number of HA mutations is important, along with the nature and location of key mutations, and likely plays a significant role in antigenic drift. In previous work, amino acid substitutions were seen in five antigenic regions of HA1 from 2007–2011 Guangdong (for short GD) isolates; in particular, regions B (N160K) and D (K174R/N). The K189E/N/Q and T228A mutations within the receptor-binding site (RBS) were present in the 2010 strains, affecting the antigenicity of HA1 . The antigenicity of epidemic H3N2 isolates in 2011 differed from that of the A/Perth/16/2009 strain.
The evolution and epidemiology of influenza H3N2 viruses is partially related to global migration dynamics . The evolution of H3N2 influenza over the past 10 years reflects the dynamics of a global metapopulation. According to the Center for Public Health Surveillance and Information Service of China , the number of influenza cases was 2.35-fold greater for January–June 2012 (74,151 cases) than for January–June 2011 (31,551 cases). Approximately 87% (365/418) of isolates were of the H3N2 subtype; these isolates were sourced from local epidemics and sporadic cases reported in Guangdong Province from March–June 2012. We analyzed HA and NA gene sequences of H3N2 viruses isolated from July 2011 to July 2012.
Materials and methods
Viruses and genes
Influenza A H3N2 viruses used in this study
Collection date (month/day/year)
Geographical location of isolate
Cell line isolated in
GenBank sequence accession no.
Molecular detection of viral genes
Primers specific for the HA and NA genes of human H3N2 isolated from 1968–2010 were designed using Primer Premier 5.0 (Premier Biosoft International, Palo Alto, CA, USA) and synthesized by Life Technologies (Shanghai, China). These primers have been used in a previous study . Viral RNA was extracted using a Qiagen QIAamp Viral RNA mini Kit, and reverse transcription polymerase chain reactions (RT-PCRs) conducted with Qiagen Sensiscript Reverse Transcriptase and Takara PyroBest Taq. Amplicons were purified with a Qiagen Gel Extraction Kit and sequenced with an ABI PRISM BigDye Terminator v3.0 Ready Reaction Cycle Sequence Kit on an ABI PRISM 3100 Genetic Analyzer. The obtained sequences were analyzed using Lasergene 7.1.
Phylogenetic analysis and evolutionary dynamics
Virus population dynamics over time were estimated using the BEAST package v1.70  and applying Bayesian Markov chain Monte Carlo (MCMC) methods. Phylogenetic trees were generated, and the reliability of trees tested by bootstrap analysis using 1000 replicates. The MCMC analysis was run for 50,000,000 generations, with stationarity and mixing efficiency examined using Tracer.
Evolutionary selection and N-glycosylation
Single likelihood ancestor counting (SLAC) was appropriate for large alignments but possibly underestimates the number of positively selected sites . The fixed-effects likelihood (FEL) and the internal fixed-effects likelihood (IFEL) took the synonymous and non-synonymous substitutions into account and could be efficiently parallelized. Potential N-glycosylation sites were predicted using NetNGlyc 1.0 server . This software application predicts N-glycosylation sites in human proteins using artificial neural networks that examine the N-X-S/T (Asn-Xaa-Ser/Thr) amino acid sequence .
Antigenic analysis, epitope region and three-dimensional (3D) structure
Cross-reactivity of isolated viruses was investigated using hemagglutination inhibition (HI) assays . The variation sites in H3N2 virus HA and NA gene sequences isolated from 2011 and 2012 were compared with previously identified epitopes of HA and NA proteins . The 3D structures of HA and NA proteins were established using the I-TASSER server . Models were modified with UCSF Chimera1.5.3 .
Phylogenic diversity and variation
Using the five vaccine strains as a reference, amino acid substitution sites in the NA genes of Guangdong and other strains isolated from 2007–2012 were apparent (Figure 2b). Since 2009, variation sites in strains from around the world have been seen: L81P for 2011 Guangdong isolates, and Russian and Georgian isolates from 2012; D93G for A/GD/222/2011, and Russian and Guangdong isolates from 2012; S367N and K369T for Guangdong isolates from 2010–2012, Danish isolates from 2011, and Russian and Georgian isolates from 2012; N402D for 2011–2012 Guangdong isolates, Japanese and Swedish isolates from 2010, and Russian and Georgian isolates from 2012. In addition, the K249E substitution occurred in four 2007 Guangdong strains, G401D was seen in three 2008 Guangdong strains, and D127N was found in A/GD/472/2009 and three other strains.
Evolutionary selection and N-glycosylation
Only codon 173 (L/S) of the 566 codons in the HA gene was a positively selected site with a dN-dS of 3.08, and a P-value of 0.0944 for SLAC. Codons 237 (P/S/T/L) and 239 (V/I) were positively selected sites with dN-dS of 7.56 × 108 and 6.96 × 108, respectively, and P-values of 0.0318 and 0.0794, respectively for FEL. Codons 174 (K/N/R/G) and 209 (S/F) were positively selected sites with normalized dN-dS of 16.7 and 18.6 respectively, and P-values of 0.0502 and 0.0315 respectively for IFEL.
Codon 338 (L/R/F), of the 469 codons for NA was a positively selected site with a dN-dS of 9.39 × 108 and a P-value of 0.0498 for FEL. Codons 338 and 372 (S/L/F) were positively selected sites with dN-dS of 1.15 × 109 and 1.23 × 109 respectively and P-values of 0.0828 and 0.0887 respectively for IFEL.
We identified 13 potential glycosylation sites in the HA gene of the A/GD/1154/2012 isolate: NST024–026, NGT038–040, NAT054–056, NSS061–063, NCT079–081, NES138–140, NWT142–144, NGT149–151, NNS160–162, NVT181–183, NST262–264, NGS301–303 and NGT499–501. For the A/Perth/16/2009 HA gene, 11 potential glycosylation sites were present, with only NSS061–063 and NNS160–162 missing. Using the A/Perth/16/2009 HA gene as a reference, S61N was seen in 2011–2012 Guangdong and other strains, resulting in a gain of the NSS061–063 glycosylation site. The N161S substitution occurred in the 2011–2012 Guangdong strains, resulting in an increase of the NNS160–162 glycosylation sites. The S140C substitution occurred in A/Tbilisi/GNCDC0557/2012, resulting in deletion of the NES138–140 glycosylation site.
There were eight potential glycosylation sites in the NA gene of the A/GD/1154/2012 isolate: NIT61–63, NTT70–72, NWS86–88, NDT146–148, NAT200–202, NGT234–236, NDS329–331, and NRS402–404. The N402D substitution occurred in the 2011–2012 Guangdong isolates, the 2010 Japanese and Swedish isolates, and the 2012 Russian and Georgian isolates. This resulted in deletion of the NRS402–404 glycosylation site.
Antigenic analysis, epitope and 3-D structure
Hemagglutination inhibition results for the H3N2 strains examined in this study
Reciprocal HI titer for antiserum
Sites of variation associated with antigenic epitopess
D → N
129460, 130279, 128706
Y → H
128649, 128428, 129924
I → V
129617, 97393, 129652, 97243, 129656
E → A/T
A → S
V → I
129617, 97393, 129652
N → S
129763, 129381. 129136
L → P
D → G
Mutated amino acids around human leukocyte antigen (HLA)-associated sites, especially those that are typically conserved, suggest that cooperative interactions act to preserve the local structural stability and protein function when mutations occur that confer evasion of cytotoxic T lymphocytes (CTLs) . For the HA genes of isolated strains, the D69N, Y110H, I246V and E296A/T substitutions occurred around 2010 , while A214S, V239I and N328S have been present since 2011. Accumulation of HA (D69N, Y110H, I246V, E296A/T, A214S, V239I and N328S), along with NA (L81P and D93G) mutations, drove antigenic drift possibly giving rise to the H3N2 influenza epidemic in Guangdong and Southern China in 2012. The HA-A214S, -V239I and -N328S mutations occurred in the second half of 2011 (Table 1 and Figure 2). As the 2012 influenza season approached, these variations in circulating strains became prevalent, resulting in local epidemics throughout Southern China.
According to a previous study that focused on isolates from 2011 , HA1 genes were sub-divided among the P (A/Perth/16/2009-Clade) and V (A/Victoria/208/2009-Clade) clades. The former included subgroups1 and 2, while the latter included subgroups 3–6. When we compared our findings with the results from Klimov  and Huang , mutations A/N144T/D and K/N145N/S were only observed in the previous studies; the mutations we have reported here were not evident in Klimov’s work. The mutations in the 2012 Guangdong isolates appear to have given rise to antigenic drift. Generally, five HA1 antigenic regions are analyzed and associated with antigenic drift that results in an epidemic. However, for the particular epidemic we investigated, these five regions were not seemingly associated with antigenic drift. Only the A214S substitution occurred in epitope D and the RBS. Beside the five HA1 antigenic regions, many amino acid sites were unclassified but still considered important . Using the H3N2 HA and NA protein sequences of A/New York/348/ 2005 as a reference, all mutated amino acid sites were related to previously identified epitopes, which might involve in antigenic presentation / recognition / response. This indicated that mutations of B cell epitope / T cell epitope regions in the Immune Database influenced antigenicity to some degree . These B cell epitopes / T cell epitopes span 17–18 amino acids in the A/New York/348/2005 strain. Further research is required to confirm the relationship between antigenicity and mutations in HA and NA genes.
As was observed for of major capsid L1 protein of HPV-16 and −18, amino acids at positions 174 (K/N/R/G), 209 (S/F), 237 (P/S/T/L) and 239 (V/I) in HA, and at positions 338 (L/R/F) and 372 (S/L/F) in NA were positively selected sites according to SLAC, FEL and IFEL (P < 0.1), indicating that these sites avoided immunological pressure for their continued persistence .
Only 76% of positive scored sequons according to NetNGlyc are modified by N-Glycans with a bias towards Thr-containing sequons . Compared with A/Perth/16/2009, two glycosylation sites (NSS061–063 and NNS160–162) in HA genes were prevalent in 2011–2012 Guangdong strains, with the NA glycosylation site N402D deleted. There are five possible functions for glycosite migration in human influenza viruses : to more effectively mask the antigenic sites; to more effectively protect enzymatic cleavage sites of NA; to stabilize polymeric structures; to regulate receptor binding and catalytic activities; and to balance the binding activity of HA with the release activity of NA. Five HA glycosites and two NA glycosites were positive sites in this study. Gains in N-glycosylation sites were likely to be positively selected for shielding antigenic sites from immune responses . The acquisition of glycosylation at residues 144 of pH1N1 was associated with viral replication, virulence and transmissibility and provided insights into the evolution dynamics of influenza viruses with implications in vaccine immunogenicity . The loss of a glycosylation site might also be important in antibody recognition , where a deletion of the NA-NRS402–404 occurred in the 2011–2012 Guangdong and other isolates in this study. With respect to HA/NA N-glycosylation changing, further analysis is required to determine any relationship between their structure mutation and function influence. At the antigenic sites of HA protein, positive selection appeared to have effected radical and conservation substitution in term of the charge of the amino acids, suggesting that antigenic drift is not a byproduct of HA evolution in H3N2 viruses .
Kenyan H3N2 viruses isolated during 2006–2007 revealed unique genetic variations, with several amino acid substitutions located at immunodominant epitopes of the HA1 protein . These mutations included V112I at site E, K173 E at site D and N278K at site C. These mutations possibly result in a conformational change to the HA molecule, thereby exposing novel epitopes and thus abrogating the binding of pre-existing antibodies at these sites. A Canadian study in 2011 involving antigenic and molecular characterization of H3N2 viruses over three seasons revealed that the number of HA mutations was important, along with the nature and location of key mutations, and played a significant role in antigenic drift . From our findings, presented in this report, the 2012 strains had evolved genetically and antigenically from the A/Perth/16/2009 vaccine-like strains. The A/GD/1154/2012 strain was antigenically distinct from the A/GD/1105/2009 strain, suggesting that it may be the parental strain responsible for the 2012 H3N2 influenza virus epidemics in Southern China. Although influenza H3N2 viruses varied genetically and antigenically during 2009–2011 , we concluded that antigenic drift in 2012 along with the accumulation of evolutionary mutations in viruses isolated in 2012 played an important role in the influenza epidemics of Guangdong and Southern China during that period.
This work was financially supported by the national transfer payments funding for the 2011 influenza survey and the National Natural Science Foundation of China (Grants no. 30972757). The authors would like to thank Li H, Yang X and Ke C from the Guangdong Provincial Center for Disease Control and Prevention, and Chen Q at the School of Public Health and Tropic Medicine, Southern Medical University.
- Schweiger B, Zadow I, Heckler R: Antigenic drift and variability of influenza viruses. Med Microbiol Immunol 2002,191(3–4):133-138.PubMedGoogle Scholar
- Stray SJ, Pittman LB: Subtype- and antigenic site-specific differences in biophysical influences on evolution of influenza virus hemagglutinin. Virol J 2012, 9: 91. 10.1186/1743-422X-9-91PubMedPubMed CentralView ArticleGoogle Scholar
- Zhong J, Huang P, Wen MH, Liang LJ, Zhang X, Tan SN: Antigenic epitope peptides of influenza H3N2 virus neuraminidase gene based on experiments. Chin Sci Bull 2012,57(22):2908-2913. 10.1007/s11434-012-5266-0View ArticleGoogle Scholar
- Rickert JD, Herpy TA, Zhong H, Lu G: Computational analysis of adaptive antigenic mutations of the human influenza hemagglutinin for vaccine strain selection. Int J Bioinform Res Appl 2012,8(1–2):81-98.PubMedView ArticleGoogle Scholar
- Ann J, Papenburg J, Bouhy X, Rhéaume C, Hamelin MÈ, Boivin G: Molecular and antigenic evolution of human influenza A/H3N2 viruses in Quebec, Canada, 2009–2011. J Clin Virol 2012,53(1):88-92. 10.1016/j.jcv.2011.09.016PubMedView ArticleGoogle Scholar
- Huang P, Liang LJ, Hou NM, Zhang X, Su WZ, Yu SY: Phylogenetic, molecular and drug-sensitivity analysis of hemagglutinin and neuraminidase of human H3N2 influenza A viruses in Guangdong, China, 2007–2011. Epidemiol Infect 2013,141(5):1061-1069. 10.1017/S0950268812001318PubMedView ArticleGoogle Scholar
- Bedford T, Cobey S, Beerli P, Pascual M: Global migration dynamics underlie evolution and persistence of human influenza a (H3N2). PLoS Pathog 2010,6(5):e1000918. 10.1371/journal.ppat.1000918PubMedPubMed CentralView ArticleGoogle Scholar
- Center for Public Health Surveillance and Information Service (CPHSIS): National data of class A, B and C communicable diseases in June 2012. Dis Surveill 2012,27(7):505. ChineseGoogle Scholar
- Drummond AJ, Suchard MA, Xie D, Rambaut A: Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol 2012,29(8):1969-1973. 10.1093/molbev/mss075PubMedPubMed CentralView ArticleGoogle Scholar
- Delport W, Poon AF, Frost SD, Kosakovsky Pond SL: Datamonkey 2010: a suite of phylogenetic analysis tools for evolutionary biology. Bioinformatics 2010,26(19):2455-2457. 10.1093/bioinformatics/btq429PubMedPubMed CentralView ArticleGoogle Scholar
- Tomar AK, Sooch BS, Yadav S: Computational analysis of concanavalin a binding glycoproteins of human seminal plasma. Bioinformation 2011,7(2):69-75.PubMedPubMed CentralView ArticleGoogle Scholar
- Abro SH, Ullman K, Belák S, Baule C: Bioinformatics and evolutionary insight on the spike glycoprotein gene of QX-like and Massachusetts strains of infectious bronchitis virus. Virol J 2012, 9: 211. 10.1186/1743-422X-9-211PubMedPubMed CentralView ArticleGoogle Scholar
- Zou W, Ke J, Zhu J, Zhou H, Jin M: The antigenic property of the H5N1 avian influenza viruses isolated in central China. Virol J 2012, 9: 148. 10.1186/1743-422X-9-148PubMedPubMed CentralView ArticleGoogle Scholar
- Vita R, Zarebski L, Greenbaum JA, Emami H, Hoof I, Salimi N: The immune epitope database 2.0. Nucleic Acids Res 2010, 38: D854-D862. 10.1093/nar/gkp1004PubMedPubMed CentralView ArticleGoogle Scholar
- Roy A, Kucukural A, Zhang Y: I-TASSER: a unified platform for automated protein structure and function prediction. Nat Protoc 2010,5(4):725-738. 10.1038/nprot.2010.5PubMedPubMed CentralView ArticleGoogle Scholar
- Huang P, Yu S, Wu C, Liang L: Highly conserved antigenic epitope regions of hemagglutinin and neuraminidase genes between 2009 H1N1 and seasonal H1N1 influenza: vaccine considerations. J Transl Med 2013,11(1):47. 10.1186/1479-5876-11-47PubMedPubMed CentralView ArticleGoogle Scholar
- Rolland M, Carlson JM, Manocheewa S, Swain JV, Lanxon-Cookson E, Deng W: Amino-acid co-variation in HIV-1 Gag subtype C: HLA-mediated selection pressure and compensatory dynamics. PLoS One 2010,5(9):e12463. 10.1371/journal.pone.0012463PubMedPubMed CentralView ArticleGoogle Scholar
- Klimov AI, Garten R, Russell C, Barr IG, Besselaar TG, Daniels R: WHO recommendations for the viruses to be used in the 2012 Southern Hemisphere Influenza Vaccine: epidemiology, antigenic and genetic characteristics of influenza A(H1N1)pdm09, A(H3N2) and B influenza viruses collected from February to September 2011. Vaccine 2012,30(45):6461-6671. 10.1016/j.vaccine.2012.07.089PubMedView ArticleGoogle Scholar
- Frati E, Bianchi S, Colzani D, Zappa A, Orlando G, Tanzi E: Genetic variability in the major capsid L1 protein of human papillomavirus type 16 (HPV-16) and 18 (HPV-18). Infect Genet Evol 2011,11(8):2119-2124. 10.1016/j.meegid.2011.06.014PubMedView ArticleGoogle Scholar
- Das SR, Puigbò P, Hensley SE, Hurt DE, Bennink JR, Yewdell JW: Glycosylation focuses sequence variation in the influenza a virus H1 hemagglutinin globular site. PLoS Pathog 2010,6(11):e1001211. 10.1371/journal.ppat.1001211PubMedPubMed CentralView ArticleGoogle Scholar
- Sun S, Wang Q, Zhao F, Chen W, Li Z: Prediction of biological functions on glycosylation site migrations in human influenza H1N1 viruses. PLoS One 2012,7(2):e32119. 10.1371/journal.pone.0032119PubMedPubMed CentralView ArticleGoogle Scholar
- Suzuki Y: Positive selection for gains of N-linked glycosylation sites in hemagglutinin during evolution of H3N2 human influenza A virus. Genes Genet Syst 2011,86(5):287-294. 10.1266/ggs.86.287PubMedView ArticleGoogle Scholar
- Kim JI, Lee I, Park S, Hwang MW, Bae JY, Lee S, Heo J, Park MS, García-Sastre A, Park MS: Genetic requirement for hemagglutinin glycosylation and its implications for influenza A H1N1 virus evolution. J Virol 2013,87(13):7539-7549. 10.1128/JVI.00373-13PubMedPubMed CentralView ArticleGoogle Scholar
- Bulimo WD, Garner JL, Schnabel DC, Bedno SA, Njenga MK, Ochieng WO: Genetic analysis of H3N2 influenza a viruses isolated in 2006–2007 in Nairobi, Kenya. Influenza Other Respi Viruses 2008,2(3):107-113. 10.1111/j.1750-2659.2008.00044.xView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. 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.