Molecular clocks in reptiles: life history influences rate of molecular evolution

Life history has been implicated as a determinant of variation in rate of molecular evolution amongst vertebrate species because of a negative correlation between body size and substitution rate for many molecular data sets. Both the generality and the cause of the negative body size trend have been debated, and the validity of key studies has been questioned (particularly concerning the failure to account for phylogenetic bias). In this study, a comparative method has been used to test for an association between a range of life-history variables-such as body size, age at maturity, and clutch size-and DNA substitution rate for three genes (NADH4, cytochrome b, and c-mos). A negative relationship between body size and rate of molecular evolution was found for phylogenetically independent pairs of reptile species spanning turtles, lizards, snakes, crocodile, and tuatara. Although this study was limited by the number of comparisons for which both sequence and life-history data were available, the results suggest that a negative body size trend in rate of molecular evolution may be a general feature of reptile molecular evolution, consistent with similar studies of mammals and birds. This observation has important implications for uncovering the mechanisms of molecular evolution and warns against assuming that related lineages will share the same substitution rate (a local molecular clock) in order to date evolutionary divergences from DNA sequences.     

Body size effects and rates of cytochrome b evolution in tube-nosed seabirds

Variation in rates of molecular evolution now appears to be widespread. The demonstration that body size is correlated with rates of molecular evolution suggests that physiological and ecological factors may be involved in molecular rate variation, but large-scale comparative studies are still lacking. Here, we use complete cytochrome b sequences from 85 species of tube-nosed seabirds (order Procellariiformes) and 5 outgroup species of penguins (order Sphenisciformes) to test for an association between body mass and rates of molecular evolution within the former avian order. Cladistic analysis of the 90 sequences estimates a phylogeny largely consistent with the traditional taxonomy of the Procellariiformes. The Diomedeidae, Procellariidae, and Pelecanoididae are monophyletic, while the Hydrobatidae are basal and paraphyletic. However, the two subfamilies within the Hydrobatidae (Hydrobatinae and Oceanitinae) are monophyletic. A likelihood ratio test detects significant deviation from clocklike evolution in our data. Using a sign test for an association between body mass and branch length in the seabird phylogeny, we find that larger taxa tend to have shorter terminal branch lengths than smaller taxa. This observation suggests that rates of mitochondrial DNA evolution are slower for larger taxa. Rate calibrations based on the fossil record reveal concordant body size effects. We interpret these results as evidence for a metabolic rate effect, as the species in this order exhibit large differences in metabolic rates, which are known to be highly correlated with body mass in this group. Our results support previous findings of body size effects and show that this effect can be significant even within a single avian order. This suggests that even lineage-specific molecular clocks may not be tenable if calibrations involve taxa with different metabolic rates.      

The Relationship Between the Rate of Molecular Evolution and the Rate of Genome Rearrangement in Animal Mitochondrial Genomes

Evolution of mitochondrial genes is far from clock-like. The substitution rate varies considerably between species, and there are many species that have a significantly increased rate with respect to their close relatives. There is also considerable variation among species in the rate of gene order rearrangement. Using a set of 55 complete arthropod mitochondrial genomes, we estimate the evolutionary distance from the common ancestor to each species using protein sequences, tRNA sequences, and breakpoint distances (a measure of the degree of genome rearrangement). All these distance measures are correlated. We use relative rate tests to compare pairs of related species in several animal phyla. In the majority of cases, the species with the more highly rearranged genome also has a significantly higher rate of sequence evolution. Species with higher amino acid substitution rates in mitochondria also have more variable amino acid composition in response to mutation pressure. We discuss the possible causes of variation in rates of sequence evolution and gene rearrangement among species and the possible reasons for the observed correlation between the two rates. [Reviewing Editor: Dr. David Pollock]     

The determinants of the molecular substitution process in turtles

Neutral rates of molecular evolution vary across species, and this variation has been shown to be related to biological traits. One of the first patterns to be observed in vertebrates has been an inverse relationship between body mass (BM) and substitution rates. The effects of three major life‐history traits (LHT) that covary with BM – metabolic rate, generation time and longevity (LON) – have been invoked to explain this relationship. However, most of the theoretical and empirical evidence supporting this relationship comes from endothermic vertebrates, that is, mammals and birds, in which the environmental conditions, especially temperature, do not have a direct impact on cellular and molecular biology. We analysed the variations in mitochondrial and nuclear rates of synonymous substitution across 224 turtle species and examined their correlation with two LHT (LON and BM) and two environmental variables [latitude (LAT) and habitat]. Our analyses indicate that in turtles, neutral rates of molecular evolution are hardly correlated with LON or BM. Rather, both the mitochondrial and nuclear substitution rates are significantly correlated with LAT – faster evolution in the tropics – and especially so for aquatic species. These results question the generality of the relationships reported in mammals and birds and suggest that environmental factors might be the strongest determinants of the mutation rate in ectotherms.     

Consistent variation in amino-acid substitution rate, despite uniformity of mutation rate: protein evolution in mammals is not neutral

Variation in mutation rate, attributed to differences in both generation time and in metabolic rate, has been invoked under the neutral theory of molecular evolution to account for differences in substitution rate among mammalian lineages. We show that substitution rates at fourfold-degenerate sites and at sites in noncoding regions do not vary between the primate and rodent lineages, implying mutation- rate uniformity. In contrast, the substitution rates at nondegenerate sites vary both within and between lineages. This difference in substitution-rate pattern between the two types of site is incompatible with neutral theory but may result from substitutions occurring by fixation of slightly deleterious mutations. Variation in the rate of protein evolution among mammalian lineages appears to be due more to differences in population fixation rates than to biochemical or physiological differences affecting mutation rates.      

Unequal evolutionary rates between annual and perennial lineages of checker mallows (Sidalcea, Malvaceae): evidence from 18S-26S rDNA internal and external transcribed spacers

Heterogeneous DNA substitution rates were found in the 18S-26S nuclear ribosomal DNA internal transcribed spacer (ITS) and external transcribed spacer (ETS) regions of Sidalcea (Malvaceae), a putatively young genus of annuals and perennials. The majority of comparisons revealed that the annual species had significantly higher molecular evolutionary rates than the perennials, whereas rates were consistently homogenous between obligate annual species. These findings led us to conclude that generation time or possibly another biological factor distinguishing annuals and perennials has influenced rates of molecular evolution in SIDALCEA: The congruence of relative-rate test results across both spacer regions reinforced the association between life history and rate of rDNA evolution across lineages of checker mallows. Evolutionary rate variation within perennials mainly involved three basally divergent lineages. The faster rate in one lineage, Sidalcea stipularis, compared with other perennials may be the result of genetic drift in the only known, small, population. The other two basally divergent lineages had slower evolutionary rates compared with the remaining perennials; possible explanations for these differences include rate-reducing effects of a suffrutescent (rather than herbaceous) habit and seed dormancy.     

Maximum-likelihood models for combined analyses of multiple sequence data

Models of nucleotide substitution were constructed for combined analyses of heterogeneous sequence data (such as those of multiple genes) from the same set of species. The models account for different aspects of the heterogeneity in the evolutionary process of different genes, such as differences in nucleotide frequencies, in substitution rate bias (for example, the transition/transversion rate bias), and in the extent of rate variation across sites. Model parameters were estimated by maximum likelihood and the likelihood ratio test was used to test hypotheses concerning sequence evolution, such as rate constancy among lineages (the assumption of a molecular clock) and proportionality of branch lengths for different genes. The example data from a segment of the mitochondrial genome of six hominoid species (human, common and pygmy chimpanzees, gorilla, orangutan, and siamang) were analyzed. Nucleotides at the three codon positions in the protein-coding regions and from the tRNA-coding regions were considered heterogeneous data sets. Statistical tests showed that the amount of evolution in the sequence data reflected in the estimated branch lengths can be explained by the codon-position effect and lineage effect of substitution rates. The assumption of a molecular clock could not be rejected when the data were analyzed separately or when the rate variation among sites was ignored. However, significant differences in substitution rate among lineages were found when the data sets were combined and when the rate variation among sites was accounted for in the models. Under the assumption that the orangutan and African apes diverged 13 million years ago, the combined analysis of the sequence data estimated the times for the human-chimpanzee separation and for the separation of the gorilla as 4.3 and 6.8 million years ago, respectively.     

You are where you live: parasitic nematode mitochondrial genome size is associated with the thermal environment generated by hosts

There exists remarkable interspecific variation in mitochondrial sequence evolution rates and in mitochondrial genome sizes. A number of hypotheses based on the forces of mutation and selection have been proposed to explain this variation. Among such hypotheses, we test three: 1) the ‘longevity‐dependent selection’, 2) the ‘functional constraints’ and 3) the ‘race for replication’ hypotheses, using published mtDNA genomic sequences of 47 Nematoda species. We did not find any relationship between body size (used as a proxy for longevity) and genome size or the substitution rate of protein sequences, providing little evidence for the first hypothesis. Parasitic species from different thermal habitats, as determined by their definitive host type (ectothermal vs. endothermal), did not differ in their rates of protein evolution. Therefore, little support was obtained for the second hypothesis. However, we revealed that mitogenomes of parasites of endotherms were significantly smaller than those of parasites of ectotherms, supporting the race for replication hypothesis. As mitochondrial genomes of endothermal animals are usually more compact than those of ectothermal animals, intriguingly, nematode parasites of endotherms and ectotherms exhibit similar patterns of mtDNA length variation to their hosts.     

Multiple nuclear insertions of mitochondrial cytochrome b sequences in callitrichine primates

We report the presence of four nuclear paralogs of a 380-bp segment of cytochrome b in callitrichine primates (marmosets and tamarins). The mitochondrial cytochrome b sequence and each nuclear paralog were obtained from several species, allowing multiple comparisons of rates and patterns of substitution both between mitochondrial and nuclear sequences and among nuclear sequences. The mitochondrial DNA had high overall rates of molecular evolution and a strong bias toward substitutions at third codon positions. Rates of molecular evolution among the nuclear sequences were low and constant, and there were small differences in substitution patterns among the nuclear clades which were probably attributable to the small number of sites involved. A novel method of phylogenetic reconstruction based on the large difference in rates of evolution at different codon positions among mitochondrial and nuclear clades was used to determine whether different nuclear paralogs represent independent transposition events or duplications following a single insertion. This method is generally applicable in cases where differences in pattern of molecular evolution are known, and it showed that at least three of the four nuclear clades represent independent transposition events. The insertion events giving rise to two of the nuclear clades predate the divergence of the callitrichines, whereas those leading to the other two nuclear clades may have occurred in the common ancestor of marmosets.     

Molecular phylogenetics and ecological diversification of the transisthmian fish genus Centropomus (Perciformes: Centropomidae).

Phylogenetic relationships among the 12 recognized fish species in the New World genus Centropomus (Pisces, Centropomidae) were analyzed using allozyme electrophoresis and 618 bp of the mitochondrial DNA 16S ribosomal RNA (rRNA) gene. Molecular phylogenetic trees were generally consistent with previously published partial hypotheses based on morphological evidence. However, previously undefined sister group relationships between major species groups were resolved using molecular data, and phylogenetic hypotheses for Centropomus based on 16S rRNA sequences were better supported than were allozyme-based hypotheses. The high level of congruence among the trees inferred from the nuclear and mitochondrial characters provided a firm phylogenetic basis for analysis of ecological diversification and molecular evolution in the genus. Compared to basal Centropomus species, members of the most nested species group were significantly larger in body size and occupied a marine niche only peripherally utilized by their congeners. We also observed substitution rate heterogeneity among 16S rRNA lineages; in contrast to expectations based on "metabolic rate" and "generation interval" models, relative substitution rates were faster than expected for the group of large-bodied snooks. Using the Pliocene rise of the Central American isthmian marine barrier to calibrate rates of 16S ribosomal gene evolution in Centropomus, we found that the rates for the genus were similar to those reported for higher vertebrates. Analysis of the three sets of transisthmian geminate taxa in Centropomus indicated that two of the pairs were probably formed during the Pliocene rise of the isthmus while the third pair diverged several million years earlier.     

An Evolutionary Footprint of Age-Related Natural Selection in Mitochondrial DNA

By comparing mtDNA sequences between different orders of mammals, we show that both longevity and generation time are significantly correlated with the nucleotide content of the mtDNA. Specifically, there is a positive correlation between generation time and mt GC content. This correlation is repeated, at a finer evolutionary scale, within the primates. Moreover, a comparison of human and chimpanzee mtDNAs shows that the effect has been very pronounced during the short evolutionary period since the divergence of these two species, with human mtDNA showing a GC-biased pattern of substitution at the variable sites. In addition to these DNA sequence patterns, comparisons between the human and the chimp mt protein sequences also revealed a surprisingly high substitution rate for threonine residues, resulting in a reduction of threonine in the human mt proteome. These patterns of both DNA and protein evolution can be explained by a balance between AT-biased mutational pressure and age-related purifying selection.     

Ascomycota has a faster evolutionary rate and higher species diversity than Basidiomycota

Differences in rates of nucleotide or amino acid substitutions among major groups of organisms are repeatedly found and well documented. A growing body of evidence suggests a link between the rate of neutral molecular change within populations and the evolution of species diversity. More than 98% of terrestrial fungi belong to the phyla Ascomycota or Basidiomycota. The former is considerably richer in number of species than the latter. We obtained DNA sequences of 21 protein-coding genes from the lichenized fungus Rhizoplaca chrysoleuca and used them together with sequences from GenBank for subsequent analyses. Three datasets were used to test rate discrepancies between Ascomycota and Basidiomycota and that within Ascomycota: (i) 13 taxa including 105 protein-coding genes, (ii) nine taxa including 21 protein-coding genes, and (iii) nuclear LSU rDNA of 299 fungal species. Based on analyses of the 105 protein-coding genes and nuclear LSU rDNA datasets, we found that the evolutionary rate was higher in Ascomycota than in Basidiomycota. The differences in substitution rates between Ascomycota and Basidiomycota were significant. Within Ascomycota, the species-rich Sordariomycetes has the fastest evolutionary rate, while Leotiomycetes has the slowest. Our results indicate that the main contribution to the higher substitution rates in Ascomycota does not come from mutualism, ecological conditions, sterility, metabolic rate or shorter generation time, but is possibly caused by the founder effect. This is another example of the correlation between species number and evolutionary rates, which is consistent with the hypothesis that the founder effect is responsible for accelerated substitution rates in diverse clades.     

Scombroid Fishes Provide Novel Insights into the Trait/Rate Associations of Molecular Evolution

The study of which life history traits primarily affect molecular evolutionary rates is often confounded by the covariance of these traits. Scombroid fishes (billfishes, tunas, barracudas, and their relatives) are unusual in that their mass-specific metabolic rate is positively associated with body size. This study exploits this atypical pattern of trait variation, which allows for direct tests of whether mass-specific metabolic rate or body size is the more important factor of molecular evolutionary rates. We inferred a phylogeny for scombroids from a supermatrix of molecular and morphological characters and used new phylogenetic comparative approaches to assess the associations of body size and mass-specific metabolic rate with substitution rate. As predicted by the body size hypothesis, there is a negative correlation between body size and substitution rate. However, unexpectedly, we also find a negative association between mass-specific metabolic and substitution rates. These relationships are supported by analyses of the total molecular data, separate mitochondrial and nuclear genes, and individual loci, and they are robust to phylogenetic uncertainty. The molecular evolutionary rates of scombroids are primarily tied to body size. This study demonstrates that groups with novel patterns of trait variation can be particularly informative for identifying which life history traits are the primary factors of molecular evolutionary rates.     

Rates of molecular evolution and the fraction of nucleotide positions free to vary

Summary Selective constraints on DNA sequence change were incorporated into a model of DNA divergence by restricting substitutions to a subset of nucleotide positions. A simple model showed that both mutation rate and the fraction of nucleotide positions free to vary are strong determinants of DNA divergence over time.When divergence between two species approaches the fraction of positions free to vary, standard methods that correct for multiple mutations yield severe underestimates of the number of substitutions per site. A modified method appropriate for use with DNA sequence, restriction site, or thermal renaturation data is derived taking this fraction into account. The model also showed that the ratio of divergence in two gene classes (e.g., nuclear and mitochondrial) may vary widely over time even if the ratio of mutation rates remains constant.DNA sequence divergence data are used increasingly to detect differences in rates of molecular evolution. Often, variation in divergence rate is assumed to represent variation in mutation rate. The present model suggests that differing divergence rates among comparisons (either among gene classes or taxa) should be interpreted cautiously. Differences in the fraction of nucleotide positions free to vary can serve as an important alternative hypothesis to explain differences in DNA divergence rates.     

Exploring the correlations between sequence evolution rate and phenotypic divergence across the Mammalian tree provides insights into adaptive evolution

Sequence evolution behaves in a relatively consistent manner, leading to one of the fundamental paradigms in biology, the existence of a ??molecular clock??. The molecular clock can be distilled to the concept of accumulation of substitutions, through time yielding a stable rate from which we can estimate lineage divergence. Over the last 50?years, evolutionary biologists have obtained an in-depth understanding of this clock??s nuances. It has been fine-tuned by taking into account the vast heterogeneity in rates across lineages and genes, leading to ??relaxed?? molecular clock methods for timetree reconstruction. Sequence rate varies with life history traits including body size, generation time and metabolic rate, and we review recent studies on this topic. However, few studies have explicitly examined correlates between molecular evolution and morphological evolution. The patterns observed across diverse lineages suggest that rates of molecular and morphological evolution are largely decoupled. We discuss how identifying the molecular mechanisms behind rapid functional radiations are central to understanding evolution. The vast functional divergence within mammalian lineages that have relatively ??slow?? sequence evolution refutes the hypotheses that pulses in diversification yielding major phenotypic change are the result of steady accumulation of substitutions. Patterns rather suggest phenotypic divergence is likely caused by regulatory alterations mediated through mechanisms such as insertions/deletions in functional regions. These can rapidly arise and sweep to fixation faster than predicted from a lineage??s sequence neutral substitution rate, enabling species to leapfrog between phenotypic ??islands??. We suggest research directions that could illuminate mechanisms behind the functional diversity we see today.     

Evolution of the mitochondrial rps3 intron in perennial and annual angiosperms and homology to nad5 intron 1.

The plant mitochondrial rps3 intron was analyzed for substitution and indel rate variation among 15 monocot and dicot angiosperms from 10 genera, including perennial and annual taxa. Overall, the intron sequence was very conserved among angiosperms. Based on length polymorphism, 10 different alleles were identified among the 10 genera. These allelic differences were mainly attributable to large indels. An insertion of 133 nucleotides, observed in the Alnus intron was partially or completely absent in the other lineages of the family Betulaceae. This insertion was located within domain IV of the secondary-structure model of this group IIA intron. A mobile element of 47 nucleotides that showed homology to sequences located in rice rps3 intron and in intergenic plant mitochondrial genomes was found within this insertion. Both substitution and indel rates were low among the Betulaceae sequences, but substitution rates were increasingly larger than indel rates in comparisons involving more distantly related taxa. From a secondary-structure model, regions involved in helical structures were shown to be well preserved from indels as compared to substitutions, but compensatory changes were not observed among the angiosperm sequences analyzed. Using approximate divergence times based on the fossil record, substitution and indel rate heterogeneity was observed between different pairs of annual and perennial taxa. In particular, the annual petunia and primrose evolved more than 15 and 10 times faster, for substitution and indel rates respectively, than the perennial birch and alder. This is the first demonstration of an evolutionary rate difference between perennial and annual forms in noncoding DNA, lending support to neutral causes such as the generation time, population size, and speciation rate effects to explain such rate heterogeneity. Surprisingly, the sequence from the rps3 intron had a high identity with the sequence of intron 1 from the angiosperm mitochondrial nad5 gene, suggesting a common origin of these two group IIA introns.     

Fast molecular evolution associated with high active metabolic rates in poison frogs

Molecular evolution is simultaneously paced by mutation rate, genetic drift, and natural selection. Life history traits also affect the speed of accumulation of nucleotide changes. For instance, small body size, rapid generation time, production of reactive oxygen species (ROS), and high resting metabolic rate (RMR) are suggested to be associated with faster rates of molecular evolution. However, phylogenetic correlation analyses failed to support a relationship between RMR and molecular evolution in ectotherms. In addition, RMR might underestimate the metabolic budget (e.g., digestion, reproduction, or escaping predation). An alternative is to test other metabolic rates, such as active metabolic rate (AMR), and their association with molecular evolution. Here, I present comparative analyses of the associations between life history traits (i.e., AMR, RMR, body mass, and fecundity) with rates of molecular evolution of and mitochondrial loci from a large ectotherm clade, the poison frogs (Dendrobatidae). My results support a strong positive association between mass-specific AMR and rates of molecular evolution for both mitochondrial and nuclear loci. In addition, I found weaker and genome-specific covariates such as body mass and fecundity for mitochondrial and nuclear loci, respectively. No direct association was found between mass-specific RMR and rates of molecular evolution. Thus, I provide a mechanistic hypothesis of the link between AMRs and the rate of molecular evolution based on an increase in ROS within germ line cells during periodic bouts of hypoxia/hyperoxia related to aerobic exercise. Finally, I propose a multifactorial model that includes AMR as a predictor of the rate of molecular evolution in ectothermic lineages.     

The genome as a life-history character: why rate of molecular evolution varies between mammal species

DNA sequences evolve at different rates in different species. This rate variation has been most closely examined in mammals, revealing a large number of characteristics that can shape the rate of molecular evolution. Many of these traits are part of the mammalian life-history continuum: species with small body size, rapid generation turnover, high fecundity and short lifespans tend to have faster rates of molecular evolution. In addition, rate of molecular evolution in mammals might be influenced by behaviour (such as mating system), ecological factors (such as range restriction) and evolutionary history (such as diversification rate). I discuss the evidence for these patterns of rate variation, and the possible explanations of these correlations. I also consider the impact of these systematic patterns of rate variation on the reliability of the molecular date estimates that have been used to suggest a Cretaceous radiation of modern mammals, before the final extinction of the dinosaurs.     

Why do species vary in their rate of molecular evolution?

Despite hopes that the processes of molecular evolution would be simple, clock-like and essentially universal, variation in the rate of molecular evolution is manifest at all levels of biological organization. Furthermore, it has become clear that rate variation has a systematic component: rate of molecular evolution can vary consistently with species body size, population dynamics, lifestyle and location. This suggests that the rate of molecular evolution should be considered part of life-history variation between species, which must be taken into account when interpreting DNA sequence differences between lineages. Uncovering the causes and correlates of rate variation may allow the development of new biologically motivated models of molecular evolution that may improve bioinformatic and phylogenetic analyses.     

Rates of molecular evolution in tree ferns are associated with body size,environmental temperature,and biological productivity

Variation in rates of molecular evolution (heterotachy) is a common phenomenon among plants. Although multiple theoretical models have been proposed, fundamental questions remain regarding the combined effects of ecological and morphological traits on rate heterogeneity. Here, we used tree ferns to explore the correlation between rates of molecular evolution in chloroplast DNA sequences and several morphological and environmental factors within a Bayesian framework. We revealed direct and indirect effects of body size, biological productivity, and temperature on substitution rates, where smaller tree ferns living in warmer and less productive environments tend to have faster rates of molecular evolution. In addition, we found that variation in the ratio of nonsynonymous to synonymous substitution rates (dN/dS) in the chloroplast rbcL gene was significantly correlated with ecological and morphological variables. Heterotachy in tree ferns may be influenced by effective population size associated with variation in body size and productivity. Macroevolutionary hypotheses should go beyond explaining heterotachy in terms of mutation rates and instead, should integrate population‐level factors to better understand the processes affecting the tempo of evolution at the molecular level.