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 |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
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 2000
) 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)
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)
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 |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
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 1990
), 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 1994
) and the ED editor (MUST package; Philippe 1993
). The program AFAS (MUST package; Philippe 1993
) 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)
and Hassanin, Lecointre, and Tillier (1998)
. Details for such an analysis are provided by Michaux and Catzeflis (2000)
.
Phylogenetic Reconstructions
The aligned sequences were treated with different approaches: the stationary Markov model (Saccone et al. 1990
) (also called general time reversible [GTR]) in PAUP 4.0b2 (Swofford 1998
), the Tamura and Nei (1993)
model, and the LogDet (Lockhaert et al. 1994
) 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 1996
) and the proportion of invariant sites (I) were estimated with the maximum-likelihood (ML) method assuming the GTRMarkov 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
(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)
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. 1995
). One thousand randomizations were performed on variable sites only (Cunningham 1997
).
Relative-Rate Test
Relative-rate tests were conducted both with RRTree, version 1.0 (Robinson et al. 1998
) (which improves the test of Wu and Li [1985]
), by taking into account the taxonomic representativity and its phylogenetic relationships, and with the likelihood approach of Muse and Gaut (1994)
. 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 |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
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 ).
|
GC3 Content
The percentages of GC at third codon positions (GC3) spanned a much larger scatter in vWF (64.189.1) than in LCAT (54.867.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 1997
), 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)
, 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)
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.
|
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 130 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)
. 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.
|
|
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).
|
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 analysesthe "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 1996
) among 945 alternative trees constructed using MOLPHY 2.3b3 (Adachi and Hasegawa 1996
). 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 1989
) 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 1896
; Miller and Gidley 1918
; Misonne 1971
; Chaline, Mein, and Petter 1977
; de Graaff 1981
; Ameur 1984
; Carleton 1984
; Flynn, Jacobs, and Lindsay 1985
; Dickerman 1992
; Musser and Carleton 1993
; Dubois, Catzeflis, and Beintema 1999
; Michaux and Catzeflis 2000
; Chevret, Catzeflis, and Michaux 2001
).
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. 1998
), we also applied the likelihood approach of Muse and Gaut (1994)
. 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)
.
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 1986
; Jacobs, Winkler, and Murry 1989
; Jacobs et al. 1990
; Jacobs and Downs 1994
; Muse and Gaut 1994
) and (2) The Gerbillus/Tatera separation, set at 8 MYA (Tong 1989
). The estimated time of divergence (41 MYA) between Mus and Rattus obtained from molecular data (Kumar and Hedges 1998
) was not used in this study because recent publications (Adkins et al. 2001
; Murphy et al. 2001
) 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.920.8 MYA for the separation between Murinae, Gerbillinae, and Acomyinae or 11.413.3 MYA between the two dendromurines Steatomys and Dendromus.
|
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
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 1997
) and vWF (Huchon, Catzeflis, and Douzery 1999
). 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. 1993
; Robinson, Gautier, and Mouchiroud 1997
). 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. 2000
).
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 1996
). 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 2000
), 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 2000
).
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 1999
; Michaux and Catzeflis 2000
)
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 1896
; Sen 1977
), we are tempted to consider spalacines and rhizomyines as the two living components of a muroid family Spalacidae, keeping the taxon name Muridae (Illiger 1811
) 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)
: 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)
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 1979
) 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]
to relate mouselike hamsters with sigmodontine Reithrodontomys), its relationships have remained obscure. Even chromosomal data were ambiguous (Matthey 1961
), with autosomal chromosomes similar to those of Eurasian cricetines, and sexual chromosomes reminiscent of North American peromyscines. Following Fahlbusch (1969)
and other paleontologists, Carleton (1984)
and Carleton and Musser (1984)
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)
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. 1992
; Chevret et al. 1994
; Hänni et al. 1995
; Dubois, Catzeflis, and Beintema 1999
). Hänni et al. (1995)
named "acomyines" the clade containing Acomys, Uranomys, and Lophuromys, to which Deomys was subsequently added (Denys et al. 1995
; Verheyen, Colyn, and Verheyen 1996
; Chevret, Catzeflis, and Michaux 2001
).
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)
suggested an acomyinegerbilline relationship, more recent experiments with additional taxa have been equivocal (Chevret, Catzeflis, and Michaux 2001
). Sequence data have provided contradictory and weakly supported results (Hänni et al. 1995
; Dubois, Catzeflis, and Beintema 1999
; Chevret, Catzeflis, and Michaux 2001
), 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.920.8 MYA. This approximation is slightly older than those obtained with the previous paleontological (16 MYA; Tong and Jaeger 1993
) and molecular (15.517 MYA; Michaux and Catzeflis 2000
) 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.518.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)
and Aguilar, Escarguel, and Michaux (1999)
, 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)
and Carleton and Musser (1984)
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 1966
). Lavocat (1973)
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 (1973
, p. 237), and were left incertae sedis.
Our present knowledge allows us to exclude Otomyinae and Tachyoryctinae from the Nesomyidae of Lavocat (1973)
, 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 2000
) about a nesting position of Myospalacinae within Cricetinae, particularly with the Asiatic species Phodopus. Thus, following a previous hypothesis based on morphological characters (Simpson 1945
), 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 1899
; Miller and Gidley 1918
; Simpson 1945
; Carleton 1984
).
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)
, the divergence time estimation based on our molecular data suggests that the separation between Myospalax and the other Cricetinae appeared earlier (Late MioceneEarly Pliocene: 56 MYA) than the Pleistocene. Further paleontological investigations on this genus should be very useful to test our hypotheses.
Moreover, the value of 8.09.3 (standard error = 0.6) MYA for the split between the Arvicolinae Clethrionomys and Dicrostonyx is much older than the estimate of 46.0 MYA suggested by Chaline and Graf (1988)
. Our estimate is similar to the one calculated by Michaux and Catzeflis (2000)
, obtained with the nuclear LCAT gene only, but is at odds with other molecular studies based on DNA/DNA hybridization (Catzeflis et al. 1987
) or on the nuclear ribonuclease gene sequences (Dubois, Catzeflis, and Beintema 1999
). 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 |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
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 2000
), the Cricetomyinae, the Arvicolinae, the Sigmodontinae, the Cricetinae, the Gerbillinae, the Acomyinae (as defined in Michaux and Catzeflis 2000
), 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)
(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 1993
) 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.
|
![]() |
Acknowledgements |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Footnotes |
---|
1 Keywords: muroids
phylogeny
LCAT
vWF
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
.
![]() |
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
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. 255265 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. 289379 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
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. 155183 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. 589616 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. 149156 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. 573586 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. 6989 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
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. 139 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
Musser G. G., M. D. Carleton, 1993 Family Muridae Pp. 501755 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. 423452 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
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]