Suprafamilial relationships among Rodentia and the phylogenetic effect of removing fast-evolving nucleotides in mitochondrial, exon and intron fragments

Evolutionary history and systematics of European blind mole rats (Rodentia: Spalacidae: Nannospalax): Multilocus phylogeny and species delimitation in a puzzling group

Species delimitation is a powerful approach to assist taxonomic decisions in challenging taxa where species boundaries are hard to establish. European taxa of the blind mole rats (genus Nannospalax) display small morphological differences and complex chromosomal evolution at a shallow evolutionary divergence level. Previous analyses led to the recognition of 25 'forms' in their distribution area. We provide a comprehensive framework to improve knowledge on the evolutionary history and revise the taxonomy of European blind mole rats based on samples from all but three of the 25 forms. We sequenced two nuclear-encoded genetic regions and the whole mitochondrial cytochrome b gene for phylogenetic tree reconstructions using concatenation and coalescence-based species-tree estimations. The phylogenetic analyses confirmed that Aegean N. insularis belongs to N. superspecies xanthodon, and that it represents the second known species of this superspecies in Europe. Mainland taxa reached Europe from Asia Minor in two colonisation events corresponding to two superspecies-level taxa: N. superspecies monticola (taxon established herewith) reached Europe c. 2.1 million years ago (Mya) and was followed by N. superspecies leucodon (re-defined herewith) c. 1.5 Mya. Species delimitation allowed the clarification of the taxonomic contents of the above superspecies. N. superspecies monticola contains three species geographically confined to the western periphery of the distribution of blind mole rats, whereas N. superspecies leucodon is more speciose with six species and several additional subspecies. The observed geographic pattern hints at a robust peripatric speciation process and rapid chromosomal evolution. The present treatment is thus regarded as the minimum taxonomic content of each lineage, which can be further refined based on other sources of information such as karyological traits, crossbreeding experiments, etc. The species delimitation models also allowed the recognition of a hitherto unnamed blind mole rat taxon from Albania, described here as a new subspecies.

Improved mammalian family phylogeny using gap-rare multiple sequence alignment: A timetree of extant placentals and marsupials

The timing of mammalian diversification in relation to the Cretaceous-Paleogene (KPg) mass extinction continues to be a subject of substantial debate. Previous studies have either focused on limited taxonomic samples with available whole-genome data or relied on short sequence alignments coupled with extensive species samples. In the present study, we improved an existing dataset from the landmark study of Meredith et al. (2011) by filling in missing fragments and further generated another dataset containing 120 taxa and 98 exonic markers. Using these two datasets, we then constructed phylogenies for extant mammalian families, providing improved resolution of many conflicting relationships. Moreover, the timetrees generated, which were calibrated using appropriate molecular clock models and multiple fossil records, indicated that the interordinal diversification of placental mammals initiated before the Late Cretaceous period. Additionally, intraordinal diversification of both extant placental and marsupial lineages accelerated after the KPg boundary, supporting the hypothesis that the availability of numerous vacant ecological niches subsequent to the mass extinction event facilitated rapid diversification. Thus, our results support a scenario of placental radiation characterized by both basal cladogenesis and active interordinal divergences spanning from the Late Cretaceous into the Paleogene.

Exploring objective feature sets in constructing the evolution relationship of animal genome sequences

BACKGROUND

Exploring evolution regularities of genome sequences and constructing more objective species evolution relationships at the genomic level are high-profile topics. Based on the evolution mechanism of genome sequences proposed in our previous research, we found that only the 8-mers containing CG or TA dinucleotides correlate directly with the evolution of genome sequences, and the relative frequency rather than the actual frequency of these 8-mers is more suitable to characterize the evolution of genome sequences.

RESULT

Therefore, two types of feature sets were obtained, they are the relative frequency sets of CG1 + CG2 8-mers and TA1 + TA2 8-mers. The evolution relationships of mammals and reptiles were constructed by the relative frequency set of CG1 + CG2 8-mers, and two types of evolution relationships of insects were constructed by the relative frequency sets of CG1 + CG2 8-mers and TA1 + TA2 8-mers respectively. Through comparison and analysis, we found that evolution relationships are consistent with the known conclusions. According to the evolution mechanism, we considered that the evolution relationship constructed by CG1 + CG2 8-mers reflects the evolution state of genome sequences in current time, and the evolution relationship constructed by TA1 + TA2 8-mers reflects the evolution state in the early stage.

CONCLUSION

Our study provides objective feature sets in constructing evolution relationships at the genomic level.

A rodent anchored hybrid enrichment probe set for a range of phylogenetic utility: From order to species

Rodents are the largest order of mammals and contain several model organisms important to scientific research in a variety of fields, yet no large set of genomic markers have been designed for this group to date, hindering evolutionary studies into relationships of the group as a whole. Here we present a genomic probe set designed and optimized for rodents with a protocol that is easy to replicate with little laboratory investment. This design utilizes an anchored hybrid enrichment approach specifically targeting rodents to generate longer loci with a higher substitution rate than existing vertebrate probes to provide utility at various taxonomic levels. Using a test set of rodents from all five suborders, we successfully obtained alignments for 416 of the 418 target loci with an average of 1379 bp per locus and a total alignment of more than half a million base pairs. This genomic data set performed well in all phylogenetic analyses, especially in recent phylogenetic splits, with ample parsimony-informative sites within genera and even within species, showing more than four times as many single nucleotide polymorphisms per locus than a recent vertebrate ultraconserved elements study. Additional support is provided in resolving deeper clades in Rodentia. By providing this probe design, we hope that more laboratories can easily generate data for answering questions in rodents from species delimitation to understanding relationships among families in rapid radiations.

Horizontal Cell Feedback to Cone Photoreceptors in Mammalian Retina: Novel Insights From the GABA-pH Hybrid Model

How neurons in the eye feed signals back to photoreceptors to optimize sensitivity to patterns of light appears to be mediated by one or more unconventional mechanisms. Via these mechanisms, horizontal cells control photoreceptor synaptic gain and enhance key aspects of temporal and spatial center-surround receptive field antagonism. After the transduction of light energy into an electrical signal in photoreceptors, the next key task in visual processing is the transmission of an optimized signal to the follower neurons in the retina. For this to happen, the release of the excitatory neurotransmitter glutamate from photoreceptors is carefully regulated via horizontal cell feedback, which acts as a thermostat to keep the synaptic transmission in an optimal range during changes to light patterns and intensities. Novel findings of a recently described model that casts a classical neurotransmitter system together with ion transport mechanisms to adjust the alkaline milieu outside the synapse are reviewed. This novel inter-neuronal messaging system carries feedback signals using two separate, but interwoven regulated systems. The complex interplay between these two signaling modalities, creating synaptic modulation-at-a-distance, has obscured it’s being defined. The foundations of our understanding of the feedback mechanism from horizontal cells to photoreceptors have been long established: Horizontal cells have broad receptive fields, suitable for providing surround inhibition, their membrane potential, a function of stimulus intensity and size, regulates inhibition of photoreceptor voltage-gated Ca²⁺ channels, and strong artificial pH buffering eliminates this action. This review compares and contrasts models of how these foundations are linked, focusing on a recent report in mammals that shows tonic horizontal cell release of GABA activating Cl⁻ and HCO₃⁻ permeable GABA autoreceptors. The membrane potential of horizontal cells provides the driving force for GABAR-mediated HCO₃⁻ efflux, alkalinizing the cleft when horizontal cells are hyperpolarized by light or adding to their depolarization in darkness and contributing to cleft acidification via NHE-mediated H⁺ efflux. This model challenges interpretations of earlier studies that were considered to rule out a role for GABA in feedback to cones.

The Dynamics, Causes, and Impacts of Mammalian Evolutionary Rates Revealed by the Analyses of Capybara Draft Genome Sequences

Capybara (Hydrochoerus hydrochaeri) is the largest species among the extant rodents. The draft genome of capybara was sequenced with the estimated genome size of 2.6 Gb. Although capybara is about 60 times larger than guinea pig, comparative analyses revealed that the neutral evolutionary rates of the two species were not substantially different. However, analyses of 39 mammalian genomes revealed very heterogeneous evolutionary rates. The highest evolutionary rate, 8.5 times higher than the human rate, was found in the Cricetidae-Muridae common ancestor after the divergence of Spalacidae. Muridae, the family with the highest number of species among mammals, emerged after the rate acceleration. Factors responsible for the evolutionary rate heterogeneity were investigated through correlations between the evolutionary rate and longevity, gestation length, litter frequency, litter size, body weight, generation interval, age at maturity, and taxonomic order. The regression analysis of these factors showed that the model with three factors (taxonomic order, generation interval, and litter size) had the highest predictive power (R2 = 0.74). These three factors determine the number of meiosis per unit time. We also conducted transcriptome analysis and found that the evolutionary rate dynamics affects the evolution of gene expression patterns.

A phylogenomic rodent tree reveals the repeated evolution of masseter architectures

Understanding the number of times a trait has evolved is a necessary foundation for comprehending its potential relationships with selective regimes, developmental constraints and evolutionary diversification. Rodents make up over 40% of extant mammalian species, and their ecological and evolutionary success has been partially attributed to the increase in biting efficiency that resulted from a forward shift of one or two portions of the masseter muscle from the zygomatic arch onto the rostrum. This forward shift has occurred in three discrete ways, but the number of times it has occurred has never been explicitly quantified. We estimated an ultrametric phylogeny, the first to include all rodent families, using thousands of ultraconserved elements. We examined support for evolutionary relationships among the five rodent suborders and then incorporated relevant fossils, fitted models of character evolution, and used stochastic character mapping to determine that a portion of the masseter muscle has moved forward onto the rostrum at least seven times (with one reversal) during the approximately 70 Myr history of rodents. Combined, the repeated evolution of this key innovation, its increasing prevalence through time, and the species diversity of clades with this character underscores the adaptive value of improved biting efficiency and the relative ease with which some advantageous traits arise.

Evolutionary history of two cryptic species of northern African jerboas

BACKGROUND

Climatic variation and geologic change both play significant roles in shaping species distributions, thus affecting their evolutionary history. In Sahara-Sahel, climatic oscillations shifted the desert extent during the Pliocene-Pleistocene interval, triggering the diversification of several species. Here, we investigated how these biogeographical and ecological events have shaped patterns of genetic diversity and divergence in African Jerboas, desert specialist rodents. We focused on two sister and cryptic species, Jaculus jaculus and J. hirtipes, where we (1) evaluated their genetic differentiation, (2) reconstructed their evolutionary and demographic history; (3) tested the level of gene flow between them, and (4) assessed their ecological niche divergence.

RESULTS

The analyses based on 231 individuals sampled throughout North Africa, 8 sequence fragments (one mitochondrial and seven single copy nuclear DNA, including two candidate genes for fur coloration: MC1R and Agouti), 6 microsatellite markers and ecological modelling revealed: (1) two distinct genetic lineages with overlapping distributions, in agreement with their classification as different species, J. jaculus and J. hirtipes, with (2) low levels of gene flow and strong species divergence, (3) high haplotypic diversity without evident geographic structure within species, and (4) a low level of large-scale ecological divergence between the two taxa, suggesting species micro-habitat specialization.

CONCLUSIONS

Overall, our results suggest a speciation event that occurred during the Pliocene-Pleistocene transition. The contemporary distribution of genetic variation suggests ongoing population expansions. Despite the largely overlapping distributions at a macrogeographic scale, our genetic results suggest that the two species remain reproductively isolated, as only negligible levels of gene flow were observed. The overlapping ecological preferences at a macro-geographic scale and the ecological divergence at the micro-habitat scale suggest that local adaptation may have played a crucial role in the speciation process of these species.

Emergence of hystricognathous rodents: Palaeogene fossil record, phylogeny, dental evolution and historical biogeography

Although phylogenetic trees imply Asia as the ancestral homeland of the Hystricognathi clade (Rodentia: Ctenohystrica), curiously the oldest known fossil occurrences of hystricognathous rodents are not from Asia, but from Africa and South America, where they appear suddenly in the fossil record of both landmasses by the Late Middle Eocene. Here we performed cladistic and Bayesian (standard and tip-dating analyses) assessments of the dental evidence documenting early ctenohystricans, including several Asian ‘ctenodactyloids’, virtually all Palaeogene Asian and African hystricognaths known thus far and two representatives of the earliest known South American hystricognaths. Our results provide a phylogenetic context of early hystricognaths (with implications on systematics) and suggest that some Eocene Asian ‘ctenodactyloids’ could be considered as stem hystricognaths and pre-hystricognaths, although they were not recognized as such originally. However, this view does not fill the gap of the Eocene Asian hystricognath record, as the proposed results imply many ghost lineages extending back to the Middle Eocene for several Asian and African taxa. They also imply a complex early historical biogeography of the group, involving multiple dispersal events from Asia to Africa (and possibly from Africa back to Asia) and then to South America sometime during the Middle Eocene. Based on these phylogenetic considerations, we discuss the emergence of hystricognathous rodents from a morpho-anatomical perspective by analysing the differentiation of their masticatory apparatus and chewing movements, notably through the evolution of their dental patterns.

Novel hybrid action of GABA mediates inhibitory feedback in the mammalian retina

The stream of visual information sent from photoreceptors to second-order bipolar cells is intercepted by laterally interacting horizontal cells that generate feedback to optimize and improve the efficiency of signal transmission. The mechanisms underlying the regulation of graded photoreceptor synaptic output in this nonspiking network have remained elusive. Here, we analyze with patch clamp recording the novel mechanisms by which horizontal cells control pH in the synaptic cleft to modulate photoreceptor neurotransmitter release. First, we show that mammalian horizontal cells respond to their own GABA release and that the results of this autaptic action affect cone voltage-gated Ca2+ channel (CaV channel) gating through changes in pH. As a proof-of-principle, we demonstrate that chemogenetic manipulation of horizontal cells with exogenous anion channel expression mimics GABA-mediated cone CaV channel inhibition. Activation of these GABA receptor anion channels can depolarize horizontal cells and increase cleft acidity via Na+/H+ exchanger (NHE) proton extrusion, which results in inhibition of cone CaV channels. This action is effectively counteracted when horizontal cells are sufficiently hyperpolarized by increased GABA receptor (GABAR)-mediated HCO3- efflux, alkalinizing the cleft and disinhibiting cone CaV channels. This demonstrates how hybrid actions of GABA operate in parallel to effect voltage-dependent pH changes, a novel mechanism for regulating synaptic output.

De novo characterization of placental transcriptome in the Eurasian beaver (Castor fiber L.)

Our pioneering data provide the first comprehensive view of placental transcriptome of the beaver during single and multiple gestation. RNA-Seq and a de novo approach allowed global pattern identification of C. fiber placental transcriptome. Non-redundant beaver transcriptome comprised 211,802,336 nt of placental transcripts, grouped into 128,459 contigs and clustered into 83,951 unigenes. An Ensembl database search revealed 14,487, 14,994, 15,004, 15,267 and 15,892 non-redundant homologs for Ictidomys tridecemlineatus, Rattus norvegicus, Mus musculus, Homo sapiens and Castor canadensis, respectively. Due to expression levels, the identified transcripts were divided into two sets: non-redundant and highly expressed (FPKM > 2 in at least three examined samples), analysed simultaneously. Among 17,009 highly expressed transcripts, 12,147 had BLASTx hits. GO annotations (175,882) were found for 4301 transcripts that were assigned to biological process (16,386), cellular component (9149) and molecular function (8338) categories; 666 unigenes were also classified into 122 KEGG pathways. Comprehensive analyses were performed for 411 and 3078 highly expressed transcripts annotated with a list of processes linked to ‘placenta’ (31 GO terms) or ‘embryo’ (324 GO terms), respectively. Among transcripts with entire CDS annotation, 281 (placenta) and 34 (embryo) alternative splicing events were identified. A total of 8499 putative SNVs (~ 6.2 SNV/transcript and 1.7 SNV/1 kb) were predicted with 0.1 minimum frequency and maximum variant quality (p value 10e⁻⁹). Our results provide a broad-based characterization of the global expression pattern of the beaver placental transcriptome. Enhancement of transcriptomic resources for C. fiber should improve understanding of crucial pathways relevant to proper placenta development and successful reproduction.

Oldest skeleton of a fossil flying squirrel casts new light on the phylogeny of the group

Flying squirrels are the only group of gliding mammals with a remarkable diversity and wide geographical range. However, their evolutionary story is not well known. Thus far, identification of extinct flying squirrels has been exclusively based on dental features, which, contrary to certain postcranial characters, are not unique to them. Therefore, fossils attributed to this clade may indeed belong to other squirrel groups. Here we report the oldest fossil skeleton of a flying squirrel (11.6 Ma) that displays the gliding-related diagnostic features shared by extant forms and allows for a recalibration of the divergence time between tree and flying squirrels. Our phylogenetic analyses combining morphological and molecular data generally support older dates than previous molecular estimates (~23 Ma), being congruent with the inclusion of some of the earliest fossils (~36 Ma) into this clade. They also show that flying squirrels experienced little morphological change for almost 12 million years.

Mammals can walk, hop, swim and fly; a few, like marsupial sugar gliders or colugos, can even glide. With 52 species scattered across the Northern hemisphere, flying squirrels are by far the most successful group that adopted this way of going airborne. To drift from tree to tree, these small animals pack their own ‘parachute’: a membrane draping between their lower limbs and the long cartilage rods that extend from their wrists. Tiny specialized wrist bones, which are unique to flying squirrels, help to support the cartilaginous extensions.

The origin of flying squirrels is a point of contention: while most genetic studies point towards the group splitting from tree squirrels about 23 million years ago, the oldest remains – mostly cheek teeth – suggest the animals were already soaring through forests 36 million years ago. However, recent studies show that the dental features used to distinguish between gliding and non-gliding squirrels may actually be shared by the two groups.

In 2002, the digging of a dump site in Barcelona unearthed a peculiar skeleton: first a tail and two thigh bones, big enough that the researchers thought it could be the fossil of a small primate. In fact, and much to the disappointment of paleoprimatologists, further excavating revealed that it was a rodent. As the specimen – nearly an entire skeleton – was being prepared, paleontologists insisted that all the ‘dirt’ attached to the bones had to be carefully screen-washed. From the mud emerged the minuscule specialized wrist bones: the primate-turned-rodent was in fact Miopetaurista neogrivensis, an extinct flying squirrel.

Here, Casanovas-Vilar et al. describe the 11.6 million years old fossil, the oldest ever found. The wrist bones reveal that the animal belongs to the group of flying squirrels that have large sizes. Evolutionary analyses that combined molecular and paleontological data demonstrated that flying squirrels evolved from tree squirrels as far back as 31 to 25 million years ago, and possibly even earlier. In addition, the results show that Miopetaurista is closely related to Petaurista, a modern group of giant flying squirrels. In fact, their skeletons are so similar that the large species that currently inhabit the tropical and subtropical forests of Asia could be considered living fossils.

Molecular and paleontological data are often at odds, but this fossil shows that they can be reconciled and combined to retrace history. Discovering older fossils, or even transitional forms, could help to retrace how flying squirrels took a leap from the rest of their evolutionary tree.

Fossils know it best: Using a new set of fossil calibrations to improve the temporal phylogenetic framework of murid rodents (Rodentia: Muridae)

Murid rodents (Rodentia: Muridae) represent the most diverse and abundant mammalian family. In this study, we provide a refined set of fossil calibrations which is used to reconstruct a dated phylogeny of the family using a multilocus dataset (six nuclear and nine mitochondrial gene fragments) encompassing 161 species representing 82 murid genera from four extant subfamilies (Deomyinae, Gerbillinae, Lophiomyinae and Murinae). In comparison with previous studies on murid or muroid rodents, our work stands out for the implementation of nine robust fossil constraints within the Muridae thanks to a thorough review of the fossil record. Before being assigned to specific nodes of the phylogeny, all potential fossil constraints were carefully assessed; they were also subjected to several cross-validation analyses. The resulting phylogeny is consistent with previous phylogenetic studies on murids, and recovers the monophyly of all sampled murid subfamilies and tribes. Based on nine controlled fossil calibrations, our inferred temporal timeframe indicates that the murid family likely originated in the course of the Early Miocene, 22.0-17.0 million years ago (Ma), and that most major lineages (i.e. tribes) started diversifying ca. 10 Ma. Historical biogeography analyses support the tropical origin for the family, with an initial internal split (vicariance event) between Afrotropical and Oriental (Indomalaya and Philippines) lineages. During the course of their diversification, the biogeographic pattern of murids is marked by several dispersal events toward the Australasian and the Palearctic regions. The Afrotropical region was also secondarily colonized at least three times from the Indomalaya, indicating that the latter region has acted as a major centre of diversification for the family.

Identification of Placental Aspartic Proteinase in the Eurasian Beaver (Castor fiber L.)

Aspartic proteinases (AP) form a multigenic group widely distributed in various organisms and includes pepsins (pep), cathepsins D and E, pregnancy associated glycoproteins (PAGs) as well as plant, fungal, and retroviral proteinases. This study describes the transcript identification and expression localization of the AP within the discoid placenta of the Castor fiber. We identified 1257 bp of the AP cDNA sequence, encoding 391 amino acids (aa) of the polypeptide precursor composed of 16 aa signal peptide, 46 aa pro-piece, and 329 aa of the mature protein. Within the AP precursor, one site of potential N-glycosylation (NPS^119–121) and two Asp residues (D) specific for the catalytic cleft of AP were identified (VLFDTGSSNLWV^91–102 and GIVDTGTSLLTV^277–288). The highest homology of the identified placental AP nucleotide and aa sequence was to mouse pepsinogen C (75.8% and 70.1%, respectively). Identified AP also shared high homology with other superfamily members: PAGs, cathepsins, and napsins. The AP identified in this study was named as pepsinogen/PAG-Like (pep/PAG-L). Diversified pep/PAG-L protein profiles with a dominant 58 kDa isoform were identified. Immune reactive signals of the pep/PAG-L were localized within the trophectodermal cells of the beaver placenta. This is the first report describing the placental AP (pep/PAG-L) in the C. fiber.

Mammalian Comparative Genomics Reveals Genetic and Epigenetic Features Associated with Genome Reshuffling in Rodentia

Understanding how mammalian genomes have been reshuffled through structural changes is fundamental to the dynamics of its composition, evolutionary relationships between species and, in the long run, speciation. In this work, we reveal the evolutionary genomic landscape in Rodentia, the most diverse and speciose mammalian order, by whole-genome comparisons of six rodent species and six representative outgroup mammalian species. The reconstruction of the evolutionary breakpoint regions across rodent phylogeny shows an increased rate of genome reshuffling that is approximately two orders of magnitude greater than in other mammalian species here considered. We identified novel lineage and clade-specific breakpoint regions within Rodentia and analyzed their gene content, recombination rates and their relationship with constitutive lamina genomic associated domains, DNase I hypersensitivity sites and chromatin modifications. We detected an accumulation of protein-coding genes in evolutionary breakpoint regions, especially genes implicated in reproduction and pheromone detection and mating. Moreover, we found an association of the evolutionary breakpoint regions with active chromatin state landscapes, most probably related to gene enrichment. Our results have two important implications for understanding the mechanisms that govern and constrain mammalian genome evolution. The first is that the presence of genes related to species-specific phenotypes in evolutionary breakpoint regions reinforces the adaptive value of genome reshuffling. Second, that chromatin conformation, an aspect that has been often overlooked in comparative genomic studies, might play a role in modeling the genomic distribution of evolutionary breakpoints.

Early embryonic development and spatiotemporal localization of mammalian primordial germ cell-associated proteins in the basal rodent Lagostomus maximus

The gene network controlling primordial germ cell (PGC) specification in eutherian mammals has been exhaustively investigated in mice. The egg-cylinder morphology of the mouse embryo is the key event enabling inductive signals from the extra-embryonic ectoderm (ExE) to specify epiblast cells as PGCs early on. We investigated the embryonic development and the spatiotemporal localization of PGC-associated proteins in the basal Hystricognathi rodent Lagostomus maximus. L. maximus develops through a flat-disc epiblast far apart from the ExE. In the primitive streak stage, OCT4-positive cells are detected in the posterior pole of the embryo disc in the mesoderm of the proximal epiblast. In the neural plate stage, a reduced 8 to 12 OCT4-positive cell population transiently expresses FRAGILIS, STELLA and SOX17 in the posterior streak. Soon after translocation to the hindgut, pluripotent OCT4 cells start expressing VASA, and then, STELLA and FRAGILIS are turned on during migration toward the genital ridge. L. maximus shows a spatiotemporal pattern of PGC-associated markers divergent from the early PGC restriction model seen in mice. This pattern conforms to alternative models that are based on a pluripotent population in the embryonic axis, where PGCs are specified later during development.

Gene structure of the pregnancy-associated glycoprotein-like (PAG-L) in the Eurasian beaver (Castor fiber L.)

The pregnancy-associated glycoprotein-like family (PAG-L) is a large group of chorionic products, expressed in the pre-placental trophoblast and later in the post-implantational chorionic epithelium, and are involved in proper placenta development and embryo-maternal interaction in eutherians. This study describes identification of the PAG-L family in the genome of the Eurasian beaver (Castor fiber L.), named CfPAG-L. We identified 7657 bp of the CfPAG-L gDNA sequence (Acc. No. KX377932), encompassing nine exons (1-9) and eight introns (A-H). The length of the CfPAG-L exons (59-200 bp) was equivalently similar to the only known counterparts of bPAG1, bPAG2, and pPAG2. The length of the CfPAG-L introns ranged 288-1937 bp and was completely different from previously known PAG introns. The exonic CfPAG-L regions revealed 50.3-72.9% homology with equivalent segments of bPAG1 and pPAG2 structure. The intronic CfPAG-L regions alignments revealed a lack of homology. Within the entire CfPAG-L gene, 31 potential single nucleotide variants (SNV: 7 transversions and 24 transitions) were predicted. The identified exonic polymorphic loci did not affect the amino acid sequence of the CfPAG-L polypeptide precursor. This is the first report describing the CfPAG-L gene sequence, structural organization, and SNVs in the Eurasian beaver, one of the largest rodents.

The Beaver's Phylogenetic Lineage Illuminated by Retroposon Reads

Solving problematic phylogenetic relationships often requires high quality genome data. However, for many organisms such data are still not available. Among rodents, the phylogenetic position of the beaver has always attracted special interest. The arrangement of the beaver's masseter (jaw-closer) muscle once suggested a strong affinity to some sciurid rodents (e.g., squirrels), placing them in the Sciuromorpha suborder. Modern molecular data, however, suggested a closer relationship of beaver to the representatives of the mouse-related clade, but significant data from virtually homoplasy-free markers (for example retroposon insertions) for the exact position of the beaver have not been available. We derived a gross genome assembly from deposited genomic Illumina paired-end reads and extracted thousands of potential phylogenetically informative retroposon markers using the new bioinformatics coordinate extractor fastCOEX, enabling us to evaluate different hypotheses for the phylogenetic position of the beaver. Comparative results provided significant support for a clear relationship between beavers (Castoridae) and kangaroo rat-related species (Geomyoidea) (p < 0.0015, six markers, no conflicting data) within a significantly supported mouse-related clade (including Myodonta, Anomaluromorpha, and Castorimorpha) (p < 0.0015, six markers, no conflicting data).

Bony labyrinth morphometry indicates locomotor adaptations in the squirrel-related clade (Rodentia, Mammalia)

The semicircular canals (SCs) of the inner ear detect angular acceleration and are located in the bony labyrinth of the petrosal bone. Based on high-resolution computed tomography, we created a size-independent database of the bony labyrinth of 50 mammalian species especially rodents of the squirrel-related clade comprising taxa with fossorial, arboreal and gliding adaptations. Our sampling also includes gliding marsupials, actively flying bats, the arboreal tree shrew and subterranean species. The morphometric anatomy of the SCs was correlated to the locomotion mode. Even if the phylogenetic signal cannot entirely be excluded, the main significance for functional morphological studies has been found in the diameter of the SCs, whereas the radius of curvature is of minor interest. Additionally, we found clear differences in the bias angle of the canals between subterranean and gliding taxa, but also between sciurids and glirids. The sensitivity of the inner ear correlates with the locomotion mode, with a higher sensitivity of the SCs in fossorial species than in flying taxa. We conclude that the inner ear of flying and gliding mammals is less sensitive due to the large information flow into this sense organ during locomotion.

Genetic adaptations of the plateau zokor in high-elevation burrows

The plateau zokor (Myospalax baileyi) spends its entire life underground in sealed burrows. Confronting limited oxygen and high carbon dioxide concentrations, and complete darkness, they epitomize a successful physiological adaptation. Here, we employ transcriptome sequencing to explore the genetic underpinnings of their adaptations to this unique habitat. Compared to Rattus norvegicus, genes belonging to GO categories related to energy metabolism (e.g. mitochondrion and fatty acid beta-oxidation) underwent accelerated evolution in the plateau zokor. Furthermore, the numbers of positively selected genes were significantly enriched in the gene categories involved in ATPase activity, blood vessel development and respiratory gaseous exchange, functional categories that are relevant to adaptation to high altitudes. Among the 787 genes with evidence of parallel evolution, and thus identified as candidate genes, several GO categories (e.g. response to hypoxia, oxygen homeostasis and erythrocyte homeostasis) are significantly enriched, are two genes, EPAS1 and AJUBA, involved in the response to hypoxia, where the parallel evolved sites are at positions that are highly conserved in sequence alignments from multiple species. Thus, accelerated evolution of GO categories, positive selection and parallel evolution at the molecular level provide evidences to parse the genetic adaptations of the plateau zokor for living in high-elevation burrows.

Capture of syncytin-Mar1, a fusogenic endogenous retroviral envelope gene involved in placentation in the Rodentia squirrel-related clade

Syncytin genes are fusogenic envelope protein (env) genes of retroviral origin that have been captured for a function in placentation. Within rodents, two such genes have previously been identified in the mouse-related clade, allowing a demonstration of their essential role via knockout mice. Here, we searched for similar genes in a second major clade of the Rodentia order, the squirrel-related clade, taking advantage of the complete sequencing of the ground squirrel Ictidomys tridecemlineatus genome. In silico search for env genes with full coding capacity identified several candidate genes with one displaying placenta-specific expression, as revealed by quantitative reverse transcription-PCR analysis of a large panel of tissues. This gene belongs to a degenerate endogenous retroviral element, with recognizable hallmarks of an integrated provirus. Cloning of the gene in an expression vector for ex vivo cell-cell fusion and pseudotype assays demonstrated fusogenicity on a large panel of mammalian cells. In situ hybridization on placenta sections showed specific expression in domains where trophoblast cells fuse into a syncytiotrophoblast at the fetomaternal interface, consistent with a role in syncytium formation. Finally, we show that the gene is conserved among the tribe Marmotini, thus dating its capture back to about at least 25 million years ago, with evidence for purifying selection and conservation of fusogenic activity. This gene that we named syncytin-Mar1 is distinct from all seven Syncytin genes identified to date in eutherian mammals and is likely to be a major effector of placentation in its related clade. Importance: Syncytin genes are fusogenic envelope genes of retroviral origin, ancestrally captured for a function in placentation. Within rodents, two such genes had been previously identified in the mouse-related clade. Here, in the squirrel-related rodent clade, we identified the envelope gene of an endogenous retrovirus with all the features of a Syncytin: it is specifically expressed in the placenta of the woodchuck Marmota monax, at the level of cells fusing into a syncytium; it can trigger cell-cell and virus-cell fusion ex vivo; and it has been conserved for >25 million years of evolution, suggesting an essential role in its host physiology. Remarkably, syncytin-Mar1 is unrelated to all other Syncytin genes identified thus far in mammals (primates, muroids, carnivores, and ruminants). These results extend the range of retroviral envelope gene "domestication" in mammals and show that these events occurred independently, on multiple occasions during evolution to improve placental development in a process of convergent evolution.

The evolution of micro-cursoriality in mammals

In this study we report on the evolution of micro-cursoriality, a unique case of cursoriality in mammals smaller than 1 kg. We obtained new running speed and limb morphology data for two species of elephant-shrews (Elephantulus spp., Macroscelidae) from Namaqualand, South Africa, which we compared with published data for other mammals. Elephantulus maximum running speeds were higher than most mammals smaller than 1 kg. Elephantulus also possess exceptionally high metatarsal:femur ratios (1.07) that are typically associated with fast unguligrade cursors. Cursoriality evolved in the Artiodactyla, Perissodactyla, and Carnivora coincident with global cooling and the replacement of forests with open landscapes in the Oligocene and Miocene. The majority of mammal species, though, remained non-cursorial, plantigrade, and small (&lt; 1 kg). The extraordinary running speed and digitigrady of elephant-shrews was established in the Early Eocene in the earliest macroscelid Prodiacodon, but was probably inherited from Paleocene, Holarctic stem macroscelids. Micro-cursoriality in macroscelids evolved from the plesiomorphic plantigrade foot of the possum-like ancestral mammal earlier than in other mammalian crown groups. Micro-cursoriality evolved first in forests, presumably in response to selection for rapid running speeds facilitated by local knowledge, in order to avoid predators. During the Miocene, micro-cursoriality was pre-adaptive to open, arid habitats, and became more derived in the newly-evolved Elephantulus and Macroscelides elephant-shrews with trail running.

Morphological diversity and evolution of the spermatozoon in the mouse-related clade of rodents

Most species in the three highly speciose families of the mouse-related clade of rodents, the Muridae, Cricetidae, and Nesomyidae (superfamily Muroidea), have a highly complex sperm head in which there is an apical hook but there are few data available for the other related families of these rodents. In the current study, using light and electron microscopies, we investigated the structure of the spermatozoon in representative species of four other families within the mouse-related clade, the Dipodidae, Spalacidae, Pedetidae, and Heteromyidae, that diverged at or near the base of the muroid lineage. Our results indicate that a diverse array of sperm head shapes and tail lengths occurs but none of the species in the families Spalacidae, Dipodidae, or Pedetidae has a sperm head with an apical hook. By contrast, a rostrally extending apical hook is present in spermatozoa of members of the Family Heteromyidae which also invariably have comparatively long sperm tails. These findings suggest that the hook-shaped sperm head in the murid, cricetid, and nesomyid rodents evolved after divergence of this lineage from its common ancestor with the other families of the mouse-related clade, and that separate, and independent, convergent evolution of a similar sperm head form, and long sperm tail, occurred in the Heteromyidae.

Jumping and gliding rodents: mitogenomic affinities of Pedetidae and Anomaluridae deduced from an RNA-Seq approach

An RNA-Seq strategy was used to obtain the complete set of protein-coding mitochondrial genes from two rodent taxa. Thanks to the next generation sequencing (NGS) 454 approach, we determined the complete mitochondrial DNA genome from Graphiurus kelleni (Mammalia: Rodentia: Gliridae) and partial mitogenome from Pedetes capensis (Pedetidae), and compared them with published rodent and outgroup mitogenomes. We finished the mitogenome sequencing by a series of amplicons using conserved PCR primers to fill the gaps corresponding to tRNA, rRNA and control regions. Phylogenetic analyses of the mitogenomes suggest a well-supported rodent phylogeny in agreement with nuclear gene trees. Pedetes groups with Anomalurus into the clade Anomaluromorpha, while Graphiurus branches within the squirrel-related clade. Moreover, Pedetes+Anomalurus branch with Castor into the mouse-related clade. Our study demonstrates the utility of NGS for obtaining new mitochondrial genomes as well as the importance of choosing adequate models of sequence evolution to infer the phylogeny of rodents.

PuF, an antimetastatic and developmental signaling protein, interacts with the Alzheimer's amyloid-β precursor protein via a tissue-specific proximal regulatory element (PRE)

BACKGROUND

Alzheimer's disease (AD) is intimately tied to amyloid-β (Aβ) peptide. Extraneuronal brain plaques consisting primarily of Aβ aggregates are a hallmark of AD. Intraneuronal Aβ subunits are strongly implicated in disease progression. Protein sequence mutations of the Aβ precursor protein (APP) account for a small proportion of AD cases, suggesting that regulation of the associated gene (APP) may play a more important role in AD etiology. The APP promoter possesses a novel 30 nucleotide sequence, or "proximal regulatory element" (PRE), at -76/-47, from the +1 transcription start site that confers cell type specificity. This PRE contains sequences that make it vulnerable to epigenetic modification and may present a viable target for drug studies. We examined PRE-nuclear protein interaction by gel electrophoretic mobility shift assay (EMSA) and PRE mutant EMSA. This was followed by functional studies of PRE mutant/reporter gene fusion clones.

RESULTS

EMSA probed with the PRE showed DNA-protein interaction in multiple nuclear extracts and in human brain tissue nuclear extract in a tissue-type specific manner. We identified transcription factors that are likely to bind the PRE, using competition gel shift and gel supershift: Activator protein 2 (AP2), nm23 nucleoside diphosphate kinase/metastatic inhibitory protein (PuF), and specificity protein 1 (SP1). These sites crossed a known single nucleotide polymorphism (SNP). EMSA with PRE mutants and promoter/reporter clone transfection analysis further implicated PuF in cells and extracts. Functional assays of mutant/reporter clone transfections were evaluated by ELISA of reporter protein levels. EMSA and ELISA results correlated by meta-analysis.

CONCLUSIONS

We propose that PuF may regulate the APP gene promoter and that AD risk may be increased by interference with PuF regulation at the PRE. PuF is targeted by calcium/calmodulin-dependent protein kinase II inhibitor 1, which also interacts with the integrins. These proteins are connected to vital cellular and neurological functions. In addition, the transcription factor PuF is a known inhibitor of metastasis and regulates cell growth during development. Given that APP is a known cell adhesion protein and ferroxidase, this suggests biochemical links among cell signaling, the cell cycle, iron metabolism in cancer, and AD in the context of overall aging.

Complete mitochondrial genome of the Eurasian flying squirrel Pteromys volans (Sciuromorpha, Sciuridae) and revision of rodent phylogeny

In this study, the complete mitochondrial genome of the Eurasian flying squirrel Pteromys volans (Rodentia, Sciuromorpha, Sciuridae) was sequenced and characterized in detail. The entire mitochondrial genome of P. volans consisted of 16,513 bp and contained 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and two non-coding regions. Its gene arrangement pattern was consistent with the mammalian ground pattern. The overall base composition and AT contents were similar to those of other rodent mitochondrial genomes. The light-strand origin generally identified between tRNA ( Asn ) and tRNA ( Cys ) consisted of a secondary structure with an 11-bp stem and an 11-bp loop. The large control region was constructed of three characteristic domains, ETAS, CD, and CSB without any repeat sequences. Each domain contained ETAS1, subsequences A, B, and C, and CSB1, respectively. In order to examine phylogenetic contentious issues of the monophyly of rodents and phylogenetic relationships among five rodent suborders, here, phylogenetic analyses based on nucleotide sequence data of the 35 rodent and 3 lagomorph mitochondrial genomes were performed using the Bayesian inference and maximum likelihood method. The result strongly supported the rodent monophyly with high node confidence values (BP 100 % in ML and BPP 1.00 in BI) and also monophylies of four rodent suborders (BP 85-100 % in ML and BPP 1.00 in BI), except for Anomalumorpha in which only one species was examined here. Also, phylogenetic relationships among the five rodent suborders were suggested and discussed in detail.

Major radiations in the evolution of Caviid rodents: reconciling fossils, ghost lineages, and relaxed molecular clocks

BACKGROUND

Caviidae is a diverse group of caviomorph rodents that is broadly distributed in South America and is divided into three highly divergent extant lineages: Caviinae (cavies), Dolichotinae (maras), and Hydrochoerinae (capybaras). The fossil record of Caviidae is only abundant and diverse since the late Miocene. Caviids belongs to Cavioidea sensu stricto (Cavioidea s.s.) that also includes a diverse assemblage of extinct taxa recorded from the late Oligocene to the middle Miocene of South America ("eocardiids").

RESULTS

A phylogenetic analysis combining morphological and molecular data is presented here, evaluating the time of diversification of selected nodes based on the calibration of phylogenetic trees with fossil taxa and the use of relaxed molecular clocks. This analysis reveals three major phases of diversification in the evolutionary history of Cavioidea s.s. The first two phases involve two successive radiations of extinct lineages that occurred during the late Oligocene and the early Miocene. The third phase consists of the diversification of Caviidae. The initial split of caviids is dated as middle Miocene by the fossil record. This date falls within the 95% higher probability distribution estimated by the relaxed Bayesian molecular clock, although the mean age estimate ages are 3.5 to 7 Myr older. The initial split of caviids is followed by an obscure period of poor fossil record (referred here as the Mayoan gap) and then by the appearance of highly differentiated modern lineages of caviids, which evidentially occurred at the late Miocene as indicated by both the fossil record and molecular clock estimates.

CONCLUSIONS

The integrated approach used here allowed us identifying the agreements and discrepancies of the fossil record and molecular clock estimates on the timing of the major events in cavioid evolution, revealing evolutionary patterns that would not have been possible to gather using only molecular or paleontological data alone.

Small mammal investigation in spotted fever focus with DNA-barcoding and taxonomic implications on rodents species from Hainan of China

Although mammals are a well-studied group of animals, making accurate field identification of small mammals is still complex because of morphological variation across developmental stages, color variation of pelages, and often damaged osteological and dental characteristics. In 2008, small mammals were collected for an epidemiological study of a spotted fever outbreak in Hainan, China. Ten species of small mammals were identified by morphological characters in the field, most using pelage color characters only. The study is extended here, in order to assess whether DNA barcoding would be suitable as an identification tool in these small mammals. Barcode clusters showed some incongruence with morphospecies, especially for some species of Rattus and Niviventer, so molecular delineation was carried out with an expanded dataset of combined cytochrome b (Cyt-b) and cytochrome c oxidase subunit I (COI) sequences. COI sequences were successfully amplified from 83% of collected mammals, but failed in all specimens of Suncus murinus, which were thus excluded in DNA barcoding analysis. Of note, ten molecular taxonomic units were found from samples of nine morphologically identified species. Accordingly, 11 species of small mammals were present in the investigated areas, including four Rattus species, three Niviventer species, Callosciurus erythraeus, Neohylomys hainanensis, Tupaia belangeri, and Suncus murinus. Based on the results of the phylogenetic and molecular delineation analyses, the systematic status of some rodent species should be redefined. R. rattus hainanicus and R. rattus sladeni are synonyms of R. andamanensis. R. losea from China and Southeast Asia comprises two independent species: R. losea and R. sakeratensis. Finally, the taxonomic status of three putative species of Niviventer should be further confirmed according to morphological, molecular and ecological characters.

Small mammal investigation in spotted fever focus with DNA-barcoding and taxonomic implications on rodents species from Hainan of China

Although mammals are a well-studied group of animals, making accurate field identification of small mammals is still complex because of morphological variation across developmental stages, color variation of pelages, and often damaged osteological and dental characteristics. In 2008, small mammals were collected for an epidemiological study of a spotted fever outbreak in Hainan, China. Ten species of small mammals were identified by morphological characters in the field, most using pelage color characters only. The study is extended here, in order to assess whether DNA barcoding would be suitable as an identification tool in these small mammals. Barcode clusters showed some incongruence with morphospecies, especially for some species of Rattus and Niviventer, so molecular delineation was carried out with an expanded dataset of combined cytochrome b (Cyt-b) and cytochrome c oxidase subunit I (COI) sequences. COI sequences were successfully amplified from 83% of collected mammals, but failed in all specimens of Suncus murinus, which were thus excluded in DNA barcoding analysis. Of note, ten molecular taxonomic units were found from samples of nine morphologically identified species. Accordingly, 11 species of small mammals were present in the investigated areas, including four Rattus species, three Niviventer species, Callosciurus erythraeus, Neohylomys hainanensis, Tupaia belangeri, and Suncus murinus. Based on the results of the phylogenetic and molecular delineation analyses, the systematic status of some rodent species should be redefined. R. rattus hainanicus and R. rattus sladeni are synonyms of R. andamanensis. R. losea from China and Southeast Asia comprises two independent species: R. losea and R. sakeratensis. Finally, the taxonomic status of three putative species of Niviventer should be further confirmed according to morphological, molecular and ecological characters.

Low, complex and probably reticulated chromosome evolution of Sciuromorpha (Rodentia) and Lagomorpha

Lagomorpha (rabbits and pikas) and Sciuromorpha (squirrels) are grouped in the Glires superorder. Their chromosome diversification, since their separation from the eutherian mammalian common ancestor, was characterized by a low rate of chromosome rearrangements. Consequently, the structure of some chromosomes was either conserved or only slightly modified, making their comparison easy at the genus, family and even order level. Interspecific in situ hybridization (Zoo-FISH) largely corroborates classical cytogenetic data but provides much more reliability in comparisons, especially for distant species. We reconstructed common ancestral karyotypes for Glires, Lagomorpha, Sciuromorpha, and Sciuridae species, and then, determined the chromosome changes separating these ancestors from their common eutherian ancestor. We propose that reticulated evolution occurred during the diversification of Glires, which implies that several pericentric inversions and Robertsonian translocations were conserved in the heterozygous status for an extensive period. Finally, among Lagomorpha and Sciuromorpha, we focused on Leporidae and Sciuridae chromosome evolution. In the various attempts to establish dichotomic evolutionary schemes, it was necessary to admit that multiple homoplasies (convergent and reverse rearrangements) occurred in Sciuridae and in a lesser degree, in Leporidae. In Leporidae, additional rearrangements were sufficient to propose a resolved phylogeny. However, a resolved phylogeny was not possible for Sciuridae because most of the rearrangements occurred in terminal branches. We conclude that a reticulated evolution took place early during the evolution of both families and lasted longer in Sciuridae than in Leporidae. In Sciuridae, most chromosome rearrangements were pericentric inversions involving short fragments. Such rearrangements have only mild meiotic consequences, which may explain the long persistence of the heterozygous status characterizing reticulated evolution.

A glimpse on the pattern of rodent diversification: a phylogenetic approach

Background

Development of phylogenetic methods that do not rely on fossils for the study of evolutionary processes through time have revolutionized the field of evolutionary biology and resulted in an unprecedented expansion of our knowledge about the tree of life. These methods have helped to shed light on the macroevolution of many taxonomic groups such as the placentals (Mammalia). However, despite the increase of studies addressing the diversification patterns of organisms, no synthesis has addressed the case of the most diversified mammalian clade: the Rodentia.

Results

Here we present a rodent maximum likelihood phylogeny inferred from a molecular supermatrix. It is based on 11 mitochondrial and nuclear genes that covers 1,265 species, i.e., respectively 56% and 81% of the known specific and generic rodent diversity. The inferred topology recovered all Rodentia clades proposed by recent molecular works. A relaxed molecular clock dating approach provided a time framework for speciation events. We found that the Myomorpha clade shows a greater degree of variation in diversification rates than Sciuroidea, Caviomorpha, Castorimorpha and Anomaluromorpha. We identified a number of shifts in diversification rates within the major clades: two in Castorimorpha, three in Ctenohystrica, 6 within the squirrel-related clade and 24 in the Myomorpha clade. The majority of these shifts occurred within the most recent familial rodent radiations: the Cricetidae and Muridae clades. Using the topological imbalances and the time line we discuss the potential role of different diversification factors that might have shaped the rodents radiation.

Conclusions

The present glimpse on the diversification pattern of rodents can be used for further comparative meta-analyses. Muroid lineages have a greater degree of variation in their diversification rates than any other rodent group. Different topological signatures suggest distinct diversification processes among rodent lineages. In particular, Muroidea and Sciuroidea display widespread distribution and have undergone evolutionary and adaptive radiation on most of the continents. Our results show that rodents experienced shifts in diversification rate regularly through the Tertiary, but at different periods for each clade. A comparison between the rodent fossil record and our results suggest that extinction led to the loss of diversification signal for most of the Paleogene nodes.

Chromosomal evolution in Rodentia

Rodentia is the most species-rich mammalian order and includes several important laboratory model species. The amount of new information on karyotypic and phylogenetic relations within and among rodent taxa is rapidly increasing, but a synthesis of these data is currently lacking. Here, we have integrated information drawn from conventional banding studies, recent comparative painting investigations and molecular phylogenetic reconstructions of different rodent taxa. This permitted a revision of several ancestral karyotypic reconstructions, and a more accurate depiction of rodent chromosomal evolution.

Expensive Brains: "Brainy" Rodents have Higher Metabolic Rate

Brains are the centers of the nervous system of animals, controlling the organ systems of the body and coordinating responses to changes in the ecological and social environment. The evolution of traits that correlate with cognitive ability, such as relative brain size is thus of broad interest. Brain mass relative to body mass (BM) varies among mammals, and diverse factors have been proposed to explain this variation. A recent study provided evidence that energetics play an important role in brain evolution (Isler and van Schaik, 2006). Using composite phylogenies and data drawn from multiple sources, these authors showed that basal metabolic rate (BMR) correlates with brain mass across mammals. However, no such relationship was found within rodents. Here we re-examined the relationship between BMR and brain mass within Rodentia using a novel species-level phylogeny. Our results are sensitive to parameter evaluation; in particular how species mass is estimated. We detect no pattern when applying an approach used by previous studies, where each species BM is represented by two different numbers, one being the individual that happened to be used for BMR estimates of that species. However, this approach may compromise the analysis. When using a single value of BM for each species, whether representing a single individual, or available species mean, our findings provide evidence that brain mass (independent of BM) and BMR are correlated. These findings are thus consistent with the hypothesis that large brains evolve when the payoff for increased brain mass is greater than the energetic cost they incur.