Mammal madness: is the mammal tree of life not yet resolved?

Genomics, the diversification of mammals, and the evolution of placentation

When and why did variations in placental structure and function evolve? Such questions cannot be addressed without a reliable version of mammalian phylogeny. Twenty-five years ago, the mammalian tree was reshaped by molecular phylogenetics. Soon it was shown, in contrast to prevailing theories, that the common ancestor of placental mammals had invasive placentation. Subsequently, evolution of many other features of extraembryonic membranes was addressed. This endeavour stimulated research to fill gaps in our knowledge of placental morphology. Last year the mammalian tree was again revised based on a large set of genomic data. With that in mind, this review provides an update on placentation in the nineteen orders of placental mammals, incorporating much recent data. The principal features such as shape, interdigitation, the interhaemal barrier and the yolk sac are summarized in synoptic tables. The evolution of placental traits and its timing is then explored by reference to the revised mammalian tree. Examples are the early appearance of epitheliochorial placentation in the common ancestor of artiodactyls, perissodactyls, pangolins and carnivores (with reversion to invasive forms in the latter) and later refinements such as the binucleate trophoblast cells and placentomes of ruminants. In primates, the intervillous space gradually evolved from the more basic labyrinth whereas trophoblast invasion of the decidua was a late development in humans and great apes. Only seldom can we glimpse the "why" of placental evolution. The best examples concern placental hormones, including some striking examples of convergent evolution such as the chorionic gonadotropins of primates and equids. In concluding, I review current ideas about what drives placental evolution and identify significant gaps in our knowledge of placentation, including several relevant to the evolution of placentation in primates.

Chromosome-level genome and population genomics of the intermediate horseshoe bat ( Rhinolophus affinis) reveal the molecular basis of virus tolerance in Rhinolophus and echolocation call frequency variation

Horseshoe bats (genus Rhinolophus, family Rhinolophidae) represent an important group within chiropteran phylogeny due to their distinctive traits, including constant high-frequency echolocation, rapid karyotype evolution, and unique immune system. Advances in evolutionary biology, supported by high-quality reference genomes and comprehensive whole-genome data, have significantly enhanced our understanding of species origins, speciation mechanisms, adaptive evolutionary processes, and phenotypic diversity. However, genomic research and understanding of the evolutionary patterns of Rhinolophus are severely constrained by limited data, with only a single published genome of R. ferrumequinum currently available. In this study, we constructed a high-quality chromosome-level reference genome for the intermediate horseshoe bat ( R. affinis). Comparative genomic analyses revealed potential genetic characteristics associated with virus tolerance in Rhinolophidae. Notably, we observed expansions in several immune-related gene families and identified various genes functionally associated with the SARS-CoV-2 signaling pathway, DNA repair, and apoptosis, which displayed signs of rapid evolution. In addition, we observed an expansion of the major histocompatibility complex class II (MHC-II) region and a higher copy number of the HLA- DQB2 gene in horseshoe bats compared to other chiropteran species. Based on whole-genome resequencing and population genomic analyses, we identified multiple candidate loci (e.g., GLI3) associated with variations in echolocation call frequency across R. affinis subspecies. This research not only expands our understanding of the genetic characteristics of the Rhinolophus genus but also establishes a valuable foundation for future research.

Phylogenomics of Afrotherian mammals and improved resolution of extant Paenungulata

Molecular investigations have gathered a diverse set of mammals-predominantly African natives like elephants, hyraxes, and aardvarks-into a clade known as Afrotheria. Nevertheless, the precise phylogenetic relationships among these species remain contentious. Here, we sourced orthologous markers and ultraconserved elements to discern the interordinal connections among Afrotherian mammals. Our phylogenetic analyses bolster the common origin of Afroinsectiphilia and Paenungulata, and propose Afrosoricida as the closer relative to Macroscelidea rather than Tubulidentata, while also challenging the notion of Sirenia and Hyracoidea as sister taxa. The approximately unbiased test and the gene concordance factor uniformly recognized the alliance of Proboscidea with Hyracoidea as the dominant topology within Paenungulata. Investigation into sites with extremly high phylogenetic signal unveiled their potential to intensify conflicts in the Paenungulata topology. Subsequent exploration suggested that incomplete lineage sorting was predominantly responsible for the observed contentious relationships, whereas introgression exerted a subsidiary influence. The divergence times estimated in our study hint at the Cretaceous-Paleogene (K-Pg) extinction event as a catalyst for Afrotherian diversification. Overall, our findings deliver a tentative but insightful overview of Afrotheria phylogeny and divergence, elucidating these relationships through the lens of phylogenomics.

SMA20/PMIS2 Is a Rapidly Evolving Sperm Membrane Alloantigen with Possible Species-Divergent Function in Fertilization

Immunodominant alloantigens in pig sperm membranes include 15 known gene products and a previously undiscovered Mr 20,000 sperm membrane-specific protein (SMA20). Here we characterize SMA20 and identify it as the unannotated pig ortholog of PMIS2. A composite SMA20 cDNA encoded a 126 amino acid polypeptide comprising two predicted transmembrane segments and an N-terminal alanine- and proline (AP)-rich region with no apparent signal peptide. The Northern blots showed that the composite SMA20 cDNA was derived from a 1.1 kb testis-specific transcript. A BLASTp search retrieved no SMA20 match from the pig genome, but it did retrieve a 99% match to the Pmis2 gene product in warthog. Sequence identity to predicted PMIS2 orthologs from other placental mammals ranged from no more than 80% overall in Cetartiodactyla to less than 60% in Primates, with the AP-rich region showing the highest divergence, including, in the extreme, its absence in most rodents, including the mouse. SMA20 immunoreactivity localized to the acrosome/apical head of methanol-fixed boar spermatozoa but not live, motile cells. Ultrastructurally, the SMA20 AP-rich domain immunolocalized to the inner leaflet of the plasma membrane, the outer acrosomal membrane, and the acrosomal contents of ejaculated spermatozoa. Gene name search failed to retrieve annotated Pmis2 from most mammalian genomes. Nevertheless, individual pairwise interrogation of loci spanning Atp4a-Haus5 identified Pmis2 in all placental mammals, but not in marsupials or monotremes. We conclude that the gene encoding sperm-specific SMA20/PMIS2 arose de novo in Eutheria after divergence from Metatheria, whereupon rapid molecular evolution likely drove the acquisition of a species-divergent function unique to fertilization in placental mammals.

Translating the Timing of Developmental Benchmarks in Short-Tailed Opossums (Monodelphisdomestica) to Facilitate Comparisons with Commonly Used Rodent Models

INTRODUCTION

The gray short-tailed opossum, Monodelhis domestica (M. domestica), is a widely used marsupial model species that presents unique advantages for neurodevelopmental studies. Notably their extremely altricial birth allows manipulation of postnatal pups at timepoints equivalent to embryonic stages of placental mammals. A robust literature exists on the development of short-tailed opossums, but many researchers working in the more conventional model species of mice and rats may find it daunting to identify the appropriate age at which to conduct experiments.

METHODS

Here, we present detailed staging diagrams taken from photographic observations of 40 individual pups, in 6 litters, over 25 timepoints across postnatal development. We also present a comparative neurodevelopmental timeline of short-tailed opossums (M. domestica), the house mouse (Mus musculus), and the laboratory rat (Rattus norvegicus) during embryonic as well as postnatal development, using timepoints taken from this study and a review of existing literature, and use this dataset to present statistical models comparing the opossum to the rat and mouse.

RESULTS

One aim of this research was to aid in testing the generalizability of results found in rodents to other mammalian brains, such as the more distantly related metatherians. However, this broad dataset also allows the identification of potential heterochronies in opossum development compared to rats and mice. In contrast to previous work, we found broad similarity between the pace of opossum neural development with that of rats and mice. We also found that development of some systems was accelerated in the opossum, such as the forelimb motor plant, oral motor control, and some aspects of the olfactory system, while the development of the cortex, some aspects of the retina, and other aspects of the olfactory system are delayed compared to the rat and mouse.

DISCUSSION

The pace of opossum development is broadly similar to that of mice and rats, which underscores the usefulness of this species as a compliment to the more commonly used rodents. Many features that differ the most between opossums and rats and mice were either clustered around the day of birth and were features that have functional importance for the pup immediately after or during birth, or were features that have reduced functional importance for the pup until later in postnatal development, given that it is initially attached to the mother.

The evolution of menopause in toothed whales

Understanding how and why menopause has evolved is a long-standing challenge across disciplines. Females can typically maximize their reproductive success by reproducing for the whole of their adult life. In humans, however, women cease reproduction several decades before the end of their natural lifespan1,2. Although progress has been made in understanding the adaptive value of menopause in humans3,4, the generality of these findings remains unclear. Toothed whales are the only mammal taxon in which menopause has evolved several times5, providing a unique opportunity to test the theories of how and why menopause evolves in a comparative context. Here, we assemble and analyse a comparative database to test competing evolutionary hypotheses. We find that menopause evolved in toothed whales by females extending their lifespan without increasing their reproductive lifespan, as predicted by the 'live-long' hypotheses. We further show that menopause results in females increasing their opportunity for intergenerational help by increasing their lifespan overlap with their grandoffspring and offspring without increasing their reproductive overlap with their daughters. Our results provide an informative comparison for the evolution of human life history and demonstrate that the same pathway that led to menopause in humans can also explain the evolution of menopause in toothed whales.

Comparative Genome Microsynteny Illuminates the Fast Evolution of Nuclear Mitochondrial Segments (NUMTs) in Mammals

The escape of DNA from mitochondria into the nuclear genome (nuclear mitochondrial DNA, NUMT) is an ongoing process. Although pervasively observed in eukaryotic genomes, their evolutionary trajectories in a mammal-wide context are poorly understood. The main challenge lies in the orthology assignment of NUMTs across species due to their fast evolution and chromosomal rearrangements over the past 200 million years. To address this issue, we systematically investigated the characteristics of NUMT insertions in 45 mammalian genomes and established a novel, synteny-based method to accurately predict orthologous NUMTs and ascertain their evolution across mammals. With a series of comparative analyses across taxa, we revealed that NUMTs may originate from nonrandom regions in mtDNA, are likely found in transposon-rich and intergenic regions, and unlikely code for functional proteins. Using our synteny-based approach, we leveraged 630 pairwise comparisons of genome-wide microsynteny and predicted the NUMT orthology relationships across 36 mammals. With the phylogenetic patterns of NUMT presence-and-absence across taxa, we constructed the ancestral state of NUMTs given the mammal tree using a coalescent method. We found support on the ancestral node of Fereuungulata within Laurasiatheria, whose subordinal relationships are still controversial. This study broadens our knowledge on NUMT insertion and evolution in mammalian genomes and highlights the merit of NUMTs as alternative genetic markers in phylogenetic inference.

Geologically calibrated mammalian tree and its correlation with global events, including the emergence of humans

A robust timetree for Mammalia was constructed using the time calibration function of BEAST v1.10.4 and MEGA 11. The analysis involved the application of times of the most recent common ancestors, including a total of 19 mammalian fossil calibration ages following Benton et al. (Palaeontologia Electronica, 2015, 1-106) for their minimum ages. Additionally, fossil calibration ages for Gorilla, Pan, and a geologic event calibration age for otters were incorporated. Using these calibration ages, I constructed a geologically calibrated tree that estimates the age of the Homo and Pan splitting to be 5.69 Ma. The tree carries several significant implications. First, after the initial rifting at 120 Ma, the Atlantic Ocean expanded by over 500 km around Chron 34 (84 Ma), and vicariant speciation between Afrotheria (Africa) and Xenarthra (South America) appears to have commenced around 70 Ma. Additionally, ordinal level differentiations began immediately following the K-Pg boundary (66.0 Ma), supporting previous hypothesis that mammalian radiation rapidly filled ecological niches left vacant by non-avian dinosaurs. I constructed a diagram depicting the relationship between base substitution rate and age using an additional function in BEAST v1.10.4. The diagram reveals an exponential increase in the base substitution rate approaching recent times. This increased base substitution rate during the Neogene period may have contributed to the expansion of biodiversity, including the extensive adaptive radiation that led to the evolution of Homo sapiens. One significant driving factor behind this radiation could be attributed to the emergence and proliferation of C4 grasses since 20 Ma. These grasses have played a role in increasing carbon fixation, reducing atmospheric CO2 concentration, inducing global cooling, and initiating Quaternary glacial-interglacial cycles, thereby causing significant climatic changes.

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.

Convergent evolution of dim light vision in owls and deep-diving whales

Animals with enhanced dim-light sensitivity are at higher risk of light-induced retinal degeneration when exposed to bright light conditions.1,2,3,4 This trade-off is mediated by the rod photoreceptor sensory protein, rhodopsin (RHO), and its toxic vitamin A chromophore by-product, all-trans retinal.5,6,7,8 Rod arrestin (Arr-1) binds to RHO and promotes sequestration of excess all-trans retinal,9,10 which has recently been suggested as a protective mechanism against photoreceptor cell death.2,11 We investigated Arr-1 evolution in animals at high risk of retinal damage due to periodic bright-light exposure of rod-dominated retinas. Here, we find the convergent evolution of enhanced Arr-1/RHO all-trans-retinal sequestration in owls and deep-diving whales. Statistical analyses reveal a parallel acceleration of Arr-1 evolutionary rates in these lineages, which is associated with the introduction of a rare Arr-1 mutation (Q69R) into the RHO-Arr-1 binding interface. Using in vitro assays, we find that this single mutation significantly enhances RHO-all-trans-retinal sequestration by ∼30%. This functional convergence across 300 million years of evolutionary divergence suggests that Arr-1 and RHO may play an underappreciated role in the photoprotection of the eye, with potentially vast clinical significance.

Homoplasy of Retrotransposon Insertions in Toothed Whales

Retrotransposon insertion patterns facilitate a virtually homoplasy-free picture of phylogenetic history. Still, a few most likely random parallel insertions or deletions result in rare cases of homoplasy in primates. The following question arises: how frequent is retrotransposon homoplasy in other phylogenetic clades? Here, we derived genome insertion data of toothed whales to evaluate the extension of homoplasy in a representative laurasiatherian group. Among more than a thousand extracted and aligned retrotransposon loci, we detected 37 cases of precise parallel insertions in species that are separated by over more than 10 million years, a time frame which minimizes the effects of incomplete lineage sorting. We compared the phylogenetic signal of insertions with the flanking sequences of these loci to further exclude potential polymorphic loci derived by incomplete lineage sorting. We found that the phylogenetic signals of retrotransposon insertion patterns exhibiting true homoplasy differ from the signals of their flanking sequences. In toothed whales, precise parallel insertions account for around 0.18-0.29% of insertion cases, which is about 12.5 times the frequency of such insertions among Alus in primates. We also detected five specific deletions of retrotransposons on various lineages of toothed whale evolution, a frequency of 0.003%, which is slightly higher than such occurrences in primates. Overall, the level of retrotransposon homoplasy in toothed whales is still marginal compared to the phylogenetic diagnostic retrotransposon presence/absence signal.

Comparative analysis of the myoglobin gene in whales and humans reveals evolutionary changes in regulatory elements and expression levels

Cetacea and other diving mammals have undergone numerous adaptations to their aquatic environment, among them high levels of the oxygen-carrying intracellular hemoprotein myoglobin in skeletal muscles. Hypotheses regarding the mechanisms leading to these high myoglobin levels often invoke the induction of gene expression by exercise, hypoxia, and other physiological gene regulatory pathways. Here we explore an alternative hypothesis: that cetacean myoglobin genes have evolved high levels of transcription driven by the intrinsic developmental mechanisms that drive muscle cell differentiation. We have used luciferase assays in differentiated C2C12 cells to test this hypothesis. Contrary to our hypothesis, we find that the myoglobin gene from the minke whale, Balaenoptera acutorostrata, shows a low level of expression, only about 8% that of humans. This low expression level is broadly shared among cetaceans and artiodactylans. Previous work on regulation of the human gene has identified a core muscle-specific enhancer comprised of two regions, the "AT element" and a C-rich sequence 5' of the AT element termed the "CCAC-box". Analysis of the minke whale gene supports the importance of the AT element, but the minke whale CCAC-box ortholog has little effect. Instead, critical positive input has been identified in a G-rich region 3' of the AT element. Also, a conserved E-box in exon 1 positively affects expression, despite having been assigned a repressive role in the human gene. Last, a novel region 5' of the core enhancer has been identified, which we hypothesize may function as a boundary element. These results illustrate regulatory flexibility during evolution. We discuss the possibility that low transcription levels are actually beneficial, and that evolution of the myoglobin protein toward enhanced stability is a critical factor in the accumulation of high myoglobin levels in adult cetacean muscle tissue.

A timescale for placental mammal diversification based on Bayesian modeling of the fossil record

The timing of the placental mammal radiation has been the focus of debate over the efficacy of competing methods for establishing evolutionary timescales. Molecular clock analyses estimate that placental mammals originated before the Cretaceous-Paleogene (K-Pg) mass extinction, anywhere from the Late Cretaceous to the Jurassic. However, the absence of definitive fossils of placentals before the K-Pg boundary is compatible with a post-Cretaceous origin. Nevertheless, lineage divergence must occur before it can be manifest phenotypically in descendent lineages. This, combined with the non-uniformity of the rock and fossil records, requires the fossil record to be interpreted rather than read literally. To achieve this, we introduce an extended Bayesian Brownian bridge model that estimates the age of origination and, where applicable, extinction through a probabilistic interpretation of the fossil record. The model estimates the origination of placentals in the Late Cretaceous, with ordinal crown groups originating at or after the K-Pg boundary. The results reduce the plausible interval for placental mammal origination to the younger range of molecular clock estimates. Our findings support both the Long Fuse and Soft Explosive models of placental mammal diversification, indicating that the placentals originated shortly prior to the K-Pg mass extinction. The origination of many modern mammal lineages overlapped with and followed the K-Pg mass extinction.

Evolutionary constraint and innovation across hundreds of placental mammals

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Integrating gene annotation with orthology inference at scale

Annotating coding genes and inferring orthologs are two classical challenges in genomics and evolutionary biology that have traditionally been approached separately, limiting scalability. We present TOGA (Tool to infer Orthologs from Genome Alignments), a method that integrates structural gene annotation and orthology inference. TOGA implements a different paradigm to infer orthologous loci, improves ortholog detection and annotation of conserved genes compared with state-of-the-art methods, and handles even highly fragmented assemblies. TOGA scales to hundreds of genomes, which we demonstrate by applying it to 488 placental mammal and 501 bird assemblies, creating the largest comparative gene resources so far. Additionally, TOGA detects gene losses, enables selection screens, and automatically provides a superior measure of mammalian genome quality. TOGA is a powerful and scalable method to annotate and compare genes in the genomic era.

A genomic timescale for placental mammal evolution

The precise pattern and timing of speciation events that gave rise to all living placental mammals remain controversial. We provide a comprehensive phylogenetic analysis of genetic variation across an alignment of 241 placental mammal genome assemblies, addressing prior concerns regarding limited genomic sampling across species. We compared neutral genome-wide phylogenomic signals using concatenation and coalescent-based approaches, interrogated phylogenetic variation across chromosomes, and analyzed extensive catalogs of structural variants. Interordinal relationships exhibit relatively low rates of phylogenomic conflict across diverse datasets and analytical methods. Conversely, X-chromosome versus autosome conflicts characterize multiple independent clades that radiated during the Cenozoic. Genomic time trees reveal an accumulation of cladogenic events before and immediately after the Cretaceous-Paleogene (K-Pg) boundary, implying important roles for Cretaceous continental vicariance and the K-Pg extinction in the placental radiation.

NGSEP 4: Efficient and accurate identification of orthogroups and whole-genome alignment

Whole-genome alignment allows researchers to understand the genomic structure and variation among genomes. Approaches based on direct pairwise comparisons of DNA sequences require large computational capacities. As a consequence, pipelines combining tools for orthologous gene identification and synteny have been developed. In this manuscript, we present the latest functionalities implemented in NGSEP 4, to identify orthogroups and perform whole genome alignments. NGSEP implements functionalities for identification of clusters of homologus genes, synteny analysis and whole genome alignment. Our results showed that the NGSEP algorithm for orthogroups identification has competitive accuracy and efficiency in comparison to commonly used tools. The implementation also includes a visualization of the whole genome alignment based on synteny of the orthogroups that were identified, and a reconstruction of the pangenome based on frequencies of the orthogroups among the genomes. NGSEP 4 also includes a new graphical user interface based on the JavaFX technology. We expect that these new developments will be very useful for several studies in evolutionary biology and population genomics.

Cholinergic, catecholaminergic, serotonergic, and orexinergic neuronal populations in the brain of the lesser hedgehog tenrec (Echinops telfairi)

The current study provides an analysis of the cholinergic, catecholaminergic, serotonergic, and orexinergic neuronal populations, or nuclei, in the brain of the lesser hedgehog tenrec, as revealed with immunohistochemical techniques. For all four of these neuromodulatory systems, the nuclear organization was very similar to that observed in other Afrotherian species and is broadly similar to that observed in other mammals. The cholinergic system shows the most variation, with the lesser hedgehog tenrec exhibiting palely immunopositive cholinergic neurons in the ventral portion of the lateral septal nucleus, and the possible absence of cholinergic neurons in the parabigeminal nucleus and the medullary tegmental field. The nuclear complement of the catecholaminergic, serotonergic and orexinergic systems showed no specific variances in the lesser hedgehog tenrec when compared to other Afrotherians, or broadly with other mammals. A striking feature of the lesser hedgehog tenrec brain is a significant mesencephalic flexure that is observed in most members of the Tenrecoidea, as well as the closely related Chrysochlorinae (golden moles), but is not present in the greater otter shrew, a species of the Potomogalidae lineage currently incorporated into the Tenrecoidea. In addition, the cholinergic neurons of the ventral portion of the lateral septal nucleus are observed in the golden moles, but not in the greater otter shrew. This indicates that either complex parallel evolution of these features occurred in the Tenrecoidea and Chrysochlorinae lineages, or that the placement of the Potomogalidae within the Tenrecoidea needs to be re-examined.

The Biology and Evolution of Fierce Females (Moles and Hyenas)

Talpid moles and spotted hyenas have become the paradigms of anatomical and behavioral female masculinization. Females of many mole species develop ovotestes that produce testosterone, show external genitalia that resemble that of males, and close their vaginal orifice after every estrus, and female spotted hyenas lack an external vaginal orifice and develop a pseudoscrotum and a large pseudopenis through which they urinate, mate, and give birth. We review current knowledge about several significant aspects of the biology and evolution of these females, including (a) their specific study methods; (b) their unique anatomical features, and how these peculiarities influence certain physiological functions; and (c) the role that steroid hormones as well as genetic and environmental factors may have in urogenital system development, aggressive behavior, and social dominance. Nevertheless, both mole and hyena females are exceptionally efficient mothers, so their peculiar genitalia should not call into question their femininity.

The current state of carnivore cognition

The field of animal cognition has advanced rapidly in the last 25 years. Through careful and creative studies of animals in captivity and in the wild, we have gained critical insights into the evolution of intelligence, the cognitive capacities of a diverse array of taxa, and the importance of ecological and social environments, as well as individual variation, in the expression of cognitive abilities. The field of animal cognition, however, is still being influenced by some historical tendencies. For example, primates and birds are still the majority of study species in the field of animal cognition. Studies of diverse taxa improve the generalizability of our results, are critical for testing evolutionary hypotheses, and open new paths for understanding cognition in species with vastly different morphologies. In this paper, we review the current state of knowledge of cognition in mammalian carnivores. We discuss the advantages of studying cognition in Carnivorans and the immense progress that has been made across many cognitive domains in both lab and field studies of carnivores. We also discuss the current constraints that are associated with studying carnivores. Finally, we explore new directions for future research in studies of carnivore cognition.

Attenuated evolution of mammals through the Cenozoic

The Cenozoic diversification of placental mammals is the archetypal adaptive radiation. Yet, discrepancies between molecular divergence estimates and the fossil record fuel ongoing debate around the timing, tempo, and drivers of this radiation. Analysis of a three-dimensional skull dataset for living and extinct placental mammals demonstrates that evolutionary rates peak early and attenuate quickly. This long-term decline in tempo is punctuated by bursts of innovation that decreased in amplitude over the past 66 million years. Social, precocial, aquatic, and herbivorous species evolve fastest, especially whales, elephants, sirenians, and extinct ungulates. Slow rates in rodents and bats indicate dissociation of taxonomic and morphological diversification. Frustratingly, highly similar ancestral shape estimates for placental mammal superorders suggest that their earliest representatives may continue to elude unequivocal identification.

The Inference of the Evolution of Immune Traits as Constrained by Phylogeny: Insight into the Immune System of the Basal Diapsid

Simple Summary

In light of emerging pathogenic threats affecting wildlife, it is important to broaden the current understanding of immune system function, development, and evolution. The relation of descent of immune traits is critical to understand the ability of organisms to handle pathogens. Here, we explore the evolution of toll-like receptors (TLRs), a series of receptors crucial to the initial immune response in reptiles. Our analysis revealed that the common ancestor may have had an immune system that lacked two receptors: TLR 15, a receptor uniquely present in Reptilia, and TLR 13, a receptor important in the recognition of pathogens. Additionally, our analysis showed a dynamic evolution for various TLRs, likely attributed to redundancies in function.

Abstract

Among vertebrates, some of the most vulnerable taxa to emergent fungal pathogens are members of Reptilia. In light of the growing threat of emergent fungal pathogens affecting wildlife, it is important to broaden the current understanding of immune system function, development, and evolution. The homologous condition of a trait is necessary in order to study its evolution, as such, homology is necessary in the study of immunological evolution. Here, we explore the evolution of toll-like receptors (TLRs), a series of homologous receptors crucial to the initial immune response. The homologous condition of TLR genes provides a unique system in which to explore the evolution of the TLR; using a Reptilian phylogeny, we elucidate the immune condition of the basal diapsid. Our analysis revealed that the basal diapsid may have had an immune system that lacked two receptors: TLR 15, a receptor uniquely present in Reptilia, and TLR 13, a receptor important in the recognition of nucleic acid motifs. Additionally, our analysis showed multiple losses and convergences for various TLRs, likely attributed to redundancies in receptor function. Further exploration into the immune condition of extinct taxa may shed light on the evolution of the reptilian immune system.

Comparative transcriptome analysis reveals molecular adaptations underlying distinct immunity and inverted resting posture in bats

Understanding how natural selection shapes unique traits in mammals is a central topic in evolutionary biology. The mammalian order Chiroptera (bats) is attractive for biologists as well as the general public due to their specific traits of extraordinary immunity and inverted resting posture. However, genomic resources for bats that occupy key phylogenetic positions are not sufficient, which hinders comprehensive investigation of the molecular mechanisms underpinning the origin of specific traits in bats. Here, we sequenced the transcriptomes of five bats that are phylogenetically divergent and occupy key positions in the phylogenetic tree of bats. In combination with the available genomes of 19 bats and 21 other mammals, we built a database consisting of 10,918 one-to-one ortholog genes and reconstructed phylogenetic relationships of these mammals. We found that genes related to immunity, bone remodeling and cardiovascular system are targets of natural selection along the ancestral branch of bats. Further analyses revealed that the T cell receptor signaling pathway involved in immune adaptation is specifically enriched in bats. Moreover, molecular adaptations of bone remodeling, cardiovascular system, and balance sensing may help to explain the reverted resting posture in bats. Our study provides valuable transcriptome resources, enabling us to tentatively identify genetic changes associated with bat-specific traits. This work is among the first to advance our understanding of molecular underpinnings of inverted resting posture in bats, which could provide insight into healthcare applications such as hypertension in humans. This article is protected by copyright. All rights reserved.

Rapid divergence of a gamete recognition gene promoted macroevolution of Eutheria

BACKGROUND

Speciation genes contribute disproportionately to species divergence, but few examples exist, especially in vertebrates. Here we test whether Zan, which encodes the sperm acrosomal protein zonadhesin that mediates species-specific adhesion to the egg's zona pellucida, is a speciation gene in placental mammals.

RESULTS

Genomic ontogeny reveals that Zan arose by repurposing of a stem vertebrate gene that was lost in multiple lineages but retained in Eutheria on acquiring a function in egg recognition. A 112-species Zan sequence phylogeny, representing 17 of 19 placental Orders, resolves all species into monophyletic groups corresponding to recognized Orders and Suborders, with <5% unsupported nodes. Three other rapidly evolving germ cell genes (Adam2, Zp2, and Prm1), a paralogous somatic cell gene (TectA), and a mitochondrial gene commonly used for phylogenetic analyses (Cytb) all yield trees with poorer resolution than the Zan tree and inferior topologies relative to a widely accepted mammalian supertree. Zan divergence by intense positive selection produces dramatic species differences in the protein's properties, with ordinal divergence rates generally reflecting species richness of placental Orders consistent with expectations for a speciation gene that acts across a wide range of taxa. Furthermore, Zan's combined phylogenetic utility and divergence exceeds those of all other genes known to have evolved in Eutheria by positive selection, including the only other mammalian speciation gene, Prdm9.

CONCLUSIONS

Species-specific egg recognition conferred by Zan's functional divergence served as a mode of prezygotic reproductive isolation that promoted the extraordinary adaptive radiation and success of Eutheria.

Ancient whale rhodopsin reconstructs dim-light vision over a major evolutionary transition: Implications for ancestral diving behavior

Cetaceans are fully aquatic mammals that descended from terrestrial ancestors, an iconic evolutionary transition characterized by adaptations for underwater foraging via breath-hold diving. Although the evolutionary history of this specialized behavior is challenging to reconstruct, coevolving sensory systems may offer valuable clues. The dim-light visual pigment, rhodopsin, which initiates phototransduction in the rod photoreceptors of the eye, has provided insight into the visual ecology of depth in several aquatic vertebrate lineages. Here, we use ancestral sequence reconstruction and protein resurrection experiments to quantify light-activation metrics in rhodopsin pigments from ancestors bracketing the cetacean terrestrial-to-aquatic transition. By comparing multiple reconstruction methods on a broadly sampled cetartiodactyl species tree, we generated highly robust ancestral sequence estimates. Our experimental results provide direct support for a blue-shift in spectral sensitivity along the branch separating cetaceans from terrestrial relatives. This blue-shift was 14 nm, resulting in a deep-sea signature (λ_max = 486 nm) similar to many mesopelagic-dwelling fish. We also discovered that the decay rates of light-activated rhodopsin increased in ancestral cetaceans, which may indicate an accelerated dark adaptation response typical of deeper-diving mammals. Because slow decay rates are thought to help sequester cytotoxic photoproducts, this surprising result could reflect an ecological trade-off between rod photoprotection and dark adaptation. Taken together, these ancestral shifts in rhodopsin function suggest that some of the first fully aquatic cetaceans could dive into the mesopelagic zone (>200 m). Moreover, our reconstructions indicate that this behavior arose before the divergence of toothed and baleen whales.

The brain of the tree pangolin (Manis tricuspis). IX. The pallial telencephalon

A cyto‐, myelo‐, and chemoarchitectonic analysis of the pallial telencephalon of the tree pangolin is provided. As certain portions of the pallial telencephalon have been described previously (olfactory pallium, hippocampal formation, and amygdaloid complex), we focus on the claustrum and endopiriform nuclear complex, the white matter and white matter interstitial cells, and the areal organization of the cerebral cortex. Our analysis indicates that the organization of the pallial telencephalon of the tree pangolin is similar to that observed in many other mammals, and specifically quite similar to the closely related carnivores. The claustrum of the tree pangolin exhibits a combination of insular and laminar architecture, while the endopiriform nuclear complex contains three nuclei, both reminiscent of observations made in other mammals. The population of white matter interstitial cells resembles that observed in other mammals, while a distinct laminated organization of the intracortical white matter was revealed with parvalbumin immunostaining. The cerebral cortex of the tree pangolin presented with indistinct laminar boundaries as well as pyramidalization of the neurons in both layers 2 and 4. All cortical regions typically found in mammals were present, with the cortical areas within these regions often corresponding to what has been reported in carnivores. Given the similarity of the organization of the pallial telencephalon of the tree pangolin to that observed in other mammals, especially carnivores, it would be reasonable to assume that the neural processing afforded the tree pangolin by these structures does not differ dramatically to that of other mammals.

The brain of the tree pangolin (Manis tricuspis). VII. The amygdaloid body

Here, we describe the cytoarchitecture and chemoarchitecture of the amygdaloid body of the tree pangolin. Our definition of the amygdaloid body includes the pallial portions of the amygdala, and the centromedial group that is a derivative of the subpallium and part of the extended amygdala. The remainder of the extended amygdala is not described herein. Within the amygdaloid body of the tree pangolin, we identified the basolateral group (composed of the lateral, basal, and accessory basal amygdaloid nuclei), the superficial, or cortical nuclei (the anterior and posterior cortical nuclei, the periamygdaloid cortex, and nuclei of the olfactory tract), the centromedial group (the central amygdaloid nucleus and the medial nuclear cluster), and other amygdaloid nuclei (the anterior amygdaloid area, the amygdalohippocampal area, the intramedullary group, and intercalated islands). The location within and relative to each other within the amygdaloid body and the internal subdivisions of these groups were very similar to that reported in other mammalian species, with no clearly derived features specific to the tree pangolin. The only variation was the lack of an insular appearance of the intercalated islands, which in the tree pangolin were observed as a continuous band of neurons located dorsomedial to the basolateral group similar in appearance to and almost continuous with the intramedullary group. In carnivores, the closest relatives of the pangolins, and laboratory rats, a similar appearance of portions of the intercalated islands has been noted.

Contradictory Phylogenetic Signals in the Laurasiatheria Anomaly Zone

Relationships among laurasiatherian clades represent one of the most highly disputed topics in mammalian phylogeny. In this study, we attempt to disentangle laurasiatherian interordinal relationships using two independent genome-level approaches: (1) quantifying retrotransposon presence/absence patterns, and (2) comparisons of exon datasets at the levels of nucleotides and amino acids. The two approaches revealed contradictory phylogenetic signals, possibly due to a high level of ancestral incomplete lineage sorting. The positions of Eulipotyphla and Chiroptera as the first and second earliest divergences were consistent across the approaches. However, the phylogenetic relationships of Perissodactyla, Cetartiodactyla, and Ferae, were contradictory. While retrotransposon insertion analyses suggest a clade with Cetartiodactyla and Ferae, the exon dataset favoured Cetartiodactyla and Perissodactyla. Future analyses of hitherto unsampled laurasiatherian lineages and synergistic analyses of retrotransposon insertions, exon and conserved intron/intergenic sequences might unravel the conflicting patterns of relationships in this major mammalian clade.

Petrosal and bony labyrinth morphology of the stem paenungulate mammal (Paenungulatomorpha) Ocepeia daouiensis from the Paleocene of Morocco

Based on high-resolution computed tomography, we describe in detail the petrosal and inner ear anatomy of one of the few known African stem paenungulates (Paenungulatomorpha), Ocepeia daouiensis from the Selandian of the Ouled Abdoun phosphate basin (Morocco). The petrosal of Ocepeia displays some remarkable, probably derived features (among eutherians) such as relatively small pars cochlearis, pars canalicularis labyrinth (including small semicircular canals), a large wing-like pars mastoidea, a large and inflated tegmen tympani, and the dorsoventral orientation of the large canal for the ramus superior. The presence of small semicircular canals in Ocepeia is an interesting shared trait with tenrecoidean afrotherians. Otherwise, and consistent with a general primitive skull morphology, the middle ear and labyrinth of Ocepeia daouiensis is characterised by many plesiomorphic traits close to the eutherian generalised plan. This adds to the rather generalised morphology of the earliest crown paenungulates such as Eritherium, Phosphatherium and Seggeurius to support an ancestral paenungulatomorph morphotype poorly derived from the eutherian pattern. As a result, Ocepeia provides key morphological and fossil data to test phylogenetic relationships of the Afrotheria (including Paenungulatomorpha) at the placental root mostly inferred from molecular studies.

Champagne: Automated Whole-Genome Phylogenomic Character Matrix Method Using Large Genomic Indels for Homoplasy-Free Inference

We present Champagne, a whole-genome method for generating character matrices for phylogenomic analysis using large genomic indel events. By rigorously picking orthologous genes and locating large insertion and deletion events, Champagne delivers a character matrix that considerably reduces homoplasy compared with morphological and nucleotide-based matrices, on both established phylogenies and difficult-to-resolve nodes in the mammalian tree. Champagne provides ample evidence in the form of genomic structural variation to support incomplete lineage sorting and possible introgression in Paenungulata and human-chimp-gorilla which were previously inferred primarily through matrices composed of aligned single-nucleotide characters. Champagne also offers further evidence for Myomorpha as sister to Sciuridae and Hystricomorpha in the rodent tree. Champagne harbors distinct theoretical advantages as an automated method that produces nearly homoplasy-free character matrices on the whole-genome scale.

The performance of outgroup-free rooting under evolutionary radiations

Tree rooting implies a temporal dimension to phylogenies. Only after defining the position of the root node is that the ancestral-descendant relationship between branches can be fully deduced. Rooting has been usually carried out by employing evolutionarily close outgroup lineages, which is a drawback when these lineages are unavailable or unknown. Alternatively, outgroup-free rooting methods were proposed, which rely on the constancy of evolutionary rates to varying degrees. In this work we analyzed the performance of two of these methods, the midpoint rooting (MPR) and the minimal ancestor deviation (MAD), in rooting topologies evolved under challenging scenarios of fast evolutionary radiations derived from empirical data, characterized by short internal branches near the crown node. Considering all branch length combinations investigated, both methods exhibited average success rates below 50%, although MAD slightly outperformed MPR. Moreover, tree balance significantly impacted the relative performance of the methods. We found that, in four-taxa unrooted trees, the outcome of whether both methodologies will correctly root the tree can be roughly predicted by two simple dimensionless metrics: the coefficient of variation of the external branch lengths, and the ratio between the internal branch length to the total sum of branch lengths, which were employed to devise a general linear model that allowed calculating the probability of correct placing the root node for any four-taxa tree. We predicted that the performance of both outgroup-free rooting methods on loci representing the placental mammal radiation ranged between 50% and 75%.

Sperm bauplan and function and underlying processes of sperm formation and selection

The spermatozoon is a highly differentiated and polarized cell, with two main structures: the head, containing a haploid nucleus and the acrosomal exocytotic granule, and the flagellum, which generates energy and propels the cell; both structures are connected by the neck. The sperm's main aim is to participate in fertilization, thus activating development. Despite this common bauplan and function, there is an enormous diversity in structure and performance of sperm cells. For example, mammalian spermatozoa may exhibit several head patterns and overall sperm lengths ranging from ∼30 to 350 µm. Mechanisms of transport in the female tract, preparation for fertilization, and recognition of and interaction with the oocyte also show considerable variation. There has been much interest in understanding the origin of this diversity, both in evolutionary terms and in relation to mechanisms underlying sperm differentiation in the testis. Here, relationships between sperm bauplan and function are examined at two levels: first, by analyzing the selective forces that drive changes in sperm structure and physiology to understand the adaptive values of this variation and impact on male reproductive success and second, by examining cellular and molecular mechanisms of sperm formation in the testis that may explain how differentiation can give rise to such a wide array of sperm forms and functions.

Evolutionary pathways to SARS-CoV-2 resistance are opened and closed by epistasis acting on ACE2

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infects a broader range of mammalian species than previously predicted, binding a diversity of angiotensin converting enzyme 2 (ACE2) orthologs despite extensive sequence divergence. Within this sequence degeneracy, we identify a rare sequence combination capable of conferring SARS-CoV-2 resistance. We demonstrate that this sequence was likely unattainable during human evolution due to deleterious effects on ACE2 carboxypeptidase activity, which has vasodilatory and cardioprotective functions in vivo. Across the 25 ACE2 sites implicated in viral binding, we identify 6 amino acid substitutions unique to mouse-one of the only known mammalian species resistant to SARS-CoV-2. Substituting human variants at these positions is sufficient to confer binding of the SARS-CoV-2 S protein to mouse ACE2, facilitating cellular infection. Conversely, substituting mouse variants into either human or dog ACE2 abolishes viral binding, diminishing cellular infection. However, these same substitutions decrease human ACE2 activity by 50% and are predicted as pathogenic, consistent with the extreme rarity of human polymorphisms at these sites. This trade-off can be avoided, however, depending on genetic background; if substituted simultaneously, these same mutations have no deleterious effect on dog ACE2 nor that of the rodent ancestor estimated to exist 70 million years ago. This genetic contingency (epistasis) may have therefore opened the road to resistance for some species, while making humans susceptible to viruses that use these ACE2 surfaces for binding, as does SARS-CoV-2.

Evolutionary history of the vertebrate Piwi gene family

PIWIs are regulatory proteins that belong to the Argonaute family. Piwis are primarily expressed in gonads and protect the germline against the mobilization and propagation of transposable elements (TEs) through transcriptional gene silencing. Vertebrate genomes encode up to four Piwi genes: Piwil1, Piwil2, Piwil3 and Piwil4, but their duplication history is unresolved. We leveraged phylogenetics, synteny and expression analyses to address this void. Our phylogenetic analysis suggests Piwil1 and Piwil2 were retained in all vertebrate members. Piwil4 was the result of Piwil1 duplication in the ancestor of gnathostomes, but was independently lost in ray-finned fishes and birds. Further, Piwil3 was derived from a tandem Piwil1 duplication in the common ancestor of marsupial and placental mammals, but was secondarily lost in Atlantogenata (Xenarthra and Afrotheria) and some rodents. The evolutionary rate of Piwil3 is considerably faster than any Piwi among all lineages, but an explanation is lacking. Our expression analyses suggest Piwi expression has mostly been constrained to gonads throughout vertebrate evolution. Vertebrate evolution is marked by two early rounds of whole genome duplication and many multigene families are linked to these events. However, our analyses suggest Piwi expansion was independent of whole genome duplications.

The complete mitochondrial genome of Tylonycteris fulvida (Chiroptera: Vespertilionidae) in South China

The complete mitochondrial genome of Tylonycteris fulvida (Peters, 1872) was obtained using high-throughput sequencing technology. The genome is a circular molecule of 16,621 bp length, containing 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and a control region. A phylogenetic tree of 13 protein-coding genes was constructed using IQ-TREE. Our result suggests that T. fulvida cluster within Chiroptera and Fereuungulata. The complete mitochondrial genome sequence of T. fulvida will be helpful for future taxonomic and phylogenetic studies on Chiroptera.

Large-scale k-mer-based analysis of the informational properties of genomes, comparative genomics and taxonomy

Information theoretic approaches are ubiquitous and effective in a wide variety of bioinformatics applications. In comparative genomics, alignment-free methods, based on short DNA words, or k-mers, are particularly powerful. We evaluated the utility of varying k-mer lengths for genome comparisons by analyzing their sequence space coverage of 5805 genomes in the KEGG GENOME database. In subsequent analyses on four k-mer lengths spanning the relevant range (11, 21, 31, 41), hierarchical clustering of 1634 genus-level representative genomes using pairwise 21- and 31-mer Jaccard similarities best recapitulated a phylogenetic/taxonomic tree of life with clear boundaries for superkingdom domains and high subtree similarity for named taxons at lower levels (family through phylum). By analyzing ~14.2M prokaryotic genome comparisons by their lowest-common-ancestor taxon levels, we detected many potential misclassification errors in a curated database, further demonstrating the need for wide-scale adoption of quantitative taxonomic classifications based on whole-genome similarity.

Natural co-infection of divergent hepatitis B and C virus homologues in carnivores

In humans, co-infection of hepatitis B and C viruses (HBV, HCV) is common and aggravates disease outcome. Infection-mediated disease aggravation is poorly understood, partly due to lack of suitable animal models. Carnivores are understudied for hepatitis virus homologues. We investigated Mexican carnivores (ringtails, Bassariscus astutus) for HBV and HCV homologues. Three out of eight animals were infected with a divergent HBV termed ringtail HBV (RtHBV) at high viral loads of 5 × 10⁹ -1.4 × 10¹⁰ copies/ml serum. Two of the RtHBV-infected animals were co-infected with a divergent hepacivirus termed ringtail hepacivirus (RtHV) at 4 × 10⁶ -7.5 × 10⁷ copies/ml in strain-specific qRT-PCR assays. Immunofluorescence assays relying on HBV core and RtHV NS3/4a proteins indicated that none of the animals had detectable hepadnavirus core-specific antibodies, whereas one RtHV-infected animal had concomitant RtHV-specific antibodies at 1:800 end-point titre. RtHBV and RtHV complete genomes showed typical HBV and HCV structure and length. All RtHBV genomes were identical, whereas RtHV genomes showed four amino acid substitutions located predominantly in the E1/E2-encoding genomic regions. Both RtHBV (>28% genomic nucleotide sequence distance) and RtHV (>30% partial NS3/NS5B amino acid sequence distance) formed new species within their virus families. Evolutionary analyses showed that RtHBV grouped with HBV homologues from different laurasiatherian hosts (carnivores, bats, and ungulates), whereas RtHV grouped predominantly with rodent-borne viruses. Ancestral state reconstructions showed that RtHV, but not RtHBV, likely emerged via a non-recent host switch involving rodent-borne hepacivirus ancestors. Conserved hepatitis virus infection patterns in naturally infected ringtails indicate that carnivores may be promising animal models to understand HBV/HCV co-infection.

Nuclear organization of catecholaminergic neurons in the brains of a lar gibbon and a chimpanzee

Using tyrosine hydroxylase immunohistochemistry, we describe the nuclear parcellation of the catecholaminergic system in the brains of a lar gibbon (Hylobates lar) and a chimpanzee (Pan troglodytes). The parcellation of catecholaminergic nuclei in the brains of both apes is virtually identical to that observed in humans and shows very strong similarities to that observed in mammals more generally, particularly other primates. Specific variations of this system in the apes studied include an unusual high-density cluster of A10dc neurons, an enlarged retrorubral nucleus (A8), and an expanded distribution of the neurons forming the dorsolateral division of the locus coeruleus (A4). The additional A10dc neurons may improve dopaminergic modulation of the extended amygdala, the enlarged A8 nucleus may be related to the increased use of communicative facial expressions in the hominoids compared to other primates, while the expansion of the A4 nucleus appears to be related to accelerated evolution of the cerebellum in the hominoids compared to other primates. In addition, we report the presence of a compact division of the locus coeruleus proper (A6c), as seen in other primates, that is not present in other mammals apart from megachiropteran bats. The presence of this nucleus in primates and megachiropteran bats may reflect homology or homoplasy, depending on the evolutionary scenario adopted. The fact that the complement of homologous catecholaminergic nuclei is mostly consistent across mammals, including primates, is advantageous for the selection of model animals for the study of specific dysfunctions of the catecholaminergic system in humans.

Odorant-binding proteins in canine anal sac glands indicate an evolutionarily conserved role in mammalian chemical communication

Background

Chemical communication is an important aspect of the behavioural ecology of a wide range of mammals. In dogs and other carnivores, anal sac glands are thought to convey information to conspecifics by secreting a pallet of small volatile molecules produced by symbiotic bacteria. Because these glands are unique to carnivores, it is unclear how their secretions relate to those of other placental mammals that make use of different tissues and secretions for chemical communication. Here we analyse the anal sac glands of domestic dogs to verify the secretion of proteins and infer their evolutionary relationship to those involved in the chemical communication of non-carnivoran mammals.

Results

Proteomic analysis of anal sac gland secretions of 17 dogs revealed the consistently abundant presence of three related proteins. Homology searches against online databases indicate that these proteins are evolutionary related to ‘odorant binding proteins’ (OBPs) found in a wide range of mammalian secretions and known to contribute to chemical communication. Screening of the dog’s genome sequence show that the newly discovered OBPs are encoded by a single cluster of three genes in the pseudoautosomal region of the X-chromosome. Comparative genomic screening indicates that the same locus is shared by a wide range of placental mammals and that it originated at least before the radiation of extant placental orders. Phylogenetic analyses suggest a dynamic evolution of gene duplication and loss, resulting in large gene clusters in some placental taxa and recurrent loss of this locus in others. The homology of OBPs in canid anal sac glands and those found in other mammalian secretions implies that these proteins maintained a function in chemical communication throughout mammalian evolutionary history by multiple shifts in expression between secretory tissues involved in signal release and nasal mucosa involved in signal reception.

Conclusions

Our study elucidates a poorly understood part of the biology of a species that lives in close association with humans. In addition, it shows that the protein repertoire underlying chemical communication in mammals is more evolutionarily stable than the variation of involved glands and tissues would suggest.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01910-w.

The effects of fossil taxa, hypothetical predicted ancestors, and a molecular scaffold on pseudoextinction analyses of extant placental orders

Pseudoextinction analyses, which simulate extinction in extant taxa, use molecular phylogenetics to assess the accuracy of morphological phylogenetics. Previous pseudoextinction analyses have shown a failure of morphological phylogenetics to place some individual placental orders in the correct superordinal clade. Recent work suggests that the inclusion of hypothetical ancestors of extant placental clades, estimated by ancestral state reconstructions of morphological characters, may increase the accuracy of morphological phylogenetic analyses. However, these studies reconstructed direct hypothetical ancestors for each extant taxon based on a well-corroborated molecular phylogeny, which is not possible for extinct taxa that lack molecular data. It remains to be determined if pseudoextinct taxa, and by proxy extinct taxa, can be accurately placed when their immediate hypothetical ancestors are unknown. To investigate this, we employed molecular scaffolds with the largest available morphological data set for placental mammals. Each placental order was sequentially treated as pseudoextinct by exempting it from the molecular scaffold and recoding soft morphological characters as missing for all its constituent species. For each pseudoextinct data set, we omitted the pseudoextinct taxon and performed a parsimony ancestral state reconstruction to obtain hypothetical predicted ancestors. Each pseudoextinct order was then evaluated in seven parsimony analyses that employed combinations of fossil taxa, hypothetical predicted ancestors, and a molecular scaffold. In treatments that included fossils, hypothetical predicted ancestors, and a molecular scaffold, only 8 of 19 pseudoextinct placental orders (42%) retained the same interordinal placement as on the molecular scaffold. In treatments that included hypothetical predicted ancestors but not fossils or a scaffold, only four placental orders (21%) were recovered in positions that are congruent with the scaffold. These results indicate that hypothetical predicted ancestors do not increase the accuracy of pseudoextinct taxon placement when the immediate hypothetical ancestor of the taxon is unknown. Hypothetical predicted ancestors are not a panacea for morphological phylogenetics.

Phylogenetic Signal and Bias in Paleontology

An unprecedented amount of evidence now illuminates the phylogeny of living mammals and birds on the Tree of Life. We use this tree to measure phylogenetic value of data typically used in paleontology (bones and teeth) from six datasets derived from five published studies. We ask three interrelated questions: 1) Can these data adequately reconstruct known parts of the Tree of Life? 2) Is accuracy generally similar for studies using morphology, or do some morphological datasets perform better than others? 3) Does the loss of non-fossilizable data cause taxa to occur in misleadingly basal positions? Adding morphology to DNA datasets usually increases congruence of resulting topologies to the well corroborated tree, but this varies among morphological datasets. Extant taxa with a high proportion of missing morphological characters can greatly reduce phylogenetic resolution when analyzed together with fossils. Attempts to ameliorate this by deleting extant taxa missing morphology are prone to decreased accuracy due to long-branch artefacts. We find no evidence that fossilization causes extinct taxa to incorrectly appear at or near topologically basal branches. Morphology comprises the evidence held in common by living taxa and fossils, and phylogenetic analysis of fossils greatly benefits from inclusion of molecular and morphological data sampled for living taxa, whatever methods are used for phylogeny estimation.

Horizontal transfer and subsequent explosive expansion of a DNA transposon in sea kraits (Laticauda)

Transposable elements (TEs) are self-replicating genetic sequences and are often described as important 'drivers of evolution'. This driving force is because TEs promote genomic novelty by enabling rearrangement, and through exaptation as coding and regulatory elements. However, most TE insertions potentially lead to neutral or harmful outcomes, therefore host genomes have evolved machinery to suppress TE expansion. Through horizontal transposon transfer (HTT) TEs can colonize new genomes, and since new hosts may not be able to regulate subsequent replication, these TEs may proliferate rapidly. Here, we describe HTT of the Harbinger-Snek DNA transposon into sea kraits (Laticauda), and its subsequent explosive expansion within Laticauda genomes. This HTT occurred following the divergence of Laticauda from terrestrial Australian elapids approximately 15-25 Mya. This has resulted in numerous insertions into introns and regulatory regions, with some insertions into exons which appear to have altered UTRs or added sequence to coding exons. Harbinger-Snek has rapidly expanded to make up 8-12% of Laticauda spp. genomes; this is the fastest known expansion of TEs in amniotes following HTT. Genomic changes caused by this rapid expansion may have contributed to adaptation to the amphibious-marine habitat.

First Description of a Satellite DNA in Manatees' Centromeric Regions

Trichechus manatus and Trichechus inunguis are the two Sirenia species that occur in the Americas. Despite their increasing extinction risk, many aspects of their biology remain understudied, including the repetitive DNA fraction of their genomes. Here we used the sequenced genome of T. manatus and TAREAN to identify satellite DNAs (satDNAs) in this species. We report the first description of TMAsat, a satDNA comprising ~0.87% of the genome, with ~684bp monomers and centromeric localization. In T. inunguis, TMAsat showed similar monomer length, chromosome localization and conserved CENP-B box-like motifs as in T. manatus. We also detected this satDNA in the Dugong dugon and in the now extinct Hydrodamalis gigas genomes. The neighbor-joining tree shows that TMAsat sequences from T. manatus, T. inunguis, D. dugon, and H. gigas lack species-specific clusters, which disagrees with the predictions of concerted evolution. We detected a divergent TMAsat-like homologous sequence in elephants and hyraxes, but not in other mammals, suggesting this sequence was already present in the common ancestor of Paenungulata, and later became a satDNA in the Sirenians. This is the first description of a centromeric satDNA in manatees and will facilitate the inclusion of Sirenia in future studies of centromeres and satDNA biology.

Non-model systems in mammalian forelimb evo-devo

Mammal forelimbs are highly diverse, ranging from the elongated wing of a bat to the stout limb of the mole. The mammal forelimb has been a long-standing system for the study of early developmental patterning, proportional variation, shape change, and the reduction of elements. However, most of this work has been performed in mice, which neglects the wide variation present across mammal forelimbs. This review emphasizes the critical role of non-model systems in limb evo-devo and highlights new emerging models and their potential. We discuss the role of gene networks in limb evolution, and touch on functional analyses that lay the groundwork for further developmental studies. Mammal limb evo-devo is a rich field, and here we aim to synthesize the findings of key recent works and the questions to which they lead.

Phylogenomic systematics of the spotted skunks (Carnivora, Mephitidae, Spilogale): Additional species diversity and Pleistocene climate change as a major driver of diversification

Four species of spotted skunks (Carnivora, Mephitidae, Spilogale) are currently recognized: Spilogale angustifrons, S. gracilis, S. putorius, and S. pygmaea. Understanding species boundaries within this group is critical for effective conservation given that regional populations or subspecies (e.g., S. p. interrupta) have experienced significant population declines. Further, there may be currently unrecognized diversity within this genus as some taxa (e.g., S. angustifrons) and geographic regions (e.g., Central America) never have been assessed using DNA sequence data. We analyzed species limits and diversification patterns in spotted skunks using multilocus nuclear (ultraconserved elements) and mitochondrial (whole mitogenomes and single gene analysis) data sets from broad geographic sampling representing all currently recognized species and subspecies. We found a high degree of genetic divergence among Spilogale that reflects seven distinct species and eight unique mitochondrial lineages. Initial divergence between S. pygmaea and all other Spilogale occurred in the Early Pliocene (∼ 5.0 million years ago). Subsequent diversification of the remaining Spilogale into an "eastern" and a "western" lineage occurred during the Early Pleistocene (∼1.5 million years ago). These two lineages experienced temporally coincident patterns of diversification at ∼0.66 and ∼0.35 million years ago into two and ultimately three distinct evolutionary units, respectively. Diversification was confined almost entirely within the Pleistocene during a timeframe characterized by alternating glacial-interglacial cycles, with the origin of this diversity occurring in northeastern Mexico and the southwestern United States of America. Mitochondrial-nuclear discordance was recovered across three lineages in geographic regions consistent with secondary contact, including a distinct mitochondrial lineage confined to the Sonoran Desert. Our results have direct consequences for conservation of threatened populations, or species, as well as for our understanding of the evolution of delayed implantation in this enigmatic group of small carnivores.

Novel de Novo Genome of Cynopterus brachyotis Reveals Evolutionarily Abrupt Shifts in Gene Family Composition across Fruit Bats

Major novel physiological or phenotypic adaptations often require accompanying modifications at the genic level. Conversely, the detection of considerable contractions and/or expansions of gene families can be an indicator of fundamental but unrecognized physiological change. We sequenced a novel fruit bat genome (Cynopterus brachyotis) and adopted a comparative approach to reconstruct the evolution of fruit bats, mapping contractions and expansions of gene families along their evolutionary history. Despite a radical change in life history as compared with other bats (e.g., loss of echolocation, large size, and frugivory), fruit bats have undergone surprisingly limited change in their genic composition, perhaps apart from a potentially novel gene family expansion relating to telomere protection and longevity. In sharp contrast, within fruit bats, the new Cynopterus genome bears the signal of unusual gene loss and gene family contraction, despite its similar morphology and lifestyle to two other major fruit bat lineages. Most missing genes are regulatory, immune-related, and olfactory in nature, illustrating the diversity of genomic strategies employed by bats to contend with responses to viral infection and olfactory requirements. Our results underscore that significant fluctuations in gene family composition are not always associated with obvious examples of novel physiological and phenotypic adaptations but may often relate to less-obvious shifts in immune strategies.