The Wilhelmine E. Key 1999 Invitational lecture. Predicting the evolution of human influenza A

We studied the evolution of the HA1 domain of the H3 hemagglutinin gene from human influenza virus type A. The phylogeny of these genes showed a single dominant lineage persisting over time. We tested the hypothesis that the progenitors of this single evolutionarily successful lineage were viruses carrying mutations at codons at which prior mutations had helped the virus to avoid human immune surveillance. We found evidence that eighteen hemagglutinin codons appeared to have been under positive selection to change the amino acid they encoded in the past. Retrospective tests show that viral lineages undergoing the greatest number of mutations in the positively selected codons were the progenitors of future H3 lineages in nine of eleven recent influenza seasons. Codons under positive selection were associated with antibody combining sites A or B or the sialic acid receptor binding site. However, not all codons in these sites had predictive value. Monitoring new H3 isolates for additional changes in positively selected codons might help identify the most fit extant viral strains that arise during antigenic drift.     

Selection pressure on Haemagglutinin genes of H9N2 influenza viruses from different hosts

Positive selection and differential selective pressure analyses were carried out to study Haemagglutinin (HA) genes of H9N2 influenza viruses from different hosts in this paper. Results showed that, although most positions in HAs were under neutral or purifying evolution, a few positions located in the antigenic regions and receptor binding sites were subject to positive selection and some of them were even positively selected at the population level. In addition, there were always some positions differentially selected for viruses from different hosts. Both selection pressure working on HA codons and positions differentially selected might account for the extension of the host range and adaptations to different hosts of H9N2 influenza viruses.     

Molecular evolution of influenza A (H3N2) viruses circulated in Fujian Province, China during the 1996–2004 period

We studied the genetic and epidemic characteristics of influenza A (H3N2) viruses circulated in human in Fujian Province, south of China from 1996 to 2004. Phylogenetic analysis was carried out for genes encoding hemagglutinin1 (HA1) of influenza A virus (14 new and 11 previously reported reference se-quences). Our studies revealed that in the 8 flu seasons, the mutations of HA1 genes occurred from time to time, which were responsible for about four times of antigenic drift of influenza H3N2 viruses in Fujian, China. The data demonstrated that amino acid changes were limited to some key codons at or near antibody binding sites A through E on the HA1 molecule. The changes at the antibody binding site B or A or sialic acid receptor binding site 226 were critical for antigenic drift. But the antigenic sites might change and the key codons for antigenic drift might change as influenza viruses evolve. It seems important to monitor new H3 isolates for mutations in the positively selected codons of HA1 gene in south of Asia.     

Heterogeneous Selective Pressure Acting on Influenza B Victoria- and Yamagata-Like Hemagglutinins

As a consequence of immune pressure, influenza virus hemagglutinin presents some of its amino acids under positive selection. Several authors have reported the existence of influenza A hemagglutinin codons under positive selective pressure (PSP). In this framework, the present work objectives were to demonstrate the presence of PSP and evaluate its effects on Victoria- and Yamagata-like influenza B viruses. Methodology adopted consisted in estimating the acceptance rate of nonsynonymous substitutions (ω = d_N/d_S) that describe the strength of selective pressure and identifying codons that may be positively selected, applying a set of continuous-time Markov chain codon-substitution models. Two groups of HA1 sequences (140 from Yamagata and 60 from Victoria lineage) were used. All the model maximum-likelihood estimates were obtained using codeml software application (PAML 3.15). The hypothesis of no existence of sites under PSP was rejected for both lineages (p < 0.001), using likelihood ratio tests. These results demonstrate the presence of positive selection acting on hemagglutinin of both Yamagata- and Victoria-like influenza B viruses. Several different sites were identified to be under PSP on Yamagata and Victoria hemagglutinins. Sites found with a posterior probability > 0.95 were codons 197 and 199 in both lineages, codon 75 in the Yamagata lineage, and codon 129 in the Victoria lineage. The detected amino acids are located at or near antigenic sites in influenza A virus H3 hemagglutinin. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Positive selection operates continuously on hemagglutinin during evolution of H3N2 human influenza A virus

It has been proposed that antigenic evolution of hemagglutinin 1 (HA1) for H3N2 human influenza A virus was punctuated. In the population genetic analysis, however, it was controversial whether positive selection operated on HA1 in a punctuated manner for the branches of the phylogenetic tree where transitions to new antigenic clusters occurred (C branches), or continuously. In the molecular evolutionary analysis, positive selection was detected for the trunk (T) branches but the relationship between antigenic evolution and positive selection was unclear. Here molecular evolutionary analysis was conducted to examine natural selection operating on HA1 of H3N2 human influenza A virus by dividing HA1 into epitopes A-E and other sites, as well as dividing the phylogenetic tree into the C branches overlapping with the T branches (C-T branches), those not overlapping with the T branches (C-NT branches), the T branches not overlapping with the C branches (NC-T branches), and other branches (NC-NT branches). Positive selection was detected for C, T, and NC-T branches, whereas evolution for the NC-NT branches appeared to be mainly neutral. Positive selection appeared to have operated throughout the trunk, which covered the entire time period of the phylogenetic tree, suggesting that positive selection operated continuously on HA1 during evolution of H3N2 human influenza A virus.     

Evolutionary analyses on the HA gene ofpandemic H1N1/09: early findings

The HA protein is responsible for influenza virus attachment and the subsequent fusion of viral and cellular membranes. Antigenic drift is driven by an accumulation of point mutations in the HA. And, the receptor-binding specificity of HA is responsible for the host range restriction of the virus. In April 2009, large outbreaks of novel H1N1 influenza in human population were reported from North America. The pandemic H1N1 virus originated from swine influenza virus. Evolutionary process of the pandemic virus after its introduction to human population remains to be clarified. We conducted phylogenetic analyses constructing a phylogenetic tree for and calculating site-by-site selective pressures in the HA gene. Phylogenetic tree showed that pandemic viruses were not clustered clearly by their geographical location or isolation time in the phylogenetic tree. The virus has been circulating the globe extensively with multiple introductions into most geographical areas. We found 3 sites positively selected in the HA gene for pandemic H1N1 virus. Among them, position 206 is located in an antigenic site. We did not find significant negative selection on any of the receptor binding sites. The virus has been evolving under unique selective pressure.     

Positive Evolutionary Selection On the RIG-I-Like Receptor Genes in Mammals

The mammalian RIG-I-like receptors, RIG-I, MDA5 and LGP2, are a family of DExD/H box RNA helicases responsible for the cytoplasmic detection of viral RNA. These receptors detect a variety of RNA viruses, or DNA viruses that express unusual RNA species, many of which are responsible for a great number of severe and lethal diseases. Host innate sentinel proteins involved in pathogen recognition must rapidly evolve in a dynamic arms race with pathogens, and thus are subjected to long-term positive selection pressures to avoid potential infections. Using six codon-based Maximum Likelihood methods, we were able to identify specific codons under positive selection in each of these three genes. The highest number of positively selected codons was detected in MDA5, but a great percentage of these codons were located outside of the currently defined protein domains for MDA5, which likely reflects the imposition of both functional and structural constraints. Additionally, our results support LGP2 as being the least prone to evolutionary change, since the lowest number of codons under selection was observed for this gene. On the other hand, the preponderance of positively selected codons for RIG-I were detected in known protein functional domains, suggesting that pressure has been imposed by the vast number of viruses that are recognized by this RNA helicase. Furthermore, the RIG-I repressor domain, the region responsible for recognizing and binding to its RNA substrates, exhibited the strongest evidence of selective pressures. Branch-site analyses were performed and several species branches on the three receptor gene trees showed evidence of episodic positive selection. In conclusion, by looking for evidence of positive evolutionary selection on mammalian RIG-I-like receptor genes, we propose that a multitude of viruses have crafted the receptors biological function in host defense, specifically for the RIG-I gene, contributing to the innate species-specific resistance/susceptibility to diverse viral pathogens.     

Detection of site-specific positive Darwinian selection on pandemic influenza A/H1N1 virus genome: integrative approaches

In the twenty-first century, the first pandemic novel human influenza A/H1N1virus (NIV) outbreak was reported at Mexico and USA on March and early April, 2009 respectively. The outbreak occurred among human populations due to the presence of meager or no immune response against newly emerged viruses. The success of vaccines and drugs depends on their low susceptibility to the formation of escape mutants in virus. Identification of excess, non-synonymous substitutions over synonymous ones is a main indicator of positive Darwinian selection in protein-coding genes of NIVs. The positive Darwinian selection operating on each site of proteins were inferred by computing ω, the ratio of the non-synonymous/synonymous substitutions [dN/dS (or) K_a/K_s], which was calculated by three different methods in terms of codon-based maximum likelihood, branch-site and empirical Bayesian methods under various models. Totally, nine sites from PB2, PB1, HA, M2 and NS1 are inferred as positively selected. The function for amino acid sites of NIVs proteins under positive selection are inferred by comparing the sites with experimentally determined functionally known amino acid sites. Completely 4 positively selected sites of PB1, HA and M2 are found to be involved in B-cell epitopes (BCEs). Interestingly, most of these sites are also involving in T-cell epitopes (TCEs). However, more sites under positive selection forces are involved in TCEs than those of BCEs. Amino acid sites engaged in both BCEs and TCEs should be measured as highly suitable targets, because these sites could induce the strong humoral and cellular immune responses against targets.     

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.     

1995～2005年中国H3N2亚型人流感病毒血凝素基因变异与流行相关性研究

为了探讨中国H3N2亚型人流感病毒血凝素基因变异与流行关系,对1995～2005年中国分离的550株H3N2流感病毒的HA1序列进行分析.HA1基因进化树显示出以很长的主干和很短的侧枝为特征.HA1的氨基酸位点变异主要位于5个已知的抗原位点及其附近,同时其它位点也有改变.通过对HA1序列资料分析发现这期间导致H3N2流行的序列改变的三种可能,第一种是同时出现多位点变化,第二种是位点变化逐渐发生累积到多个位点变化,第三种是单个抗原位点和受体结合位点同时改变,均可以引起H3N2流行.     

The role of selection in the evolution of human mitochondrial genomes

High mutation rate in mammalian mitochondrial DNA generates a highly divergent pool of alleles even within species that have dispersed and expanded in size recently. Phylogenetic analysis of 277 human mitochondrial genomes revealed a significant (P < 0.01) excess of rRNA and nonsynonymous base substitutions among hotspots of recurrent mutation. Most hotspots involved transitions from guanine to adenine that, with thymine-to-cytosine transitions, illustrate the asymmetric bias in codon usage at synonymous sites on the heavy-strand DNA. The mitochondrion-encoded tRNAThr varied significantly more than any other tRNA gene. Threonine and valine codons were involved in 259 of the 414 amino acid replacements observed. The ratio of nonsynonymous changes from and to threonine and valine differed significantly (P = 0.003) between populations with neutral (22/58) and populations with significantly negative Tajima's D values (70/76), independent of their geographic location. In contrast to a recent suggestion that the excess of nonsilent mutations is characteristic of Arctic populations, implying their role in cold adaptation, we demonstrate that the surplus of nonsynonymous mutations is a general feature of the young branches of the phylogenetic tree, affecting also those that are found only in Africa. We introduce a new calibration method of the mutation rate of synonymous transitions to estimate the coalescent times of mtDNA haplogroups.     

Inference of genotype-phenotype relationships in the antigenic evolution of human influenza A (H3N2) viruses

Distinguishing mutations that determine an organism's phenotype from (near-) neutral 'hitchhikers' is a fundamental challenge in genome research, and is relevant for numerous medical and biotechnological applications. For human influenza viruses, recognizing changes in the antigenic phenotype and a strains' capability to evade pre-existing host immunity is important for the production of efficient vaccines. We have developed a method for inferring 'antigenic trees' for the major viral surface protein hemagglutinin. In the antigenic tree, antigenic weights are assigned to all tree branches, which allows us to resolve the antigenic impact of the associated amino acid changes. Our technique predicted antigenic distances with comparable accuracy to antigenic cartography. Additionally, it identified both known and novel sites, and amino acid changes with antigenic impact in the evolution of influenza A (H3N2) viruses from 1968 to 2003. The technique can also be applied for inference of 'phenotype trees' and genotype-phenotype relationships from other types of pairwise phenotype distances.     

Detecting patches of protein sites of influenza A viruses under positive selection

Influenza A viruses are single-stranded RNA viruses capable of evolving rapidly to adapt to environmental conditions. Examples include the establishment of a virus in a novel host or an adaptation to increasing immunity within the host population due to prior infection or vaccination against a circulating strain. Knowledge of the viral protein regions under positive selection is therefore crucial for surveillance. We have developed a method for detecting positively selected patches of sites on the surface of viral proteins, which we assume to be relevant for adaptive evolution. We measure positive selection based on dN/dS ratios of genetic changes inferred by considering the phylogenetic structure of the data and suggest a graph-cut algorithm to identify such regions. Our algorithm searches for dense and spatially distinct clusters of sites under positive selection on the protein surface. For the hemagglutinin protein of human influenza A viruses of the subtypes H3N2 and H1N1, our predicted sites significantly overlap with known antigenic and receptor-binding sites. From the structure and sequence data of the 2009 swine-origin influenza A/H1N1 hemagglutinin and PB2 protein, we identified regions that provide evidence of evolution under positive selection since introduction of the virus into the human population. The changes in PB2 overlap with sites reported to be associated with mammalian adaptation of the influenza A virus. Application of our technique to the protein structures of viruses of yet unknown adaptive behavior could identify further candidate regions that are important for host-virus interaction.     

The variable codons of H5N1 avian influenza A virus haemagglutinin genes

We investigated the selection pressures on the haemagglutinin genes of H5N1 avian influenza viruses using fixed effects likelihood models. We found evidence of positive selection in the sequences from isolates from 1997 to 2007, except viruses from 2000. The haemagglutinin sequences of viruses from southeast Asia, Hong Kong and mainland China were the most polymorphic and had similar nonsynonymous profiles. Some sites were positively selected in viruses from most regions and a few of these sites displayed different amino acid patterns. Selection appeared to produce different outcomes in viruses from Europe, Africa and Russia and from different host types. One position was found to be positively selected for human isolates only. Although the functions of some positively selected positions are unknown, our analysis provided evidence of different temporal, spatial and host adaptations for H5N1 avian influenza viruses.     

The variable codons of H5N1 avian influenza A virus haemagglutinin genes

We investigated the selection pressures on the haemagglutinin genes of H5N1 avian influenza viruses using fixed effects likelihood models. We found evidence of positive selection in the sequences from isolates from 1997 to 2007, except viruses from 2000. The haemagglutinin sequences of viruses from southeast Asia, Hong Kong and mainland China were the most polymorphic and had similar nonsyn-onymous profiles. Some sites were positively selected in viruses from most regions and a few of these sites displayed different amino acid patterns. Selection appeared to produce different outcomes in vi-ruses from Europe, Africa and Russia and from different host types. One position was found to be positively selected for human isolates only. Although the functions of some positively selected posi-tions are unknown, our analysis provided evidence of different temporal, spatial and host adaptations for H5N1 avian influenza viruses.     

Influenza A viruses lacking sialidase activity can undergo multiple cycles of replication in cell culture, eggs, or mice

Influenza A viruses possess both hemagglutinin (HA), which is responsible for binding to the terminal sialic acid of sialyloligosaccharides on the cell surface, and neuraminidase (NA), which contains sialidase activity that removes sialic acid from sialyloligosaccharides. Interplay between HA receptor-binding and NA receptor-destroying sialidase activity appears to be important for replication of the virus. Previous studies by others have shown that influenza A viruses lacking sialidase activity can undergo multiple cycles of replication if sialidase activity is provided exogenously. To investigate the sialidase requirement of influenza viruses further, we generated a series of sialidase-deficient mutants. Although their growth was less efficient than that of the parental NA-dependent virus, these viruses underwent multiple cycles of replication in cell culture, eggs, and mice. To understand the molecular basis of this viral growth adaptation in the absence of sialidase activity, we investigated changes in the HA receptor-binding affinity of the sialidase-deficient mutants. The results show that mutations around the HA receptor-binding pocket reduce the virus's affinity for cellular receptors, compensating for the loss of sialidase. Thus, sialidase activity is not absolutely required in the influenza A virus life cycle but appears to be necessary for efficient virus replication.     

Fast, accurate and simulation-free stochastic mapping

Mapping evolutionary trajectories of discrete traits onto phylogenies receives considerable attention in evolutionary biology. Given the trait observations at the tips of a phylogenetic tree, researchers are often interested where on the tree the trait changes its state and whether some changes are preferential in certain parts of the tree. In a model-based phylogenetic framework, such questions translate into characterizing probabilistic properties of evolutionary trajectories. Current methods of assessing these properties rely on computationally expensive simulations. In this paper, we present an efficient, simulation-free algorithm for computing two important and ubiquitous evolutionary trajectory properties. The first is the mean number of trait changes, where changes can be divided into classes of interest (e.g. synonymous/non-synonymous mutations). The mean evolutionary reward, accrued proportionally to the time a trait occupies each of its states, is the second property. To illustrate the usefulness of our results, we first employ our simulation-free stochastic mapping to execute a posterior predictive test of correlation between two evolutionary traits. We conclude by mapping synonymous and non-synonymous mutations onto branches of an HIV intrahost phylogenetic tree and comparing selection pressure on terminal and internal tree branches.     

在疫苗免疫选择压力下H9N2亚型禽流行性感冒病毒HA基因的遗传变异

为了解H9N2亚型禽流行性感冒(流感)病毒在同亚型灭活疫苗的选择压力下的遗传变异情况,对某鸡场的感染鸡群进行连续4年的跟踪监测,对使用疫苗前和持续使用疫苗后不同时段分离到的H9N2亚型禽流感病毒的HA基因进行全序列分析.结果表明,在使用第一次分离的病毒株制备的疫苗后8个月分离到的病毒株,其HA基因仅发生一个氨基酸的差异;但在继续使用该疫苗的第二个和第三个年头分离的病毒株,它们的HA基因则一直在发生较大的变化.这一发现对进一步研究禽流感病毒在不断使用疫苗的选择压力下发生变异的规律,指导制定正确的禽流感防制对策具有重要意义.     

Inferring Weak Selection from Patterns of Polymorphism and Divergence at ``silent' Sites in Drosophila DNA

Patterns of codon usage and ``silent'''' DNA divergence suggest that natural selection discriminates among synonymous codons in Drosophila. ``Preferred'''' codons are consistently found in higher frequencies within their synonymous families in Drosophila melanogaster genes. This suggests a simple model of silent DNA evolution where natural selection favors mutations from unpreferred to preferred codons (preferred changes). Changes in the opposite direction, from preferred to unpreferred synonymous codons (unpreferred changes), are selected against. Here, selection on synonymous DNA mutations is investigated by comparing the evolutionary dynamics of these two categories of silent DNA changes. Sequences from outgroups are used to determine the direction of synonymous DNA changes within and between D. melanogaster and Drosophila simulans for five genes. Population genetics theory shows that differences in the fitness effect of mutations can be inferred from the comparison of ratios of polymorphism to divergence. Unpreferred changes show a significantly higher ratio of polymorphism to divergence than preferred changes in the D. simulans lineage, confirming the action of selection at silent sites. An excess of unpreferred fixations in 28 genes suggests a relaxation of selection on synonymous mutations in D. melanogaster. Estimates of selection coefficients for synonymous mutations (3.6 <|N(e)s| < 1.3) in D. simulans are consistent with the reduced efficacy of natural selection (|N(e)s| < 1) in the three- to sixfold smaller effective population size of D. melanogaster. Synonymous DNA changes appear to be a prevalent class of weakly selected mutations in Drosophila.