key: cord-0907659-eh1k2c33 authors: Wang, Sheng‐Fan; Lee, Yuan‐Ming; Chan, Yu‐Jiun; Liu, Hsin‐Fu; Yen, Yung‐Fong; Liu, Wu‐Tse; Huang, Jason C.; Chen, Yi‐Ming Arthur title: Influenza A virus in Taiwan, 1980–2006: Phylogenetic and antigenic characteristics of the hemagglutinin gene date: 2009-06-23 journal: J Med Virol DOI: 10.1002/jmv.21531 sha: ee1859efc506500f2f4b82530938ed41220417c1 doc_id: 907659 cord_uid: eh1k2c33 Limited amount of information is available in Taiwan on the genetic or antigenic characteristics of influenza A virus prior to the establishment of a Taiwan surveillance network in 2000. Isolates of H1N1 and H3N2 viruses in Taiwan between 1980 and 2006 were studied, and part of the hemagglutinin gene was analyzed due to its importance in terms of viral infection and antibody neutralization. Results from a phylogenetic analysis indicate continuous evolutionary topology in H3N2 isolates, and two distinct H1N1 lineages. Many genetic relationships between vaccine strains and epidemic isolates appearing in Taiwan before other global locations were also observed and recorded in addition to a gradual increase in the number of N‐linked glycosylation sites on partial HA1 proteins since 1980. The results from pairwise comparisons of HA1 nucleotide and deduced amino acid sequences indicate shared identities within groups organized according to their bootstrap and P‐values of approximately 95.5–100% and 95.7–100% in H1N1 and 94.5–100% and 93.2–100% in H3N2 viruses, respectively. Comparisons of amino acid substitutions in the five antigenic regions reveal highly non‐synonymous changes occurring in the Sb region of H1N1 and in the B region of H3N2. The results of an antigenic analysis using a hemagglutinin inhibition (HI) test indicate the presence of some epidemic strains 1–2 years earlier in Taiwan than in other parts of the world, as well as higher vaccine mismatch rates. This information supports the need for continuous surveillance of emerging influenza viruses in Taiwan, which will be useful for making global vaccine decisions. J. Med. Virol. 81:1457–1470, 2009. © 2009 Wiley‐Liss, Inc. Limited amount of information is available in Taiwan on the genetic or antigenic characteristics of influenza A virus prior to the establishment of a Taiwan surveillance network in 2000. Isolates of H1N1 and H3N2 viruses in Taiwan between 1980 and 2006 were studied, and part of the hemagglutinin gene was analyzed due to its importance in terms of viral infection and antibody neutralization. Results from a phylogenetic analysis indicate continuous evolutionary topology in H3N2 isolates, and two distinct H1N1 lineages. Many genetic relationships between vaccine strains and epidemic isolates appearing in Taiwan before other global locations were also observed and recorded in addition to a gradual increase in the number of N-linked glycosylation sites on partial HA1 proteins since 1980. The results from pairwise comparisons of HA1 nucleotide and deduced amino acid sequences indicate shared identities within groups organized according to their bootstrap and P-values of approximately 95.5-100% and 95.7-100% in H1N1 and 94.5-100% and 93.2-100% in H3N2 viruses, respectively. Comparisons of amino acid substitutions in the five antigenic regions reveal highly non-synonymous changes occurring in the Sb region of H1N1 and in the B region of H3N2. The results of an antigenic analysis using a hemagglutinin inhibition (HI) test indicate the presence of some epidemic strains 1-2 years earlier in Taiwan than in other parts of the world, as well as higher vaccine mismatch rates. This information supports the need for continuous surveillance of emerging influenza viruses in Taiwan, which will be useful for making global vaccine decisions. J. Med. Virol. 81:1457 -1470 , 2009 . ß 2009 INTRODUCTION Influenza viruses cause symptomatic diseases in about 20% of children and 5% of adults every year [Glezen et al., 1997; Redlberger et al., 2007] . Influenza viruses belong to the Orthomyxoviridae family of RNA viruses, which is characterized by minus-stranded RNA and a segmented genome. Influenza viruses are divided into subtypes according to the antigenic properties of surface glycoprotein, hemagglutinin (HA), and neuraminidase (NA) [Webster et al., 1992; Saito, 2006 ]. Influenza viruses have caused epidemics and occasional pandemics. The 1918 Spanish, 1957 -1958 Asia, and 1968 -1969 Hong Kong influenza strains were at the center of three 20th century pandemics [Silverstein, 2006] . Any new influenza A virus subtype poses a threat to human health with significant morbidity and mortality. Changes in the antigenicity of influenza virus are the results of antigenic drifts and shifts. Antigenic drifts are caused by the low fidelity of viral RNA polymerase and the accumulation of point mutations in the HA and NA genes. Antigenic shifts result from the reassortment of gene segments between human, avian, and swine influenza viruses [Nicholson et al., 2003; Pontoriero et al., 2003; Carrat and Flahault, 2007] . The H1N1 and H3N2 influenza viruses have been co-circulating worldwide since 1977. Due to the high mutation frequency that is characteristic of the influenza virus genome, the continuous monitoring of viral genetic and antigenic changes is considered essential to public health. Influenza virus HA protein is the primary target of neutralizing antibodies. Influenza viruses can evade host immune systems by mutations on specific antigenic regions that are no longer recognized by neutralizing antibodies. Such mutations are usually generated on the HA1 globular domain, near receptor binding sites [Underwood et al., 1987; Ruigrok et al., 1988] . Results from sequencing analyses of HA1 indicate that amino acid substitutions resulting in antigenic drift primarily accumulate on five antigenic sites in H1N1 (Sa, Sb, Ca1, Ca2, Cb) and H3N2 (A, B, C, D, E) [Gerhard et al., 1981; Wiley et al., 1981] . Previous reports have suggested that more severe epidemics have resulted from antigenic changes in the HA proteins of new antigenic variants [Shih et al., 2005] . Global influenza virus activity is monitored by a surveillance network of laboratories coordinated by the World Health Organization (WHO). Each year WHO makes recommendations for the strain compositions of trivalent influenza vaccines for the northern and southern hemispheres. Suitable vaccine strains are chosen based on antigenic and genetic data for virus strains isolated in previous years [Hsieh et al., 2005; Carrat and Flahault, 2007] . Influenza vaccine efficacy depends on whether currently circulating strains are related antigenically to the vaccine strains. The present study is based on clinical specimens to gather isolates of epidemic H1N1 and H3N2 influenza strains that appeared in Taiwan between 1980 and 2006, and to carry out phylogenetic analyses and antigenic site comparisons to define genetic and antigenic characteristics of the HA gene. Hemagglutinin inhibition (HI) assays were carried out to evaluate antigenicity between the collected isolates and vaccine strains recommended by WHO. Throat swab specimens containing epidemic influenza A strains were collected by physicians and researchers at the Veterans General Hospital-Taipei between 1980 and 2006. Virus subtypes were determined by immunofluorescent antibody (IFA) staining and real-time RT-PCR. Genetic and phylogenetic analyses were carried out on 28 H3N2 and 35 H1N1 influenza A isolates. All virus isolate data were confirmed by the US Centers for Disease Control (US CDC). Viral RNA was extracted from 140 ml of virus-infected MDCK culture fluid with a QIAmp Viral RNA Mini Kit (Qiagen, Hilden, Germany). Reverse transcription was performed using 1 mM of Unil2 primer (5 0 -AGCAAAAGCAGG-3 0 ) with SuperScript II reverse transcriptase (Gibco-BRL, Burlington, ON, Canada) as described by Chan et al. [2006] and Hoffmann et al. [2001] . Partial HA1 segments containing major antigenic sites of H1N1 and H3N2 viruses were amplified by PCR using Taq polymerase with the following primers: HA-H1F, 5 0 -GATGCAGACACAA TATGTAGAGG-3 0 ; HA-H1R, 5 0 -CNCTACAGAGACA TAAGGATTT-3 0 ; HA-H3F, 5 0 -CTGTTACCCTTATGATGTGC-3 0 ; and HA-H3R, 5 0 -GCTYCCATTTGGAGTRATGC-3. Cycling consisted of 958C for 3 min (1 cycle); 948C for 1 min, 558C for 1 min, and 728C for 2 min (35 cycles); and 728C for 7 min (1 cycle). Resulting amplicons were analyzed by 1% gel electrophoresis and purified using a QIAquick Gel Extraction Kit (Qiagen, Valencia, CA). DNA sequencing was performed with a Big Dye Terminator V.3.0 Cycle Sequencing Ready Reaction Kit (ABI, Foster City, CA) and ABI-310 Genetic Analyzer (PerkinElmer, Norwalk, CT). Sequences were analyzed using Chromas Lite version 2.01 (Technelysium Pty. Ltd., Helensvale, Australia) and Bioedit Sequence Alignment Editor V.7.0.5.3. All nucleotide sequences were submitted to GenBank. HI tests were carried out using whole viruses and panels of post-infection ferret sera with representative vaccine strains obtained from the US CDC. Some results were confirmed by the CDC. Any isolate reacting with a reference antiserum at a titer equal to that at which it reacted with the reference homologous virus (AE1-to 2-fold dilution) was said to resemble closely that virus [Kendal et al., 1981; Cox et al., 1983] . This criterion was used to analyze all matches between circulating and vaccine strains [Pontoriero et al., 2003; Hsieh et al., 2005] . Relationships between individual isolates and vaccine strains of the H1N1 and H3N2 influenza A viruses for each year were elucidated by analyzing their amino acid signature patterns. Five major antigenic sites were identified in the HA1 domain of H1N1 (Sa, Sb, Ca1, Ca2, Cb) and five in the same domain of H3N2 (A, B, C, D, E). Deduced amino acid sequences were aligned using BioEdit Sequence Alignment Editor V.7.0.5.3. The amino acid sequences of A/PR/8/34 and A/HK/1/68 were used as references. Phylogenetic trees were constructed using a neighborjoining (NJ) method based on Kimura's two-parameter distance matrix with 1,000 bootstrap replicates, using the MEGA (version 3.0) and PHYLIP (version 3.5) software packages [Kumar et al., 2004; Lin et al., 2007] . MEGA was used to construct NJ trees and PHYLIP was used with parsimony and maximum-likelihood (ML) methods to verify the topology of taxas shown in the trees. All trees presented in this article are NJ; labels on tree branching points indicate ML-calculated P-values. Predictions of N-linked glycosylation on H1N1 or H3N2 HA proteins were made using the NetNGlyc 1.0 Server (http://www.cbs.dtu.dk/services/NetNGlyc/). Amino acid sequences from yearly isolates were uploaded to the website and total N-linked glycosylation sites were predicted according to the Asn-X-Ser/Thr rule. For partial HA1 segments, sites were predicted within a range of a.a 103-280 for H3N2 and a.a 12-197 for H1N1. Subtypes of 63 human influenza A isolates determined by IFA staining and real-time RT-PCR are shown in Table I . Nucleotide sequences of partial HA segments were used to construct phylogenetic trees in order to characterize genetic variations and to determine relationship between vaccine strains and the collected human H1N1/H3N2 influenza viruses. H3N2 phylogenetic analysis results reveal occurrences of antigenic drift every 1-2 years. The tree topologies indicate continuous change patterns. The resulting phylogenetic tree can be divided into 12 groups according to bootstrap values (NJ method) and P-values (ML method) (Fig. 1A ). In contrast, the H1N1 phylogenetic analysis shows antigenic drift occurring every 1-4 years, with some early isolated strains reemerging more than once (Fig. 1B) . For example, the 2005-2006 isolates TW030110-06, TW120202-05, TW031612-06, and TW021011-06 clustered with some 2004 and 2005 isolates and had individual sections that were closer in structure to 1999 Taiwan isolates. The H1N1 phylogenetic tree can be divided into 11 groups and two lineages originating in 1992, with one lineage evolving into 1995-1996 isolates and the other into 1998-2006 isolates (Fig. 1B) . Note that the H1N1 and H3N2 strains recommended by WHO for constructing vaccinations clustered frequently with earlier Taiwanese isolates (Fig. 1A,B) . HA1 nucleotide and deduced amino acid sequences were compared within and between groups of isolates and other Taiwanese reference strains. Each group was identified according to its bootstrap and P-values. The results indicate that nucleotide and amino acid identities within each group were 95.5-100% and 95.7-100% in H1N1 and 94.5-100% and 93.2-100% in H3N2, respectively. The identities between each group were 88.5-100% and 84.2-100% in H1N1 and 84.2-100% and 79.7-100% in H3N2, respectively (Tables IIA and IIB ). Note that amino acid identities were sometimes lower than nucleotide identities within or between the two groups. Hemagglutinin proteins constitute receptor-binding and membrane fusion glycoproteins in influenza viruses. The antigenic residues of H1N1 and H3N2 (after aligning deduced amino acid sequences containing major targets for neutralizing antibodies) are shown, respectively, in Figure 2A ,B. The greatest variation was noted in Sb (H1N1) and B (H3N2) antigenic sites ( Fig. 2A,B) . Results from a comparison of HA1 antigenic sites for H1N1 indicate that most 1983-1997 isolates contained Ser(S)73, Phe(F)75, Glu(E)156, and Ile(I)194, and that these changed to Leu(L)73, Ile(I)75, Gly(G)156, and Leu(L)194 in the 1998-2006 isolates. Furthermore, sites 189, 192, 193, 196, and 197 (Fig. 2B) . HI assays using reference antisera were carried out to identify the antigenic characteristics of H1N1 and H3N2 isolates. As shown in Table III, the H3H2 results show similarities between the 1980-1981, 1984-1985, 1987, 1988, 1990, 1992, 2002, 2004, 1983, 1985, 1988-1990, 1991-1992, 1995-1996, the 2002-2006 isolates and A/Brazil/11/78, A/Chile/1/83, A/Singapore/ 6/86, A/Texas/36/91, A/Bayern/7/95, A/New Caledonia/ 20/99, and A/Solomon Island/3/06 strains. The data in Table III also show that antigenic relationships between isolates and WHO-recommended vaccine strains were stronger in H1N1 (8/11; 72.7%) than in H3N2 (6/12; 50%), and that N-linked glycosylation sites on partial HA1 segments increased gradually from 1-2 in 1980 to 6-7 in 2006. After combining these results with existing data from the Clinical Virology Laboratory of the Taipei Bootstrap and P-values were used to divide H3N2 influenza virus isolates into 12 groups (G1-G12) and H1N1 isolates into 11 groups (G1-G11). Changes in amino acids on antigenic sites in different periods are labeled, with comparisons of earlier and later periods. Minor amino acid substitutions are displayed within brackets. Squares indicate vaccine strains, and triangles indicate our isolates. *P < 0.05; **P < 0.01. Shih et al. [2005] , comparisons were made with lists of vaccine strains recommended by WHO, world epidemic strains, and Taiwanese isolates. The vaccine match rate for H3N2 (19/28, 67.8%) was lower than for H1N1 (23/28, 82.1%). The data also indicate that some world epidemic strains appeared in Taiwan before appearing in other parts of the world (Table IV) . The Taiwan Centers for Disease Control (CDC) were founded in 1999, and the surveillance network coordinated by the Taiwan CDC was established the following year [Shih et al., 2005; Lin et al., 2008] . The Veterans General Hospital-Taipei is a major teaching hospital and medical center. Prior to 2000, the hospital's Department of Virology served as a national reference laboratory for influenza virus isolation and characterization. After 2000, the hospital's specimen collection activity was restricted to the northern region of Taiwan [Shih et al., 2005] . In response to this limitation, influenza A viruses collected and/or examined by other Taiwanese laboratories were included in this analysis. As a result, the viral sequences described in this study represent influenza virus A strains prevalent in Taiwan between 1980 and 2006. The three most commonly used algorithms in molecular phylogenetic analyses are NJ, ML, and Parsimony (Pars); all three were used to obtain the tree topologies discussed in this article. The phylogenetic tree can be segregated into 12 (H3N2) and 11 (H1N1) groups according to bootstrap values obtained from the NJ J. Med. Virol. DOI 10.1002/jmv 4 -83.7 86.5 -88.4 81.4-83.7 86.3-88.4 81.4-84.2 G3 (1986-1992) 94.5 -99.2 93.2-99.4 91.7-94.3 87.6 -90.4 90.8 -93.6 85.3-89.3 90.4-93.0 84.8 -88.7 89.5 -92.7 83.1-87.6 89.7-92.7 83.7 -87.6 88.5 -90.8 82.5-84.8 87.4-89.7 82.0 -84.2 87.4 -89.7 82.0-84.2 87.1-90.0 80.8-84.8 G4 (1995-1997) method and P-values from the ML method. In terms of robustness, a bootstrap value of 70% is cited frequently as a cutoff for cluster reliability [Hillis, 1997 [Hillis, , 1998 Lan et al., 2005] . When using the ML method, phylogenetic tree topology is considered significant when P < 0.05. The results of the H3N2 phylogenetic analysis reveal a single evolutionary lineage and a pattern of continuous nucleotide change, identical to recent epidemic and pandemic histories of influenza A viruses resulting from continuous H3N2 antigenic drift [Abed et al., 2002; Redlberger et al., 2007] . Previous researchers have reported H3 HA gene evolution as a single lineage following the introduction of H3N2 viruses into human populations in 1968; H1N1 viruses in particular represent distinct co-circulating lineages [Hay et al., 2001; Bragstad et al., 2008] . However, the previous data indicate the disappearance and reemergence of H1N1 circulating initially from 1918 to 1957 and reappearing in 1977, with the new strains sharing considerable similarities with those circulating in the 1950s [Hsieh Fig. 2 . Amino acid substitution patterns on major antigenic sites of the human influenza A virus in Taiwan from 1980 to 2006. A: Synonymous and non-synonymous amino acid changes on antigenic sites Sa, Sb, Ca1, Ca2, and Cb of the HA1 gene of the H1N1 influenza virus, corresponding to residues 12-197 of A/PR/8/34. Periods (labeled with bold dashed lines) were identified based on amino acid signature patterns and phylogenetic tree results. Some Taiwanese reference strains (non-bold dashed lines) lack information for partial Sb regions. B: Synonymous and non-synonymous amino acid changes on antigenic sites A, B, C, D, and E of the HA1 gene of the H3N2 influenza virus, corresponding to residues 103-280 of A/HongKong/1/86. Periods (labeled with bold dashed lines) were identified based on amino acid signature patterns and phylogenetic tree results. et al., 2005] . The present results also show two separate H1N1 lineages originating in 1992, supporting the findings of Shih et al. [2005] that H1N1 isolates in Taiwan from 1995 to 2003 can be divided into two major groups. The results from the phylogenetic analyses and antigenic site comparisons indicate higher evolutionary rates in the five H3N2 antigenic regions compared to their H1N1 counterparts. An average of 1-3 amino acid changes at H3N2 and H1N1 antigenic sites were found over 1-to 2-and 1-to 4-year intervals, respectively. According to Smith et al. [2002] , new drift variants of epidemiologic importance generally have four or more amino acid substitutions located in two or more antigenic sites on HA molecules. The present results also indicate that the identities of the amino acid sequences in the HA1 genes of Taiwanese isolates were at times lower than those of the nucleotide sequences. These findings indicate the dominance of non-synonymous change in the antigenic sites of human influenza A viruses, especially in the H3N2 subtype. Jian et al. [2008] reported significantly more non-synonymous than synonymous nucleotide substitutions at these antigenic sites. Plotkin and Dushoff [2003] , who studied influenza A gene sequence evolution for the past two decades, state that HA codons, especially those in antibody-binding regions, express a significant bias toward substitution mutations compared to other internal viral genes. The results from Hsieh et al. [2005] and Tseng et al. [1996] also suggest that among Taiwan epidemic strains, H3N2 shows the highest antigenic variation, followed by type B and H1N1. These findings indicate that non-synonymous changes in the HA1 gene, especially in antigenic regions, have been more frequent in H3N2 than in H1N1 over the past 25 years in Taiwan, perhaps triggering outbreaks and/or epidemics in this country. HI assays are used most commonly for measuring antigenic differences in epidemic strains, an important aspect of influenza virus surveillance [Sheth and Jhala, 1976; Ikic et al., 1977] . The HI test results of the present study indicate that epidemic strains and vaccine strains J. Med. Virol. DOI 10.1002/jmv Fig. 2. (Continued) were less similar in terms of antigenicity, and that vaccine-like strains often circulated in Taiwan 1-2 years before other global locations (Tables III and IV) . The second finding was confirmed by phylogenetic analyses and comparisons of amino acid substitution patterns at antigenic sites, that is, the Taiwanese epidemic strains not only clustered with later vaccine strains, but also displayed similar amino acid signature patterns. In addition, Hsieh et al. [2005] described the initial appearance of strains in Taiwan, including the A/ Beijing/262/95(H1N1)-like strain (used for the 1998-1999 vaccine) that was found in France, Senegal, South Africa, and the United States; this strain was isolated in Taiwan in 1996 . Shih et al. [2005 also reported on the isolation of the A/Fujian/411/2002(H3N2)-like strain in Taiwan in [2003] [2004] 1 year prior to its proposed use as a vaccine strain. In the present study, a comparison of WHO-recommended and Taiwanese epidemic strains for H1N1 and H3N2 produced matching rates of 82.1% and 67.8%, respectively. Compare these with H1N1 vaccine match rates from 1987/1988 to 2003/2004 of 82% for Taiwan and 98.6% for Canada, and with H3N2 vaccine match rates of 53% for Taiwan and 64.9% for Canada [Hsieh et al., 2005 ]. Other researchers have described a 1-to 2-year time lag between WHO-recommended vaccine strains for northern hemisphere nations and circulating strains in Taiwan and western countries [Hsieh et al., 2005; Shih et al., 2005] . Mismatches between Taiwanese epidemic and vaccine strains may be due to the unique location of Taiwan on the boundary between north and south Asia and its proximity to China, which is considered the epicenter of influenza viruses. New antigenic drift variants may spread into Taiwan by cross-border travel or migrating birds. A second possible explanation is the 6-to 9-month lag time between the selection of vaccine strain and manufacture, a phenomenon also observed in European and North American countries [de Jong et al., 2000; Carrat and Flahault, 2007; Skowronski et al., 2007] . Compared to the results from phylogenetic analyses and antigenic subtypes, the data described in this article suggest a strong association between HA genetic variation and antigenic changes in influenza isolates used to link epidemic and vaccine strains. Due to the limited availability of HI testing equipment in clinical virology laboratories in Taiwan, there may be a greater reliance at present on HA sequencing for evaluating and A/India/6263/80 a comparing differences between local epidemic isolates and vaccine strains. Surface protein glycosylation plays an important role in the biology of many viruses, including the severe acute respiratory syndrome coronavirus (SARS-CoV) ], hepatitis C virus [Goffard et al., 2005] , West Nile virus [Hanna et al., 2005] , HIV-1 [Bashirova et al., 2001] , and influenza viruses [Roberts et al., 1993] . Surface glycosylation may help these viruses enter target cells and evade host immune responses [Fournillier et al., 2001; Abe et al., 2004] . HIV-1 and influenza A viruses have been described as using glycosylated gp120 and HA proteins to interact with and mediate entry via DC-SIGN and mannose molecules on dendritic cells [Reading et al., 2000; Bashirova et al., 2001; Wang et al., 2008] . However, N-linked glycosylation on the HA proteins of influenza viruses has been reported as serving four important biological functions: (a) adding carbohydrates to antigenic sites, which interferes with neutralizing antibody access [Skehel et al., 1984] ; (b) influencing virus pathogenicity when it occurs near cleavage sites [Ohuchi et al., 1989] ; (c) sialylating oligosaccharide chains near receptor-binding sites, which interferes with viral binding [Ohuchi et al., 1995] ; and (d) promoting the proper folding and maintenance of protein conformation or stability [Roberts et al., 1993] . According to the present results, N-linked glycosylation sites on partial HA1 proteins in both H1N1 and H3N2 have gradually increased over time. When describing this same phenomenon, Vigerust et al. [2007] and Zhang et al. [2004] also described human influenza A HA H3 and H1 HA1 accumulating sequons over time, and a gradual increase in N-linked glycosylation in and around the globular head of H3N2 viruses during the past four decades. Vigerust et al. [2007] suggest that the increased glycosylation of influenza viruses may result in decreased virulence. The ongoing surveillance of influenza virus activity is essential for preventing human pandemics or severe epidemics. With this goal in mind, the Taiwan Centers for Disease Control and Prevention established a national influenza surveillance network in 2000 [Shih et al., 2005; Jian et al., 2008] . However, its work has been hindered by limited information on epidemic strains prior to 2000. This article described the genetic and antigenic characteristics of the HA gene in influenza A isolates from 1980 to 2006 and addresses Taiwan's vaccine mismatch phenomenon. Other gene segments of these isolates are being studied. A/New Caledonia/20/99 a A/New Caledonia/20/99 a A/Panama/2007/99 a,b (Continued) Mismatch between WHO-recommended vaccine and epidemic strains in Taiwan. b Epidemic strain circulated in Taiwan earlier than in other countries and prior to WHO vaccine recommendation Effect of the addition of oligosaccharides on the biological activities and antigenicity of influenza A/H3N2 virus hemagglutinin Divergent evolution of hemagglutinin and neuraminidase genes in recent influenza A:H3N2 viruses isolated in Canada A dendritic cell-specific intercellular adhesion molecule 3-grabbing nonintegrin (DC-SIGN)-related protein is highly expressed on human liver sinusoidal endothelial cells and promotes HIV-1 infection The evolution of human influenza A viruses from 1999 to 2006: A complete genome study Influenza vaccine: The challenge of antigenic drift Amplification of the entire genome of influenza A virus H1N1 and H3N2 subtypes by reverse-transcription polymerase chain reaction Laboratory-based surveillance of influenza A(H1N1) and A(H3N2) viruses in 1980-81: Antigenic and genomic analyses Mismatch between the 1997/1998 influenza vaccine and the major epidemic A(H3N2) virus strain as the cause of an inadequate vaccine-induced antibody response to this strain in the elderly Induction of hepatitis C virus E1 envelope protein-specific immune response can be enhanced by mutation of N-glycosylation sites Antigenic structure of influenza virus haemagglutinin defined by hybridoma antibodies Influenza virus infections in infants Role of N-linked glycans in the functions of hepatitis C virus envelope glycoproteins N-linked glycosylation of west nile virus envelope proteins influences particle assembly and infectivity The evolution of human influenza viruses Phylogenetic analysis Taxonomic sampling, phylogenetic accuracy, and investigator bias Universal primer set for the full-length amplification of all influenza A viruses Influenza in Taiwan: Seasonality and vaccine strain match HI antibody response in volunteers vaccinated with live influenza A/New Jersey/76 vaccine Genetic and epidemiological analysis of influenza virus epidemics in Taiwan during Effect of test system on the ability of monoclonal antibodies to detect antigenic drift in influenza A(H1N1) virus haemagglutinins MEGA3: Integrated software for molecular evolutionary genetics analysis and sequence alignment Molecular epidemiology of severe acute respiratory syndrome-associated coronavirus infections in Taiwan Molecular epidemiology of HIV-1 infection and full-length genomic analysis of circulating recombinant form 07_BC strains from injection drug users in Taiwan Genetic and antigenic analysis of epidemic influenza viruses isolated during 2006-2007 season in Taiwan Mutations at the cleavage site of the hemagglutinin after the pathogenicity of influenza virus A/chick/Penn/83 (H5N2) Neuraminidase is essential for fowl plague virus hemagglutinin to show hemagglutinating activity Codon bias and frequency-dependent selection on the hemagglutinin epitopes of influenza A virus Antigenic and genomic relation between human influenza viruses that circulated in Argentina in the period 1995-1999 and the corresponding vaccine components Involvement of the mannose receptor in infection of macrophages by influenza virus Dynamics of antigenic and genetic changes in the hemagglutinins of influenza A/H3N2 viruses of three consecutive seasons Role of conserved glycosylation sites in maturation and transport of influenza A virus hemagglutinin Studies on the structure of the influenza virus haemagglutinin at the pH of membrane fusion Review on replication cycle of influenza virus Seropatterns of influenza HI antibodies Laboratory-based surveillance and molecular epidemiology of influenza virus in Taiwan Identifying epitopes responsible for neutralizing antibody and DC-SIGN binding on the spike glycoprotein of the severe acute respiratory syndrome coronavirus Preparing for pandemic influenza A carbohydrate side chain on hemagglutinins of Hong Kong influenza viruses inhibits recognition by a monoclonal antibody Estimating vaccine effectiveness against laboratory-confirmed influenza using a sentinel physician network: Results from the 2005-2006 season of dual A and B vaccine mismatch in Canada Molecular epidemiology of influenza A(H3N2) virus reinfections Influenza virus infections in Taiwan from 1979 to 1995 Receptor-binding characteristics of monoclonal antibody-selected antigenic variants of influenza virus N-linked glycosylation attenuates H3N2 influenza viruses DC-SIGN mediates avian H5N1 influenza virus infection in cis and in trans Evolution and ecology of influenza A viruses Structural identification of the antibody-binding sites of Hong Kong influenza haemagglutinin and their involvement in antigenic variation Tracking global patterns of N-linked glycosylation site variation in highly variable viral glycoproteins: HIV, SIV, and HCV envelopes and influenza hemagglutinin We thank our colleagues from the AIDS Prevention and Research Center of National Yang-Ming University and the Division of Clinical Virology, Department of Pathology and Laboratory Medicine of Taipei Veterans General Hospital for their helpful feedback, Ms. Diana Kou and Mr. Jon Lindemann for editing assistance.