Journal of Clinical Microbiology, February 2008, p. 399-405, Vol. 46, No. 2
0095-1137/08/$08.00+0 doi:10.1128/JCM.01549-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Department of Public Health, Graduate School of Medical and Dental Sciences, Niigata University, Niigata, Japan,1 National Institute of Hygiene and Epidemiology, Hanoi, Vietnam,2 Niigata Prefectural Institute of Public Health and Environmental Sciences, Niigata City, Niigata Prefecture, Japan3
Received 3 August 2007/ Returned for modification 13 September 2007/ Accepted 10 October 2007
|
|
|---|
|
|
|---|
The effectiveness of annually applied trivalent influenza vaccines depends on the selection of component strains that offer optimal immunity from the numerous variants in the global influenza virus circulation. Studies based on sequencing analyses of viruses can be utilized as surveillance tools and can contribute to the vaccine selection process when they are combined with classical serological antigenic analysis (10). Continuous monitoring of viral genetic changes throughout the year is necessary for us to develop our ability to precisely define variation in influenza virus.
In the tropical zone, influenza virus circulates throughout the year, as reported in southern China, Indonesia, and Thailand (2, 18, 20). A laboratory-based influenza virus surveillance system has been in place in Hanoi, Vietnam, since 2001, and we have reported on the presence of influenza virus throughout the year, albeit with summer and winter peaks (21). In order to elucidate the evolutionary patterns for influenza viruses in Vietnam, in this study we undertook a genetic analysis of the influenza A viruses circulating from 2001 to 2006, focusing on two external genes (the HA gene and the NA gene) and one internal regulatory gene (the M gene, which codes for the M1 and the M2 channel proteins).
|
|
|---|
RNA extraction and PCR. One hundred-microliter aliquots of the supernatants after the third culture passage were used for viral RNA extraction with an Extragen II kit (Kainos, Tokyo, Japan), according to the manufacturer's instructions. RNA was transcribed to cDNA with the influenza A virus universal primer Uni12, as described elsewhere (24). The HA genes (segment 4), the NA genes (segment 6), and the M genes (segment 7) of H1N1 and H3N2 viruses were amplified with segment-specific primers as described elsewhere (3, 24).
Nucleotide sequencing and phylogenetic analysis. The PCR products were purified with a MicroSpin S-300 HR column PCR purification kit (Amersham Bioscience, Buckinghamshire, United Kingdom); labeled by use of a BigDye Terminator (version 3.1) cycle sequencing kit (Applied Biosystems, Foster, CA), according to the manufacturer's instructions; and then analyzed on an ABI 3100 automatic DNA sequencer. The sequences were assembled by use of the MEGA (version 3.1) program (26), and multiple-sequence alignment was conducted with the Clustal W program for the major coding regions of the three segments: HA1 (906 bp), NA (1,363 bp), and M (with regions of the overlapping reading frames of M1-M2; 923 bp) for H3N2 isolates and HA1 (829 bp), NA (1,363 bp), and M (938 bp) for H1N1 isolates. Phylogenetic trees were constructed by using a neighbor-joining and bootstrap analysis (n = 1,000) program to determine the best fits for the HA, NA, and M genes. Major branches with bootstrap values of >70% were identified as distinct groups. The phylogenetic grouping was not consistent for the three genes. The genomic sequences of the vaccine strains and the other strains used in this study were obtained from the Influenza Sequences Database (http://www.flu.lanl.gov).
Nucleotide sequence accession numbers. The nucleotide sequence data from this study were deposited in the DDBJ (DNA Data Bank of Japan), with the accession numbers listed in Table 1.
|
View this table: [in a new window] |
TABLE 1. Sequence data of influenza isolates used in this study
|
|
|
|---|
A/H3N2 influenza virus. For analysis of the sequences of the HA, NA, and M genes of the influenza A/H3N2 viruses, we used the consensus sequence of A/Moscow/10/99 as the phylogenetic root for the HA gene. Other isolates were also used as reference strains for each year's epidemic.
In the phylogenetic tree of the HA genes of the A/H3N2 isolates, the viruses formed a monophyletic group which could be divided into four major subgroups with a bootstrap value of >70% (Fig. 1A). These four groups were concurrent with the influenza seasons: group I in 2002; group II in 2003; group III in 2003, 2004, and 2005; and group IV in 2005. Six isolates circulating during the 2002 season in Hanoi were in positions between the A/Panama/2007/99 and the A/Fujian/411/02 groups of viruses and were clustered in group I. One 2003 isolate and the vaccine strain (A/Fujian/411/02) were clustered in group II. Two isolates circulating in 2003, nine isolates circulating in 2004, and strain A/Wellington/1/04 were classified as group IIIa; and three isolates circulating in 2005 and strain A/California/7/04 were classified as group IIIb. The remaining nine isolates circulating in 2005 and the Northern Hemisphere vaccine strain A/Wisconsin/67/05 clustered in group IV.
![]() View larger version (48K): [in a new window] |
FIG. 1. Phylogenetic analysis of the HA1 domain of HA gene nucleotide sequences (906 bp), NA gene nucleotide sequences (1,363 bp), and M gene (with regions of the overlapping reading frames of M1-M2) nucleotide sequences (923 bp) of influenza A/H3N2 viruses circulating in Vietnam from 2002 to 2005. Reference strains shown in black were obtained from the genetic database. The isolates were assigned to group I (orange), group II (green), group III (purple), or group IV (red) in the phylogenetic tree of the HA gene (A). The grouping of the NA and the M genes was made in accordance with the branching in each tree. The colors of the strains in the phylogenetic trees of the NA gene (B) and the M gene (C) comply with those for the HA gene. Sequence data for the reference strains were obtained from the GenBank database. The asterisks denote the amantadine-resistant strain with the amino acid change S31N in the M2 protein. Bootstrap values of >70% are shown for the main groups.
|
In the phylogenetic tree of the M gene (with regions of the overlapping reading frames of M1-M2), the viruses were divided into two groups. Isolates circulating in 2002, 2003, and 2005 were clustered in group I (group I, group II, and group IV for HA) (Fig. 1C). Nine isolates circulating in 2004 and three isolates circulating in 2005 were clustered in group II (group III for HA).
The amino acid sequences encoded by the HA (HA1 subunit) genes are shown in Table 2. For group I, all six isolates had amino acid changes, at positions S21P, R50G, E83K, N145K, S186G, V202I, W222R, and G225D, in comparison with the sequence of A/Panama/2007/99 strain. Four of the six isolates had additional double changes at positions N144D and G275D, and two of the six had an additional amino acid difference at position A131T. The single isolate circulating in 2003 and some other reference strains circulating in the same year were in group II and had four additional amino acid changes, L25I, H75Q, H155T, and Q156H, compared to the sequence of the isolates in group I. Among the isolates in group IIIa, 11 had additional changes at Y159F, S189N, and S227P and 4 of the 11 had another change at V226I, while among the isolates in group IIIb, three isolates circulating in 2005 further changed at K145N. All the remaining nine isolates in 2005, which were in group IV, had double amino acid changes at S193F and D225N. All isolates in this group were amantadine resistant and had the S31N amino acid change in the M2 gene.
|
View this table: [in a new window] |
TABLE 2. Comparison of amino acid sequences for the HA1 subunit of H3N2 viruses collected from 2002 to 2005
|
![]() View larger version (35K): [in a new window] |
FIG. 2. Phylogenetic analysis of the HA1 domain of HA gene nucleotide sequences (829 bp), NA gene nucleotide sequences (1,363 bp), and M gene (with regions of the overlapping reading frames of M1-M2) nucleotide sequences (938 bp) of influenza A/H1N1 viruses circulating in Vietnam from 2001 to 2006. Reference strains shown in black were obtained from the genetic database. The isolates were assigned to group I (orange), group II (dark red), group III (purple), or group IV (red) in the phylogenetic tree of HA gene (Fig. 1A). The grouping of the NA and M genes was made in accordance with the branching in each tree. The colors of the strains in the phylogenetic trees of the NA gene (B) and the M gene (C) comply with those for the HA gene. The asterisks denote the amantadine-resistant strain with the amino acid change S31N in the M2 protein. Bootstrap values of >70% are shown for the main groups.
|
In the phylogenetic analysis of the M genes (with regions of the overlapping reading frames of M1-M2) of the viruses, all except two isolates from 2001 to 2006 clustered in group I (Fig. 2C); isolates A/Hanoi/BM344/06 and A/Hanoi/TX09/06 clustered in group II with a bootstrap value of 89%.
In comparison to vaccine strain A/New Caledonia/20/99, the HA genes of the A/H1N1 isolates in groups I, II, and III from 2002 to 2006 had mutations at positions 158, 169, 190, 255, and 256. In particular, two amino acid changes at V169A and W255R consistently occurred in the samples (Table 3). At position 158, isolates in group II from 2002 and 2003 demonstrated an N158S change, while an isolate in group III from 2005 had a change at N158K. At position 190, isolates in groups I, II, and III demonstrated N190D and N190V changes. At position 256, isolates in group III from 2005 and 2006 featured Y256F. Isolate A/Hanoi/BM344/06 of group IV, on the other hand, exhibited 12 amino acid changes in comparison with the sequence of A/New Calidonia/20/99; these were E72D, K77R, T86K, Y98H, V132T, K144T, R149K, V169A, V190D, R212K, W255R, and T270N.
|
View this table: [in a new window] |
TABLE 3. Comparison of amino acid sequences for the HA1 subunit of H1N1 isolates collected from 2002 to 2006
|
|
|
|---|
During the 5 years, influenza A/H3N2 virus evolved from a strain intermediate between A/Panama/2007/99 and A/Fujian/411/02 to an A/Wiscosin/67/05-like strain. The HA and NA genes evolved independently in a sequential way, whereas the M gene demonstrated a nonlinear pattern not linked to the evolution of the HA and NA genes (29). Furthermore, the HA gene could be divided into four groups with shifts from groups I to IV. The NA and M genes were divided into two groups, but the grouping of the M genes was not in accordance with that of the NA or the HA gene. The results also suggest that the HA gene evolved more rapidly than both the NA and the M genes and that these three genes of influenza A/H3N2 virus change in a nonsynchronous manner.
The genetic drift of the HA genes of the influenza A/H3N2 viruses, located at positions believed to have functional or antigenic significance, was found to have occurred during all epidemic periods. The HA genes of viruses isolated in 2002 were clustered between strains A/Panama/2007/99 and A/Fujian/411/2002 in the phylogeny, with some strains having five amino acid changes within the antigenic sites (38, 39), namely, R50G (site C), E83K (site E), A131T (site A), N144K (site A), and S186G (site D). The same variant began to circulate in other parts of Asia between the end of 2001 and the beginning of 2002 (9). One isolate that was collected in August 2003 showed two key amino acid changes at residues H155T and Q156T (site B). This result indicated invasion of the A/Fujian/411/02 like strain into Vietnam, which was roughly coincident with our antigenic findings from influenza virus surveys in Hanoi (21). Between September 2003 and February 2004, the predominant strain circulating in Asia, the Americas, and Europe was an A/Fujian/411/02-like strain (40). However, two Vietnamese isolates circulating in September 2003 had drifted from the A/Fujian/411/02-like strain and showed additional changes at residues 159 (site B), 189 (site B), and 227 (site D). They had a higher degree of genetic homology to strain A/Wellington/1/04, which was isolated in New Zealand in January 2004 and which also circulated in Vietnam in the winter season of 2003-2004. Moreover, strain A/California/7/04 like, characterized by further changes at residues 145 (site A) and 226 (11, 30, 31), was isolated in Vietnam in May 2005 in this study.
High proportions of amantadine-resistant influenza A/H3N2 viruses with specific amino acid changes at residues 193 and 225 (named clade N) were reported from the 2005-2006 season in Japan and the United States (4, 12, 34, 35). Residue 193 is located within antigenic site B (37), and residue 225 is located within the receptor-binding site. In this study, clade N viruses were found from April to May of 2005 in Vietnam and the surrounding areas (1). However, it was in September and October of 2005 and the subsequent winter in Japan and North America when the strain caused a large community outbreak (5, 12, 34, 35).
The HA genes of the A/H1N1 isolates in groups I, II, and III showed only five amino acid changes from the sequence of the A/New Caledonia/20/99 vaccine strain. The HA1 subunit includes the globular head and contains five major antibody-binding sites, namely, Sa, Sb, Ca1, Ca2, and Cb (7). A substitution at residue 169 (V169A) occurred in the Ca1 antigenic site in all isolates, while substitution at residue 158, located in the Sb site, changed irregularly with the strain. The WHO Influenza Center reported that the majority of H1N1 viruses circulating in 2005 were antigenically closely related to current vaccine strain A/New Caledonia/20/99 (22). One of the strains analyzed at the WHO center, A/Virginia/4/05 (22), is antigenically related to strain A/New Caledonia/20/99 and grouped together with our isolates in group III that were circulating in Vietnam in 2005 and 2006. Therefore, we conclude that the most of the H1N1 viruses circulated in Vietnam from 2001 to 2006 had antigenicity similar to that of strain A/New Caledonia/20/99 and underwent changes more slowly than the influenza A/H3N2 virus subtype, as reported elsewhere (15).
An A/Hanoi/BM344/06 strain collected in May 2006 showed 12 additional amino acid changes in the HA gene in comparison with the sequence of strain A/New Caledonia/20/99. Two of the 12 substitutions, K81R and V169A, were located in antigenic sites Cb and Ca1, respectively. In addition, substitutions at V133T and K144E were located in antibody-binding sites (38). Also, strain A/Solomon Islands/3/06, which was isolated in August 2006 and which was later recommended by WHO as the vaccine strain for the 2007-2008 season (41), clustered with strain A/Hanoi/BM344/06 in the phylogenetic tree of the HA gene. Therefore, it is likely that we obtained an early insight into the arrival of this novel strain of influenza A/H1N1 virus.
In general, reassortment events are reported more often in A/H3N2 viruses than in A/H1N1 viruses (25). Our analysis indicated that reassortment occurred for both H1N1 and H3N2 during this 6-year period. Isolate A/Hanoi/184/02 (H3N2) contained an HA gene intermediate between those of A/Panama/2007/99 and A/Fujian/411/02 and an NA gene similar to that of A/Wyoming/3/03, whereas the M gene proved similar to that of A/Fujian/411/02. For H1N1, isolate A/Hanoi/TX09/06 (H1N1) represents the result of a rarely reported reassortment. The HA gene of the virus was closely related to those in group III, which were antigenically related to the A/New Caledonia/20/99-like strain, but the NA and M genes were closely related to the A/Solomon Islands/3/06 strain. An intersubtype reassortant, A/H1N2, was isolated from Japan as well as from many other areas (13, 19, 32, 42) after 2001, but this was not detected in our present analysis.
Continuous monitoring of viral genetic changes throughout the year is warranted to monitor the variations of influenza viruses. As influenza viruses circulate throughout the year in the tropical and subtropical zones (18, 20), the results for viruses from those areas are very suitable for monitoring purposes. This applies especially in Asia, where the emergence of new strains may be detected earlier than the emergence of new strains in temperate zones, where there is only a single peak of activity in each year.
We declare that none of the authors have any conflict of interest.
This work was supported by Grants-in-Aid for Scientific Research from Monbu Kagakusho (Ministry of Education, Culture, Sports, Science and Technology, Japan) and by the Acute Respiratory Infections Panels, United States-Japan Cooperative Medical Science Program (which is supported by the U.S. Department of Health and Human Services and the U.S. Department of State in the United States and the Ministry of Foreign Affairs; the Ministry of Health, Labor, and Welfare; and the Ministry of Education, Culture, Sports, Science and Technology in Japan).
Published ahead of print on 17 October 2007. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»