*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;
Department of Biochemistry, University of Nijmegen, The Netherlands;
Biology and Biochemistry, Queen's University of Belfast, UK
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
The relationships between rodent families are also disputed. Morphological approaches have been frustrated by convergent evolution of characters (e.g., Jaeger 1988
), 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 1985
; Nedbal, Honeycutt, and Schlitter 1996
; Huchon, Catzeflis, and Douzery 2000
; Adkins et al. 2001
). However, the monophyly of most rodent families seems well established (Hartenberger 1985
). 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. 1995
; Nedbal, Honeycutt, and Schlitter 1996
; Matthee and Robinson 1997
; Robinson et al. 1997
; Huchon, Catzeflis, and Douzery 1999
, 2000
; Robinson-Rechavi, Ponger, and Mouchiroud 2000
; Adkins et al. 2001
; DeBry and Sagel 2001
; Murphy et al. 2001a;
Montgelard et al. 2002b
). For example, singular families like Ctenodactylidae or Anomaluridae are rarely sampled. In the present study we include for the first time all major rodent lineagesnamely all sciurognath families and hystricognath superfamiliesthe 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 1998
), whereas molecular clocks based on a few rodent species tend to support a Cretaceous origin and diversification of rodents, 89125 MYA for the divergence of murids or hystricognaths (Janke, Xu, and Arnason 1997
; Kumar and Hedges 1998
; Cao et al. 2000
). 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 2000
; Adkins et al. 2001
). 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 1998
).
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 1996
; Springer et al. 1997a,
1997b;
DeBry and Sagel 2001
; Delsuc et al. 2001
; Madsen et al. 2001
). 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. 2001
). 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 |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
DNA extractions and PCR reactions were conducted as described before (Huchon, Catzeflis, and Douzery 1999
), with slight modifications: 1 M betaine was included in the PCR mixture (Henke et al. 1997
), 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)
, Springer et al. (1997b
), and Huchon, Catzeflis, and Douzery (1999,
2000)
, 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 [
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 (AJ427226AJ427270) 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 (8).
Neighbor-Joining (NJ) analyses with weighted average (WAVE) maximum likelihood distances were performed with PHYLIP 3.573 (Felsenstein 1995
) and WAVEBOOT 1.2 (Krajewski et al. 1999
). 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 1999
), version 4, releases beta 4 and 8; and TREE-PUZZLE 4.0.2 (Strimmer and von Haeseler 1996
). 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 1998
).
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 () and the parameter (
) of the Gamma distribution of the HKY +
8 model. For amino acid sequences, ML reconstructions under the JTT +
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 2001
). 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-WAVEwhere each codon position had its own rateand MP analyses.
Robustness of the nodes of the phylogenetic trees was assessed by bootstrap (Felsenstein 1985
) 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 +
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 + 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 +
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 1997
), version 3.0d. In the latter case, all ML parametersi.e., 6 transition-transversion rate parameters, 6
parameters, and 6 x 81 branch parameterswere reestimated by PAML for each evaluated topology. Partitioned log-likelihoods were then compared using the Kishino and Hasegawa (1989)
test with the Shimodaira and Hasegawa (1999)
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 + 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)
test with the Shimodaira and Hasegawa (1999)
correction was performed.
Molecular Dating
Divergence times were estimated using quartet dating (Rambaut and Bromham 1998
), 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
.
|
|
![]() |
Results and Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
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).
|
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 + 8 model, the log-likelihood of the best topology was lnL = -26,415.56. Assuming six independent HKY +
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 bywith increasing relative rateA2AB (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.
|
|
|
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 2000
) contains gundi and Guinea pigrelated 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 1985
; Meng 1990
) and molecular data (Huchon, Catzeflis, and Douzery 1999
, 2000
; Adkins et al. 2001
; DeBry and Sagel 2001
). The increase in taxonomic sampling did not reduce the support for this clade (BPML = 100; RPPROT = 85). In agreement with DeBry and Sagel (2001)
, 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 1985
; Meng 1990
) and mitochondrial (Reyes et al. 2000
) studies; however, similar to other molecular studies (e.g., Huchon, Catzeflis, and Douzery 1999
; Adkins et al. 2001
; DeBry and Sagel 2001
; Montgelard et al. 2002b
), 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 1855
; Tullberg 1899
). However, some morphological studies could not confirm this relationship (Bugge 1985
; Lavocat and Parent 1985
; Meng 1990
), and it has even been suggested that Castoridae might be more closely related to Muridae than to Sciuridae (Meng 1990
). Anomaluridae and Pedetidae were considered as enigmatic families, possibly related to Ctenodactylidae or Hystricognathi (e.g., Luckett and Hartenberger 1985
; Jaeger 1988
). Geomyoidea has, less ambiguously, been associated with the Muridae (e.g., Wahlert 1985
; Ryan 1989
).
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. 2001
; Murphy et al. 2001a
). 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 1985
; Hartenberger 1998
). 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 inor related toEarly- and Middle-Eocene families: Cricetidae, Zapodidae, Zegdoumyidae, Eomyidae, and Eutypomyidae, respectively (Hartenberger 1998
).
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 1985
; Nedbal, Honeycutt, and Schlitter 1996
; Huchon, Catzeflis, and Douzery 1999
; DeBry and Sagel 2001
), and (2) Anomaluromorpha, i.e., scaly-tailed flying squirrels and springhares (Bugge 1985
; Lavocat and Parent 1985
; Montgelard et al. 2002a
). The third subclade unexpectedly clusters beavers with pocket mice and pocket gophers (Adkins et al. 2001
). 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. 2001
), on a combined data set of 15 genes (Murphy et al. 2001a
), and on 12S rRNA and cytochrome b sequences (Montgelard et al. 2002b
). In contrast, DeBry and Sagel (2001)
found different results with the IRBP gene: Geomyoidea are the sister clade of Ctenohystrica and Castoridae clusters with Sciuroideaan 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 2000
) and the GHR gene (Adkins et al. 2001
). 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)
. 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)
, rodents have been divided into two subordersSciurognathi 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 1855
). The contents of these groups changed according to the authors, and we here follow the classification of McKenna and Bell (1997)
. 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 2000
). 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.353.2 Myr (vWF) against 66.9137.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 1998
), 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).
|
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.463.7; fig. 4
). Paleontological studies place the rodent radiation at 55 MYA at the Paleocene-Eocene boundary (Hartenberger 1998
). The split between rodents and lagomorphs is estimated to be 64.5 MYA (57.373.3), just at the Cretaceous-Tertiary boundary. These results indicate that rodents and lagomorphs diverged only at the beginning of the Tertiary (fig. 4
).
|
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. 1996
; Alroy 1999
; Eizirik, Murphy, and O'Brien 2001
; Madsen et al. 2001
). 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 |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
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 20002002 (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 |
---|
1 Present address: Molecular Evolution Laboratory, Faculty of Bioscience and Biotechnology, Tokyo Institute of Technology, Nagatsuta-cho, Midori-ku, Yokohama, Japan
2 Present address: Bioinformatics, GlaxoSmithKline, UP1345, Collegeville Pennsylvania
Keywords: Rodentia
Glires
Eutheria
vWF
IRBP
A2AB
phylogeny
molecular dating
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
![]() |
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
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. 381402 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
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. 133 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. 1933 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
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
Jacobs L. L., W. R. Downs, 1994 The evolution of murine rodents in Asia Pp. 149156 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. 177199 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
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. 333354 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. 227276 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
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
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
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
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
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
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
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
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. 277309 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. 392409 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]