Protein evolution in different cellular environments: cytochrome b in sharks and mammals

DNA sequences for the mitochondrial cytochrome b gene were determined for 13 species of sharks. Rates and patterns of amino acid replacement are compared for sharks and mammals. Absolute rates of cytochrome b evolution are six times slower in sharks than in mammals. Bivariate plots of the number of nonsynonymous and silent transversions are indistinguishable in the two groups, however, suggesting that the differences in amino acid replacement rates are due primarily to differences in DNA substitution rates. Patterns of amino acid replacement are also similar in the two groups. Conserved and variable regions occur in the same parts of the cytochrome b gene, and there is little evidence that the types of amino acid changes are significantly different between the groups. Similarity in the relative rates and patterns of protein change between the two groups prevails despite dramatic differences in the cellular environments of sharks and mammals. Poor penetrance of physiological differences through to rates of protein evolution provides support for the neutral theory and suggests that, for cytochrome b, patterns of evolution have been relatively constant throughout much of vertebrate history.      

Metabolic rate and directional nucleotide substitution in animal mitochondrial DNA

There is marked heterogeneity of nucleotide composition in mitochondrial DNA across divergent animals. Differences in nucleotide composition presumably reflect differences in directional nucleotide substitution for A+T or G+C nucleotides. In mitochondrial DNA, there is A+T directional nucleotide substitution in most (if not all) animals surveyed, and the magnitude of directional A+T nucleotide substitution differs greatly within and among groups. Differences in directional nucleotide substitution among lineages of mammals can be explained by changes in metabolic physiology. This relationship is thought to be mediated by the effect of oxygen radicals because these toxic compounds are by-products of aerobic metabolism and are known mutagens. Association between metabolism and nucleotide composition provides additional evidence in favor of the hypothesis that rates and patterns of nucleotide substitution in mitochondrial DNA can be influenced by factors that impinge on rates of endogenous DNA damage.      

Estimating hybridization in the presence of coalescence using phylogenetic intraspecific sampling

Background

Recent research has indicated a positive association between rates of molecular evolution and diversification in a number of taxa. However debate continues concerning the universality and cause of this relationship. Here, we present the first systematic investigation of this relationship within the mammals. We use phylogenetically independent sister-pair comparisons to test for a relationship between substitution rates and clade size at a number of taxonomic levels. Total, non-synonymous and synonymous substitution rates were estimated from mitochondrial and nuclear DNA sequences.

Results

We found no evidence for an association between clade size and substitution rates in mammals, for either the nuclear or the mitochondrial sequences. We found significant associations between body size and substitution rates, as previously reported.

Conclusions

Our results present a contrast to previous research, which has reported significant positive associations between substitution rates and diversification for birds, angiosperms and reptiles. There are three possible reasons for the differences between the observed results in mammals versus other clades. First, there may be no link between substitution rates and diversification in mammals. Second, this link may exist, but may be much weaker in mammals than in other clades. Third, the link between substitution rates and diversification may exist in mammals, but may be confounded by other variables.     

Sociality and the rate of molecular evolution

The molecular clock does not tick at a uniform rate in all taxa but may be influenced by species characteristics. Eusocial species (those with reproductive division of labor) have been predicted to have faster rates of molecular evolution than their nonsocial relatives because of greatly reduced effective population size; if most individuals in a population are nonreproductive and only one or few queens produce all the offspring, then eusocial animals could have much lower effective population sizes than their solitary relatives, which should increase the rate of substitution of "nearly neutral" mutations. An earlier study reported faster rates in eusocial honeybees and vespid wasps but failed to correct for phylogenetic nonindependence or to distinguish between potential causes of rate variation. Because sociality has evolved independently in many different lineages, it is possible to conduct a more wide-ranging study to test the generality of the relationship. We have conducted a comparative analysis of 25 phylogenetically independent pairs of social lineages and their nonsocial relatives, including bees, wasps, ants, termites, shrimps, and mole rats, using a range of available DNA sequences (mitochondrial and nuclear DNA coding for proteins and RNAs, and nontranslated sequences). By including a wide range of social taxa, we were able to test whether there is a general influence of sociality on rates of molecular evolution and to test specific predictions of the hypothesis: (1) that social species have faster rates because they have reduced effective population sizes; (2) that mitochondrial genes would show a greater effect of sociality than nuclear genes; and (3) that rates of molecular evolution should be correlated with the degree of sociality. We find no consistent pattern in rates of molecular evolution between social and nonsocial lineages and no evidence that mitochondrial genes show faster rates in social taxa. However, we show that the most highly eusocial Hymenoptera do have faster rates than their nonsocial relatives. We also find that social parasites (that utilize the workers from related species to produce their own offspring) have faster rates than their social relatives, which is consistent with an effect of lower effective population size on rate of molecular evolution. Our results illustrate the importance of allowing for phylogenetic nonindependence when conducting investigations of determinants of variation in rate of molecular evolution.     

Phylogeny of elasmobranchs based on LSU and SSU ribosomal RNA genes

The dominant view of the phylogeny of living elasmobranchs, based on morphological characters, is that batoids (skates and rays) are derived sharks, joined with saw sharks, and angel sharks in the clade Hypnosqualea [S. Shirai, Squalean Phylogeny: A New Framework of 'Squaloid' Sharks and Related Taxa, Hokkaido University Press, Sapporo, 1992]. By contrast, a recent molecular-phylogenetic study based on mitochondrial genes for 12S and 16S rRNA and tRNA valine [C.J. Douady et al., Mol. Phylogenet. Evol., 26 (2003) 215-221] supported the older view that batoids and sharks are separate lineages. Here, we tested these two different views using combined, nuclear large-subunit and small-subunit rRNA gene sequences ( approximately 5.3kb) from 22 elasmobranchs, two chimeras, and two bony fishes. We used maximum likelihood, maximum parsimony, minimum evolution, and Bayesian inference for tree reconstruction, and found the large-subunit rRNA gene to contain far more signal than the small-subunit gene for resolving this mostly Mesozoic radiation. Our findings matched those of in separating batoids from sharks and in statistically rejecting Hypnosqualea. The angel shark (Squatina) was the sister group to squaliforms (dogfish sharks), and our findings are consistent with the idea that "orbitostylic" sharks form a monophyletic group (squaliforms+the hexanchiform Chlamydoselachus+Squatina+Pristiophorus). In the galeomorph sharks, however, lamniforms grouped with orectolobiforms, opposing the widely accepted 'lamniform+carcharhiniform' grouping. A tree based on the mitochondrial gene for cytochrome b also supported a separation of sharks and batoids, in contrast to Hypnosqualea. Among elasmobranchs, variation in the evolutionary rates of the nuclear rRNA genes was higher than that of cytochrome b genes, mainly due to the relatively rapid evolution of rRNA in some carcharhiniforms. In conclusion, several different molecular studies now refute the Hypnosqualea hypothesis of elasmobranch interrelationships.     

Accelerated evolution of mitochondrial but not nuclear genomes of Hymenoptera: new evidence from crabronid wasps

Mitochondrial genes in animals are especially useful as molecular markers for the reconstruction of phylogenies among closely related taxa, due to the generally high substitution rates. Several insect orders, notably Hymenoptera and Phthiraptera, show exceptionally high rates of mitochondrial molecular evolution, which has been attributed to the parasitic lifestyle of current or ancestral members of these taxa. Parasitism has been hypothesized to entail frequent population bottlenecks that increase rates of molecular evolution by reducing the efficiency of purifying selection. This effect should result in elevated substitution rates of both nuclear and mitochondrial genes, but to date no extensive comparative study has tested this hypothesis in insects. Here we report the mitochondrial genome of a crabronid wasp, the European beewolf (Philanthus triangulum, Hymenoptera, Crabronidae), and we use it to compare evolutionary rates among the four largest holometabolous insect orders (Coleoptera, Diptera, Hymenoptera, Lepidoptera) based on phylogenies reconstructed with whole mitochondrial genomes as well as four single-copy nuclear genes (18S rRNA, arginine kinase, wingless, phosphoenolpyruvate carboxykinase). The mt-genome of P. triangulum is 16,029 bp in size with a mean A+T content of 83.6%, and it encodes the 37 genes typically found in arthropod mt genomes (13 protein-coding, 22 tRNA, and two rRNA genes). Five translocations of tRNA genes were discovered relative to the putative ancestral genome arrangement in insects, and the unusual start codon TTG was predicted for cox2. Phylogenetic analyses revealed significantly longer branches leading to the apocritan Hymenoptera as well as the Orussoidea, to a lesser extent the Cephoidea, and, possibly, the Tenthredinoidea than any of the other holometabolous insect orders for all mitochondrial but none of the four nuclear genes tested. Thus, our results suggest that the ancestral parasitic lifestyle of Apocrita is unlikely to be the major cause for the elevated substitution rates observed in hymenopteran mitochondrial genomes.     

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.     

Mitochondrial genome sequence evolution in Chlamydomonas

The mitochondrial genomes of the Chlorophyta exhibit significant diversity with respect to gene content and genome compactness; however, quantitative data on the rates of nucleotide substitution in mitochondrial DNA, which might help explain the origin of this diversity, are lacking. To gain insight into the evolutionary forces responsible for mitochondrial genome diversification, we sequenced to near completion the mitochondrial genome of the chlorophyte Chlamydomonas incerta, estimated the evolutionary divergence between Chlamydomonas reinhardtii and C. incerta mitochondrial protein-coding genes and rRNA-coding regions, and compared the relative evolutionary rates in mitochondrial and nuclear genes. Synonymous and nonsynonymous substitution rates do not differ significantly between the mitochondrial and nuclear protein-coding genes. The mitochondrial rRNA-coding regions, however, are evolving much faster than their nuclear counterparts, and this difference might be explained by relaxed functional constraints on the mitochondrial translational apparatus due to the small number of proteins synthesized in Chlamydomonas mitochondria. Substitution rates at synonymous sites in a nonstandard mitochondrial gene (rtl) and at intronic and synonymous sites in nuclear genes expressed at low levels suggest that the mutation rate is similar in these two genetic compartments. Potential evolutionary forces shaping mitochondrial genome evolution in Chlamydomonas are discussed.     

Maximum-Likelihood Analysis of the Tethythere Hypothesis Based on a Multigene Data Set and a Comparison of Different Models of Sequence Evolution

Hyracoids have been allied with either perissodactyls or tethytheres (i.e., Proboscidea + Sirenia) based on morphological data. The latter hypothesis, termed Paenungulata, is corroborated by numerous molecular studies. However, molecular studies have failed to support Tethytheria, a group that is supported by morphological data. We examined relationships among living paenungulate orders using a multigene data set that included sequences from four mitochondrial genes (12S rRNA, tRNA valine, 16S rRNA, cytochrome b) and four nuclear genes (aquaporin, A2AB, IRBP, vWF). Nineteen maximum-likelihood models were employed, including models with process partitions for base composition and substitution parameterizations. With the inclusion of partitions with a heterogeneous base composition, 18 of 19 models favored Hyracoidea + Sirenia. All 19 models favored Hyracoidea + Sirenia after excluding heterogeneous base composition partitions. Most of the support for Hyracoidea + Sirenia derived from the mitochondrial genes (bootstrap support ranged from 51 to 99%); Tethytheria, in turn, received 0 to 19% support in different analyses. Bootstrap support deriving from the nuclear genes was more evenly split among the competing hypotheses (3 to 45% for Tethytheria; 17.5 to 62% for Hyracoidea + Sirenia). Lineage-specific rate variation among both mitochondrial and nuclear genes may contribute to the different results that were obtained with mitochondrial versus nuclear data. Whether Tethytheria or a competing hypothesis is correct, short internodes on the molecular phylogenies suggest that paenungulate orders diverged from each other over a 5- to 8-million-year time window extending from the late Paleocene into the early Eocene. We also used likelihood-ratio tests to compare different models of sequence evolution. A gamma distribution of rates results in a greater improvement in likelihood scores than does an allowance for invariant sites. Twenty-one rate partitions corresponding to stems, loops, and codon positions of different genes result in higher likelihood scores than a gamma distribution of rates and/or an allowance for invariant sites. Process partitions of the data that incorporate base composition and substitution parameterizations result in significant improvements in likelihood scores in comparison to models that allow only for relative rate differences among partitions.     

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.     

Likelihood analysis of asymmetrical mutation bias gradients in vertebrate mitochondrial genomes

Protein-coding genes in mitochondrial genomes have varying degrees of asymmetric skew in base frequencies at the third codon position. The variation in skew among genes appears to be caused by varying durations of time that the heavy strand spends in the mutagenic single-strand state during replication (D(ssH)). The primary data used to study skew have been the gene-by-gene base frequencies in individual taxa, which provide little information on exactly what kinds of mutations are responsible for the base frequency skew. To assess the contribution of individual mutation components to the ancestral vertebrate substitution pattern, here we analyze a large data set of complete vertebrate mitochondrial genomes in a phylogeny-based likelihood context. This also allows us to evaluate the change in skew continuously along the mitochondrial genome and to directly estimate relative substitution rates. Our results indicate that different types of mutation respond differently to the D(ssH) gradient. A primary role for hydrolytic deamination of cytosines in creating variance in skew among genes was not supported, but rather linearly increasing rates of mutation from adenine to hypoxanthine with D(ssH) appear to drive regional differences in skew. Substitutions due to hydrolytic deamination of cytosines, although common, appear to quickly saturate, possibly due to stabilization by the mitochondrial DNA single-strand-binding protein. These results should form the basis of more realistic models of DNA and protein evolution in mitochondria.     

Is GC bias in the nuclear genome of the carnivorous plant Utricularia?driven by ROS-based mutation and biased gene conversion?

At less than 90 Mbp, the tiny nuclear genome of the carnivorous bladderwort plant Utricularia is an attractive model system for studying molecular evolutionary processes leading to genome miniaturization. Recently, we reported that expression of genes encoding DNA repair and reactive oxygen species (ROS) detoxification enzymes is highest in Utricularia traps, and we argued that ROS mutagenic action correlates with the high nucleotide substitution rates observed in the Utricularia plastid, mitochondrial, and nuclear genomes. Here, we extend our analysis of 100 nuclear genes from Utricularia and related asterid eudicots to examine nucleotide substitution biases and their potential correlation with ROS-induced DNA lesions. We discovered an unusual bias toward GC nucleotides, most prominently in transition substitutions at the third position of codons, which are presumably silent with respect to adaptation. Given the general tendency of biased gene conversion to drive GC bias, and of ROS to induce double strand breaks requiring recombinational repair, we propose that some of the unusual features of the bladderwort and its genome may be more reflective of these nonadaptive processes than of natural selection.     

Evidence for a slowed rate of molecular evolution in the order acipenseriformes

A test of the hypothesis that the members of the order Acipenseriformes (sturgeons and paddlefishes) possess a slowed rate of molecular evolution was carried out by conducting relative-rate comparisons with representatives of four groups of teleost fishes (Cypriniformes, Elopomorpha, Salmonidae, and Percomorpha) using 21 nuclear or mitochondrial protein loci and the nuclear and mitochondrial small subunit rRNA genes, obtained from the literature or our own research. In 70 out of 81 comparisons between individual taxa (86%), acipenseriform sequences showed slower rates of change than the homologous teleost loci examined. When teleost sequences are considered together, 21 of the 23 loci show slower rates of substitution in the acipenseriform lineage. Teleost proteins show 1.85 times as many unique amino acid differences as acipenseriform proteins, when both are compared with outlier sequences. These results support a hypothesis of slowed molecular evolutionary rate in the Acipenseriformes.     

Housekeeping genes for phylogenetic analysis of eutherian relationships

The molecular relationship of placental mammals has attracted great interest in recent years. However, 2 crucial and conflicting hypotheses remain, one with respect to the position of the root of the eutherian tree and the other the relationship between the orders Rodentia, Lagomorpha (rabbits, hares), and Primates. Although most mitochondrial (mt) analyses have suggested that rodents have a basal position in the eutherian tree, some nuclear data in combination with mt-rRNA genes have placed the root on the so-called African clade or on a branch that includes this clade and the Xenarthra (e.g., anteater and armadillo). In order to generate a new and independent set of molecular data for phylogenetic analysis, we have established cDNA sequences from different tissues of various mammalian species. With this in mind, we have identified and sequenced 8 housekeeping genes with moderately fast rate of evolution from 22 placental mammals, representing 11 orders. In order to determine the root of the eutherian tree, the same genes were also sequenced for 3 marsupial species, which were used as outgroup. Inconsistent with the analyses of nuclear + mt-rRNA gene data, the current data set did not favor a basal position of the African clade or Xenarthra in the eutherian tree. Similarly, by joining rodents and lagomorphs on the same basal branch (Glires hypothesis), the data set is also inconsistent with the tree commonly favored in mtDNA analyses. The analyses of the currently established sequences have helped examination of problematic parts in the eutherian tree at the same time as they caution against suggestions that have claimed that basal eutherian relationships have been conclusively settled.     

Rates of nuclear and cytoplasmic mitochondrial DNA sequence divergence in mammals

Differential rates of nucleotide substitution among different gene segments and between distinct evolutionary lineages is well documented among mitochondrial genes and is likely a consequence of locus-specific selective constraints that delimit mutational divergence over evolutionary time. We compared sequence variation of 18 homologous loci (15 coding genes and 3 parts of the control region) among 10 mammalian mitochondrial DNA genomes which allowed us to describe different mitochondrial evolutionary patterns and to produce an estimation of the relative order of gene divergence. The relative rates of divergence of mitochondrial DNA genes in the family Felidae were estimated by comparing their divergence from homologous counterpart genes included in nuclear mitochondrial DNA (Numt, pronounced "new might"), a genomic fossil that represents an ancient transfer of 7.9 kb of mitochondrial DNA to the nuclear genome of an ancestral species of the domestic cat (Felis catus). Phylogenetic analyses of mitochondrial (mtDNA) sequences with multiple outgroup species were conducted to date the ancestral node common to the Numt and the cytoplasmic (Cymt) mtDNA genes and to calibrate the rate of sequence divergence of mitochondrial genes relative to nuclear homologous counterparts. By setting the fastest substitution rate as strictly mutational, an empirical "selective retardation index" is computed to quantify the sum of all constraints, selective and otherwise, that limit sequence divergence of mitochondrial gene sequences over time.      

Evolutionary rates of commonly used nuclear and organelle markers of Arabidopsis relatives (Brassicaceae)

Recovering the genetic divergence between species is one of the major interests in the evolutionary biology. It requires accurate estimation of the neutral substitution rates. Arabidopsis thaliana, the first whole-genome sequenced plant, and its out-crossing relatives provide an ideal model for examining the split between sister species. In the study, rates of molecular evolution at markers frequently used for systematics and population genetics, including 14 nuclear genes spanning most chromosomes, three noncoding regions of chloroplast genome, and one intron of mitochondrial genome, between A. thaliana and four relatives were estimated. No deviation from neutrality was detected in the genes examined. Based on the known divergence between A. thaliana and its sisters about 8.0-17.6 MYA, evolutionary rates of the eighteen genes were estimated. Accordingly, the ratio of rates of synonymous substitutions among mitochondrial, chloroplast and nuclear genes was calculated with an average and 95% confidence interval of 1 (0.25-1.75): 15.77 (7.48-114.09): 74.79 (36.27-534.61). Molecular evolutionary rates of nuclear genes varied, with a range of 0.383-0.856×10(-8) for synonymous substitutions per site per year and 0.036-0.081×10(-9) for nonsynonymous substitutions per site per year. Compared with orthologs in Populus, a long life-span tree, genes in Arabidopsis evolved faster in an order of magnitude at the gene level, agreeing with a generation time hypothesis. The estimated substitution rates of these genes can be used as a reference for molecular dating.     

Evolutionary rates of mitochondrial genomes correspond to diversification rates and to contemporary species richness in birds and reptiles

Rates of biological diversification should ultimately correspond to rates of genome evolution. Recent studies have compared diversification rates with phylogenetic branch lengths, but incomplete phylogenies hamper such analyses for many taxa. Herein, we use pairwise comparisons of confamilial sauropsid (bird and reptile) mitochondrial DNA (mtDNA) genome sequences to estimate substitution rates. These molecular evolutionary rates are considered in light of the age and species richness of each taxonomic family, using a random-walk speciation–extinction process to estimate rates of diversification. We find the molecular clock ticks at disparate rates in different families and at different genes. For example, evolutionary rates are relatively fast in snakes and lizards, intermediate in crocodilians and slow in turtles and birds. There was also rate variation across genes, where non-synonymous substitution rates were fastest at ATP8 and slowest at CO3. Family-by-gene interactions were significant, indicating that local clocks vary substantially among sauropsids. Most importantly, we find evidence that mitochondrial genome evolutionary rates are positively correlated with speciation rates and with contemporary species richness. Nuclear sequences are poorly represented among reptiles, but the correlation between rates of molecular evolution and species diversification also extends to 18 avian nuclear genes we tested. Thus, the nuclear data buttress our mtDNA findings.     

Spatial covariation of mutation and nonsynonymous substitution rates in vertebrate mitochondrial genomes

Mitochondrial genomes encode fundamental subunits of the basic energy producing machinery of eukaryotic cells that are under strong functional constraint. Paradoxically, these genes evolve rapidly in general, and there is substantial variation in evolutionary rates among genes within genomes. In order to investigate spatial variation in selection intensity, we conducted tests of neutrality using ratios of synonymous to nonsynonymous substitutions (dN/dS = omega) on numerous protein gene segments from fishes and mammals. Values of omega were very low for nearly all genomic regions. However, values of both omega and dN varied in a clinal pattern with increasing distance from the light-strand origin of replication. Spatial heterogeneity of nonsynonymous substitution rates exhibits a significantly positive correlation with variation in mutation rates that are related to the mode of mitochondrial DNA replication. The finding that nonsynonymous substitution rates are proportional to mutation rates is expected if a majority of substitutions are selectively neutral or slightly deleterious. Spatial patterns of among-gene variation in nonsynonymous rates were highly similar between fishes and mammals, suggesting that forces governing mitochondrial gene evolution have remained relatively constant over 450 Myr of vertebrate evolution. Conservation of substitution patterns despite major shifts in thermal habit and metabolic demands among taxa implicates a conserved replication mechanism controlling relative mutation rates as a major determinant of mitochondrial protein evolution.     

Different rates of substitution may produce different phylogenies of the eutherian mammals

Summary In an attempt to resolve some points of branching order in the phylogeny of the eutherian mammals, a phylogenetic analysis of 26 nuclear and 6 mitochondrial genes was undertaken using a maximum likelihood method on a constant rate stochastic model of molecular evolution. Seventeen of the nuclear genes gave a primates/artiodactyls grouping highest support whereas three of the mitochondrial genes found a rodents/artiodactyls grouping to be best supported. The primates/rodents grouping was never the best supported. On the assumption that rodents are indeed an outgroup to primates and artiodactyls and that the latter taxa diverged 70 million years ago, an estimation was made, for each gene, of the time of divergence of the rodent lineage. In most cases such estimates were beyond the limits set by present interpretations of the paleontological record as were many estimates of the divergence time of mouse and rat. These results suggest that, although there is locus variation, the divergent position of the rodent lineage may be an artifact of an elevated rate of nucleotide substitution in this order.