Genetic analysis and evaluation of the reassortment of influenza B viruses isolated in Taiwan during the 2004-2005 and 2006-2007 epidemics

Host-range shift of H3N8 canine influenza virus: a phylodynamic analysis of its origin and adaptation from equine to canine host

Prior to the emergence of H3N8 canine influenza virus (CIV) and the latest avian-origin H3N2 CIV, there was no evidence of a circulating canine-specific influenza virus. Molecular and epidemiological evidence suggest that H3N8 CIV emerged from H3N8 equine influenza virus (EIV). This host-range shift of EIV from equine to canine hosts and its subsequent establishment as an enzootic CIV is unique because this host-range shift was from one mammalian host to another. To further understand this host-range shift, we conducted a comprehensive phylodynamic analysis using all the available whole-genome sequences of H3N8 CIV. We found that (1) the emergence of H3N8 CIV from H3N8 EIV occurred in approximately 2002; (2) this interspecies transmission was by a reassortant virus of the circulating Florida-1 clade H3N8 EIV; (3) once in the canine species, H3N8 CIV spread efficiently and remained an enzootic virus; (4) H3N8 CIV evolved and diverged into multiple clades or sublineages, with intra and inter-lineage reassortment. Our results provide a framework to understand the molecular basis of host-range shifts of influenza viruses and that dogs are potential “mixing vessels” for the establishment of novel influenza viruses.

Seasonal dynamics of influenza viruses and age distribution of infected individuals across nine seasons covering 2009-2018 in Taiwan

BACKGROUND/PURPOSE

A swine-origin influenza A/H1N1 virus (termed A/H1N1pdm) caused a pandemic in 2009 and has continuously circulated in the human population. To investigate its possible ecological effects on circulating influenza strains, the seasonal patterns of influenza viruses and the respective age distribution of infected patients were studies.

METHODS

The data obtained from national influenza surveillance systems in Taiwan from July 2009 to June 2018 were analyzed.

RESULTS

The A/H1N1pdm and A/H3N2 strains usually caused a higher ratio of severe to mild cases than influenza B. New variants of A/H1N1pdm and A/H3N2 emerged accompanied by a large epidemic peak. However, the new influenza B variants intended to circulate for several seasons before causing a large epidemic. The major group of outpatients affected by A/H1N1pdm were aged 13-23 years in the pandemic wave, and the age range of infected individuals gradually shifted to 24-49 and 0-6 years across seasons; A/H1N1pdm-infected inpatients were aged 24-49 years in 2009-2011, and the age range gradually switched to older groups aged 50-65 and >65 years. Individuals aged 0-6 or 24-49 years accounted for the majority of A/H3N2-infected outpatients across seasons, whereas most of the inpatients affected by A/H3N2 were aged >65 years.

CONCLUSION

Understanding the effects of new variants and changes in dominant circulating viral strains on the age distribution of the affected human population, disease severity and epidemic levels is useful for the establishment of fine-tuned strategies for further improvement of influenza control.

Unreliable usage of a single influenza virus IgM antibody assay in influenza-like illness: A retrospective study of the 2016-2018 flu epidemic

We retrospectively analyzed serum IgM antibodies (Abs) to influenza viruses from two tertiary hospitals in Beijing from December 2016 to February 2018. Samples from 36,792 patients, aged 0–98 years, were collected and tested. Among the patients, 923 children from two winter flu seasons were assayed with both antigens and IgM Abs to Flu A and Flu B and assigned as paired groups. Another 2,340 adults and 1,978 children with only antigen tested in the 2016 and 2017 winter flu seasons were named as unpaired groups. IgM Abs-positivity rates in children were 0.80% and 36.57% for Flu A and Flu B, respectively, peaking at 4–5 years of age. For adults, the Flu A and Flu B IgM Abs-positivity rates were 10.34% and 21.49%, respectively, peaking at 18–35 years of age. The trend of temporal distribution between the children and the adults was significantly correlated for IgM Abs to Flu B, but not for Flu A. Compared with unpaired groups, the detection rate of Flu A antigen was significantly higher than IgM Abs in children, whereas frequencies of IgM Abs were higher than antigen in adults. Incidence of Flu B antigen was sharply increased in 2017 winter than in the 2016 winter in both children and adults, but no concomitant increase was observed in IgM Abs to Flu B. For paired children groups, incidence of Flu B antigen in the 2017 flu season was significantly higher than that in the 2016 flu season; in contrast, positive rates of IgM Abs in the 2017 flu season were even lower than those in 2016. Considering antigen detection may reflect the Flu A/Flu B epidemic, our results indicate single-assayed IgM Abs were less effective in the diagnosis of acute influenza virus infection, and the use of this assay for epidemiology evaluations was not supported by these findings.

Broadened immunity against influenza by vaccination with computationally designed influenza virus N1 neuraminidase constructs

Split inactivated influenza vaccines remain one of the primary preventative strategies against severe influenza disease in the population. However, current vaccines are only effective against a limited number of matched strains. The need for broadly protective vaccines is acute due to the high mutational rate of influenza viruses and multiple strain variants in circulation at any one time. The neuraminidase (NA) glycoprotein expressed on the influenza virion surface has recently regained recognition as a valuable vaccine candidate. We sought to broaden the protection provided by NA within the N1 subtype by computationally engineering consensus NA sequences. Three NA antigens (NA5200, NA7900, NA9100) were designed based on sequence clusters encompassing three major groupings of NA sequence space; (i) H1N1 2009 pandemic and Swine H1N1, (ii) historical seasonal H1N1 and (iii) H1N1 viruses ranging from 1933 till current times. Recombinant NA proteins were produced as a vaccine and used in a mouse challenge model. The design of the protein dictated the protection provided against the challenge strains. NA5200 protected against H1N1 pdm09, a Swine isolate from 1998 and NIBRG-14 (H5N1). NA7900 protected against all seasonal H1N1 viruses tested, and NA9100 showed the broadest range of protection covering all N1 viruses tested. By passive transfer studies and serological assays, the protection provided by the cluster-based consensus (CBC) designs correlated to antibodies capable of mediating NA inhibition. Importantly, sera raised to the consensus NAs displayed a broader pattern of reactivity and protection than naturally occurring NAs, potentially supporting a predictive approach to antigen design.

The high variability of the influenza virus — arising from its high mutation rate and wide range of strains — limits the effectiveness of influenza vaccines unless they induce a broad immune response, a difficult task when relying on natural viral antigens. Here, Xavier Saelens, Thorsten Vogel, Ray Oomen and colleagues applied a ‘cluster-based’ consensus computational approach to design three consensus sequences of the viral protein neuroaminidase (NA) subtype 1 that induce broadly protective immune responses in vaccinated mice. NA9100, a consensus NA sequence based on H1N1 virus strains collected from 1933 to today, was protective against all N1 viruses tested. By using a computational method to integrate multiple sequences of viral proteins into one consensus protein, the researchers provide a strategy that can be applied to develop broadly protective vaccine formulations for influenza virus.

Significant circulation of influenza B viruses mismatching the recommended vaccine-lineage in South Korea, 2007-2014

We aimed to characterize the lineages of influenza B viruses obtained from clinical specimens during the 2007-2014 seasons in South Korea. RT-PCR for the partial hemagglutinin gene of influenza B virus was performed on laboratory-confirmed influenza B samples from the 2007-2008 season to 2013-2014 season. A phylogenetic tree was generated, and current influenza vaccine strains for the Northern Hemisphere were used as representative strains of Victoria and Yamagata lineages. A total of 571 influenza B virus sequences were analyzed. During the 2009-2010 season, most of the circulating influenza B viruses matched the vaccine strain; 91.0% (91/100) of viruses belonged to the Victoria lineage. In the 2007-2008, 2011-2012, and 2013-2014 seasons, co-circulation of each influenza B lineage was found with a match ratio to the vaccine strain of 53.2% (42/79), 40.9% (63/154), and 58.3% (134/230), respectively. Overall, 41.7% (238/571) of the circulating influenza B viruses belonged to the lineage mismatching the vaccine strain. During the seven influenza seasons, influenza B epidemics were substantial in four seasons in South Korea. Significant mismatches of the vaccine and lineage of the circulating influenza B viruses were found. The current trivalent influenza vaccine may not be fully suitable for effective protection against influenza B.

Literature review of the epidemiology of influenza B disease in 15 countries in the Asia-Pacific region

Influenza control strategies focus on the use of trivalent influenza vaccines containing two influenza A virus subtypes and one of the two circulating influenza type B lineages (Yamagata or Victoria). Mismatches between the vaccine B lineage and the circulating lineage have been regularly documented in many countries, including those in the Asia‐Pacific region. We conducted a literature review with the aim of understanding the relative circulation of influenza B viruses in Asia‐Pacific countries. PubMed and Western Pacific Region Index Medicus were searched for relevant articles on influenza type B published since 1990 in English language for 15 Asia‐Pacific countries. Gray literature was also accessed. From 4834 articles identified, 121 full‐text articles were analyzed. Influenza was reported as an important cause of morbidity in the Asia‐Pacific region, affecting all age groups. In all 15 countries, influenza B was identified and associated with between 0% and 92% of laboratory‐confirmed influenza cases in any one season/year. Influenza type B appeared to cause more illness in children aged between 1 and 10 years than in other age groups. Epidemiological data for the two circulating influenza type B lineages remain limited in several countries in the Asia‐Pacific, although the co‐circulation of both lineages was seen in countries where strain surveillance data were available. Mismatches between circulating B lineages and vaccine strains were observed in all countries with available data. The data suggest that a shift from trivalent to quadrivalent seasonal influenza vaccines could provide additional benefits by providing broader protection.

Drug-Resistant and Genetic Evolutionary Analysis of Influenza Virus from Patients During the 2013 and 2014 Influenza Season in Beijing

The study aimed to analyze drug resistance and mutations and genetic evolution of influenza A and influenza B viruses during the 2013 and 2014 influenza season in Beijing, China. RNA was extracted from pharyngeal or nasal swabs of 28 patients, and determination of influenza genotypes was performed by using real-time reverse-transcription polymerase chain reaction. Influenza A virus samples were sequenced with the neuraminidase (NA) gene and M2 matrix protein gene to determine the NA inhibitor (NAI) resistance and amantadine resistance mutations, and influenza B virus samples were sequenced with the NA gene and hemagglutinin (HA) gene to analyze NAI resistance mutations. As a result, the enrolled subjects consisted of 19 patients with the A(H1N1)pdm09 subtype, four with A(H3N2) subtype and five with influenza B virus. All of the 23 samples with influenza A viruses harbored amantadine resistance mutation S31N in M2 matrix protein. V241I, a compensatory NAI resistance mutation, was detected in all of the 19 A(H1N1)pdm09 viruses. No other NAI resistance mutation was observed in both influenza A and B viruses. The NA gene of the five influenza B virus strains was classified as B-Victoria lineage, while the HA gene of five strains was classified as B-Yamagata lineage. In summary, all influenza A viruses from patients in Beijing in the 2013-2014 season were resistant to amantadine agent. Both influenza A and B viruses kept sensitive to NAIs. Lineage recombination was detected in influenza B virus strains and may impair the efficacy of influenza vaccination.

Circulating pattern and genomic characteristics of influenza B viruses in Taiwan from 2003 to 2014

BACKGROUND/PURPOSE

Influenza B viruses are antigenically classified into Yamagata and Victoria lineages according to their hemagglutinin (HA) proteins. These two lineages are known to either appear sequentially or cocirculate in Taiwan.

METHODS

Taiwanese influenza B viral HA and neuraminidase (NA) sequences between 2003 and 2014 were determined and analyzed. A time-scaled phylogenetic tree was constructed to decipher the evolutionary trends of these sequences, and the reassortment between the two lineages. Positively selected amino acids were predicted, demonstrating the adaptive mutations of the circulating pattern.

RESULTS

The HA phylogenetic tree revealed that the Victoria lineage evolved into a ladder-like pattern, whereas the Yamagata lineage exhibited complex topology with several independently evolved clades on which viruses from different influenza seasons interlaced. For several seasons, HA sequences were found to be dominated by strains of the same lineage as the corresponding vaccine strain. Inspecting these sequences revealed that frequent mutations occurred in neutralizing epitopes and glycosylation sites. Amino acid positions 212 and 214 of N-glycosylation sites, which are known to be critical determinants of receptor-binding specificity, were found to be subject to positive selection. No drug-resistant sites were noticed in the NA sequences. In addition, we identified several cases of NA reassortment with an overall incidence rate of 6% for the investigated Taiwan strains.

CONCLUSION

We highlighted the interplay between mutations in the glycosylation sites and epitope during HA evolution. These are crucial molecular signatures to be monitored for influenza B epidemics in the future.

Patterns of influenza B circulation in Brazil and its relevance to seasonal vaccine composition

Data on the burden of disease and circulation patterns of influenza B lineages for Brazil are limited. This review aims to describe the pattern of influenza B occurrence in Brazil to have a better understanding of its epidemiology and its relevance when considering seasonal influenza vaccine composition. A review of the data including analysis of international and local surveillance data as well as information from online search of databases using Medical Subject Headings terms in conjunction with screening of abstracts from scientific events was performed. Based on international epidemiologic surveillance data, moderate levels of influenza B disease (19%; 2006-2014) were observed. Of these nine years, it was possible to compare data from three years (2007, 2008 and 2013) which have information on the circulating influenza B lineage. Co-circulation of influenza B lineages was observed in all these three influenza seasons, of which, during one season, a high degree of mismatch between the vaccine lineage and the predominant circulating lineage (91.4% [2013]) was observed. Local surveillance data reveal a distinct and dynamic distribution of respiratory viruses over the years. Data from published literature and abstracts show that influenza B is a significant cause of disease with an unpredictable circulation pattern and showing trends indicating reemergence of the B/Victoria lineage. The abstracts report notable levels of co-circulation of both influenza B lineages (2000-2013). Mismatch between the Southern hemisphere vaccine and the most prevalent circulating viruses in Brazil were observed in five influenza seasons. The evidence on co-circulation of two influenza B lineages and mismatched seasons in Brazil indicates the benefit of quadrivalent influenza vaccines in conferring broader seasonal influenza protection. Additionally, improving influenza surveillance platforms in Brazil is important for monitoring disease trends and the impact of introducing seasonal influenza vaccination.

Epidemiology of pandemic influenza A/H1N1 virus during 2009-2010 in Taiwan

Outbreak of swine-origin influenza A/H1N1 virus (pdmH1N1) occurred in 2009. Taiwanese authorities implemented nationwide vaccinations with pdmH1N1-specific inactivated vaccine as of November 2009. This study evaluates prevalence, HA phylogenetic relationship, and transmission dynamic of influenza A and B viruses in Taiwan in 2009-2010. Respiratory tract specimens were analyzed for influenza A and B viruses. The pdmH1N1 peaked in November 2009, was predominant from August 2009 to January 2010, then sharply dropped in February 2010. Significant prevalence peaks of influenza B in April-June of 2010 and H3N2 virus in July and August were observed. Highest percentage of pdmH1N1- and H3N2-positive cases appeared among 11-15-year-olds; influenza B-positive cases were dominant among those 6-10 years old. Maximum likelihood phylogenetic trees showed 11 unique clusters of pdmH1N1, seasonal H3N2 influenza A and B viruses, as well as transmission clusters and mixed infections of influenza strains in Taiwan. The 2009 pdmH1N1 virus was predominant in Taiwan from August 2009 to January 2010; seasonal H3N2 influenza A and B viruses exhibited small prevalence peaks after nationwide vaccinations. Phylogenetic evidence indicated transmission clusters and multiple independent clades of co-circulating influenza A and B strains in Taiwan.

Surveillance and vaccine effectiveness of an influenza epidemic predominated by vaccine-mismatched influenza B/Yamagata-lineage viruses in Taiwan, 2011-12 season

Introduction

The 2011−12 trivalent influenza vaccine contains a strain of influenza B/Victoria-lineage viruses. Despite free provision of influenza vaccine among target populations, an epidemic predominated by influenza B/Yamagata-lineage viruses occurred during the 2011−12 season in Taiwan. We characterized this vaccine-mismatched epidemic and estimated influenza vaccine effectiveness (VE).

Methods

Influenza activity was monitored through sentinel viral surveillance, emergency department (ED) and outpatient influenza-like illness (ILI) syndromic surveillance, and case-based surveillance of influenza with complications and deaths. VE against laboratory-confirmed influenza was evaluated through a case-control study on ILI patients enrolled into sentinel viral surveillance. Logistic regression was used to estimate VE adjusted for confounding factors.

Results

During July 2011−June 2012, influenza B accounted for 2,382 (72.5%) of 3,285 influenza-positive respiratory specimens. Of 329 influenza B viral isolates with antigen characterization, 287 (87.2%) were B/Yamagata-lineage viruses. Proportions of ED and outpatient visits being ILI-related increased from November 2011 to January 2012. Of 1,704 confirmed cases of influenza with complications, including 154 (9.0%) deaths, influenza B accounted for 1,034 (60.7%) of the confirmed cases and 103 (66.9%) of the deaths. Reporting rates of confirmed influenza with complications and deaths were 73.5 and 6.6 per 1,000,000, respectively, highest among those aged ≥65 years, 50−64 years, 3−6 years, and 0−2 years. Adjusted VE was −31% (95% CI: −80, 4) against all influenza, 54% (95% CI: 3, 78) against influenza A, and −66% (95% CI: −132, −18) against influenza B.

Conclusions

This influenza epidemic in Taiwan was predominated by B/Yamagata-lineage viruses unprotected by the 2011−12 trivalent vaccine. The morbidity and mortality of this vaccine-mismatched epidemic warrants careful consideration of introducing a quadrivalent influenza vaccine that includes strains of both B lineages.

Phylogenetic and evolutionary history of influenza B viruses, which caused a large epidemic in 2011-2012, Taiwan

The annual recurrence of the influenza epidemic is considered to be primarily associated with immune escape due to changes to the virus. In 2011-2012, the influenza B epidemic in Taiwan was unusually large, and influenza B was predominant for a long time. To investigate the genetic dynamics of influenza B viruses during the 2011-2012 epidemic, we analyzed the sequences of 4,386 influenza B viruses collected in Taiwan from 2004 to 2012. The data provided detailed insight into the flux patterns of multiple genotypes. We found that a re-emergent TW08-I virus, which was the major genotype and had co-circulated with the two others, TW08-II and TW08-III, from 2007 to 2009 in Taiwan, successively overtook TW08-II in March and then underwent a lineage switch in July 2011. This lineage switch was followed by the large epidemic in Taiwan. The whole-genome compositions and phylogenetic relationships of the representative viruses of various genotypes were compared to determine the viral evolutionary histories. We demonstrated that the large influenza B epidemic of 2011-2012 was caused by Yamagata lineage TW08-I viruses that were derived from TW04-II viruses in 2004-2005 through genetic drifts without detectable reassortments. The TW08-I viruses isolated in both 2011-2012 and 2007-2009 were antigenically similar, indicating that an influenza B virus have persisted for 5 years in antigenic stasis before causing a large epidemic. The results suggest that in addition to the emergence of new variants with mutations or reassortments, other factors, including the interference of multi-types or lineages of influenza viruses and the accumulation of susceptible hosts, can also affect the scale and time of an influenza B epidemic.

A new antigenic variant of human influenza A (H3N2) virus isolated from airport and community surveillance in Taiwan in early 2009

A new variant of influenza A H3N2 virus emerged in January 2009 and became the dominant strain in Taiwan in April 2009. The variant was also detected in imported cases from various regions, including East and Southeast Asia and North America, indicating that it has circulated globally. Compared to the 2009-2010 vaccine strain, A/Brisbane/10/2007, the hemagglutinin gene of this variant exhibited five substitutions, E62K, N144K, K158N, K173Q and N189K, which are located in the antigenic sites E, A, B, D and B respectively, and it was antigenically distinct from A/Brisbane/10/2007 with more than eight-fold titer reduction in the hemagglutination inhibition reaction. The A/Perth/16/2009 (H3N2)-like virus recommended by World Health Organization for use in the 2010 southern hemisphere and 2010-2011 northern influenza seasons exhibited the same substitutions like this new variant. In addition to regional or community influenza surveillance, the imported cases or airport fever screening surveillance may be a good resource to monitor the evolution of the virus and benefit the real-time information of global influenza circulation.

Probabilistic indoor transmission modeling for influenza (sub)type viruses

OBJECTIVES

To use a probability based transmission modeling approach to examine the influenza risk of infection virus in indoor environments. This was based on 10 years of data gathered from influenza-like illness sentinel physician and laboratory surveillance, and experimental viral shedding data in Taiwan.

METHODS

We integrated sentinel physician-reported cases and positive rates of influenza A (H1N1), A (H3N2), influenza B, and respiratory syncytial virus in Taiwan using the Wells-Riley mathematical model. This model incorporates environmental factors such as room ventilation and breathing rates. We also linked vaccine match rate with related transmission estimations to predict the controllable potential using a control model characterized by basic reproduction number (R(0)) and proportion of asymptomatic infections (theta).

RESULTS

A quantitative framework was developed to better understand the infection risk and R(0) estimates of A (H1N1), A (H3N2), and B viruses. The viral concentration in human fluid was linked successfully with quantum generation rates to estimate virus-specific infection risks. Our results revealed that A (H3N2) virus had a higher transmissibility and uncontrollable potential than the A (H1N1) and B viruses.

CONCLUSIONS

Probabilistic transmission model can incorporate virus-specific data on experimental viral shedding, long-term sentinel physician and laboratory surveillance to predict virus-specific infection risks in Taiwan.

Virological investigation of four outbreaks of influenza B reassortants in the northern region of Taiwan from October 2006 to February 2007

Background

From October 2006 to February 2007, clinical specimens from 452 patients with symptoms related to respiratory tract infection in the northern region of Taiwan were collected. Real-time PCR and direct immunofluorescent antibody tests showed that 145 (32%) patients had influenza B virus infections. Subsequently, nucleotide sequence analyses of both hemagglutinin (HA) and neuraminidase (NA) genes of 39 isolates were performed. Isolated viruses were antigenically characterized using hemagglutinin inhibition (HI) test.

Findings

Phylogenetic tree analysis showed that all the isolates belonged to the B reassortant lineage with HA gene belonged to the B/Victoria/2/87 lineage and the NA gene belonged to the B/Yamagata/16/88 lineage. In addition, a group of children aged between 6 to 8 years old resided in Yilan county were infected with a variant strain. Hemagglutinin inhibition (HI) tests confirmed that all the reassortant influenza B viruses were B/Malaysia/2506/04-like viruses. Pre- and post-immunized serum samples from 4 normal volunteers inoculated with 2007 influenza vaccine were evaluated for their HI activity on 6 reassortant B isolates including two variants that we found in the Yilan county. The results demonstrated that after vaccination, all four vaccinees had at least 4-fold increases of their HI titers.

Conclusion

The results indicate that the 2006–2007 seasonal influenza vaccine was effective in stimulating protective immunity against the influenza B variants identified in Yilan county. Continuous surveillance of emerging influenza B variants in the northern region of Taiwan is important for the selection of proper vaccine candidate in the future.

Epidemiologic study of influenza infection in Okinawa, Japan, from 2001 to 2007: changing patterns of seasonality and prevalence of amantadine-resistant influenza A virus

To clarify seasonal influenza patterns and the prevalence of amantadine-resistant influenza A viruses in Okinawa, located at the southern extremity of Japan in a subtropical climate, we conducted a laboratory-based study of influenza virus infections from 2001 to 2007. The annual outbreaks tended to show two peaks in Okinawa, in summer and winter, although the main islands of Japan, located in a temperate climate area, showed only winter influenza activity. Epidemic types and subtypes in Okinawa mostly matched those on the main islands of Japan in winter and those in Taiwan in summer. Rates of amantadine resistance dramatically increased, from 7.3% in the November 2002-to-March 2003 season to 90.0% in summer 2005, and a similarly high rate of resistance continued for the rest of the study period. Phylogenetic analysis of the hemagglutinin gene of A/H3N2 isolates collected from 2002 to 2007 revealed a monophyletic lineage that was divided into four period groups. Each group included amantadine-sensitive and -resistant viruses within independent clusters. In the November 2005-to-March 2006 season, all of the amantadine-resistant viruses were clustered in clade N, with dual (position 193 and 225) amino acid mutations in their HA1 subunits. In 2005, clade N amantadine-resistant viruses existed in Okinawa several months before the circulation of this clade on the main islands of Japan. In conclusion, surveillance in Okinawa to monitor influenza virus circulation is important for elucidating the dynamics of virus transmission in a border area between temperate and subtropical areas, as Okinawa is one of the best sentinel points in Japan.

Genetic and epidemiological analysis of influenza virus epidemics in Taiwan during 2003 to 2006

The genetic characterization of Taiwanese influenza A and B viruses on the basis of analyses of pairwise amino acid variations, genetic clustering, and phylogenetics was performed. A total of 548, 2,123, and 1,336 sequences of the HA1 genes of influenza A virus subtypes H1 and H3 and influenza B virus, respectively, collected during 2003 to 2006 from an island-wide surveillance network were determined. Influenza A virus H3 showed activity during all periods, although it was dominant only in the winters of 2002-2003 and 2003-2004. Instead, influenza B virus and influenza A virus H1 were dominant in the winters of 2004-2005 and 2005-2006, respectively. Additionally, two influenza A virus H3 peaks were found in the summers of 2004 and 2005. From clustering analysis, similar characteristics of high sequence diversity and short life spans for the influenza A virus H1 and H3 clusters were observed, despite their distinct seasonal patterns. In contrast, clusters with longer life spans and fewer but larger clusters were found among the influenza B viruses. We also noticed that more amino acid changes at antigenic sites, especially at sites B and D in the H3 viruses, were found in 2003 and 2004 than in the following 2 years. The only epidemic of the H1 viruses, which occurred in the winter of 2005-2006, was caused by two genetically distinct lineages, and neither of them showed apparent antigenic changes compared with the antigens of the vaccine strain. For the influenza B viruses, the multiple dominant lineages of Yamagata-like strains with large genetic variations observed reflected the evolutionary pressure caused by the Yamagata-like vaccine strain. On the other hand, only one dominant lineage of Victoria-like strains circulated from 2004 to 2006.