Evolutionary History of the Most Speciose Mammals: Molecular Phylogeny of Muroid Rodents

Johan Michaux, Aurelio Reyes and François Catzeflis

Laboratoire de Paléontologie, Institut des Sciences de l'Evolution, UMR 5554 Centre National de la Recherche Scientifique, Université de Montpellier 2, Montpellier, France;
Centro di Studio sui Mitocondri e Metabolismo Energetico, Consiglio Nazionale delle Ricerche, Bari, Italy


    Abstract
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions
 Acknowledgements
 References
 
Phylogenetic relationships between 32 species of rodents representing 14 subfamilies of Muridae and four subfamilies of Dipodidae were studied using sequences of the nuclear protein-coding genes Lecithin Cholesterol Acyl Transferase (LCAT) and von Willebrand Factor (vWF). An examination of some evolutionary properties of each data matrix indicates that the two genes are rather complementary, with lower rates of nonsynonymous substitutions for LCAT. Both markers exhibit a wide range of GC3 percentages (55%–89%), with several taxa above 70% GC3 for vWF, which indicates that those exonic regions might belong to the richest class of isochores. The primary sequence data apparently harbor few saturations, except for transitions on third codon positions for vWF, as indicated by comparisons of observed and expected pairwise values of substitutions. Phylogenetic trees based on 1,962 nucleotidic sites from the two genes indicate that the 14 Muridae subfamilies are organized into five major lineages. An early isolation leads to the clade uniting the fossorial Spalacinae and semifossorial Rhizomyinae with a strong robustness. The second lineage includes a series of African taxa representing nesomyines, dendromurines, cricetomyines, and the sole living member of mystromyines. The third one comprises only the mouselike hamster Calomyscus. The fourth clade represents the cricetines, myospalacines, sigmodontines, and arvicolines, whereas the fifth one comprises four "traditional" subfamilies (Gerbillinae, Murinae, Otomyinae, and Acomyinae). Within these groups, we confirm the monophyly of almost all studied subfamilies, namely, Spalacinae, Rhizomyinae, Nesomyinae, Cricetomyinae, Arvicolinae, Sigmodontinae, Cricetinae, Gerbillinae, Acomyinae, and Murinae. Finally, we present evidence that the sister group of Acomyinae is Gerbillinae, and we confirm a nested position of Myospalacinae within Cricetinae and Otomyinae within Murinae. From a biogeographical point of view, the five main lineages spread and radiated from Asia with different degrees of success: the first three groups are now represented by a limited number of species and genera localized in some regions, whereas the last two groups radiated in a large variety of species and genera dispersed all over the world.


    Introduction
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions
 Acknowledgements
 References
 
The rodents of the family Muridae are the most diverse group of mammals, encompassing at least 1,326 species spanning more than 281 genera (Musser and Carleton 1993Citation ). The evolutionary systematics of this family have been very difficult, and in spite of many attempts (Miller and Gidley 1918Citation ; Simpson 1945Citation ; Hooper and Musser 1964Citation ; Chaline, Mein, and Petter 1977Citation ; Carleton and Musser 1984Citation ), several uncertainties, confusions, and conflicting views have persisted for these animals. For this reason, in their recent review, Musser and Carleton (1993)Citation decided to keep a prudent state of uncertainty about the hierarchical pattern of muroid suprageneric groups and to divide the family Muridae into 17 subfamilies considered at the same taxonomic level. These "major lineages" within murids are Arvicolinae (26 genera/143 species), Calomyscinae (1/6), Cricetinae (7/18), Cricetomyinae (3/6), Dendromurinae (8/23), Gerbillinae (14/110), Lophiomyinae (1/1), Murinae (122/529), Myospalacinae (1/7), Mystromyinae (1/1), Nesomyinae (7/14), Otomyinae (2/14), Petromyscinae (2/5), Platacanthomyinae (2/3), Rhizomyinae (3/15), Sigmodontinae (79/423), and Spalacinae (2/8).

However, there remains a strong need to better define the taxonomic boundaries of these subfamilies and especially the relationships existing between them. Many important questions have not yet been adequately answered concerning the evolutionary origins of most of the 17 subfamilies or of their rates of evolution.

In a recent publication based on 40 rodent sequences of the nuclear protein-coding gene Lecithin Cholesterol Acyl Transferase (LCAT) (Michaux and Catzeflis 2000Citation ) representing 13 Muridae subfamilies, we proposed a first step toward a better knowledge of this speciose group. The results showed that within Muridae, a first branching leads to the fossorial Spalacinae and the semifossorial Rhyzomyinae. The remaining components of Muridae appeared as a polytomy from which Sigmodontinae, Calomyscinae, Arvicolinae, Cricetinae, Mystromyinae, Nesomyinae and some Dendromurinae (Steatomys and Dendromus) were issued. This phylogeny was interpreted by Michaux and Catzeflis (2000)Citation as the result of a bushlike radiation at the end of the early Miocene, leading to emergence of most living subfamilies. The separation between three additional taxa, Murinae, Gerbillinae, and "Acomyinae" (which comprises the genera Acomys, Deomys, Uranomys, and Lophuromys), occurred more recently, from a common ancestor issued from the main basal radiation. As previously shown by other molecular studies, the vlei rats Otomyinae are nested within Old World Murinae. In the same way, the zokors Myospalacinae appear strongly nested within the hamsters Cricetinae. Finally, Michaux and Catzeflis (2000)Citation proposed a sister group relationship between Malagasy Nesomyinae and South African Mystromyinae.

In order to verify these hypotheses and to know if the majority of the modern Muridae subfamilies indeed appeared as a bushlike radiation, we here tested the relationships between 14 Muridae subfamilies using another nuclear gene: the von Willebrand Factor (vWF). Then, we combined these new sequences with those of the LCAT gene.


    Materials and Methods
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions
 Acknowledgements
 References
 
DNA used to sequence exon 28 of the vWF gene was extracted and purified from ethanol-preserved tissues taken from the collection of mammal tissues housed at Montpellier (Catzeflis 1991Citation ). Whenever possible, we selected two species for each studied subfamily (see table 1 ). This biological sampling was aimed at getting an equilibrated representation of each murid lineage diminishing a possible "long-branch attraction effect." These species were also chosen according to the sampling already performed on the nuclear LCAT gene (Michaux and Catzeflis 2000Citation ) in order to combine the sequences of both genes.


View this table:
[in this window]
[in a new window]
 
Table 1 References for Rodent Tissues Used in the Experiments

 
DNA Sequencing of vWF Gene Exon 28
DNA extraction from the 95% ethanol-preserved tissues was performed according to Sambrook, Fritsch, and Maniatis (1989)Citation . The main part of exon 28 from the vWF gene (1,265 bp) was amplified using the PCR primers V1 (5'-TGTCAACCTCACCTGTGAAGCCTG-3') and W1 (5'-TGCAGGACCAGGTCAGGAGCCTCTC-3') previously designed by Huchon, Catzeflis, and Douzery (1999)Citation . All PCRs used 5 min at 94°C; 33 cycles of 45 s at 94°C, 30 s at 52°C, and 1 min at 72°C; and 10 min at 72°C in an Appligen Crocodile 3 or a Labover PTC100 thermal cycler. The total reaction volume was 50 µl. PCR products were purified using the Ultra-free DNA Amicon kit (Millipore) and directly sequenced. Sequencing on both strands was done using a dye terminator or Big dye terminator (Perkin Elmer) sequencing kit and ABI 373 and ABI 310 (Perkin Elmer) automatic sequencers. The external (V1 and W1) and internal (V2 and W2) primers designed by Huchon, Catzeflis, and Douzery (1999)Citation were used for sequencing.

Technical difficulties prevented us from obtaining a high-quality sequence of Tachyoryctes for the LCAT gene, whereas we obtained only a partial sequence (624 nt) for Rhizomys exon 28 of the vWF gene. Thus, our representation of the taxon Rhizomyinae is chimeric: complete sequence of vWF for Tachyoryctes and of LCAT for Rhizomys. We nevertheless checked that Tachyoryctes and Rhizomys were sister taxa (100% bootstrap support; data not shown) when the only 624 available sites of vWF were compared for 33 Myodonta taxa. Thus, we feel confident that, as documented by comparative morphology of fossil and living taxa (Flynn 1990Citation ), Rhizomys and Tachyoryctes are sister genera among all other rodents of our study.

Sequence Alignment and Saturation Analysis
Previously known sequences for vWF and LCAT genes were extracted from GenBank and aligned with the new sequences using CLUSTAL W (Thompson, Higgins, and Gibson 1994Citation ) and the ED editor (MUST package; Philippe 1993Citation ). The program AFAS (MUST package; Philippe 1993Citation ) was used to combine the aligned matrices of vWF and LCAT.

We performed a saturation analysis for vWF as described by Philippe and Douzery (1994)Citation and Hassanin, Lecointre, and Tillier (1998)Citation . Details for such an analysis are provided by Michaux and Catzeflis (2000)Citation .

Phylogenetic Reconstructions
The aligned sequences were treated with different approaches: the stationary Markov model (Saccone et al. 1990Citation ) (also called general time reversible [GTR]) in PAUP 4.0b2 (Swofford 1998Citation ), the Tamura and Nei (1993)Citation model, and the LogDet (Lockhaert et al. 1994Citation ) estimator were used for the calculation of genetic distances. The last estimator was used to take into account the differences in GC composition between species. The GTR analyses were also performed assuming a gamma distribution for substitution rates across sites, where the parameter alpha (Yang 1996Citation ) and the proportion of invariant sites (I) were estimated with the maximum-likelihood (ML) method assuming the GTR–Markov evolution (ME) phylogeny using PAUP 4.0b2. Maximum-parsimony (MP; heuristic search, TBR branch swapping option) and ML (GTR model of sequence evolution) analyses were also conducted using PAUP 4.0b2 (Swofford 1998).

The robustness of inferences was assessed through bootstrap resampling (BP) under ML after 100 replicates (with neighbor-joining [NJ] starting trees, NNI branch swapping, and model parameters fixed to values estimated from the original data) and under MP and NJ or ME after 1,000 replicates. Bremer's support index Citation (BSI, Bremer 1988) was also calculated on the most parsimonious tree with enforcement of topological constraints. Alternative topologies were evaluated by the Kishino and Hasegawa (1989)Citation test implemented in PAUP 4.0b2.

The level of incongruence between the two genes was tested using PAUP4 (option Hompart; Swofford 1998). This approach used the incongruence length difference (ILD) test with the parsimony criterion (Farris et al. 1995Citation ). One thousand randomizations were performed on variable sites only (Cunningham 1997Citation ).

Relative-Rate Test
Relative-rate tests were conducted both with RRTree, version 1.0 (Robinson et al. 1998Citation ) (which improves the test of Wu and Li [1985]Citation ), by taking into account the taxonomic representativity and its phylogenetic relationships, and with the likelihood approach of Muse and Gaut (1994)Citation . In the latter case, statistical significance of the test was assessed by means of the chi-square test (P < 0.05). Quantifications of the rate differences were performed either on the proportions of synonymous (Ks) and nonsynonymous (Ka) substitutions for RRTree, or on all three codon positions and GTR model for the likelihood approach. Relative-rate tests were performed among rodents at intrafamilial levels; Dipodidae were chosen as outgroups.


    Results
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions
 Acknowledgements
 References
 
Sequenced Species
The 27 new rodent sequences of vWF gene exon 28 have been deposited in the EMBL gene bank under accession numbers AJ297764, AJ297765, and AJ402693–AJ402717. The newly determined sequences were compared with four rodent sequences determined by Huchon, Catzeflis, and Douzery (1999)Citation , as well as with a Mus musculus sequence available in GenBank (U27810).

Evolutionary Properties of LCAT and Exon 28 vWF Genes
The data matrix for 32 Myodonta taxa had 1,962 nucleotidic sites, of which 711 (about 36%) were parsimony-informative. Exons of LCAT comprised 804 sites; about half of them were variable, and 239 (30%) were parsimony-informative; exon 28 of the vWF gene included 1,158 positions, of which 710 (61%) were variable and 472 (41%) were parsimony-informative. Thus, there was more absolute and relative variation in vWF than in LCAT. The average ratio of transitions to transversions was 2.00 for LCAT, ranging from 1.1 (Uranomys ruddi/M. musculus comparison) to 5.00 (Allactaga elater/Dipus sagitta). For vWF, this ratio was 1.94, ranging from 1.23 (Saccostomus campestris/D. sagitta and S. campestris/Jaculus jaculus comparisons) to 3.52 (Mesocricetus auratus/Myospalax sp.).

From comparing the 496 pairs of sequences, it appears that LCAT exons have evolved slower than exon 28 of vWF in muroids and dipodoids. The LCAT/vWF ratios for nucleotidic and amino acid percentages of differences were 0.78% ± 0.16% and 0.63% ± 0.21%, respectively. Nonsynonymous changes (Ka values) were much less frequent in LCAT (values of up to 0.08), on average 59% ± 20% of the values observed in vWF (values of up to 0.15); synonymous changes were as frequent in LCAT as in vWF (values of up to 0.45 and 0.55, respectively), and this might well be due to saturation of both genes for this kind of change (fig. 1 ).



View larger version (24K):
[in this window]
[in a new window]
 
Fig. 1.—Synonymous (Ks) and nonsynonymous (Ka) changes for 496 pairs of rodent sequences in the LCAT gene and exon 28 of the vWF gene. In the two curves, data are organized by increasing values in order to illustrate (A) the similar rates of synonymous change in both genes and (B) the faster rate of nonsynonymous change in vWF sequences

 
The nucleotidic and amino acid compositions were different for each gene (P < 0.05; nonparametric Mann-Whitney t-tests) at each codon position (nucleotides) and for each amino acid except proline (P = 0.47). LCAT had much more T in first codon positions and A and T in third positions, whereas vWF contained much more A in first positions and G in third positions. The most pronounced differences in frequencies of amino acids concerned an excess of C, N, D, and Y as opposed to a deficiency of K, I, S, A, and E for LCAT. Exon 28 of vWF had a very low representation of amino acids C (0.51% ± 0.06%) and W (0.03% ± 0.09%) (data not shown).

GC3 Content
The percentages of GC at third codon positions (GC3) spanned a much larger scatter in vWF (64.1–89.1) than in LCAT (54.8–67.6) (data available on request). On average, the four dipodoids had a higher GC3 (81.7 ± 8.3 for vWF; 64.4 ± 3.1 for LCAT) than the 28 muroids (68.6 ± 2.9 for vWF; 59.7 ± 3.1 for LCAT) (unpaired Mann-Whitney t-tests, P < 0.009). These differences held also for vWF (P < 0.005) in comparisons of dipodoids and several subsets of muroids, such as Gerbillinae + Murinae + Acomyinae, or Arvicolinae + Sigmodontinae + Cricetinae, or Nesomyinae + Cricetomyinae + Dendromurinae. For LCAT, such comparisons also indicated that advanced muroids had a lower GC content, although the differences were marginally significant (P values of 0.070, 0.019, and 0.052, respectively). We here confirm, based on a larger data set, that the murids have a lower GC3 value for LCAT (Robinson, Gautier, and Mouchiroud 1997Citation ), and demonstrate that for vWF the average GC3 difference is still larger for murids and nonmurids.

Saturation Analyses
The levels of saturation for LCAT had already been assessed by Michaux and Catzeflis (2000)Citation , showing that the C-T and A-G transitions exhibit a lower slope (average of 0.33 for A-G and 0.35 for C-T) with regard to the transversions, whatever the codon position. Consequently, Michaux and Catzeflis (2000)Citation performed a weighted analysis according to the slopes of each substitution. However, although the slopes were significantly lower than those observed for vWF (see below), the results of the weighted analysis did not show any interesting differences with regard to the unweighted parsimony analysis. Thus, all substitutions were retained for phylogenetic reconstruction performed with the LCAT gene.

The levels of saturation for vWF were moderate on each kind of nucleotide substitution and position (slopes of regression analysis of 0.51 for TS1 and 0.54 for TS2 and of 0.86, 0.87, and 0.87, respectively, for TV1, TV2 and TV3), with the exception of the TS at the third position (fig. 2 ). However, the graph shows that intra-Muridae pairwise comparisons also appear to be moderately subjected to saturation (slope of the regression analysis S = 0.50). Thus, our major interest being in intra-Muridae relationships, all substitutions were retained for phylogenetic reconstructions derived from the vWF gene.



View larger version (22K):
[in this window]
[in a new window]
 
Fig. 2.—Graphic estimation of saturation in vWF for transitions at third codon positions. The inferred number of changes, derived from the patristic distance matrix taken from a most-parsimonious tree, are plotted against the pairwise numbers of observed differences. Equations of the linear regression (straight lines) and correlation coefficients are given for Dipodidae-Muridae and intra-Muridae pairwise comparisons

 
Phylogenetic Reconstructions
According to the results obtained with the LCAT nuclear gene (Michaux and Catzeflis 2000Citation ) and other nuclear genes (GHR [Adkins et al. 2001Citation ], IRBP [DeBry and Sagel 2001Citation ], A2AB [D. Huchon, personal communication]) which showed a high robustness for a sister group relationship between Dipodidae and Muridae and for the monophyly of each superfamily, we used four species of Dipodidae as an outgroup for the Muridae.

LCAT Gene
The first set of analyses considered the 32 rodent species sequenced for the LCAT gene. The MP reconstruction yielded one most-parsimonious tree (L = 1,080 steps; consistency index [CI] = 0.41 and retention index [RI] = 0.52). Bootstrapping and BSI values are indicated in table 2 for the ancestral segments labeled 1–30 in figure 3 . The majority of the results of this analysis and those of the ML and distance approaches have been presented in Michaux and Catzeflis (2000)Citation . To summarize, the first dichotomy isolates the Spalacinae and the Rhizomyinae. The remaining 13 Muridae subfamilies are clustered together in a strongly supported clade. The first branching event of this clade is a large polytomy leading to eight lineages, of which only four encompass more than one subfamily: Mystromyinae and Nesomyinae (node 10), Cricetomyinae and Dendromurinae (node 12), Cricetinae and Myospalacinae (node 19), and Murinae, Otomyinae, Gerbillinae, and Acomyinae (node 21). Thus, with this gene, most of the "advanced" murid subfamilies appear to be of a polytomous origin, suggesting the phenomenon of a spectacular bushlike radiation having led to the majority of them.


View this table:
[in this window]
[in a new window]
 
Table 2 Indices of Robustness for the Nodes (labeled 1–30) of the Phylogenetic Trees Represented in Figures 4 and 5 Using Maximum-Parsimony (MP), Distance (NJ), Maximum-Likelihood (ML), and Markov Evolution Analyses on the Three Codon Positions (ME123) (only for the combined analysis)

 


View larger version (44K):
[in this window]
[in a new window]
 
Fig. 3.—Synthetic trees summarizing the results derived from four approaches on 32 mammalian DNA sequences of the LCAT and vWF genes. The robustness of each node (labeled 1–30) is described in table 2 for maximum-parsimony (bootstrap percentage [BP] and Bremer's support index), distance (BP), maximum-likelihood (BP), and Markov evolution analyses (BP). The trees were rooted by the four Dipodidae sequences

 
Exon 28 of vWF
The second set of analyses considered the 32 rodent species sequenced for vWF. The MP reconstruction yielded one most-parsimonious tree (L = 2,043 steps; CI = 0.42 and RI = 0.51). Bootstrapping and BSI values are indicated in table 2 for the ancestral segments labeled 1–30 in figure 3 . As already observed with LCAT, the first dichotomy isolates the Spalacinae and the Rhizomyinae. The remaining Muridae built up a strongly supported clade (node 4: BP = 100%, °BSI = +17). However, in contrast to the LCAT gene results, these "modern" subfamilies do not appear to be of a polytomous origin. Several clades appear clearly, i.e., an "African" clade (node 6) clustering the Nesomyinae, Mystromyinae, Cricetomyinae, and Dendromurinae (as defined in Michaux and Catzeflis 2000Citation ) or a "Cricetoid" clade (node 15) uniting the Arvicolinae, the Cricetinae, the Sigmodontinae, and the Myospalacinae.

However, some terminal nodes (i.e., the monophyly of the Cricetomyinae, the position of Otomys within the Murinae) appear less robust with the vWF gene than with LCAT. Thus, as these two genes seem complementary, we decided to combine the two data matrices in order to perform a new combined analysis.

Combined Analysis of the LCAT and vWF Genes
The ILD test revealed that there was some incongruence between the two nuclear genes (P = 0.002 < 0.05). This result might be explained by the extensive heterogeneity of the two data matrices (large differences in GC contents, some species are evolving very quickly or slowly with regard to the other taxa, etc.) and by the fact that these genes seem to evolve at different rates (see above). According to our results, LCAT yields a better (i.e., more robust) resolution for terminal nodes, whereas vWF seems to perform better for the deeper ancestral segments.

Thus, these genes, although they are suggested to be incongruent by the ILD test, also appear rather complementary. For this reason, the two data sets were concatenated.

The MP analysis yielded one most-parsimonious tree (L = 3,145 steps; CI = 0.41 and RI = 0.51) identical to the ML and NJ and ME trees (fig. 4 and table 2 ). Once again, the first dichotomy isolates Spalacinae and Rhizomyinae, and the remaining Muridae built up a strongly supported clade (node 4: BP = 100%, °BSI = +24).



View larger version (49K):
[in this window]
[in a new window]
 
Fig. 4.—Synthetic tree summarizing the results from the combined analysis using LCAT and vWF genes. The robustness of each node (labeled 1–30) is described in table 2 for maximum-parsimony (bootstrap percentage [BP] and Bremer's support index), distance (BP), maximum-likelihood (reliability percentages), and Markov evolution analyses (BP). The tree was rooted by four Dipodidae sequences

 
The monophyly of almost all subfamilies represented by at least two genera is robust (nodes 7, 9, 11, 16, 17, 19, 23, 24, 28 supported by BP of 97%–100% and by BSI of +10–+40): Nesomyinae (node 7), Dendromurinae (as previously defined in Michaux and Catzeflis 2000Citation ) (node 9), Cricetomyinae (node 11), Arvicolinae (node 16), Sigmodontinae (node 17), Cricetinae (node 19), Gerbillinae (node 23), Acomyinae (as previously defined in Michaux and Catzeflis 2000Citation ) (node 24), and Murinae (node 28). These different subfamilies are divided into four clades. The first (node 6), highly supported (BP = 97%, BSI = +9), includes the Dendromurinae, Mystromys, the Cricetomyinae, and the Nesomyinae. Within it, the first three subfamilies seem related, although with a lower robustness (BP = 75%, BSI = +3) (node 8). This result is probably linked to a conflict between the two genes. Indeed, LCAT associates Mystromys with the Nesomyinae (node 10, fig. 3 ). The second major lineage includes Calomyscus only. The third clade associates the Arvicolinae, the Cricetinae, the Sigmodontinae, and the Myospalacinae (node 15: BP = 89%, BSI = +5). Finally, the last group includes the Murinae, the Gerbillinae, and the Acomyinae (BP = 98%, BSI = +13) (node 21). Within it, the ancestral fragment uniting the Gerbillinae with Acomyinae is strongly supported (BP = 96%, BSI = +8) (node 22). A relationship between the third and the fourth lineages is also observed (node 14: BP = 86%, BSI = +5).

The NJ (GTR and LogDet), ML, and ME analyses yielded results similar to those obtained with the parsimony approach. However, some ancestral segments merit further comment, as they appear to be more robustly supported by some analyses—the "Holarctic" group uniting the third and fourth clades (node 14) (BP for ML = 97%); sister group relationships between the "Holarctic" group and Calomyscus (node 13) (BP for ML = 91%) and between Sigmodontinae and Cricetinae (node 18) (BP for ME = 94%); within Acomyinae, an external position of Lophuromys (node 25) (BP for NJ and ME = 87% and 88%, respectively); and the node 26, uniting Uranomys and Acomys (BP for NJ and ME = 80% and 86%, respectively).

Likelihood Alternatives to the Best Tree
The highest-likelihood tree (ln L = -18,177.01) for 32 Myodonta species was identified with PUZZLE (Strimmer and von Haeseler 1996Citation ) among 945 alternative trees constructed using MOLPHY 2.3b3 (Adachi and Hasegawa 1996Citation ). This tree has the same topology previously obtained by the parsimony and distance criteria and was used as the reference topology to apply Kishino-Hasegawa tests (Kishino and Hasegawa 1989Citation ) for assessing the following clades: the existence of the five Muridae lineages (Spalacinae + Rhizomyinae and the four modern groups); the existence of a "Holarctic" group; an isolated position of Calomyscus; a closer relationship between Calomyscus and the "Holarctic" subfamilies than with the "African" clade (Nesomyinae, Dendromurinae, Cricetomyinae, Mystromyinae); the sister group relationships between (a) Mystromys, Cricetomyinae, and Dendromurinae, (b) Gerbillinae + Acomyinae, and (c) Sigmodontinae and Cricetinae; within the Acomyinae, (a) an early offshoot of Lophuromys and (b) a closer relationship between Uranomys and Acomys; the monophyly of (a) the Sigmodontinae, (b) the Nesomyinae, and (c) the Acomyinae; and the nested position of (a) Otomyinae within Murinae and (b) Myospalacinae within Cricetinae.

For doing so, we tested 19 different alternative topologies derived from previous morphological or molecular studies (Thomas 1896Citation ; Miller and Gidley 1918Citation ; Misonne 1971Citation ; Chaline, Mein, and Petter 1977Citation ; de Graaff 1981Citation ; Ameur 1984Citation ; Carleton 1984Citation ; Flynn, Jacobs, and Lindsay 1985Citation ; Dickerman 1992Citation ; Musser and Carleton 1993Citation ; Dubois, Catzeflis, and Beintema 1999Citation ; Michaux and Catzeflis 2000Citation ; Chevret, Catzeflis, and Michaux 2001Citation ).

Almost all the alternative topologies exhibited significantly worse log likelihoods at 5% probability. However, the alternative hypotheses against a Holarctic group and the sister group relationships between Sigmodontinae and Cricetinae were only significantly worse at 8% probability. Moreover, the topologies uniting Mystromys to the Nesomyinae and Acomys to Deomys and proposing Uranomys as an early offshoot of the Acomyinae were not significantly worse than the best tree.

Relative-Rate Tests
To determine if differences of rates of vWF change existed between the different Muridae subfamilies, relative-rate tests were conducted with each of these subfamilies against the remaining lineages. The Dipodidae were used as the outgroup. As previously, synonymous (Ks) changes did not show significant differences between the different Muridae subfamilies. However, Ka comparisons showed that Spalacinae, Nesomyinae, Calomyscinae (Calomyscus), and Sigmodontinae were slowly evolving (P < 0.00001, 0.01, 0.04, and 0.002, respectively) and that Cricetomyinae and Rhizomyinae were quickly evolving (P < 0.0001 and 0.00001). Within Malagasy rodents and the New World Sigmodontinae, further analyses showed that only Macrotarsomys and Neotoma, respectively, had a lower rate of evolution (P < 0.02). Within the African Cricetomyinae, only Saccostomus was significantly more quickly evolving (P < 0.0001).

As the RRTree test is rather sensitive to the rate of the species used as outgroups (Robinson et al. 1998Citation ), we also applied the likelihood approach of Muse and Gaut (1994)Citation . This confirmed the results derived from the above RRTree tests based on Ka values.

The different behaviors obtained between Ks and Ka can be explained by the fact that, for the muroids, synonymous substitutions saturate when the phylogenetic level is higher than within the subfamily, as already suggested in Michaux and Catzeflis (2000)Citation .

Consequently, in order to apply a molecular clock and date the separation between the murid genera and subfamilies, we performed another ML analysis with Dipodidae as the outgroup and all the Muridae except the slowest- and fastest-evolving species (Spalax, Calomyscus, Rhizomys, Macrotarsomys, Neotoma, and Saccostomus). The inferred ML distances were the basis for estimating separation times. Two calibration points derived from paleontological data were chosen: (1) the Mus/Rattus dichotomy, set at 12 MYA (Jaeger, Tong, and Buffetaut 1986Citation ; Jacobs, Winkler, and Murry 1989Citation ; Jacobs et al. 1990Citation ; Jacobs and Downs 1994Citation ; Muse and Gaut 1994Citation ) and (2) The Gerbillus/Tatera separation, set at 8 MYA (Tong 1989Citation ). The estimated time of divergence (41 MYA) between Mus and Rattus obtained from molecular data (Kumar and Hedges 1998Citation ) was not used in this study because recent publications (Adkins et al. 2001Citation ; Murphy et al. 2001Citation ) have confirmed that it is an overestimate. The ML distance between Mus and Rattus is 0.049, whereas that between Tatera and Gerbillus is 0.028. These values give a rate of 0.0041 (Mus/Rattus) or 0.0035 ML distance per million years. When these rates are applied to the different lineages of Muridae, molecular datings can be suggested for several dichotomies of interest (table 3 ), such as 17.9–20.8 MYA for the separation between Murinae, Gerbillinae, and Acomyinae or 11.4–13.3 MYA between the two dendromurines Steatomys and Dendromus.


View this table:
[in this window]
[in a new window]
 
Table 3 Estimates of the Separation Times of Different Events Within Muridae on the Basis of the Molecular Data

 

    Discussion
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions
 Acknowledgements
 References
 
Base Composition
The GC3 compositions of vWF and LCAT (64.1%–89.1% and 54.8%–67.6%, respectively) indicate that these genes belong to the richest class of isochores, as already observed in Huchon, Catzeflis, and Douzery (1999)Citation and Robinson, Gautier, and Mouchiroud (1997)Citation . Also, other nuclear genes sequenced in different rodent species show similar GC3 contents, such as the nuclear ribonuclease (67.7%–82.3%; Dubois, Catzeflis, and Bientema 1999Citation ) or the p53 gene (54.3%–61.3%; D'Erchia et al. 1999Citation ). Accordingly, most of the genes would be distributed in the GC-rich fractions of DNA, which seems to be a general trend in mammals. Indeed, for humans, it has been shown that about 90% of genes are distributed in the two most GC-rich isochores, H2 and H3 (Zoubak, Clay, and Bernardi 1996Citation ).

Nevertheless, large differences between Muridae (vWF: 68.6% ± 2.9%; LCAT: 59.7% ± 3.1%), Dipodidae (vWF: 81.7% ± 8.3%; LCAT: 64.4% ± 3.1%), and human (vWF: 82%; LCAT: 77%) base compositions have been observed, being more pronounced for GC3. These results are in agreement with previous surveys of both LCAT (Robinson, Gautier, and Mouchiroud 1997Citation ) and vWF (Huchon, Catzeflis, and Douzery 1999Citation ). The observed pattern for murid LCAT and vWF has proved to follow a general trend described for the whole genome of murids and known as the murid pattern (Sabeur et al. 1993Citation ; Robinson, Gautier, and Mouchiroud 1997Citation ). Accordingly, murids would be characterized by a shift in base composition toward GC-poorer DNA when compared with the pattern observed in nonmurid rodents and all other mammals (general pattern). For Dipodidae, we found that GC3 values were between human and Muridae values. Similar results were obtained from CsCl profiles, in which Dipodidae showed a skewness value higher than those of murids but slightly lower than those of nonmurid rodents (Douady et al. 2000Citation ).

Even though these two genes are characterized by large GC3 values, the LCAT gene shows lower values than vWF, regardless of the species, Muridae, Dipodidae, or humans. For humans, this difference may be due to the fact that these two genes belong to different GC-rich isochores, namely, H2 and H3 (Zoubak, Clay, and Bernardi 1996Citation ). The same thing can be postulated for Muridae and Diopodidae, but taking into account the general shift toward GC-poorer DNA.

From the taxonomic point of view, no apparent relationship between taxonomy (or systematics) and GC3 content is readily observable for either LCAT or vWF genes within Muridae (data available on request). Average GC3 percentages for selected taxa are as follows for vWF and LCAT, respectively: Acomyinae (69.3% and 60.0%), Arvicolinae (68.7% and 62.7%), Cricetinae (66.0% and 56.0%), Cricetomyinae (69.5% and 62.5%), Dendromurinae (69.2% and 60.5%), Gerbillinae (71.0% and 61.6%), Murinae (69.1% and 60.2%), Nesomyinae (69.9% and 60.5%), and Sigmodontinae (66.2% and 59.8%).

Rate of Evolution
As we have stressed for base composition, LCAT and vWF show different properties that can be also detected at the level of rate of evolution and, subsequently, at the level of saturation. Indeed, LCAT is more conserved at both the amino acid and the nucleotide levels than is vWF and presents a lower asynonymous substitution rate value than the latter (fig. 1 ). This variation is thought to reflect differences in functional constraints in both genes, mainly in the proportion of the sequences that are critical to the function of the proteins. Nucleotide substitution rate, and thus the degree of sequence conservation, may be influenced by different factors affecting either mutation or fixation rates:

Regarding synonymous rates of substitution, LCAT and vWF genes show similar substitution rates, as shown in figure 1 . Since it seems there is no significant variation of synonymous rates between tissue-specific genes (Duret and Mouchiroud 2000Citation ), the observed difference between both genes may be due to differences in mutation rates. Nevertheless, the observed difference is most likely an underestimation of the real values because of a certain level of saturation in the transitions of both the vWF (fig. 2 ) and the LCAT genes (Michaux and Catzeflis 2000Citation ).

Because of the different evolutionary properties of vWF and LCAT, their combination for resolving phylogenetic relationships among muroid rodents has proven to give results with more resolution than those obtained using each single gene (Huchon, Catzeflis, and Douzery 1999Citation ; Michaux and Catzeflis 2000Citation )

Fourteen Subfamilies of Muroids Are Organized into Five Major Lineages
An examination of figure 4 , with the associated robustness values indicated in table 2 , indicates that the combined LCAT and vWF genes suggest a muroid branching pattern with five clades comprising one or more traditional subfamilies.

The most ancient speciation event leads to ancestral segments 4 and 5, a dichotomy separating spalacines and rhizomyines from the rest of the muroids. Both clades are robust (bootstrap values from 91 to 100, Bremer's decay indices of 24 and 10). Moreover, the alternative topologies against this first isolation were always significantly worse. Following other authors (Thomas 1896Citation ; Sen 1977Citation ), we are tempted to consider spalacines and rhizomyines as the two living components of a muroid family Spalacidae, keeping the taxon name Muridae (Illiger 1811Citation ) for all the remaining muroid taxa.

Thus, the ancestral segment 5 would define Muridae, for which our data set includes representatives of 12 subfamilies. The reliability of our Muridae/Spalacidae split will be settled after the inclusion of the three missing subfamilies listed by Musser and Carleton (1993)Citation : Lophiomyinae (one genus), Petromyscinae (two genera), and Platacanthomyinae (two genera).

Among the four major lineages of murids, a first clade (ancestral segment 6, bootstrap values from 95% to 99%) unites African taxa representing nesomyines, dendromurines, cricetomyines, and the sole living member of mystromyines. This cluster of genera is reminiscent of Lavocat's (1973, 1978)Citation concept of Nesomyidae, a taxon comprising archaic African "cricetids" which could all be derived from the fossil Afrocricetodontinae subfamily. Until additional taxa are examined for their molecular relationships (Malagasy Eliurus and Gymnuromys, African Petromyscus and Beamys, additional dendromurines), we refrain from suggesting a name for this clade.

The second lineage in murids includes Calomyscus only (see below). The third (ancestral segment 15) and fourth (ancestral segment 21) major murid clades are related, although not convincingly so (bootstrap values from <50% to 97%; a Bremer's decay index of 5; alternative topologies against this relationship were only significantly worse at 8% probability). All taxa representing the cricetines, myospalacines, sigmodontines, and arvicolines are united by a strong ancestral segment (15) (BP 89%–99%).

Finally, four "traditional" subfamilies (Gerbillinae, Murinae, Otomyinae, and Acomyines) cluster in a robust clade (ancestral segment 21: bootstrap values from 94% to 100%; a Bremer's decay index of 13).

An Isolated and Ancient Muroid: Calomyscus
The mouselike hamster Calomyscus (six living species) was only recently given a subfamilial rank (Calomyscinae) (Vorontsov and Potapova 1979Citation ) with regard to sharp differences in several morphological characters. Although previous studies had suggested affinities of Calomyscus with different muroid taxa, sometimes with compelling evidence (similarities in the auditory ossicles led Pavlinov [1980]Citation to relate mouselike hamsters with sigmodontine Reithrodontomys), its relationships have remained obscure. Even chromosomal data were ambiguous (Matthey 1961Citation ), with autosomal chromosomes similar to those of Eurasian cricetines, and sexual chromosomes reminiscent of North American peromyscines. Following Fahlbusch (1969)Citation and other paleontologists, Carleton (1984)Citation and Carleton and Musser (1984)Citation concluded that "Calomyscus could be justifiably classified among the cricetodontines, a group hitherto supposed extinct."

Both LCAT and exon 28 of the vWF gene concur for an isolated position of Calomyscus, which might be more related (node 13: bootstrap values from 64% to 91%) to Holarctic cricetids (arvicolines, sigmodontines, cricetines) and advanced murids (gerbillines, murines, acomyines) than to the African muroids (nesomyines, cricetomyines, dendromurines). Moreover, this hypothesis was confirmed with the Kishino-Hazegawa test, which showed that alternative topologies against these relationships always exhibit significantly worse log likelihood values. In any case, Calomyscus belongs to the same radiation (nodes 5, 13, and 14 combined) which has led to the three major clades (ancestral segments 6, 15, and 21) of murids.

Acomyines: A Clade Sister to the Gerbils
Much progress with regard to the systematics of Acomys has been made since Sarich's (1985)Citation immunological study, which was the first to suggest that spiny mice were not related to true murines (rats and mice). Different morphological and molecular studies have since shown evidence of a cluster of related genera, comprising Acomys, Uranomys, Lophuromys, and Deomys, which are neither murines nor dendromurines (see references cited in Denys et al. 1992Citation ; Chevret et al. 1994Citation ; Hänni et al. 1995Citation ; Dubois, Catzeflis, and Beintema 1999Citation ). Hänni et al. (1995)Citation named "acomyines" the clade containing Acomys, Uranomys, and Lophuromys, to which Deomys was subsequently added (Denys et al. 1995Citation ; Verheyen, Colyn, and Verheyen 1996Citation ; Chevret, Catzeflis, and Michaux 2001Citation ).

Although convincing biochemical and molecular data exist for sustaining this cluster, the relationships of acomyines with regard to other living muroids have remained obscure or controversial. Whereas the DNA-DNA hybridization of Chevret et al. (1993)Citation suggested an acomyine–gerbilline relationship, more recent experiments with additional taxa have been equivocal (Chevret, Catzeflis, and Michaux 2001Citation ). Sequence data have provided contradictory and weakly supported results (Hänni et al. 1995Citation ; Dubois, Catzeflis, and Beintema 1999Citation ; Chevret, Catzeflis, and Michaux 2001Citation ), with acomyines being either external to a murine + gerbilline clade or belonging to a polytomy also comprising murines and cricetines.

The combined use of two nuclear genes (LCAT and exon 28 of vWF) provides a clear and strong picture: acomyines cluster with gerbillines (ancestral segment 22; bootstrap values 96%–100%; BSI = 8, alternative topologies always significantly worse), and these two lineages are sister to murines (segment 21: BP = 94%–100%; BSI = 13). That these relationships are not biased by saturation is indicated by the robustness of deeper ancestral segments such as 14, 13, and 5 of figure 4 . Although our taxonomic sampling is rather limited, each of the clades acomyines, gerbillines, and murines is defined by high values of robustness (BP = 100%; BSI = 22, 40, and 34, respectively). According to the molecular-clock analysis, the separation between the three subfamilies appeared 17.9–20.8 MYA. This approximation is slightly older than those obtained with the previous paleontological (16 MYA; Tong and Jaeger 1993Citation ) and molecular (15.5–17 MYA; Michaux and Catzeflis 2000Citation ) data. Thus, on the basis of these results, the first radiation event within the modern muroids should have appeared earlier than previously supposed. Moreover, a second event of separation seems to have occurred 16.5–18.5 MYA. This led to the separation between Acomyinae and Gerbillinae, but also to the radiation of the "Cricetid" group (Arvicolinae, Sigmodontinae, Cricetinae, Myospalacinae) and the "African" clade (Nesomyinae, Mystromyinae, Cricetomyinae, Dendromurinae) (see table 3 ). According to Aguilar et al. (1996)Citation and Aguilar, Escarguel, and Michaux (1999)Citation , this period (end of early Miocene) was characterized by changes in climate which favored the spread in Europe, Africa, and America of allochtonous rodent groups probably coming from Asia.

The South African Mystromys Clusters with African Murids
Carleton (1984)Citation and Carleton and Musser (1984)Citation summarize perfectly the uncertain classification of the African white-tailed hamster Mystromys, which has been either allied to Holarctic cricetines or placed in its own subfamily (Mystromyinae) (Vorontsov 1966Citation ). Lavocat (1973)Citation suggested that Miocene Cricetodontidae found in African (Ethiopia) deposits were ancestral to fossil and recent Mystromys and proposed a classification of African muroids into two families: the family Nesomyidae, including not only the Malagasy Nesomyinae and the Mystromyinae, but also Cricetomyinae, Lophiomyinae, Tachyoryctinae, and Otomyinae; and the family Muridae, comprising Murinae, Rhizomyinae, and Dendromurinae. Gerbillinae were not allocated by Lavocat (1973Citation , p. 237), and were left incertae sedis.

Our present knowledge allows us to exclude Otomyinae and Tachyoryctinae from the Nesomyidae of Lavocat (1973)Citation , whose emended taxon thus corresponds to the clade defined by ancestral segment 6 in figure 4 . Thus, Mystromys belongs to a group of archaic African muroids, together with nesomyines, cricetomyines, and some dendromurines. Its relationships among this cluster are not resolved in this study, although exon 28 of the vWF gene suggests the tritomy Mystromys-cricetomyines-dendromurines. Additional taxa (such as Petromyscus or Beamys, Lophiomys) and gene sequences will be necessary to confirm these observations.

A Nesting Position of Myospalacinae Within Cricetinae and of Otomyinae Within Murinae
The analyses performed on the basis of the vWF sequences confirm those obtained with the LCAT gene (Michaux and Catzeflis 2000Citation ) about a nesting position of Myospalacinae within Cricetinae, particularly with the Asiatic species Phodopus. Thus, following a previous hypothesis based on morphological characters (Simpson 1945Citation ), we propose invalidation of the Myospalacinae subfamily and consideration of the genus Myospalax as defining a tribe among the subfamily Cricetinae. In the same way, the Otomyinae must be definitively considered a tribe of the Murinae subfamily, as already suggested by several morphological studies (Tullberg 1899Citation ; Miller and Gidley 1918Citation ; Simpson 1945Citation ; Carleton 1984Citation ).

Differences Between Datings Estimated by Fossils and Molecules
Although showing a trend toward overvaluation, most of the separation times estimated on the basis of the molecular data are in good agreement with those obtained with the fossil records (table 3 ). However, some of them seem more unclear. As already observed in Michaux and Catzeflis (2000)Citation , the divergence time estimation based on our molecular data suggests that the separation between Myospalax and the other Cricetinae appeared earlier (Late Miocene–Early Pliocene: 5–6 MYA) than the Pleistocene. Further paleontological investigations on this genus should be very useful to test our hypotheses.

Moreover, the value of 8.0–9.3 (standard error = 0.6) MYA for the split between the Arvicolinae Clethrionomys and Dicrostonyx is much older than the estimate of 4–6.0 MYA suggested by Chaline and Graf (1988)Citation . Our estimate is similar to the one calculated by Michaux and Catzeflis (2000)Citation , obtained with the nuclear LCAT gene only, but is at odds with other molecular studies based on DNA/DNA hybridization (Catzeflis et al. 1987Citation ) or on the nuclear ribonuclease gene sequences (Dubois, Catzeflis, and Beintema 1999Citation ). The relative-rate tests for the combined molecular data showed that the two Arvicolinae taxa do not evolve at a particular rate of evolution with regard to the other Muridae. Thus, the molecular clock of the LCAT and vWF genes seems also valid for these taxa. New paleontological investigations are needed to confirm this older Dicrostonyx/Clethrionomys divergence.


    Conclusions
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions
 Acknowledgements
 References
 
The present molecular study was performed on two nuclear genes sequenced for representatives of 14 among 17 Muridae subfamilies and four of the seven Dipodidae subfamilies.

This taxonomic sampling led to evidence that from the base composition and rate-of-evolution points of view, muroids would be characterized by a shift in base composition toward GC-poorer DNA when compared with the pattern observed in nonmurid rodents and all other mammals (general pattern). The LCAT gene shows lower GC3 values than vWF, regardless of the taxa under consideration (Muridae, Dipodidae, or humans), and is more conserved at both the amino acid and the nucleotide levels than is vWF. No apparent relation between taxonomy (or systematics) and GC3 content is readily observable.

The sampling also led to evidence that the studied subfamilies of muroids are organized into five major lineages: a first isolation unites the Spalacinae and Rhizomyinae with a strong robustness; the second one includes the African taxa representing nesomyines, dendromurines, cricetomyines, and the sole living member of mystromyines; the third one comprises only the genus Calomyscus; the fourth one represents the cricetines, myospalacine, sigmodontines, and arvicolines; and the fifth one comprises four "traditional" subfamilies (Gerbillinae, Murinae, Otomyinae, and Acomyines).

Also evidenced by the present study were the monophyly of almost all studied subfamilies, including the Spalacinae, the Rhizomyinae, the Nesomyinae, the Dendromurinae (as understood by Michaux and Catzeflis 2000Citation ), the Cricetomyinae, the Arvicolinae, the Sigmodontinae, the Cricetinae, the Gerbillinae, the Acomyinae (as defined in Michaux and Catzeflis 2000Citation ), and the Murinae; an isolated position of Calomyscus with regard to the other modern subfamilies; a sister group relationship of the Acomyinae and the Gerbillinae; confirmation of a nested position of Myospalacinae within Cricetinae and Otomyinae within Murinae; and a close relationship between the South African Mystromys and the African murids (Nesomyinae, Cricetomyinae, Dendromurinae).

Future studies including members of the remaining subfamilies listed by Musser and Carleton (1993)Citation (Lophiomyinae, Petromyscinae, and Platacanthomyinae), as well as some important missing genera for several subfamilies considered in our study (i.e., Murinae, Arvicolinae, Sigmodontinae), will confirm definitively the reliability of these results.

From the biogeographical point of view, an examination of the recent distribution of members of those five major lineages (fig. 5 ) allows the following comments, assuming that all muroids originated, spread, and radiated from Asia with different degrees of success: the first three groups radiated in a limited number of species (23, 43, and one, respectively) and genera (five, 18, and one, respectively) (Wilson and Reeder 1993Citation ) localized in some regions, whereas the last two groups literally exploded in a wide variety of species (591 and 654) and genera (109 and 139) dispersed all over the world.



View larger version (48K):
[in this window]
[in a new window]
 
Fig. 5.—Approximate geographic distribution (from data in Wilson and Reeder 1993) of living species in each major lineage within Muridae. The numbers of species are those derived from Musser and Carleton (1993)Citation assuming that each subfamily is monophyletic

 


    Acknowledgements
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Conclusions
 Acknowledgements
 References
 
This research was funded by European Community Project ERB-FMRX-CT98-0221. We thank Dr. Michael Carleton for his very helpful comments. We thank all people who provided tissue samples of rodent taxa which they collected, often in remote countries: Marina Baskevitch, Gary Bronner, Giancarlo Contrafatto, D. Bellars, Jean-Marc Duplantier, Piotr Gambarian, Jean-Claude Gautun, Laurent Granjon, I. Hansmann, Rolf Anker Ims, Petros Lymberakis, Daniel Rakotondravony, Marco Salvioni, Majad Zadé, and Vitaly Volobouev. Thanks also go to Emmanuel Douzery, Dorothée Huchon, and Pascale Chevret for their helpful discussions about analyses of saturation, relative-rate tests, the phylogeny of Muridae, and the vWF sequencing (D. Huchon) of tissues T-662 and T-552. Special thanks go to Marie-Claire Beckers and Corinne Lautier for much advice on automatic sequencing and to Professor F. Gricorescu (Institut Universitaire de Recherche Clinique, Montpellier) for the use of his automatic sequencer. This is contribution ISEM 2001-66 of Institut des Sciences de l'Evolution de Montpellier (UMR 5554 Centre National de la Recherche Scientifique, Université Montpellier II).


    Footnotes
 
Manolo Gouy, Reviewing Editor

1 Keywords: muroids phylogeny LCAT vWF Back

2 Address for correspondence and reprints: Johan Michaux, Laboratoire de Paléontologie, Institut des Sciences de l'Evolution, UMR 5554 Centre National de la Recherche Scientifique, Université de Montpellier 2, 34095 Montpellier, France. johan{at}isem.univ-montp2.fr . Back


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

    Adachi J., M. Hasegawa, 1996 MOLPHY: programs for molecular phylogenetics. Version 2.3 Institute of Mathematics, Tokyo, Japan

    Adkins R. M., E. L. Gelke, D. Rowe, R. 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]

    Aguilar J. P., G. Clauzon, A. de Goer de Herve, H. Maluski, J. Michaux, J. L. Welcomme, 1996 The MN3 fossil mammal-bearing locality of Beaulieu (France): biochronology, radiometric dating, and lower age limit of the Early Neogene renewal of the mammalian fauna in Europe Newsl. Stratigr 34:171-191

    Aguilar J. P., G. Escarguel, J. Michaux, 1999 A succession of Miocene rodent assemblages from fissure fillings in southern France: palaeoenvironmental interpretation and comparison with Spain Palaeogeog. Palaeoclimatatol. Palaeoecol 145:215-230[ISI]

    Ameur R., 1984 Découverte de nouveaux rongeurs dans la formation miocène de Bou Hanifia (Algérie Occidentale) Géobios 17:165-175

    Bremer K., 1988 The limits of amino acid sequence data in angiosperm phylogenetic reconstruction Evolution 42:795-803[ISI]

    Carleton M. D., 1984 Introduction to rodents Pp. 255–265 in S. Anderson and J. K. Jones, eds. Orders and families of Recent mammals of the world. John Wiley & Sons, New York

    Carleton M. D., G. G. Musser, 1984 Muroid rodents Pp. 289–379 in S. Anderson and J. K. Jones, eds. Orders and families of Recent mammals of the world. John Wiley & Sons, New York

    Catzeflis F., 1991 Animal tissue collections for molecular genetics and systematics Trends Ecol. Evol 6:168[ISI]

    Catzeflis F. M., F. H. Sheldon, J. E. Ahlquist, C. G. Sibley, 1987 DNA-DNA hybridization evidence of the rapid rate of muroid rodent DNA evolution Mol. Biol. Evol 4:242-253[Abstract]

    Chaline J., J.-D. Graf, 1988 Phylogeny of the Arvicolidae (Rodentia): biochemical and paleontological evidence J. Mamm 69:22-33[ISI]

    Chaline J., P. Mein, F. Petter, 1977 Les grandes lignes d'une classification évolutive des Muroidea Mammalia 41:245-252[ISI]

    Chevret P., R. M. Catzeflis, J. Michaux, 2001 Acomyinae: new molecular evidences ORSTOM eds. (special review) (in press)

    Chevret P., C. Denys, J.-J. Jaeger, J. Michaux, F. M. Catzeflis, 1993 Molecular evidence that the spiny mouse (Acomys) is more closely related to gerbils (Gerbillinae) than to true mice (Murinae) Proc. Natl. Acad. Sci. USA 90:3433-3436[Abstract]

    Chevret P., L. Granjon, J.-M. Duplantier, C. Denys, F. M. Catzeflis, 1994 Molecular phylogeny of the Praomys complex (Rodentia, Murinae): a study based on DNA/DNA hybridization experiments Zool. J. Linn. Soc 112:425-442[ISI]

    Conroy G. C., M. Pickford, B. Senut, J. Van Couvering, P. Mein, 1992 Otavipithecus namibiensis, first Miocene hominoid from southern Africa Nature 356:144-148[ISI][Medline]

    Cunningham C. W., 1997 Can three incongruence tests predict when data should be combined? Mol. Biol. Evol 14:733-740[Abstract]

    de Graaff G., 1981 The rodents of Southern Africa Butterworths, Durban, South Africa

    DeBry R., R. Sagel, 2001 Phylogeny of Rodentia (Mammalia) inferred from the nuclear-encoded gene IRBP Mol. Phylogenet. Evol 19:209-301

    Denys C., J. Michaux, F. Catzeflis, S. Ducrocq, P. Chevret, 1995 Morphological and molecular data against the monophyly of Dendromurinae (Muridae: Rodentia) Bonn. Zool. Beitr 45:173-190

    Denys C., J. Michaux, F. Petter, J.-P. Aguilar, J.-J. Jaeger, 1992 Molar morphology as a clue to the phylogenetic relationship of Acomys to the Murinae Isr. J. Zool 38:253-262[ISI]

    D'Erchia A. M., G. Pesole, A. Tullo, C. Saccone, E. Sbisa, 1999 Guinea pig p53 mRNA: identification of new elements in coding and untranslated regions and their functional and evolutionary implications Genomics 58:50-64[ISI][Medline]

    Dickerman A. W., 1992 Molecular systematics of some New World muroid rodents PhD dissertation, University of Wisconsin, Madison

    Douady C., N. Carels, O. Clay, F. Catzeflis, G. Bernardi, 2000 Diversity and phylogenetic implications of CsCl profiles from rodent DNAs Mol. Phylogenet. Evol 17:219-230[ISI][Medline]

    Dubois J.-Y., F. M. Catzeflis, J. Beintema, 1999 The phylogenetic position of "Acomyinae" (Rodentia, Mammalia) as sister-group of a Murinae + Gerbillinae clade: evidence from the nuclear ribonuclease gene Mol. Phylogenet. Evol 13:181-192[ISI][Medline]

    Duret L., D. Mouchiroud, 2000 Determinants of substitution rates in mammalian genes: expression pattern affects selection but not mutation rate Mol. Biol. Evol 17:68-74[Abstract/Free Full Text]

    Fahlbusch V., 1969 Pliozäne Pleistozäne Cricetinae (Rodentia, Mammalia) aus Polen Acta Zool. Cracoviensia 14:99-138

    Farris J. S., M. K, A. G. Kluge, C. Bult, 1995 Testing significance of incongruence Cladistics 10:315-319[ISI]

    Flynn L. J., 1990 The natural history of rhizomyid rodents Pp. 155–183 in E. Nevo and O. A. Reig, eds. Evolution of subterranean mammals at the organismal and molecular levels. Alan R. Liss, New York

    Flynn L. J., L. L. Jacobs, E. H. Lindsay, 1985 Problems in muroid phylogeny: relationships to other rodents and origin of major groups Pp. 589–616 in W. P. Luckett and J.-L. Hartenberger, eds. Evolutionary relationships among rodents. Plenum Press, New York

    Hänni C., V. Laudet, V. Barriel, F. M. Catzeflis, 1995 Evolutionary relationships of Acomys and other murids (Rodentia, Mammalia) based on complete 12S rRNA mitochondrial gene sequences Isr. J. Zool 41:131-146[ISI]

    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 321:611-620[Medline]

    Hooper E. T., M. D. Musser, 1964 The glans penis in Neotropical cricetines (family Muridae) with comments on classification of muroid rodents Misc. Publ. Univ. Mich 123:1-57

    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]

    Hughes A. L., 1997 Rapid evolution of immunoglobulin superfamily C2 domains expressed in immune system cells Mol. Biol. Evol 14:1-5[Abstract]

    Hurst L. D., N. G. C. Smith, 1999 Do essential genes evolve slowly? Curr. Biol 9:747-750[ISI][Medline]

    Illiger C., 1811 Prodromus systematis mammalium et avium additis terminis zoographicis utriusque classis C. Salfeld, Berlin

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

    Jacobs L. L., L. J. Flynn, W. R. Downs, J. C. Barry, 1990 Quo vadis, Antemus? The Siwalik muroid record. Pp. 573–586 in E. H. Linday, V. Fahlbusch, and P. Mein, eds. European Neogene mammal chronology. Plenum Press, New York

    Jacobs L. L., D. A. Winkler, P. A. Murry, 1989 Modern mammal origins: evolutionary grades in the Early Cretaceous of North America Proc. Natl. Acad. Sci. USA 86:4992-4995[Abstract]

    Jaeger J.-J., H. Tong, E. Buffetaut, 1986 The age of MusRattus divergence: paleontological data compared with the molecular clock C. R. Acad. Sci. Paris 302:917-922

    Kishino H., M. Hasegawa, 1989 Examination 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]

    Kuma K., N. Iwabe, T. Miyata, 1995 Functional constraints against variations on molecules from the tissue level: slowly evolving brain-specific genes demonstrated by protein kinase and immunoglobulin supergene families Mol. Biol. Evol 12:123-130[Abstract]

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

    Lavocat R., 1973 Les rongeurs du Miocene d'Afrique Orientale Mem. Trav. EPHE 1:1-284

    ———. 1978 Rodentia and Lagomorpha Pp. 69–89 in V. J. Maglio and H. B. S. Cooke, eds. Evolution of African mammals. Harvard University Press, Cambridge, Mass

    Lockhaert P. J., M. A. Steel, M. D. Hendy, D. Denny, 1994 Recovering evolutionary trees under a more realistic model of sequence evolution Mol. Biol. Evol 11:605-612[Free Full Text]

    Matassi G., P. M. Sharp, C. Gautier, 1999 Chromosomal location effects on gene sequence evolution in mammals Curr. Biol 9:786-791[ISI][Medline]

    Matthey R., 1961 Cytologie comparée des Cricetinae paléarctiques et américains Les chromosomes de Calomyscus bailwardi Thomas et la position du genre Calomyscus dans la systématique des Cricetinae. Rev. Suisse Zool 68:41-61

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

    Michaux J., F. Catzeflis, 2000 The bushlike radiation of muroid rodents is exemplified by the molecular phylogeny of the LCAT nuclear gene Mol. Phylogenet. Evol 17:280-293[ISI][Medline]

    Miller G. S., J. W. Gidley, 1918 Synopsis of the super-generic groups of rodents J. Wash. Acad. Sci 8:431-448

    Misonne X., 1971 Order Rodentia. Part 6 Pp. 1–39 in J. Meester and H. W. Setzer, eds. The mammals of Africa: an identification manual. Smithsonian Institution Press, Washington, D.C

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

    Muse S. V., B. S. Gaut, 1994 A likelihood approach for comparing synonymous and nonsynonymous nucleotide substitution rates, with application to the chloroplast genome Mol. Biol. Evol 11:715-724[Abstract/Free Full Text]

    Musser G. G., M. D. Carleton, 1993 Family Muridae Pp. 501–755 in D. E. Wilson and D. M. Reeder, eds. Mammal species of the world. A taxonomic and geographic reference. Smithsonian Institution Press, Washington, D.C., and London

    Pavlinov I. Y., 1980 Taxonomic status of Calomyscus Thomas (Rodentia, Cricetidae) on the basis of structure of auditory ossicles Zool. Zh 59:312-316[ISI]

    Perry J., A. Ashworth, 1999 Evolutionary rate of a gene affected by chromosomal position Curr. Biol 9:987-989[ISI][Medline]

    Philippe H., 1993 MUST, a computer package for management utilities for sequences and trees Nucleic Acids Res 21:5264-5272[Abstract]

    Philippe H., E. Douzery, 1994 The pitfalls of molecular phylogeny based on four species as illustrated by the Cetacea/Artiodactyla relationships J. Mamm. Evol 2:133-152

    Robinson M., C. Gautier, D. Mouchiroud, 1997 Evolution of isochores in rodents Mol. Biol. Evol 14:823-828[Abstract]

    Robinson M., M. Gouy, C. Gautier, D. Mouchiroud, 1998 Sensibility of the relative-rate test to taxonomic sampling Mol. Biol. Evol 15:1091-1098[Abstract]

    Sabeur G., G. Macaya, F. Kadi, G. Bernardi, 1993 The isochore patterns of mammalian genomes and their phylogenetic implications J. Mol. Evol 37:93-108[ISI][Medline]

    Saccone C., C. Lanave, G. Pesole, G. Preparata, 1990 Influence of base composition on quantitative estimates of gene evolution Meth. Enzymol 183:570-583[ISI][Medline]

    Saitou N., M. Nei, 1987 The neighbor-joining method: a new method for reconstructing phylogenetic trees Mol. Biol. Evol 4:406-425[Abstract]

    Sambrook J., E. F. Fritsch, T. Maniatis, 1989 Molecular cloning: a laboratory manual Cold Spring Harbor Laboratory, Cold Spring Harbor, N.Y

    Sarich V. M., 1985 Rodent macromolecular systematics Pp. 423–452 in W. P. Luckett and J.-L. Hartenberger, eds. Evolutionary relationships among rodents. Plenum Press, New York

    Sen S., 1977 La faune de rongeurs Pliocènes de Calta (Ankara, Turquie) Bull. Mus. Natl. Hist. Nat. Sci. Terre 465:1-171

    Simpson G. G., 1945 The principles of classification and a classification of mammals Bull. Am. Mus. Nat. Hist 85:1-350[ISI]

    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]

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

    Tamura K., M. Nei, 1993 Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees Mol. Biol. Evol 10:512-529[Abstract]

    Thomas O., 1896 On the genera of rodents: an attempt to bring up to date the current arrangement of the order Proc. Zool. Soc. Lond 1896:1012-1028

    Thompson J. D., D. G. Higgins, T. J. Gibson, 1994 CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, positions-specific gap penalties and weight matrix choice Nucleic Acids Res 22:4673-4680[Abstract]

    Tong H., 1989 Origine et évolution des Gerbillidae (Mammalia, Rodentia) en Afrique du Nord Mem. Soc. Geol. Fr 155:120

    Tong H., J.-J. Jaeger, 1993 Muroid rodents from the Middle Miocene Fort Ternan locality (Kenya) and their contribution to the phylogeny of muroids Paleontol. Abt. A 229:51-73

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

    Verheyen E., M. Colyn, W. Verheyen, 1996 A mitochondrial cytochrome b phylogeny confirms the paraphyly of the Dendromurinae Alston, 1896 (Muridae, Rodentia) Mammalia 60:780-785[ISI]

    Vorontsov N., 1966 Taxonomic position and a survey of the hamsters of the genus Mystromys Wagn. (Mammalia, Glires) Zool. Zh 45:436-446

    Vorontsov N. N., E. G. Potapova, 1979 Taxonomy of the genus Calomyscus (Cricetidae). 2. Status of Calomyscus in the system of Cricetinae Zool. Zh 58:1391-1397

    Wilson D. E., D. M. Reeder, 1993 Mammals species of the world. A taxonomic and geographic reference Smithsonian Institution Press, Washington

    Wu C.-I., W. H. Li, 1985 Evidence for higher rates of nucleotide substitutions in rodents than in man Proc. Natl. Acad. Sci. USA 82:1741-1745[Abstract]

    Yang Z., 1996 Among-site rate variation and its impact on phylogenetic analyses Trends Ecol. Evol 11:367-372[ISI]

    Zoubak S., O. Clay, G. Bernardi, 1996 The gene distribution of the human genome Gene 174:95-102[ISI][Medline]

Accepted for publication July 2, 2001.