Rodent Lineages Research Articles

A new Late Miocene stem Octodontoid (Rodentia: Caviomorpha: Octochinchilloi) from the Solimões Formation (Brazil): paleobiogeographic implications

ABSTRACT Caviomorphs are the neotropical hystricognathous rodents (Ctenohystrica: Hystricognathi: Caviomorpha). Octodontoidea are today the richest and most speciose clade of caviomorphs. Their fossil record also reveals a wide array of extinct species, including Acarechimys, a taxon with a widespread South American distribution and extensive temporal range. This genus comprises at least five species documented in deposits ranging from the late Oligocene to the Middle Miocene. Here, we report the discovery of a new species of Acarechimys from Upper Miocene deposits of the Envira River (PRE06 outcrop), situated in the Acre State, western Brazilian Amazonia. The new species differs from the other species of the genus in having lower molars with a strong posterior arm of the metaconid, the presence of a mesostylid (usually prominent) and a neomesolophid forwardly directed and connected to the metalophulid I on m1 but not on m2–3, and the absence of posterior arm of the protoconid. This taxon represents the youngest known species of the Acarechimys lineage to date. Its occurrence in Western Amazonia strengthens support for the hypothesis that South American rodent lineages survived longer during the Neogene at low latitudes than in the southern regions of the continent.

Genomic insights into adaptation to bipedal saltation and desert-like habitats of jerboas.

Jerboas is a lineage of small rodents displaying atypical mouse-like morphology with elongated strong hindlimbs and short forelimbs. They have evolved obligate bipedal saltation and acute senses, and been well-adapted to vast desert-like habitats. Using a newly sequenced chromosome-scale genome of the Mongolian five-toed jerboa (Orientallactaga sibirica), our comparative genomic analyses and in vitro functional assays showed that the genetic innovations in both protein-coding and non-coding regions played an important role in jerboa morphological and physiological adaptation. Jerboa-specific amino acid substitutions, and segment insertions/deletions (indels) in conserved non-coding elements (CNEs) were found in components of proteoglycan biosynthesis pathway (XYLT1 and CHSY1), which plays an important role in limb development. Meanwhile, we found specific evolutionary changes functionally associated with energy or water metabolism (e.g., specific amino acid substitutions in ND5 and indels in CNEs physically near ROR2) and senses (e.g., expansion of vomeronasal receptors and the FAM136A gene family) in jerboas. Further dual-luciferase reporter assay verified that some of the CNEs with jerboa-specific segment indels exerted a significantly different influence on luciferase activity, suggesting changes in their regulatory function in jerboas. Our results revealed the potential molecular mechanisms underlying jerboa adaptation since the divergence from the Eocene-Oligocene transition, and provided more resources and new insights to enhance our understanding of the molecular basis underlying the phenotypic diversity and the environmental adaptation of mammals.

Phylogenetic relationships and systematics of tapeworms of the family Davaineidae (Cestoda, Cyclophyllidea), with emphasis on species in rodents.

The present study aims at clarifying the poorly known phylogenetic relationships and systematics of cestodes of the family Davaineidae Braun, 1900 (Cyclophyllidea), primarily the genus Raillietina Fuhrmann, 1920 and of the subfamily Inermicapsiferinae (Anoplocephalidae) from mammals (mostly rodents, 31 new isolates) and birds (eight new isolates). Phylogenetic analyses are based on sequences of the large subunit ribosomal RNA gene (28S) and mitochondrial NADH dehydrogenase subunit 1 gene (nad1). The main phylogenetic pattern emerging from the present analysis is the presence of three independent lineages within the main clade of the subfamily Davaineinae, one of which is almost entirely confined to species from rodents and the other two show a mixture of species from birds and mammals. It is suggested that the major diversification of the main clade took place in birds, possibly in galliforms. The subsequent diversification included repeated host shifts from birds to mammals and to other birds, and from rodents to other mammals, showing that colonisation of new host lineages has been the main driver in the diversification of davaineine cestodes. It is also shown that all isolates of Inermicapsifer Janicki, 1910, mainly from rodents, form a monophyletic group positioned among Raillietina spp. in the "rodent lineage", indicating that the genus Inermicapsifer is a member of the family Davaineidae. This means that the subfamily Inermicapsiferinae and the family Inermicapsiferidae should be treated as synonyms of the Davaineidae, specifically the subfamily Davaineinae. Three additional genera generally included in the Inermicapsiferinae, i.e. Metacapsifer Spasskii, 1951, Pericapsifer Spasskii, 1951 and Thysanotaenia Beddard, 1911, are also assigned here to the Davaineidae (subfamily Davaineinae). Raillietina spp. were present in all three main lineages and appeared as multiple independent sublineages from bird and mammalian hosts, verifying the non-monophyly of the genus Raillietina and suggesting a presence of multiple new species and genera.

The structure of the TH/INS locus and the parental allele expressed are not conserved between mammals

Parent-of-origin-specific expression of imprinted genes is critical for successful mammalian growth and development. Insulin, coded by the INS gene, is an important growth factor expressed from the paternal allele in the yolk sac placenta of therian mammals. The tyrosine hydroxylase gene TH encodes an enzyme involved in dopamine synthesis. TH and INS are closely associated in most vertebrates, but the mouse orthologues, Th and Ins2, are separated by repeated DNA. In mice, Th is expressed from the maternal allele, but the parental origin of expression is not known for any other mammal so it is unclear whether the maternal expression observed in the mouse represents an evolutionary divergence or an ancestral condition. We compared the length of the DNA segment between TH and INS across species and show that separation of these genes occurred in the rodent lineage with an accumulation of repeated DNA. We found that the region containing TH and INS in the tammar wallaby produces at least five distinct RNA transcripts: TH, TH-INS1, TH-INS2, lncINS and INS. Using allele-specific expression analysis, we show that the TH/INS locus is expressed from the paternal allele in pre- and postnatal tammar wallaby tissues. Determining the imprinting pattern of TH/INS in other mammals might clarify if paternal expression is the ancestral condition which has been flipped to maternal expression in rodents by the accumulation of repeat sequences.

Does climate shape geographical morphological variation of spines in Neotropical spiny rats (Rodentia, Echimyidae, Trinomys)?

Abstract Flattened and flexible spiny hairs, the aristiforms, evolved independently in some lineages of rodents, but unlikely capable of causing severe damage to predators. Their biological roles are not yet clearly known and some different thermoregulatory functions have been proposed, often associated with severe climatic conditions but not yet tested. In this study we investigated whether climate may be a major factor shaping geographic variation of length and width of aristiform hairs of 2 species of Neotropical spiny rats of the genus Trinomys (Echimyidae) occurring in eastern Brazil: T. albispinus and T. setosus. Results revealed clinal variation in these hairs. While bioclimatic factors and geographic distance together accounted for up to a third of the morphometric variation, the exclusive impact of climate was limited, explaining only up to 3% of the observed variation. However, considering the close association between climate and geographic distance, the role of climate in shaping aristiform hair cannot be discounted. While our results did not fully support the hypothesis that well-developed, wide spiny hairs are selected for greater protection against rain for underhairs, they marginally supported their role in facilitating heat loss in warm environments by facilitating water condensation. More consistently, results show that narrow spiny hairs tended to evolve in cold environments, likely enhancing thermal insulation. The length of aristiform hair in T. albispinus also appeared to respond to this function, becoming longer in cold and high-altitude regions. However, we cannot rule out the possibility of an evolutionary response correlated with selection on underhairs. Future experimental studies should test the thermoregulatory function of flexible spiny hairs and investigate the geographic variation of these structures in other rodent species.

Insertion-Deletion Events Are Depleted in Protein Regions with Predicted Secondary Structure.

A fundamental goal in evolutionary biology and population genetics is to understand how selection shapes the fate of new mutations. Here, we test the null hypothesis that insertion-deletion (indel) events in protein-coding regions occur randomly with respect to secondary structures. We identified indels across 11,444 sequence alignments in mouse, rat, human, chimp, and dog genomes and then quantified their overlap with four different types of secondary structure-alpha helices, beta strands, protein bends, and protein turns-predicted by deep-learning methods of AlphaFold2. Indels overlapped secondary structures 54% as much as expected and were especially underrepresented over beta strands, which tend to form internal, stable regions of proteins. In contrast, indels were enriched by 155% over regions without any predicted secondary structures. These skews were stronger in the rodent lineages compared to the primate lineages, consistent with population genetic theory predicting that natural selection will be more efficient in species with larger effective population sizes. Nonsynonymous substitutions were also less common in regions of protein secondary structure, although not as strongly reduced as in indels. In a complementary analysis of thousands of human genomes, we showed that indels overlapping secondary structure segregated at significantly lower frequency than indels outside of secondary structure. Taken together, our study shows that indels are selected against if they overlap secondary structure, presumably because they disrupt the tertiary structure and function of a protein.

Where Are the Formerly Y-linked Genes in the Ryukyu Spiny Rat that has Lost its Y Chromosome?

It has been predicted that the highly degenerate mammalian Y chromosome will be lost eventually. Indeed, Y was lost in the Ryukyu spiny rat Tokudaia osimensis, but the fate of the formerly Y-linked genes is not completely known. We looked for all 12 ancestrally Y-linked genes in a draft T. osimensis genome sequence. Zfy1, Zfy2, Kdm5d, Eif2s3y, Usp9y, Uty, and Ddx3y are putatively functional and are now located on the X chromosome, whereas Rbmy, Uba1y, Ssty1, Ssty2, and Sry are missing or pseudogenized. Tissue expressions of the mouse orthologs of the retained genes are significantly broader/higher than those of the lost genes, suggesting that the destinies of the formerly Y-linked genes are related to their original expressions. Interestingly, patterns of gene retention/loss are significantly more similar than by chance across four rodent lineages where Y has been independently lost, indicating a level of certainty in the fate of Y-linked genes even when the chromosome is gone.

Proportional variation and scaling in the hindlimbs of hopping mammals, including convergent evolution in argyrolagids and jerboas

Bipedal hopping is a mode of locomotion seen today in four rodent lineages and one clade of marsupials. The Argyrolagidae, marsupials from the Oligocene to Pliocene of South America, have also been considered to be hoppers. These lineages all convergently evolved similar general morphologies, with elongated hindlimbs, reduced forelimbs, and elongated tails, and their similarities and variations may be informative in understanding the evolution of hopping in mammals. This study uses principal components analysis and log-log regressions to investigate variation in the hindlimb proportions of these hopping mammals and how this relates to body mass. We find that the distribution of hopping mammal masses is bimodal, divided at roughly 500 g. These two domains among hopping mammals may reflect optimisation for different forms of hopping locomotion; species under 500 g tend to have more elongated metatarsals relative to the rest of their hindlimbs, perhaps to facilitate rapid vertical jumps for predator evasion, a behaviour not seen in larger hoppers. Despite this bimodal distribution in body mass, hindlimb proportions cluster more by clade than mass, with some similarities among clades being especially noteworthy. The jerboas (Dipodidae, Rodentia) and Argyrolagidae share a particularly extreme degree of metatarsal elongation. The drivers of this convergence are unclear, but we hypothesise that the elongation may be related to the reduction/fusion of metatarsals in these groups, or a greater reliance on bipedality at slow speeds, as jerboas are known to utilise multiple bipedal gaits in addition to hopping.

Irreversible evolutionary loss of chitin-degrading ability in the chitinase-like protein Ym1 under positive selection in rodents.

Ym1 (chitinase-like 3, Chil3) expressed in mice is a nonenzymatic chitinase-like protein, which shows 67% identity with mouse acidic chitinase (Chia). Similar to Chia, Ym1 is overexpressed in asthma and parasitic infections in mouse lungs. Due to the lack of chitin-degrading activity, the biomedical role of Ym1 under these pathophysiological conditions remains to be determined. In this study, we investigated what region and amino acid changes in Ym1 resulted in the loss of enzymatic activity. Replacing two amino acids at the catalytic motif to obtain a Chia-like sequence (N136D and Q140E; MT-Ym1) did not activate the protein. We conducted a comparative study of Ym1 and Chia. We found that three protein segments-(i) the catalytic motif residues, (ii) exons 6 and 7, and (iii) exon 10-are responsible for chitinase activity loss in Ym1. We show that replacing each of these three segments in Chia that are also involved in substrate recognition and binding by the Ym1 sequence can fully abolish the enzymatic activity. In addition, we show that there have been extensive gene duplication events at the Ym1 locus specific to the rodent lineages. Consistent with this result, Ym1 orthologs from the rodent genome were under positive selection when analyzed through the CODEML program. These data suggest that numerous amino acid substitutions in the regions involved in the chitin recognition, binding, and degradation ability of the ancestor Ym1 molecule lead to the irreversible inactivation of the protein.

A phylogenomic analysis of Limosilactobacillus reuteri reveals ancient and stable evolutionary relationships with rodents and birds and zoonotic transmission to humans

BackgroundGut microbes play crucial roles in the development and health of their animal hosts. However, the evolutionary relationships of gut microbes with vertebrate hosts, and the consequences that arise for the ecology and lifestyle of the microbes are still insufficiently understood. Specifically, the mechanisms by which strain-level diversity evolved, the degree by which lineages remain stably associated with hosts, and how their evolutionary history influences their ecological performance remain a critical gap in our understanding of vertebrate-microbe symbiosis.ResultsThis study presents the characterization of an extended collection of strains of Limosilactobacillus reuteri and closely related species from a wide variety of hosts by phylogenomic and comparative genomic analyses combined with colonization experiments in mice to gain insight into the long-term evolutionary relationship of a bacterial symbiont with vertebrates. The phylogenetic analysis of L. reuteri revealed early-branching lineages that primarily consist of isolates from rodents (four lineages) and birds (one lineage), while lineages dominated by strains from herbivores, humans, pigs, and primates arose more recently and were less host specific. Strains from rodent lineages, despite their phylogenetic divergence, showed tight clustering in gene-content-based analyses. These L. reuteri strains but not those ones from non-rodent lineages efficiently colonize the forestomach epithelium of germ-free mice. The findings support a long-term evolutionary relationships of L. reuteri lineages with rodents and a stable host switch to birds. Associations of L. reuteri with other host species are likely more dynamic and transient. Interestingly, human isolates of L. reuteri cluster phylogenetically closely with strains from domesticated animals, such as chickens and herbivores, suggesting zoonotic transmissions.ConclusionsOverall, this study demonstrates that the evolutionary relationship of a vertebrate gut symbiont can be stable in particular hosts over time scales that allow major adaptations and specialization, but also emphasizes the diversity of symbiont lifestyles even within a single bacterial species. For L. reuteri, symbiont lifestyles ranged from autochthonous, likely based on vertical transmission and stably aligned to rodents and birds over evolutionary time, to allochthonous possibly reliant on zoonotic transmission in humans. Such information contributes to our ability to use these microbes in microbial-based therapeutics.

DNA barcodes highlight genetic diversity patterns in rodents from lowland desert and Andean areas in Argentina

Rodents are an important component of South America fauna. Their high diversity has motivated researchers to continually review their taxonomy, genetic diversity, species limits, and phylogenetic relationships. Here, we applied DNA-barcodes for assessing the taxonomic and genetic diversity in the two major lineages of South American rodents: caviomorphs and sigmodontines. We analysed 335 COI barcodes in 34morphologically determined species from 39localities along central Andes and arid lands of Argentina. Neighbour-joining and maximum likelihood reconstruction provided clear separation between species. The Barcode Index number and Bayesian Poisson tree processes were used to confirm concordance between sequence clusters and species designations by taxonomy. We found deep divergence within the Phyllotis xanthopygus species complex, with distances up to 13.0% between geographically separated lineages. Minor divergences (3.30% and 2.52%) were found within Abrothrix hirta, and Tympanoctomys barrerae, respectively, with differentiation in their genetic lineages. Also, we documented geographically separated clusters for Akodon spegazzinii and A.oenos with up to 2.3% divergence, but clustering methods failed to distinguish them as different species. Sequence results show a clear barcode gap with a mean intraspecific divergence (0.56%) versus a minimum nearest-neighbour distance averaging (10.1%). Distances between congeneric species varied from 4.1 to 14%, with the exception of two related forms within Euneomys and the sister species Akodon spegazzinii and A.oenos. This study constitutes a substantial contribution to the global barcode reference library. It provides insights into the complex phylogeographic patterns and speciation scenarios in rodents, while highlighting areas that require in-depth taxonomic and integrative research.

Secondary loss of miR-3607 reduced cortical progenitor amplification during rodent evolution.

The evolutionary expansion and folding of the mammalian cerebral cortex resulted from amplification of progenitor cells during embryonic development. This process was reversed in the rodent lineage after splitting from primates, leading to smaller and smooth brains. Genetic mechanisms underlying this secondary loss in rodent evolution remain unknown. We show that microRNA miR-3607 is expressed embryonically in the large cortex of primates and ferret, distant from the primate-rodent lineage, but not in mouse. Experimental expression of miR-3607 in embryonic mouse cortex led to increased Wnt/β-catenin signaling, amplification of radial glia cells (RGCs), and expansion of the ventricular zone (VZ), via blocking the β-catenin inhibitor APC (adenomatous polyposis coli). Accordingly, loss of endogenous miR-3607 in ferret reduced RGC proliferation, while overexpression in human cerebral organoids promoted VZ expansion. Our results identify a gene selected for secondary loss during mammalian evolution to limit RGC amplification and, potentially, cortex size in rodents.

Comparative morphology of the dormouse skull and the influence of size and ecology.

Dormice are widely dispersed across various ecosystems in Eurasia and Africa and among the oldest extant rodent lineages. Despite distinct morphological variation to be present between groups, comprehensive morphometrical studies on the dormouse skull are limited. Here, the form of eight out of the nine extant dormouse genera was analysed using 3D geometric morphometrics and linear biomechanical measurements, providing a better understanding of the overall morphological variation present within Gliridae. Species-, genus- and family-specific morphological trends in both the size and shape of the cranium and mandible are linked with specific habitats and feeding strategies. Smaller dormice show adaptations to a more arboreal lifestyle such as a relatively enlarged braincase and an inferiorly reoriented foramen magnum. Larger dormice show cranial modifications, including clear flattening of the skull and a more posteriorly positioned foramen magnum, hinting towards a more rupicolous lifestyle. Furthermore, specimens inhabiting arid areas appear to have more inflated auditory bullae, whereas other variable features, such as the length of the incisive foramen, were not associated with either size changes or climatic variables. Lastly, more robust and horizontally orientated zygomatic arches as well as increased robusticity of the molar row appear to be linked with herbivory in dormice, whereas thinner arches and small concave molars are seen in more insectivorous species. This study reveals clear convergent adaptations between dormouse species and results in a better understanding of ecological drivers underpinning the morphological divergence present within Gliridae.

The Evolutionary History of the Chymase Locus -a Locus Encoding Several of the Major Hematopoietic Serine Proteases.

Several hematopoietic cells of the immune system store large amounts of proteases in cytoplasmic granules. The absolute majority of these proteases belong to the large family of chymotrypsin-related serine proteases. The chymase locus is one of four loci encoding these granule-associated serine proteases in mammals. The chymase locus encodes only four genes in primates, (1) the gene for a mast-cell-specific chymotryptic enzyme, the chymase; (2) a T-cell-expressed asp-ase, granzyme B; (3) a neutrophil-expressed chymotryptic enzyme, cathepsin G; and (4) a T-cell-expressed chymotryptic enzyme named granzyme H. Interestingly, this locus has experienced a number of quite dramatic expansions during mammalian evolution. This is illustrated by the very large number of functional protease genes found in the chymase locus of mice (15 genes) and rats (18 genes). A separate expansion has also occurred in ruminants, where we find a new class of protease genes, the duodenases, which are expressed in the intestinal region. In contrast, the opossum has only two functional genes in this locus, the mast cell (MC) chymase and granzyme B. This low number of genes may be the result of an inversion, which may have hindered unequal crossing over, a mechanism which may have been a major factor in the expansion within the rodent lineage. The chymase locus can be traced back to early tetrapods as genes that cluster with the mammalian genes in phylogenetic trees can be found in frogs, alligators and turtles, but appear to have been lost in birds. We here present the collected data concerning the evolution of this rapidly evolving locus, and how these changes in gene numbers and specificities may have affected the immune functions in the various tetrapod species.

Widespread loss of mammalian lineage and dietary diversity in the early Oligocene of Afro-Arabia

Diverse lines of geological and geochemical evidence indicate that the Eocene-Oligocene transition (EOT) marked the onset of a global cooling phase, rapid growth of the Antarctic ice sheet, and a worldwide drop in sea level. Paleontologists have established that shifts in mammalian community structure in Europe and Asia were broadly coincident with these events, but the potential impact of early Oligocene climate change on the mammalian communities of Afro-Arabia has long been unclear. Here we employ dated phylogenies of multiple endemic Afro-Arabian mammal clades (anomaluroid and hystricognath rodents, anthropoid and strepsirrhine primates, and carnivorous hyaenodonts) to investigate lineage diversification and loss since the early Eocene. These analyses provide evidence for widespread mammalian extinction in the early Oligocene of Afro-Arabia, with almost two-thirds of peak late Eocene diversity lost in these clades by ~30 Ma. Using homology-free dental topographic metrics, we further demonstrate that the loss of Afro-Arabian rodent and primate lineages was associated with a major reduction in molar occlusal topographic disparity, suggesting a correlated loss of dietary diversity. These results raise new questions about the relative importance of global versus local influences in shaping the evolutionary trajectories of Afro-Arabia’s endemic mammals during the Oligocene.

Two Functional Epithelial Sodium Channel Isoforms Are Present in Rodents despite Pronounced Evolutionary Pseudogenization and Exon Fusion.

The epithelial sodium channel (ENaC) plays a key role in salt and water homeostasis in tetrapod vertebrates. There are four ENaC subunits (α, β, γ, δ), forming heterotrimeric αβγ- or δβγ-ENaCs. Although the physiology of αβγ-ENaC is well understood, for decades the field has stalled with respect to δβγ-ENaC due to the lack of mammalian model organisms. The SCNN1D gene coding for δ-ENaC was previously believed to be absent in rodents, hindering studies using standard laboratory animals. We analyzed all currently available rodent genomes and discovered that SCNN1D is present in rodents but was independently lost in five rodent lineages, including the Muridae (mice and rats). The independent loss of SCNN1D in rodent lineages may be constrained by phylogeny and taxon-specific adaptation to dry habitats, however habitat aridity does not provide a selection pressure for maintenance of SCNN1D across Rodentia. A fusion of two exons coding for a structurally flexible region in the extracellular domain of δ-ENaC appeared in the Hystricognathi (a group that includes guinea pigs). This conserved pattern evolved at least 41 Ma and represents a new autapomorphic feature for this clade. Exon fusion does not impair functionality of guinea pig (Cavia porcellus) δβγ-ENaC expressed in Xenopus oocytes. Electrophysiological characterization at the whole-cell and single-channel level revealed conserved biophysical features and mechanisms controlling guinea pig αβγ- and δβγ-ENaC function as compared with human orthologs. Guinea pigs therefore represent commercially available mammalian model animals that will help shed light on the physiological function of δ-ENaC.

Compromised Function of the Pancreatic Transcription Factor PDX1 in a Lineage of Desert Rodents

Gerbils are a subfamily of rodents living in arid regions of Asia and Africa. Recent studies have shown that several gerbil species have unusual amino acid changes in the PDX1 protein, a homeodomain transcription factor essential for pancreatic development and β-cell function. These changes were linked to strong GC-bias in the genome that may be caused by GC-biased gene conversion, and it has been hypothesized that this caused accumulation of deleterious changes. Here we use two approaches to examine if the unusual changes are adaptive or deleterious. First, we compare PDX1 protein sequences between 38 rodents to test for association with habitat. We show the PDX1 homeodomain is almost totally conserved in rodents, apart from gerbils, regardless of habitat. Second, we use ectopic gene overexpression and gene editing in cell culture to compare functional properties of PDX1 proteins. We show that the divergent gerbil PDX1 protein inefficiently binds an insulin gene promoter and ineffectively regulates insulin expression in response to high glucose in rat cells. The protein has, however, retained the ability to regulate some other β-cell genes. We suggest that during the evolution of gerbils, the selection-blind process of biased gene conversion pushed fixation of mutations adversely affecting function of a normally conserved homeodomain protein. We argue these changes were not entirely adaptive and may be associated with metabolic disorders in gerbil species on high carbohydrate diets. This unusual pattern of molecular evolution could have had a constraining effect on habitat and diet choice in the gerbil lineage.

Rapid Evolution of HERC6 and Duplication of a Chimeric HERC5/6 Gene in Rodents and Bats Suggest an Overlooked Role of HERCs in Mammalian Immunity

Studying the evolutionary diversification of mammalian antiviral defenses is of main importance to better understand our innate immune repertoire. The small HERC proteins are part of a multigene family, including HERC5 and HERC6, which have probably diversified through complex evolutionary history in mammals. Here, we performed mammalian-wide phylogenetic and genomic analyses of HERC5 and HERC6, using 83 orthologous sequences from bats, rodents, primates, artiodactyls, and carnivores—the top five representative groups of mammalian evolution. We found that HERC5 has been under weak and differential positive selection in mammals, with only primate HERC5 showing evidences of pathogen-driven selection. In contrast, HERC6 has been under strong and recurrent adaptive evolution in mammals, suggesting past and widespread genetic arms-races with viral pathogens. Importantly, the rapid evolution of mammalian HERC6 spacer domain suggests that it might be a host-pathogen interface, targeting viral proteins and/or being the target of virus antagonists. Finally, we identified a HERC5/6 chimeric gene that arose from independent duplication in rodent and bat lineages and encodes for a conserved HERC5 N-terminal domain and divergent HERC6 spacer and HECT domains. This duplicated chimeric gene highlights adaptations that potentially contribute to rodent and bat immunity. Our findings open new research avenues on the functions of HERC6 and HERC5/6 in mammals, and on their implication in antiviral innate immunity.

A Gliriform Tooth from the Eocene of the Erlian Basin (Nei Mongol, China) and the Premolar Morphology of Anagalidan Mammals at a Crossroads

The middle Eocene in Nei Mongol (China) was an interval of profound faunal changes as regards the basal Glires and gliriform mammals in general. A major diversification of rodent lineages (ctenodactyloids) and more modern small-sized lagomorphs was accompanied by a decline of mimotonids (Gomphos and Mimolagus) and anagalids. The latter was an enigmatic group of basal Euarchontoglires endemic to China and Mongolia. Here, we describe the first anagalid tooth (a P4) from the Huheboerhe classic site in the Erlian Basin. The tooth, characterized by its unique morphology intermediate between mimotonids and anagalids is semihypsodont, has a single buccal root typical of mimotonids, a large paracone located anteriorly, and a nascent hypocone, characteristic of advanced anagalids. The new finding of neither an abundant nor speciose group suggests a greater diversity of anagalids in the Eocene of China. This discovery is important because it demonstrates the convergent adaptations in anagalids, possibly of ecological significance.

Functional signatures of evolutionarily young CTCF binding sites

BackgroundThe introduction of novel CTCF binding sites in gene regulatory regions in the rodent lineage is partly the effect of transposable element expansion, particularly in the murine lineage. The exact mechanism and functional impact of evolutionarily novel CTCF binding sites are not yet fully understood. We investigated the impact of novel subspecies-specific CTCF binding sites in two Mus genus subspecies, Mus musculus domesticus and Mus musculus castaneus, that diverged 0.5 million years ago.ResultsCTCF binding site evolution is influenced by the action of the B2-B4 family of transposable elements independently in both lineages, leading to the proliferation of novel CTCF binding sites. A subset of evolutionarily young sites may harbour transcriptional functionality as evidenced by the stability of their binding across multiple tissues in M. musculus domesticus (BL6), while overall the distance of subspecies-specific CTCF binding to the nearest transcription start sites and/or topologically associated domains (TADs) is largely similar to musculus-common CTCF sites. Remarkably, we discovered a recurrent regulatory architecture consisting of a CTCF binding site and an interferon gene that appears to have been tandemly duplicated to create a 15-gene cluster on chromosome 4, thus forming a novel BL6 specific immune locus in which CTCF may play a regulatory role.ConclusionsOur results demonstrate that thousands of CTCF binding sites show multiple functional signatures rapidly after incorporation into the genome.