Stochastic processes are key determinants of short-term evolution in influenza a virus

Understanding the evolutionary dynamics of influenza A virus is central to its surveillance and control. While immune-driven antigenic drift is a key determinant of viral evolution across epidemic seasons, the evolutionary processes shaping influenza virus diversity within seasons are less clear. Here we show with a phylogenetic analysis of 413 complete genomes of human H3N2 influenza A viruses collected between 1997 and 2005 from New York State, United States, that genetic diversity is both abundant and largely generated through the seasonal importation of multiple divergent clades of the same subtype. These clades cocirculated within New York State, allowing frequent reassortment and generating genome-wide diversity. However, relatively low levels of positive selection and genetic diversity were observed at amino acid sites considered important in antigenic drift. These results indicate that adaptive evolution occurs only sporadically in influenza A virus; rather, the stochastic processes of viral migration and clade reassortment play a vital role in shaping short-term evolutionary dynamics. Thus, predicting future patterns of influenza virus evolution for vaccine strain selection is inherently complex and requires intensive surveillance, whole-genome sequencing, and phenotypic analysis.     

Exploring the Molecular Epidemiology and Evolutionary Dynamics of Influenza A Virus in Taiwan

The evolution and population dynamics of human influenza in Taiwan is a microcosm of the viruses circulating worldwide, which has not yet been studied in detail. We collected 343 representative full genome sequences of human influenza A viruses isolated in Taiwan between 1979 and 2009. Phylogenetic and antigenic data analysis revealed that H1N1 and H3N2 viruses consistently co-circulated in Taiwan, although they were characterized by different temporal dynamics and degrees of genetic diversity. Moreover, influenza A viruses of both subtypes underwent internal gene reassortment involving all eight segments of the viral genome, some of which also occurred during non-epidemic periods. The patterns of gene reassortment were different in the two subtypes. The internal genes of H1N1 viruses moved as a unit, separately from the co-evolving HA and NA genes. On the other hand, the HA and NA genes of H3N2 viruses tended to segregate consistently with different sets of internal gene segments. In particular, as reassortment occurred, H3HA always segregated as a group with the PB1, PA and M genes, while N2NA consistently segregated with PB2 and NP. Finally, the analysis showed that new phylogenetic lineages and antigenic variants emerging in summer were likely to be the progenitors of the epidemic strains in the following season. The synchronized seasonal patterns and high genetic diversity of influenza A viruses observed in Taiwan make possible to capture the evolutionary dynamic and epidemiological rules governing antigenic drift and reassortment and may serve as a “warning” system that recapitulates the global epidemic.     

Allele dynamics plots for the study of evolutionary dynamics in viral populations

Phylodynamic techniques combine epidemiological and genetic information to analyze the evolutionary and spatiotemporal dynamics of rapidly evolving pathogens, such as influenza A or human immunodeficiency viruses. We introduce 'allele dynamics plots' (AD plots) as a method for visualizing the evolutionary dynamics of a gene in a population. Using AD plots, we propose how to identify the alleles that are likely to be subject to directional selection. We analyze the method's merits with a detailed study of the evolutionary dynamics of seasonal influenza A viruses. AD plots for the major surface protein of seasonal influenza A (H3N2) and the 2009 swine-origin influenza A (H1N1) viruses show the succession of substitutions that became fixed in the evolution of the two viral populations. They also allow the early identification of those viral strains that later rise to predominance, which is important for the problem of vaccine strain selection. In summary, we describe a technique that reveals the evolutionary dynamics of a rapidly evolving population and allows us to identify alleles and associated genetic changes that might be under directional selection. The method can be applied for the study of influenza A viruses and other rapidly evolving species or viruses.     

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.     

Multiple reassortment events in the evolutionary history of H1N1 influenza A virus since 1918

The H1N1 subtype of influenza A virus has caused substantial morbidity and mortality in humans, first documented in the global pandemic of 1918 and continuing to the present day. Despite this disease burden, the evolutionary history of the A/H1N1 virus is not well understood, particularly whether there is a virological basis for several notable epidemics of unusual severity in the 1940s and 1950s. Using a data set of 71 representative complete genome sequences sampled between 1918 and 2006, we show that segmental reassortment has played an important role in the genomic evolution of A/H1N1 since 1918. Specifically, we demonstrate that an A/H1N1 isolate from the 1947 epidemic acquired novel PB2 and HA genes through intra-subtype reassortment, which may explain the abrupt antigenic evolution of this virus. Similarly, the 1951 influenza epidemic may also have been associated with reassortant A/H1N1 viruses. Intra-subtype reassortment therefore appears to be a more important process in the evolution and epidemiology of H1N1 influenza A virus than previously realized.     

Reassortment between seasonal H1N1 and pandemic (H1N1) 2009 influenza viruses is restricted by limited compatibility among polymerase subunits

Reassortment is important for influenza virus evolution and the generation of novel viruses with pandemic potential; however, the factors influencing reassortment are still poorly understood. Here, using reverse genetics and a replicon assay, we demonstrated that a mixed polymerase complex containing a pandemic (H1N1) 2009 influenza virus PB2 on a seasonal H1N1 virus background has reduced polymerase activity, leading to impaired virus viability. Adaptation of viruses containing the mixed polymerase complex resulted in compensatory mutations in PB1. Taken together, our results identify the cooperation between PB2 and PB1 as an important restricting factor for reassortment of influenza viruses.     

The Evolutionary Dynamics of Human Influenza B Virus

Despite their close phylogenetic relationship, type A and B influenza viruses exhibit major epidemiological differences in humans, with the latter both less common and less often associated with severe disease. However, it is unclear what processes determine the evolutionary dynamics of influenza B virus, and how influenza viruses A and B interact at the evolutionary scale. To address these questions we inferred the phylogenetic history of human influenza B virus using complete genome sequences for which the date (day) of isolation was available. By comparing the phylogenetic patterns of all eight viral segments we determined the occurrence of segment reassortment over a 30-year sampling period. An analysis of rates of nucleotide substitution and selection pressures revealed sporadic occurrences of adaptive evolution, most notably in the viral hemagglutinin and compatible with the action of antigenic drift, yet lower rates of overall and nonsynonymous nucleotide substitution compared to influenza A virus. Overall, these results led us to propose a model in which evolutionary changes within and between the antigenically distinct 'Yam88' and 'Vic87' lineages of influenza B virus are the result of changes in herd immunity, with reassortment continuously generating novel genetic variation. Additionally, we suggest that the interaction with influenza A virus may be central in shaping the evolutionary dynamics of influenza B virus, facilitating the shift of dominance between the Vic87 and the Yam88 lineages.     

Revelation of Influencing Factors in Overall Codon Usage Bias of Equine Influenza Viruses

Equine influenza viruses (EIVs) of H3N8 subtype are culprits of severe acute respiratory infections in horses, and are still responsible for significant outbreaks worldwide. Adaptability of influenza viruses to a particular host is significantly influenced by their codon usage preference, due to an absolute dependence on the host cellular machinery for their replication. In the present study, we analyzed genome-wide codon usage patterns in 92 EIV strains, including both H3N8 and H7N7 subtypes by computing several codon usage indices and applying multivariate statistical methods. Relative synonymous codon usage (RSCU) analysis disclosed bias of preferred synonymous codons towards A/U-ended codons. The overall codon usage bias in EIVs was slightly lower, and mainly affected by the nucleotide compositional constraints as inferred from the RSCU and effective number of codon (ENc) analysis. Our data suggested that codon usage pattern in EIVs is governed by the interplay of mutation pressure, natural selection from its hosts and undefined factors. The H7N7 subtype was found less fit to its host (horse) in comparison to H3N8, by possessing higher codon bias, lower mutation pressure and much less adaptation to tRNA pool of equine cells. To the best of our knowledge, this is the first report describing the codon usage analysis of the complete genomes of EIVs. The outcome of our study is likely to enhance our understanding of factors involved in viral adaptation, evolution, and fitness towards their hosts.     

Molecular epidemiology of A/H3N2 and A/H1N1 influenza virus during a single epidemic season in the United States

To determine the spatial and temporal dynamics of influenza A virus during a single epidemic, we examined whole-genome sequences of 284 A/H1N1 and 69 A/H3N2 viruses collected across the continental United States during the 2006-2007 influenza season, representing the largest study of its kind undertaken to date. A phylogenetic analysis revealed that multiple clades of both A/H1N1 and A/H3N2 entered and co-circulated in the United States during this season, even in localities that are distant from major metropolitan areas, and with no clear pattern of spatial spread. In addition, co-circulating clades of the same subtype exchanged genome segments through reassortment, producing both a minor clade of A/H3N2 viruses that appears to have re-acquired sensitivity to the adamantane class of antiviral drugs, as well as a likely antigenically distinct A/H1N1 clade that became globally dominant following this season. Overall, the co-circulation of multiple viral clades during the 2006-2007 epidemic season revealed patterns of spatial spread that are far more complex than observed previously, and suggests a major role for both migration and reassortment in shaping the epidemiological dynamics of human influenza A virus.     

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.     

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.     

Evolutionary complexities of swine flu H1N1 gene sequences of 2009

A recently emerged novel influenza A (H1N1) virus continues to spread globally. The pandemic caused by this new H1N1 swine influenza virus presents an opportunity to analyze the evolutionary significance of the origin of the new strain of swine flu. Our study clearly suggests that strong purifying selection is responsible for the evolution of the novel influenza A (H1N1) virus among human. We observed that the 2009 viral sequences are evolutionarily widely different from the past few years’ sequences. Rather, the 2009 sequences are evolutionarily more similar to the most ancient sequence reported in the NCBI Influenza Virus Resource Database collected in 1918. Analysis of evolutionary rates also supports the view that all the genes in the pandemic strain of 2009 except NA and M genes are derived from triple reassorted swine viruses. Our study demonstrates the importance of using complete-genome approach as more sequences will become available to investigate the evolutionary origin of the 1918 influenza A (H1N1) swine flu strain and the possibility of future reassortment events.     

Influenza A virus neuraminidase limits viral superinfection

Enveloped viruses use multiple mechanisms to inhibit infection of a target cell by more than one virion. These mechanisms may be of particular importance for the evolution of segmented viruses, because superinfection exclusion may limit the frequency of reassortment of viral genes. Here, we show that cellular expression of influenza A virus neuraminidase (NA), but not hemagglutinin (HA) or the M2 proton pump, inhibits entry of HA-pseudotyped retroviruses. Cells infected with H1N1 or H3N2 influenza A virus were similarly refractory to HA-mediated infection and to superinfection with a second influenza A virus. Both HA-mediated entry and viral superinfection were rescued by the neuraminidase inhibitors oseltamivir carboxylate and zanamivir. These inhibitors also prevented the removal of alpha-2,3- and alpha-2,6-linked sialic acid observed in cells expressing NA or infected with influenza A viruses. Our data indicate that NA alone among viral proteins limits influenza A virus superinfection.     

1968−2014年中国甲型H3N2流感病毒PB1基因进化分析

【目的】分析季节性H3N2流感病毒PB1基因序列的变异情况,揭示H3N2流感病毒PB1基因的分子特征与进化趋势。【方法】对1968?2014年中国地区82株人H3N2毒株、2012?2014年江苏省分离的81株甲型H3N2流感病毒、6株SIV和4株AIV H3N2亚型PB1、PB1-F2基因进行分子进化分析。【结果】1968?2014年中国H3N2流感毒株PB1核苷酸和氨基酸相似性分别为90.91%?100%和96.91%?100%。系统进化树分析,1968?2014年共173株H3N2流感病毒总体上分为4个分支,2002?2014年分离毒株位于第IV分支上,1968?1994年分离毒株位于第II和III分支;猪源H3N2亚型分布于第I、II、IV分支上;分子特征显示PB1氨基酸52、113、179、216、576、586、619、621、709位在2002年以后发生适应性改变,替换了原来的氨基酸;PB1-F2基因编码截断型蛋白长度有52、34、25、24、11 aa (猪源),PB1-F2蛋白毒力关键位点上未出现高致病性特征突变。【结论】自1968年起H3N2亚型PB1基因变异逐步趋于稳定,且PB1-F2截断型毒株正逐渐成为一类新的进化特征,但PB1基因与其他亚型之间发生重配以及关键毒力位点的变异仍应是监测的重点。     

Prospective of Genomics in Revealing Transmission, Reassortment and Evolution of Wildlife-Borne Avian Influenza A (H5N1) Viruses

The outbreak of highly pathogenic avian influenza (HPAI) H5N1 disease has led to significant loss of poultry and wild life and case fatality rates in humans of 60%. Wild birds are natural hosts for all avian influenza virus subtypes and over120 bird species have been reported with evidence of H5N1 infection. Influenza A viruses possess a segmented RNA genome and are characterized by frequently occurring genetic reassortment events, which play a very important role in virus evolution and the spread of novel gene constellations in immunologically naïve human and animal populations. Phylogenetic analysis of whole genome or sub-genomic sequences is a standard means for delineating genetic variation, novel reassortment events, and surveillance to trace the global transmission pathways. In this paper, special emphasis is given to the transmission and circulation of H5N1 among wild life populations, and to the reassortment events that are associated with inter-host transmission of the H5N1 viruses when they infect different hosts, such as birds, pigs and humans. In addition, we review the inter-subtype reassortment of the viral segments encoding inner proteins between the H5N1 viruses and viruses of other subtypes, such as H9N2 and H6N1. Finally, we highlight the usefulness of genomic sequences in molecular epidemiological analysis of HPAI H5N1 and the technical limitations in existing analytical methods that hinder them from playing a greater role in virological research.     

Genetically Diverse Low Pathogenicity Avian Influenza A Virus Subtypes Co-Circulate among Poultry in Bangladesh

Influenza virus surveillance, poultry outbreak investigations and genomic sequencing were assessed to understand the ecology and evolution of low pathogenicity avian influenza (LPAI) A viruses in Bangladesh from 2007 to 2013. We analyzed 506 avian specimens collected from poultry in live bird markets and backyard flocks to identify influenza A viruses. Virus isolation-positive specimens (n = 50) were subtyped and their coding-complete genomes were sequenced. The most frequently identified subtypes among LPAI isolates were H9N2, H11N3, H4N6, and H1N1. Less frequently detected subtypes included H1N3, H2N4, H3N2, H3N6, H3N8, H4N2, H5N2, H6N1, H6N7, and H7N9. Gene sequences were compared to publicly available sequences using phylogenetic inference approaches. Among the 14 subtypes identified, the majority of viral gene segments were most closely related to poultry or wild bird viruses commonly found in Southeast Asia, Europe, and/or northern Africa. LPAI subtypes were distributed over several geographic locations in Bangladesh, and surface and internal protein gene segments clustered phylogenetically with a diverse number of viral subtypes suggesting extensive reassortment among these LPAI viruses. H9N2 subtype viruses differed from other LPAI subtypes because genes from these viruses consistently clustered together, indicating this subtype is enzootic in Bangladesh. The H9N2 strains identified in Bangladesh were phylogenetically and antigenically related to previous human-derived H9N2 viruses detected in Bangladesh representing a potential source for human infection. In contrast, the circulating LPAI H5N2 and H7N9 viruses were both phylogenetically and antigenically unrelated to H5 viruses identified previously in humans in Bangladesh and H7N9 strains isolated from humans in China. In Bangladesh, domestic poultry sold in live bird markets carried a wide range of LPAI virus subtypes and a high diversity of genotypes. These findings, combined with the seven year timeframe of sampling, indicate a continuous circulation of these viruses in the country.     

PCR-based molecular characterization and insilico analysis of food-borne trematode parasites Paragonimus westermani, Fasciolopsis buski and Fasciola gigantica from Northeast India using ITS2 rDNA

In early 2009, new swine-origin influenza A (H1N1) virus emerged in Mexico and the United States. The emerging influenza virus had made global influenza pandemic for nearly one year. To every emerging pathogen, exploring the origin sources is vital for viral control and clearance. Influenza virus is different from other virus in that it has 8 segments, making the segment reassortment a main drive in virus evolution. In exploring reassortment evolution origins of a newly emerging influenza virus, integrated comparing of the origin sources of all the segments is necessary. If some segments have high homologous with one parental strain, lower homologous with another parental strain, while other segments are reverse, can we proposed that this emerging influenza virus may re-assort from the two parental strains. Here we try to explore the multilevel reassortment evolution origins of 2009 H1N1 influenza virus using this method. By further validating the fidelity of this strategy, this method might be useful in judging the reassortment origins of newly emerging influenza virus.     

Virulence and genetic compatibility of polymerase reassortant viruses derived from the pandemic (H1N1) 2009 influenza virus and circulating influenza A viruses

Gene mutations and reassortment are key mechanisms by which influenza A virus acquires virulence factors. To evaluate the role of the viral polymerase replication machinery in producing virulent pandemic (H1N1) 2009 influenza viruses, we generated various polymerase point mutants (PB2, 627K/701N; PB1, expression of PB1-F2 protein; and PA, 97I) and reassortant viruses with various sources of influenza viruses by reverse genetics. Although the point mutations produced no significant change in pathogenicity, reassortment between the pandemic A/California/04/09 (CA04, H1N1) and current human and animal influenza viruses produced variants possessing a broad spectrum of pathogenicity in the mouse model. Although most polymerase reassortants had attenuated pathogenicity (including those containing seasonal human H3N2 and high-pathogenicity H5N1 virus segments) compared to that of the parental CA04 (H1N1) virus, some recombinants had significantly enhanced virulence. Unexpectedly, one of the five highly virulent reassortants contained a A/Swine/Korea/JNS06/04(H3N2)-like PB2 gene with no known virulence factors; the other four had mammalian-passaged avian-like genes encoding PB2 featuring 627K, PA featuring 97I, or both. Overall, the reassorted polymerase complexes were only moderately compatible for virus rescue, probably because of disrupted molecular interactions involving viral or host proteins. Although we observed close cooperation between PB2 and PB1 from similar virus origins, we found that PA appears to be crucial in maintaining viral gene functions in the context of the CA04 (H1N1) virus. These observations provide helpful insights into the pathogenic potential of reassortant influenza viruses composed of the pandemic (H1N1) 2009 influenza virus and prevailing human or animal influenza viruses that could emerge in the future.     

Evolution of novel reassortant A/H3N2 influenza viruses in North American swine and humans, 2009-2011

Novel H3N2 influenza viruses (H3N2v) containing seven genome segments from swine lineage triple-reassortant H3N2 viruses and a 2009 pandemic H1N1 (H1N1pdm09) matrix protein segment (pM) were isolated from 12 humans in the United States between August and December 2011. To understand the evolution of these novel H3N2 viruses in swine and humans, we undertook a phylogenetic analysis of 674 M sequences and 388 HA and NA sequences from influenza viruses isolated from North American swine during 2009-2011, as well as HA, NA, and M sequences from eight H3N2v viruses isolated from humans. We identified 34 swine influenza viruses (termed rH3N2p) with the same combination of H3, N2, and pM segments as the H3N2v viruses isolated from humans. Notably, these rH3N2p viruses were generated in swine via reassortment events between H3N2 viruses and the pM segment approximately 4 to 10 times since 2009. The pM segment has also reassorted with multiple distinct lineages of H1 virus, especially H1δ viruses. Importantly, the N2 segment of all H3N2v viruses isolated from humans is derived from a genetically distinct N2 lineage that has circulated in swine since being acquired by reassortment with seasonal human H3N2 viruses in 2001-2002, rather than from the N2 that is associated with the 1998 H3N2 swine lineage. The identification of this N2 variant may have implications for influenza vaccine design and the potential pandemic threat of H3N2v to human age groups with differing levels of prior exposure and immunity.     

Genome Sequence of a Natural Reassortant H5N2 Avian Influenza Virus from Domestic Mallard Ducks in Eastern China

Here, we report the genomic sequence of a Chinese reassortant H5N2 avian influenza virus which possessed the polybasic motif PLREKRRK-R/GL at the hemagglutinin cleavage site. Phylogenetic analysis showed that all eight genes were of the Eurasian lineage, five of which were highly homologous to the endemic clade 2.3.4 H5N1 viruses and their H5N5 reassortant descendants. These data suggested that novel multisubtypic NA reassortants bearing the H5N1 backbone could be generated through genetic reassortment in H5N1 circulating regions, which will help in understanding the evolution and segment reassortment mechanism of H5 subtype avian influenza viruses.