- Research article
- Open Access
Amino acid substitutions in hemagglutinin of the 2009 pandemic influenza A(H1N1) viruses that might affect the viral antigenicity
BMC Research Notes volume 7, Article number: 951 (2014)
During 2009 to 2012, Thailand had encountered 4 distinctive waves of the 2009 pandemic influenza A(H1N1) (H1N1pdm) outbreaks. Considering the RNA nature of the influenza viral genome, a mutation in hemagglutinin (HA) gene which led to change in antigenicity of the strains circulating during those epidemic periods is anticipated. It is also uncertain whether the A/California/07/2009 (H1N1) (CA/07) vaccine strain still confers protective immunity against those evolved viruses, the causative agents of the later epidemic waves.
HA gene segments of 10 H1N1pdm isolates obtained during 2009 to 2012 were sequenced and phylogenetically analysed using ClustalW and MEGA5 programs. A total of 124 convalescent serum samples collected from patients naturally infected during 3 epidemic waves were employed as tools to investigate for antigenic change in HA of these 10 circulating H1N1pdm viruses by hemagglutination inhibition (HI) assay.
A phylogenetic analysis showed that the 10 virus isolates were grouped into 4 clusters corresponding to the time of 4 consecutive outbreaks. An accumulation of amino acid substitutions in HA was observed in viruses derived from the late epidemic waves. Significantly lower antibody titers were observed when CA/07 was tested against convalescent sera collected from the 3 waves (p < 0.05) compared to most of Thai isolates; and significantly lower antibody titers were also obtained when virus isolates, retrieved from the third epidemic wave were tested against convalescent sera collected during the first and second wave. These results were suggestive of change in antigenicity of the evolved viruses. Our results also showed some mutation position residing outside the previously reported antigenic site that may involve in an alteration of the viral antigenicity.
Our study demonstrated that convalescent sera collected from individuals naturally infected with H1N1pdm virus were successfully used to reveal a statistically significant change in antibody titers against the currently evolved H1N1pdm viruses as determined by HI assay. Nevertheless, the antibody titers of individual serum against various viruses were less than 4-folded difference as compared to that against the CA/07 vaccine strain. Therefore, CA/07 is still a potent vaccine strain for those evolved H1N1pdm viruses.
In mid-April 2009, the emergence of a novel influenza A virus was first noticed in Mexico . The virus is a quadruple reassortant, in which the RNA genome is originated from swine, avian, and human influenza viruses [2, 3]. The virus capable of infecting humans, at the time of its emergence, was antigenically new to the world’s population and subsequently spread and caused uncontained influenza outbreaks among humans worldwide. The World Health Organization (WHO) declared the influenza pandemic period between June 2009 and August 2010 [4, 5]. It involved over 18,449 laboratory-confirmed deaths reported to WHO , which probably under-represents the total number. On the other hand, an indirect estimation using statistical modeling suggested 201,200 respiratory deaths associated with the H1N1pdm during the first year of the virus circulation .
Thailand was among the firstly countries in Southeast Asia that was attacked by the 2009 pandemic influenza. The first documented case was recorded in May 2009 , followed by three subsequent dominant waves of the epidemics which lasted for 18 months. The Bureau of Epidemiology, Department of Disease Control, Ministry of Public Health, reported the first wave between May-October 2009, followed by the second in November 2009-April 2010, and the third in May-October 2010. Thailand was heavily attacked by H1N1pdm before the WHO announcement of the pandemic phase. In addition, the actual subsidence of the third wave in Thailand was noted some time after the WHO had announced the end of the pandemic. In January-June 2010, a monovalent pandemic H1N1 vaccine derived from A/California/07/2009 (H1N1) virus (CA/07) was administered to groups of people at risk in Thailand, including healthcare workers . Since then, the H1N1pdm virus has replaced the previous A(H1N1) virus and become the circulating strain in Thailand and worldwide. At present, CA/07 remains a component of the trivalent seasonal influenza vaccine together with influenza A(H3N2) and influenza B strains .
Considering the RNA nature of the influenza viral genome, a high rate of mutation resulting in antigenic drift is anticipated due to the lack of the proof-reading capacity of viral RNA polymerase. Influenza hemagglutinin (HA), the surface glycoprotein of a virion, comprises HA1 and HA2 domains. HA1 domain is highly variable; while HA2 is more conserved. The HA1 domain constitutes the HA globular head, which consists of Sa, Sb, Ca, and Cb antigenic sites . HA1 binds sialic acid receptor on the host cytoplasmic membrane, while HA2 mediates fusion of the viral envelope to the endosomal membrane. Thus, the HA molecule determines host cell specificity and induces a protective antibody response [11, 12]. Mutation in the HA gene sometimes results in the HA antigenic change and may result in changing the strain components of a seasonal trivalent vaccine. The hemagglutination inhibition (HI) assay is used to measure HA antibody titer for serodiagnosis or vaccine evaluation. HI antibody titers equal or greater than 40, either as the result of natural infection or vaccination, are considered immune to re-infection or protective against infection by homologous influenza strains . During the post-pandemic period, the World Health Organization has recommended monitoring of the virus for important genetic and antigenic changes . Knowing the changes will be useful for prediction of potential future outbreaks, and maintaining a matching between circulating influenza strains and vaccine strains.
In this study, 10 H1N1pdm strains isolated from Thai patients during the 4 epidemic waves (2009-2012) were analysed for genetic and antigenic variation in comparison with the CA/07 vaccine strain. A phylogenetic tree was constructed from their HA gene sequences. These 10 virus isolates were allowed to react with convalescent sera collected from patients infected during the first 3 epidemic waves; and a possible change in the virus antigenicity was analysed based on the levels of the HI antibody titers obtained. Although a panel of immune ferret sera immunized with particular strains is usually employed for determining the change in virus antigenicity by HI assay, our study suggested that convalescent sera collected during different epidemic waves from individuals naturally infected with H1N1pdm virus might serve as an alternative tool to reveal an antigenic variation of the evolved viruses. The 3-dimensional structure of the viral HA, with mutation, was also predicted by homology modeling.
Analysis of the HA sequences of the Thailand isolates
Detailed information of the 10 H1N1pdm isolates used in this study is shown in Table 1. Eight strains were isolated during the three epidemic waves which lasted from May 2009 to August 2010; and two strains were isolated during post-pandemic period in July 2012 (defined here as the fourth wave of the Thailand epidemic). Full-length HA genes of the study viruses were analyzed against CA/07 (vaccine strain) and other 21 H1N1pdm strains isolated in Thailand and neighboring countries during the years 2009-2013, using sequencing data retrieved from the NCBI influenza database (http://www.ncbi.nlm.nih.gov/genomes/FLU/). Multiple sequence alignment showed >97% identity to the HA genes of the 32 viruses analyzed. The Thailand viruses isolated from each epidemic wave were accordingly clustered together, corresponding with the times of epidemics (Figure 1). Although the viruses from wave 1 and wave 2 were not separated into distinctive clusters, a clear separation could be observed with the viruses from wave 3 and wave 4.
Amino acid sequences of HA protein, including the HA1 globular head (amino acid positions 42-275 which spanned the Sa, Sb, Ca, and Cb antigenic sites) and the HA2 domain of the study viruses were analyzed. The amino acid changes in those regions were compared with those of the vaccine strain, CA/07, as the reference sequence (Table 2). The result showed a higher number of amino acid substitutions in the viruses collected from the late epidemic waves than those from the early epidemic waves. As expected, mutations were more common in the HA1 domain than the HA2 domain.Considering the antigenic site, Sa, Sb, Ca, and Cb on the HA1 globular head, no difference was found in the amino acid compositions between Thai/104 (the imported virus from Mexico) and CA/07. However, amino acid substitution S203T on Ca was noted in all other Thai strains, while the mutation in Sb was detected in viruses isolated during the fourth epidemic wave. Surprisingly, no amino acid change was observed in the Sa and Cb sites. In addition, the BKK4 strain from the fourth wave contained amino acid substitution K119N, which contributes to the acquisition of a new N-linked glycosylation site. Therefore, it is predicted that BKK4 contains 7 glycosylation sites on the HA1 globular head instead of 6 sites as found in the other viruses. Amino acid substitutions that may affect antigenicity and N-linked glycosylation sites are demonstrated in the homology model of HA protein in Figure 2.
HI assay of sera from the H1N1pdm patients against various virus isolates
In this study, convalescent sera collected from laboratory-confirmed cases of H1N1pdm infection during the three epidemic waves were tested by HI assay against 11 virus isolates, including CA/07 virus and 10 Thailand isolates originating from the four epidemic waves. Using the CA/07 virus as the test antigen, the HI antibody titers ≥40 were found in 81.4% (48 of 59), 71.9% (23 of 32), and 78.8% (26 of 33) of the convalescent serum samples collected from the first, second and third epidemic waves, respectively. Variation in the number of sera with HI antibody titers ≥40 was found when these panels of sera were assayed against the Thailand viruses belonging to different epidemic waves (Table 3). The geometric mean titer (GMT) of HI antibody in sera collected from the 3 epidemic waves, as assayed against CA/07 and 10 Thailand isolates are shown in Figure 3. The GMT of sera from the first and second epidemic waves against CA/07 virus was significantly lower than GMTs against the Thailand viruses originating from the first and second epidemic waves (p < 0.05), but not against viruses from the third wave, which contained a higher number of amino acid substitution positions. Regarding serum samples from the third epidemic wave, the GMT against all viruses were higher than CA/07. Furthermore, BKK4 virus of the fourth wave showed a significantly lower GMT compared with the BKK5 virus of the same wave, when reacted with sera from all 3 waves.
HA is a viral protein that induces protective antibody response. With a pressure from neutralizing antibodies generated upon natural infection or vaccination, changes in the HA antigenic determinant are driven in order that the virus can evade the host immune response . It has been demonstrated that the percentages of amino acid divergence in HA are highest, when compared with other genome segments . In this study, HA sequences of 10 Thailand H1N1pdm viruses isolated in 2009-2012 were analysed genetically against sequencing data from the NCBI Influenza Virus Resource database. Phylogenetic analysis showed that the HA genes of our isolates formed 4 clusters, corresponding to the time of virus isolation, i.e., the first, second, third and fourth epidemic waves in Thailand. The accumulation of nucleotide substitutions in the late viruses indicated that the H1N1pdm viruses gradually evolved at the nucleic acid and amino acid levels during the past 5 years.
We sought to elucidate whether changes in the HA gene of our viruses affected the ability of HA protein to react with the anti-HA antibodies naturally raised in Thai patients. Convalescent sera from H1N1pdm patients infected during the 3 epidemic waves were assayed against 10 H1N1pdm viruses derived from the 4 epidemic waves by HI assay, the method with a potential to detect changes on the HA globular head. It was noticed that CA/07, the vaccine strain, yielded a statistically significant lower HI titers than most of the Thailand viruses originating from different epidemic waves; nevertheless, the antibody titer against CA/07 in individual serum was less than 4-folded difference as compared to those Thailand isolates (data not shown). Even though the cross reactive HI titers of ≥ 40 were found in 33% of adults of age more than 60 years prior to H1N1pdm vaccination , 2 (6.1%) of 33 of the acute sera from our subjects with unknown age possessed this antibody level (data not shown). Therefore, an issue on cross reactive antibody present in old aged people that might affect our laboratory interpretation is less likely in this study. P83S amino acid substitution in the HA1 globular head has been observed in all Thailand viruses, but not CA/07, suggesting that this point mutation may be associated with an increase in HI titer as investigated with the Thailand viruses. Nevertheless, P83S substitution is not located in the previously reported antigenic sites (Sa, Sb, Ca, and Cb). The sera from all 3 epidemic waves reacted with the viruses of the first and second waves at comparable titers, indicating that viruses of the first and the second waves are not antigenically distinct. This concurs with a report by the Ministry of Public Health, Thailand which showed that the attack rate of H1N1pdm in the second wave was lower than that of the first wave . This could be due to the protective immunity developed in Thai population after the first H1N1pdm epidemic wave. Based on amino acid analysis, the viruses from the third wave which accumulated a high number of amino acid substitutions, showed relatively poor activity when reacted with sera from the first and second waves. Observation of the HA1 globular head has demonstrated amino acid substitutions I216V and R205K in viruses from the third wave only, and this might relate to decreased HI titers when these viruses were used as test antigens. Nevertheless, the substitutions I216V are located outside previously reported antigenic sites; while R205K is located on antigenic site Ca. No substitution was found in the Sa and Cb regions of any virus studied. It is likely that amino acid residues outside the previously reported antigenic sites also play a role as additional antigenic determinants of viruses. Our finding was concordant with that reported by Deem et al. who employed an entropy-based method and demonstrated that more amino acid positions outside the previously defined antigenic site of H1 HA may involve in the HA antigenicity.
D97N, S185T, S203T, D222G, and E374K, observed in this study, have been reported in viruses from Thailand and other regions [16, 19, 20]. S203T present in all viruses except Thai/104 and S185T present in the viruses of the fourth wave, were located in the Ca and Sb sites, respectively. However, their involvement in antigenicity changes is uncertain. HA protein modeling suggested that the amino acid position 203 is not located on the surface of HA trimer. D222G, the mutation found in a virus isolate from a dead case (CBI10), is previously shown to be associated with the more virulent phenotype in humans and mice [19, 21, 22]. E374K substitution on HA2 stalk region has been observed in most of the H1N1pdm strains isolated after the year 2010 . This position, in particular, has been found in viruses isolated from vaccinees who received the monovalent 2009 H1N1 influenza vaccine, and from a number of fatal cases . The E374K is suggestive for an increase in a structural and acid stability of the HA trimers involving in membrane fusion, and potentially affects the viral antigenicity and adaptation [23, 25, 26]. Among 10 H1N1pdm viruses in this study, none of our 3 strains in the first epidemic wave contained this E374K substitution, while it was detected in all 7 strains isolated in the later epidemic wave. Four of these 7 strains were isolated from cases of vaccine failure (BKK1, BKK2, BKK4, and BKK5), but no case was fatal. The simultaneous occurrences of E374K with D97N and S203T, and of E374K with S185T among our Thailand viruses were previously observed among the strains isolated from the other geographical regions [19, 27]. K119N, a potential glycosylation site on the HA globular head, was found in BKK4, but not BKK5 virus derived from the fourth epidemic wave. This substitution may be responsible for a decrease in GMT when BKK4 was compared with the BKK5 virus. This K119N glycosylation site is situated close to the ligand binding site; hence, it could mask the receptor binding epitope on the HA surface . It was noted that K119N was found among the candidate vaccine strains with improved virus yields in eggs [28, 29].
Collectively, our results showed a gradual change in HA proteins of the H1N1pdm viruses isolated in Thailand during 2009-2012. This variation has resulted in the viral antigenic changes as a significant difference in HI antibody titer could be clearly demonstrated when convalescent sera were cross-checked with the vaccine strain and virus isolates derived from various epidemic waves. However, the HI titers found in this study are not more than 4-folded difference, indicating that the current CA/07 vaccine strain should still be effective.
Subjects and samples
A total of 124 archival samples of convalescent sera comprised 59 samples collected from male conscripts in July 2009 during the first epidemic wave, 32 samples (15 male conscripts and 17 general population of age between 2-42 years old) collected between December 2009 and March 2010 during the second epidemic wave, and 33 samples collected from male inmates in September 2010 during the third epidemic wave. These sera had been kept at -20°C until tested. As suggested by occupation, the ages of male conscripts were approximately 21 years old. Unfortunately, ages of some subjects including the male inmates were not known. Infection status of the subjects was confirmed by detection of the viral genome in respiratory samples by real time reverse transcription-polymerase chain reaction (real time RT-PCR) using the protocol established by the U.S. Centers for Disease Control and Prevention (CDC) . Collection of respiratory and blood samples is the standard practice for investigation of influenza outbreak by the Bureau of Epidemiology, Ministry of Public Health. This study which employed these blood samples was approved by Siriraj Institutional Review Board, Faculty of Medicine Siriraj Hospital, Mahidol University, Thailand.
Eleven strains of H1N1pdm, isolated during the years 2009-2012, were employed as test antigens in HI assay. Except for CA/07, which was obtained from the National Institute for Biological Standards and Control (NIBSC), UK, all viruses were isolated in Thailand. The HA gene of these viruses were subjected to nucleotide sequencing, and the complete HA sequences were submitted to the NCBI Influenza Virus Resource database (http://www.ncbi.nlm.nih.gov/genomes/FLU/). Descriptions of the viruses are shown in Table 1. The viruses were propagated in Madin-Darby Canine Kidney (MDCK) cells, obtained from the American Type Culture Collection (CCL-34), maintained in Eagle's minimum essential medium (Invitrogen, Grand Island, NY) supplemented with 2 μg/ml of L-1-tosylamido-2-phenylethyl chloromethyl ketone (TPCK)-treated trypsin (Sigma, St. Louis, MO). The culture supernatant was harvested, centrifuged at 4°C for 15 minutes, aliquoted, and kept frozen as virus stock at -80°C until used.
Nucleotide sequencing and analysis
Total RNA was extracted from the virus stock by QIAamp® Viral RNA Mini Kit (QIAGEN Inc., Valencia, CA) using the procedure described in the manufacturer’s instructions. The HA gene was amplified by QIAGEN® OneStep RT-PCR kit (QIAGEN Inc.) using 2 overlapping HA primer sets: HA1F; 5’ATACGACTAGCAAAAGCAGGGG3’ and HA943R; 5’GAAAKGGGAGRCTGGTGTTTA3’ (product size of 943 bp) and HA736F; 5’AGRATGRACTATTACTGGAC3’ and HA1778R; 5’GTGTCAGTAGAAACAAGGGTGTTT3’ (product size 1042 bp). The protocol for RT-PCR comprised the steps of reverse transcription at 50°C for 45 min and PCR activation at 95°C for 15 min, followed by 5 cycles of denaturation at 95°C for 30 sec, annealing at 50°C for 30 sec, and extension at 72°C for 1 min, then another 30 cycles of denaturation at 95°C for 30 sec, annealing at 55°C for 30 sec and extension at 72°C for 1 min, with the final extension step at 72°C for 10 min. The PCR products were gel purified and sent to BioDesign Co., Ltd., Thailand for nucleotide sequencing. Nucleotide sequences of the 2 partial HA fragments were assembled using the BioEdit program to yield the full length HA sequence. A phylogenetic tree was constructed using the MEGA 5 program, by applying the neighbor-joining algorithms with Kimura’s two-parameter distance model and 1,000 bootstrap replicates [16, 31]. Prediction of the potential glycosylation site employed NetNGlyc 1.0 (http://www.cbs.dtu.dk/services/NetNGlyc/). The glycosylation site was defined as Asn-X-Ser/Thr, where X represents any amino acid, except proline .
Hemagglutination inhibition (HI) assay
The HI assay was performed as previously described [32, 33]. The test sera were mixed with receptor-destroying enzyme (RDE) (Denka Seiken, Tokyo, Japan) at a serum: RDE ratio of 1: 3, to remove non-specific inhibitor. The mixture was incubated at 37°C in a water bath for 16-18 h, followed by heat inactivation at 56°C for 30 min. The RDE-treated serum was added with 5 volumes of 0.9% NaCl and 1 volume of 50% erythrocyte followed by incubation at 4°C for 1 h, to remove non-specific agglutinator. The reaction tube was spun and the treated serum at a dilution of 1: 10 was obtained. For the HI assay, the test serum was 2-fold serially diluted from a starting dilution of 1: 10 to 1: 1280 in a 25 μl volume, and 25 μl of the virus suspension were added to a final concentration of 4 HA units. After incubation at ambient temperature for 30 min, 50 μl of 0.5% goose erythrocytes were added and the mixture was incubated at 4°C for 30 min. HI antibody titers were determined and defined as the reciprocal of the highest serum dilution that completely inhibited the hemagglutination reaction. Each run included positive control serum and back titration of viral antigen. The experiments were run in duplicate and the results were presented as geometric mean titer (GMT). The difference in HI titers of the test sera as assayed against different virus strains was statistically analyzed by Wilcoxon Signed Ranks Test. The statistical significance was set at p < 0.05.
A structural model of the influenza HA molecule was constructed using SWISS-MODEL, based on homology modeling using the crystal structure of the HA molecule of H1N1pdm (PDB ID: 3LZG) as template . The model obtained was validated by PROCHECK .
Centers for Disease Control and Prevention: US CDC (2009) Update: infections with a swine-origin influenza A (H1N1) virus-United States and other countries, April 28, 2009. MMWR. 2009, 58: 431-433.
Itoh Y, Shinya K, Kiso M, Watanabe T, Sakoda Y, Hatta M, Muramoto Y, Tamura D, Sakai-Tagawa Y, Noda T, Sakabe S, Imai M, Hatta Y, Watanabe S, Li C, Yamada S, Fujii K, Murakami S, Imai H, Kakugawa S, Ito M, Takano R, Iwatsuki-Horimoto K, Shimojima M, Horimoto T, Goto H, Takahashi K, Makino A, Ishigaki H, Nakayama M: In vitro and in vivo characterization of new swine-origin H1N1 influenza viruses. Nature. 2009, 460: 1021-1025.
Neumann G, Noda T, Kawaoka Y: Emergence and pandemic potential of swine-origin H1N1 influenza virus. Nature. 2009, 459: 931-939. 10.1038/nature08157.
World Health Organization: Statement to the press by WHO Director-General Dr Margaret Chan: World now at the start of. 2009, [http://www.who.int/mediacentre/news/statements/2009/h1n1_pandemic_phase6_20090611/en/], influenza pandemic
World Health Organization: Director-General's opening statement at virtual press conference: H1N1 in post-pandemic period. [http://www.who.int/mediacentre/news/statements/2010/h1n1_vpc_20100810/en/]
World Health Organization: Pandemic (H1N1) 2009 - update 112. [http://www.who.int/csr/don/2010_08_06/en/]
Dawood FS, Iuliano AD, Reed C, Meltzer MI, Shay DK, Cheng PY, Bandaranayake D, Breiman RF, Brooks WA, Buchy P, Feikin DR, Fowler KB, Gordon A, Hien NT, Horby P, Huang QS, Katz MA, Krishnan A, Lal R, Montgomery JM, Mølbak K, Pebody R, Presanis AM, Razuri H, Steens A, Tinoco YO, Wallinga J, Yu H, Vong S, Bresee J: Estimated global mortality associated with the first 12 months of 2009 pandemic influenza A H1N1 virus circulation: a modelling study. Lancet Infect Dis. 2012, 12: 687-695. 10.1016/S1473-3099(12)70121-4. Erratum in: Lancet Infect Dis 2012, 12:655
Iamsirithaworn S, Akarasewi P, Yingyong T, Suthachana S, Pittayawonganon C, Ungchusak K: Three waves of the 2009 H1N1 influenza pandemic in Thailand. Siriraj Med J. 2011, 63: 64-67.
World Health Organization: WHO recommendations on the composition of influenza virus vaccines. [http://www.who.int/influenza/vaccines/virus/recommendations/en/]
Igarashi M, Ito K, Yoshida R, Tomabechi D, Kida H, Takada A: Predicting the antigenic structure of the pandemic (H1N1) 2009 influenza virus hemagglutinin. PLoS One. 2010, 5: e8553-10.1371/journal.pone.0008553.
Skehel JJ, Wiley DC: Receptor binding and membrane fusion in virus entry: the influenza hemagglutinin. Annu Rev Biochem. 2000, 69: 531-569. 10.1146/annurev.biochem.69.1.531.
Maines TR, Szretter KJ, Perrone L, Belser JA, Bright RA, Zeng H, Tumpey TM, Katz JM: Pathogenesis of emerging avian influenza viruses in mammals and the host innate immune response. Immunol Rev. 2008, 225: 68-84. 10.1111/j.1600-065X.2008.00690.x.
de Jong JC, Palache AM, Beyer WE, Rimmelzwaan GF, Boon AC, Osterhaus AD: Haemagglutination-inhibiting antibody to influenza virus. Dev Biol (Basel). 2003, 115: 63-73.
World Health Organization: WHO recommendations for the post-pandemic period. [http://www.who.int/csr/disease/swineflu/notes/briefing_20100810/en/]
Schmolke M, García-Sastre A: Evasion of innate and adaptive immune responses by influenza A virus. Cell Microbiol. 2010, 12: 873-880. 10.1111/j.1462-5822.2010.01475.x.
Makkoch J, Suwannakarn K, Payungporn S, Prachayangprecha S, Cheiocharnsin T, Linsuwanon P, Linsuwanon PTheamboonlers A, Poovorawan Y: Whole genome characterization, phylogenetic and genome signature analysis of human pandemic H1N1 virus in Thailand, 2009-2012. PLoS One. 2012, 7: e51275-10.1371/journal.pone.0051275.
Centers for Disease Control and Prevention (CDC): Serum cross-reactive antibody response to a novel influenza A (H1N1) virus after vaccination with seasonal influenza vaccine. MMWR Morb Mortal Wkly Rep. 2009, 58: 521-524.
Deem MW, Pan K: The epitope regions of H1-subtype influenza A, with application to vaccine efficacy. Protein Eng Des Sel. 2009, 22: 543-546. 10.1093/protein/gzp027.
Ikonen N, Haanpää M, Rönkkö E, Lyytikäinen O, Kuusi M, Ruutu P, Kallio-Kokko H, Mannonen L, Lappalainen M, Ziegler T, Julkunen I: Genetic diversity of the 2009 pandemic influenza A(H1N1) viruses in Finland. PLoS One. 2010, 5: e13329-10.1371/journal.pone.0013329.
Strengell M, Ikonen N, Ziegler T, Julkunen I: Minor changes in the hemagglutinin of influenza A(H1N1)2009 virus alter its antigenic properties. PLoS One. 2011, 6: e25848-10.1371/journal.pone.0025848.
Kilander A, Rykkvin R, Dudman SG, Hungnes O: Observed association between the HA1 mutation D222G in the 2009 pandemic influenza A(H1N1) virus and severe clinical outcome, Norway 2009-2010. Euro Surveill. 2010, 15: 19498-
Mak GC, Au KW, Tai LS, Chuang KC, Cheng KC, Shiu TC, Lim W: Association of D222G substitution in haemagglutinin of 2009 pandemic influenza A (H1N1) with severe disease. Euro Surveill. 2010, 15: 19534-
Yang H, Chang JC, Guo Z, Carney PJ, Shore DA, Donis RO, Cox NJ, Villanueva JM, Klimov AI, Stevens J: Structural stability of influenza A(H1N1)pdm09 virus hemagglutinins. J Virol. 2014, 88: 4828-4838. 10.1128/JVI.02278-13.
Barr IG, Cui L, Komadina N, Lee RT, Lin RT, Deng Y, Caldwell N, Shaw R, Maurer-Stroh S: A new pandemic influenza A(H1N1) genetic variant predominated in the winter 2010 influenza season in Australia, New Zealand and Singapore. Euro Surveill. 2010, 15: 19692-
Maurer-Stroh S, Lee RT, Eisenhaber F, Cui L, Phuah SP, Lin RT: A new common mutation in the hemagglutinin of the 2009 (H1N1) influenza A virus. PLoS Curr. 2010, 2: RRN1162-
Cotter CR, Jin H, Chen Z: A single amino acid in the stalk region of the H1N1pdm influenza virus HA protein affects viral fusion, stability and infectivity. PLoS Pathog. 2014, 10: e1003831-10.1371/journal.ppat.1003831.
Mak GC, Leung CK, Cheng KC, Wong KY, Lim W: Evolution of the haemagglutinin gene of the influenza A(H1N1)2009 virus isolated in Hong Kong, 2009-2011. Euro Surveill. 2011, 16: 19807-
Robertson JS, Nicolson C, Harvey R, Johnson R, Major D, Guilfoyle K, Roseby S, Newman R, Collin R, Wallis C, Engelhardt OG, Wood JM, Le J, Manojkumar R, Pokorny BA, Silverman J, Devis R, Bucher D, Verity E, Agius C, Camuglia S, Ong C, Rockman S, Curtis A, Schoofs P, Zoueva O, Xie H, Li X, Lin Z, Ye Z: The development of vaccine viruses against pandemic A(H1N1) influenza. Vaccine. 2011, 29: 1836-1843. 10.1016/j.vaccine.2010.12.044.
Nicolson C, Harvey R, Johnson R, Guilfoyle K, Engelhardt OG, Robertson JS: An additional oligosaccharide moiety in the HA of a pandemic influenza H1N1 candidate vaccine virus confers increased antigen yield in eggs. Vaccine. 2012, 30: 745-751. 10.1016/j.vaccine.2011.11.081.
World Health Organization: CDC protocol of realtime RTPCR for influenza A(H1N1). [http://www.who.int/csr/resources/publications/swineflu/CDCRealtimeRTPCR_SwineH1Assay-2009_20090430.pdf?ua=1]
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.
World Health Organization: WHO manual on animal influenza diagnosis and surveillance. [http://www.who.int/csr/resources/publications/influenza/en/whocdscsrncs20025rev.pdf]
Lerdsamran H, Pittayawonganon C, Pooruk P, Mungaomklang A, Iamsirithaworn S, Thongcharoen P, Kositanont U, Auewarakul P, Chokephaibulkit K, Oota S, Pongkankham W, Silaporn P, Komolsiri S, Noisumdaeng P, Chotpitayasunondh T, Sangsajja C, Wiriyarat W, Louisirirotchanakul S, Puthavathana P: Serological response to the 2009 pandemic influenza A (H1N1) virus for disease diagnosis and estimating the infection rate in Thai population. PLoS One. 2011, 6: e16164-10.1371/journal.pone.0016164.
Schwede T, Kopp J, Guex N, Peitsch MC: SWISS-MODEL: an automated protein homology-modeling server. Nucleic Acids Res. 2003, 31: 3381-3385. 10.1093/nar/gkg520.
Laskowski RA, MacArthur MW, Moss DS, Thornton JM: PROCHECK: a program to check the stereochemical quality of protein structures. J Appl Cryst. 1993, 26: 283-291. 10.1107/S0021889892009944.
This work was supported by the Thailand Research Fund for Senior Research Scholar through PP, the Mahidol University Post-doctoral Scholarship, and the Office of the Higher Education Commission and Mahidol University under the National Research Universities Initiative. The authors thank Ms. Rungnapa Bunruang and Ms. Kannika Nateerom for assistance with the HI assay. The manuscript was grammatically edited by Mr. Glad Rotaru, the Office of Research Services, Faculty of Tropical Medicine, Mahidol University.
The authors declare that they have no competing interests.
NK carried out the HA phylogenetic study and hemagglutination inhibition assay, analysed data and wrote the manuscript. UB performed the protein homology modeling. PhP preformed the virus isolation. SI, AM and KC collected samples. PA coordinated the study. PP conceived of the study, designed the study and wrote the manuscript. All authors read and approved the final manuscript.