Natural selection on the influenza virus genome

Influenza viruses are the etiological agents of influenza. Although vaccines and drugs are available for the prophylaxis and treatment of influenza virus infections, the generation of escape mutants has been reported. To develop vaccines and drugs that are less susceptible to the generation of escape mutants, it is important to understand the evolutionary mechanisms of the viruses. Here natural selection operating on all the proteins encoded by the H3N2 human influenza A virus genome was inferred by comparing the numbers of synonymous (d(S) [D(S)]) and nonsynonymous (d(N) [D(N)]) substitutions per site. Natural selection was also inferred for the groups of functional amino acid sites involved in B-cell epitopes (BCEs), T-cell epitopes (TCEs), drug resistance, and growth in eggs. The entire region of PB1-F2 was positively selected, and positive selection also appeared to operate on BCEs, TCEs, and growth in eggs. The frequency of escape mutant generation appeared to be positively correlated with the d(N)/d(S) (D(N)/D(S)) values for the targets of vaccines and drugs, suggesting that the amino acid sites under strong functional constraint are suitable targets. In particular, TCEs may represent candidate targets because the d(N)/d(S) (D(N)/D(S)) values were small and negative selection was inferred for many of them.     

Adaptive evolution after duplication of penaeidin antimicrobial peptides

Penaeidin antimicrobial peptides in penaeid shrimps are an important component of their innate immune system that provides immunity against infection caused by several gram-positive bacteria and filamentous fungal species. Despite the knowledge on the identification and characterization of these peptides in penaeid shrimps, little is known about the evolutionary pattern of these peptides and the underlying genetic mechanisms that maintain high sequence diversities in the penaeidin gene family. Based on the phylogenetic analyses and maximum likelihood-based codon substitution analyses, here we present the convincing evidence that multiple copies of penaeidins have evolved by gene duplication, and positive Darwinian selection (adaptive evolution) is the likely cause of accelerated rate of amino acid substitutions among these duplicated genes. While the average ratio of non-synonymous to synonymous substitutions (omega) for the entire coding region of both active domains is 0.9805, few codon sites showed significantly higher omega (3.73). The likelihood ratio tests that compare models incorporating positive selection (omega>1) at certain codon sites with models not incorporating positive selection (omega<1), failed to reject (p=0) the evidence of positive Darwinian selection. The rapid adaptive evolution of this gene family might be directed by the pathogens and the faster rate of amino acid substitutions in the N-terminal proline-rich and C-terminal cysteine-rich domains could be due to their direct involvement in the protection against pathogens. When the host expose to different habitats/environment an accelerated rate of amino acid substitutions in both the active domains may also be expected.     

Codon-based tests of positive selection,branch lengths,and the evolution of mammalian immune system genes

Using basic probability theory, we show that there is a substantial likelihood that even in the presence of strong purifying selection, there will be a number of codons in which the number of synonymous nucleotide substitutions per site (d (S)) exceeds the number of non-synonymous nucleotide substitutions per site (d (N)). In an empirical study, we examined the numbers of synonymous (b (S)) and non-synonymous substitutions (b (N)) along branches of the phylogenies of 69 single-copy orthologous genes from seven species of mammals. A pattern of b (N) > b (S) was most commonly seen in the shortest branches of the tree and was associated with a high coefficient of variation in both b (N) and b (S), suggesting that high stochastic error in b (N) and b (S) on short branches, rather than positive Darwinian selection, is the explanation of most cases where b (N) is greater than b (S) on a given branch. The branch-site method of Zhang et al. (Zhang, Nielsen, Yang, Mol Biol Evol, 22:2472-2479, 2005) identified 117 codons on 35 branches as "positively selected," but a majority of these codons lacked synonymous substitutions, while in the others, synonymous and non-synonymous differences per site occurred in approximately equal frequencies. Thus, it was impossible to rule out the hypothesis that chance variation in the pattern of mutation across sites, rather than positive selection, accounted for the observed pattern. Our results showed that b (N)/b (S) was consistently elevated in immune system genes, but neither the search for branches with b (N) > b (S) nor the branch-site method revealed this trend.     

Within-host evolution of CD8+-TL epitopes encoded by overlapping and non-overlapping reading frames of simian immunodeficiency virus

In order to understand the impact of overlapping reading frames on natural selection by host CD8+ T lymphocytes (CD8(+)-TL), we analyzed the pattern of nucleotide substitution in simian immunodeficiency virus (SIV) genomes sampled from populations at time of death in 35 rhesus monkeys. Both the mean number of nonsynonymous nucleotide substitutions per nonsynonymous site (d(N)) and the mean number of synonymous nucleotide substitutions per synonymous site (d(S)) were elevated in overlap regions in comparison to non-overlap regions. Mean d(N) exceeded mean d(S) in CD8(+)-TL epitopes restricted by the host's class I major histocompatibility complex molecules. This pattern, which is indicative of positive Darwinian selection favoring amino acid changes in these epitopes, was seen in both overlap and non-overlap regions; but mean d(N) was particularly elevated in restricted CD8(+)-TL epitopes encoded in overlap regions. Amino acid changes from the inoculum were defined as parallel if the same amino acid change occurred at the same site independently in two or more monkeys, and a surprisingly high proportion (71.9%) of observed amino acid changes throughout the SIV genome occurred in parallel in different monkeys. The proportion of parallel changes in restricted epitopes encoded by overlapping reading frames was still higher (80%), supporting the hypothesis that the interaction of positive selection and overlapping reading frames enhances the probability of convergent or parallel amino acid change.     

Overdispersed molecular clock at the major histocompatibility complex loci.

The extent of amino acid differences of major histocompatibility complex molecules within species is unusually high, consistent with the finding that some pairs of alleles have persisted for more than ten million years and the view that the polymorphism has been maintained by natural selection. The disparity between synonymous and non-synonymous substitutions in the antigen recognition site, however, suggests that some non-synonymous sites have undergone a number of substitutions whereas others have little or none. To describe statistically such an overdispersed underlying process, commonly used Poisson processes are inadequate. An alternative process leads to the surprising conclusion that each non-synonymous site has accumulated as many as 2.6 substitutions, on the average, in the two lineages leading to humans and mice. The standard deviation is also very large (6.6) and the dispersion index (the ratio of the variance to the mean) is at least 17. The substitution process thus inferred qualitatively agrees with the disposition (a boomerang pattern) of substitutions between HLA-A2 and Aw68 alleles, and quantitatively agrees well with that expected where the evolution of major histocompatibility complex molecules has long been driven mostly by balancing selection.     

Reevaluation of amino acid variability of the human immunodeficiency virus type 1 gp120 envelope glycoprotein and prediction of new discontinuous epitopes

To elucidate the evolutionary mechanisms of the human immunodeficiency virus type 1 gp120 envelope glycoprotein at the single-site level, the degree of amino acid variation and the numbers of synonymous and nonsynonymous substitutions were examined in 186 nucleotide sequences for gp120 (subtype B). Analyses of amino acid variabilities showed that the level of variability was very different from site to site in both conserved (C1 to C5) and variable (V1 to V5) regions previously assigned. To examine the relative importance of positive and negative selection for each amino acid position, the numbers of synonymous and nonsynonymous substitutions that occurred at each codon position were estimated by taking phylogenetic relationships into account. Among the 414 codon positions examined, we identified 33 positions where nonsynonymous substitutions were significantly predominant. These positions where positive selection may be operating, which we call putative positive selection (PS) sites, were found not only in the variable loops but also in the conserved regions (C1 to C4). In particular, we found seven PS sites at the surface positions of the alpha-helix (positions 335 to 347 in the C3 region) in the opposite face for CD4 binding. Furthermore, two PS sites in the C2 region and four PS sites in the C4 region were detected in the same face of the protein. The PS sites found in the C2, C3, and C4 regions were separated in the amino acid sequence but close together in the three-dimensional structure. This observation suggests the existence of discontinuous epitopes in the protein's surface including this alpha-helix, although the antigenicity of this area has not been reported yet.     

Adaptive molecular evolution of virulence genes of avian influenza - A virus subtype H5N1: An analysis of host radiation

The phenomenon of host radiation is strongly influenced by the rates of mutation of their virulence genes. We have studied the molecular evolution of virulence genes (HA, NS, PB2) of the Avian Influenza Virus H5N1 from avian to human hosts. We used a site-specific comparison of synonymous (silent) and non-synonymous (amino acid altering) nucleotide substitutions for the three chosen genes in parasite populations from different hosts. Analyses were made using Maximum Likelihood (ML) genealogies for the null and alternate hypothesis based on differential gamma distribution rates. The null hypothesis had a higher rate of substitution and was found to be more suitable for all the studied genes by Likelihood Ratio Test (LRT). The study showed the NS gene to be having the fastest rate of evolution.     

藏鸡主要组织相容性复合体B-LBⅡ基因第二外显子序列遗传变异分析

根据鸡主要组织相容性复合体B-LBⅡ基因序列设计特异性引物,在藏鸡基因组中扩增了一个包括其第二外显子和第二内含子在内长度为374 bp的片段,并通过克隆和PCR直接测序获得了该片段的核苷酸序列。发现了15个B-LBⅡ新等位基因。对18个B-LBⅡ等位基因核苷酸序列和其所编码的MHCB-LBⅡ分子β1结构域的氨基酸序列分析显示,第二外显子核苷酸序列遗传多态性异常丰富,存在着62个多态变异位点(共包括80个突变),其中41个为简约性多态位点;衡量该序列遗传多样性的π值为0.0718;反映其群体内遗传变异度的平均遗传距离为0.056±0.008,低于在5个外来品种所估算的平均遗传距离。该编码区核苷酸相对异义替换率(15.61±2.69%)显著高于其同义替换率(3.25±0.94%),进一步分析表明,基因重组和平衡选择机制可能是引起B-LBⅡ基因序列变异的主要因素。在β1结构域氨基酸序列中,存在11个同义替换和27个异义替换;在24个肽结合位点中有12个变异位点;与其他6个中国地方鸡品种和一个外来品种比较发现,有11个异义氨基酸替换仅出现在藏鸡群体中,并被认为与藏鸡的免疫特异性有关,可为鸡的抗病力研究提供分子依据。     

The Site-Wise Log-Likelihood Score is a Good Predictor of Genes under Positive Selection

The strength and direction of selection on the identity of an amino acid residue in a protein is typically measured by the ratio of the rate of non-synonymous substitutions to the rate of synonymous substitutions. In attempting to predict positively selected sites from amino acid alignments, we made the unexpected observation that the site likelihood of an alignment column for a given tree tends to be negatively correlated with the posterior probability that site is in the positive selection class under widely-used codon models. This is likely because positively selected sites tend to be more variable and display more “radical” amino acid changes; both of these features are expected to result in low site log-likelihoods. We explored the efficacy of using the site log-likelihood (SLL) score as a predictor for positive selection. Through simulation we show that a SLL-based test has a low false positive rate and comparable power as the codon models. In one case where the simulated data violated the assumption that synonymous substitution rates were constant across the sites, the codon models were not able to detect positive selection in the data while the SLL test did. We applied the new method to ten empirical datasets and found that it made similar predictions as the codon models in eight of them. For the tax gene dataset the SLL test seemed to produce more reasonable results. The SLL methods are a valuable complement to codon models, especially for some cases where the assumptions of codon models are likely violated.     

New Methods for Detecting Positive Selection at Single Amino Acid Sites

Inferring positive selection at single amino acid sites is of particular importance for studying evolutionary mechanisms of a protein. For this purpose, Suzuki and Gojobori (1999) developed a method (SG method) for comparing the rates of synonymous and nonsynonymous substitutions at each codon site in a protein-coding nucleotide sequence, using ancestral codons at interior nodes of the phylogenetic tree as inferred by the maximum parsimony method. In the SG method, however, selective neutrality of nucleotide substitutions cannot be tested at codon sites, where only termination codons are inferred at any interior node or the number of equally parsimonious inferences of ancestral codons at all interior nodes exceeds 10,000. Here I present a modified SG method which is free from these problems. Specifically, I use the distance-based Bayesian method for inferring the single most likely ancestral codon from 61 sense codons at each interior node. In the computer simulation and real data analysis, the modified SG method showed a higher overall efficiency of detecting positive selection than the original SG method, particularly at highly polymorphic codon sites. These results indicate that the modified SG method is useful for inferring positive selection at codon sites where neutrality cannot be tested by the original SG method. I also discuss that the p-distance is preferable to the number of synonymous substitutions for inferring the phylogenetic tree in the SG method, and present a maximum likelihood method for detecting positive selection at single amino acid sites, which produced reasonable results in the real data analysis.     

Positively selected amino acid sites in the entire coding region of hepatitis C virus subtype 1b

To predict the amino acid sites important for the clearance of hepatitis C virus (HCV) subtype 1b in vivo, positively selected amino acid sites were detected by analyzing the sequence data collected from the international DNA databank. The rate of nonsynonymous substitutions per nonsynonymous site was compared with that of synonymous substitutions per synonymous site for each codon site in the entire coding region. As a result, 13 out of 3010 amino acid sites were found to be positively selected. Among the 13 positively selected amino acid sites, eight were located in the structural proteins and five were in the nonstructural proteins. Moreover, eight were located in B-cell epitopes and two were in T-cell epitopes. These observations suggest that both the antibody and the cytotoxic T lymphocyte are involved in the clearance of HCV subtype 1b in vivo. These positively selected amino acid sites represent candidate vaccination targets for HCV subtype 1b.     

Codon-substitution models for heterogeneous selection pressure at amino acid sites

Comparison of relative fixation rates of synonymous (silent) and nonsynonymous (amino acid-altering) mutations provides a means for understanding the mechanisms of molecular sequence evolution. The nonsynonymous/synonymous rate ratio (omega = d(N)d(S)) is an important indicator of selective pressure at the protein level, with omega = 1 meaning neutral mutations, omega < 1 purifying selection, and omega > 1 diversifying positive selection. Amino acid sites in a protein are expected to be under different selective pressures and have different underlying omega ratios. We develop models that account for heterogeneous omega ratios among amino acid sites and apply them to phylogenetic analyses of protein-coding DNA sequences. These models are useful for testing for adaptive molecular evolution and identifying amino acid sites under diversifying selection. Ten data sets of genes from nuclear, mitochondrial, and viral genomes are analyzed to estimate the distributions of omega among sites. In all data sets analyzed, the selective pressure indicated by the omega ratio is found to be highly heterogeneous among sites. Previously unsuspected Darwinian selection is detected in several genes in which the average omega ratio across sites is <1, but in which some sites are clearly under diversifying selection with omega > 1. Genes undergoing positive selection include the beta-globin gene from vertebrates, mitochondrial protein-coding genes from hominoids, the hemagglutinin (HA) gene from human influenza virus A, and HIV-1 env, vif, and pol genes. Tests for the presence of positively selected sites and their subsequent identification appear quite robust to the specific distributional form assumed for omega and can be achieved using any of several models we implement. However, we encountered difficulties in estimating the precise distribution of omega among sites from real data sets.     

Positive selection in AvrP4 avirulence gene homologues across the genus Melampsora

Pathogen genes involved in interactions with their plant hosts are expected to evolve under positive Darwinian selection or balancing selection. In this study a single copy avirulence gene, AvrP4, in the plant pathogen Melampsora lini, was used to investigate the evolution of such a gene across species. Partial translation elongation factor 1-alpha sequences were obtained to establish phylogenetic relationships among the Melampsora species. We amplified AvrP4 homologues from species pathogenic on hosts from different plant families and orders, across the inferred phylogeny. Translations of the AvrP4 sequences revealed a predicted signal peptide and towards the C-terminus of the protein, six identically spaced cysteines were identified in all sequences. Maximum likelihood analysis of synonymous versus non-synonymous substitution rates indicated that positive selection played a role in the evolution of the gene during the diversification of the genus. Fourteen codons under significant positive selection reside in the C-terminal 28 amino acid region, suggesting that this region interacts with host molecules in most sequenced accessions. Selection pressures on the gene may be either due to the pathogenicity or avirulence function of the gene or both.     

Natural selection in the <Emphasis Type="Italic">TLR-</Emphasis>related genes in the course of primate evolution

The innate immune system constitutes the front line of host defense against pathogens. Toll-like receptors (TLRs) recognize molecules derived from pathogens and play crucial roles in the innate immune system. Here, we provide evidence that the TLR-related genes have come under natural selection pressure in the course of primate evolution. We compared the nucleotide sequences of 16 TLR-related genes, including TLRs (TLR1–10), MYD88, TILAP, TICAM1, TICAM2, MD2, and CD14, among seven primate species. Analysis of the non-synonymous/synonymous substitution ratio revealed the presence of both strictly conserved and rapidly evolving regions in the TLR-related genes. The genomic segments encoding the intracellular Toll/interleukin 1 receptor domains, which exhibited lower rates of non-synonymous substitution, have undergone purifying selection. In contrast, TLR4, which carried a high proportion of non-synonymous substitutions in the part of extracellular domain spanning 200 amino acids, was found to have been the suggestive target of positive Darwinian selection in primate evolution. However, sequence analyses from 25 primate species, including eight hominoids, six Old World monkeys, eight New World monkeys, and three prosimians, showed no evidence that the pressure of positive Darwinian selection has shaped the pattern of sequence variations in TLR4 among New World monkeys and prosimians. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Further Examples of Evolution by Gene Duplication Revealed through DNA Sequence Comparisons

To test the theory that evolution by gene duplication occurs as a result of positive Darwinian selection that accompanies the acceleration of mutant substitutions, DNA sequences of recent duplication were analyzed by estimating the numbers of synonymous and nonsynonymous substitutions. For the troponin C family, at the period of differentiation of the fast and slow isoforms, amino acid substitutions were shown to have been accelerated relative to synonymous substitutions. Comparison of the first exon of α-actin genes revealed that amino acid substitutions were accelerated when the smooth muscle, skeletal and cardiac isoforms differentiated. Analysis of members of the heat shock protein 70 gene family of mammals indicates that heat shock responsive genes including duplicated copies are evolving rapidly, contrary to the cognitive genes which have been evolutionarily conservative. For the α(1)-antitrypsin reactive center, the acceleration of amino acid substitution has been found for gene pairs of recent duplication.     

Positive selection for gains of N-linked glycosylation sites in hemagglutinin during evolution of H3N2 human influenza A virus

The number of N-linked glycosylation sites in the globular head of hemagglutinin (HA) has increased during evolution of H3N2 human influenza A virus. Here natural selection operating on the gains of N-linked glycosylation sites was examined by using the single-site analysis and the single-substitution analysis. In the single-site analysis, positive selection was not inferred at the amino acid sites where the substitutions generating N-linked glycosylation sites were observed, but was detected at antigenic sites. In contrast, in the single-substitution analysis, positive selection was detected for the amino acid substitutions generating N-linked glycosylation sites. The single-site analysis and the single-substitution analysis appeared to be suitable for detecting recurrent and episodic natural selection, respectively. The gains of N-linked glycosylation sites were likely to be positively selected for the function of shielding antigenic sites from immune responses. At the antigenic sites, positive selection appeared to have operated not only on the radical substitution but also on the conservative substitution in terms of the charge of amino acids, suggesting that the antigenic drift is not a by-product of the evolution of receptor binding avidity in HA of human H3N2 virus.     

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.     

Dengue Virus 1 in Buenos Aires from 1999 to 2010: Towards Local Spread

Dengue virus (DENV) is a public health problem representing the most important arthropod-borne viral disease in humans. In Argentina, Northern provinces have reported autochthonous cases since 1997, though these outbreaks have originated in bordering countries, where co-circulation of more than one serotype has been reported. In the last decade, imported dengue cases have been reported in Buenos Aires, the urban area of Argentina with the highest population density. In 2009, a dengue outbreak affected Buenos Aires and, for the first time, local transmission was detected. All cases of this outbreak were caused by DENV-1. In this report, we present the full-length sequences of 27 DENV-1 isolates, corresponding to imported cases of 1999–2000, as well as local and imported cases of the 2009 and 2010 outbreaks. We analyzed their phylogenetic and phylodynamic relationships and their global and local spread. Additionally, we characterized their genomic and phenotypic features. All cases belonged to DENV-1 genotype V. The most recent ancestor for this genotype was dated ∼1934, whereas that for the 2009 outbreak was dated ∼2007. The mean rates of nucleotide substitution were 4.98E-4 and 8.53E-4 subs./site/yr, respectively. We inferred an introduction from Paraguay in 1999–2000 and mainly from Venezuela during 2009–2010. Overall, the number of synonymous substitutions per synonymous site significantly exceeded the number of non-synonymous substitutions per site and 12 positively selected sites were detected. These analyses could contribute to a better understanding regarding spread and evolution of this pathogen in the Southern Cone of South America.     

Maximum Likelihood Estimation on Large Phylogenies and Analysis of Adaptive Evolution in Human Influenza Virus A

Algorithmic details to obtain maximum likelihood estimates of parameters on a large phylogeny are discussed. On a large tree, an efficient approach is to optimize branch lengths one at a time while updating parameters in the substitution model simultaneously. Codon substitution models that allow for variable nonsynonymous/synonymous rate ratios (ω=d _N/d _S) among sites are used to analyze a data set of human influenza virus type A hemagglutinin (HA) genes. The data set has 349 sequences. Methods for obtaining approximate estimates of branch lengths for codon models are explored, and the estimates are used to test for positive selection and to identify sites under selection. Compared with results obtained from the exact method estimating all parameters by maximum likelihood, the approximate methods produced reliable results. The analysis identified a number of sites in the viral gene under diversifying Darwinian selection and demonstrated the importance of including many sequences in the data in detecting positive selection at individual sites. Received: 25 April 2000 / Accepted: 24 July 2000