Widespread positive selection in synonymous sites of mammalian genes

Evolution of protein sequences is largely governed by purifying selection, with a small fraction of proteins evolving under positive selection. The evolution at synonymous positions in protein-coding genes is not nearly as well understood, with the extent and types of selection remaining, largely, unclear. A statistical test to identify purifying and positive selection at synonymous sites in protein-coding genes was developed. The method compares the rate of evolution at synonymous sites (Ks) to that in intron sequences of the same gene after sampling the aligned intron sequences to mimic the statistical properties of coding sequences. We detected purifying selection at synonymous sites in approximately 28% of the 1,562 analyzed orthologous genes from mouse and rat, and positive selection in approximately 12% of the genes. Thus, the fraction of genes with readily detectable positive selection at synonymous sites is much greater than the fraction of genes with comparable positive selection at nonsynonymous sites, i.e., at the level of the protein sequence. Unlike other genes, the genes with positive selection at synonymous sites showed no correlation between Ks and the rate of evolution in nonsynonymous sites (Ka), indicating that evolution of synonymous sites under positive selection is decoupled from protein evolution. The genes with purifying selection at synonymous sites showed significant anticorrelation between Ks and expression level and breadth, indicating that highly expressed genes evolve slowly. The genes with positive selection at synonymous sites showed the opposite trend, i.e., highly expressed genes had, on average, higher Ks. For the genes with positive selection at synonymous sites, a significantly lower mRNA stability is predicted compared to the genes with negative selection. Thus, mRNA destabilization could be an important factor driving positive selection in nonsynonymous sites, probably, through regulation of expression at the level of mRNA degradation and, possibly, also translation rate. So, unexpectedly, we found that positive selection at synonymous sites of mammalian genes is substantially more common than positive selection at the level of protein sequences. Positive selection at synonymous sites might act through mRNA destabilization affecting mRNA levels and translation.     

Rapid evolution of a primate sperm protein: relaxation of functional constraint or positive Darwinian selection?

Protamines are arginine-rich proteins that replace histones and bind sperm DNA during spermatogenesis in vertebrates. Previous studies have shown that protamine exons evolve faster than does the protamine intron. It has been suggested that this is a result of a relaxation of functional constraint. However, a more likely explanation is that the evolutionary rate of exons has been accelerated by positive Darwinian selection, because introns are generally believed to evolve in a neutral fashion. Therefore, we examined the possibility that positive selection has been acting on the protamine genes of three groups of placental mammals: primates (hominoids and Old World monkeys), rodents (mice, rats, and guinea pigs), and pecoran ruminants (deer and bovids). We found that the nucleotide substitution rate at nonsynonymous sites is significantly higher than the rate at synonymous and intron sites for protamine P1 of hominoids and Old World monkeys. This result suggests that positive selection has been operating on protamine P1 of these species. In contrast, no clear-cut evidence of positive selection was found for protamine P1 of ruminants and rodents or protamine P2 of primates. The agent of positive selection on primate protamine P1 remains unknown, though sperm competition is a possibility. Further investigations on the function and intraspecific polymorphism of this protein are needed in order to identify the selection agent.     

Hypermutable Non-Synonymous Sites Are under Stronger Negative Selection

Mutation rate varies greatly between nucleotide sites of the human genome and depends both on the global genomic location and the local sequence context of a site. In particular, CpG context elevates the mutation rate by an order of magnitude. Mutations also vary widely in their effect on the molecular function, phenotype, and fitness. Independence of the probability of occurrence of a new mutation''s effect has been a fundamental premise in genetics. However, highly mutable contexts may be preserved by negative selection at important sites but destroyed by mutation at sites under no selection. Thus, there may be a positive correlation between the rate of mutations at a nucleotide site and the magnitude of their effect on fitness. We studied the impact of CpG context on the rate of human–chimpanzee divergence and on intrahuman nucleotide diversity at non-synonymous coding sites. We compared nucleotides that occupy identical positions within codons of identical amino acids and only differ by being within versus outside CpG context. Nucleotides within CpG context are under a stronger negative selection, as revealed by their lower, proportionally to the mutation rate, rate of evolution and nucleotide diversity. In particular, the probability of fixation of a non-synonymous transition at a CpG site is two times lower than at a CpG site. Thus, sites with different mutation rates are not necessarily selectively equivalent. This suggests that the mutation rate may complement sequence conservation as a characteristic predictive of functional importance of nucleotide sites.     

Variation in synonymous codon use and DNA polymorphism within the Drosophila genome

A strong negative correlation between the rate of amino-acid substitution and codon usage bias in Drosophila has been attributed to interference between positive selection at nonsynonymous sites and weak selection on codon usage. To further explore this possibility we have investigated polymorphism and divergence at three kinds of sites: synonymous, nonsynonymous and intronic in relation to codon bias in D. melanogaster and D. simulans. We confirmed that protein evolution is one of the main explicative parameters for interlocus codon bias variation (r(2) approximately 40%). However, intron or synonymous diversities, which could have been expected to be good indicators of local interference [here defined as the additional increase of drift due to selection on tightly linked sites, also called 'genetic draft' by Gillespie (2000)] did not covary significantly with codon bias or with protein evolution. Concurrently, levels of polymorphism were reduced in regions of low recombination rates whereas codon bias was not. Finally, while nonsynonymous diversities were very well correlated between species, neither synonymous nor intron diversities observed in D. melanogaster were correlated with those observed in D. simulans. All together, our results suggest that the selective constraint on the protein is a stable component of gene evolution while local interference is not. The pattern of variation in genetic draft along the genome therefore seems to be instable through evolutionary times and should therefore be considered as a minor determinant of codon bias variance. We argue that selective constraints for optimal codon usage are likely to be correlated with selective constraints on the protein, both between codons within a gene, as previously suggested, and also between genes within a genome.     

Selection, recombination and demographic history in Drosophila miranda

Selection, recombination, and the demographic history of a species can all have profound effects on genomewide patterns of variability. To assess the impact of these forces in the genome of Drosophila miranda, we examine polymorphism and divergence patterns at 62 loci scattered across the genome. In accordance with recent findings in D. melanogaster, we find that noncoding DNA generally evolves more slowly than synonymous sites, that the distribution of polymorphism frequencies in noncoding DNA is significantly skewed toward rare variants relative to synonymous sites, and that long introns evolve significantly slower than short introns or synonymous sites. These observations suggest that most noncoding DNA is functionally constrained and evolving under purifying selection. However, in contrast to findings in the D. melanogaster species group, we find little evidence of adaptive evolution acting on either coding or noncoding sequences in D. miranda. Levels of linkage disequilibrium (LD) in D. miranda are comparable to those observed in D. melanogaster, but vary considerably among chromosomes. These patterns suggest a significantly lower rate of recombination on autosomes, possibly due to the presence of polymorphic autosomal inversions and/or differences in chromosome sizes. All chromosomes show significant departures from the standard neutral model, including too much heterogeneity in synonymous site polymorphism relative to divergence among loci and a general excess of rare synonymous polymorphisms. These departures from neutral equilibrium expectations are discussed in the context of nonequilibrium models of demography and selection.     

Evidence for purifying selection against synonymous mutations in mammalian exonic splicing enhancers

Silent sites in mammals have classically been assumed to be free from selective pressures. Consequently, the synonymous substitution rate (Ks) is often used as a proxy for the mutation rate. Although accumulating evidence demonstrates that the assumption is not valid, the mechanism by which selection acts remain unclear. Recent work has revealed that the presence of exonic splicing enhancers (ESEs) in coding sequence might influence synonymous evolution. ESEs are predominantly located near intron-exon junctions, which may explain the reduced single-nucleotide polymorphism (SNP) density in these regions. Here we show that synonymous sites in putative ESEs evolve more slowly than the remaining exonic sequence. Differential mutabilities of ESEs do not appear to explain this difference. We observe that substitution frequency at fourfold synonymous sites decreases as one approaches the ends of exons, consistent with the existing SNP data. This gradient is at least in part explained by ESEs being more abundant near junctions. Between-gene variation in Ks is hence partly explained by the proportion of the gene that acts as an ESE. Given the relative abundance of ESEs and the reduced rates of synonymous divergence within them, we estimate that constraints on synonymous evolution within ESEs causes the true mutation rate to be underestimated by not more than approximately 8%. We also find that Ks outside of ESEs is much lower in alternatively spliced exons than in constitutive exons, implying that other causes of selection on synonymous mutations exist. Additionally, selection on ESEs appears to affect nonsynonymous sites and may explain why amino acid usage near intron-exon junctions is nonrandom.     

The interaction between the first and last intron nucleotides in the second step of pre-mRNA splicing is independent of other conserved intron nucleotides.

Virtually all pre-mRNA introns begin with the sequence /GU and end with AG/ (where / indicates a border between an exon and an intron). We have previously shown that the G residues at the first and last positions of the yeast actin intron interact during the second step of splicing. In this work, we ask if other highly conserved intron nucleotides also take part in this /G-G/ interaction. Of special interest is the penultimate intron nucleotide (AG/), which is important for the second step of splicing and is in proximity to other conserved intron nucleotides. Therefore, we tested interactions of the penultimate intron nucleotide with the second intron nucleotide (/GU) and with the branch site nucleotide. We also tested two models that predict interactions between sets of three conserved intron nucleotides. In addition, we used random mutagenesis and genetic selection to search for interactions between nucleotides in the pre-mRNA. We find no evidence for other interactions between intron nucleotides besides the interaction between the first and last intron nucleotides.     

A theoretical method for evaluating the relative importance of positive selection and neutral drift from observed base changes

To evaluate the relative importance of positive selection and neutral drift from the nucleotide base changes observed in the homologous alignment of genes, a theoretical equation of base changes is formulated by including both the influence of selection and the base substitutions due to mutations. Under the assumption that the average rate of base substitutions estimated from synonymous changes is the ``true' mutation rate applicable at all positions, this method is applied to the vertebrate globin gene family, and evaluates the departures of base change rates from the ``true' mutation rate at the first and second codon positions as a consequence of preferential selection for the conservation of important function. In addition to the strong effect of selection on the amino acid residues in the internal region mostly common to myoglobin and hemoglobin chains, the distinctive directions of selective parameter values are seen at sites on the globin surface, distinguishing the subunit contact residues of hemoglobins from the polar residues on the surface of myoglobins. Moreover, this effect of selection distinguishing between the myoglobin and hemoglobin chain genes becomes weaker in cold-blooded vertebrates, especially in fish, strongly suggesting the possibility that the clear distinction between these globins is a result of selection out of the changes regarded as neutral ones in an ancestor of vertebrates. Thus, the present method may also serve to investigate the homology of many other proteins from the aspect of molecular evolution, mainly focusing on the evolution of their biological functions. Received: 2 January 1996 / Accepted: 20 February 1997     

Mammalian gene evolution: Nucleotide sequence divergence between mouse and rat

As a paradigm of mammalian gene evolution, the nature and extent of DNA sequence divergence between homologous protein-coding genes from mouse and rat have been investigated. The data set examined includes 363 genes totalling 411 kilobases, making this by far the largest comparison conducted between a single pair of species. Mouse and rat genes are on average 93.4% identical in nucleotide sequence and 93.9% identical in amino acid sequence. Individual genes vary substantially in the extent of nonsynonymous nucleotide substitution, as expected from protein evolution studies; here the variation is characterized. The extent of synonymous (or silent) substitution also varies considerably among genes, though the coefficient of variation is about four times smaller than for nonsynonymous substitutions. A small number of genes mapped to the X-chromosome have a slower rate of molecular evolution than average, as predicted if molecular evolution is male-driven. Base composition at silent sites varies from 33% to 95% G + C in different genes; mouse and rat homologues differ on average by only 1.7% in silent-site G + C, but it is shown that this is not necessarily due to any selective constraint on their base composition. Synonymous substitution rates and silent site base composition appear to be related (genes at intermediate G + C have on average higher rates), but the relationship is not as strong as in our earlier analyses. Rates of synonymous and nonsynonymous substitution are correlated, apparently because of an excess of substitutions involving adjacent pairs of nucleotides. Several factors suggest that synonymous codon usage in rodent genes is not subject to selection.     

Within- and between-species DNA sequence variation and the 'footprint' of natural selection

Extensive DNA data emerging from genome-sequencing projects have revitalized interest in the mechanisms of molecular evolution. Although the contribution of natural selection at the molecular level has been debated for over 30 years, the relevant data and appropriate statistical methods to address this issue have only begun to emerge. This paper will first present the predominant models of neutral, nearly neutral, and adaptive molecular evolution. Then, a method to identify the role of natural selection in molecular evolution by comparing within- and between-species DNA sequence variation will be presented. Computer simulations show that such methods are powerful for detecting even very weak selection. Examination of DNA variation data within and between Drosophila species suggests that 'silent' sites evolve under a balance between weak selection and genetic drift. Simulated data also show that sequence comparisons are a powerful method to detect adaptive protein evolution, even when selection is weak or affects a small fraction of nucleotide sites. In the Drosophila data examined, positive selection appears to be a predominant force in protein evolution.     

Molecular Population Genetics of an Electrophoretically Monomorphic Protein in the Alcohol Dehydrogenase Region of Drosophila Pseudoobscura

Nucleotide sequence data from the alcohol dehydrogenase (Adh) region of 18 isochromosomal strains of Drosophila pseudoobscura were used to determine whether the lack of amino acid polymorphism in ADH results from a low neutral mutation rate or a recent directional selection event. We estimated the neutral mutation parameter, 4Nmu, in synonymous sites for 17 subregions of Adh. The nucleotide diversity data were tested for departures from an equilibrium neutral model with two statistical tests. The Tajima test and the Hudson, Kreitman and Aguade test each failed to reject a neutral model. These results suggest that the ADH enzyme of D. pseudoobscura lacks amino acid polymorphisms because the neutral mutation rate of nonsynonymous sites is low. The neutral mutation parameter for synonymous sites is heterogeneous between domains of the Adh region. These data indicate that selective constrains on synonymous sites can vary between functional domains.     

Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene.

Several codon-based models for the evolution of protein-coding DNA sequences are developed that account for varying selection intensity among amino acid sites. The "neutral model" assumes two categories of sites at which amino acid replacements are either neutral or deleterious. The "positive-selection model" assumes an additional category of positively selected sites at which nonsynonymous substitutions occur at a higher rate than synonymous ones. This model is also used to identify target sites for positive selection. The models are applied to a data set of the V3 region of the HIV-1 envelope gene, sequenced at different years after the infection of one patient. The results provide strong support for variable selection intensity among amino acid sites The neutral model is rejected in favor of the positive-selection model, indicating the operation of positive selection in the region. Positively selected sites are found in both the V3 region and the flanking regions.     

Uncorrected nucleotide bias in mtDNA can mimic the effects of positive Darwinian selection

The relative rates of nucleotide substitution at synonymous and nonsynonymous sites within protein-coding regions have been widely used to infer the action of natural selection from comparative sequence data. It is known, however, that mutational and repair biases can affect rates of evolution at both synonymous and nonsynonymous sites. More importantly, it is also known that synonymous sites are particularly prone to the effects of nucleotide bias. This means that nucleotide biases may affect the calculated ratio of substitution rates at synonymous and nonsynonymous sites. Using a large data set of animal mitochondrial sequences, we demonstrate that this is, in fact, the case. Highly biased nucleotide sequences are characterized by significantly elevated dN/dS ratios, but only when the nucleotide frequencies are not taken into account. When the analysis is repeated taking the nucleotide frequencies at each codon position into account, such elevated ratios disappear. These results suggest that the recently reported differences in dN/dS ratios between vertebrate and invertebrate mitochondrial sequences could be explained by variations in mitochondrial nucleotide frequencies rather than the effects of positive Darwinian selection.     

Evolution of Mhc–DRB Introns: Implications for the Origin of Primates

Introns are generally believed to evolve too rapidly and too erratically to be of much use in phylogenetic reconstructions. Few phylogenetically informative intron sequences are available, however, to ascertain the validity of this supposition. In the present study the supposition was tested on the example of the mammalian class II major histocompatibility complex (Mhc) genes of the DRB family. Since the Mhc genes evolve under balancing selection and are believed to recombine or rearrange frequently, the evolution of their introns could be expected to be particularly rapid and subject to scrambling. Sequences of intron 4 and 5 DRB genes were obtained from polymerase chain reaction-amplified fragments of genomic DNA from representatives of six eutherian orders—Primates, Scandentia, Chiroptera, Dermoptera, Lagomorpha, and Insectivora. Although short stretches of the introns have indeed proved to be unalignable, the bulk of the intron sequences from all six orders, spanning >85 million years (my) of evolution, could be aligned and used in a study of the tempo and mode of intron evolution. The analysis has revealed the Mhc introns to evolve at a rate similar to that of other genes and of synonymous sites of non-Mhc genes. No evidence of homogenization or large-scale scrambling of the intron sequences could be found. The Mhc introns apparently evolve largely by point mutations and insertions/deletions. The phylogenetic signals contained in the intron sequences could be used to identify Scandentia as the sister group of Primates, to support the existence of the Archonta superorder, and to confirm the monophyly of the Chiroptera. Received: 26 October 1998 / Accepted: 21 December 1998     

The cDNA and protein sequences of mouse lactate dehydrogenase B. Molecular evolution of vertebrate lactate dehydrogenase genes A (muscle), B (heart) and C (testis)

Mouse lactate dehydrogenase-B cDNAs were isolated from cDNA libraries of macrophage (ICR strain) and thymus (F1 hybrid of C57BL/6 and CBA strains), and their nucleotide sequences determined. The lactate dehydrogenase-B cDNA insert of thymus clone mB188 consists of the protein-coding sequence (1002 nucleotides), the 5' (46 nucleotides) and 3' (190 nucleotides) non-coding regions, and poly(A) tail (19 nucleotides), while macrophage clone mB168 contains a partial lactate dehydrogenase cDNA insert from codon no. 55 to the poly(A) tail. Seven silent nucleotide substitutions at codon no. 142, 143, 186, 187, 241, 285 and 292, as well as a single nucleotide change in the 3' non-coding region, were found between these different strains of mice. The predicted sequence of 333 amino acids, excluding initiation methionine, was confirmed by sequencing and/or compositional analyses of a total of 103 (31%) amino acids from tryptic peptides of mouse lactate dehydrogenase-B protein. The nucleotide sequence of the mouse coding region for lactate dehydrogenase B shows 86% identity with that of the human isoenzyme, and only eight of the 139 nucleotide differences resulted in amino acid substitutions at residues 10, 13, 14, 17, 52, 132, 236 and 317. The rates of nucleotide substitutions at synonymous and nonsynonymous sites in the mammalian lactate dehydrogenase genes are calculated. The rates of synonymous substitutions for lactate dehydrogenase genes A (muscle) and B (heart) are considerably higher than the average rate computed from human and rodent genes. The rates of nonsynonymous substitutions for lactate dehydrogenase genes A (muscle) and B (heart), particularly the latter, are highly conservative. The rates of synonymous and nonsynonymous substitutions for the lactate dehydrogenase-C gene are about the same as the average rates for mammalian genes. A phylogenetic tree of vertebrate lactate dehydrogenase protein sequences is constructed. In agreement with the previous results, this analysis further indicates that lactate dehydrogenase-C gene branched off earlier than did lactate dehydrogenase-A and lactate dehydrogenase-B genes.     

蕨类植物rpoC1内含子缺失及其分子进化速率

rpoC1基因编码RNA聚合酶β°亚基蛋白, 在转录过程中与DNA模板结合, 与β亚基形成的β-β°亚基复合体构成RNA合成的催化中心。以rpoC1基因为研究对象, 在贝叶斯因子大于20的条件下, 用HyPhy软件位点模型检测到3个正选择位点和541个负选择位点; 用PAML软件位点模型检测到10个正选择位点, 其中3个位点的后验概率超过99%。此外, 基于最大似然法构建64种蕨类植物的系统发育树, 结合HyPhy软件分析rpoC1基因的转换率、颠换率、转换率/颠换率、同义替换率、非同义替换率以及同义替换率/非同义替换率, 探讨rpoC1基因内含子丢失与分子进化速率的关系。结果表明, rpoC1基因内含子缺失对转换率、颠换率以及非同义替换率有一定影响。     

Nucleic acid composition,codon usage,and the rate of synonymous substitution in protein-coding genes

Summary Based on the rates of synonymous substitution in 42 protein-codin gene pairs from rat and human, a correlation is shown to exist between the frequency of the nucleotides in all positions of the codon and the synonymous substitution rate. The correlation coefficients were positive for A and T and negative for C and G. This means that AT-rich genes accumulate more synonymous substitutions than GC-rich genes. Biased patterns of mutation could not account for this phenomenon. Thus, the variation in synonymous substitution rates and the resulting unequal codon usage must be the consequence of selection against A and T in synonymous positions. Most of the varition in rates of synonymous substitution can be explained by the nucleotide composition in synonymous positions. Codon-anticodon interactions, dinucleotide frequencies, and contextual factors influence neither the rates of synonymous substitution nor codon usage. Interestingly, the nucleotide in the second position of codons (always a nonsynonymous position) was found to affect the rate of synonymous substitution. This finding links the rate of nonsynonymous substitution with the synonymous rate. Consequently, highly conservative proteins are expected to be encoded by genes that evolve slowly in terms of synonymous substitutions, and are consequently highly biased in their codon usage.