Determinants of rate variation in mammalian DNA sequence evolution

Attempts to analyze variation in the rates of molecular evolution among mammalian lineages have been hampered by paucity of data and by nonindependent comparisons. Using phylogenetically independent comparisons, we test three explanations for rate variation which predict correlations between rate variation and generation time, metabolic rate, and body size. Mitochondrial and nuclear genes, protein coding, rRNA, and nontranslated sequences from 61 mammal species representing 14 orders are used to compare the relative rates of sequence evolution. Correlation analyses performed on differences in genetic distance since common origin of each pair against differences in body mass, generation time, and metabolic rate reveal that substitution rate at fourfold degenerate sites in two out of three protein sequences is negatively correlated with generation time. In addition, there is a relationship between the rate of molecular evolution and body size for two nuclear-encoded sequences. No evidence is found for an effect of metabolic rate on rate of sequence evolution. Possible causes of variation in substitution rate between species are discussed.     

Mammalian DNA replication: mutation biases and the mutation rate

Experimental studies have shown that the fidelity of DNA replication can be affected by the concentrations of free deoxyribonucleotides present in the cell. Replication of mammalian chromosomes is achieved using pools of newly-synthesized deoxyribonucleotides which fluctuate during the cell cycle. Since regions of mammalian chromosomes are replicated sequentially, there is the potential for differences among mammalian loci in both the relative and absolute frequencies of the various transitional and transversional mutations which may occur. Where these mutations are effectively neutral, at silent sites in genes and in non-coding sequences, this may result in different rates of evolution and in different base compositions, as have been observed in data from mammalian genes. A simple model of the DNA replication process is developed to describe how the mutation rate could be affected by the G + C contents of the deoxyribonucleotide pools and of the replicating DNA. Mutation rates are predicted to vary from locus to locus; only in the particular case of identical G + C contents in the DNA locus and the deoxyribonucleotide pools, and no proofreading, will the mutation rate be uniform over all loci.     

Lineage-specific evolutionary rate in mammalian mtDNA

The existence of a lineage-specific nucleotide substitution rate in mammalian mtDNA has been investigated by analyzing the mtDNA of all available species, that is, 35 complete mitochondrial genomes from 14 mammalian orders. A detailed study of their evolutionary dynamics has been carried out on both ribosomal RNA and first and second codon positions (P12) of H-strand protein-coding genes by using two different types of relative-rate tests. Results are quite congruent between ribosomal and P12 sites. Significant rate variations have been observed among orders and among species of the same order. However, rate variation does not exceed 1.8-fold between the fastest (Proboscidea and Primates) and the slowest (Perissodactyla) evolving orders. Thus, the observed mitochondrial rate variations among taxa do not invalidate the suitability of mtDNA for drawing mammalian phylogeny. Dependence of evolutionary rate differences on variations in mutation and/or fixation rates was examined. Body size, generation time, and metabolic rate were tested, and no significant correlation was observed between them and the taxon-specific evolutionary rates, most likely because the latter might be influenced by multiple overlapping variable constraints.     

Rates of nucleotide substitution in Cornaceae (Cornales)-Pattern of variation and underlying causal factors

Identifying causes of genetic divergence is a central goal in evolutionary biology. Although rates of nucleotide substitution vary among taxa and among genes, the causes of this variation tend to be poorly understood. In the present study, we examined the rate and pattern of molecular evolution for five DNA regions over a phylogeny of Cornus, the single genus of Cornaceae. To identify evolutionary mechanisms underlying the molecular variation, we employed Bayesian methods to estimate divergence times and to infer how absolute rates of synonymous and nonsynonymous substitutions and their ratios change over time. We found that the rates vary among genes, lineages, and through time, and differences in mutation rates, selection type and intensity, and possibly genetic drift all contributed to the variation of substitution rates observed among the major lineages of Cornus. We applied independent contrast analysis to explore whether speciation rates are linked to rates of molecular evolution. The results showed no relationships for individual genes, but suggested a possible localized link between species richness and rate of nonsynonymous nucleotide substitution for the combined cpDNA regions. Furthermore, we detected a positive correlation between rates of molecular evolution and morphological change in Cornus. This was particularly pronounced in the dwarf dogwood lineage, in which genome-wide acceleration in both molecular and morphological evolution has likely occurred.     

Population size and molecular evolution on islands

The nearly neutral theory predicts that the rate and pattern of molecular evolution will be influenced by effective population size (Ne), because in small populations more slightly deleterious mutations are expected to drift to fixation. This important prediction has not been widely empirically tested, largely because of the difficulty of comparing rates of molecular evolution in sufficient numbers of independent lineages which differ only in Ne. Island endemic species provide an ideal test of the effect of Ne on molecular evolution because species restricted to islands frequently have smaller Ne than closely related mainland species, and island endemics have arisen from mainland lineages many times in a wide range of taxa. We collated a dataset of 70 phylogenetically independent comparisons between island and mainland taxa, including vertebrates, invertebrates and plants, from 19 different island groups. The rate of molecular evolution in these lineages was estimated by maximum likelihood using two measures: overall substitution rate and the ratio of non-synonymous to synonymous substitution rates. We show that island lineages have significantly higher ratios of non-synonymous to synonymous substitution rates than mainland lineages, as predicted by the nearly neutral theory, although overall substitution rates do not differ significantly.     

Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages

The nonsynonymous (amino acid-altering) to synonymous (silent) substitution rate ratio (omega = d(N)/d(S)) provides a measure of natural selection at the protein level, with omega = 1, >1, and <1, indicating neutral evolution, purifying selection, and positive selection, respectively. Previous studies that used this measure to detect positive selection have often taken an approach of pairwise comparison, estimating substitution rates by averaging over all sites in the protein. As most amino acids in a functional protein are under structural and functional constraints and adaptive evolution probably affects only a few sites at a few time points, this approach of averaging rates over sites and over time has little power. Previously, we developed codon-based substitution models that allow the omega ratio to vary either among lineages or among sites. In this paper we extend previous models to allow the omega ratio to vary both among sites and among lineages and implement the new models in the likelihood framework. These models may be useful for identifying positive selection along prespecified lineages that affects only a few sites in the protein. We apply those branch-site models as well as previous branch- and site-specific models to three data sets: the lysozyme genes from primates, the tumor suppressor BRCA1 genes from primates, and the phytochrome (PHY) gene family in angiosperms. Positive selection is detected in the lysozyme and BRCA genes by both the new and the old models. However, only the new models detected positive selection acting on lineages after gene duplication in the PHY gene family. Additional tests on several data sets suggest that the new models may be useful in detecting positive selection after gene duplication in gene family evolution.     

Nucleotide sequence analysis of the lemur beta-globin gene family: evidence for major rate fluctuations in globin polypeptide evolution

Lemur beta-related globin genes have been isolated and sequenced. Orthology of prosimian and human epsilon-, gamma-, and beta-related globin genes was established by dot-matrix analysis. All of these lemur globin genes potentially encode functional beta-related globin polypeptides, though precisely when the gamma-globin gene is expressed remains unknown. The organization of the 18-kb brown lemur beta-globin gene cluster (5' epsilon-gamma-[psi eta-delta]-beta 3') is consistent with its evolution by contraction via unequal crossing-over from the putative ancestral mammalian beta-globin gene cluster (5' epsilon-gamma- eta-delta-beta 3'). The dwarf lemur nonadult globin genes are arranged as in the brown lemur. Similar levels of synonymous (silent) nucleotide substitutions and noncoding DNA sequence differences have accumulated between species in all of these genes, suggesting a uniform rate of noncoding DNA divergence throughout primate beta-globin gene clusters. These differences are comparable with those observed in the nonfunctional psi eta pseudogene and have therefore accumulated at the presumably maximal neutral rate. In contrast, nonsynonymous (replacement) nucleotide substitutions show a significant heterogeneity in distribution for both the same gene in different lineages and different genes in the same lineage. These major fluctuations in replacement but not silent substitution rates cannot be attributed to changes in mutation rate, suggesting that changes in the rate of globin polypeptide evolution in primates is not governed solely by variable mutation rates.      

On the constancy of the evolutionary rate of cistrons

Summary The variations of evolutionary rates in hemoglobins and cytochrome c among various lines of vertebrates are analysed by estimating the variance. The observed variances appear to be larger than expected purely by chance.If the amino acid substitutions in evolution are the result of random fixation of selectively neutral or nearly neutral mutations, the evolutionary rate of cistrons can be represented by the integral of the product of mutation rate and fixation probability in terms of selective values around the neutral point. This integral is called the effective neutral mutation rate.The influence of effective population number and generation time on the effective neutral mutation rate is discussed. It is concluded that the uniformity of the rate of amino acid substitutions over diverse lines is compatible with random fixation of neutral or very slightly deleterious mutations which have some chance of being selected against during the course of substitution. On the other hand, definitely advantageous mutations will introduce significant variation in the substitution rate among lines. Approximately 10% of the amino acid substitutions of average cistrons might be adaptive and create slight but significant variations in evolutionary rate among vertebrate lines, although the uniformity of evolutionary rate is still valid as a first approximation.Contribution No. 813 from the National Institute of Genetics, Mishima, Shizuokaken 411 Japan. Aided in part by a grant-in-aid from the Ministry of Education, Japan.     

Evolutionary rate variation in Old World monkeys

We analysed over 8 million base pairs of bacterial artificial chromosome-based sequence alignments of four Old World monkeys and the human genome. Our findings are as follows. (i) Genomic divergences among several Old World monkeys mirror those between well-studied hominoids. (ii) The X-chromosome evolves slower than autosomes, in accord with ‘male-driven evolution’. However, the degree of male mutation bias is lower in Old World monkeys than in hominoids. (iii) Evolutionary rates vary significantly between lineages. The baboon branch shows a particularly slow molecular evolution. Thus, lineage-specific evolutionary rate variation is a common theme of primate genome evolution. (iv) In contrast to the overall pattern, mutations originating from DNA methylation exhibit little variation between lineages. Our study illustrates the potential of primates as a model system to investigate genome evolution, in particular to elucidate molecular mechanisms of substitution rate variation.     

Misconceptions and false expectations in neutral evolution

Neutral evolution results from random recurrent mutation and genetic drift. A small part of random evolution, that which is related to protein or DNA polymorphisms, is the subject of the Neutral Theory of Evolution. One of the foundations of this theory is the demonstration that the mutation rate (m) is equal to the substitution rate. Since both rates are independent of population size, they are independent of drift, which is dependent upon population size. Neutralists have erroneously equated the substitution rate with the fixation rate, despite the fact that they are antithetical conceptions. The neutralists then applied the random walk stochastic model to justify alleles or bases that were fixated or eliminated. In this model, once the allele or base frequencies reach the monomorphic states (values of 1.0 or 0.0), the absorbing barriers, they can no longer return to the polymorphic state. This operates in a pure mathematical model. If recurrent mutation occurs (as in biotic real systems) fixation and elimination are impossible. A population of bacteria in which m = 10(-8) base mutation (or substitution)/site/generation and the reproduction rate is 1000 cell cycle/year should replace all its genome bases in approximately 100,000 years. The expected situation for all sites is polymorphism for the four bases rather than monomorphism at 1.0 or 0.0 frequencies. If fixation and elimination of a base for more than 500,000 years are impossible, then most of the neutral theory is untenable. A new complete neutral model, which allows for recurrent substitutions, is proposed here based on recurrent mutation or substitution and drift alone. The model fits a binomial or Poisson distribution and not a geometric one, as does neutral theory.     

Synonymous and nonsynonymous rate variation in nuclear genes of mammals

A maximum likelihood approach was used to estimate the synonymous and nonsynonymous substitution rates in 48 nuclear genes from primates, artiodactyls, and rodents. A codon-substitution model was assumed, which accounts for the genetic code structure, transition/transversion bias, and base frequency biases at codon positions. Likelihood ratio tests were applied to test the constancy of nonsynonymous to synonymous rate ratios among branches (evolutionary lineages). It is found that at 22 of the 48 nuclear loci examined, the nonsynonymous/synonymous rate ratio varies significantly across branches of the tree. The result provides strong evidence against a strictly neutral model of molecular evolution. Our likelihood estimates of synonymous and nonsynonymous rates differ considerably from previous results obtained from approximate pairwise sequence comparisons. The differences between the methods are explored by detailed analyses of data from several genes. Transition/transversion rate bias and codon frequency biases are found to have significant effects on the estimation of synonymous and nonsynonymous rates, and approximate methods do not adequately account for those factors. The likelihood approach is preferable, even for pairwise sequence comparison, because more-realistic models about the mutation and substitution processes can be incorporated in the analysis. Received: 17 May 1997 / Accepted: 28 September 1997     

Substitution rates in hepatitis delta virus

Substitution rates were estimated for the coding and noncoding regions of the hepatitis delta virus (HDV). The estimated rates of synonymous substitution in HDV were lower than the rates of substitution at nonsynonymous sites and in the noncoding region. HDV has lower synonymous substitution rates than the hepatitis C virus, though both are RNA viruses. The relatively low rate of synonymous substitution in HDV may be due to a strong preference of G and C nucleotides at third codon positions. Variation in substitution rate among HDV lineages may be correlated with the clinical development of the HDV-induced hepatitis. The phylogenetic tree inferred for 24 HDV strains reveals similarities between lineages isolated from the same geographic region. Correspondence to: W.-H. Li     

Inverse relationship between longevity and evolutionary rate of mitochondrial proteins in mammals and birds

Recently, an unexpected, positive correlation between the rate of evolution of mitochondrial proteins and longevity was reported. Here we re-analyze this relationship in various mammalian lineages using a bayesian phylogenetic analysis of amino-acid sequences, allowing for variable evolutionary rates across sites and species. A negative relationship between protein evolutionary rate and species longevity is reported for all oxidative phosphorylation complexes. A detailed analysis of the cytochrome b in 528 mammals reinforced this result, which contradicts previous publications. Reconducting the analysis in birds yielded similar results. We explain the discrepancy between this and previous reports by our improved taxon sampling and more appropriate methodology: unlike distance-based methods, the tree-based bayesian approach can take into account the high variation of substitution rate across amino-acid sites, and the resulting multiple substitution events. We discuss how our analysis contradicts Rottenberg’s rationale, but does not dismiss his proposal of a longevity-dependent selective pressure on mitochondrial mutation rate in mammals and birds. This is because his interpretation invokes adaptation as the single evolutionary force at work, disregarding the effects of mutation, genetic drift, and purifying selection.     

Correlation between the substitution rate and rate variation among sites in protein evolution.

It is well known that the rate of amino acid substitution varies among different proteins and among different sites of a protein. It is, however, unclear whether the extent of rate variation among sites of a protein and the mean substitution rate of the protein are correlated. We used two approaches to analyze orthologous protein sequences of 51 nuclear genes of vertebrates and 13 mitochondrial genes of mammals. In the first approach, no assumptions of the distribution of the rate variation among sites were made, and in the second approach, the gamma distribution was assumed. Through both approaches, we found a negative correlation between the extent of among-site rate variation and the average substitution rate of a protein. That is, slowly evolving proteins tend to have a high level of rate variation among sites, and vice versa. We found this observation consistent with a simple model of the neutral theory where most sites are either invariable or neutral. We conclude that the correlation is a general feature of protein evolution and discuss its implications in statistical tests of positive Darwinian selection and molecular time estimation of deep divergences.     

Correlates of substitution rate variation in mammalian protein-coding sequences

Background  

Rates of molecular evolution in different lineages can vary widely, and some of this variation might be predictable from aspects of species' biology. Investigating such predictable rate variation can help us to understand the causes of molecular evolution, and could also help to improve molecular dating methods. Here we present a comprehensive study of the life history correlates of substitution rate variation across the mammals, comparing results for mitochondrial and nuclear loci, and for synonymous and non-synonymous sites. We use phylogenetic comparative methods, refined to take into account the special nature of substitution rate data. Particular attention is paid to the widespread correlations between the components of mammalian life history, which can complicate the interpretation of results.     

Variation in mutation rate and polymorphism among mitochondrial genes of Silene vulgaris

The prevailing wisdom of the plant mitochondrial genome is that it has very low substitution rates, thus it is generally assumed that nucleotide diversity within species will also be low. However, recent evidence suggests plant mitochondrial genes may harbor variable and sometimes high levels of within-species polymorphism, a result attributed to variance in the influence of selection. However, insufficient attention has been paid to the effect of among-gene variation in mutation rate on varying levels of polymorphism across loci. We measured levels of polymorphism in seven mitochondrial gene regions across a geographically wide sample of the plant Silene vulgaris to investigate whether individual mitochondrial genes accumulate polymorphisms equally. We found that genes vary significantly in polymorphism. Tests based on coalescence theory show that the genes vary significantly in their scaled mutation rate, which, in the absence of differences among genes in effective population size, suggests these genes vary in their underlying mutation rate. Further evidence that among-gene variance in polymorphism is due to variation in the underlying mutation rate comes from a significant positive relationship between the number of segregating sites and silent site divergence from an outgroup. Contrary to recent studies, we found unconvincing evidence of recombination in the mitochondrial genome, and generally confirm the standard model of plant mitochondria characterized by low substitution rates and no recombination. We also show no evidence of significant variation in the strength or direction of selection among genes; this result may be expected if there is no recombination. The present study provides some of the most thorough data on plant mitochondrial polymorphism, and provides compelling evidence for mutation rate variation among genes. The study also demonstrates the difficulty in establishing a null model of mitochondrial genome polymorphism, and thus the difficulty, in the absence of a comparative approach, in testing the assumption that low substitution rates in plant mitochondria lead to low polymorphism.     

The causes of synonymous rate variation in the rodent genome. Can substitution rates be used to estimate the sex bias in mutation rate?

Miyata et al. have suggested that the male-to-female mutation rate ratio (alpha) can be estimated by comparing the neutral substitution rates of X-linked (X), Y-linked (Y), and autosomal (A) genes. Rodent silent site X/A comparisons provide very different estimates from X/Y comparisons. We examine three explanations for this discrepancy: (1) statistical biases and artifacts, (2) nonneutral evolution, and (3) differences in mutation rate per germline replication. By estimating errors and using a variety of methodologies, we tentatively reject explanation 1. Our analyses of patterns of codon usage, synonymous rates, and nonsynonymous rates suggest that silent sites in rodents are evolving neutrally, and we can therefore reject explanation 2. We find both base composition and methylation differences between the different sets of chromosomes, a result consistent with explanation 3, but these differences do not appear to explain the observed discrepancies in estimates of alpha. Our finding of significantly low synonymous substitution rates in genomically imprinted genes suggests a link between hemizygous expression and an adaptive reduction in the mutation rate, which is consistent with explanation 3. Therefore our results provide circumstantial evidence in favor of the hypothesis that the discrepancies in estimates of alpha are due to differences in the mutation rate per germline replication between different parts of the genome. This explanation violates a critical assumption of the method of Miyata et al., and hence we suggest that estimates of alpha, obtained using this method, need to be treated with caution.     

Lineage-specific variation in slow- and fast-X evolution in primates

Theories predict that the evolutionary rates of X-linked regions can differ from those of autosomal regions. The male-biased mutation theory predicts a slower rate of neutral substitution on the X chromosome (slow-X evolution), as the X spends less time in male germlines, where more mutations originate per generation than in female germlines. The fast-X theory, however, predicts a faster rate of adaptive substitution on the X chromosome when newly arising beneficial mutations are, on average, partially recessive (fast-X evolution), as the X enjoys a greater efficacy of positive selection. The slow- and fast-X processes are expected to interact as the degree of male-biased mutation can in turn influence the relative rate of adaptive evolution on the X. Here, we investigate lineage-specific variation in, and the interaction of, slow- and fast-X processes using genomic data from four primates. We find consistent evidence for slow-X evolution in all lineages. In contrast, evidence for fast-X evolution exists in only a subset of lineages. In particular, the marmoset lineage, which shows the strongest evidence of fast-X, exhibits the lowest male mutation bias. We discuss the possible interaction between slow- and fast-X evolution and other factors that influence the degrees of slow- and fast-X evolution.     

Rates of molecular evolution in bacteria are relatively constant despite spore dormancy

Rates of molecular evolution are known to vary considerably among lineages, partially due to differences in life-history traits such as generation time. The generation-time effect has been well documented in some eukaryotes, but its prevalence in prokaryotes is unknown. "Because many species of Firmicute bacteria spend long periods of time as metabolically dormant spores, which could result in fewer DNA substitutions per unit time, they present an excellent system for testing predictions of the molecular clock hypothesis." To test whether spore-forming bacteria evolve more slowly than their non-spore-forming relatives, I used phylogenetic methods to determine if there were differences in rates of amino acid substitution between spore-forming and non-spore-forming lineages of Firmicute bacteria. Although rates of evolution do vary among lineages, I find no evidence for an effect of spore-formation on evolutionary rate and, furthermore, evolutionary rates are similar to those calculated for enteric bacteria. These results support the notion that variation in generation time does not affect evolutionary rates in bacterial lineages.