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.     

Close relationship between the 2009 H1N1 virus and South Dakota AIV strains

Although previous publications suggest the 2009 pandemic influenza A (H1N1) virus was reassorted from swine viruses of North America and Eurasia, the immediate ancestry still remains elusive due to the big evolutionary distance between the 2009 H1N1 virus and the previously isolated strains. Since the unveiling of the 2009 H1N1 influenza, great deal of interest has been drawn to influenza, consequently a large number of influenza virus sequences have been deposited into the public sequence databases. Blast analysis demonstrated that the recently submitted 2007 South Dakota avian influenza virus strains and other North American avian strains contained genetic segments very closely related to the 2009 H1N1 virus, which suggests these avian influenza viruses are very close relatives of the 2009 H1N1 virus. Phylogenetic analyses also indicate that the 2009 H1N1 viruses are associated with both avian and swine influenza viruses circulating in North America. Since the migrating wild birds are preferable to pigs as the carrier to spread the influenza viruses across vast distances, it is very likely that birds played an important role in the inter-continental evolution of the 2009 H1N1 virus. It is essential to understand the evolutionary route of the emerging influenza virus in order to find a way to prevent further emerging cases. This study suggests the close relationship between 2009 pandemic virus and the North America avian viruses and underscores enhanced surveillance of influenza in birds for understanding the evolution of the 2009 pandemic influenza.     

The evolution of epidemic influenza

Recent developments in complete-genome sequencing, antigenic mapping and epidemiological modelling are greatly improving our knowledge of the evolution of human influenza virus at the epidemiological scale. In particular, recent studies have revealed a more complex relationship between antigenic evolution, natural selection and reassortment than previously realized. Despite these advances, there is much that remains to be understood about the epidemiology of influenza virus, particularly the processes that determine the virus's strong seasonality. We argue that a complete understanding of the evolutionary biology of this important human pathogen will require a genomic view of genetic diversity, including the acquisition of polymorphism data from within individual hosts and from geographical regions, particularly the tropics, which have been poorly surveyed to date.     

2017–2019年我国南方地区高致病性H5N6亚型禽流感病毒血凝素蛋白分子特征分析

【背景】自2014年以来,H5N6禽流感病毒在我国家禽和活禽市场持续进化,成为人类和动物健康的重大威胁。【目的】对2017–2019年中国南方地区93株高致病性H5N6禽流感病毒的HA基因进行分子进化分析。【方法】接种9–11日龄鸡胚分离核酸检测阳性的H5N6标本,运用下一代测序平台对病毒分离物进行全基因组测序,从NCBI和GISAID数据库下载参考序列,利用BLAST、MEGA6.1及Clustal X等软件进行序列分析。【结果】2017–2019年,从189份江苏省H5亚型禽类/环境标本和1名H5N6患者咽拭子标本中共分离到43株病毒,完成了33株H5N6病毒的全基因组测序。下载网上同时期中国其他地区流行的H5N6毒株序列,对总计93株H5N6病毒的HA基因进行分子进化分析。93株H5N6病毒中有78株属于Clade 2.3.4.4h,9株病毒属于Clade 2.3.4.4e,4株H5N6病毒属于Clade 2.3.4.4b,1株属于Clade 2.3.4.4f,1株属于Clade 2.3.4.4g。所有93株病毒HA蛋白的裂解位点含有多个碱性氨基酸,表明它们都属于高致病性禽流感病...     

A complete analysis of HA and NA genes of influenza A viruses

Background

More and more nucleotide sequences of type A influenza virus are available in public databases. Although these sequences have been the focus of many molecular epidemiological and phylogenetic analyses, most studies only deal with a few representative sequences. In this paper, we present a complete analysis of all Haemagglutinin (HA) and Neuraminidase (NA) gene sequences available to allow large scale analyses of the evolution and epidemiology of type A influenza.

Methodology/Principal Findings

This paper describes an analysis and complete classification of all HA and NA gene sequences available in public databases using multivariate and phylogenetic methods.

Conclusions/Significance

We analyzed 18975 HA sequences and divided them into 280 subgroups according to multivariate and phylogenetic analyses. Similarly, we divided 11362 NA sequences into 202 subgroups. Compared to previous analyses, this work is more detailed and comprehensive, especially for the bigger datasets. Therefore, it can be used to show the full and complex phylogenetic diversity and provides a framework for studying the molecular evolution and epidemiology of type A influenza virus. For more than 85% of type A influenza HA and NA sequences into GenBank, they are categorized in one unambiguous and unique group. Therefore, our results are a kind of genetic and phylogenetic annotation for influenza HA and NA sequences. In addition, sequences of swine influenza viruses come from 56 HA and 45 NA subgroups. Most of these subgroups also include viruses from other hosts indicating cross species transmission of the viruses between pigs and other hosts. Furthermore, the phylogenetic diversity of swine influenza viruses from Eurasia is greater than that of North American strains and both of them are becoming more diverse. Apart from viruses from human, pigs, birds and horses, viruses from other species show very low phylogenetic diversity. This might indicate that viruses have not become established in these species. Based on current evidence, there is no simple pattern of inter-hemisphere transmission of avian influenza viruses and it appears to happen sporadically. However, for H6 subtype avian influenza viruses, such transmissions might have happened very frequently and multiple and bidirectional transmission events might exist.     

Influenza virus sequence feature variant type analysis: evidence of a role for NS1 in influenza virus host range restriction

Genetic drift of influenza virus genomic sequences occurs through the combined effects of sequence alterations introduced by a low-fidelity polymerase and the varying selective pressures experienced as the virus migrates through different host environments. While traditional phylogenetic analysis is useful in tracking the evolutionary heritage of these viruses, the specific genetic determinants that dictate important phenotypic characteristics are often difficult to discern within the complex genetic background arising through evolution. Here we describe a novel influenza virus sequence feature variant type (Flu-SFVT) approach, made available through the public Influenza Research Database resource (www.fludb.org), in which variant types (VTs) identified in defined influenza virus protein sequence features (SFs) are used for genotype-phenotype association studies. Since SFs have been defined for all influenza virus proteins based on known structural, functional, and immune epitope recognition properties, the Flu-SFVT approach allows the rapid identification of the molecular genetic determinants of important influenza virus characteristics and their connection to underlying biological functions. We demonstrate the use of the SFVT approach to obtain statistical evidence for effects of NS1 protein sequence variations in dictating influenza virus host range restriction.     

Teaching an Old Virus New Tricks: A Review on New Approaches to Study Age-Old Questions in Influenza Biology

Influenza viruses have been studied for over 80 years, yet much about the basic viral lifecycle remain unknown. However, new imaging, biochemical, and sequencing techniques have revealed significant insight into many age-old questions of influenza virus biology. In this review, we will cover the role of imaging techniques to describe unique aspects of influenza virus assembly, biochemical techniques to study viral genomic organization, and next-generation sequencing to explore influenza genomic evolution. Our goal is to provide a brief overview of how emerging techniques are being used to answer basic questions about influenza viruses. This is not a comprehensive list of emerging techniques, rather ones that we feel will continue to make significant contributions to field of influenza biology.     

Intraseasonal Dynamics and Dominant Sequences in H3N2 Influenza

Long-term influenza evolution has been well studied, but the patterns of sequence diversity within seasons are less clear. H3N2 influenza genomes sampled from New York State over ten years indicated intraseasonal changes in evolutionary dynamics. Using the mean Hamming distance of a set of amino acid or nucleotide sequences as an indicator of its diversity, we found that influenza sequence diversity was significantly higher during the early epidemic period than later in the influenza season. Diversity was lowest during the peak of the epidemic, most likely due to the high prevalence of a single dominant amino acid sequence or very few dominant sequences during the peak epidemic period, corresponding with rapid expansion of the viral population. The frequency and duration of dominant sequences varied by influenza protein, but all proteins had an abundance of one distinct sequence during the peak epidemic period. In New York State from 1995 to 2005, high sequence diversity during the early epidemic suggested that seasonal antigenic drift could have occurred primarily in this period, followed by a clonal expansion of typically one clade during the peak of the epidemic, possibly indicating a shift to neutral drift or purifying selection.     

Sequence analysis of the 2009 pandemic influenza A H1N1 virus haemagglutinin gene from 2009-2010 Brazilian clinical samples

In this paper, we analysed the haemagglutinin (HA) gene identified by polymerase chain reaction from 90 influenza A H1N1 virus strains that circulated in Brazil from April 2009-June 2010. A World Health Organization sequencing protocol allowed us to identify amino acid mutations in the HA protein at positions S220T (71%), D239G/N/S (20%), Y247H (4.5%), E252K (3.3%), M274V (2.2%), Q310H (26.7%) and E391K (12%). A fatal outcome was associated with the D239G mutation (p < 0.0001). Brazilian HA genetic diversity, in comparison to a reference strain from California, highlights the role of influenza virus surveillance for study of viral evolution, in addition to monitoring the spread of the virus worldwide.     

H7N9亚型禽流感病毒非结构蛋白1(NS1)基因序列分析及2013上海分离株NS1的原核表达

目的:对2013年3月发生的感染人的新型H7N9亚型禽流感病毒的非结构蛋白1(NS1)基因序列进行同源性分析,构建NS1重组质粒并表达。方法:从GenBank获得2006~2013年不同来源的H7N9亚型病毒NS1序列,并进行同源性比较;利用PCR方法从H7N9亚型禽流感病毒株A/Shanghai/4664T/2013(H7N9)基因组cDNA中扩增得到全长NS1基因,并将该片段定向克隆到原核表达载体pET28a上,构建重组质粒pET28a-NS1,经酶切鉴定,将重组质粒转化大肠杆菌BL21(DE3)感受态细胞后,IPTG诱导表达,且进行Western印迹分析。结果:经序列分析,2013年暴发的H7N9型禽流感病毒的NS1基因核苷酸序列同源性为95%~100%,与之前暴发的H7N9型流感病毒NS1基因序列的同源性为86.4%~90.7%,表明2次暴发的该型流感分离株属于不同的进化分支;PCR扩增得到约680 bp的NS1基因序列,所克隆的NS1基因在原核细胞中的表达产物主要以包涵体形式存在,SDS-PAGE检测结果表明重组蛋白相对分子质量为25×103,Western印迹分析证实表达产物为H7N9禽流感病毒NS1蛋白。结论:为进一步研究H7N9亚型流感病毒NS1蛋白功能及基于NS1蛋白的抗病毒药物奠定了基础。     

Triplet entropy analysis of hemagglutinin and neuraminidase sequences measures influenza virus phylodynamics

The influenza virus has been a challenge to science due to its ability to withstand new environmental conditions. Taking into account the development of virus sequence databases, computational approaches can be helpful to understand virus behavior over time. Furthermore, they can suggest new directions to deal with influenza. This work presents triplet entropy analysis as a potential phylodynamic tool to quantify nucleotide organization of viral sequences. The application of this measure to segments of hemagglutinin (HA) and neuraminidase (NA) of H1N1 and H3N2 virus subtypes has shown some variability effects along timeline, inferring about virus evolution. Sequences were divided by year and compared for virus subtype (H1N1 and H3N2). The nonparametric Mann–Whitney test was used for comparison between groups. Results show that differentiation in entropy precedes differentiation in GC content for both groups. Considering the HA fragment, both triplet entropy as well as GC concentration show intersection in 2009, year of the recent pandemic. Some conclusions about possible flu evolutionary lines were drawn.     

Large-scale sequencing of plant genomes

The large number of ESTs generated for Arabidopsis and rice in recent years now act as an important complement to whole genome sequencing projects. The Arabidopsis Genome Initiative has begun a coordinated effort to sequence the entire genome and, as a result, increasing numbers of large sequence entries can be found in the public databases. In addition, the mitochondrial genome of Arabidopsis has been completely sequenced. Genome sequencing studies and the public sequence databases have begun to influence the direction of diverse areas of research from physiology to evolution.     

Analysis of the evolution and variation of the human influenza A virus nucleoprotein gene from 1933 to 1990.

This study examined the evolution and variation of the human influenza virus nucleoprotein gene from the earliest isolates to the present. Phylogenetic reconstruction of the most parsimonious evolutionary path connecting 49 nucleoprotein sequences yielded a single lineage. The average calculated rate of mutation was 3.6 nucleotide substitutions per year (2.3 x 10(-3) substitutions per site per year). Thirty-two percent of these mutations resulted in amino acid substitutions, and the remainder were silent mutations. Analysis of virus isolates from China and elsewhere showed no significant differences in their rate of evolution, genetic diversity, or mean survival time. The nearly constant rate of change was maintained through the two antigenic shifts, and there were no obvious changes in the number or types of mutations associated with the changes in the surface proteins. A detailed comparison of the changes that have occurred on the main evolutionary path with those that have occurred on the side branches of the phylogenetic tree was made. This showed that while 35% of the mutations on the side branches resulted in amino acid changes, only 21% of those on the main path affected the protein sequence. These results suggest that although the rate of change of the human influenza virus nucleoprotein is much higher than that previously described for avian influenza viruses, there are measurable constraints on the evolution of the surviving virus lineage. Comparison of the nucleoproteins of virus isolates adapted to chicken embryos with the nucleoproteins of those grown only in MDCK cells revealed no consistent differences between the virus pairs. Thus, although the nucleoprotein is known to be critical for host specificity, its adaptation to growth in eggs apparently involves no immediate selective pressures, such as are found with hemagglutinin.     

Analysis of the genetic diversity of influenza A viruses using next-generation DNA sequencing

Background

Influenza viruses exist as a large group of closely related viral genomes, also called quasispecies. The composition of this influenza viral quasispecies can be determined by an accurate and sensitive sequencing technique and data analysis pipeline. We compared the suitability of two benchtop next-generation sequencers for whole genome influenza A quasispecies analysis: the Illumina MiSeq sequencing-by-synthesis and the Ion Torrent PGM semiconductor sequencing technique.

Results

We first compared the accuracy and sensitivity of both sequencers using plasmid DNA and different ratios of wild type and mutant plasmid. Illumina MiSeq sequencing reads were one and a half times more accurate than those of the Ion Torrent PGM. The majority of sequencing errors were substitutions on the Illumina MiSeq and insertions and deletions, mostly in homopolymer regions, on the Ion Torrent PGM. To evaluate the suitability of the two techniques for determining the genome diversity of influenza A virus, we generated plasmid-derived PR8 virus and grew this virus in vitro. We also optimized an RT-PCR protocol to obtain uniform coverage of all eight genomic RNA segments. The sequencing reads obtained with both sequencers could successfully be assembled de novo into the segmented influenza virus genome. After mapping of the reads to the reference genome, we found that the detection limit for reliable recognition of variants in the viral genome required a frequency of 0.5% or higher. This threshold exceeds the background error rate resulting from the RT-PCR reaction and the sequencing method. Most of the variants in the PR8 virus genome were present in hemagglutinin, and these mutations were detected by both sequencers.

Conclusions

Our approach underlines the power and limitations of two commonly used next-generation sequencers for the analysis of influenza virus gene diversity. We conclude that the Illumina MiSeq platform is better suited for detecting variant sequences whereas the Ion Torrent PGM platform has a shorter turnaround time. The data analysis pipeline that we propose here will also help to standardize variant calling in small RNA genomes based on next-generation sequencing data.     

Antigenic structure of the influenza B virus hemagglutinin: nucleotide sequence analysis of antigenic variants selected with monoclonal antibodies.

We report here the complete nucleotide sequence of the hemagglutinin (HA) gene of influenza B virus B/Oregon/5/80 and, through comparative sequence analysis, identify amino acid substitutions in the HA1 polypeptide responsible for the antigenic alterations in laboratory-selected antigenic variants of this virus. The complete nucleotide sequence of the B/Oregon/5/80 HA gene was established by a combination of chemical sequencing of a full-length cDNA clone and dideoxy sequencing of the virion RNA. The nucleotide sequence is very similar to previously reported influenza B virus HA gene sequences and differs at only nine nucleotide positions from the B/Singapore/222/79 HA gene (Verhoeyen et al., Nucleic Acids Res. 11:4703-4712, 1983). The nucleotide sequences of the HA1 portions of the HA genes of 18 laboratory-selected antigenic variants were determined by the dideoxy method. Comparison of the deduced amino acid sequences of the parental and variant HA1 polypeptides revealed 16 different amino acid substitutions at nine positions. All amino acid substitutions resulted from single-point mutations, and no double mutants were detected, demonstrating that as in the influenza A viruses, single amino acid substitutions are sufficient to alter the antigenicity of the HA molecule. Many of the amino acid substitutions in the variants occurred at positions also observed to change in natural drift strains. The substitutions appear to identify at least two immunodominant regions which correspond to proposed antigenic sites A and B on the influenza A virus H3 HA.     

A Novel Nonhuman Primate Model for Influenza Transmission

Studies of influenza transmission are necessary to predict the pandemic potential of emerging influenza viruses. Currently, both ferrets and guinea pigs are used in such studies, but these species are distantly related to humans. Nonhuman primates (NHP) share a close phylogenetic relationship with humans and may provide an enhanced means to model the virological and immunological events in influenza virus transmission. Here, for the first time, it was demonstrated that a human influenza virus isolate can productively infect and be transmitted between common marmosets (Callithrix jacchus), a New World monkey species. We inoculated four marmosets with the 2009 pandemic virus A/California/07/2009 (H1N1pdm) and housed each together with a naïve cage mate. We collected bronchoalveolar lavage and nasal wash samples from all animals at regular intervals for three weeks post-inoculation to track virus replication and sequence evolution. The unadapted 2009 H1N1pdm virus replicated to high titers in all four index animals by 1 day post-infection. Infected animals seroconverted and presented human-like symptoms including sneezing, nasal discharge, labored breathing, and lung damage. Transmission occurred in one cohabitating pair. Deep sequencing detected relatively few genetic changes in H1N1pdm viruses replicating in any infected animal. Together our data suggest that human H1N1pdm viruses require little adaptation to replicate and cause disease in marmosets, and that these viruses can be transmitted between animals. Marmosets may therefore be a viable model for studying influenza virus transmission.     

Multiple genotypes of influenza B virus circulated between 1979 and 2003

The segmented genome of influenza B virus allows exchange of gene segments between cocirculating strains. Through this process of reassortment, diversity is generated by the mixing of genes between viruses that differ in one or more gene segments. Phylogenetic and evolutionary analyses of all 11 genes of 31 influenza B viruses isolated from 1979 to 2003 were used to study the evolution of whole genomes. All 11 genes diverged into two new lineages prior to 1987. All genes except the NS1 gene were undergoing linear evolution, although the rate of evolution and the degree to which nucleotide changes translated into amino acid changes varied between lineages and by gene. Frequent reassortment generated 14 different genotypes distinct from the gene constellation of viruses circulating prior to 1979. Multiple genotypes cocirculated in some locations, and a sequence of reassortment events over time could not be established. The surprising diversity of the viruses, unrestricted mixing of lineages, and lack of evidence for coevolution of gene segments do not support the hypothesis that the reassortment process is driven by selection for functional differences.     

Evolutionarily conserved protein sequences of influenza a viruses, avian and human, as vaccine targets

Background

Influenza A viruses generate an extreme genetic diversity through point mutation and gene segment exchange, resulting in many new strains that emerge from the animal reservoirs, among which was the recent highly pathogenic H5N1 virus. This genetic diversity also endows these viruses with a dynamic adaptability to their habitats, one result being the rapid selection of genomic variants that resist the immune responses of infected hosts. With the possibility of an influenza A pandemic, a critical need is a vaccine that will recognize and protect against any influenza A pathogen. One feasible approach is a vaccine containing conserved immunogenic protein sequences that represent the genotypic diversity of all current and future avian and human influenza viruses as an alternative to current vaccines that address only the known circulating virus strains.

Methodology/Principal Findings

Methodologies for large-scale analysis of the evolutionary variability of the influenza A virus proteins recorded in public databases were developed and used to elucidate the amino acid sequence diversity and conservation of 36,343 sequences of the 11 viral proteins of the recorded virus isolates of the past 30 years. Technologies were also applied to identify the conserved amino acid sequences from isolates of the past decade, and to evaluate the predicted human lymphocyte antigen (HLA) supertype-restricted class I and II T-cell epitopes of the conserved sequences. Fifty-five (55) sequences of 9 or more amino acids of the polymerases (PB2, PB1, and PA), nucleoprotein (NP), and matrix 1 (M1) proteins were completely conserved in at least 80%, many in 95 to 100%, of the avian and human influenza A virus isolates despite the marked evolutionary variability of the viruses. Almost all (50) of these conserved sequences contained putative supertype HLA class I or class II epitopes as predicted by 4 peptide-HLA binding algorithms. Additionally, data of the Immune Epitope Database (IEDB) include 29 experimentally identified HLA class I and II T-cell epitopes present in 14 of the conserved sequences.

Conclusions/Significance

This study of all reported influenza A virus protein sequences, avian and human, has identified 55 highly conserved sequences, most of which are predicted to have immune relevance as T-cell epitopes. This is a necessary first step in the design and analysis of a polyepitope, pan-influenza A vaccine. In addition to the application described herein, these technologies can be applied to other pathogens and to other therapeutic modalities designed to attack DNA, RNA, or protein sequences critical to pathogen function.