A new method for estimating synonymous and nonsynonymous rates of nucleotide substitution considering the relative likelihood of nucleotide and codon changes

A new method is proposed for estimating the number of synonymous and nonsynonymous nucleotide substitutions between homologous genes. In this method, a nucleotide site is classified as nondegenerate, twofold degenerate, or fourfold degenerate, depending on how often nucleotide substitutions will result in amino acid replacement; nucleotide changes are classified as either transitional or transversional, and changes between codons are assumed to occur with different probabilities, which are determined by their relative frequencies among more than 3,000 changes in mammalian genes. The method is applied to a large number of mammalian genes. The rate of nonsynonymous substitution is extremely variable among genes; it ranges from 0.004 X 10(-9) (histone H4) to 2.80 X 10(-9) (interferon gamma), with a mean of 0.88 X 10(-9) substitutions per nonsynonymous site per year. The rate of synonymous substitution is also variable among genes; the highest rate is three to four times higher than the lowest one, with a mean of 4.7 X 10(-9) substitutions per synonymous site per year. The rate of nucleotide substitution is lowest at nondegenerate sites (the average being 0.94 X 10(-9), intermediate at twofold degenerate sites (2.26 X 10(-9)). and highest at fourfold degenerate sites (4.2 X 10(-9)). The implication of our results for the mechanisms of DNA evolution and that of the relative likelihood of codon interchanges in parsimonious phylogenetic reconstruction are discussed.     

Low Nucleotide Diversity in Man

The nucleotide diversity (pi) in humans is studied by using published cDNA and genomic sequences that have been carefully checked for sequencing accuracy. This measure of genetic variability is defined as the number of nucleotide differences per site between two randomly chosen sequences from a population. A total of more than 75,000 base pairs from 49 loci are compared. The DNA regions studied are the 5' and 3' untranslated regions and the amino acid coding regions. The coding regions are divided into nondegenerate sites (i.e., sites at which all possible changes are nonsynonymous), twofold degenerate sites (i.e., sites at each of which one of the three possible changes is synonymous) and fourfold degenerate sites (i.e., sites at which all three possible changes are synonymous). The pi values estimated are, respectively, 0.03 and 0.04% for the 5' and 3' UT regions, and 0.03, 0.06 and 0.11% for nondegenerate, twofold degenerate and fourfold degenerate sites. Since the highest pi value is only 0.11%, which is about one order of magnitude lower than those in Drosophila populations, the nucleotide diversity in humans is very low. The low diversity is probably due to a relatively small long-term effective population size rather than any severe bottleneck during human evolution.     

Low occurrence of gene transposition events during the evolution of the genus Drosophila

Abstract.— The role played by gene transpositions during the evolution of eukaryotic genomes is still poorly understood and indeed has been analyzed in detail only in nematodes. In Drosophila , a limited number of transpositions have been detected by comparing the chromosomal location of genes between different species. The relative importance of gene transposition versus other types of chromosomal rearrangements, for example, inversions, has not yet been evaluated. Here, we use physical mapping to perform an extensive search for long-distance gene transpositions and assess their impact during the evolution of the Drosophila genome. We compare the relative order of 297 molecular markers that cover 60% of the euchromatic fraction of the genome between two related Drosophila species and conclude that the frequency of gene transpositions is very low, namely one order of magnitude lower than that of nematodes. In addition, gene transpositions seem to be events almost exclusively associated with genes of repetitive nature such as the Histone gene complex ( HIS-C ).     

Nonrandom frequency patterns of synonymous substitutions in homologous mammalian genes

All 69 homologous coding sequences that are currently available in four mammalian orders were aligned and the synonymous positions of quartet and duet (fourfold and twofold degenerate) codons were divided into three classes (that will be called conserved, intermediate, and variable) according to whether they show no change, one change, or more than one change, respectively. We observed (1) that the frequencies of conserved, intermediate, and variable positions of quartet and duet codons are different in different genes; (2) that the frequencies of the three classes are significantly different from expectations based on a random substitution process in the majority of genes (especially for GC-rich genes) for quartet codons and in a minority of genes for doublet codons; and (3) that the frequencies of the three classes of positions of quartet codons are correlated with those of duet codons, the conserved positions of quartet and duet codons being, in addition, correlated with the degree of amino acid conservation. Our main conclusions are that synonymous substitution frequencies: (1) are gene-specific; (2) are not simply the result of a stochastic process in which nucleotide substitutions accumulate at random, over time; and (3) are correlated in quartet and duet codons.     

Contrasting levels of DNA polymorphism at the autosomal and X-linked visual color pigment loci in humans and squirrel monkeys

The X-linked color pigment (opsin) locus is known to be highly polymorphic in the squirrel monkey and other New World monkeys. To see whether this is also the case for the autosomal (blue) opsin locus, we obtained 32 squirrel monkey and 30 human blue opsin gene sequences. No amino acid polymorphism was found in either the squirrel monkey sample or the human sample, contrary to the situation at the X-linked opsin locus. This sharp contrast in the level of polymorphism might be due to differences in gene expression between the autosomal and the X-linked loci. At the X-linked locus, heterozygote advantage can occur because, owing to X-inactivation, the two alleles in a heterozygote are expressed in different cone cells, producing two types of cone cell, whereas at the autosomal locus, heterozygote advantage cannot occur because the two alleles in a heterozygote are expressed in the same cone cells, producing only one type of cone cell (i.e., phenotypically a homozygote). From the sequence data, the levels of nucleotide diversity (pi, i.e., the number of nucleotide differences per site) are estimated: for the human sample, pi = 0.00% per nondegenerate site, 0.00% per twofold degenerate site, and 0.04% per fourfold degenerate site in the coding regions and 0.01% per site in intron 4; for the squirrel monkey sample, pi = 0.00% per nondegenerate site, 0.00% per twofold degenerate site, and 0.15% per fourfold degenerate site in the coding regions and 0.17% per site in intron 4. The blue opsin genes from the common and pygmy chimpanzees, the gorilla, the capuchin, and the howler monkey were also sequenced. Features critical to the function of the opsin are well conserved in all known mammalian sequences. However, the interhelical loops are, on average, actually more conservative than the transmembrane helical regions. In addition, these sequence data and those from some other genes indicate that the common and pygmy chimpanzees are not closely related, their divergence data being from one third to one half the date of the human-chimpanzee divergence.      

Transposition of mobile dispersed genes (MDG) after substitution of various chromosome pairs in inbred Drosophila melanogaster strains

Transpositions of MDG-1, MDG-3 and copia were detected as a result of crosses of the inbred maladaptive LA stock with laboratory stocks, in order to construct the genomes carrying different combinations of the LA or non-La chromosomal pairs. Changes of the mobile gene distributions were revealed in chromosomes of hybrid genotypes, as compared to parental chromosomal pairs. A trivial source of variability of chromosomal molecular structure ensured by crossing over was excluded by inversions which serve as suppressors of crossing over in corresponding crosses. Multiple transpositions of mobile genes in definite chromosomal sites were detected in genotypes carrying chromosomal pair 2 originated from the LA stock. No such transpositions were observed, when the pair 2 was substituted by the chromosome 2 originated from the Swedish-b line or in control crosses, where the LA stock was not involved. Both LA chromosomes 2 and 3 were shown to be the targets of transpositions. Comparison of hot spot transposition sites of MDG-1, as a result of crosses, with the earlier described rare events of spontaneous transpositions in the LA stock, coupled with its fitness increase, revealed that the hot spot sites were shared in both series of experiments. The data obtained show that transpositions of mobile genetic elements may change the genetic and molecular structure of the chromosome involved in crosses, in spite of suppression of crossing over by inversions usually suggested as a tool for keeping chromosomal genetic structure intact.     

Rates of Transition and Transversion in Coding Sequences since the Human-Rodent Divergence

Protein-coding sequences of 337 human genes were compared with those of homologous genes from rodent (mouse or rat). A composite alignment containing 477,189 nucleotide positions was constructed, and 21,570 amino acid replacements were inferred. The rates of transitional and transversional silent substitutions in fourfold degenerate sites are estimated as 1.71 × 10^-9 and 1.22 × 10^-9 site^-1 year^-1, respectively. Rates of substitutions in replacement sites, subject to selective constraints mediated by the genetic code, are lower, but also reflect a transitional bias. The amino acid exchange rejected least often during evolution is Asp/Glu, which is fixed at 30% the rate of transversions in silent sites. The most mutable amino acids in this survey are threonine and serine; serine coded by AGY is more mutable than serine coded by TCN. A scoring matrix for evaluating amino acid similarity was derived from this study.     

Unbiased estimation of the rates of synonymous and nonsynonymous substitution

Summary The current convention in estimating the number of substitutions per synonymous site (K _S ) and per nonsynonymous site (K _A ) between two protein-coding genes is to count each twofold degenerate site as one-third synonymous and two-thirds nonsynonymous because one of the three possible changes at such a site is synonymous and the other two are nonsynonymous. This counting rule can considerably overestimate theK _S value because transitional mutations tend to occur more often than transversional mutations and because most transitional mutations at twofold degenerate sites are synonymous. A new method that gives unbiased estimates is proposed. An application of the new and the old method to 14 pairs of mouse and rat genes shows that the new method gives aK _S value very close to the number of substitutions per fourfold degenerate site whereas the old method gives a value 30% higher. Both methods give aK _A value close to the number of substitutions per nondegenerate site.     

Accelerated Evolution of Cytochrome b in Simian Primates: Adaptive Evolution in Concert with Other Mitochondrial Proteins?

We have sequenced the cytochrome b gene of Horsfield's tarsier, Tarsius bancanus, to complete a data set of sequences for this gene from representatives of each primate infraorder. These primate cytochrome b sequences were combined with those from representatives of three other mammalian orders (cat, whale, and rat) in an analysis of relative evolutionary rates. The nonsynonymous nucleotide substitution rate of the cytochrome b gene has increased approximately twofold along lineages leading to simian primates compared to that of the tarsier and other primate and nonprimate mammalian species. However, the rate of transversional substitutions at fourfold degenerate sites has remained uniform among all lineages. This increase in the evolutionary rate of cytochrome b is similar in character and magnitude to that described previously for the cytochrome c oxidase subunit II gene. We propose that the evolutionary rate increase observed for cytochrome b and cytochrome c oxidase subunit II may underlie an episode of coadaptive evolution of these two proteins in the mitochondria of simian primates. Received: 15 December 1997 / Accepted: 24 February 1998     

Three patterns of mitochondrial DNA nucleotide divergence in the meadow vole,Microtus pennsylvanicus

Summary The DNA sequence was determined for the cytochrome c oxidase II (COII), tRNA^Lys, and ATPase 8 genes from the mitochondrial genome of the meadow vole, Microtus pennsylvanicus. When compared to other rodents, three different patterns of evolutionary divergence were found. Nucleotide variation in tRNA^Lys is concentrated in the TC loop. Nucleotide variation in the COII gene in three genera of rodents (Microtus, Mus, Rattus) consists predominantly of transitions in the third base positions of codons. The predicted amino acid sequence in highly conserved (>92% similarity). Analysis of the ATPase 8 gene among four genera (Microtus, Cricetulus, Mus, Rattus) revealed more detectable transversions than transitions, many fixed first and second position mutations, and considerable amino acid divergence. The rate of nucleotide substitution at nonsynonymous sites in the ATPase 8 gene is 10 times the rate in the COII gene. In contrast, the estimated absolute mutation rate as determined by analysis of nucleotide substitutions at fourfold degenerate sites probably is the same for the two genes. The primary sequences of the ATPase 8 and COII peptides are constrained differently, but each peptide is conserved in terms of predicted secondary-level configuration.     

The Effects of Mutation and Natural Selection on Codon Bias in the Genes of Drosophila

Codon bias varies widely among the loci of Drosophila melanogaster, and some of this diversity has been explained by variation in the strength of natural selection. A study of correlations between intron and coding region base composition shows that variation in mutation pattern also contributes to codon bias variation. This finding is corroborated by an analysis of variance (ANOVA), which shows a tendency for introns from the same gene to be similar in base composition. The strength of base composition correlations between introns and codon third positions is greater for genes with low codon bias than for genes with high codon bias. This pattern can be explained by an overwhelming effect of natural selection, relative to mutation, in highly biased loci. In particular, this correlation is absent when examining fourfold degenerate sites of highly biased genes. In general, it appears that selection acts more strongly in choosing among fourfold degenerate codons than among twofold degenerate codons. Although the results indicate regional variation in mutational bias, no evidence is found for large scale regions of compositional homogeneity.     

Rates of DNA change and phylogeny from the DNA sequences of the alcohol dehydrogenase gene for five closely related species of Hawaiian Drosophila

The sequence of 1.6 kb of DNA surrounding the alcohol dehydrogenase (Adh) gene from five species of the Planitibia subgroup of the Hawaiian picture-winged Drosophila, with estimated divergence times of 0.4-5.1 Myr, has been determined. The gene trees which were found by using the sequence divergence from different regions of the sequences are generally in accord with the phylogeny proposed for these species when chromosomal inversions and island of origin are used. One of the species (D. picticornis) appears to be more distant from the other species in this group than they are from a member of the Grimshawi group (D. affinidisjuncta) which is chromosomally more distant. Two of the species (D. differens and D. plantibia) show heterogeneity in the nucleotide changes in the Adh coding region, heterogeneity which is interpreted to be due to a gene conversion or recombination after hybridization between the two species. The minimal rate of nucleotide substitution of synonymous nucleotides and of nontranscribed nucleotides downstream from the coding region is estimated as 1.5 x 10(-8) and 1.1 x 10(-8) substitutions/nucleotide/year, respectively. This rate is two to three times the maximal rate estimated for mammalian synonymous substitutions.     

Distance-Based Genome Rearrangement Phylogeny

Evolution operates on whole genomes through direct rearrangements of genes, such as inversions, transpositions, and inverted transpositions, as well as through operations, such as duplications, losses, and transfers, that also affect the gene content of the genomes. Because these events are rare relative to nucleotide substitutions, gene order data offer the possibility of resolving ancient branches in the tree of life; the combination of gene order data with sequence data also has the potential to provide more robust phylogenetic reconstructions, since each can elucidate evolution at different time scales. Distance corrections greatly improve the accuracy of phylogeny reconstructions from DNA sequences, enabling distance-based methods to approach the accuracy of the more elaborate methods based on parsimony or likelihood at a fraction of the computational cost. This paper focuses on developing distance correction methods for phylogeny reconstruction from whole genomes. The main question we investigate is how to estimate evolutionary histories from whole genomes with equal gene content, and we present a technique, the empirically derived estimator (EDE), that we have developed for this purpose. We study the use of EDE on whole genomes with identical gene content, and we explore the accuracy of phylogenies inferred using EDE with the neighbor joining and minimum evolution methods under a wide range of model conditions. Our study shows that tree reconstruction under these two methods is much more accurate when based on EDE distances than when based on other distances previously suggested for whole genomes. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Martin Kreitman]     

Speciation in bacteria: comparison of the 16S rRNA gene for closely related Enterococcus species

Sequencing of the 16S rRNA genes from enterococcal strains used as starters suggested the existence of specialized taxa of lactic acid enterococci within the species Enterococcus durans and E. faecium and a new species, E. lactis. Comparisons showed that the 16S rRNA genes of closely related species have the same sets of variable positions with different combinations of nucleotides. The presence of identical combinations of nucleotide substitutions in different species was assumed to result from a transfer of genetic information via gene conversion between different rRNA operons. Such events were presumably associated with speciation in bacteria.     

DNA序列进化过程中核苷酸替代的非独立性研究

本文评述了DNA序列间核苷酸替代数的估计方法,并通过对七个物种中组蛋白基因的比较对DNA进化的模型进行了考察。发现H2A基因第三位点上的碱基组成在物种间变异很大,并且跟H2A基因第一位点、H4基因第一、三位点及H2A上游,下游序列中的碱基组成有强正相关,提示DNA序列进化过程中存在着物种特异的区域性约束力。可能的原因是高等真核生物中GC含量升高,或者是染色体重组使这些同源序列位于不同的等质区段,从而受到不同的选择突变压。密码内各位点上核苷酸替代的相关性分析表明不同位点的替代是非独立的,其原因可能是一次替代事件引起多个位点的变化。文中讨论了这些结果对进化树推断的意义。     

Molecular evolution and mosaicism of leptospiral outer membrane proteins involves horizontal DNA transfer

Leptospires belong to a genus of parasitic bacterial spirochetes that have adapted to a broad range of mammalian hosts. Mechanisms of leptospiral molecular evolution were explored by sequence analysis of four genes shared by 38 strains belonging to the core group of pathogenic Leptospira species: L. interrogans, L. kirschneri, L. noguchii, L. borgpetersenii, L. santarosai, and L. weilii. The 16S rRNA and lipL32 genes were highly conserved, and the lipL41 and ompL1 genes were significantly more variable. Synonymous substitutions are distributed throughout the ompL1 gene, whereas nonsynonymous substitutions are clustered in four variable regions encoding surface loops. While phylogenetic trees for the 16S, lipL32, and lipL41 genes were relatively stable, 8 of 38 (20%) ompL1 sequences had mosaic compositions consistent with horizontal transfer of DNA between related bacterial species. A novel Bayesian multiple change point model was used to identify the most likely sites of recombination and to determine the phylogenetic relatedness of the segments of the mosaic ompL1 genes. Segments of the mosaic ompL1 genes encoding two of the surface-exposed loops were likely acquired by horizontal transfer from a peregrine allele of unknown ancestry. Identification of the most likely sites of recombination with the Bayesian multiple change point model, an approach which has not previously been applied to prokaryotic gene sequence analysis, serves as a model for future studies of recombination in molecular evolution of genes.     

On transition bias in mitochondrial genes of pocket gophers

The relative contribution of mutation and purifying selection to transition bias has not been quantitatively assessed in mitochondrial protein genes. The observed transition/transversion (s/v) ratio is (μ _s P _s)/(μ _v P _v), where μ _s and μ _v denote mutation rate of transitions and transversions, respectively, andP _s andP _v denote fixation probabilities of transitions and transversions, respectively. Because selection against synonymous transitions can be assumed to be roughly equal to that against synonymous transversions,P _s/P_v ≈ 1 at fourfold degenerate sites, so that thes/v ratio at fourfold degenerate sites is approximately μ _s /μ _v , which is a measure of mutational contribution to transition bias. Similarly, thes/v ratio at nondegenerate sites is also an estimate of μ _s /μ _v if we assume that selection against nonsynonymous transitions is roughly equal to that against nonsynonymous transversions. In two mitochondrial genes, cytochrome oxidase subunit I (COI) and cytochromeb (cyt-b) in pocket gophers, thes/v ratio is about two at nondegenerate and fourfold degenerate sites for both the COI and the cyt-b genes. This implies that mutation contribution to transition bias is relatively small. In contrast, thes/v ratio is much greater at twofold degenerate sites, being 48 for COI and 40 for cyt-b. Given that the μ _s /μ _v ratio is about 2, theP _s/P_v ratio at twofold degenerate sites must be on the order of 20 or greater. This suggests a great effect of purifying selection on transition bias in mitochondrial protein genes because transitions are synonymous and transversions are nonsynonymous at twofold degenerate sites in mammalian mitochondrial genes. We also found that nonsynonymous mutations at twofold degenerate sites are more neutral than nonsynonymous mutations at nondegenerate sites, and that the COI gene is subject to stronger purifying selection than is the cyt-b gene. A model is presented to integrate the effect of purifying selection, codon bias, DNA repair and GC content ons/v ratio of protein-coding genes. Correspondence to: X. Xia     

突脐孢属Brnl基因核苷酸序列比较及系统发育研究

对所有供试突脐孢菌株的Brnl基因(1,3,8-三羟基萘还原酶基因)扩增均获得PCR产物。序列比较表明：在种内各菌株间没有核苷酸序列长度变化,存在核苷酸序列简单代换;在种间核苷酸序列长度有变化,核苷酸的缺失或插入发生在内含子区;所有菌株编码区核苷酸序列长度相同;在种内水平氨基酸序列没有差别,显示出高度的保守性。利用最大简约法(Maximum Parsimony)和邻近结合法(Neighbor-joining)构建系统发育树,两个系统发育树的拓扑结构相似,不同种在不同的分支上。Brnl基因适合突脐孢属种级水平的分子系统学研究。