- Research article
- Open Access
Correlates of substitution rate variation in mammalian protein-coding sequences
BMC Evolutionary Biology volume 8, Article number: 53 (2008)
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.
We find that mitochondrial synonymous substitution rates, estimated from the 9 longest mitochondrial genes, show strong negative correlations with body mass and with maximum recorded lifespan. But lifespan is the sole variable to remain after multiple regression and model simplification. Nuclear synonymous substitution rates, estimated from 6 genes, show strong negative correlations with body mass and generation time, and a strong positive correlation with fecundity. In contrast to the mitochondrial results, the same trends are evident in rates of nonsynonymous substitution.
A substantial proportion of variation in mammalian substitution rates can be explained by aspects of their life history, implying that molecular and life history evolution are closely interlinked in this group. The strength and consistency of the nuclear body mass effect suggests that molecular dating studies may have been systematically misled, but also that methods could be improved by incorporating the finding as a priori information. Mitochondrial synonymous rates also show the body mass effect, but for apparently quite different reasons, and the strength of the relationship with maximum lifespan provides support for the hypothesis that mtDNA damage is causally linked to aging.
There is now a great deal of evidence that rates of DNA substitution can vary widely between closely related lineages [1–3], and while some of this variation is erratic and locus-specific, trends also apply consistently across many loci [2, 4–7]. There is also increasing evidence, particularly in vertebrates, that some of this lineage-specific rate variation may be predictable from aspects of a species' biology [8–17]. Uncovering this predictable rate variation is an important part of understanding the causes of molecular evolution . It may also bring practical benefits, because variation in the rate of substitution complicates the production of dated molecular phylogenies, which are increasingly relied upon in diverse areas of biology. If reliable correlates of rate variation can be identified, this information could be exploited to improve molecular dating methods [15, 16, 19, 20].
The present study investigates the causes and correlates of lineage-specific rate variation in mammalian protein coding sequences. Mammals were chosen because of the unrivalled availability of relevant data: mitochondrial and nuclear DNA sequences, life-history records and comprehensive phylogenetic information [21, 22]. In addition, the timings of ordinal and higher-level radiations of the mammals have been controversial, and systematic change in substitution rates has been proposed as an explanation of disagreements between molecular and palaeontological dating approaches [23–27].
Many previous studies have tested for correlates of mammalian rate variation [5, 6, 14, 15, 17, 28–32], but no clear consensus has emerged. One problem is that negative findings have been difficult to interpret because of issues surrounding statistical power in the comparative study of substitution rates. In some cases, sample size was inflated artificially by a failure to correct for phylogenetic non-independence. This means that single instances of evolutionary change can be counted multiple times, and this can increase both false positive and false negative error rates [33, 34]. Even studies that did take shared ancestry into account, e.g. , did not always use appropriate statistical methods or diagnostic tests of the parametric assumptions, and this too can yield misleading results [35, 36]. The question of statistical power is particularly vexed in the study of substitution rates, because rates cannot be measured directly, but must be inferred from DNA sequence data. For short sequences, or short periods of divergence, a small number of substitutions can make a large difference to the rate inferred, making the estimates very noisy [16, 37]. Accordingly, there is often a trade-off between the number of data points in an analysis, and the accuracy with which each data point estimates a change in rate. In the present study, we address these problems by using large multi-gene alignments, and the method of phylogenetically independent contrasts [33, 34], combined with new procedures for establishing a minimum comparison depth .
Another difficulty, particularly acute for mammals, stems from the widespread correlations between all aspects of their life histories and wider biology; these include not just strong and ubiquitous allometries, but also strong correlations between traits after correction for body size [22, 38–40]. These colinearities make regression analyses particularly difficult to interpret, because a significant result can always be plausibly attributed to some absent variable covarying with the predictor. For example, many comparative studies of mammalian aging have been criticised for failing to take into account the covariation of lifespan with body mass [13, 41, 42]. Here, the problem of covarying life history traits is mitigated by including, in a single multiple regression analysis, most of the correlates of vertebrate rate variation that have been identified or hypothesised in the literature. Specifically, the variables included are body mass [13, 14, 20, 28], which has often been used as a proxy for basal metabolic rate [11, 15, 17, 43–45], organismal generation time [8, 9, 14, 28], fecundity [10, 46], and maximum recorded lifespan [13, 32, 41]. These variables are defined in detail below, and an assessment of the various causal hypotheses with which they have been associated follows in the discussion.
The single variable regressions, summarised in Table 1, show that of the four variables tested, body mass and maximum lifespan are individually significant predictors of mitochondrial synonymous substitution rates: Species with greater mass, or longer lifespan tend to have slower rates of synonymous substitution. The significant regressions are plotted in Figure 1a–d, together with their associated raw cross-species plots, which show closely similar trends. A multiple regression analysis, including pairs with measurements for all four predictor traits, shows that only maximum lifespan remains significant (Table 2). Furthermore, model simplification shows that lifespan alone explains almost as much variation as does the four-predictor model (Table 2). Table 3 contains results for the larger subset of points for which lifespan and body mass measurements were available (i.e., including those pairs lacking fecundity or generation time data). Again, maximum lifespan was found to be the sole significant predictor. Diagnostic tests identified the sperm whale pair, Kogia breviceps-Physeter catodon, as a weak outlier, and excluding this pair increases the r2to levels matching those in Table 2. (It is possible that the outlying nature of this point reflects the recognised difficulties obtaining accurate measurements of maximum lifespan for cetaceans: .)
To check for consistency across the major superordinal groups, we carried out separate analyses for the three clades with sufficient comparisons, namely Metatheria (marsupials), Laurasiatheria (Cetartiodactyla, Carnivora, Perissodactyla, Chiroptera, and some former Insectivora) and Euarchontoglires (Rodentia, Lagomorpha, Scandentia, Dermoptera and Primates) [see Additional file 1]. For both Metatheria and Laurasiatheria, results for all traits are closely consistent with each other, and with the complete data set (for example, for body mass we have Metatheria: n = 11, slope = -0.06, r2= 0.22; Laurasiatheria: n = 21, slope = -0.11, r2= 0.24). For the Euarchontoglires, by contrast, lifespan, fecundity and generation time had slopes close to zero, and no visible trend in the raw cross-species plots. Furthermore, in this group synonymous rates actually appear to increase with body mass (Figure 1e–f), albeit with a non-significant regression (Euarchontoglires: n = 10, slope = 0.10, r2= 0.25, p = 0.12). Accordingly, we repeated the full analyses with the Euarchontoglires excluded. For the two-variable regression (Table 3), results were not robust, with neither predictor reaching significance for the 27 remaining pairs, but both reaching significance when the outlying sperm whale pair was excluded (Table 3). However, for the subset of pairs with all variables measured, maximum lifespan was again identified as the sole significant predictor (Table 2). Excluding Euarchontoglires also led to r2 values that were much higher in both cases (Tables 2, 3).
In stark contrast to results for synonymous substitutions, mitochondrial nonsynonymous rates show no trend with any of the predictor variables (Tables 1 and 2). This applied equally to the independent comparisons, and to the raw cross-species plots (not shown).
Plots for the nuclear synonymous data set are shown in Figure 2. Table 1 shows that body mass, generation time and fecundity are all significant predictors of nuclear synonymous rates. Lifespan, by contrast, shows no trend, but this is attributable to two clear outliers, appearing as such in both the diagnostic tests of the regressions (Fig. 2c) and in the raw cross-species scatter-plots (Fig. 2d). (The two small values are for Cynocephalus variegatus, the malayan flying lemur, and Amblysomus hottentotus, the hottentot golden mole). Removing these two outlying points makes maximum lifespan, too, a significant predictor of synonymous rates (Table 1). Of course, these outliers could represent true biological variation, but problems with maximum recorded lifespan as a statistic make measurement error particularly likely in this case [13, 22, 47, 48]. (For example, recorded lifespan can only increase with new observations, and so is strongly dependent on the number of animals sampled; a species' suitability for captivity will also have an important influence.) If outliers are excluded, then for each of the predictors, r2 values are much higher than for the mitochondrial data set, even when corrections for reduced sample size are made. Furthermore, and again in contrast to the mitochondrial results, all effects apply consistently to Laurasiatheria and Euarchontoglires (the nuclear data set contains no Metatheria).
The multiple regression results for the nuclear data (Table 2) are unfortunately difficult to interpret. Model simplification leads to a model with generation time and fecundity remaining. But neither coefficient is individually significant, and model fit is poor in several respects (for example, the Shapiro-Wilks test shows significant departures from normality). More telling are the two-variable regressions with body mass and lifespan (Table 3). For nuclear rates, body mass is found to be the superior predictor, whether or not the outlying lifespan values are excluded, and, in either case, lifespan drops out in model simplification. This is the opposite result to that obtained for the mitochondrial data set (Table 3).
Finally, and again in direct contrast to the mitochondrial results, for the nuclear data set the same trends are observed for both synonymous and non-synonymous rates (Tables 1, 2, 3). While the non-synonymous data are noisier, with generation time and longevity effects not reaching individual significance, the slopes of the regressions are very similar (Table 1), and body mass is again favoured over maximum lifespan in multiple regression (Table 3).
It is often possible to attribute an apparent change in the rate of substitution to some other factor, such as misspecification of the substitution model, inaccurate divergence dates, changes in base composition, or a transient and locus-specific burst of adaptive changes. But the identification of significant predictors of rates, within the framework of phylogenetically independent comparisons (Table 1; Figures 1, 2), is difficult to reconcile with any of these explanations, and argues strongly that changes of substitution rate are frequent and substantial within the mammals. Conversely, the various sources of error inherent in the estimation of rates makes it probable that the r2 values estimated here – high as they are – are systematic underestimates of the true relationship between rate and trait. (Specifically, rate differs from traits such as body mass, in that it must be estimated from the number of substitutions accrued stochastically over a period of time, which must in turn be estimated from comparison of sequence data, and all of these factors can obscure the true relationship between trait and rate .)
Implications for molecular dating
From the perspective of molecular phylogenetics, the strength and consistency of the body size effect across our nuclear data set is of particular interest: If such an effect applies, concerted changes in the average body size of a clade over time can create systematic biases in the molecular date estimates for that clade, even if variable-rate dating methods are used [16, 19, 49]. The same biases can apply in the absence of concerted body size change, if larger or smaller taxa are over-represented in the data sampled. For these reasons, if there was a concerted increase in mammalian size between the late Cretaceous and early Tertiary, or if larger mammals have been disproportionately sampled, then molecular date estimates of the ordinal level radiation of mammals could have been systematically misled [16, 23–26]. While the widespread rate variation presents a problem for molecular dating, the identification of robust predictors of rate gives grounds for optimism. In particular, body size is both widely sampled for extant species, and can be estimated for extinct or ancestral taxa from their fossil record. Consequently, this information could be used to improve future molecular dating methods, allowing the development of "corrected molecular clocks", or empirically informed priors for Bayesian approaches [15, 19, 20].
Mitochondrial rate variation, and mitochondrial theories of aging
For the mitochondrial data set, maximum lifespan was found to be the strongest predictor of synonymous rate (Table 1), and was the sole significant predictor in multiple regressions (Tables 2, 3). This implies that the mitochondrial body size effect, e.g.,  can be attributed solely to the correlation of mammalian body mass and lifespan [13, 39–41]. In contrast, for the nuclear data set, maximum lifespan dropped out of the model in favour of body mass – even when outlying lifespan values were removed (Tables 1, 2, 3). Together, these findings support the longstanding, though still controversial, theory that mitochondrial DNA damage is causally linked to aging [17, 50–54]. The possibility that this link is causal is strengthened, ironically, by the many problems with maximum recorded lifespan as a statistic [13, 47, 48]. It has been shown many times that maximum lifespan is a poor proxy for typical longevity in the wild, with the two quantities differing substantially and unpredictably in many cases [55, 56]. As such, maximum recorded lifespan relates only weakly to the realised life histories on which selection can act [39, 40, 57].
A link between mtDNA damage and mammalian aging also has some direct experimental support. Premature aging has been reported in mice expressing defective mitochondrial DNA polymerase [58–60], and an extension of youthful lifespan shown to result from the over-expression of mitochondrially-targeted catalase, but not of nuclear-targeted catalase . However, there are also reports that mitochondrial point mutations have no effect on mouse lifespan [62, 63]. These latter results might be consistent with recent modifications of the mitochondrial theory of aging, implicating mtDNA deletions rather than point mutations [51, 52, 54], and consistent with our results if the two types of DNA damage are highly correlated in nature.
Although the precise mechanisms linking mtDNA damage to senescence are much disputed, most theories implicate Reactive Oxygen Species (ROS): Mitochondria are the major site of ROS production, and mtDNA a major site of oxidative damage [50–52]. ROS production is an inevitable by-product of aerobic metabolism, and this raises the possibility that the true determinant of substitution rates is basal metabolic rate (BMR) [11, 15, 43, 64]. BMR was not included in the present study, but there are several reasons to believe that its inclusion would not have altered our conclusions. First, the correlation of BMR with body mass is particularly strong in mammals [44, 65], making it very unlikely that lifespan would be identified as a better predictor than body size if BMR were the true causal factor. Second, the mechanistic basis of the metabolic rate hypothesis is doubtful, both because the basal rate is a poor measure of total energy metabolism , and because mammalian metabolism can be decoupled from ROS production by various means [45, 65, 66]. Nuclear DNA, where the body mass effect is strongest (Tables 1 and 3), also appears to be protected from mitochondrially-generated ROS . Third, there is some evidence that mtDNA damage can accelerate aging without increasing oxidative stress [59, 60]. Finally, previous analyses of substitution rates have not identified BMR as a significant predictor when body size was also included in the model [12, 14, 45].
Particularly revealing case studies come from species whose ecologies have led to departures from the mammalian norm. For example, the naked mole rat, particularly long-lived and metabolically slow for its size, shows peculiarly low levels of anti-oxidative defences, and correspondingly high levels of oxidative damage in its nuclear DNA , but this might be attributable, paradoxically, to reduced levels of ROS production in the mitochondria, implying reduced levels mtDNA damage [68, 69]. Just such a situation appears to explain the greater longevity of birds compared to mammals . It would also be of great interest to compare patterns of rate variation in Chiroptera, a group not well represented in the present study [17, 47, 48, 57, 70].
A final puzzling aspect of the mitochondrial results is the anomalous patterns shown by the Euarchontoglires – the group that contains both rodents and primates, and so all of the most highly studied mammalian models. For this superordinal group, synonymous rate appears to increase with body mass (Fig. 1e–f), and shows no trend with lifespan. One possible explanation is that these results reflect the confounding influence of weak purifying selection. Theory has shown that, if certain assumptions hold, we can expect a negative correlation between the mutation rate per year, and the fixation probability of weakly deleterious mutants [18, 28, 71]. For example, if species with higher mutation rates also have larger populations, then the increased efficacy of purifying selection could lead to a reduction in the rate of substitution due to genetic drift. It is possible that weak purifying selection acts on mitochondrial synonymous sites in mammals , and there is also some evidence that the positive relationship between mutation rate and population size does hold . This effect is likely to be dampened by the complete linkage of mitochondrial DNA, which weakens the dependence of substitution rate on population size [71, 74], but typical mammalian populations may be small enough for drift-based effects to remain important . Furthermore, when linkage is tight, a negative correlation between mutation rate and the fixation probability could still hold, if populations with higher mutation rates per year, also have higher rates of adaptive substitution per generation .
The hypothetical explanation given above remains incomplete, however, because it remains unclear why effects counteracting the variation in mutation rates should be so much stronger in the Euarchontoglires than in the other major mammalian groups.
Nuclear rate variation and mammalian life history
Compared to the mitochondrial data, our nuclear data set is smaller in every respect, and particularly in terms of the fraction of the genome sampled. However, for the loci studied, the influence of life history on substitution rate appears to be even more pervasive.
For example, in common with previous studies, no successful predictor was found for mitochondrial non-synonymous changes [5, 14, 30], but our nuclear results were similar for both classes of site (Tables 1, 2, 3). The nuclear results contrast with some previous empirical studies [14, 29], and also contradict theoretical predictions, mentioned above, that the evolution of selected sites should be relatively immune to variation in life history . We cannot know whether our nuclear data set is representative, in this respect, of nuclear loci in general, but examination of the dN/dS ratios yields no evidence that levels of selective constraint are anomalous (values for the 32 lineages in the main data set range from 0.046 to 0.200 [see Additional file 1], which can be compared to results from complete genomes of primates and rodents ). Furthermore, the theoretical predictions rely on assumptions that are unlikely to hold at all loci (e.g., the influence of population size on levels of selective constraint may be slight for highly leptokurtic distributions of selective effects ). Our results are also consistent with observations that a substantial component of nuclear nonsynonymous rate variation in mammals is lineage-specific , and with the strong correlation between nuclear dS and dN  (for contrasting patterns in mitochondrial sequences, see [2, 5, 72]).
Results for our nuclear loci also contrast with the mitochondrial patterns in terms of the relative success the individual predictors. A strong effect of lifespan on nuclear rates can be discounted, but no single other factor was unambiguously favoured (Table 2). Nevertheless, the contrast between the two sets of results allows us to draw some tentative conclusions about various causal hypotheses.
For example, the argument for body mass as a true causal factor stems from evidence that larger bodied mammals suffer increased risk of cancer [13, 42], and that this might select for increased DNA repair [13, 41, 77]. But this hypothesis is difficult to reconcile with the comparative weakness of the body size effect in mtDNA, where links between mutation and carcinogenesis are well established [62, 77]. In contrast, the success of generation time as a predictor of nuclear, but not mitochondrial rates, might reflect differences in the biology of the two genomes. In particular, germline mosaicism , where a large fraction of nuclear mutations appear in just one or two meiotic divisions, is certain to strengthen the correlation between generation time and mutation rate per year; while the replication of mtDNA is potentially decoupled from cell division , weakening such a dependency.
A more perplexing result is the apparent increase in nuclear rates with fecundity. Non-significant for mitochondrial rates, this predictor has an r2 of around 50% for both classes of nuclear site (Table 1). It has been hypothesised that increases in fecundity might entail an increased mutation rate , or reflect increased variance in offspring viability, and so lower effective population size . But given the limited variation in mammalian fecundity, it is unlikely that either of these explanations could explain our results [22, 80]. A third intriguing possibility views germline mutation rate as an integral part of life history strategy [22, 39, 40], suggesting that species with fewer offspring should invest more heavily in each, implying a selective pressure to lower germline mutation rates .
We have presented a comprehensive study of substitution rate variation in protein-coding sequences across the mammals, and shown that for mitochondrial synonymous sites, and for both synonymous and nonsynonymous sites in 6 nuclear loci, a substantial fraction of the between-lineage rate variation can be explained by aspects of life history.
The results imply that molecular dating studies of mammalian evolution might have been misled, particularly if there was a systematic change in the life history of the clade. However, future methods could exploit results presented here, incorporating measurements of body size as a priori information about the substitution rates to be inferred.
While both mitochondrial and nuclear rates show an inverse correlation with body mass, the results differ in important respects, implying that the causal mechanisms are quite different in the two sets of loci. The success of maximum lifespan as a predictor of mitochondrial rates, and its comparative failure to predict nuclear rates, provides support for theories linking mtDNA damage to aging. Causal interpretation of the nuclear results is more difficult, but one conclusion is clear. Molecular change may be decoupled from phenotypic change in the sense that many changes in germline DNA may lead to vanishingly small changes in phenotype [18, 28], but mammalian molecular evolution is nevertheless a part of life history evolution, and the treatment of germline mutation rates as a component of life history deserves further attention.
All genetic sequence data were obtained from GenBank , and aligned by eye using Se-Al . Accession Numbers are listed in the Data Supplement [Additional file 1], and alignments are available on request.
For the mitochondrial dataset, the 9 longest protein-coding genes (ATP6, COI, COII, COIII, CYTB, ND1, ND2, ND4, ND5) were aligned for 160 mammalian species. (The four short genes, less comprehensively sampled, were excluded to increase taxonomic coverage.) The resulting data set allowed us to choose at least one comparison pair from the monotremes, from 5 of the 7 marsupial orders, and from 13 of the 19 'molecular consensus' placental orders . The absent orders were small, containing just 100 species in total, of which half are Afrosoricida.
For the nuclear data set, there was a clear tradeoff between alignment length and taxonomic coverage. We chose a data set consisting of partial coding sequences from 6 nuclear loci, which were available for 58 species. This data set allowed us to choose pairs only for the Eutheria, but 16 of their 19 orders were represented. Genes and approximate sequence lengths are ADRB2 ~830 bp; ATP7A ~680 bp; BDNF ~590 bp; CNR1 ~1000 bp; EDG1 ~980 bp; RAG2 ~740 bp. In addition to annotated sequences, BLAST searches were carried out on genomic contigs from Echinops telfairi (lesser hedgehog tenrec) and Myotis lucifugus (little brown bat; Lindblad-Toh, K., J. L. Chang, S. Gnerre, M. Clamp, and E. S. Lander, unpublished), confirming orthology by reciprocal BLAST.
Life history data
In most cases, the life history data were taken from the PanTHERIA database of mammalian life history and ecology (K. E. Jones, J. Bielby, A. Purvis, D. Orme, A. Teacher, J. L. Gittleman, R. Grenyer, et al. unpublished; ). Body mass measurements were the unique median of adult (or age unspecified) mass, of males and females, based on the GLM equation ln(body mass) = species + sex. For two Chiropterans in the mitochondrial data set (Rhinolophus monoceros and Chalinolobus tuberculatus) this value was extrapolated from head-body length. Our measure of fecundity was the product of litter size, and litters per year. Litter size is the full median of offspring number born per litter per female, counted before birth, at birth, or after birth, based on the GLM equation ln(litter size) = species + litter size definition. Litters per year is the full median of the number of litters per female per year for non-captive individuals. Maximum lifespan was simply the maximum recorded adult age from a wild or captive individual. To supplement the PanTHERIA data, we obtained additional measurements from the database AnAge (Build 9, Feb. 2006; ), which contains values from multiple recent compilations from the literature. Our preferred proxy of generation time was the unique median of age at first birth. Values from the PanTHERIA database were supplemented by equivalent data from the compilation of Wooton . Because coverage remained insufficient for a robust analysis, we also included some measurements where generation time was defined as age at sexual maturity plus gestation period , with values obtained from the AnAge database (Build 9, Feb. 2006; ).
Choice of phylogenetically independent comparison pairs
For the comparative analysis, phylogenetically-independent sister pairs were chosen from the recent species-level supertree of all mammals [21, 26]. Comparisons between reconstructed states at internodes  were not included, as these are problematic in the comparative study of substitution rates [see Additional file 2]. For three comparisons, no suitable outgroup was available, and so molecular branch lengths were estimated from a split along one of the lineages, rather than from the pair's common ancestor, making the pair paraphyletic with respect to another pair. In each case the distance of the chosen node from the common ancestor of the pair was small compared to the total divergence between the pair. Nevertheless, in all three cases, rate estimates were corrected to take into account the different time periods represented by the two molecular branch lengths. Repeating analyses with these pairs excluded was found to make little difference to the results (not shown).
Branch length estimation
Nonsynonymous and synonymous molecular branch lengths for each pair were estimated via maximum likelihood, using the codon-based substitution model of Goldman and Yang [85, 86]. Both the overall substitution rate and the ratio of nonsynonymous to synonymous changes were allowed to take branch-specific values, and so estimating branch lengths for the complete tree risks overparameterisation. Furthermore, whole-tree estimation can be unreliable if nuisance parameters, such as base composition, vary across groups. For these reasons, the data were split into small groups of 3–8 related species, and branch lengths estimated separately for these small subtrees. For both nuclear and mitochondrial data sets, results reported here are for concatenated multi-locus alignments. Analyses were also carried out for individual loci, and for a lineage-effect component of rate variation, estimated via the method of Smith and Eyre-Walker [2, 4], but as results did not differ qualitatively from those obtained by the simpler method of concatenating the sequences, they are not reported further. In a few cases, one or more branches of a subtree showed signs of saturation, but the relevant pairs fit the same trends as the remainder of the data set, implying that signal was present in the sequence data, and so they were not excluded from the analyses.
Full details of the comparison pairs, the molecular branch lengths, and life history measurements are included in the Data Supplement [Additional file 1].
The independent comparisons were analysed with multiple regression, forcing all regression lines through the origin [34, 35]. For results to be valid, it is important to meet the assumptions of the parametric test, and this typically involves transforming the trait measurements, and standardising comparisons to account for their different periods of divergence . For this purpose, estimated divergence dates were taken from . (An advantage of the use of sister pairs is that the dates were required solely for this purpose, and have no effect on the estimated differences in log rates.) To assess the standardisations and transformations used, we employed standard regression diagnostics, and customised tests [35, 36]. These methods were extended to solve the special problems posed by substitution rate data . The extended methods use patterns in the contrast variance to define a minimum depth below which comparison pairs were excluded [see Additional file 2]. In this way, we hoped to prevent the strength of any effect being masked by stochastic noise in the substitution process [1, 4, 18]. All statistical tests were implemented in R .
We also produced raw cross-species plots, for which the absolute rates along each lineage were calculated using the estimated divergence dates. These plots were solely illustrative, but they allowed us to include taxa excluded from the independent contrasts analysis, either because of missing life-history data for one member of the pair, or because of the diagnostic tests.
Bromham L, Penny D: The modern molecular clock. Nat Rev Genet. 2003, 4: 216-224. 10.1038/nrg1020.
Smith NGC, Eyre-Walker A: Partitioning the variation in mammalian substitution rates. Mol Biol Evol. 2003, 20: 10-17. 10.1093/molbev/msg003.
Bininda-Emonds ORP: Fast genes and slow clades: comparative rates of molecular evolution in mammals. Evolutionary Bioinformatics. 2007, 3: 59-85.
Gillespie JH: Lineage effects and the index of dispersion of molecular evolution. Mol Biol Evol. 1989, 6: 636-647.
Gissi C, Reyes A, Pesole G, Saccone C: Lineage-specific evolutionary rate in mammalian mtDNA. Mol Biol Evol. 2000, 17: 1022-1031.
Li W-H, Gouy M, Sharp P, O'Huigin C, Yang Y-W: Molecular phylogeny of Rodentia, Lagomorpha, Artiodactyla and Carnivora and molecular clocks. Proc Natl Acad Sci USA. 1990, 87: 6703-6707. 10.1073/pnas.87.17.6703.
Williams EJB, Hurst LD: Is the synonymous substitution rate in mammals gene-specific?. Mol Biol Evol. 2002, 19: 1395-1398.
Laird CD, McConaughy BL, McCarthy BJ: Rate of fixation of nucleotide substitutions in evolution. Nature. 1969, 224: 149-154. 10.1038/224149a0.
Sarich VM, Wilson AC: Generation time and genomic evolution in primates. Science. 1973, 179: 1144-1147. 10.1126/science.179.4078.1144.
Britten RJ: Rates of DNA sequence evolution differ between taxonomic groups. Science. 1986, 231: 1393-1398. 10.1126/science.3082006.
Martin AP, Palumbi SR: Body size, metabolic rate, generation time and the molecular clock. Proc Natl Acad Sci USA. 1993, 90: 4087-4091. 10.1073/pnas.90.9.4087.
Mooers AØ, Harvey PH: Metabolic rate, generation time and the rate of molecular evolution in birds. Mol Phylogenet Evol. 1994, 3: 344-350. 10.1006/mpev.1994.1040.
Promislow DE: DNA repair and the evolution of longevity: a critical analysis. J Theor Biol. 1994, 170: 291-300. 10.1006/jtbi.1994.1190.
Bromham L, Rambaut A, Harvey PH: Determinants of rate variation in mammalian DNA sequence evolution. J Mol Evol. 1996, 43: 610-621. 10.1007/BF02202109.
Gillooly JF, Allen AP, West GB, Brown JH: The rate of DNA evolution: Effects of body size and temperature on the molecular clock. Proc Natl Acad Sci USA. 2005, 102: 140-145. 10.1073/pnas.0407735101.
Fontanillas E, Welch JJ, Thomas JA, Bromham L: The influence of body size and diversification rate on molecular evolution during the Cambrian Explosion of animal phyla. BMC Evol Biol. 2007, 7: 95-10.1186/1471-2148-7-95.
Nabholz B, Glémin S, Galtier N: Strong variations of mitochondrial mutation rate across mammals – the longevity hypothesis. Mol Biol Evol. 2008, 25: 120-130. 10.1093/molbev/msm248.
Gillespie JH: The causes of molecular evolution. 1991, Oxford, UK: Oxford University Press
Welch JJ, Bromham L: Molecular dating when rates vary. Trends Ecol Evol. 2005, 20: 320-327. 10.1016/j.tree.2005.02.007.
Thomas JA, Welch JJ, Woolfit M, Bromham L: There is no universal molecular clock for invertebrates, but rate variation does not scale with body size. Proc Natl Acad Sci USA. 2006, 103: 7366-7371. 10.1073/pnas.0510251103.
Beck RMD, Bininda-Emonds ORP, Cardillo M, Liu FR, Purvis A: A higher-level MRP supertree of placental mammals. BMC Evol Biol. 2006, 6: 93-10.1186/1471-2148-6-93.
Bielby J, Mace GM, Bininda-Emonds ORP, Cardillo M, Gittleman JL, Jones KE, Orme CDL, Purvis A: The fast-slow continuum in mammalian life history: an empirical reevaluation. Am Nat. 2007, 169: 748-757. 10.1086/516847.
Bromham LD, Phillips MJ, Penny D: Growing up with dinosaurs: molecular dates and the mammalian radiation. Trends Ecol Evol. 1999, 14: 113-118. 10.1016/S0169-5347(98)01507-9.
Douzery EJ, Delsuc F, Stanhope MJ, Huchon D: Local molecular clocks in three nuclear genes: divergence times for rodents and other mammals and incompatibility among fossil calibrations. J Mol Evol. 2003, 57: S201-S213. 10.1007/s00239-003-0028-x.
Brochu CA, Sumrall CD, Theodor JM: When clocks (and communities) collide: Estimating divergence time from molecules and the fossil record. Journal of Palaeontology. 2004, 78: 1-6. 10.1666/0022-3360(2004)078<0001:WCACCE>2.0.CO;2.
Bininda-Emonds ORP, Cardillo M, Jones KE, MacPhee RDE, Beck RMD, Greyner R, Price SA, Vos RA, Gittleman JA, Purvis A: The delayed rise of the present-day mammals. Nature. 2007, 446: 507-512. 10.1038/nature05634.
Wible JR, Rougier GW, Novacek MJ, Asher RJ: Cretaceous eutherians and Laurasian origin for placental mammals near the K/T boundary. Nature. 2007, 447: 1003-1006. 10.1038/nature05854.
Ohta T: An examination of the generation time effect on molecular evolution. Proc Natl Acad Sci USA. 1993, 90: 10676-10680. 10.1073/pnas.90.22.10676.
Huchon D, Catzeflis FM, Douzery EJP: Molecular evolution of the nuclear von Willebrand Factor gene in mammals and the phylogeny of rodents. Mol Biol Evol. 1999, 16: 577-589.
Spradling TA, Hafner MS, Demastes JW: Differences in rate of Cytochrome-b evolution among species of rodents. Journal of Mammology. 2001, 82: 65-80. 10.1644/1545-1542(2001)082<0065:DIROCB>2.0.CO;2.
Rowe DL, Honeycutt RL: Phylogenetic relationships, ecological correlates, and molecular evolution within the Cavoidea (Mammalia, Rodentia). Mol Biol Evol. 2002, 19: 263-277.
Rottenberg H: Coevolution of exceptional longevity, exceptionally high metabolic rates, and mitochondrial DNA-coded proteins in mammals. Exp Gerontol. 2007, 42: 364-373. 10.1016/j.exger.2006.10.016.
Felsenstein J: Phylogenies and the comparative method. Am Nat. 1985, 125: 1-15. 10.1086/284325.
Harvey PH, Pagel MD: The comparative method in evolutionary biology. 1991, Oxford, UK: Oxford University Press
Garland TJ, Harvey PH, Ives AR: Procedures for the analysis of comparative data using phylogenetically independent contrasts. Syst Biol. 1992, 41: 18-32. 10.2307/2992503.
Freckleton RP: Phylogenetic tests of ecological and evolutionary hypotheses: checking for phylogenetic independence. Functional Ecology. 2000, 14: 129-134. 10.1046/j.1365-2435.2000.00400.x.
Welch JJ, Waxman D: Calculating independent contrasts for the comparative study of substitution rates. J Theor Biol. 2007, Epub ahead of print
Millar JS: Adaptive features of mammalian reproduction. Evolution. 1977, 31: 370-386. 10.2307/2407759.
Stearns SC: The evolution of life histories. 1992, Oxford, UK: Oxford University Press
Harvey PH, Purvis A: Understanding the ecological and evolutionary reasons for life history variation: mammals as a case study. Advanced ecological theory. Edited by: McGlade J. 1999, Oxford, UK: Blackwell Science Ltd, 232-248.
Speakman JR: Correlations between physiology and lifespan – two widely ignored problems with comparative studies. Aging Cell. 2005, 4: 167-175. 10.1111/j.1474-9726.2005.00162.x.
Seluanov A, Chen Z, Hine C, Sasahara THC, Ribeiro AACM, Catania KC, Presgraves DC, Gorbunova V: Telomerase activity coevolves with body mass, not lifespan. Aging cell. 2007, 6: 45-52. 10.1111/j.1474-9726.2006.00262.x.
Brown WM, Prager EM, Wang A, Wilson AC: Mitochondrial DNA sequences of primates: tempo and mode of evolution. J Mol Evol. 1982, 18: 225-239. 10.1007/BF01734101.
Savage VM, Gillooly JF, Woodruff WH, West GB, Allen AP, J EB, Brown JH: The predominance of quarter power scaling in biology. Functional Ecology. 2004, 18: 257-282. 10.1111/j.0269-8463.2004.00856.x.
Lanfear R, Thomas JA, Welch JJ, Bromham L: Metabolic Rate Does Not Calibrate the Molecular Clock. Proc Natl Acad Sci USA. 2007, 104: 15388-15393. 10.1073/pnas.0703359104.
Bromham L: Molecular Clocks in Reptiles: Life History Influences Rate of Molecular Evolution. Mol Biol Evol. 2002, 19: 302-309.
de Magalhães JP, Costa J, Church GM: An analysis of the relationship between metabolism, developmental schedules, and longevity using phylogenetic independent contrasts. J Gerontol Biol Sci Med Sci. 2007, 62 (2): 149-160.
Austad SN, Fischer KE: Mammalian aging, metabolism, and ecology: evidence from the bats and marsupials. Journal of Gerontology. 1991, 46: B47-B53.
Felsenstein J: Inferring Phylogenies. 2004, Sunderland Mass.: Sinauer Associates
Harman D: The biologic clock: the mitochondria?. J Am Geriatr Soc. 1972, 20: 145-147.
de Grey ADNJ: Mitochondria in homeotherm aging: Will detailed mechanisms consistent with the evidence now receive attention?. Aging Cell. 2004, 3: 77-10.1111/j.1474-9728.2004.00091.x.
Skulachev VP: Mitochondria, reactive oxygen species and longevity: some lessons from the Barja group. Aging Cell. 2004, 3: 17-19. 10.1111/j.1474-9728.2003.00076.x.
Fridovich I: Mitochondria: are they the seat of senescence?. Aging Cell. 2004, 3: 13-16. 10.1046/j.1474-9728.2003.00075.x.
Lee H-C, Wei Y-H: Oxidative stress, mitochondrial DNA mutation, and apoptosis in aging. Experimental biology and medicine. 2007, 232 (5): 592-606.
Krementz DG, Sauer JR, Nichols JD: Model-based estimates of annual survival rate are preferable to observed maximum lifespan statistics for use in comparative life history studies. Oikos. 1989, 56: 203-208. 10.2307/3565337.
Ricklefs RE, Scheuerlein A: Comparison of aging-related mortality among birds and mammals. Exp Gerontol. 2001, 36: 845-857. 10.1016/S0531-5565(00)00245-X.
Prothero J: Adult life span as a function of age at maturity. Exp Gerontol. 1993, 28: 529-536. 10.1016/0531-5565(93)90041-B.
Trifunovic A, Wredenberg A, Falkenberg M, Spelbrink JN, Rovio AT, Bruder CE, Bohooly YM, Gidlof S, Oldfors A, Wibom R, Törnell J, Jacobs HT, Larsson NG: Premature aging in mice expressing defective mitochondrial DNA polymerase. Nature. 2004, 429: 417-423. 10.1038/nature02517.
Kujoth GC, Hiona A, Pugh TD, Someya S, Panzer K, Wohlgemuth SE, Hofer T, Seo AY, Sullivan R, Jobling WA, Morrow JD, Van Remmen H, Sedivy JM, Yamasoba T, Tanokura M, Weindruch R, Leeuwenburgh C, Prolla TA: Mitochondrial DNA mutations, oxidative stress, and apoptosis in mammalian aging. Science. 2005, 309: 481-484. 10.1126/science.1112125.
Trifunovic A, Hansson A, Wredenberg A, Rovio AT, Durfour E, Khvorostov I, Spelbrink JN, Wibom R, Jacobs HT, Larsson NG: Somatic mtDNA mutations cause aging phenotypes without affecting reactive oxygen species production. Proc Natl Acad Sci USA. 2005, 102: 17993-17998. 10.1073/pnas.0508886102.
Schriner SE, Linford NJ, Martin GM, Treuting P, Ogburn CE, Emond M, Coskun PE, Ladiges W, Wolf N, Van Remmen H, Wallace DC, Rabinovitch PS: Extension of murine life span by overexpression of catalase targeted to mitochondria. Science. 2005, 308: 1909-1911. 10.1126/science.1106653.
Van Remmen H, Ikeno Y, Hamilton M, Pahlavani M, Wolf N, Thorpe SR, Alderson NL, Baynes JW, Epstein CJ, Huang TT, Nelson J, Strong R, Richardson A: Life-long reduction in MnSOD activity results in increased DNA damage and higher incidence of cancer but does not accelerate aging. Physiological Genomics. 2003, 16: 29-37. 10.1152/physiolgenomics.00122.2003.
Vermulst M, Bielas JH, Kujoth GC, Ladiges WC, Rabinovitch PS, Prolla TA, Loeb LA: Mitochondrial point mutations do not limit the natural lifespan of mice. Nat Genet. 2007, 39: 540-543. 10.1038/ng1988.
Rand DM: Thermal habit, metabolic rate and the evolution of mitochondrial DNA. Trends Ecol Evol. 1994, 9: 125-131. 10.1016/0169-5347(94)90176-7.
Speakman JR: Body size, energy metabolism and lifespan. J Exp Biol. 2005, 208: 1717-1730. 10.1242/jeb.01556.
Speakman JR, Talbot DA, Selman C, Snart S, McLaren JS, Redman P, Krol E, Jackson DM, Johnson MS, Brand MD: Uncoupled and surviving: individual mice with high metabolism have greater mitochondrial uncoupling and live longe. Aging Cell. 2004, 3: 87-95. 10.1111/j.1474-9728.2004.00097.x.
Hoffmann S, Spitkovsky D, Radicella JP, Epe B, Wiesner RJ: Reactive oxygen species derived from the mitochondrial respiratory chain are not responsible for the basal levels of oxidative base modifications observed in nuclear DNA of Mammalian cells. Free Radical Biology and Medicine. 2004, 36: 765-773. 10.1016/j.freeradbiomed.2003.12.019.
Andziak B, O'Connor TP, Qi W, DeWaal EM, Pierce A, Chaudhuri AR, Van Remmen H, Buffenstein R: High oxidative damage levels in the longest-living rodent, the naked mole-rat. Aging Cell. 2006, 5: 463-471. 10.1111/j.1474-9726.2006.00237.x.
Hulbert AJ, Faulks SC, Buffenstein R: Oxidation-resistant membrane phospholipids can explain longevity differences among the longest-living rodents and similarly-sized mice. J Geront Biol Sci Med Sci. 2006, 61 (10): 1009-1018.
Jones KE, MacLarnon A: Bat life histories: Testing models of mammalian life-history evolution. Evolutionary Ecology Research. 2001, 3: 465-476.
Gillespie JH: Is the population size of a species relevant to its evolution?. Evolution. 2001, 55 (11): 2161-2169.
Ballard JWO, Whitlock MC: The incomplete natural history of mitochondria. Mol Ecol. 2004, 13: 729-744. 10.1046/j.1365-294X.2003.02063.x.
Chao L, Carr DE: The molecular clock and the relationship between population size and generation time. Evolution. 1993, 47: 688-690. 10.2307/2410082.
Bazin E, Glémin S, Galtier N: Population size does not influence mitochondrial genetic diversity in animals. Science. 2006, 312: 570-571. 10.1126/science.1122033.
Mulligan CJ, Kitchen A, Miyamoto MM: Comment on "Population Size Does Not Influence Mitochondrial Genetic Diversity in Animals". Science. 2006, 314: 1390-10.1126/science.1132585.
Rhesus Macaque Genome Sequencing and Analysis Consortium: Evolutionary and biomedical insights from the rhesus macaque genome. Science. 2007, 316: 222-234. 10.1126/science.1139247.
Andressoo J-O, Hoeijmakers JHJ, Mitchell JR: Nucleotide excision repair disorders and the balance between cancer and aging. Cell Cycle. 2006, 5: 2886-2888.
Woodruff RC, Hual H, Thompson JN: Clusters of identical new mutation in the evolutionary landscape. Genetica. 1996, 98: 149-160. 10.1007/BF00121363.
Hedrick P: Large variance in reproductive success and the Ne/N ratio. Evolution. 2005, 59: 1596-1599.
Nunney L: The influence of variation in female fecundity on effective population size. Biological Journal of the Linnean Society. 1996, 59: 411-425.
National Center for Biotechnology Information. [http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/]
Rambaut A: Se-Al: Sequence Alignment Editor. 1996, [http://tree.bio.ed.ac.uk/software/seal]
de Magalhães JP, Costa J, Toussaint O: HAGR: the Human Ageing Genomic Resources. Nucleic Acids Research. 2005, 33: D537-D543. 10.1093/nar/gki017.
Wooton JT: The effects of body mass, phylogeny, habitat, and trophic level on mammalian age at first reproduction. Evolution. 1987, 41: 732-749. 10.2307/2408884.
Goldman N, Yang Z: A codon-based model of nucleotide substitution for protein-coding DNA sequences. Mol Biol Evol. 1994, 11: 725-736.
Yang Z: PAML: a program package for phylogenetic analysis by maximum likelihood. Computer Applications in BioSciences. 1997, 13 (5): 555-556.
R Development Core Team: R: a language and environment for statistical computing. V 2.4. 2006, [http://www.r-project.org]
We thank Marcel Cardillo, Andy Purvis, Kate Jones, and all members of the PanTHERIA Project for sharing their unpublished mammalian life-history data. Thanks are also due to Rob Lanfear, Fraser Lewis, Andrew Rambaut, Jess Thomas, David Waxman and an anonymous reviewer, for helpful advice and comments on the manuscript. JJW is supported by BBSRC grant DO17750 awarded to Andrew Rambaut. OBE was supported by a Heisenberg Fellowship of the DFG (BI 825/2-1).
LB and JW designed research project. JW and OBE compiled data. JW performed analyses. All authors wrote and approved the final manuscript.
Electronic supplementary material
About this article
Cite this article
Welch, J.J., Bininda-Emonds, O.R. & Bromham, L. Correlates of substitution rate variation in mammalian protein-coding sequences. BMC Evol Biol 8, 53 (2008) doi:10.1186/1471-2148-8-53
- Substitution Rate
- Basal Metabolic Rate
- Maximum Lifespan
- Synonymous Rate
- Substitution Rate Variation