Phylogenetic diversity and genotypical complexity of H9N2 influenza A viruses revealed by genomic sequence analysis

H9N2 influenza A viruses have become established worldwide in terrestrial poultry and wild birds, and are occasionally transmitted to mammals including humans and pigs. To comprehensively elucidate the genetic and evolutionary characteristics of H9N2 influenza viruses, we performed a large-scale sequence analysis of 571 viral genomes from the NCBI Influenza Virus Resource Database, representing the spectrum of H9N2 influenza viruses isolated from 1966 to 2009. Our study provides a panoramic framework for better understanding the genesis and evolution of H9N2 influenza viruses, and for describing the history of H9N2 viruses circulating in diverse hosts. Panorama phylogenetic analysis of the eight viral gene segments revealed the complexity and diversity of H9N2 influenza viruses. The 571 H9N2 viral genomes were classified into 74 separate lineages, which had marked host and geographical differences in phylogeny. Panorama genotypical analysis also revealed that H9N2 viruses include at least 98 genotypes, which were further divided according to their HA lineages into seven series (A-G). Phylogenetic analysis of the internal genes showed that H9N2 viruses are closely related to H3, H4, H5, H7, H10, and H14 subtype influenza viruses. Our results indicate that H9N2 viruses have undergone extensive reassortments to generate multiple reassortants and genotypes, suggesting that the continued circulation of multiple genotypical H9N2 viruses throughout the world in diverse hosts has the potential to cause future influenza outbreaks in poultry and epidemics in humans. We propose a nomenclature system for identifying and unifying all lineages and genotypes of H9N2 influenza viruses in order to facilitate international communication on the evolution, ecology and epidemiology of H9N2 influenza viruses.     

Rapid evolution of H5N1 influenza viruses in chickens in Hong Kong

The H5N1 avian influenza virus that killed 6 of 18 persons infected in Hong Kong in 1997 was transmitted directly from poultry to humans. Viral isolates from this outbreak may provide molecular clues to zoonotic transfer. Here we demonstrate that the H5N1 viruses circulating in poultry comprised two distinguishable phylogenetic lineages in all genes that were in very rapid evolution. When introduced into new hosts, influenza viruses usually undergo rapid alteration of their surface glycoproteins, especially in the hemagglutinin (HA). Surprisingly, these H5N1 isolates had a large proportion of amino acid changes in all gene products except in the HA. These viruses maybe reassortants each of whose HA gene is well adapted to domestic poultry while the rest of the genome arises from a different source. The consensus amino acid sequences of "internal" virion proteins reveal amino acids previously found in human strains. These human-specific amino acids may be important factors in zoonotic transmission.     

Novel Genotypes of H9N2 Influenza A Viruses Isolated from Poultry in Pakistan Containing NS Genes Similar to Highly Pathogenic H7N3 and H5N1 Viruses

The impact of avian influenza caused by H9N2 viruses in Pakistan is now significantly more severe than in previous years. Since all gene segments contribute towards the virulence of avian influenza virus, it was imperative to investigate the molecular features and genetic relationships of H9N2 viruses prevalent in this region. Analysis of the gene sequences of all eight RNA segments from 12 viruses isolated between 2005 and 2008 was undertaken. The hemagglutinin (HA) sequences of all isolates were closely related to H9N2 viruses isolated from Iran between 2004 and 2007 and contained leucine instead of glutamine at position 226 in the receptor binding pocket, a recognised marker for the recognition of sialic acids linked α2–6 to galactose. The neuraminidase (NA) of two isolates contained a unique five residue deletion in the stalk (from residues 80 to 84), a possible indication of greater adaptation of these viruses to the chicken host. The HA, NA, nucleoprotein (NP), and matrix (M) genes showed close identity with H9N2 viruses isolated during 1999 in Pakistan and clustered in the A/Quail/Hong Kong/G1/97 virus lineage. In contrast, the polymerase genes clustered with H9N2 viruses from India, Iran and Dubai. The NS gene segment showed greater genetic diversity and shared a high level of similarity with NS genes from either H5 or H7 subtypes rather than with established H9N2 Eurasian lineages. These results indicate that during recent years the H9N2 viruses have undergone extensive genetic reassortment which has led to the generation of H9N2 viruses of novel genotypes in the Indian sub-continent. The novel genotypes of H9N2 viruses may play a role in the increased problems observed by H9N2 to poultry and reinforce the continued need to monitor H9N2 infections for their zoonotic potential.     

中国H5N1禽流感病毒HA基因的分子进化分析

为了阐明中国H5N1禽流感病毒血凝素(HA)基因的分子进化规律,从GenBank下载中国近年来分离的27株禽流感H5N1病毒的HA基因,利用LaserGene软件包中的EditSeq程序剪切共同序列并翻译成氨基酸序列,与DNAMAN软件进行比对,并以LaserGene软件包中的MegAlign程序构建分子进化树,结合背景资料分析其传播流行规律。结果显示:①在HA裂解位点上游毗邻的位置,gdch04、gddu04、gxch0405、gxdu05、gxqu0502、hkcph05、stch05、stqu05、ynch0502株都缺失了三个碱基;②27株毒株在HA裂解位点上游毗邻的位置存在多个连续的碱性氨基酸;③来自广西的gxch0405、gxdu05、gxqu0502和来自香港的hkcph05、hkgh04的序列同源性很高,很可能与苍鹭的迁移有关;④来自广西的gxqu0502和来自香港的hkcph05亲源关系很近,很可能来自同一祖先。     

H9N2 influenza virus in China: a cause of concern

The recent human infection with avian influenza virus revealed that H9N2 influenza virus is the gene donor for H7N9 and H10N8 viruses infecting humans. The crucial role of H9N2 viruses at the animal-human interface might be due to the wide host range, adaptation in both poultry and mammalian, and extensive gene reassortment. As the most prevalent subtype of influenza viruses in chickens in China, H9N2 also causes a great economic loss for the poultry industry, even under the long-term vaccination programs. The history, epidemiology, biological characteristics, and molecular determinants of H9N2 influenza virus are reviewed in this paper. The contribution of H9N2 genes, especially RNP genes, to the infection of humans needs to be investigated in the future.     

Origin and molecular characterization of the human-infecting H6N1 influenza virus in Taiwan

In June 2013, the first human H6N1 influenza virus infection was confirmed in Taiwan. However, the origin and molecular characterization of this virus, A/Taiwan/2/2013 (H6N1), have not been well studied thus far. In the present report, we performed phylogenetic and coalescent analyses of this virus and compared its molecular profile/characteristics with other closely related strains. Molecular characterization of H6N1 revealed that it is a typical avian influenza virus of low pathogenicity, which might not replicate and propagate well in the upper airway in mammals. Phylogenetic analysis revealed that the virus clusters with A/chicken/Taiwan/A2837/2013 (H6N1) in seven genes, except PB1. For the PB1 gene, A/Taiwan/2/2013 was clustered with a different H6N1 lineage from A/chicken/Taiwan/A2837/2013. Although a previous study demonstrated that the PB2, PA, and M genes of A/Taiwan/2/2013 might be derived from the H5N2 viruses, coalescent analyses revealed that these H5N2 viruses were derived from more recent strains than that of the ancestor of A/Taiwan/2/2013. Therefore, we propose that A/Taiwan/2/2013 is a reassortant from different H6N1 lineages circulating in chickens in Taiwan. Furthermore, compared to avian isolates, a single P186L (H3 numbering) substitution in the hemagglutinin H6 of the human isolate might increase the mammalian receptor binding and, hence, this strain’s pathogenicity in humans. Overall, human infection with this virus seems an accidental event and is unlikely to cause an influenza pandemic. However, its co-circulation and potential reassortment with other influenza subtypes are still worthy of attention.     

Molecular Evolution of H6 Influenza Viruses from Poultry in Southeastern China: Prevalence of H6N1 Influenza Viruses Possessing Seven A/Hong Kong/156/97 (H5N1)-Like Genes in Poultry

The A/teal/Hong Kong/W312/97 (H6N1) influenza virus and the human H5N1 and H9N2 influenza viruses possess similar genes encoding internal proteins, suggesting that H6N1 viruses could become novel human pathogens. The molecular epidemiology and evolution of H6 influenza viruses were characterized by antigenic and genetic analyses of 29 H6 influenza viruses isolated from 1975 to 1981 and 1997 to 2000. Two distinct groups were identified on the basis of their antigenic characteristics. Phylogenetic analysis revealed that all H6N1 viruses isolated from terrestrial poultry in 1999 and 2000 are closely related to A/teal/Hong Kong/W312/97 (H6N1), and the nucleotide sequences of these viruses and of A/Hong Kong/156/97 (H5N1) were more than 96% homologous. The hemagglutinin (HA) of the 1999 and 2000 terrestrial viruses does not have multiple basic amino acids at the site of cleavage of HA1 to HA2; however, a unique insertion of aspartic acid in HA1 between positions 144 and 145 (H3 numbering) was found. The neuraminidase of these terrestrial H6N1 viruses has a deletion of 19 amino acids characteristic of A/Hong Kong/156/97 (H5N1). Evolutionary analysis suggested that these H6N1 viruses coevolved with A/quail/Hong Kong/G1/97-like H9N2 viruses and became more adapted to terrestrial poultry. These terrestrial 1999 and 2000 A/teal/Hong Kong/W312/97 (H6N1)-like viruses, along with the H9N2 viruses, could have been involved in the genesis of the pathogenic H5N1 influenza viruses of 1997. The presence of H6N1 viruses in poultry markets in Hong Kong that possess seven of the eight genes of the A/Hong Kong/156/97 (H5N1) virus raises the following fundamental questions relevant to influenza pandemic preparedness: could the pathogenic H5N1 virus reemerge and could the H6N1 viruses directly cross the species barrier to mammals?     

Molecular Epidemiology of Influenza A/H3N2 Viruses Circulating in Mexico from 2003 to 2012

In this work, nineteen influenza A/H3N2 viruses isolated in Mexico between 2003 and 2012 were studied. Our findings show that different human A/H3N2 viral lineages co-circulate within a same season and can also persist locally in between different influenza seasons, increasing the chance for genetic reassortment events. A novel minor cluster was also identified, named here as Korea, that circulated worldwide during 2003. Frequently, phylogenetic characterization did not correlate with the determined antigenic identity, supporting the need for the use of molecular evolutionary tools additionally to antigenic data for the surveillance and characterization of viral diversity during each flu season. This work represents the first long-term molecular epidemiology study of influenza A/H3N2 viruses in Mexico based on the complete genomic sequences and contributes to the monitoring of evolutionary trends of A/H3N2 influenza viruses within North and Central America.     

Whole-genome analysis of human influenza A virus reveals multiple persistent lineages and reassortment among recent H3N2 viruses

Understanding the evolution of influenza A viruses in humans is important for surveillance and vaccine strain selection. We performed a phylogenetic analysis of 156 complete genomes of human H3N2 influenza A viruses collected between 1999 and 2004 from New York State, United States, and observed multiple co-circulating clades with different population frequencies. Strikingly, phylogenies inferred for individual gene segments revealed that multiple reassortment events had occurred among these clades, such that one clade of H3N2 viruses present at least since 2000 had provided the hemagglutinin gene for all those H3N2 viruses sampled after the 2002–2003 influenza season. This reassortment event was the likely progenitor of the antigenically variant influenza strains that caused the A/Fujian/411/2002-like epidemic of the 2003–2004 influenza season. However, despite sharing the same hemagglutinin, these phylogenetically distinct lineages of viruses continue to co-circulate in the same population. These data, derived from the first large-scale analysis of H3N2 viruses, convincingly demonstrate that multiple lineages can co-circulate, persist, and reassort in epidemiologically significant ways, and underscore the importance of genomic analyses for future influenza surveillance.     

Envelope Proteins Pertain with Evolution and Adaptive Mechanism of the Novel Influenza A/H1N1 in Humans

The novel swine-origin influenza A/H1N1 virus (S-OIV) first detected in April 2009 has been identified to transmit from human to human directly and is the cause of currently emerged pandemic. In this study, nucleotide and deduced amino acid sequences of hemagglutinin (HA) and neuraminidase (NA) of the S-OIV and other influenza A viruses were analyzed through bioinformatic tools for phylogenetic analysis, genetic recombination and point mutation to investigate the emergence and adaptation of the S-OIV in human. The phylogenetic analysis showed that the HA comes from triple reassortant influenza A/H1N2 and the NA from Eurasian swine influenza A/H1N1 indicating HA and NA to descend from different lineages during the genesis of the S-OIV. Recombination analysis nullified the possibility of occurrence of recombination in HA and NA denoting the role of reassortment in the outbreak. Several conservative mutations are observed in the amino acid sequences of the HA and NA and this mutated residues are identical in the S-OIV. The results reported herein suggested the notion that the recent pandemic is the result of reassortment of different genes from different lineages of two envelope proteins, HA and NA which are responsible for antigenic activity of virus. This study further suggests that the adaptive capability of the S-OIV in human is acquired by the unique mutations generated during emergence.     

甲型H1N1流感病毒血凝素单克隆抗体轻链和重链可变区基因的巢式PCR扩增及序列分析

目的：采用巢式PCR对甲型H1N1流感病毒血凝素单克隆抗体的轻链和重链基因进行扩增,对获得的基因进行序列分析,并找出克隆鼠Igκ轻链和重链可变区基因的通用方法。方法：设计22对扩增鼠Igκ轻链可变区和重链可变区基因的引物,对6株鼠抗人甲型H1N1流感病毒血凝素单克隆抗体的轻链和重链可变区基因进行克隆并测序,与NCBI公布的鼠免疫球蛋白序列比对分析。结果：巢式PCR方法可以有效避免单克隆抗体克隆过程的假基因,并且得到的单克隆抗体的氨基酸序列均符合鼠免疫球蛋白可变区特征。结论：建立了克隆鼠免疫球蛋白轻链和重链可变区基因的通用方法,为后期克隆鼠源性单克隆抗体的可变区基因提供了基础,并为研究甲型H1N1流感病毒血凝素与抗体的结合位点提供了实验数据。     

Production and characterization of monoclonal antibodies against the hemagglutinin of H5N1 and antigenic investigation of avian influenza H5N1 viruses isolated from China

Eight monoclonal antibodies against hemagglutinin of influenza A virus A/Chicken/Henan/01/2004(H5N1) were produced by a DNA prime and inactivated virions-boost immunization strategy. Among the monoclonal antibodies, 3 (H50, H56, and H57) exhibited hemagglutination inhibition activity. Western blot analyses revealed that all the monoclonal antibodies reacted to the prokaryotically expressed HA1 of A/Chicken/Henan/01/2004(H5N1). The monoclonal antibodies were then used to characterize 10 avian influenza H5N1 viruses isolated from China during 2004 to 2007, by using the hemagglutination inhibition test and the antigen-capture enzyme-linked immunosorbent assay. The isolates could be divided into 4 different antigenic groups according to their responses to the monoclonal antibodies. The antigenic grouping of these 10 H5N1 isolates, using these antibodies, did not completely match their phylogenetic classification based on the hemagglutinin sequences. The results showed there were antigenic variations within the subclade 2.3.4 of H5N1, which is predominant in China.     

感染CDV的犬病料N、H、F基因的克隆与序列分析

目的了解病料中犬瘟热病毒（CDV）的变异情况,为犬瘟热的预防与治疗提供理论依据。方法采用RT-PCR方法对犬病料中CDV的血凝素蛋白（H）、融合蛋白（F）和核衣壳蛋白（N）基因进行扩增,将扩增产物克隆到pGEM-T-easy载体中测序,并进行序列分析。结果测定一株CDV的H、F和N基因大小分别为1863 bp、1653 bp和2235 bp。同源性分析表明,未发现碱基的插入和缺失。该病毒的H、F和N基因分别与国内强毒株BJ080127株的H基因、GN株的F基因和HL株的N基因亲缘关系最近,氨基酸同源性分别为97.3%、97.7%和99.3%。与国外标准野毒株A75-17在核苷酸水平上同源性依次为94.4%、93.8%和97.2%,与疫苗株CDV3在核苷酸水平上同源性分别为88.5%、88.8%和93.9%。系统进化树表明该病毒与国内强毒株在同一谱系。糖基化分析显示,该病毒在F基因3～5位氨基酸处多出1个糖基化位点。结论本研究从犬病料中成功克隆了犬瘟热病毒血凝素蛋白、融合蛋白和核衣壳蛋白,在F基因中新增了一个糖基化位点,为犬瘟热的遗传变异和流行学研究提供了分子生物学依据。     

Continuing evolution of H9N2 influenza viruses in Southeastern China

H9N2 influenza viruses are panzootic in domestic poultry in Eurasia and since 1999 have caused transient infections in humans and pigs. To investigate the zoonotic potential of H9N2 viruses, we studied the evolution of the viruses in live-poultry markets in Hong Kong in 2003. H9N2 was the most prevalent influenza virus subtype in the live-poultry markets between 2001 and 2003. Antigenic and phylogenetic analysis of hemagglutinin (HA) showed that all of the 19 isolates found except one belonged to the lineage represented by A/Duck/Hong Kong/Y280/97 (H9N2). The exception was A/Guinea fowl/NT184/03 (H9N2), whose HA is most closely related to that of the human isolate A/Guangzhou/333/99 (H9N2), a virus belonging to the A/Chicken/Beijing/1/94-like (H9N2) lineage. At least six different genotypes were recognized. The majority of the viruses had nonstructural (and HA) genes derived from the A/Duck/Hong Kong/Y280/97-like virus lineage but had other genes of mixed avian virus origin, including genes similar to those of H5N1 viruses isolated in 2001. Viruses of all six genotypes of H9N2 found were able to replicate in chickens and mice without adaptation. The infected chickens showed no signs of disease, but representatives of two viral genotypes were lethal to mice. Three genotypes of virus replicated in the respiratory tracts of swine, which shed virus for at least 5 days. These results show an increasing genetic and biologic diversity of H9N2 viruses in Hong Kong and support their potential role as pandemic influenza agents.     

Study of highly pathogenic H5N1 influenza virus isolated from sick and dead birds in Western Siberia

The mass destruction of domestic birds has been registered in July, 2005 in Novosibirsk region. Influenza virus H5N1 have been isolated from bodies of the lost birds on developing chicken embryos and identified by serological and molecular biological methods. M-gene and genes coding hemagglutinin and neurominidase were in part sequening. The phylogenetic analysis of hemagglutinin gene has shown, that the isolated viruses are forming a claster with strains, isolated from birds during outbreak of the bird's flu on lake Tsinghai (Qinghai) in China in 2005, in Japan in 2004, and also with H5N1 strains, isolated from the person and birds in the countries of Southeast Asia during ipizootia in 2003-2004. The site of the restriction which associated with pathogenicity of isolated avian influenza viruses H5 serogroup, corresponds to sequence of high pathogenic strains, circulating in the countries of Southeast Asia. The test for pathogenicity with use of chickens has confirmed, that researched strains were high pathogenic for birds.     

Comparisons of Highly Virulent H5N1 Influenza A Viruses Isolated from Humans and Chickens from Hong Kong

Genes of an influenza A (H5N1) virus from a human in Hong Kong isolated in May 1997 were sequenced and found to be all avian-like (K. Subbarao et al., Science 279:393–395, 1998). Gene sequences of this human isolate were compared to those of a highly pathogenic chicken H5N1 influenza virus isolated from Hong Kong in April 1997. Sequence comparisons of all eight RNA segments from the two viruses show greater than 99% sequence identity between them. However, neither isolate’s gene sequence was closely (>95% sequence identity) related to any other gene sequences found in the GenBank database. Phylogenetic analysis demonstrated that the nucleotide sequences of at least four of the eight RNA segments clustered with Eurasian origin avian influenza viruses. The hemagglutinin gene phylogenetic analysis also included the sequences from an additional three human and two chicken H5N1 virus isolates from Hong Kong, and the isolates separated into two closely related groups. However, no single amino acid change separated the chicken origin and human origin isolates, but they all contained multiple basic amino acids at the hemagglutinin cleavage site, which is associated with a highly pathogenic phenotype in poultry. In experimental intravenous inoculation studies with chickens, all seven viruses were highly pathogenic, killing most birds within 24 h. All infected chickens had virtually identical pathologic lesions, including moderate to severe diffuse edema and interstitial pneumonitis. Viral nucleoprotein was most frequently demonstrated in vascular endothelium, macrophages, heterophils, and cardiac myocytes. Asphyxiation from pulmonary edema and generalized cardiovascular collapse were the most likely pathogenic mechanisms responsible for illness and death. In summary, a small number of changes in hemagglutinin gene sequences defined two closely related subgroups, with both subgroups having human and chicken members, among the seven viruses examined from Hong Kong, and all seven viruses were highly pathogenic in chickens and caused similar lesions in experimental inoculations.     

Comparative study of the hemagglutinin and neuraminidase genes of influenza A virus H3N2, H9N2, and H5N1 subtypes using bioinformatics techniques

To investigate the genomic patterns of influenza A virus subtypes, such as H3N2, H9N2, and H5N1, we collected 1842 sequences of the hemagglutinin and neuraminidase genes from the NCBI database and parsed them into 7 categories: accession number, host species, sampling year, country, subtype, gene name, and sequence. The sequences that were isolated from the human, avian, and swine populations were extracted and stored in a MySQL database for intensive analysis. The GC content and relative synonymous codon usage (RSCU) values were calculated using JAVA codes. As a result, correspondence analysis of the RSCU values yielded the unique codon usage pattern (CUP) of each subtype and revealed no extreme differences among the human, avian, and swine isolates. H5N1 subtype viruses exhibited little variation in CUPs compared with other subtypes, suggesting that the H5N1 CUP has not yet undergone significant changes within each host species. Moreover, some observations may be relevant to CUP variation that has occurred over time among the H3N2 subtype viruses isolated from humans. All the sequences were divided into 3 groups over time, and each group seemed to have preferred synonymous codon patterns for each amino acid, especially for arginine, glycine, leucine, and valine. The bioinformatics technique we introduce in this study may be useful in predicting the evolutionary patterns of pandemic viruses.     

Genetic variability of isolates of pandemic influenza A virus H1N1 isolated in Russia in 2009

Complete nucleotide sequence of the genome segments encoding the surface glycoproteins, hemagglutinin, and neuraminidase of influenza A virus H1N1 derived from the patients with influenza in the context of pandemic (H1N1) 2009 was determined out of 14 isolates of pandemic influenza. The philogenetic analysis of these sequences demonstrated their genetic similarity to the corresponding genes of the pandemic influenza virus A (H1N1) 2009 isolates obtained in other countries; each gene homology was greater than 99%. Neuraminidase mutations causing virus resistance to oseltamivir and other neuraminidase inhibitors, known from the literature, were not detected. The hemagglutinin gene mutation D222G was found in 4 isolates from autopsy material. In the hemagglutinin of pandemic A/Salekhard/01/2009(H1N1) isolate a mutation G155E leading to the increase in viral replication in developing chick embryos was detected. The nature and frequency of nucleotides substitutions within HA and NA genes were determined in the current research.     

Chinese and Global Distribution of H9 Subtype Avian Influenza Viruses

H9 subtype avian influenza viruses (AIVs) are of significance in poultry and public health, but epidemiological studies about the viruses are scarce. In this study, phylogenetic relationships of the viruses were analyzed based on 1233 previously reported sequences and 745 novel sequences of the viral hemagglutinin gene. The novel sequences were obtained through large-scale surveys conducted in 2008-2011 in China. The results revealed distinct distributions of H9 subtype AIVs in different hosts, sites and regions in China and in the world: (1) the dominant lineage of H9 subtype AIVs in China in recent years is lineage h9.4.2.5 represented by A/chicken/Guangxi/55/2005; (2) the newly emerging lineage h9.4.2.6, represented by A/chicken/Guangdong/FZH/2011, has also become prevalent in China; (3) lineages h9.3.3, h9.4.1 and h9.4.2, represented by A/duck/Hokkaido/26/99, A/quail/Hong Kong/G1/97 and A/chicken/Hong Kong/G9/97, respectively, have become globally dominant in recent years; (4) lineages h9.4.1 and h9.4.2 are likely of more risk to public health than others; (5) different lineages have different transmission features and host tropisms. This study also provided novel experimental data which indicated that the Leu-234 (H9 numbering) motif in the viral hemagglutinin gene is an important but not unique determinant in receptor-binding preference. This report provides a detailed and updated panoramic view of the epidemiological distributions of H9 subtype AIVs globally and in China, and sheds new insights for the prevention of infection in poultry and preparedness for a potential pandemic caused by the viruses.     

Comparative analysis of evolutionary mechanisms of the hemagglutinin and three internal protein genes of influenza B virus: multiple cocirculating lineages and frequent reassortment of the NP, M, and NS genes

Phylogenetic profiles of the genes coding for the hemagglutinin (HA) protein, nucleoprotein (NP), matrix (M) protein, and nonstructural (NS) proteins of influenza B viruses isolated from 1940 to 1998 were analyzed in a parallel manner in order to understand the evolutionary mechanisms of these viruses. Unlike human influenza A (H3N2) viruses, the evolutionary pathways of all four genes of recent influenza B viruses revealed similar patterns of genetic divergence into two major lineages. Although evolutionary rates of the HA, NP, M, and NS genes of influenza B viruses were estimated to be generally lower than those of human influenza A viruses, genes of influenza B viruses demonstrated complex phylogenetic patterns, indicating alternative mechanisms for generation of virus variability. Topologies of the evolutionary trees of each gene were determined to be quite distinct from one another, showing that these genes were evolving in an independent manner. Furthermore, variable topologies were apparently the result of frequent genetic exchange among cocirculating epidemic viruses. Evolutionary analysis done in the present study provided further evidence for cocirculation of multiple lineages as well as sequestering and reemergence of phylogenetic lineages of the internal genes. In addition, comparison of deduced amino acid sequences revealed a novel amino acid deletion in the HA1 domain of the HA protein of recent isolates from 1998 belonging to the B/Yamagata/16/88-like lineage. It thus became apparent that, despite lower evolutionary rates, influenza B viruses were able to generate genetic diversity among circulating viruses through a combination of evolutionary mechanisms involving cocirculating lineages and genetic reassortment by which new variants with distinct gene constellations emerged.