Rodent Phylogeny and a Timescale for the Evolution of Glires: Evidence from an Extensive Taxon Sampling Using Three Nuclear Genes

Dorothée Huchon*,1, Ole Madsen{dagger}, Mark J. J. B. Sibbald*, Kai Ament{dagger}, Michael J. Stanhope{ddagger},2, François Catzeflis*, Wilfried W. de Jong{dagger} and Emmanuel J. P. Douzery*

*Laboratoire de Paléontologie, Paléobiologie et Phylogénie-CC064, Institut des Sciences de l'Evolution UMR 5554/CNRS, Université Montpellier II, Place E. Bataillon, Montpellier Cedex 05, France;
{dagger}Department of Biochemistry, University of Nijmegen, The Netherlands;
{ddagger}Biology and Biochemistry, Queen's University of Belfast, UK


    Abstract
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 
Rodentia is the largest order of placental mammals, with approximately 2,050 species divided into 28 families. It is also one of the most controversial with respect to its monophyly, relationships between families, and divergence dates. Here, we have analyzed and compared the performance of three nuclear genes (von Willebrand Factor, interphotoreceptor retinoid-binding protein, and Alpha 2B adrenergic receptor) for a large taxonomic sampling, covering the whole rodent and placental diversity. The phylogenetic results significantly support rodent monophyly, the association of Rodentia with Lagomorpha (the Glires clade), and a Glires + Euarchonta (Primates, Dermoptera, and Scandentia) clade. The resolution of relationships among rodents is also greatly improved. The currently recognized families are divided here into seven well-defined clades (Anomaluromorpha, Castoridae, Ctenohystrica, Geomyoidea, Gliridae, Myodonta, and Sciuroidea) that can be grouped into three major clades: Ctenohystrica, Gliridae + Sciuroidea, and a mouse-related clade (Anomaluromorpha, Castoridae + Geomyoidea, and Myodonta). Molecular datings based on these three genes suggest that the rodent radiation took place at the transition between Paleocene and Eocene. The divergence between rodents and lagomorphs is placed just at the K-T boundary and the first splits among placentals in the Late Cretaceous. Our results thus tend to reconcile molecular and morphological-paleontological insights.


    Introduction
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 
Rodents form the most abundant and diversified order of living mammals, representing about 40% of the total number of mammalian species. Their tremendous diversity has always been a challenge for those interested in their origins, ways of radiation, and times of diversification. A major debate was raised by the claim that "the guinea pig is not a rodent" (Graur, Hide, and Li 1991Citation ; D'Erchia et al. 1996Citation ). This proposal not only contradicted the conventional view of a monophyletic origin of the rodents but also conflicted with the familiar association of rodents and lagomorphs into a clade called Glires. Like Rodentia, the Glires are well supported by morphological synapomorphies (Luckett and Hartenberger 1993Citation ). As a consequence, these questions relating to rodent relationships are currently disputed between morphology and some molecular approaches. On the one hand, molecular results—mainly based on complete mitochondrial genomes—have been criticized because only a few rodent and lagomorph species were included in the analyses, i.e., mouse, rat, guinea pig (D'Erchia et al. 1996Citation ), more recently complemented with dormouse (Reyes, Pesole, and Saccone 1998Citation ), squirrel (Reyes et al. 2000Citation ), cane rat (Mouchaty et al. 2001Citation ), and rabbit (Gissi, Gullberg, and Arnason 1998Citation ). It has been suggested that the proposed paraphyly of rodents is an artifact that might be corrected by the analysis of a larger taxonomic sampling (Luckett and Hartenberger 1993Citation ; Philippe 1997Citation ) and the use of appropriate models of sequence evolution (Sullivan and Swofford 1997Citation ). On the other hand, morphological synapomorphies for rodents are mainly based on dental and cranial characters (Luckett and Hartenberger 1993Citation ) that could be the result of ecological constraints and homoplasies (e.g., Li et al. 1992Citation ). Recently, the growing number of molecular markers—in combination with a broader species sampling within the order Rodentia—indeed provided weak to moderate support for rodent monophyly, both with mitochondrial (Nedbal, Honeycutt, and Schlitter 1996Citation ) and nuclear sequences (Huchon, Catzeflis, and Douzery 1999Citation , 2000Citation ; Adkins et al. 2001Citation ; DeBry and Sagel 2001Citation ). The combination of numerous independent nuclear markers now even led to a robust support for rodent monophyly (Murphy et al. 2001aCitation ).

The relationships between rodent families are also disputed. Morphological approaches have been frustrated by convergent evolution of characters (e.g., Jaeger 1988Citation ), and the current intraorder classification is therefore largely unsatisfactory. For example, the long-standing division of rodents into Hystricomorpha, Myomorpha, and Sciuromorpha, on the basis of the insertion patterns of masseter muscles, or alternatively into Hystricognathi and Sciurognathi, on the basis of the plane of incisor insertions, have both been shown to be inadequate (Hartenberger 1985Citation ; Nedbal, Honeycutt, and Schlitter 1996Citation ; Huchon, Catzeflis, and Douzery 2000Citation ; Adkins et al. 2001Citation ). However, the monophyly of most rodent families seems well established (Hartenberger 1985Citation ). Few molecular studies have until now investigated the relationships among rodent families, and none have included all rodent lineages or the same taxonomic sampling (or both), thus limiting the comparison of the phylogenetic results (Catzeflis et al. 1995Citation ; Nedbal, Honeycutt, and Schlitter 1996Citation ; Matthee and Robinson 1997Citation ; Robinson et al. 1997Citation ; Huchon, Catzeflis, and Douzery 1999Citation , 2000Citation ; Robinson-Rechavi, Ponger, and Mouchiroud 2000Citation ; Adkins et al. 2001Citation ; DeBry and Sagel 2001Citation ; Murphy et al. 2001a;Citation Montgelard et al. 2002bCitation ). For example, singular families like Ctenodactylidae or Anomaluridae are rarely sampled. In the present study we include for the first time all major rodent lineages—namely all sciurognath families and hystricognath superfamilies—the two lagomorph families, and representatives of all other eutherian orders, in order to investigate rodent monophyly as well as the phylogenetic relationships between rodent families.

The timing of the origin of rodents is also controversial. Fossil evidence indicates a radiation of the rodents 55 MYA (Hartenberger 1998Citation ), whereas molecular clocks based on a few rodent species tend to support a Cretaceous origin and diversification of rodents, 89–125 MYA for the divergence of murids or hystricognaths (Janke, Xu, and Arnason 1997Citation ; Kumar and Hedges 1998Citation ; Cao et al. 2000Citation ). Recent analyses, including a broader rodent sampling, but based on single sequences, suggest that the radiation of rodent families is older than the K-T limit (75 Myr) (Huchon, Catzeflis, and Douzery 2000Citation ; Adkins et al. 2001Citation ). One should note that all molecular datings for rodents are based on different genes, different sampling, and different methods. We compare here the dating performance of three genes using quartet dating, a method to estimate molecular ages of divergence that allows for rate heterogeneity between lineages (Rambaut and Bromham 1998Citation ).

We analyze three nuclear genes that are not genetically linked and code for proteins that do not display any known biological interactions: the alpha 2B adrenergic receptor (A2AB), the exon 1 of interphotoreceptor retinoid-binding protein (IRBP), and the exon 28 of von Willebrand factor (vWF). The choice of these nuclear markers was based on three considerations. First, all three sequences have successfully been used to reconstruct the phylogeny of eutherians at various taxonomic levels (Porter, Goodman, and Stanhope 1996Citation ; Springer et al. 1997a,Citation 1997b;Citation DeBry and Sagel 2001Citation ; Delsuc et al. 2001Citation ; Madsen et al. 2001Citation ). Second, they display similar sizes and numbers of variable sites, which favors the comparison of their phylogenetic performance. Third, nuclear genes have been suggested to perform better than mitochondrial ones (Springer et al. 2001Citation ). The separate and combined phylogenetic analyses of A2AB, IRBP, and vWF allow to (1) define new clades among rodents, (2) strongly confirm the monophyly of Rodentia and Glires, (3) evaluate the properties of independent markers for dating purposes, and (4) suggest a Tertiary radiation of rodent families, a K-T split between rodents and lagomorphs, and a Late Cretaceous origin for placental orders.


    Materials and Methods
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 
Amplification and Sequencing
Twenty-two rodent taxa were selected to cover the current diversity of the order, with at least one representative per family or superfamily (See Supplementary Material at MBE Web site: http://www.molbiolevol.org). All 10 sciurognath lineages were sampled (Muridae, Dipodidae, Geomyidae, Heteromyidae, Gliridae, Sciuridae, Aplodontidae, Castoridae, Anomaluromorpha, Ctenodactylidae) as well as the eight lineages that represent the whole hystricognath diversity (Thryonomyidae, Petromuridae, Bathyergidae, Hystricidae, Chinchilloidea, Octodontoidea, Cavioidea, Erethizontoidea) (Huchon and Douzery 2001Citation ). The vWF exon 28 (1236 aligned positions), the 5' third of the IRBP exon 1 (1227 positions), and the A2AB gene (1170 positions) were newly determined for five (Tachyoryctes, Dipodomys, Thomomys, Castor, and Anomalurus), 21 (all but Mus), and 18 (all but Mus, Rattus, and Cavia) rodent species, respectively. Additionally, the two lagomorphs Lepus and Ochotona were sequenced for IRBP and A2AB.

DNA extractions and PCR reactions were conducted as described before (Huchon, Catzeflis, and Douzery 1999Citation ), with slight modifications: 1 M betaine was included in the PCR mixture (Henke et al. 1997Citation ), and annealing was performed at 50°C. The PCR primers for amplifying the A2AB, IRBP, and vWF sequences have been described in Stanhope et al. (1992)Citation , Springer et al. (1997bCitation ), and Huchon, Catzeflis, and Douzery (1999,Citation 2000)Citation , respectively. When the amount of amplified DNA was insufficient for direct sequencing, smaller overlapping DNA fragments were obtained by reamplification of the initial PCR product and then sequenced. Additional sequencing primers were designed when required. All vWF sequences, the A2AB sequences of Anomalurus, Aplodontia, Bathyergus, Castor, Chinchilla, Dipodomys, Echimys, Glis, Massoutiera, Petromus, Thomomys, Trichys, and the IRBP sequences of Dinomys, Lepus, Macropus, Pedetes, and Tachyoryctes were obtained manually with [{alpha}33P]ddNTP. The other IRBP sequences were determined using dye-terminator cycle-sequencing reactions and an Applied Biosystems 373A automatic sequencer. Despite several attempts, it was not possible to amplify the A2AB of Pedetes.

PCR amplifications of A2AB from Dipus, Dryomys, Erethizon, Lepus, Marmota, Tachyoryctes, and Thryonomys were as previously described (Springer et al. 1997b), except for Dryomys and Lepus where 1 M betaine and 1.3% DMSO were added to the reaction mixture. PCR products were cloned into the pGEM-T easy vector (Promega), and sequences were determined using the Big dye-terminator cycle-sequencing kit and an ABI Prism 3700 DNA Analyser (Applied Biosystems). Clones from at least two independent PCR amplifications were sequenced to detect ambiguity caused by the PCR reaction.

Accession numbers of the newly determined sequences (AJ427226–AJ427270) are given in the Supplementary Material, together with those for the other species used in this study.

Phylogenetic Analyses
A2AB, IRBP, and vWF sequences of 40 placentals and two marsupials were manually aligned. Only few gaps had to be introduced, which were coded as missing data. A glutamic acid repeat of variable length in the A2AB gene was excluded from subsequent analyses. Five nucleotide data sets were analyzed: each gene separately, vWF and IRBP combined (to assess the position of Pedetidae among rodents), and the three genes in concatenation. One protein data set comprising the concatenation of the A2AB, IRBP, and vWF amino acid sequences was also analyzed.

The models of sequence evolution used were HKY85 for nucleotides and JTT for amino acids. The HKY85 model was favored, relative to more complex models (i.e., GTR), because of computation time limitations. However, additional Bayesian analyses (data not shown) indicated that the use of either GTR or HKY85 did not have an impact on the phylogenetic conclusions. Rate heterogeneity among DNA and protein sites was described by a discrete Gamma distribution with eight categories ({Gamma}8).

Neighbor-Joining (NJ) analyses with weighted average (WAVE) maximum likelihood distances were performed with PHYLIP 3.573 (Felsenstein 1995Citation ) and WAVEBOOT 1.2 (Krajewski et al. 1999Citation ). In these analyses the nucleotide data matrix was partitioned as follows: each codon position in each of the three genes was allowed to have its own rate, base frequency, and transition-to-transversion ratio.

Maximum parsimony (MP) and maximum likelihood trees were inferred using PAUP* (Swofford 1999Citation ), version 4, releases beta 4 and 8; and TREE-PUZZLE 4.0.2 (Strimmer and von Haeseler 1996Citation ). MP nucleotide sequence analyses were conducted with equal or differential weighting of character state changes. In the latter case, the six possible nucleotide substitutions were weighted at each codon position for each of the three genes, according to their consistency index, excluding uninformative characters (Hassanin, Lecointre, and Tillier 1998Citation ).

Before running individual and combined heuristic ML PAUP* searches on nucleotide sequences, the program TREE-PUZZLE 4.0.2 was used to estimate the transition-transversion parameter ({kappa}) and the parameter ({alpha}) of the Gamma distribution of the HKY + {Gamma}8 model. For amino acid sequences, ML reconstructions under the JTT + {Gamma}8 model were obtained using the quartet-puzzling method with TREE-PUZZLE 4.0.2.

Bayesian phylogenetic analyses were performed with MrBayes 2.1 (Huelsenbeck and Ronquist 2001Citation ). The Metropolis-coupled Markov chain Monte Carlo sampling approach was used to calculate posterior probabilities. Prior probabilities for all trees were equal, starting trees were random, tree sampling was done every 20 generations, and burn-in values were determined empirically from the likelihood values. To check for consistency of results, four Markov chains were run simultaneously, twice for 200,000 and twice for 500,000 generations.

For nucleotide sequence analysis, third codon positions were excluded from ML and Bayesian analyses because of their base compositional heterogeneity between taxa (as evaluated by TREE-PUZZLE 4.0.2 and PAUP*), whereas they were kept in NJ-WAVE—where each codon position had its own rate—and MP analyses.

Robustness of the nodes of the phylogenetic trees was assessed by bootstrap (Felsenstein 1985Citation ) with PAUP* and reliability percentages (RP) with TREE-PUZZLE 4.0.2. For NJ-WAVE and MP analyses, 1,000 bootstrap replicates were computed. For ML analyses, computing time limitations forced us to estimate bootstrap percentages (BP) after only 100 replicates, with HKY + {Gamma}8 parameters set to the values estimated for the best tree, with NJ starting trees, and with 1,000 rearrangements of TBR branch swapping. RP were estimated after 10,000 puzzling steps.

Comparison of Alternative Phylogenies
The first step was to reconstruct alternative tree topologies. PAUP* heuristic searches under a single HKY + {Gamma}8 model and incorporating a topological constraint were conducted in order to identify the highest-likelihood topology that satisfied a given hypothesis (e.g., the paraphyly of rodents). Second, the alternative topologies previously identified were evaluated and compared relative to the best ML topology found for the nucleotide sequences. Statistical comparisons were conducted with partitioned maximum likelihood on the nucleotide matrix and on the protein matrix of combined data. For the nucleotide sequences, to account for differences in evolutionary substitution processes between codon positions and between the three genetically independent nuclear markers, first and second codon positions were distinguished for A2AB, IRBP, and vWF. This resulted in the definition of six partitions of sites. Six independent HKY + {Gamma}8 models were thus assumed for each of these six partitions. The topologies previously identified after PAUP* heuristic searches were evaluated and compared under the more complex 6-partitions model with PAML (Yang 1997Citation ), version 3.0d. In the latter case, all ML parameters—i.e., 6 transition-transversion rate parameters, 6 {alpha} parameters, and 6 x 81 branch parameters—were reestimated by PAML for each evaluated topology. Partitioned log-likelihoods were then compared using the Kishino and Hasegawa (1989)Citation test with the Shimodaira and Hasegawa (1999)Citation correction.

Because of computation time limitations, a similar approach could not be applied for the protein sequences. A single partition was considered, and topology comparisons were done under a JTT + {Gamma}8 model assumed for the concatenated protein data set. The 81 branch parameters were reestimated by PAML for each evaluated topology, and the Kishino and Hasegawa (1989)Citation test with the Shimodaira and Hasegawa (1999)Citation correction was performed.

Molecular Dating
Divergence times were estimated using quartet dating (Rambaut and Bromham 1998Citation ), a method allowing each lineage to have a different rate of evolution, and implemented for nucleotide sequence analysis in the program QDATE 1.1. This ML method calculates the divergence date between two pairs of calibrating lineages, each one being represented by two species for which the time of divergence is known. Four molecular datings were estimated: the first split within Rodentia ("R" in fig. 2 ), the divergence between the two Glires lineages ("G"), the divergence between Laurasiatheria and [Glires + Euarchonta] ("L/E + G"), and the first split among placental mammals ("P"). These molecular datings were derived from all combinations of, respectively (1) two time-calibrated pairs of rodents ("D1" or "D2" vs. "D3" or "D4"; fig. 2 ), (2) one pair of rodents versus one pair of lagomorphs ("D5"), (3) one pair of glires versus one pair of cetartiodactyls ("D6"), and (4) one pair of paenungulates ("D7") versus the remaining pairs of placentals. The seven pairs of calibrating taxa are detailed in table 1 .



View larger version (34K):
[in this window]
[in a new window]
 
Fig. 2.—Maximum likelihood tree (-lnL = -26,567.53) reconstructed from first and second codon positions of combined A2AB, IRBP, and vWF nucleotide sequences of 40 placentals and two marsupials. Branch lengths were computed assuming a single ML model for the three combined genes. Nucleotide substitutions were described by a HKY model with parameter {kappa} = 1.78 and rate heterogeneity among sites described by an 8-categories discrete Gamma distribution with parameter {alpha} = 0.39. The branch leading to the marsupial outgroup has been shortened three times. ML BP derived after 100 replicates are given for each node. Note that all nodes supported by more than 58% of bootstrap in ML do have Bayesian posterior probabilities ranging from 0.91 to 1.00. The three major rodent clades are indicated: Sciuroidea + Gliridae (S + G), the mouse-related clade (M), and Ctenohystrica (C). Black circles indicate the nodes that have been dated using the Quartet dating method (abbreviations: E, Euarchonta; G, Glires; L, Laurasiatheria; P, Placentalia; R, Rodentia; cf. fig. 4 ), whereas the white circles indicate the calibration points that have been used (i.e., the time-calibrated pairs D2D7; cf. table 1 )

 

View this table:
[in this window]
[in a new window]
 
Table 1 Calibration Points Used in the Quartet Dating Analyses

 
Rate constancy was tested with likelihood ratio tests (LRT) (Felsenstein 1988Citation ), and only quartets that fitted a two-rate constrained model (i.e., each dating pair of the quartet does have its own rate), relative to a free-rate model (i.e., each branch of the quartet has a different rate), were kept in the dating estimations. Because rate heterogeneity between sites can affect estimation of divergence date (Rambaut and Bromham 1998Citation ), sequence evolution was described by the HKY + {Gamma}8 model. Quartet dating results were described by the median of the distribution of the quartet date estimates, and the associated lower and upper limits of the 95% confidence interval were those of the median quartet.


    Results and Discussion
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 
Base composition homogeneity of all codon positions was evaluated at the 1% level of chi-square tests for the five data sets, A2AB, IRBP, vWF, IRBP + vWF, and A2AB + IRBP + vWF, by comparing the nucleotide composition of each sequence with the frequency distribution assumed in the ML model. These five data sets actually displayed a significant base compositional heterogeneity for 5, 4, 7, 10, and 13 species, respectively—including 2, 1, 2, 5, and 7 Glires, respectively—among the 42 mammals. A closer examination revealed that the base composition at third codon positions was responsible for this heterogeneity. For each codon position of the three genes, we computed the difference in GC content between the GC-richest sequence and the GC-poorest sequence. The difference in GC content ranges between 8.1% and 16.9% for the first codon position, 6.3% and 10.3% for the second codon position, and 24.7% and 35.9% for the third codon position, with vWF showing the most heterogeneous base composition. Some sequences presented high (e.g., Bradypus) or low (e.g., Mus, Orycteropus, Marsupialia) GC content at third codon positions for the three genes. Other sequences showed extreme GC values at third codon positions for one gene but presented average GC levels for the two other genes (e.g., Geomyidae for vWF vs. A2AB and IRBP). After exclusion of third codon positions, base composition became homogeneous for all placental taxa. Because preliminary results indicated that ML reconstructions might be affected by the base composition heterogeneity at third codon positions, all subsequent ML nucleotide analyses were conducted on first + second codon positions only. No compositional heterogeneity was detected in the protein sequences.

Phylogenies Reconstructed from Individual Genes
All tree-building methods either on nucleotide or on protein sequences (NJ-WAVE, standard and differentially weighted MP, ML, and Bayesian approach) provided the same overall phylogenies, with minor topological variations involving only weakly supported nodes. Considering ML analysis on nucleotide sequences as best representing the results of the individual genes, only those results will be detailed here. For comparative purposes, quartet puzzling reliability percentage observed on protein sequences (RPPROT) and Bayesian posterior probabilities for nucleotides (PPNUC) and proteins (PPPROT) are however indicated for the combined data set.

The trees reconstructed from first plus second codon positions of each of the three genes are given in figure 1 . All three trees suggest the monophyly of rodents and that of Sciuroidea (Marmota + Aplodontia), Geomyoidea (Dipodomys + Thomomys), Hystricognathi (Bathyergus, Thryonomys, Petromus, Trichys, Cavia, Chinchilla, Echimys, and Erethizon), and Ctenohystrica (Massoutiera + hystricognaths). Discrepancies between the trees always involve weakly supported nodes (i.e., not involving two conflicting nodes with BPML > 50), except for the position of Dipodidae and the relationships among Cetartiodactyla. According to A2AB and IRBP, Dipodidae are the sister clade of Muridae (BPML = 67 and 74, respectively), whereas vWF clusters Dipodidae with Geomyoidea (BPML = 65). The grouping of Dipodidae with Geomyidae may be an artifact resulting from similar base compositions and rapid rates of evolution of their vWF sequences (data not shown). Within Cetartiodactyla, A2AB and IRBP place Lama in the most basal position relative to other cetartiodactyls (BPML = 51 and 81, respectively), whereas vWF clusters Lama with Sus (BPML = 66).



View larger version (31K):
[in this window]
[in a new window]
 
Fig. 1.—Maximum likelihood trees and corresponding BP computed from first and second codon positions of each nuclear marker. The length of the branch leading to marsupials has been reduced three times. The following ML parameters of the HKY + {Gamma}8 model maximize the log-likelihood for A2AB (-lnL = 6,858.67): transition-transversion parameter {kappa} = 1.61, {alpha} = 0.33, for IRBP (-lnL = 8,824.71): {kappa} = 2.11, {alpha} = 0.41, and for vWF (-lnL = 10,399.48): {kappa} = 1.58, {alpha} = 0.52. Nodes labeled with asterisks do not occur in the bootstrap consensus tree

 
In spite of having similar lengths, the three nuclear genes do not contain the same phylogenetic signal. Each gene strongly supports some nodes, but these nodes might be different from one gene to another, and none is able to solve the whole rodent phylogeny. For example, only A2AB provides high support for rodent monophyly (BPML = 83), IRBP for Myodonta monophyly (i.e., Muridae + Dipodidae, BPML = 74), and vWF for Ctenohystrica monophyly (BPML = 91). The results also indicate that the resolving power of each gene is not restricted to a given taxonomic level. For example, A2AB, unlike vWF, is able to solve deep relationships like rodent monophyly as well as more recent relationships like Sciuroidea (BPML = 98) but not intermediate clades like Ctenohystrica (BPML = 41), whereas vWF does so.

Phylogenies Reconstructed from the Concatenated Genes
When codon positions 1 and 2 of the combined A2AB + IRBP + vWF genes were analyzed under a single HKY + {Gamma}8 model, the log-likelihood of the best topology was lnL = -26,415.56. Assuming six independent HKY + {Gamma}8 models for each of the six partitions yielded lnL = -25,988.88 and resulted in a significant increase of log-likelihood (LRT statistics = 853.36; df = 430; P < 0.0001). The ML parameters estimated for the six partitions are given in table 2 . All three genes exhibit similar base compositions, on first as well as on second codon positions. The slowest evolving partition is the second codon position of A2AB, followed by—with increasing relative rate—A2AB (first codon position), IRBP (second), vWF (second), IRBP (first), and vWF (first). Second codon positions are more heterogeneous than first positions in terms of substitution rates for the three markers, but they display a higher transition-transversion rate parameter for IRBP and vWF.


View this table:
[in this window]
[in a new window]
 
Table 2 General Characteristics and Maximum Likelihood Estimates of Nucleotide Substitution Parameters for the Six Partitions Used in Phylogenetic Analyses: First and Second Codon Positions of the Three Nuclear Markers (A2AB, IRBP, and vWF)

 
Combination of the Three Markers
Although the A2AB, IRBP, and vWF trees do not show any major topological incongruences, crossed Shimodaira-Hasegawa tests indicate that each nucleotide data set rejects the highest-likelihood topology of the two other data sets (table 3 ). However, none of the three genes rejects the ML topology obtained from the combined data set. Consequently, this a posteriori observation suggests that the three genes can be combined and that the combined data tree accordingly appears to be the "best provisional phylogenetic hypothesis" (Adkins et al. 2001Citation ). The combination of A2AB, IRBP, and vWF leads to a topology that stabilizes the phylogenetic position of rodents among mammals and contributes to resolve most of the relationships between rodent families.


View this table:
[in this window]
[in a new window]
 
Table 3 Maximum Likelihood Evaluation of the Congruence Between the Three Nuclear Markers

 
Position of Rodentia Among Mammals
Phylogenetic analyses based on first plus second codon positions of the three concatenated nuclear genes all indicate the monophyly of rodents, its support being the highest under the weighted MP (not shown) and the maximum likelihood approaches (BPML = 95: fig. 2 ; RPPROT = 74). Maximum likelihood tests of various phylogenetic hypotheses indicate that the alternative to the monophyly of rodents is significantly less likely (table 4 ). The Bayesian approach also provides a posterior probability of 1.00 for the monophyly of rodents for both nucleotide and protein sequences (trees not shown). Our results thus confirm statistically the monophyly of Rodentia with an extended sampling of this order—including all sciurognath families and hystricognath superfamilies—and with representatives of all other placental orders. In fact, with the recent increase of available complete mitochondrial genomes, rodent monophyly is no longer statistically rejected either (Cao et al. 2000Citation ; Mouchaty et al. 2001Citation ). It thus appears that the paraphyly of rodents will be difficult to defend with the broader taxonomic sampling within the order and with the growing number of molecular markers supporting their monophyly.


View this table:
[in this window]
[in a new window]
 
Table 4 Partitioned Maximum Likelihood Comparisons of Alternative Topologies Suggested by Morphology, Paleontology, and Molecular Data Relative to the Best Tree and Conducted Under the Combination of the Nucleotide and Protein Matrices of the Three Markers

 
In agreement with the independent nuclear markers analyzed by Murphy et al. (2001a,Citation 2001b)Citation , the concatenation of A2AB, IRBP, and vWF suggests that Rodentia is the sister group of Lagomorpha, constituting the superorder Glires (PPNUC and PPPROT = 1.00; BPML = 84; RPPROT = 37). The best alternative hypothesis to the Glires monophyly appears to be significantly less likely with protein sequences but only marginally significant with nucleotide sequences (table 4 ). We also verify that the Euarchonta (Primates + Dermoptera, and Scandentia) are the sister clade of Glires (PPNUC and PPPROT = 1.00; BPML = 66: fig. 2 ; RPPROT < 20), and the best alternative hypothesis appears to be significantly less likely with the nucleotide sequences only (table 4 ). This superordinal clade "Euarchontoglires" has recently been proposed by molecular studies (Madsen et al. 2001Citation ; Murphy et al. 2001a,Citation 2001bCitation ) but has never been suggested by morphology. Morphological studies generally cluster Glires with Macroscelidae, to form the Anagalida, and Primates, Dermoptera, and Scandentia with the Chiroptera, to form the Archonta (e.g., McKenna and Bell 1997, p. 295Citation ). However, Anagalida and Archonta are based on only a few synapomorphies, and these clades have been rejected by molecular studies (e.g., Murphy et al. 2001a;Citation Springer et al. 1997a,Citation 1999Citation ). It is interesting to note that a recent paleontological study clusters Macroscelidae with Proboscidea (Tabuce et al. 2001Citation ) and that a relationship between Glires and Primates has been suggested (McKenna 1986Citation ). Our results are thus not at odds with morphology and suggest that paleontological evidence should search for a sister clade relationship between Glires and Euarchonta. Relationships among mammals as shown in figure 2 are not further discussed because of insufficient species sampling.

Relationships Among Rodents
Our results suggest the division of Rodentia into three major infraordinal clades. The first one comprises squirrel- and dormouse-related animals: Sciuroidea (Sciuridae [squirrels] + Aplodontidae [mountain beavers]) and Gliridae (dormice). The second clade contains mouse-related rodents: Myodonta (Muridae [mouse, rats] + Dipodidae [jerboas]), Castoridae (beavers), Geomyoidea (Geomyidae [pocket gophers] + Heteromyidae [pocket mice]), and Anomaluromorpha (Anomaluridae [scaly-tailed flying squirrels] + Pedetidae [springhares]). The third clade (Ctenohystrica sensu Huchon, Catzeflis, and Douzery 2000Citation ) contains gundi and Guinea pig–related rodents: Ctenodactylidae and Hystricognathi. The interrelationships between these three clades are poorly resolved. None of the three bifurcating topologies connecting them involves significantly different log-likelihood (0.28 < PSH < 0.32), and the results are sensitive to the method of reconstruction and the data set considered. Ctenohystrica clusters with the mouse-related clade for ML and Bayesian nucleotide analyses (BPML = 40: fig. 2 ; PPNUC = 0.83). However, Bayesian protein analysis clusters Ctenohystrica with the squirrel- and dormouse-related clade (PPPROT = 0.70), and quartet puzzling protein analysis clusters the mouse-related clade with the squirrel- and dormouse-related clade (RPPROT = 42).

The Squirrel- and Dormouse-Related Clade
The monophyly of Sciuroidea (Sciuridae and Aplodontidae) has been supported by morphological (e.g., Lavocat and Parent 1985Citation ; Meng 1990Citation ) and molecular data (Huchon, Catzeflis, and Douzery 1999Citation , 2000Citation ; Adkins et al. 2001Citation ; DeBry and Sagel 2001Citation ). The increase in taxonomic sampling did not reduce the support for this clade (BPML = 100; RPPROT = 85). In agreement with DeBry and Sagel (2001)Citation , we observe that Sciuroidea are characterized by a unique insertion of three amino acids in the IRBP gene. The sister clade of Sciuroidea appears to be the Gliridae (fig. 2 ). The existence of this suprafamilial clade has been suggested by morphological (e.g., Lavocat and Parent 1985Citation ; Meng 1990Citation ) and mitochondrial (Reyes et al. 2000Citation ) studies; however, similar to other molecular studies (e.g., Huchon, Catzeflis, and Douzery 1999Citation ; Adkins et al. 2001Citation ; DeBry and Sagel 2001Citation ; Montgelard et al. 2002bCitation ), the support here remains moderate to strong. It is noteworthy that the support is higher with protein sequences (PPPROT = 1.00; RPPROT = 75; PSH < 0.05) than with nucleotide sequences (PPNUC = 1.00; BPML = 58; PSH < 0.15).

The Mouse-Related Clade
The grouping of Anomaluromorpha, Castoridae, Geomyoidea, and Myodonta (fig. 2 ) has never been suggested by morphological and paleontological observations. Castoridae has usually been related to Sciuridae because both families share the sciuromorph and sciurognath states (Brandt 1855Citation ; Tullberg 1899Citation ). However, some morphological studies could not confirm this relationship (Bugge 1985Citation ; Lavocat and Parent 1985Citation ; Meng 1990Citation ), and it has even been suggested that Castoridae might be more closely related to Muridae than to Sciuridae (Meng 1990Citation ). Anomaluridae and Pedetidae were considered as enigmatic families, possibly related to Ctenodactylidae or Hystricognathi (e.g., Luckett and Hartenberger 1985Citation ; Jaeger 1988Citation ). Geomyoidea has, less ambiguously, been associated with the Muridae (e.g., Wahlert 1985Citation ; Ryan 1989Citation ).

The mouse-related clade is moderately to strongly supported (PPNUC and PPPROT = 1.00; BPML = 65: fig. 2 ; RPPROT = 43), and alternative hypotheses cannot be significantly rejected (table 4 ). Such a molecular support and lack of morphological evidence might reflect that we are dealing with a phylogenetic artifact, but independent molecular studies identified the same node, although it was not explicitly noticed and discussed (Adkins et al. 2001Citation ; Murphy et al. 2001aCitation ). Consequently, the naturalness of this new superfamilial arrangement must be seriously considered in future studies. Another interesting aspect of this mouse-related clade is that it includes animals displaying different jaw patterns and having different geographical origins. Myodonta has been suggested to have originated from Asian hystricomorph rodents, Anomaluromorpha might have originated from African hystricomorph rodents, and Geomyidae and Castoridae are sciuromorph rodents from North America (e.g., Vianey-Liaud 1985Citation ; Hartenberger 1998Citation ). This suggests a complicated biogeographical history and an ancient origin for this rodent clade. This group might have a Paleocene or an Early-Eocene origin because Muridae, Dipodidae, Anomaluridae, Geomyoidea, and Castoridae are rooted in—or related to—Early- and Middle-Eocene families: Cricetidae, Zapodidae, Zegdoumyidae, Eomyidae, and Eutypomyidae, respectively (Hartenberger 1998Citation ).

The mouse-related clade is divided into three subclades of morphologically distinct rodents, of which the relationships are unclear. Two subclades have been recognized earlier: (1) Myodonta, i.e., murids and jerboas (Luckett and Hartenberger 1985Citation ; Nedbal, Honeycutt, and Schlitter 1996Citation ; Huchon, Catzeflis, and Douzery 1999Citation ; DeBry and Sagel 2001Citation ), and (2) Anomaluromorpha, i.e., scaly-tailed flying squirrels and springhares (Bugge 1985Citation ; Lavocat and Parent 1985Citation ; Montgelard et al. 2002aCitation ). The third subclade unexpectedly clusters beavers with pocket mice and pocket gophers (Adkins et al. 2001Citation ). Here, Myodonta is highly supported (PPNUC and PPPROT = 1.00; BPML = 81, PSH < 0.03: fig. 2 , table 4 ; RPPROT = 57). With the vWF + IRBP combined data set, Anomaluromorpha is highly supported (BPML > 86, whatever the reconstruction method used; data not shown). Alternative hypotheses do not appear significantly less likely (PSH < 0.10), but the analysis is based on a shorter data set. The Geomyoidea plus Castoridae clade is moderately to strongly supported (PPNUC and PPPROT = 1.00; BPML = 63; RPPROT = 41; PSH < 0.13). This enigmatic clade has never been clearly suggested by morphological studies but had already been evidenced in molecular studies based on the GHR gene (Adkins et al. 2001Citation ), on a combined data set of 15 genes (Murphy et al. 2001aCitation ), and on 12S rRNA and cytochrome b sequences (Montgelard et al. 2002bCitation ). In contrast, DeBry and Sagel (2001)Citation found different results with the IRBP gene: Geomyoidea are the sister clade of Ctenohystrica and Castoridae clusters with Sciuroidea—an observation which might be the consequence of a smaller rodent sampling (i.e., Anomaluromorpha are missing) because our IRBP tree clusters Geomyidae and Castoridae (BPML = 57; fig. 1 ).

The lack of resolution for the branching order of the three clades Anomaluromorpha, Castoridae + Geomyidae, and Myodonta is illustrated by the ML and Bayesian analyses. Anomalurus connects either with Castor + Geomyidae (PPNUC = 0.56 and PPPROT = 0.33; BPML = 29; fig. 2 ) or with Myodonta (PPNUC = 0.40 and PPPROT = 0.36; BPML = 37).

The Ctenohystrica
The grouping of Ctenodactylidae and Hystricognathi in a Ctenohystrica clade has been strongly supported by the vWF gene (Huchon, Catzeflis, and Douzery 2000Citation ) and the GHR gene (Adkins et al. 2001Citation ). Our results confirm these observations (PPNUC and PPPROT = 1.00; BPML = 100; RPPROT = 93), and breaking the monophyly of Ctenohystrica is a significantly worse alternative (PSH < 0.03; table 4 ), as is the case for constraining Ctenodactylidae to branch with other sciurognaths (PSH < 0.03). Among Hystricognathi, the relationships obtained agree with the conclusions of Huchon and Douzery (2001)Citation . Hystricognathi is divided into three clades: Hystricidae, Phiomorpha s. s. (i.e., Bathyergidae, Thryonomyidae, and Petromuridae) and Caviomorpha (i.e., Octodontoidea, Cavioidea, Erethizontoidea, and Chinchilloidea). It is interesting to note that the increase of sequence length favors a basal position of Hystricidae (PPNUC and PPPROT = 1.00; BPML = 87; RPPROT = 74), but alternative hypotheses are not significantly less likely (PSH < 0.14; table 4 ).

Robustness of Alternative Topologies Suggested by Morphology
Following Tullberg (1899)Citation , rodents have been divided into two suborders—Sciurognathi and Hystricognathi. According to our molecular data, Hystricognathi remains a valid clade (BPML = 100; fig. 2 , RPPROT = 93), but the monophyly of Sciurognathi is statistically rejected (see earlier; table 4 ). An alternative classification divided rodents into Myomorpha, Sciuromorpha, and Hystricomorpha (Brandt 1855Citation ). The contents of these groups changed according to the authors, and we here follow the classification of McKenna and Bell (1997)Citation . Constraining the monophyly of Myomorpha, i.e., clustering Muridae + Dipodidae with Geomyoidea and Gliridae, is significantly less likely than the best ML topology with our A2AB + IRBP + vWF nucleotide and protein data (PSH < 0.05; table 4 ). The grouping of Aplodontidae, Sciuridae, and Castoridae into Sciuromorpha is marginally rejected (PSH < 0.10). Hystricomorpha actually corresponds to hystricognaths and thus appears monophyletic. These results suggest that the current subordinal classification of rodents should be thoroughly revised.

Quartet Dating of Divergence Times
Impact of the Choice of Markers for Dating Purposes
The combination of the independent paleontological calibration points, D1D7 (table 1 and fig. 2 ), allowed the construction and evaluation of 59 quartets of species for A2AB, IRBP, vWF, and their combination. Among them, only 11 quartets fitted the two-rate constrained model for the three genes. Divergence times were estimated by quartet dating for the splits between the mouse-related clade and Ctenohystrica, the rodent superfamilies, the glires, and the Laurasiatheria versus Euarchonta + Glires (fig. 3 ). The results illustrate the impact of the choice of calibration points on the date estimates (cf. Huchon, Catzeflis, and Douzery 2000Citation ). For example, A2AB estimates the divergence between Laurasiatheria and Euarchonta + Glires at 79.5 MYA with the quartet Lama-Physeter versus Marmota-Aplodontia (fig. 3 , square in quartet 11), whereas it is 111.9 Myr with the quartet Lama-Sus versus Echimys-Cavia (fig. 3 , square in quartet 8). There also is an impact of the choice of the gene on the dating results. Two different genes can lead to very different date estimates, even for the same quartet of species. For example, Mus-Rattus versus Echimys-Erethizon gives a date of 42.2 Myr with vWF (fig. 3 , circle in quartet 1), against 94.2 Myr with A2AB (fig. 3 , square in quartet 1). In this case, even the confidence intervals do not overlap: 34.3–53.2 Myr (vWF) against 66.9–137.3 Myr (A2AB). Most of the A2AB estimates appear to give older dates and larger confidence intervals (fig. 3 ), but this might just reflect the small number of quartets considered. Even when the medians of the quartet dates are compared between genes, the dating results remain divergent. For example, the median date for the divergence between Laurasiatheria and Euarchonta-Glires is 98.9 Myr for A2AB, 81.8 Myr for IRBP, and 76.9 Myr for vWF over 15, 24, and 5 quartets, respectively (data not shown). Because the three genes have similar lengths, this indicates that they contain quite different dating information. Therefore, because dating estimation becomes more accurate with longer sequences (e.g., Rambaut and Bromham 1998Citation ), we combined the three nuclear genes. As expected, the distribution of the quartet estimations clearly becomes narrower than with the single genes (fig. 3 ). Confidence intervals are also up to four times smaller when the three genes are concatenated, as for example illustrated by the comparison between A2AB versus combination of markers (fig. 3 , square vs. diamond in quartet 1).



View larger version (24K):
[in this window]
[in a new window]
 
Fig. 3.—Influence of gene and species choice on quartet dating results. The 11 independent quartets of species that are statistically compatible with a two-rates model for all three markers are represented. For each quartet, four dating estimates with their 95% confidence intervals are given: white squares for A2AB, black triangles for IRBP, white circles for vWF, and black diamonds for the combination. The following divergence dates were computed: MUR-CTE = divergence of Muridae and Ctenohystrica (based on D1 vs. D2 comparisons: see fig. 2 ); ROD = radiation of rodent superfamilies (D1 or D2 vs. D3 or D4); GLI = divergence of Rodentia and Lagomorpha (D1, D2, D3, or D4 vs. D5); LAU-EUG = divergence of Laurasiatheria and Euarchonta + Glires (D1, D2, D3, D4, or D5 vs. D6). Note that the age of the placental radiation has not been estimated here because no quartet fitted a two-rate constrained model for each of the three genes. The dotted lines indicate boundaries between geological periods (from bottom to top: Early-Late Cretaceous, K-T boundary, Paleocene-Eocene, and Eocene-Oligocene). The detail of species content and divergence date of each individual quartet is available upon request

 
Molecular Dating Using Combined Markers
The timing of the diversification of extant placental orders is heavily debated. According to paleontology, they diverged in the Paleocene, 65–55 MYA (e.g., Alroy 1999Citation ). Molecules rather suggest that placentals were already diversified in the Cretaceous (review in Bromham, Phillips, and Penny 1999Citation ). However, various molecular studies give deviating datings. Their estimations are generally based on few calibration points and use either distance (e.g., Hedges et al. 1996Citation ; Janke, Xu, and Arnason 1997Citation ; Springer et al. 1997b;Citation Kumar and Hedges 1998Citation ; Waddell et al. 1999Citation ; review in Bromham, Phillips, and Penny 1999Citation ), or Bayesian (Cao et al. 2000Citation ), or ML (Eizirik, Murphy, and O'Brien 2001Citation ) approaches.

Among the 59 quartets of species involving calibration points D1D7, only 41 fitted the two-rates model for the concatenation of the three markers. Among them, respectively, 1, 3, 6, 7, 17, and 7 allowed to estimate the age of the split between Sciuroidea and Gliridae (not shown), Muridae and Ctenohystrica (not shown), rodent superfamilies, glires, Laurasiatheria and Euarchonta + Glires, and placentals (fig. 4 ). Our results, being based on ML estimates and several independent calibrating pairs of taxa, are congruent with paleontology about ages of rodent radiation and Glires divergence. They support a rodent diversification at the transition between Paleocene and Eocene, 55.8 MYA (median of confidence intervals: 49.4–63.7; fig. 4 ). Paleontological studies place the rodent radiation at 55 MYA at the Paleocene-Eocene boundary (Hartenberger 1998Citation ). The split between rodents and lagomorphs is estimated to be 64.5 MYA (57.3–73.3), just at the Cretaceous-Tertiary boundary. These results indicate that rodents and lagomorphs diverged only at the beginning of the Tertiary (fig. 4 ).



View larger version (28K):
[in this window]
[in a new window]
 
Fig. 4.—Maximum likelihood estimations (black circles) and their 95% confidence intervals (bars) for the divergence times between the Rodentia superfamilies, the orders Lagomorpha and Rodentia (Glires), the Laurasiatheria and Euarchonta-Glires, and the placental orders. Quartet dating results are presented in descending order, and the medians of the quartet distributions are represented by empty circles and bars, indicating the median of the 95% confidence intervals: 55.8 Myr (49.4–63.7); 64.5 Myr (57.3–73.3); 83.2 Myr (74.1–94.4); and 101.2 (88.5–116.4). The Cretaceous-Tertiary boundary (K-T) is indicated by a dashed line. The detail of species content and divergence date of each individual quartet is available upon request

 
Molecular dating indicates a Late Cretaceous divergence, 83.2 MYA (74.1–94.4; fig. 4 ), for the last common ancestor of Glires and Laurasiatheria (represented here by cetartiodactyls). This date is congruent with paleontological studies: ungulates have been linked to the 85 MYA zhelestids (Archibald 1996Citation ), and Glires have been strongly related to the zalambdalestids whose oldest fossil is estimated to be 90 MYA (Archibald, Averianov, and Ekdale 2001Citation ). This estimate also agrees with the quartet dating (63.5–95.3 Myr) of Eizirik, Murphy, and O'Brien (2001)Citation using independent calibration points. First divergences within placentals seem to occur 101.2 MYA (88.5–116.4; fig. 4 ), at the transition between Early and Late Cretaceous, but there are huge variations between ages provided by different quartets. Additional calibration points within Laurasiatheria are likely to stabilize these estimates (Eizirik, Murphy, and O'Brien 2001Citation ). These molecular dates should also be confirmed in the future using other dating methods, such as for example the Bayesian approach of Thorne, Kishino, and Painter (1998)Citation .

Our dating results on Glires agree with the hypothesis that Mesozoic placental divergences only involved the first basal clades of the placental tree, whereas the diversification of extant lineages occurred after the K-T boundary (Hedges et al. 1996Citation ; Alroy 1999Citation ; Eizirik, Murphy, and O'Brien 2001Citation ; Madsen et al. 2001Citation ). Dating involving other placental clades are required to verify whether the pattern observed for rodents (i.e., a diversification only after the K-T boundary) is valid for other placental orders too.


    Acknowledgements
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 
This work would not have been possible without the essential contribution of all collectors who provided mammalian tissues now housed in the collection of Montpellier: Marina Baskevich (for Dryomys), Dona L. Dittmann (Thomomys; Collection of Genetic Resources, Louisiana State University Museum of Natural Sciences), Jean-Marc Duplantier (Lepus), Chris G. Faulkes (Bathyergus), Piotr Gambarian (Dipus), Laurent Granjon (Thryonomys, Lepus), Jean-Claude Gautun (Anomalurus), Patrick Gouat (Massoutiera), Jack Hayes and Marilyn Banta (Dipodomys), Robert S. Hoffmann (Marmota), Reginald Hoyt and John Trupkiewicz (Petromus; Zoological Society of Philadelphia), John A. W. Kirsh (Erethizon), Vincent Laudet (Chinchilla), Conrad Matthee (Pedetes), David Nolta (Aplodontia), E. Pelé and Vitaly Volobouev (Tachyoryctes), Thierry Petit (Macropus; Zoo de La Palmyre, France), Philippe Perret (Cavia), Rob Stuebing (Trichys), and the Zoo du Lunaret in Montpellier (Castor).

We wish to thank Ron DeBry for sharing the Dipodomys spectabilis vWF sequence, Christophe Douady for providing the Ochotona IRBP sequence, Jean-Louis Hartenberger, Jean-Jacques Jaeger, and Laurent Marivaux for useful paleontological discussions and comments, Jacques Demaille and Denis Pugnère (Institut de Génétique Humaine de Montpellier, CNRS UPR 1142) for computing facilities, and Emma Teeling for bench advice with regard to the automatic sequencer.

This work was supported by ACC-SV7 (Réseau National de Biosystématique), ACC-SV3 (Réseau coordonné par D. Mouchiroud), European Community TMR Network "Mammalian phylogeny" FMRX-CT98-0221 (to M.S., F.C., and W. de J.), and the Genopole Région Montpellier Languedoc-Roussillon and the Action inter-EPST Bioinformatique 2000–2002 (to E.D.). D.H. acknowledges the financial support of a Lavoisier grant from the French Ministry of Foreign Affairs. This is contribution number 2002-008 of the Institut des Sciences de l'Evolution de Montpellier (UMR 5554-CNRS).


    Footnotes
 
Dan Graur, Reviewing Editor

1 Present address: Molecular Evolution Laboratory, Faculty of Bioscience and Biotechnology, Tokyo Institute of Technology, Nagatsuta-cho, Midori-ku, Yokohama, Japan Back

2 Present address: Bioinformatics, GlaxoSmithKline, UP1345, Collegeville Pennsylvania Back

Keywords: Rodentia Glires Eutheria vWF IRBP A2AB phylogeny molecular dating Back

Address for correspondence and reprints: Emmanuel J. P. Douzery, Laboratoire de Paléontologie, Paléobiologie et Phylogénie-CC064, Institut des Sciences de l'Evolution UMR 5554/CNRS, Université Montpellier II, Place E. Bataillon, 34 095 Montpellier Cedex 05, France. douzery{at}isem.univ-montp2.fr Back


    References
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 

    Adkins R. M., E. L. Gelke, D. Rowe, R. L. Honeycutt, 2001 Molecular phylogeny and divergence time estimates for major rodent groups: evidence from multiple genes Mol. Biol. Evol 18:777-791[Abstract/Free Full Text]

    Alroy J., 1999 The fossil record of North American mammals: evidence for a Paleocene evolutionary radiation Syst. Biol 48:107-118[ISI][Medline]

    Archibald J. D., 1996 Fossil evidence for a late Cretaceous origin of "hoofed" mammals Nature 272:1150-1153

    Archibald J. D., A. O. Averianov, E. G. Ekdale, 2001 Late Cretaceous relatives of rabbits, rodents, and other extant eutherian mammals Nature 414:62-65[ISI][Medline]

    Brandt J. F., 1855 Beiträge zur nahern Kenntniss der Säugethiere Russlands Mem. Acad. Imp. St. Petersbourg Ser 69:1-375

    Bromham L., M. J. Phillips, D. Penny, 1999 Growing up with dinosaurs: molecular dates and the mammalian radiation Trends Ecol. Evol 14:113-118[ISI][Medline]

    Bugge J., 1985 Systematic value of the carotid arterial pattern in rodents Pp. 381–402 in W. P. Luckett and J.-L. Hartenberger, eds. Evolutionary relationships among rodents: a multidisciplinary analysis. Plenum Press, New York

    Cao Y., M. Fujiwara, M. Nikaido, N. Okada, M. Hasegawa, 2000 Interordinal relationships and timescale of eutherian evolution as inferred from mitochondrial genome data Gene 259:149-158[ISI][Medline]

    Catzeflis F. M., C. Hänni, P. Sourrouille, E. Douzery, 1995 Molecular systematics of hystricognath rodents: the contribution of sciurognath mitochondrial 12S rRNA sequences Mol. Phylogenet. Evol 4:357-360[ISI][Medline]

    DeBry R. W., R. M. Sagel, 2001 Phylogeny of Rodentia (Mammalia) inferred from the nuclear-encoded gene IRBP Mol. Phylogenet. Evol 19:290-301[ISI][Medline]

    Delsuc F., F. M. Catzeflis, M. Stanhope, E. J. P. Douzery, 2001 The evolution of armadillos, anteaters and sloths depicted by nuclear and mitochondrial phylogenies: implications for the status of the enigmatic fossil Eurotamandua Proc. R. Soc. Lond. B 268:1605-1615[ISI][Medline]

    D'Erchia A. M., C. Gissi, G. Pesole, C. Saccone, U. Arnason, 1996 The guinea-pig is not a rodent Nature 381:597-600[ISI][Medline]

    Eizirik E., W. J. Murphy, S. J. O'Brien, 2001 Molecular dating and biogeography of the early placental mammal radiation J. Hered 92:212-219[Abstract/Free Full Text]

    Felsenstein J., 1985 Confidence limits on phylogenies: an approach using the bootstrap Evolution 39:783-791[ISI]

    ———. 1988 Phylogenies from molecular sequences: inference and reliability Ann. Rev. Genet 22:521-565[ISI][Medline]

    ———. 1995 PHYLIP (phylogeny inference package) Version 3.573. Distributed over the Internet (http://evolution.genetics.washington.edu)

    Gheerbrant E., J. Sudre, H. Cappetta, 1996 A Palaeocene proboscidean from Morocco Nature 383:68-70.[ISI]

    Gingerich P. D., M. D. Uhen, 1998 Likelihood estimation of the time of origin of Cetacea and the time of divergence of Cetacea and Artiodactyla Paleontol. Electronica 2:1-47

    Gissi C., A. Gullberg, U. Arnason, 1998 The complete mitochondrial DNA sequence of the rabbit, Oryctolagus cuniculus Genomics 50:161-169[ISI][Medline]

    Graur D., W. A. Hide, W.-H. Li, 1991 Is the guinea-pig a rodent? Nature 351:649-652[ISI][Medline]

    Hartenberger J.-L., 1985 The order Rodentia: major questions on their evolutionary origin, relationships and suprafamilial systematics Pp. 1–33 in W. P. Luckett and J.-L. Hartenberger, eds. Evolutionary relationships among rodents: a multidisciplinary analysis. Plenum Press, New York

    ———. 1994 The evolution of the Gliroidea Pp. 19–33 in Y. Tomida, C.-K. Li, and T. Setoguchi, eds. Rodent and lagomorph families of Asian origins and diversification. National Science Museum Monographs, Tokyo

    ———. 1998 Description de la radiation des Rodentia (Mammalia) du Paléocène supérieur au Miocène; incidences phylogénétiques C. R. Acad. Sci. Paris, Sciences de la terre et des planètes 326:439-444

    Hassanin A., G. Lecointre, S. Tillier, 1998 The "evolutionary signal" of homoplasy in protein-coding gene sequences and its consequences for a priori weighting in phylogeny C. R. Acad. Sci. Paris, Life sci 321:611-620

    Hedges S. B., P. H. Parker, C. G. Sibley, S. Kumar, 1996 Continental breakup and the ordinal diversification of birds and mammals Nature 381:226-229[ISI][Medline]

    Henke W., K. Herdel, K. Jung, D. Schnorr, S. A. Loening, 1997 Betaine improves the PCR amplification of GC-rich DNA sequences Nucleic Acids Res 25:3957-3958[Abstract/Free Full Text]

    Huchon D., F. M. Catzeflis, E. J. P. Douzery, 1999 Molecular evolution of the nuclear von Willebrand Factor gene in mammals and the phylogeny of rodents Mol. Biol. Evol 16:577-589[Abstract]

    Huchon D., F. M. Catzeflis, E. J. P. Douzery, 2000 Variance of molecular datings, evolution of rodents, and the phylogenetic affinities between Ctenodactylidae and Hystricognathi Proc. R. Soc. Lond. B 267:393-402[ISI][Medline]

    Huchon D., E. J. P. Douzery, 2001 From the Old World to the New World: a molecular chronicle of the phylogeny and biogeography of hystricognath rodents Mol. Phylogenet. Evol 20:238-251[ISI][Medline]

    Huelsenbeck J. P., F. Ronquist, 2001 MrBayes: Bayesian inference of phylogenetic trees Bioinformatics 17:754-755[Abstract/Free Full Text]

    Jacobs L. L., W. R. Downs, 1994 The evolution of murine rodents in Asia Pp. 149–156 in Y. Tomida, C.-K. Li, and T. Setoguchi, eds. Rodent and lagomorph families of Asian origins and diversification. National Science Museum Monographs, Tokyo

    Jaeger J.-J., 1988 Rodent phylogeny: new data and old problems Pp. 177–199 in M. J. Benton, ed. The phylogeny and classification of the Tetrapods. Clarendon Press, Oxford

    Janke A., X. Xu, U. Arnason, 1997 The complete mitochondrial genome of the wallaroo (Macropus robustus) and the phylogenetic relationship among Monotremata, Marsupialia, and Eutheria Proc. Natl. Acad. Sci. USA 94:1276-1281[Abstract/Free Full Text]

    Kishino H., M. Hasegawa, 1989 Evaluation of the maximum likelihood estimate of the evolutionary tree topologies from DNA sequence data, and the branching order in Hominoidea J. Mol. Evol 29:170-179[ISI][Medline]

    Krajewski C., M. G. Fain, L. Buckley, D. G. King, 1999 Dynamically heterogenous partitions and phylogenetic inference: an evaluation of analytical strategies with cytochrome b and ND6 gene sequences in cranes Mol. Phylogenet. Evol 13:302-313[ISI][Medline]

    Kumar S., S. B. Hedges, 1998 A molecular timescale for vertebrate evolution Nature 392:917-920[ISI][Medline]

    Lavocat R., J.-P. Parent, 1985 Phylogenetic analyses of middle ear features in fossil and living rodents Pp. 333–354 in W. P. Luckett, and J.-L. Hartenberger, eds. Evolutionary relationships among rodents: a multidisciplinary analysis. Plenum Press, New York

    Li W.-H., W. A. Hide, A. Zharkikh, D.-P. Ma, D. Graur, 1992 The molecular taxonomy and evolution of the guinea-pig J. Hered 83:174-181[ISI][Medline]

    Luckett W. P., J.-L. Hartenberger, 1985 Evolutionary relationships among rodents: comment and conclusions Pp. 227–276 in W. P. Luckett and J.-L. Hartenberger, eds. Evolutionary relationships among rodents: a multidisciplinary analysis. Plenum Press, New York

    ———. 1993 Monophyly or polyphyly of the order Rodentia: possible conflict between morphological and molecular interpretations J. Mamm. Evol 1:127-147

    Madsen O., M. Scally, C. J. Douady, D. J. Kao, R. W. DeBry, R. Adkins, H. Amrine, M. J. Stanhope, W. W. de Jong, M. S. Springer, 2001 Parallel adaptative radiations in two major clades of placental mammals Nature 409:610-614[ISI][Medline]

    Matthee C. A., T. J. Robinson, 1997 Molecular phylogeny of the springhare, Pedetes capensis, based on mitochondrial DNA sequences Mol. Biol. Evol 14:20-29[Abstract]

    McKenna M. C., 1986 Glirology Science 231:1666-1667

    McKenna M. C., S. K. Bell, 1997 Classification of mammals above the species level Columbia University Press, New York

    Meng J., 1990 The auditory region of Reithroparamys delicatissimus (Mammalia, Rodentia) and its systematic implications Am. Mus. Novit 2972:1-35

    Montgelard C., S. Bentz, C. Douady, J. Lauquin, F. M. Catzeflis, 2002a. Molecular phylogeny of the sciurognath rodent families Gliridae, Anomaluridae and Pedetidae: morphological and paleontological implications 8th ASM Symposium, Paris, Edition IRD, Collection colloques et séminaires (in press)

    Montgelard C., S. Bentz, C. Tirard, O. Verneau, F. M. Catzeflis, 2002b. Molecular systematics of Sciurognathi (Rodentia): the mitochondrial cytochrome b and 12S rRNA genes support the Anomaluroidea (Pedetidae and Anomaluridae) Mol. Phylogenet. Evol 22:220-233[ISI][Medline]

    Mouchaty S. K., F. M. Catzeflis, A. Janke, U. Arnason, 2001 Molecular evidence of an African Phiomorpha-South American Caviomorpha clade and support for Hystricognathi based on the complete mitochondrial genome of the cane rat (Thryonomys swinderianus) Mol. Phylogenet. Evol 18:127-135[ISI][Medline]

    Murphy W. J., E. Eizirik, W. E. Johnson, Y. P. Zhang, O. A. Ryder, S. J. O'Brien, 2001a. Molecular phylogenetics and the origins of placental mammals Nature 409:614-618[ISI][Medline]

    Murphy W. J., E. Eizirik, S. J. O'Brien, et al. (11 co-authors) 2001b. Resolution of the early placental mammal radiation using Bayesian phylogenetics Science 294:2348-2351[Abstract/Free Full Text]

    Nedbal M. A., R. L. Honeycutt, D. A. Schlitter, 1996 Higher-level systematics of rodents (Mammalia, Rodentia): evidence from the mitochondrial 12S rRNA gene J. Mamm. Evol 3:201-237

    Philippe H., 1997 Rodent monophyly: pitfalls of molecular phylogenies J. Mol. Evol 45:712-715[ISI][Medline]

    Porter A. P., M. Goodman, M. J. Stanhope, 1996 Evidence on mammalian phylogeny from sequences of exon 28 of the von Willebrand Factor gene Mol. Phylogenet. Evol 5:89-101[ISI][Medline]

    Rambaut A., L. Bromham, 1998 Estimating divergence dates from molecular sequences Mol. Biol. Evol 15:442-448[Abstract]

    Reyes A., C. Gissi, G. Pesole, F. M. Catzeflis, C. Saccone, 2000 Where do rodents fit? Evidence from the complete mitochondrial genome of Sciurus vulgaris Mol. Biol. Evol 17:979-983[Free Full Text]

    Reyes A., G. Pesole, C. Saccone, 1998 Complete mitochondrial DNA sequence of the fat dormouse, Glis glis: further evidence of rodent paraphyly Mol. Biol. Evol 15:499-505[Abstract]

    Robinson M., F. Catzeflis, J. Briolay, D. Mouchiroud, 1997 Molecular phylogeny of rodents, with special emphasis on murids: evidence from nuclear gene LCAT Mol. Phylogenet. Evol 8:423-434[ISI][Medline]

    Robinson-Rechavi M., L. Ponger, D. Mouchiroud, 2000 Nuclear gene LCAT support rodent monophyly Mol. Biol. Evol 17:1410-1412[Free Full Text]

    Ryan J. M., 1989 Comparative myology and phylogenetic systematics of the Heteromyidae (Mammalia, Rodentia) Misc. Publ. Mus. Zool., Univ. Mich 176:1-103

    Shimodaira H., M. Hasegawa, 1999 Multiple comparisons of log-likelihoods with applications to phylogenetic inference Mol. Biol. Evol 16:1114-1116[Free Full Text]

    Springer M. S., H. Amrine, A. Burk, M. Stanhope, 1999 Additional support for Afrotheria and Paenungulata, the performance of mitochondrial versus nuclear genes, and the impact of data partitions with heterogeneous base composition Syst. Biol 48:65-75[ISI][Medline]

    Springer M. S., A. Burk, J. R. Kavanagh, V. G. Waddell, M. J. Stanhope, 1997a. The Interphotoreceptor Retinoid Binding Protein gene in therian mammals: implications for higher level relationships and evidence for loss of function in the marsupial mole Proc. Natl. Acad. Sci. USA 94:13754-13759[Abstract/Free Full Text]

    Springer M. S., G. C. Cleven, O. Madsen, W. W. de Jong, V. G. Waddell, H. M. Amrine, M. J. Stanhope, 1997b. Endemic african mammals shake the phylogenetic tree Nature 388:61-64[ISI][Medline]

    Springer M. S., R. W. DeBry, C. Douady, H. M. Amrine, O. Madsen, W. W. de Jong, M. J. Stanhope, 2001 Mitochondrial versus nuclear gene sequences in deep-level mammalian phylogeny reconstruction Mol. Biol. Evol 18:132-143[Abstract/Free Full Text]

    Stanhope M. J., J. Czelusniak, J.-S. Si, J. Nickerson, M. Goodman, 1992 A molecular perspective on mammalian evolution from the gene encoding interphotoreceptor retinoid binding protein, with convincing evidence for bat monophyly Mol. Phylogenet. Evol 1:148-160[Medline]

    Strimmer K., A. von Haeseler, 1996 Quartet puzzling: a quartet maximum-likelihood method for reconstructing tree topologies Mol. Biol. Evol 13:964-969[Free Full Text]

    Sullivan J., D. L. Swofford, 1997 Are guinea pigs rodents? The importance of adequate models in molecular phylogenetics J. Mamm. Evol 4:77-86

    Swofford D. L., 1999 PAUP* Phylogenetic analysis using parsimony (* and other methods). Version 4. Sinauer Associates, Sunderland, Mass

    Tabuce R., B. Coiffait, P.-E. Coiffait, M. Mahboubi, J.-J. Jaeger, 2001 A new genus of Macroscelidea (Mammalia) from the Eocene of Algeria: a possible origin for elephant-shrews J. Vertebr. Paleontol 21:535-546[ISI]

    Thorne J. L., H. Kishino, I. S. Painter, 1998 Estimating the rate of evolution of the rate of molecular evolution Mol. Biol. Evol 15:1647-1657[Abstract/Free Full Text]

    Tullberg T., 1899 Ueber das System der Nagetiere: Eine phylogenetische Studie Nova Acta Reg. Soc. Sci. Upsala Ser. 3 18:1-514

    Vianey-Liaud M., 1985 Possible evolutionary relationships among Eocene and Lower Oligocene rodents of Asia, Europe and North America Pp. 277–309 in W. P. Luckett and J.-L. Hartenberger, eds. Evolutionary relationships among rodents: a multidisciplinary analysis. Plenum Press, New York

    Waddell P. J., Y. Cao, M. Hasegawa, D. Mindell, 1999 Assessing the Cretaceous superordinal divergence times within birds and placental mammals by using whole mitochondrial protein sequences and an extended statistical framework Syst. Biol 48:119-137[ISI][Medline]

    Wahlert J. H., 1985 Skull morphology and relationships of geomyoid rodents Am. Mus. Novit 2812:1-20

    Walton A. H., 1997 Rodents Pp. 392–409 in R. F. Kay, R. H. Madden, R. L. Ciffelli, and J. J. Flynn, eds. Vertebrate paleontology in the neotropics. The Miocene fauna of La Venta, Columbia. Smithsonian Institution Press, Washington

    Wyss A. R., J. J. Flynn, M. A. Norell, C. C. Swisher III, R. Charrier, M. J. Novacek, M. C. McKenna, 1993 South America's earliest rodent and recognition of a new interval of mammalian evolution Nature 365:434-437[ISI]

    Yang Z., 1997 PAML: a program package for phylogenetic analysis by maximum likelihood CABIOS 13:555-556[Medline]

Accepted for publication February 11, 2002.