The new framework for understanding placental mammal evolution

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.

A global catalog of whole-genome diversity from 233 primate species

The rich diversity of morphology and behavior displayed across primate species provides an informative context in which to study the impact of genomic diversity on fundamental biological processes. Analysis of that diversity provides insight into long-standing questions in evolutionary and conservation biology and is urgent given severe threats these species are facing. Here, we present high-coverage whole-genome data from 233 primate species representing 86% of genera and all 16 families. This dataset was used, together with fossil calibration, to create a nuclear DNA phylogeny and to reassess evolutionary divergence times among primate clades. We found within-species genetic diversity across families and geographic regions to be associated with climate and sociality, but not with extinction risk. Furthermore, mutation rates differ across species, potentially influenced by effective population sizes. Lastly, we identified extensive recurrence of missense mutations previously thought to be human specific. This study will open a wide range of research avenues for future primate genomic research.

Divergent vertebral formulae shape the evolution of axial complexity in mammals

Complexity, defined as the number of parts and their degree of differentiation, is a poorly explored aspect of macroevolutionary dynamics. The maximum anatomical complexity of organisms has undoubtedly increased through evolutionary time. However, it is unclear whether this increase is a purely diffusive process or whether it is at least partly driven, occurring in parallel in most or many lineages and with increases in the minima as well as the means. Highly differentiated and serially repeated structures, such as vertebrae, are useful systems with which to investigate these patterns. We focus on the serial differentiation of the vertebral column in 1,136 extant mammal species, using two indices that quantify complexity as the numerical richness and proportional distribution of vertebrae across presacral regions and a third expressing the ratio between thoracic and lumbar vertebrae. We address three questions. First, we ask whether the distribution of complexity values in major mammal groups is similar or whether clades have specific signatures associated with their ecology. Second, we ask whether changes in complexity throughout the phylogeny are biased towards increases and whether there is evidence of driven trends. Third, we ask whether evolutionary shifts in complexity depart from a uniform Brownian motion model. Vertebral counts, but not complexity indices, differ significantly between major groups and exhibit greater within-group variation than recognized hitherto. We find strong evidence of a trend towards increasing complexity, where higher values propagate further increases in descendant lineages. Several increases are inferred to have coincided with major ecological or environmental shifts. We find support for multiple-rate models of evolution for all complexity metrics, suggesting that increases in complexity occurred in stepwise shifts, with evidence for widespread episodes of recent rapid divergence. Different subclades evolve more complex vertebral columns in different configurations and probably under different selective pressures and constraints, with widespread convergence on the same formulae. Further work should therefore focus on the ecological relevance of differences in complexity and a more detailed understanding of historical patterns.

Molecular phylogenies map to biogeography better than morphological ones

Phylogenetic relationships are inferred principally from two classes of data: morphological and molecular. Currently, most phylogenies of extant taxa are inferred from molecules and when morphological and molecular trees conflict the latter are often preferred. Although supported by simulations, the superiority of molecular trees has rarely been assessed empirically. Here we test phylogenetic accuracy using two independent data sources: biogeographic distributions and fossil first occurrences. For 48 pairs of morphological and molecular trees we show that, on average, molecular trees provide a better fit to biogeographic data than their morphological counterparts and that biogeographic congruence increases over research time. We find no significant differences in stratigraphic congruence between morphological and molecular trees. These results have implications for understanding the distribution of homoplasy in morphological data sets, the utility of morphology as a test of molecular hypotheses and the implications of analysing fossil groups for which molecular data are unavailable.

Using biogeographical and phylogenetic data, it is shown that molecular trees fit species geographical data better than trees inferred from morphology, and that these differences are not simply due to better tree resolution.

Ossification pattern of the unusual pisiform in two-toed (Choloepus) and three-toed sloths (Bradypus)

Two-toed (Choloepus sp.) and three-toed (Bradypus sp.) sloths possess short, rounded pisiforms that are rare among mammals and differ from other members of Xenarthra like the giant anteater (Myrmecophaga tridactyla) which retain elongated, rod-like pisiforms in common with most mammals. Using photographs, radiographs, and μCT, we assessed ossification patterns in the pisiform and the paralogous tarsal, the calcaneus, for two-toed sloths, three-toed sloths, and giant anteaters to determine the process by which pisiform reduction occurs in sloths and compare it to other previously studied examples of pisiform reduction in humans and orangutans. Both extant sloth genera achieve pisiform reduction through the loss of a secondary ossification center and the likely disruption of the associated growth plate based on an unusually porous subchondral surface. This represents a third unique mechanism of pisiform reduction among mammals, along with primary ossification center loss in humans and retention of two ossification centers with likely reduced growth periods in orangutans. Given the remarkable similarities between two-toed and three-toed sloth pisiform ossification patterns and the presence of pisiform reduction in fossil sloths, extant sloth pisiform morphology does not appear to represent a recent convergent adaptation to suspensory locomotion, but instead is likely to be an ancestral trait of Folivora that emerged early in the radiation of extant and fossil sloths.

New Middle Eocene proboscidean from Togo illuminates the early evolution of the elephantiform-like dental pattern

Africa has played a pivotal role in the evolution of early proboscideans (elephants and their extinct relatives), yet vast temporal and geographical zones remain uncharted on the continent. A long hiatus encompassing most of the Eocene (Ypresian to the Early Priabonian, around 13 Myr timespan) considerably hampers our understanding of the early evolutionary history of the group. It is notably the case with the origin of its most successful members, the Elephantiformes, i.e. all elephant-like proboscideans most closely related to modern elephants. Here, we describe a proboscidean lower molar discovered in Lutetian phosphate deposits from Togo, and name a new genus and species, Dagbatitherium tassyi. We show that Dagbatitherium displays several elephantiform dental characteristics such as a three-layered Schmelzmuster, the presence of a mesoconid, transversely enlarged buccal cusps and the individualization of a third lophid closely appressed to a minute distal cingulid. Dagbatitherium represents a stem Elephantiformes, pushing back the origin of the group by about 10 Myr, i.e. a third of its currently known evolutionary history. More importantly, Dagbatitherium potentially unlocks the puzzle of the origin of the unique elephantiform tooth crown organization by bridging a critical temporal and morphological gap between early bunodont incipiently bilophodont proboscidean taxa and more derived elephantiforms.

The squirrel is in the detail: Anatomy and morphometry of the tail in Sciuromorpha (Rodentia, Mammalia)

In mammals, the caudal vertebrae are certainly among the least studied elements of their skeleton. However, the tail plays an important role in locomotion (e.g., balance, prehensility) and behavior (e.g., signaling). Previous studies largely focused on prehensile tails in Primates and Carnivora, in which certain osteological features were selected and used to define tail regions (proximal, transitional, distal). Interestingly, the distribution pattern of these anatomical characters and the relative proportions of the tail regions were similar in both orders. In order to test if such tail regionalization can be applied to Rodentia, we investigated the caudal vertebrae of 20 Sciuridae and six Gliridae species. Furthermore, we examined relationships between tail anatomy/morphometry and locomotion. The position of selected characters along the tail was recorded and their distribution was compared statistically using Spearman rank correlation. Vertebral body length (VBL) was measured to calculate the proportions of each tail region and to perform procrustes analysis on the shape of relative vertebral body length (rVBL) progressions. Our results show that tail regionalization, as defined for Primates and Carnivora, can be applied to almost all investigated squirrels, regardless of their locomotor category. Moreover, major locomotor categories can be distinguished by rVBL progression and tail region proportions. In particular, the small flying squirrels Glaucomys volans and Hylopetes sagitta show an extremely short transitional region. Likewise, several semifossorial taxa can be distinguished by their short distal region. Moreover, among flying squirrels, Petaurista petaurista shows differences with the small flying squirrels, mirroring previous observations on locomotory adaptations based on their inner ear morphometry. Our results show furthermore that the tail region proportions of P. petaurista, phylogenetically more basal than the small flying squirrels, are similar to those of bauplan-conservative arboreal squirrels.

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.

An annotated checklist of the insectivores (Mammalia, Eulipotyphla) of Laos

An annotated checklist of the insectivorous mammals (order Eulipotyphla) of Laos is presented for the first time. The checklist is based on an extensive literature review and an examination of available museum specimens from world collections. Confirmed fauna of eulipotyphlans includes 11 species from 7 genera and 3 families. Additionally, several species which are known from neighbouring countries and could potentially occur in Laos are included in this checklist. According to these data, the Laotian fauna can include at least 23 species in 11 genera.

Mechanics of Arboreal Locomotion in Swinhoe’s Striped Squirrels: A Potential Model for Early Euarchontoglires

Differences between arboreal and terrestrial supports likely pose less contrasting functional demands on the locomotor system at a small body size. For arboreal mammals of small body size, asymmetrical gaits have been demonstrated to be advantageous to increase dynamic stability. Many of the extant arboreal squirrel-related rodents display a small body size, claws on all digits, and limited prehensility, a combination that was proposed to have characterized the earliest Euarchontoglires. Thus, motion analysis of such a modern analog could shed light onto the early locomotor evolution of eurarchontoglirans. In this study, we investigated how Swinhoe’s striped squirrels (Tamiops swinhoei; Scuiromorpha) adjust their locomotion when faced with different orientations on broad supports and simulated small branches. We simultaneously recorded high-Hz videos (501 trials) and support reaction forces (451 trials) of squirrels running on two types of instrumented trackways installed at either a 45° incline (we recorded locomotion on inclines and declines) or with a horizontal orientation. The striped squirrels almost exclusively used asymmetrical gaits with a preference for full bounds. Locomotion on simulated branches did not differ substantially from locomotion on the flat trackway. We interpreted several of the quantified adjustments on declines and inclines (in comparison to horizontal supports) as mechanisms to increase stability (e.g., by minimizing toppling moments) and as adjustments to the differential loading of fore- and hind limbs on inclined supports. Our data, in addition to published comparative data and similarities to the locomotion of other small arboreal rodents, tree shrews, and primates as well as a likely small body size at the crown-group node of Euarchontoglires, render a preference for asymmetrical gaits in early members of the clade plausible. This contributes to our understanding of the ancestral lifestyle of this mammalian ‘superclade’.

P*R*O*P: a web application to perform phylogenetic analysis considering the effect of gaps

Background

Phylogenetic analysis strongly depends on evolutionary models. Most evolutionary models for estimating genetic differences and phylogenetic relationships do not treat gap sites in the alignment of sequences. Appropriately incorporating evolutionary information of sites containing insertions and deletions into genetic difference measures will be improve the accuracy of phylogenetic estimates.

Results

We introduced a new measure for estimating genetic differences, and presented P*R*O*P, a web application for performing phylogenetic analysis based on genetic difference considering the effect of gaps. As an example of phylogenetic analysis using P*R*O*P, we used complete p53 amino acid sequences of 31 organisms and illustrated that the genetic differences with and without information on sites containing gaps result in trees with different topologies.

Conclusions

P*R*O*P is available at https://www.rs.tus.ac.jp/bioinformatics/prop and the user can perform phylogenetic analysis by uploading sequence data on the website. The most distinctive feature of P*R*O*P is its genetic difference that is estimated without eliminating gap sites for alignment sequences, which helps users detect meaningful difference in an evolutionary process. The source code is available in GitHub: https://github.com/TUS-Satolab/PROP.

Distribution, number, and certain neurochemical identities of infracortical white matter neurons in the brains of three megachiropteran bat species

A large population of infracortical white matter neurons, or white matter interstitial cells (WMICs), are found within the subcortical white matter of the mammalian telencephalon. We examined WMICs in three species of megachiropterans, Megaloglossus woermanni, Casinycteris argynnis, and Rousettus aegyptiacus, using immunohistochemical and stereological techniques. Immunostaining for neuronal nuclear marker (NeuN) revealed substantial numbers of WMICs in each species-M. woermanni 124,496 WMICs, C. argynnis 138,458 WMICs, and the larger brained R. aegyptiacus having an estimated WMIC population of 360,503. To examine the range of inhibitory neurochemical types we used antibodies against parvalbumin, calbindin, calretinin, and neural nitric oxide synthase (nNOS). The calbindin and nNOS immunostained neurons were the most commonly observed, while those immunoreactive for calretinin and parvalbumin were sparse. The proportion of WMICs exhibiting inhibitory neurochemical profiles was ~26%, similar to that observed in previously studied primates. While for the most part the WMIC population in the megachiropterans studied was similar to that observed in other mammals, the one feature that differed was the high proportion of WMICs immunoreactive to calbindin, whereas in primates (macaque monkey, lar gibbon and human) the highest proportion of inhibitory WMICs contain calretinin. Interestingly, there appears to be an allometric scaling of WMIC numbers with brain mass. Further quantitative comparative work across more mammalian species will reveal the developmental and evolutionary trends associated with this infrequently studied neuronal population.

Dental enamel structure in long-nosed armadillos (Xenarthra: Dasypus) and its evolutionary implications

Most xenarthrans have a reduced and simplified dentition that lacks enamel. However, the presence of prismatic enamel has been recorded in the Eocene armadillos Utaetus buccatus (Euphractinae) and Astegotherium dichotomus (Astegotheriini). Among extant xenarthrans, the occurrence of enamel has been recognized only in the long-nosed armadillo, Dasypus novemcinctus (Dasypodinae), but its microstructure has never been described. In this contribution, we analyse the enamel microstructure in deciduous and permanent teeth of four Dasypus species. In deciduous molariform teeth of some species, we identify an apical cap of vestigial enamel (without crystalline structure), interpreted as an amorphous ameloblastic secretion. In permanent teeth, a thin layer of true enamel is found in the apical portion of unworn molariforms. The enamel is prismatic in D. novemcinctus, but in Dasypus hybridus, Dasypus sabanicola and Dasypus punctatus it is prismless. Taking into account the Eocene species of armadillos, the ancestral condition of enamel in cingulates could have been more complex (as in other placentals) and undergone progressive reduction, as shown in the Dasypus lineage. In light of previous genetic and developmental studies, we review and briefly discuss the processes that can account for the reduction/loss of enamel in extant and extinct armadillos. The retention of enamel and the fact that this genus is the only living xenarthran with two functional generations of teeth support the early divergence of the Dasypus lineage among living cingulates. This is in agreement with morphological and molecular analyses.

Cranial endocast of the stem lagomorph Megalagus and brain structure of basal Euarchontoglires

Early lagomorphs are central to our understanding of how the brain evolved in Glires (rodents, lagomorphs and their kin) from basal members of Euarchontoglires (Glires + Euarchonta, the latter grouping primates, treeshrews, and colugos). Here, we report the first virtual endocast of the fossil lagomorph Megalagus turgidus, from the Orella Member of the Brule Formation, early Oligocene, Nebraska, USA. The specimen represents one of the oldest nearly complete lagomorph skulls known. Primitive aspects of the endocranial morphology in Megalagus include large olfactory bulbs, exposure of the midbrain, a small neocortex and a relatively low encephalization quotient. Overall, this suggests a brain morphology closer to that of other basal members of Euarchontoglires (e.g. plesiadapiforms and ischyromyid rodents) than to that of living lagomorphs. However, the well-developed petrosal lobules in Megalagus, comparable to the condition in modern lagomorphs, suggest early specialization in that order for the stabilization of eye movements necessary for accurate visual tracking. Our study sheds new light on the reconstructed morphology of the ancestral brain in Euarchontoglires and fills a critical gap in the understanding of palaeoneuroanatomy of this major group of placental mammals.

Zmat2 in mammals: conservation and diversification among genes and Pseudogenes

BACKGROUND

Recent advances in genetics and genomics present unique opportunities for enhancing our understanding of mammalian biology and evolution through detailed multi-species comparative analysis of gene organization and expression. Yet, of the more than 20,000 protein coding genes found in mammalian genomes, fewer than 10% have been examined in any detail. Here we elucidate the power of data available in publicly-accessible genomic and genetic resources by querying them to evaluate Zmat2, a minimally studied gene whose human ortholog has been implicated in spliceosome function and in keratinocyte differentiation.

RESULTS

We find extensive conservation in coding regions and overall structure of Zmat2 in 18 mammals representing 13 orders and spanning ~ 165 million years of evolutionary development, and in their encoded proteins. We identify a tandem duplication in the Zmat2 gene and locus in opossum, but not in other monotremes, marsupials, or other mammals, indicating that this event occurred subsequent to the divergence of these species from one another. We also define a collection of Zmat2 pseudogenes in half of the mammals studied, and suggest based on phylogenetic analysis that they each arose independently in the recent evolutionary past.

CONCLUSIONS

Mammalian Zmat2 genes and ZMAT2 proteins illustrate conservation of structure and sequence, along with the development and diversification of pseudogenes in a large fraction of species. Collectively, these observations also illustrate how the focused identification and interpretation of data found in public genomic and gene expression resources can be leveraged to reveal new insights of potentially high biological significance.

New insights from RADseq data on differentiation in the Hottentot golden mole species complex from South Africa

Golden moles (Family Chrysochloridae) are small subterranean mammals, endemic to sub-Saharan Africa, and many of the 21 species are listed as threatened on the IUCN Red List. Most species have highly restricted ranges; however two species, the Hottentot golden mole (Amblysomus hottentotus) and the Cape golden mole (Chrysochloris asiatica) have relatively wide ranges. We recently uncovered cryptic diversity within A. hottentotus, through a phylogeographic analysis of this taxon using two mitochondrial gene regions and a nuclear intron. To further investigate this cryptic diversity, we generated nuclear SNP data from across the genome of A. hottentotus, by means of double-digest restriction-site associated DNA sequencing (ddRADSeq), and mapped reads to the Cape golden mole genome. We conducted a phylogenetic analysis and investigated population differentiation. Our results support the distinctiveness of A. h. meesteri. Furthermore, we provide evidence from nuclear SNPs in support of our previous finding that Central coastal samples represent a unique cryptic lineage that is highly divergent from A. h. pondoliae farther south. Although mtDNA suggests that Umtata may represent a unique lineage sister to A. h. longiceps, mito-nuclear discordance from our RADseq data indicate that these samples may instead be closer to A. h. pondoliae, and therefore may not represent a distinct lineage. We stress the importance of recognizing that understudied populations, such as that of Umtata, may represent populations or ESUs under threat and in need of conservation attention. We present a high-quality filtered SNP dataset, comprising thousands of SNPs, which may serve as a useful resource for future golden mole studies. We have thus added to the growing body of research demonstrating the power and utility of RADseq to investigate population differentiation.

The insulin-like growth factor 2 gene in mammals: Organizational complexity within a conserved locus

The secreted protein, insulin-like growth factor 2 (IGF2), plays a central role in fetal and prenatal growth and development, and is regulated at the genetic level by parental imprinting, being expressed predominantly from the paternally derived chromosome in mice and humans. Here, IGF2/Igf2 and its locus has been examined in 19 mammals from 13 orders spanning ~166 million years of evolutionary development. By using human or mouse DNA segments as queries in genome analyses, and by assessing gene expression using RNA-sequencing libraries, more complexity was identified within IGF2/Igf2 than was annotated previously. Multiple potential 5’ non-coding exons were mapped in most mammals and are presumably linked to distinct IGF2/Igf2 promoters, as shown for several species by interrogating RNA-sequencing libraries. DNA similarity was highest in IGF2/Igf2 coding exons; yet, even though the mature IGF2 protein was conserved, versions of 67 or 70 residues are produced secondary to species-specific maintenance of alternative RNA splicing at a variable intron-exon junction. Adjacent H19 was more divergent than IGF2/Igf2, as expected in a gene for a noncoding RNA, and was identified in only 10/19 species. These results show that common features, including those defining IGF2/Igf2 coding and several non-coding exons, were likely present at the onset of the mammalian radiation, but that others, such as a putative imprinting control region 5’ to H19 and potential enhancer elements 3’ to H19, diversified with speciation. This study also demonstrates that careful analysis of genomic and gene expression repositories can provide new insights into gene structure and regulation.

Improvements in the fossil record may largely resolve current conflicts between morphological and molecular estimates of mammal phylogeny

Phylogenies of mammals based on morphological data continue to show several major areas of conflict with the current consensus view of their relationships, which is based largely on molecular data. This raises doubts as to whether current morphological character sets are able to accurately resolve mammal relationships. We tested this under a hypothetical 'best case scenario' by using ancestral state reconstruction (under both maximum parsimony and maximum likelihood) to infer the morphologies of fossil ancestors for all clades present in a recent comprehensive DNA sequence-based phylogeny of mammals, and then seeing what effect the subsequent inclusion of these predicted ancestors had on unconstrained phylogenetic analyses of morphological data. We found that this resulted in topologies that are highly congruent with the current consensus phylogeny, at least when the predicted ancestors are assumed to be well preserved and densely sampled. Most strikingly, several analyses recovered the monophyly of clades that have never been found in previous morphology-only studies, such as Afrotheria and Laurasiatheria. Our results suggest that, at least in principle, improvements in the fossil record-specifically the discovery of fossil taxa that preserve the ancestral or near-ancestral morphologies of the nodes in the current consensus-may be sufficient to largely reconcile morphological and molecular estimates of mammal phylogeny, even using current morphological character sets.

Bone cortical compactness in 'tree sloths' reflects convergent evolution

Bone remodeling, one of the main processes that regulate bone microstructure, consists of bone resorption followed by the deposition of secondary bone at the same location. Remodeling intensity varies among taxa, but a characteristically compact cortex is ubiquitous in the long bones of mature terrestrial mammals. A previous analysis found that cortical bone in a few 'tree sloth' (Bradypus and Choloepus) specimens is heavily remodeled and characterized by numerous immature secondary osteons, suggesting that these animals were remodeling their bones at high rate until late in their ontogeny. This study aims at testing if this remodeling is generally present in 'tree sloths', using a quantitative analysis of the humeral cortical compactness (CC) among xenarthrans. The results of the investigation of humeral diaphyseal cross-sections of 26 specimens belonging to 10 xenarthran species including specimens from both extinct and extant species indicate that in 'tree sloths' the CC is significantly lower than in the other sampled xenarthrans. No significant difference was found between the CC of the two genera of 'tree sloths'. Our results are consistent with the hypothesis that the cortical bone of 'tree sloths' in general undergoes intense and balanced remodeling that is maintained until late (possibly throughout) in their ontogeny. In the light of xenarthran phylogeny, low CC represents another convergence between the long-separated 'tree sloth' lineages. Although the exact structural and/or functional demands that are associated with this trait are hitherto unknown, several hypotheses are suggested here, including a relationship to their relatively low metabolism and to the mechanical demands imposed upon the bones by the suspensory posture and locomotion, which was independently acquired by the two genera of 'tree sloths'.

Possibilities and limitations of three-dimensional reconstruction and simulation techniques to identify patterns, rhythms and functions of apoptosis in the early developing neural tube

The now classical idea that programmed cell death (apoptosis) contributes to a plethora of developmental processes still has lost nothing of its impact. It is, therefore, important to establish effective three-dimensional (3D) reconstruction as well as simulation techniques to decipher the exact patterns and functions of such apoptotic events. The present study focuses on the question whether and how apoptosis promotes neurulation-associated processes in the spinal cord of Tupaia belangeri (Tupaiidae, Scandentia, Mammalia). Our 3D reconstructions demonstrate that at least two craniocaudal waves of apoptosis consecutively pass through the dorsal spinal cord. The first wave appears to be involved in neural fold fusion and/or in selection processes among premigratory neural crest cells. The second one seems to assist in establishing the dorsal signaling center known as the roof plate. In the hindbrain, in contrast, apoptosis among premigratory neural crest cells progresses craniocaudally but discontinuously, in a segment-specific manner. Unlike apoptosis in the spinal cord, these segment-specific apoptotic events, however, precede later ones that seemingly support neural fold fusion and/or postfusion remodeling. Arguing with Whitehead that biological patterns and rhythms differ in that biological rhythms depend "upon the differences involved in each exhibition of the pattern" (Whitehead in An enquiry concerning the principles of natural knowledge. Cambridge University Press, London, 1919, p. 198) we show that 3D reconstruction and simulation techniques can contribute to distinguish between (static) patterns and (dynamic) rhythms of apoptosis. By deciphering novel patterns and rhythms of developmental apoptosis, our reconstructions help to reconcile seemingly inconsistent earlier findings in chick and mouse embryos, and to create rules for computer simulations.

Loss of RXFP2 and INSL3 genes in Afrotheria shows that testicular descent is the ancestral condition in placental mammals

Descent of testes from a position near the kidneys into the lower abdomen or into the scrotum is an important developmental process that occurs in all placental mammals, with the exception of five afrotherian lineages. Since soft-tissue structures like testes are not preserved in the fossil record and since key parts of the placental mammal phylogeny remain controversial, it has been debated whether testicular descent is the ancestral or derived condition in placental mammals. To resolve this debate, we used genomic data of 71 mammalian species and analyzed the evolution of two key genes (relaxin/insulin-like family peptide receptor 2 [RXFP2] and insulin-like 3 [INSL3]) that induce the development of the gubernaculum, the ligament that is crucial for testicular descent. We show that both RXFP2 and INSL3 are lost or nonfunctional exclusively in four afrotherians (tenrec, cape elephant shrew, cape golden mole, and manatee) that completely lack testicular descent. The presence of remnants of once functional orthologs of both genes in these afrotherian species shows that these gene losses happened after the split from the placental mammal ancestor. These “molecular vestiges” provide strong evidence that testicular descent is the ancestral condition, irrespective of persisting phylogenetic discrepancies. Furthermore, the absence of shared gene-inactivating mutations and our estimates that the loss of RXFP2 happened at different time points strongly suggest that testicular descent was lost independently in Afrotheria. Our results provide a molecular mechanism that explains the loss of testicular descent in afrotherians and, more generally, highlight how molecular vestiges can provide insights into the evolution of soft-tissue characters.

While fossils of whales with legs demonstrate that these species evolved from legged ancestors, the ancestral state of nonfossilizing soft-tissue structures can only be indirectly inferred. This difficulty is also confounded by uncertainties in the phylogenetic relationships between the animals concerned. A prime example is the case of testicular descent, a developmental process that determines the final position of testes, which occurs in most placental mammals but is absent from several afrotherian lineages. Here, we discovered that afrotherians possess remnants of genes known to be required for testicular descent. These “molecular vestiges” show that testicular descent was already present in the placental ancestor and was subsequently lost in Afrotheria. Our study highlights the potential of molecular vestiges in resolving contradictory ancestral states of soft-tissue characters.

Genomic evidence reveals a radiation of placental mammals uninterrupted by the KPg boundary

The timing of the diversification of placental mammals relative to the Cretaceous-Paleogene (KPg) boundary mass extinction remains highly controversial. In particular, there have been seemingly irreconcilable differences in the dating of the early placental radiation not only between fossil-based and molecular datasets but also among molecular datasets. To help resolve this discrepancy, we performed genome-scale analyses using 4,388 loci from 90 taxa, including representatives of all extant placental orders and transcriptome data from flying lemurs (Dermoptera) and pangolins (Pholidota). Depending on the gene partitioning scheme, molecular clock model, and genic deviation from molecular clock assumptions, extensive sensitivity analyses recovered widely varying diversification scenarios for placental mammals from a given gene set, ranging from a deep Cretaceous origin and diversification to a scenario spanning the KPg boundary, suggesting that the use of suboptimal molecular clock markers and methodologies is a major cause of controversies regarding placental diversification timing. We demonstrate that reconciliation between molecular and paleontological estimates of placental divergence times can be achieved using the appropriate clock model and gene partitioning scheme while accounting for the degree to which individual genes violate molecular clock assumptions. A birth-death-shift analysis suggests that placental mammals underwent a continuous radiation across the KPg boundary without apparent interruption by the mass extinction, paralleling a genus-level radiation of multituberculates and ecomorphological diversification of both multituberculates and therians. These findings suggest that the KPg catastrophe evidently played a limited role in placental diversification, which, instead, was likely a delayed response to the slightly earlier radiation of angiosperms.

Macrosystematics of eutherian mammals combining HTS data to expand taxon coverage

In the last few years high-throughput sequencing technologies have permitted significant advances in mammalian phylogenetic studies from a genomic perspective. However, these studies have been restricted to a sparse number of species with available reference genomes. Thus, several issues inside the eutherian mammals phylogeny remain unresolved. This may be due in part to limited taxon sampling, as taxonomic density is known to affect phylogenetic resolution. In this context, we present a protocol to increase taxon coverage using high-throughput sequencing data (RNA or DNA) generated for other biological studies and available in public databases. Following this procedure we addressed pending or controversial issues concerning the phylogenetic position of Dermoptera, Pholidota and Chiroptera, considering multiple and independent loci. Also for Chiroptera and Arctoidea we evaluated the relationships of the lineages that compose it. Although the maximum number of genes used is moderate (95), in some cases taxon coverage doubles that of previous related studies. Globally, all coalescent-based (STAR, MP-EST and ASTRAL) and concatenated (IQ-TREE and BEAST2) methods used for species tree reconstruction were consistent to each other and most of interrogated nodes received high statistical support.

On the development of the chondrocranium and the histological anatomy of the head in perinatal stages of marsupial mammals

An overview of the literature on the chondrocranium of marsupial mammals reveals a relative conservatism in shape and structures. We document the histological cranial anatomy of individuals representing Monodelphis domestica, Dromiciops gliroides, Perameles sp. and Macropus eugenii. The marsupial chondrocranium is generally characterized by the great breadth of the lamina basalis, absence of pila metoptica and large otic capsules. Its most anterior portion (cupula nasi anterior) is robust, and anterior to it there are well-developed tactile sensory structures, functionally important in the neonate. Investigations of ossification centers at and around the nasal septum are needed to trace the presence of certain bones (e.g., mesethmoid, parasphenoid) across marsupial taxa. In many adult marsupials, the tympanic floor is formed by at least three bones: alisphenoid (alisphenoid tympanic process), ectotympanic and petrosal (rostral and caudal tympanic processes); the squamosal also contributes in some diprotodontians. The presence of an entotympanic in marsupials has not been convincingly demonstrated. The tubal element surrounding the auditory tube in most marsupials is fibrous connective tissue rather than cartilage; the latter is the case in most placentals recorded to date. However, we detected fibrocartilage in a late juvenile of Dromiciops, and a similar tissue has been reported for Tarsipes. Contradictory reports on the presence of the tegmen tympani can be found in the literature. We describe a small tegmen tympani in Macropus. Several heterochronic shifts in the timing of development of the chondocranium and associated structures (e.g., nerves, muscles) and in the ossification sequence have been interpreted as largely being influenced by functional requirements related to the altriciality of the newborn marsupial during early postnatal life. Comparative studies of chondocranial development of mammals can benefit from a solid phylogenetic framework, research on non-classical model organisms, and integration with imaging and sectional data derived from computer-tomography.

Golgi Analysis of Neuron Morphology in the Presumptive Somatosensory Cortex and Visual Cortex of the Florida Manatee (Trichechus manatus latirostris)

The current study investigates neuron morphology in presumptive primary somatosensory (S1) and primary visual (V1) cortices of the Florida manatee (Trichechus manatus latirostris) as revealed by Golgi impregnation. Sirenians, including manatees, have an aquatic lifestyle, a large body size, and a relatively large lissencephalic brain. The present study examines neuron morphology in 3 cortical areas: in S1, dorsolateral cortex area 1 (DL1) and cluster cortex area 2 (CL2) and in V1, dorsolateral cortex area 4 (DL4). Neurons exhibited a variety of morphological types, with pyramidal neurons being the most common. The large variety of neuron types present in the manatee cortex was comparable to that seen in other eutherian mammals, except for rodents and primates, where pyramid-shaped neurons predominate. A comparison between pyramidal neurons in S1 and V1 indicated relatively greater dendritic branching in S1. Across all 3 areas, the dendritic arborization pattern of pyramidal neurons was also similar to that observed previously in the afrotherian rock hyrax, cetartiodactyls, opossums, and echidnas but did not resemble the widely bifurcated dendrites seen in the large-brained African elephant. Despite adaptations for an aquatic environment, manatees did not share specific neuron types such as tritufted and star-like neurons that have been found in cetaceans. Manatees exhibit an evolutionarily primitive pattern of cortical neuron morphology shared with most other mammals and do not appear to have neuronal specializations for an aquatic niche.

Ancient collagen reveals evolutionary history of the endemic South American 'ungulates'

Since the late eighteenth century, fossils of bizarre extinct creatures have been described from the Americas, revealing a previously unimagined chapter in the history of mammals. The most bizarre of these are the ‘native’ South American ungulates thought to represent a group of mammals that evolved in relative isolation on South America, but with an uncertain affinity to any particular placental lineage. Many authors have considered them descended from Laurasian ‘condylarths’, which also includes the probable ancestors of perissodactyls and artiodactyls, whereas others have placed them either closer to the uniquely South American xenarthrans (anteaters, armadillos and sloths) or the basal afrotherians (e.g. elephants and hyraxes). These hypotheses have been debated owing to conflicting morphological characteristics and the hitherto inability to retrieve molecular information. Of the ‘native’ South American mammals, only the toxodonts and litopterns persisted until the Late Pleistocene–Early Holocene. Owing to known difficulties in retrieving ancient DNA (aDNA) from specimens from warm climates, this research presents a molecular phylogeny for both Macrauchenia patachonica (Litopterna) and Toxodon platensis (Notoungulata) recovered using proteomics-based (liquid chromatography–tandem mass spectrometry) sequencing analyses of bone collagen. The results place both taxa in a clade that is monophyletic with the perissodactyls, which today are represented by horses, rhinoceroses and tapirs.

Neuron Types in the Presumptive Primary Somatosensory Cortex of the Florida Manatee (Trichechus manatus latirostris)

Within afrotherians, sirenians are unusual due to their aquatic lifestyle, large body size and relatively large lissencephalic brain. However, little is known about the neuron type distributions of the cerebral cortex in sirenians within the context of other afrotherians and aquatic mammals. The present study investigated two cortical regions, dorsolateral cortex area 1 (DL1) and cluster cortex area 2 (CL2), in the presumptive primary somatosensory cortex (S1) in Florida manatees (Trichechus manatus latirostris) to characterize cyto- and chemoarchitecture. The mean neuron density for both cortical regions was 35,617 neurons/mm(3) and fell within the 95% prediction intervals relative to brain mass based on a reference group of afrotherians and xenarthrans. Densities of inhibitory interneuron subtypes labeled against calcium-binding proteins and neuropeptide Y were relatively low compared to afrotherians and xenarthrans and also formed a small percentage of the overall population of inhibitory interneurons as revealed by GAD67 immunoreactivity. Nonphosphorylated neurofilament protein-immunoreactive (NPNFP-ir) neurons comprised a mean of 60% of neurons in layer V across DL1 and CL2. DL1 contained a higher percentage of NPNFP-ir neurons than CL2, although CL2 had a higher variety of morphological types. The mean percentage of NPNFP-ir neurons in the two regions of the presumptive S1 were low compared to other afrotherians and xenarthrans but were within the 95% prediction intervals relative to brain mass, and their morphologies were comparable to those found in other afrotherians and xenarthrans. Although this specific pattern of neuron types and densities sets the manatee apart from other afrotherians and xenarthrans, the manatee isocortex does not appear to be explicitly adapted for an aquatic habitat. Many of the features that are shared between manatees and cetaceans are also shared with a diverse array of terrestrial mammals and likely represent highly conserved neural features. A comparative study across manatees and dugongs is necessary to determine whether these traits are specific to one or more of the manatee species, or can be generalized to all sirenians.

Ancient dates or accelerated rates? Morphological clocks and the antiquity of placental mammals

Analyses of a comprehensive morphological character matrix of mammals using 'relaxed' clock models (which simultaneously estimate topology, divergence dates and evolutionary rates), either alone or in combination with an 8.5 kb nuclear sequence dataset, retrieve implausibly ancient, Late Jurassic-Early Cretaceous estimates for the initial diversification of Placentalia (crown-group Eutheria). These dates are much older than all recent molecular and palaeontological estimates. They are recovered using two very different clock models, and regardless of whether the tree topology is freely estimated or constrained using scaffolds to match the current consensus placental phylogeny. This raises the possibility that divergence dates have been overestimated in previous analyses that have applied such clock models to morphological and total evidence datasets. Enforcing additional age constraints on selected internal divergences results in only a slight reduction of the age of Placentalia. Constraining Placentalia to less than 93.8 Ma, congruent with recent molecular estimates, does not require major changes in morphological or molecular evolutionary rates. Even constraining Placentalia to less than 66 Ma to match the 'explosive' palaeontological model results in only a 10- to 20-fold increase in maximum evolutionary rate for morphology, and fivefold for molecules. The large discrepancies between clock- and fossil-based estimates for divergence dates might therefore be attributable to relatively small changes in evolutionary rates through time, although other explanations (such as overly simplistic models of morphological evolution) need to be investigated. Conversely, dates inferred using relaxed clock models (especially with discrete morphological data and MrBayes) should be treated cautiously, as relatively minor deviations in rate patterns can generate large effects on estimated divergence dates.

Tex19 and Sectm1 concordant molecular phylogenies support co-evolution of both eutherian-specific genes

BACKGROUND

Transposable elements (TE) have attracted much attention since they shape the genome and contribute to species evolution. Organisms have evolved mechanisms to control TE activity. Testis expressed 19 (Tex19) represses TE expression in mouse testis and placenta. In the human and mouse genomes, Tex19 and Secreted and transmembrane 1 (Sectm1) are neighbors but are not homologs. Sectm1 is involved in immunity and its molecular phylogeny is unknown.

METHODS

Using multiple alignments of complete protein sequences (MACS), we inferred Tex19 and Sectm1 molecular phylogenies. Protein conserved regions were identified and folds were predicted. Finally, expression patterns were studied across tissues and species using RNA-seq public data and RT-PCR.

RESULTS

We present 2 high quality alignments of 58 Tex19 and 58 Sectm1 protein sequences from 48 organisms. First, both genes are eutherian-specific, i.e., exclusively present in mammals except monotremes (platypus) and marsupials. Second, Tex19 and Sectm1 have both duplicated in Sciurognathi and Bovidae while they have remained as single copy genes in all further placental mammals. Phylogenetic concordance between both genes was significant (p-value < 0.05) and supported co-evolution and functional relationship. At the protein level, Tex19 exhibits 3 conserved regions and 4 invariant cysteines. In particular, a CXXC motif is present in the N-terminal conserved region. Sectm1 exhibits 2 invariant cysteines and an Ig-like domain. Strikingly, Tex19 C-terminal conserved region was lost in Haplorrhini primates while a Sectm1 C-terminal extra domain was acquired. Finally, we have determined that Tex19 and Sectm1 expression levels anti-correlate across the testis of several primates (ρ = -0.72) which supports anti-regulation.

CONCLUSIONS

Tex19 and Sectm1 co-evolution and anti-regulated expressions support a strong functional relationship between both genes. Since Tex19 operates a control on TE and Sectm1 plays a role in immunity, Tex19 might suppress an immune response directed against cells that show TE activity in eutherian reproductive tissues.

Ontogenetic criteria to distinguish vertebral types on the debated xenarthran synsacrum

The presence of a synsacrum formed by the fusion of vertebrae that come into closed contact with the ilium and ischium is a feature that characterizes the clade Xenarthra. Nevertheless, the proper identity of each vertebral element that forms it is a matter of discussion. In this article, we provide ontogenetic information about skeletal ossification of the xenarthran synsacrum and define the position of the sacrocaudal limit within it. We analyzed the synsacrum of 25 specimens of nonadult and 101 adult armadillos and anteaters: Dasypus hybridus, D. novemcinctus, Chaetophractus vellerosus, C. villosus, Tamandua tetradactyla, and Myrmecophaga tridactyla. Two sets of vertebrae were identified: an anterior set, often attached to the iliac bones, in which transverse processes are originated mainly from an expansion of the base of the neural arches, and secondarily from a lateroventral ossification center. A posterior set is characterized by a series of vertebrae along which extra lateral ossifications (described here for the first time) are developed and form exclusively the transverse processes. Among armadillos, the sacrocaudal limit is set between the last vertebrae attached to the iliac bones and the first vertebrae that form the dorsal border of the sacroischial fenestra. In addition, anterior free caudals also showed extra lateral ossifications forming exclusively the transverse processes, supporting the notion that more posterior synsacrals are in fact caudal vertebrae that were incorporated to the synsacrum. In pilosans, the sacrocaudal limit is set between the first vertebrae that come into contact with the ischial bones and the immediately anterior one. However, the pattern of homologies is obscured by the low resolution in the ontogenetic sequence when compared to that of armadillos.

Placental Evolution within the Supraordinal Clades of Eutheria with the Perspective of Alternative Animal Models for Human Placentation

Here a survey of placental evolution is conducted. Placentation is a key factor for the evolution of placental mammals that had evolved an astonishing diversity. As a temporary organ that does not allow easy access, it is still not well understood. The lack of data also is a restriction for better understanding of placental development, structure, and function in the human. Animal models are essential, because experimental access to the human placenta is naturally restricted. However, there is not a single ideal model that is entirely similar to humans. It is particularly important to establish other models than the mouse, which is characterised by a short gestation period and poorly developed neonates that may provide insights only for early human pregnancy. In conclusion, current evolutionary studies have contributed essentially to providing a pool of experimental models for recent and future approaches that may also meet the requirements of a long gestation period and advanced developmental status of the newborn in the human. Suitability and limitations of taxa as alternative animal models are discussed. However, further investigations especially in wildlife taxa should be conducted in order to learn more about the full evolutionary plasticity of the placenta system.

Vγ9 and Vδ2 T cell antigen receptor genes and butyrophilin 3 (BTN3) emerged with placental mammals and are concomitantly preserved in selected species like alpaca (Vicugna pacos)

Human Vγ9Vδ2 T cells recognize phosphorylated products of isoprenoid metabolism (phosphoantigens) PAg with TCR comprising Vγ9JP γ-chains and Vδ2 δ-chains dependent on butyrophilin 3 (BTN3) expressed by antigen-presenting cells. They are massively activated in many infections and show anti-tumor activity and so far, they have been considered to exist only in higher primates. We performed a comprehensive analysis of databases and identified the three genes in species of both placental magnorders, but not in rodents. The common occurrence or loss of in silico translatable Vγ9, Vδ2, and BTN3 genes suggested their co-evolution based on a functional relationship. In the peripheral lymphocytes of alpaca (Vicugna pacos), characteristic Vγ9JP rearrangements and in-frame Vδ2 rearrangements were found and could be co-expressed in a TCR-negative mouse T cell hybridoma where they rescued CD3 expression and function. Finally, database sequence analysis of the extracellular domain of alpaca BTN3 revealed complete conservation of proposed PAg binding residues of human BTN3A1. In summary, we show emergence and preservation of Vγ9 and Vδ2 TCR genes with the gene of the putative antigen-presenting molecule BTN3 in placental mammals and lay the ground for analysis of alpaca as candidate for a first non-primate species to possess Vγ9Vδ2 T cells.

Layer 4 in primary visual cortex of the awake rabbit: contrasting properties of simple cells and putative feedforward inhibitory interneurons

Extracellular recordings were obtained from two cell classes in layer 4 of the awake rabbit primary visual cortex (V1): putative inhibitory interneurons [suspected inhibitory interneurons (SINs)] and putative excitatory cells with simple receptive fields. SINs were identified solely by their characteristic response to electrical stimulation of the lateral geniculate nucleus (LGN, 3+ spikes at >600 Hz), and simple cells were identified solely by receptive field structure, requiring spatially separate ON and/or OFF subfields. Notably, no cells met both criteria, and we studied 62 simple cells and 33 SINs. Fourteen cells met neither criterion. These layer 4 populations were markedly distinct. Thus, SINs were far less linear (F1/F0 < 1), more broadly tuned to stimulus orientation, direction, spatial and temporal frequency, more sensitive to contrast, had much higher spontaneous and stimulus-driven activity, and always had spatially overlapping ON/OFF receptive subfields. SINs responded to drifting gratings with increased firing rates (F0) for all orientations and directions. However, some SINs showed a weaker modulated (F1) response sharply tuned to orientation and/or direction. SINs responded at shorter latencies than simple cells to stationary stimuli, and the responses of both populations could be sustained or transient. Transient simple cells were more sensitive to contrast than sustained simple cells and their visual responses were more frequently suppressed by high contrasts. Finally, cross-correlation between LGN and SIN spike trains confirmed a fast and precisely timed monosynaptic connectivity, supporting the notion that SINs are well suited to provide a fast feedforward inhibition onto targeted cortical populations.

Evolutionary patterns of bone histology and bone compactness in xenarthran mammal long bones

Bone microstructure reflects physiological characteristics and has been shown to contain phylogenetic and ecological signals. Although mammalian long bone histology is receiving increasing attention, systematic examination of the main clades has not yet been performed. Here we describe the long bone microstructure of Xenarthra based on thin sections representing twenty-two species. Additionally, patterns in bone compactness of humeri and femora are investigated. The primary bone tissue of xenarthran long bones is composed of a mixture of woven, parallel-fibered and lamellar bone. The vascular canals have a longitudinal, reticular or radial orientation and are mostly arranged in an irregular manner. Concentric rows of vascular canals and laminar organization of the tissue are only found in anteater bones. The long bones of adult specimens are marked by dense Haversian bone, a feature that has been noted for most groups of mammals. In the long bones of armadillos, secondary osteons have an oblique orientation within the three-dimensional bone tissue, thus resulting in their irregular shape when the bones are sectioned transversely. Secondary remodeling is generally more extensive in large taxa than in small taxa, and this could be caused by increased loading. Lines of arrested growth are assumed to be present in all specimens, but they are restricted to the outermost layer in bones of armadillos and are often masked by secondary remodeling in large taxa. Parameters of bone compactness show a pattern in the femur that separates Cingulata and Pilosa (Folivora and Vermilingua), with cingulates having a lower compactness than pilosans. In addition, cingulates show an allometric relationship between humeral and femoral bone compactness.