Rodent phylogeny and a timescale for the evolution of Glires: evidence from an extensive taxon sampling using three nuclear genes

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

Molecular evolution of the nuclear von Willebrand factor gene in mammals and the phylogeny of rodents.

Nucleotide sequences of exon 28 of the von Willebrand Factor (vWF) were analyzed for a representative sampling of rodent families and eutherian orders, with one marsupial sequence as outgroup. The aim of this study was to test if inclusion of an increased taxonomic diversity in molecular analyses would shed light on three uncertainties concerning rodent phylogeny: (1) relationships between rodent families, (2) Rodentia monophyly, and (3) the sister group relationship of rodents and lagomorphs. The results did not give evidence of any particular rodent pattern of molecular evolution relative to a general eutherian pattern. Base compositions and rates of evolution of vWF sequences of rodents were in the range of placental variation. The 10 rodent families studied here cluster in five clades: Hystricognathi, Sciuridae and Aplodontidae (Sciuroidea), Muridae, Dipodidae, and Gliridae. Among hystricognaths, the following conclusions are drawn: a single colonization event in South America by Caviomorpha, a paraphyly of Old World and New World porcupines, and an African origin for Old World porcupines. Despite a broader taxonomic sampling diversity, we did not obtain a robust answer to the question of Rodentia monophyly, but in the absence of any other alternative, we cannot reject the hypothesis of a single origin of rodents. Moreover, the phylogenetic position of Lagomorpha remains totally unsettled.     

Variance of molecular datings, evolution of rodents and the phylogenetic affinities between Ctenodactylidae and Hystricognathi

The von Willebrand factor (vWF) gene has been used to understand the origin and timing of Rodentia evolution in the context of placental phylogeny vWF exon 28 sequences of 15 rodent families and eight non-rodent eutherian clades are analysed with two different molecular dating methods (uniform clock on a linearized tree; quartet dating). Three main conclusions are drawn from the study of this nuclear exon. First, Ctenodactylidae (gundis) and Hystricognathi (e.g. porcupines, guinea-pigs, chinchillas) robustly cluster together in a newly recognized clade, named 'Ctenohystrica'. The Sciurognathi monophyly is subsequently rejected. Pedetidae (springhares) is an independent and early diverging rodent lineage, suggesting a convergent evolution of the multiserial enamel of rodent incisors. Second, molecular date estimates are here more influenced by accuracy and choice of the palaeontological temporal references used to calibrate the molecular clock than by either characters analysed (nucleotides versus amino acids) or species sampling. The caviomorph radiation at 31 million years (Myr) and the pig porpoise split at 63 Myr appear to be reciprocally compatible dates. Third, during the radiation of Rodentia, at least three lineages (Gliridae, Sciuroidea and Ctenohystrica) emerged close to the Cretaceous-Tertiary boundary, and their common ancestor separated from other placental orders in the Late Cretaceous.     

High-level phylogeny of early Tertiary rodents: dental evidence

Major crown‐groups of rodents were well established in the early Tertiary, and fossils provide an invaluable window into their evolutionary history. The main focus of this project was to perform a cladistic assessment of the dental evidence for early Tertiary rodent cladogenesis – the masticatory apparatus and teeth are the most frequently preserved anatomical features in the fossil record. We focused on groups that existed in a period corresponding to their early history, combining fossils belonging to extinct lineages and to stem‐groups leading to modern lineages. While the monophyly of some groups is not systematically explored, our results have important implications for high‐level rodent relationships and systematics. These results are consistent with those of recent molecular phylogenies and reliably congruent with the stratigraphic record, thus enhancing the pertinence of dental characters for phylogenetic inference. Our approach provides evidence of a fundamental dichotomy in early rodent history. Two major clades have been identified: (1) the earliest ‘ctenodactyloid’ (Ctenodactylidae, Chapattimyidae, Yuomyidae, Diatomyidae) and hystricognathous (Tsaganomyidae, Baluchimyinae, ‘phiomorphs’, ‘caviomorphs’) rodents, and (2) the earliest ‘ischyromyoid’ rodents with their closest relatives (Muroidea + Dipodoidea + Geomyoidea + Anomaluroidea + Castoroidea + Sciuravidae + Gliroidea, and Sciuroidea + Aplodontoidea + Theridomorpha). This topology has led us to endorse Ctenohystrica as the first clade and propose a new taxon, Ischyromyiformes, for the second. Although minimized in our working hypothesis, the homoplasy in dental characters remains significant. However, a number of homoplasic characters reveal structuring in their internal distribution, allowing us to discern evolutionary morphological patterns, notably the pentalophodonty of molars, zygomasseteric complex and incisor enamel microstructure. © 2004 The Linnean Society of London, Zoological Journal of the Linnean Society, 2004, 142 , 105–134.     

Understanding angiosperm diversification using small and large phylogenetic trees

How will the emerging possibility of inferring ultra-large phylogenies influence our ability to identify shifts in diversification rate? For several large angiosperm clades (Angiospermae, Monocotyledonae, Orchidaceae, Poaceae, Eudicotyledonae, Fabaceae, and Asteraceae), we explore this issue by contrasting two approaches: (1) using small backbone trees with an inferred number of extant species assigned to each terminal clade and (2) using a mega-phylogeny of 55473 seed plant species represented in GenBank. The mega-phylogeny approach assumes that the sample of species in GenBank is at least roughly proportional to the actual species diversity of different lineages, as appears to be the case for many major angiosperm lineages. Using both approaches, we found that diversification rate shifts are not directly associated with the major named clades examined here, with the sole exception of Fabaceae in the GenBank mega-phylogeny. These agreements are encouraging and may support a generality about angiosperm evolution: major shifts in diversification may not be directly associated with major named clades, but rather with clades that are nested not far within these groups. An alternative explanation is that there have been increased extinction rates in early-diverging lineages within these clades. Based on our mega-phylogeny, the shifts in diversification appear to be distributed quite evenly throughout the angiosperms. Mega-phylogenetic studies of diversification hold great promise for revealing new patterns, but we will need to focus more attention on properly specifying null expectation.     

ARE RATES OF DIVERSIFICATION IN SUBTERRANEAN SOUTH AMERICAN TUCO-TUCOS (GENUS CTENOMYS,RODENTIA: OCTODONTIDAE) UNUSUALLY HIGH?

Subterranean rodents have been used frequently as examples of explosive speciation in mammals. We tested for differential rates of diversification by using information from molecular phylogenies to focus primarily on tuco-tucos (Rodentia: Octodontidae), the most speciose lineage of subterranean rodents. Tuco-tucos were not significantly more diverse than their sister taxon (octodontines); however, a lineages-through-time analysis suggests an increase in diversification at the base of the tuco-tuco clade.     

Long-branch attraction phenomenon and the impact of among-site rate variation on rodent phylogeny

The phylogenetic relationships among major lineages of rodents is one of the issues most debated by both paleontologists and molecular biologists. In the present study, we have analyzed all complete mammalian mitochondrial genomes available in the databases, including five rodent species (rat, mouse, dormouse, squirrel and guinea-pig). Phylogenetic analyses were performed on H-strand amino acid sequences by means of maximum-likelihood and on H-strand protein-coding and ribosomal genes by means of distance methods. Also, log-likelihood ratio tests were applied to different tree topologies under the assumption of rodent monophyly, paraphyly or polyphyly. The analyses significantly rejected rodent monophyly and showed the existence of two differentiated clades, one containing non-murids (dormouse, squirrel and guinea-pig) and the other containing murids (rat and mouse). Long-branch attraction between murids and the outgroups could not be responsible for the existence of two different rodent clades, as no significant differences in evolutionary rate have been observed, except in the case of the squirrel, which shows a lower rate. The impact of among-site rate variation models on the phylogeny of rodents has been evaluated using the gamma distribution model. Results have shown that relationships among rodents remained unchanged, and the general topology of the tree was not affected, even though some branches were not properly resolved, most likely due to a lack of fit between estimated and real rate heterogeneity parameters.     

Diversity and phylogenetic implications of CsCl profiles from rodent DNAs

Buoyant density profiles of high-molecular-weight DNAs sedimented in CsCl gradients, i.e., compositional distributions of 50- to 100-kb genomic fragments, have revealed a clear difference between the murids so far studied and most other mammals, including other rodents. Sequence analyses have revealed other, related, compositional differences between murids and nonmurids. In the present study, we obtained CsCl profiles of 17 rodent species representing 13 families. The modal buoyant densities obtained for rodents span the full range of values observed in other eutherians. More remarkably, the skewness (asymmetry, mean - modal buoyant density) of the rodent profiles extends to values well below those of other eutherians. Scatterplots of these and related CsCl profile parameters show groups of rodent families that agree largely with established rodent taxonomy, in particular with the monophyly of the Geomyoidea superfamily and the position of the Dipodidae family within the Myomorpha. In contrast, while confirming and extending previously reported differences between the profiles of Myomorpha and those of other rodents, the CsCl data question a traditional hypothesis positing Gliridae within Myomorpha, as does the recently sequenced mitochondrial genome of dormouse. Analysis of CsCl profiles is presented here as a rapid, robust method for exploring rodent and other vertebrate systematics.     

Tempo and mode of mandibular shape and size evolution reveal mixed support for incumbency effects in two clades of island‐endemic rodents (Muridae: Murinae)*

Existing radiations in a spatially limited system such as an oceanic island may limit the ecological opportunity experienced by later colonists, resulting in lower macroevolutionary rates for secondary radiations. Additionally, potential colonists may be competitively excluded by these incumbent (resident) species, unless they are biologically distinct (biotic filtering). The extant phenotypic diversity of secondary colonists may thus be impacted by lower rates of phenotypic evolution, exclusion from certain phenotypes, and transitions to new morphotypes to escape competition from incumbent lineages. We used geometric morphometric methods to test whether the rates and patterns of mandibular evolution of the Luzon “old endemic” rodent clades, Phloeomyini and Chrotomyini, are consistent with these predictions. Each clade occupied nearly completely separate shape space and partially separate size space. We detected limited support for decelerating and clade‐specific evolutionary rates for both shape and size, with strong evidence for a shift in evolutionary mode within Chrotomyini. Our results suggest that decelerating phenotypic evolutionary rates are not a necessary result of incumbency interactions; rather, incumbency effects may be more likely to determine which clades can become established in the system. Nonincumbent clades that pass a biotic filter can potentially exhibit relatively unfettered evolution.     

Structure and evolution of opossum,guinea pig,and porcupine cytochrome b genes

Summary We have sequenced the mitochondrial cytochrome b gene from the guinea pig, the African porcupine, and a South American opossum. A phylogenetic analysis, which includes 22 eutherian and four other vertebrate cytochrome b sequences, indicates that the guinea pig and the porcupine constitute a natural clade (Hystricomorpha) that is not a sister group to the clade of mice and rats (Myomorpha). Therefore, the hypothesis that the Rodentia is paraphyletic receives additional support. The artiodactyls, the perissodactyls, and the cetaceans form a group that is separated from the primates and the rodents. The 26 sequences are used to study the structure/function relationships in cytochrome b, whose function is electron transport. Most of the amino acid residues involved in the two reaction centers are well conserved in evolution. The four histidines that are believed to ligate the two hemes are invariant among the 26 sequences, but their nearby residues are not well conserved in evolution. The eight transmembrane domains represent some of the most divergent regions in the cytochrome b sequence. The rate of nonsynonymous substitution is considerably faster in the human and elephant lineages than in other eutherian lineages; the faster rate might be due to coevolution between cytochrome b and cytochrome c. Offprint requests to: W.-H. Li     

Pleistocene effects on the European freshwater fish fauna: double origin of the cobitid genus Sabanejewia in the Danube basin (Osteichthyes: Cobitidae)

Biogeographical hypotheses of European freshwater fishes were inferred using phylogeographic analysis of the complete cytochrome b and ATP synthase 8 and 6 mitochondrial genes (1982bp). To test the relative importance of drainage origin versus Pleistocene glaciations in the origin of primary freshwater fishes in Europe, we reconstructed the phylogenetic relationships of the genus Sabanejewia which is distributed in European waters. The phylogenetic relationships recovered for the genus Sabanejewia (n=75) provide support for the monophyly of six main evolutionary mtDNA lineages: Sabanejewia larvata, Sabanejewia romanica, Sabanejewia aurata/Sabanejewia caucasica, Sabanejewia kubanica, Sabanejewia baltica, and the Danubian-Balkanian complex. The Caucasian-Caspian mtDNA lineages, S. kubanica, S. aurata/S. caucasica, and the Northern European S. baltica represents the sister group of the Danubian-Balkanian complex mtDNAclade, supporting a Caucasian-Northern European origin of most of mtDNA lineages of the Central European freshwater fish fauna. The mtDNA divergence observed between the Danubian Sabanejewia species is too dissimilar to support their contemporary origin. Rather, the mtDNA data suggest that the Danubian Sabanejewia lineages most likely have a double origin, indicating that the European Sabanejewia lineages have experienced different historical processes for the following reasons. First, the origin of the S. larvata and S. romanica mtDNA clades predates the origin of the Danubian-Balkanian complex, and our results showed that the completion of the Alps and the origin of the Danube drainage seem to have promoted the speciation of the earliest Sabanejewia clades in the Miocene. Second, small genetic distances and the geographical pattern found within the Danubian-Balkanian complex clade indicate that the lineages included in this clade spread recently across the Danube and Greek river drainages. The inclusion of the S. balcanica species within all mtDNA lineages suggests that cyclical cold periods during the Pleistocene glaciations have favoured its rapid expansion and genetic homogenisation across Central European and Greek waters.     

MORPHOLOGICAL RATES OF ANGIOSPERM SEED SIZE EVOLUTION

The evolution of seed size among angiosperms reflects their ecological diversification in a complex fitness landscape of life‐history strategies. The lineages that have evolved seeds beyond the upper and lower boundaries that defined nonflowering seed plants since the Paleozoic are more dispersed across the angiosperm phylogeny than would be expected under a neutral model of phenotypic evolution. Morphological rates of seed size evolution estimated for 40 clades based on 17,375 species ranged from 0.001 (Garryales) to 0.207 (Malvales). Comparative phylogenetic analysis indicated that morphological rates are not associated with the clade's seed size but are negatively correlated with the clade's position in the overall distribution of angiosperm seed sizes; clades with seed sizes closer to the angiosperm mean had significantly higher morphological rates than clades with extremely small or extremely large seeds. Likewise, per‐clade taxonomic diversification rates are not associated with the seed size of the clade but with where the clade falls within the angiosperm seed size distribution. These results suggest that evolutionary rates (morphological and taxonomic) are elevated in densely occupied regions of the seed morphospace relative to lineages whose ecophenotypic innovations have moved them toward the edges.     

Phylogeography of two cryptic species of African desert jerboas (Dipodidae: Jaculus)

The lesser Egyptian jerboa Jaculus jaculus is a desert dwelling rodent that inhabits a broad Arabian–Saharan arid zone. Recently, two distant sympatric lineages were described in North‐West Africa, based on morphometric and molecular data, which may correspond to two cryptic species. In the current study, phylogenetic relationships and phylogeographical structure among those lineages and geographical populations from North Africa and the Middle East were investigated. The phylogeographical patterns and genetic diversity of the cytochrome b gene (1110 bp) were addressed on 111 jerboas from 41 localities. We found that the variation in Africa is partitioned into two divergent mitochondrial clades (10.5% divergence relating to 1.65–4.92 Mya) that corresponds to the two cryptic species: J. jaculus and J. deserti. Diversifications within those cryptic species/clades were dated to 0.23–1.13 Mya, suggesting that the Middle Pleistocene climatic change and its environmental consequences affected the evolutionary history of African jerboas. The third distant clade detected, found in the Middle East region, most likely represents a distinct evolutionary unit, independent of the two African lineages. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ??, ??–??.     

Local molecular clocks in three nuclear genes: divergence times for rodents and other mammals and incompatibility among fossil calibrations

Reconstructing the chronology of mammalian evolution is a debated issue between molecule- and fossil-based inferences. A methodological limitation of molecules is the evolutionary rate variation among lineages, precluding the application of the global molecular clock. We considered 2422 first and second codon positions of the combined ADRA2B, IRBP, and vWF nuclear genes for a well-documented set of placentals including an extensive sampling of rodents. Using seven independent calibration points and a maximum-likelihood framework, we evaluated whether molecular and paleontological estimates of mammalian divergence dates may be reconciled by the local molecular clocks approach, allowing local constancy of substitution rates with variations at larger phylogenetic scales. To handle the difficulty of choosing among all possible rate assignments for various lineages, local molecular clocks were based on the results of branch-length and two-cluster tests. Extensive lineage-specific variation of evolutionary rates was detected, even among rodents. Cross-calibrations indicated some incompatibilities between divergence dates based on different paleontological references. To decrease the impact of a single calibration point, estimates derived from independent calibrations displaying only slight reciprocal incompatibility were averaged. The divergence dates inferred for the split between mice and rats (approximately 13-19 Myr) was younger than previously published molecular estimates. The most recent common ancestors of rodents, primates and rodents, boreoeutherians, and placentals were estimated to be, respectively, approximately 60, 70, 75, and 78 Myr old. Global clocks, local clocks, and quartet dating analyses suggested a Late Cretaceous origin of the crown placental clades followed by a Tertiary radiation of some placental orders like rodents.     

Interactions within and between clades shaped the diversification of terrestrial carnivores

A longstanding debate in evolutionary biology and paleontology is whether ecological interactions such as competition impose diversity dependence on speciation and extinction rates. Here, we analyze the fossil record of terrestrial mammalian carnivores in North America and Eurasia using a Bayesian framework to assess whether their diversity dynamics were affected by diversity dependence within and between families (12 in Eurasia, 10 in North America). We found eight instances of within‐clade diversity dependence suppressing speciation rates and detected between‐clade effects increasing extinction rates in six instances. Diversity dependence often involved lineages that migrated between continents and we found that speciation was more responsive to diversity changes within the clade, whereas extinction responded to diversity of taxa in other clades. The analysis of the fossil record of Carnivora suggests that interactions within and between clades are associated with different speciation and extinction regimes, opening room for a broader theory of diversity dependence.     

A major shift in diversification rate helps explain macroevolutionary patterns in primate species diversity

Primates represent one of the most species rich, wide ranging, and ecologically diverse clades of mammals. What major macroevolutionary factors have driven their diversification and contributed to the modern distribution of primate species remains widely debated. We employed phylogenetic comparative methods to examine the role of clade age and evolutionary rate heterogeneity in the modern distribution of species diversity of Primates. Primate diversification has accelerated since its origin, with decreased extinction leading to a shift to even higher evolutionary rates in the most species rich family (Cercopithecidae). Older primate clades tended to be more diverse, however a shift in evolutionary rate was necessary to adequately explain the imbalance in species diversity. Species richness was also poorly explained by geographic distribution, especially once clade age and evolutionary rate shifts were accounted for, and may relate instead to other ecological factors. The global distribution of primate species diversity appears to have been strongly impacted by heterogeneity in evolutionary rates.     

Species and shape diversification are inversely correlated among gobies and cardinalfishes (Teleostei: Gobiiformes)

Gobies and their relatives are significant components of nearshore marine, estuarine, and freshwater fish faunas in both tropical and temperate habitats worldwide. They are remarkable for their ability to adapt to and diversify in a wide range of environments. Among gobiiform clades, species diversities vary widely, ranging from two species in Kurtidae to more than 1,000 species in Gobiidae. There is also great variation in head and body shape and in environmental preferences (fresh, brackish, or marine habitats). In this study, I used a time-calibrated molecular phylogeny, coupled with morphometric and comparative analyses, to examine evolutionary rates of both speciation and morphological diversification among gobiiform lineages. Projection of the phylogeny onto a shape-derived morphospace shows that Gobioidei is morphometrically distinct from its sister taxon Apogonoidei, but that families within Gobioidei overlap in morphospace. Analysis of species diversification rates indicates that three rate shifts have occurred over the evolutionary history of Gobiiformes. Relative to the other lineages, Kurtidae has exhibited a slowdown in speciation, whereas both Apogonidae and Gobiidae?+?Gobionellidae have experienced an increase in diversification. Comparative analyses show that in Apogonidae and Gobiidae?+?Gobionellidae, increased speciation is correlated with diminished rates of morphological diversification, differently manifested in either clade and among the various sublineages. The elevation in speciation rates and diminishment in rates of morphological change in both Apogonidae and the clade Gobiidae?+?Gobionellidae are correlated with shifts to oceanic habitats from freshwater. This pattern is the complement to that seen across the global radiation of acanthomorph fishes in which a decrease in species diversification is associated with an increase in morphological disparity.     

Higher-level systematics of rodents (Mammalia,Rodentia): Evidence from the mitochondrial 12S rRNA gene

Phylogenetic relationships among major rodent superfamilies traditionally have been difficult to establish because of the apparent high level of convergence and parallelism seen among morphological characters and/or rapid differentiation of rodent groups in the Paleocene/Eocene. Nucleotide sequence data from the mitochondrial 12S rRNA gene were used to clarify phylogenetic relationships among the major groups of rodents as defined by Brandt (1855) and Tullberg (1899). Based on the approximately 800 bp analyzed for the 12S rRNA gene in 59 mammalian species, including 25 of the 32 extant rodent families, the major rodent groups that could be defined as monophyletic clades were the Hystricognathi, the Muroidea, and the Geomyoidea. In addition, support for superfamilial sister-group relationships was found for Aplodontoidea with Sciuroidea and Dipodoidea with Muroidea.     

Phylogeographic structure and regional history of Lemniscomys striatus (Rodentia: Muridae) in tropical Africa

Aim This study aims to elucidate the phylogeography of the murid rodent Lemniscomys striatus and to evaluate the relative roles of ecological change, habitat patchiness, rivers and geological barriers in structuring patterns of diversity. Location Sub‐Saharan Africa. Methods The extent of phylogeographic patterns and molecular genetic diversity (cytochrome b gene) were addressed in a survey of 128 individuals of L. striatus from 42 localities. Using maximum parsimony, maximum likelihood, Bayesian, network and genetic structure analyses, we inferred intraspecific relationships and tested hypotheses for historical patterns of gene flow within L. striatus. Results Our results identified four major geographical clades within L. striatus: a West African clade, a Benin‐Nigeria clade, a Central African clade, and an East African clade. Several subclades were identified within these four major clades. Restricted gene flow with isolation by distance was recorded, which is congruent with the low dispersal ability of such a small murid rodent. No clear signal of population expansion was detected within clades or subclades. Main conclusions The western rift system and the Volta and Niger rivers may have acted as long‐term extrinsic barriers to gene flow, resulting in the emergence of the four main clades of L. striatus. The observed pattern of mitochondrial variation observed within each clade probably results from late Pleistocene climatic and vegetation changes: during adverse conditions (forest expansion), L. striatus may have survived only in refugia, and then experienced range expansion under favourable conditions (savanna expansion).