Origin and evolution of influenza virus hemagglutinin genes

Influenza A, B, and C viruses are the etiological agents of influenza. Hemagglutinin (HA) is the major envelope glycoprotein of influenza A and B viruses, and hemagglutinin-esterase (HE) in influenza C viruses is a protein homologous to HA. Because influenza A virus pandemics in humans appear to occur when new subtypes of HA genes are introduced from aquatic birds that are known to be the natural reservoir of the viruses, an understanding of the origin and evolution of HA genes is of particular importance. We therefore conducted a phylogenetic analysis of HA and HE genes and showed that the influenza A and B virus HA genes diverged much earlier than the divergence between different subtypes of influenza A virus HA genes. The rate of amino acid substitution for A virus HAs from duck, a natural reservoir, was estimated to be 3.19 x 10(-4) per site per year, which was slower than that for human and swine A virus HAs but similar to that for influenza B and C virus HAs (HEs). Using this substitution rate from the duck, we estimated that the divergences between different subtypes of A virus HA genes occurred from several thousand to several hundred years ago. In particular, the earliest divergence time was estimated to be about 2,000 years ago. Also, the A virus HA gene diverged from the B virus HA gene about 4,000 years ago and from the C virus HE gene about 8,000 years ago. These time estimates are much earlier than the previous ones.     

Evolution of influenza A virus PB2 genes: implications for evolution of the ribonucleoprotein complex and origin of human influenza A virus

Phylogenetic analysis of 20 influenza A virus PB2 genes showed that PB2 genes have evolved into the following four major lineages: (i) equine/Prague/56 (EQPR56); (ii and iii) two distinct avian PB2 lineages, one containing FPV/34 and H13 gull virus strains and the other containing North American avian and recent equine strains; and (iv) human virus strains joined with classic swine virus strains (i.e., H1N1 swine virus strains related to swine/Iowa/15/30). The human virus lineage showed the greatest divergence from its root relative to other lineages. The estimated nucleotide evolutionary rate for the human PB2 lineage was 1.82 x 10(-3) changes per nucleotide per year, which is within the range of published estimates for NP and NS genes of human influenza A viruses. At the amino acid level, PB2s of human viruses have accumulated 34 amino acid changes over the past 55 years. In contrast, the avian PB2 lineages showed much less evolution, e.g., recent avian PB2s showed as few as three amino acid changes relative to the avian root. The completion of evolutionary analyses of the PB1, PB2, PA and NP genes of the ribonucleoprotein (RNP) complex permits comparison of evolutionary pathways. Different patterns of evolution among the RNP genes indicate that the genes of the complex are not coevolving as a unit. Evolution of the PB1 and PB2 genes is less correlated with host-specific factors, and their proteins appear to be evolving more slowly than NP and PA. This suggests that protein functional constraints are limiting the evolutionary divergence of PB1 and PB2 genes. The parallel host-specific evolutionary pathways of the NP and PA genes suggest that these proteins are coevolving in response to host-specific factors. PB2s of human influenza A viruses share a common ancestor with classic swine virus PB2s, and the pattern of evolution suggests that the ancestor was an avian virus PB2. This same pattern of evolution appears in the other genes of the RNP complex. Antigenic studies of HA and NA proteins and sequence comparisons of NS and M genes also suggest a close ancestry for these genes in human and classic swine viruses. From our review of the evolutionary patterns of influenza A virus genes, we propose the following hypothesis: the common ancestor to current strains of human and classic swine influenza viruses predated the 1918 human pandemic virus and was recently derived from the avian host reservoir.     

Phylodynamics and molecular evolution of influenza A virus nucleoprotein genes in Taiwan between 1979 and 2009

Background

Many studies concentrate on variation in the hemagglutinin glycoprotein (HA) because of its significance in host immune response, the evolution of this virus is even more complex when other genome segments are considered. Recently, it was found that cytotoxic T lymphocytes (CTL) play an important role in immunity against influenza and most CTL epitopes of human influenza viruses were remarkably conserved. The NP gene has evolved independently in human and avian hosts after 1918 flu pandemic and it has been assigned a putative role as a determinant of host range.

Methods and Findings

Phylodynamic patterns of the genes encoding nucleoprotein (NP) of influenza A viruses isolated from 1979–2009 were analyzed by applying the Bayesian Markov Chain Monte Carlo framework to better understand the evolutionary mechanisms of these Taiwanese isolates. Phylogenetic analysis of the NP gene showed that all available H3 worldwide isolates collected so far were genetically similar and divided into two major clades after the year 2004. We compared the deduced amino acid sequences of the NP sequences from human, avian and swine hosts to investigate the emergence of potential adaptive mutations. Overall, selective pressure on the NP gene of human influenza A viruses appeared to be dominated by purifying selection with a mean d_N/d_S ratio of 0.105. Site-selection analysis of 488 codons, however, also revealed 3 positively selected sites in addition to 139 negatively selected ones.

Conclusions

The demographic history inferred by Bayesian skyline plot showed that the effective number of infections underwent a period of smooth and steady growth from 1998 to 2001, followed by a more recent rise in the rate of spread. Further understanding the correlates of interspecies transmission of influenza A virus genes from other host reservoirs to the human population may help to elucidate the mechanisms of variability among influenza A virus.     

Evolution of influenza virus genes

The nucleotide sequences of the eight different influenza A virus segments (genes) were compared among 14 different subtypes. These comparisons demonstrate the presence of molecular clocks in the viral genes; they accumulated both silent and amino acid-changing substitutions at approximately constant rates with respect to time during evolution. In addition, comparison of the rates of evolution among the eight viral genes, excluding the P2 gene, revealed a rapid and roughly equal rate of silent substitution for different genes. The P2 gene exception is explained as the result of recombination (reassortment) between distantly related strains. The rate of amino acid-changing substitution differs greatly from gene to gene. The rate of silent substitution was estimated to be 1.1 X 10(-2)/site/year on the average--that is, about 2 X 10(6) times higher than eukaryotic gene equivalents, which is remarkable. Strain A/USSR/90/77 was shown to evolve with a rate that is similar to those of other strains but to behave as if replication was frozen during a certain period (Nakajima et al. 1978). The frozen period was estimated to be 25 yr on the basis of the molecular clock. A similar analysis revealed another example of frozen replication--in this case, apparently for a period of about 9 yr- -in a duck strain, A/duck/Ontario/77.      

The genes of influenza virus

The 5′ terminal sequences of several adenovirus 2 (Ad2) mRNAs, isolated late in infection, are complementary to sequences within the Ad2 genome which are remote from the DNA from which the main coding sequence of each mRNA is transcribed. This has been observed by forming RNA displacement loops (R loops) between Ad2 DNA and unfractionated polysomal RNA from infected cells. The 5′ terminal sequences of mRNAs in R loops, variously located between positions 36 and 92, form complex secondary hybrids with single-stranded DNA from restriction endonuclease fragments containing sequences to the left of position 36 on the Ad2 genome. The structures visualized in the electron microscope show that short sequences coded at map positions 16.6, 19.6 and 26.6 on the R strand are joined to form a leader sequence of 150–200 nucleotides at the 5′ end of many late mRNAs. A late mRNA which maps to the left of position 16.6 shows a different pattern of second site hybridization. It contains sequences from 4.9?6.0 linked directly to those from 9.6?10.9. These findings imply a new mechanism for the biosynthesis of Ad2 mRNA in mammalian cells.     

Compensatory evolution of net-charge in influenza A virus hemagglutinin

The propagation of influenza A virus depends on the balance between the activities of hemagglutinin (HA) for binding to host cells and neuraminidase (NA) for releasing from infected cells (HA-NA balance). Since the host cell membrane and the sialic acid receptor are negatively charged, the amino acid substitutions increasing (charge+) and decreasing (charge-) the positive charge of HA subunit 1 (HA1) enhance and reduce, respectively, the binding avidity and affinity. The positive charge of HA1 in human influenza A virus bearing subtype H3N2 (A/H3N2 virus) was observed to have increased during evolution, but the evolutionary mechanism for this observation was unclear because this may disrupt the HA-NA balance. Here we show, from the phylogenetic analysis of HA for human A/H3N2 and A/H1N1 viruses, that the relative frequencies of charge+ and charge- substitutions were elevated on the branches where the number of N-glycosylation sites (NGS) increased and decreased, respectively, compared to those where the number of NGS did not change. On the latter branches, the net-charge of HA1 appeared to have been largely maintained to preserve its structure and function. Since the charge+ and charge- substitutions in HA1 have opposite effects to the gain and loss of NGS on the binding and release of the virus, the net-charge of HA1 may have evolved to compensate for the effect of the gain and loss of NGS, probably through changing the avidity. Apparently, the relative frequency of charge- substitutions in HA1 of A/H3N2 virus was elevated after the introduction of oseltamivir, and that of charge+ substitutions in HA1 of A/H1N1 virus was elevated after the spread of oseltamivir resistance. These observations may also be explained by the compensatory effect of the net-charge in HA1 on the NA activity for keeping the HA-NA balance.     

流感病毒基因组进化研究进展

流感病毒先后造成了1918、1957、1968和2009年等多次全球性大流感,对人类的生命健康和社会生活形成了巨大的威胁。流感病毒的基因组进化研究为揭示病毒致病机理、疫情监测、准备疫苗和研发抗病毒药物提供了巨大的帮助。文章以流感病毒基因组进化机制为核心,结合与基因组进化密切相关的抗原性和抗药性等表型进化,对流感病毒基因组进化研究的相关进展予以介绍。     

The role of RNA folding free energy in the evolution of the polymerase genes of the influenza A virus

Background  

The influenza A virus genome is composed of eight single-stranded RNA segments of negative polarity. Although the hemagglutinin and neuraminidase genes are known to play a key role in host adaptation, the polymerase genes (which encode the polymerase segments PB2, PB1, PA) and the nucleoprotein gene are also important for the efficient propagation of the virus in the host and for its adaptation to new hosts. Current efforts to understand the host-specificity of the virus have largely focused on the amino acid differences between avian and human isolates.     

Genetic diversity and evolution of the influenza C virus

The influenza C virus is spread worldwide and causes diseases of the upper and (less frequently) lower respiratory tract in human. The virus is not pandemic, but it circulates together with pandemic influenza A and B viruses during winter months and has quite similar clinical manifestations. The influenza C virus is also encountered in animals (pigs and dogs) and is known to override the interspecific barriers of transmssion. The immune system of mammals often fails to recognize new antigenic variants of influenza C virus, which invariably arise in nature, resulting in outbreaks of diseases, although the structure of antigens in influenza C virus in general is much more stable than those of influenza viruses A and B. Variability of genetic information in natural isolates of viruses is determined by mutations, reassortment, and recombination. However, recombination events very rarely occur in genomes of negative-strand RNA viruses, including those of influenza, and virtually have no effect on their evolution. Unambiguous explanations for this phenomenon have thus far not been proposed. There is no proof of recombination processes in the influenza C virus genome. On the contrary, reassortant viruses derived from different strains of influenza C virus frequently appear in vitro and are likely to be common in nature. The genome of influenza C virus comprises seven segments. Based on the comparison of sequences in one of its genes (HEF), six genetic or antigenic lineages of this virus can be distinguished (Yamagata/26/81, Aichi/1/81, Mississippi/80, Taylor/1233/47, Sao Paulo/378/82, and Kanagawa/1/76). However, the available genetic data show that all the seven segments of the influenza C virus genome evolve independently.     

Polymorphism and evolution of collagenolytic serine protease genes in crustaceans.

Two genomic DNA fragments encoding crustacean collagenolytic serine protease genes show coding fragments that span 1522-1526 base pairs and contain seven exons encoding the complete amino acid sequence of two enzymes, CHYA and CHYB. As in serine protease genes from other organisms, the region coding for the residues around the active site is split by two introns. Although the introns differ from those of other organisms in size and nucleotide sequence, their number and location are more or less the same as found in mammalian chymotrypsin or elastase genes that evolved lately, but different for trypsin genes. Meanwhile, the junction that occurs between the propeptide and the maturation site is only found in the shrimp genes. This is also the case for the junction located 13 amino acids after the active site aspartic acid in these genes. Between 40 and 50 copies of the genes are reported by Southern analysis. Seven different genes within ChyA Pv family present 0-6% base changes, whereas five different genes belonging to ChyB Pv family show changes of up to 27% in the short studied portion of exon 4. This last family presents a mosaic organization of the coding parts, which are also expressed in the hepatopancreas of the shrimp as the variant PVC5 cDNA.     

Genome-scale evolution and phylodynamics of equine H3N8 influenza A virus

Equine influenza viruses (EIVs) of the H3N8 and H7N7 subtypes are the causative agents of an important disease of horses. While EIV H7N7 apparently is extinct, H3N8 viruses have circulated for more than 50 years. Like human influenza viruses, EIV H3N8 caused a transcontinental pandemic followed by further outbreaks and epidemics, even in populations with high vaccination coverage. Recently, EIV H3N8 jumped the species barrier to infect dogs. Despite its importance as an agent of infectious disease, the mechanisms that underpin the evolutionary and epidemiological dynamics of EIV are poorly understood, particularly at a genomic scale. To determine the evolutionary history and phylodynamics of EIV H3N8, we conducted an extensive analysis of 82 complete viral genomes sampled during a 45-year span. We show that both intra- and intersubtype reassortment have played a major role in the evolution of EIV, and we suggest that intrasubtype reassortment resulted in enhanced virulence while heterosubtypic reassortment contributed to the extinction of EIV H7N7. We also show that EIV evolves at a slower rate than other influenza viruses, even though it seems to be subject to similar immune selection pressures. However, a relatively high rate of amino acid replacement is observed in the polymerase acidic (PA) segment, with some evidence for adaptive evolution. Most notably, an analysis of viral population dynamics provided evidence for a major population bottleneck of EIV H3N8 during the 1980s, which we suggest resulted from changes in herd immunity due to an increase in vaccination coverage.     

Rohlin distance and the evolution of influenza A virus: weak attractors and precursors

The evolution of the hemagglutinin amino acids sequences of Influenza A virus is studied by a method based on an informational metrics, originally introduced by Rohlin for partitions in abstract probability spaces. This metrics does not require any previous functional or syntactic knowledge about the sequences and it is sensitive to the correlated variations in the characters disposition. Its efficiency is improved by algorithmic tools, designed to enhance the detection of the novelty and to reduce the noise of useless mutations. We focus on the USA data from 1993/94 to 2010/2011 for A/H3N2 and on USA data from 2006/07 to 2010/2011 for A/H1N1. We show that the clusterization of the distance matrix gives strong evidence to a structure of domains in the sequence space, acting as weak attractors for the evolution, in very good agreement with the epidemiological history of the virus. The structure proves very robust with respect to the variations of the clusterization parameters, and extremely coherent when restricting the observation window. The results suggest an efficient strategy in the vaccine forecast, based on the presence of "precursors" (or "buds") populating the most recent attractor.     

A robust tool highlights the influence of bird migration on influenza A virus evolution

One of the fundamental unknowns in the field of influenza biology is a panoramic understanding of the role wild birds play in the global maintenance and spread of influenza A viruses. Wild aquatic birds are considered a reservoir host for all lowly pathogenic avian influenza A viruses (AIV) and thus serve as a potential source of zoonotic AIV, such as Australasian‐origin H5N1 responsible for morbidity and mortality in both poultry and humans, as well as genes that may contribute to the emergence of pandemic viruses. Years of broad, in‐depth wild bird AIV surveillance have helped to decipher key observations and ideas regarding AIV evolution and viral ecology including the trending of viral lineages, patterns of gene flow within and between migratory flyways and the role of geographic boundaries in shaping viral evolution (Bahl et al. 2009 ; Lam et al. 2012 ). While these generally ‘virus‐centric’ studies have ultimately advanced our broader understanding of AIV dynamics, recent studies have been more host‐focused, directed at determining the potential impact of host behaviour on AIV, specifically, the influence of bird migration upon AIV maintenance and transmission. A large number of surveillance studies have taken place in Alaska, United States—a region where several global flyways overlap—with the aim of detecting the introduction of novel, Australasian‐origin highly pathogenic H5N1 AIV into North America. By targeting bird species with known migration habits, long‐distance migrators were determined to be involved in the intercontinental movement of individual AIV gene segments, but not entire viruses, between the Australasian and North American flyways (Koehler et al. 2008 ; Pearce et al. 2010 ). Yet, bird movement is not solely limited to long‐distance migration, and the relationship of resident or nonmigratory and intermediate‐distance migrant populations with AIV ecology has only recently been explored by Hill et al. ( 2012 ) in this issue of Molecular Ecology. Applying a uniquely refined, multidimensional approach, Hill et al. validate the innovative use of stable isotope assays for qualifying migration status of wild mallards within the Pacific flyway. The authors reveal that AIV prevalence and diversity did not differ in wintering mallard ducks with different migration strategies, and while migrant mallards do indeed introduce AIV, these viruses do not circulate as the predominant viruses in resident birds. On the other hand, resident mallards from more temperate regions act as reservoirs, possibly contributing to the unseasonal circulation and extended transmission period of AIV. This study highlights the impact of animal behaviour on shaping viral evolution, and the unique observations made will help inform prospective AIV surveillance efforts in wild birds.     

Rapid and reliable universal cloning of influenza A virus genes by target-primed plasmid amplification

Reverse genetics has become pivotal in influenza virus research relying on rapid generation of tailored recombinant influenza viruses. They are rescued from transfected plasmids encoding the eight influenza virus gene segments, which have been cloned using restriction endonucleases and DNA ligation. However, suitable restriction cleavage sites often are not available. Here, we describe a cloning method universal for any influenza A virus strain which is independent of restriction sites. It is based on target-primed plasmid amplification in which the insert provides two megaprimers and contains termini homologous to plasmid regions adjacent to the insertion site. For improved efficiency, a cloning vector was designed containing the negative selection marker ccdB flanked by the highly conserved influenza A virus gene termini. Using this method, we generated complete sets of functional gene segments from seven influenza A strains and three haemagglutinin genes from different serotypes amounting to 59 cloned influenza genes. These results demonstrate that this approach allows rapid and reliable cloning of any segment from any influenza A strain without any information about restriction sites. In case the PCR amplicon ends are homologous to the plasmid annealing sites only, this method is suitable for cloning of any insert with conserved termini.     

Expression of host genes in influenza virus infected cells

The NS1 protein of influenza virus shuts off host gene expression by inhibiting the polyadenylation-site cleavage of host pre-mRNAs, resulting in a general decline in cellular protein synthesis. On the other hand, an activation of several host genes related to host antiviral defense such as interferon- alpha/beta, MxA, 2',5'-oligoadenylate synthetase, and Fas occures upon infection. Therefore, balance of the shut-off and the activation of cellular genes during virus growth may be crucial in determining the outcome of infection. To obtain a comprehensive view of the global effects of influenza virus infection on human respiratory epithelial cells at the cytoplasmic mRNA level, we performed oligo DNA microarray analysis using GeneChip arrays (Affymetrix). In NCl-H292 cells infected with A/Udorn/72 virus, more than 4-fold increase of expression level was observed for 164 genes at 12 h pi. Approximately 60% of the virus-stimulated genes (VSGs) were also stimulated with interferon-beta treatment and contained the genes known to possess antiviral activity. Interestingly, majority of the VSGs were stimulated before induction of interferons, suggesting that the stimulation of the VSGs during early phase of infection is not mediated by interferons, but it is triggered from within by the virus infection.     

猪流感病毒进化方式及其流行特点

摘要：猪在甲型流感病毒生态分布和遗传进化中占有重要地位。猪的呼吸道上皮同时具有禽和人流感病毒2种类型的受体,因此人和禽流感病毒都可以感染猪,猪被认为是禽、人流感病毒的中间宿主和不同来源流感病毒的基因“混合器”。猪流感病毒（Swine influenza virus, SIV）的进化方式包括基因重配、抗原漂移和宿主适应性进化,其中基因重配是主要进化方式。与人类季节性流感病毒相比,SIV在全球的流行情况各不相同,呈地方流行性,并具有明显的地区差异。全球范围内流行的SIV亚型主要有3种：H1N1、H3N2和H1N2亚型,其中各亚型内病毒基因来源又不尽相同。欧洲、北美及我国猪群中流行的流感病毒在遗传进化和基因来源方面各具特色。目前欧洲猪群中流行的主要是类禽H1N1、H1N2和H3N2病毒,其中后两者是基因重配病毒。从1998年开始,古典猪H1N1、“人-猪-禽”三源基因重配H3N2、H1N1和H1N2病毒共存于北美的猪群中,其遗传变异日趋复杂。在基因进化上,欧洲和北美基因重配的SIV是目前新的人类大流感病毒-“甲型H1N1病毒”-的母源病毒。我国猪群中流感病毒主要是古典猪H1N1和类人H3N2病毒,但近年来在我国猪群中分离到遗传上与欧洲和北美SIV高度相关的病毒,提示我国SIV的进化趋势值得关注。1970年代以来,全球已报道了50多起人感染SIV事件,表明SIV也是一种值得重视的人兽共患病,预示了SIV可能成为人类大流感毒株或为大流感毒株提供基因。鉴于SIV在甲型流感病毒生态学上的重要意义,以及对人类公共卫生的潜在威胁,建议应尽早启动我国SIV的常规监测工作,密切关注SIV的流行动态,掌握其分子遗传进化规律。同时,将SIV的监测工作纳入整个流感病毒（人和动物流感病毒）的监测网络,在信息上实现共享,从生态学的高度把握我国流感病毒的流行和进化趋势,这对保护动物健康和预防人类大流感都有重要意义。     

The impact of seasonal and year-round transmission regimes on the evolution of influenza A virus

Punctuated antigenic change is believed to be a key element in the evolution of influenza A; clusters of antigenically similar strains predominate worldwide for several years until an antigenically distant mutant emerges and instigates a selective sweep. It is thought that a region of East-Southeast Asia with year-round transmission acts as a source of antigenic diversity for influenza A and seasonal epidemics in temperate regions make little contribution to antigenic evolution. We use a mathematical model to examine how different transmission regimes affect the evolutionary dynamics of influenza over the lifespan of an antigenic cluster. Our model indicates that, in non-seasonal regions, mutants that cause significant outbreaks appear before the peak of the wild-type epidemic. A relatively large proportion of these mutants spread globally. In seasonal regions, mutants that cause significant local outbreaks appear each year before the seasonal peak of the wild-type epidemic, but only a small proportion spread globally. The potential for global spread is strongly influenced by the intensity of non-seasonal circulation and coupling between non-seasonal and seasonal regions. Results are similar if mutations are neutral, or confer a weak to moderate antigenic advantage. However, there is a threshold antigenic advantage, depending on the non-seasonal transmission intensity, beyond which mutants can escape herd immunity in the non-seasonal region and there is a global explosion in diversity. We conclude that non-seasonal transmission regions are fundamental to the generation and maintenance of influenza diversity owing to their epidemiology. More extensive sampling of viral diversity in such regions could facilitate earlier identification of antigenically novel strains and extend the critical window for vaccine development.     

Identification of hemagglutinin and neuraminidase genes of influenza B virus.

The genome of influenza B viruses was shown by electrophoresis to consist of eight RNA segments. The fifth largest segment coded for hemagglutinin and the sixth coded for neuraminidase.     

Molecular evolution of hemagglutinin genes of H1N1 swine and human influenza A viruses

Summary The hemagglutinin (HA) genes of influenza type A (H1N1) viruses isolated from swine were cloned into plasmid vectors and their nucleotide sequences were determined. A phylogenetic tree for the HA genes of swine and human influenza viruses was constructed by the neighbor-joining method. It showed that the divergence between swine and human HA genes might have occurred around 1905. The estimated rates of synonymous (silent) substitutions for swine and human influenza viruses were almost the same. For both viruses, the rate of synonymous substitution was much higher than that of nonsynonymous (amino acid altering) substitution. It is the case even for only the antigenic sites of the HA. This feature is consistent with the neutral theory of molecular evolution. The rate of nonsynonymous substitution for human influenza viruses was three times the rate for swine influenza viruses. In particular, nonsynonymous substitutions at antigenic sites occurred less frequently in swine than in humans. The difference in the rate of nonsynonymous substitution between swine and human influenza viruses can be explained by the different degrees of functional constraint operating on the amino acid sequence of the HA in both hosts.