Molecular phylogeny and evolution of the Perissodactyla

Genetic characterization of diagnostic epitopes of cardiac troponin I in African rhinoceros

African rhinoceros undergo chemical immobilization and prolonged transport during translocations for conservation purposes and, hence, experience several pathophysiologic changes, including skeletal muscle injury. Potential concurrent myocardial injury has not been investigated due to a lack of validated immunoassays. We aimed to use inferred cardiac troponin I (cTnI) amino acid sequences of southern white (Ceratotherium simum simum) and southern-central black (Diceros bicornis minor) rhinoceros to assess the potential usefulness of several commercial cTnI immunoassays for detecting cTnI in African rhinoceros. We extracted RNA from the myocardium of deceased rhinoceros (2 white, 1 black rhinoceros) followed by primer design, cDNA synthesis via RT-PCR, and Sanger sequencing. The inferred cTnI amino acid sequences were obtained from the mRNA transcript sequences. The homology of epitope binding sites recognized by capture and detection antibodies in 6 human immunoassays was visually evaluated using aligned inferred rhinoceros cTnI amino acid sequences. Percentage identity between white and black rhinoceros cDNA nucleotide sequences was 99%; inferred amino acid sequences were identical. There were 5 amino acid differences between humans and rhinoceros in the epitope binding sites of immunoassay antibodies; 5 assays contained antibodies against epitopes that were not conserved. For one assay, the single capture antibody targeted a short heterologous epitope (residue 87-91), and cross-reactivity with rhinoceros cTnI was deemed unlikely. For the other 5 assays, complete antibody-epitope homology, or the inclusion of multiple detection or capture antibodies, or targeting of long epitopes, indicated that these assays could be suitable for further investigation of cTnI measurement in African rhinoceros.

Equus mitochondrial pangenome reveals independent domestication imprints in donkeys and horses

Mitochondria are semi-autonomous organelles that play a crucial role in the energy budget of animal cells and are closely related to the locomotor abilities of animals. Equidae is renowned for including two domesticated species with distinct purposes: the endurance-oriented donkey and the power-driven horse, making it an ideal system for studying the relationship between mitochondria and locomotor abilities. In this study, to cover the genetic diversity of donkeys, we sequenced and assembled six new mitochondrial genomes from China. Meanwhile, we downloaded the published mitochondrial genomes of all species within Equus and conducted a comprehensive pan-mitochondrial genome analysis. We found that the mitochondrial genomes of Equus are highly conserved, each encoding 37 genes, including 13 protein-coding genes (PCGs). Phylogenetic analysis based on mitochondrial genomes supports previous research, indicating that the extant species in Equus are divided into three main branches: horses, donkeys, and zebras. Specifically, 761 genetic variants were identified between donkeys and horses, 68 of which were non-synonymous mutations in PCGs, potentially linked to their different locomotor abilities. Structural protein modeling indicated that despite genetic differences, the overall protein structures between donkeys and horses remain similar. This study revealed the mitochondrial genome variation patterns of domesticated animals, offering novelty perspectives on domestication imprints. Additionally, it provides reliable candidate molecular markers for the identification of donkeys and horses.

New Insights into the Molecular Evolution of Tapirus pinchaque (Tapiridae, Perissodactyla) and the Rise and Fall of Tapirus kabomani as a Full Species

Large wild mammals are extremely important in their respective ecological communities and are frequently considered to be emblematic. This is the case of the different tapir species, the largest terrestrial mammals from the Neotropics. Despite their large size and being objects of interest for many naturalists, the field still lacks critical genetics and systematics information about tapir species. In the current work, we analyzed four molecular datasets (mitogenomes, and three nuclear genes, RAG 1-2, IRBP, and BRCA1) of two South American tapirs: the Andean tapir (Tapirus pinchaque) and the alleged new species of tapir, Tapirus kabomani. We derived four main findings. (1) Our molecular phylogenetic analyses showed T. pinchaque as the youngest tapir branch in Neotropics and a sister species of Tapirus terrestris. This contradicts the traditional morphological observations of renowned zoologists and paleontologists, who considered T. pinchaque as the oldest Neotropical tapir. (2) Our data does not support that the alleged T. kabomani is a full species. Rather, it is a specific group within T. terrestris. (3) T. pinchaque is the Neotropical tapir species which yielded the lowest levels of genetic diversity (both for mitochondrial and nuclear data). (4) The spatial genetic structure for T. pinchaque shows differences depending on the type of molecular marker used. With mitogenomes, the spatial structure is relatively weak, whereas with two nuclear genes (RAG 1-2 and IRBP), the spatial structure is highly significant. Curiously, for the other nuclear gene (BRCA1), the spatial structure is practically nonexistent. In any case, the northernmost population of T. pinchaque we studied (Los Nevados National Park in Colombia) was in a peripatric situation and was the most genetically differentiated. This is important for the adequate conservation of this population. (5) T. pinchaque showed clear evidence of population expansion during the last part of the Pleistocene, a period during which the dryness and glacial cold extinguished many large mammals in the Americas. However, T. pinchaque survived and spread throughout the Northern Andes.

Deciphering the influence of evolutionary legacy and functional constraints on the patella: an example in modern rhinoceroses amongst perissodactyls

In mammals, the patella is the biggest sesamoid bone of the skeleton and is of crucial importance in posture and locomotion, ensuring the role of a pulley for leg extensors while protecting and stabilizing the knee joint. Despite its central biomechanical role, the relation between the shape of the patella and functional factors, such as body mass or locomotor habit, in the light of evolutionary legacy are poorly known. Here, we propose a morphofunctional investigation of the shape variation of the patella among modern rhinoceroses and more generally among perissodactyls, this order of ungulates displaying a broad range of body plan, body mass and locomotor habits, to understand how the shape of this sesamoid bone varies between species and relatively to these functional factors. Our investigation, relying on three dimensional geometric morphometrics and comparative analyses, reveals that, within Rhinocerotidae and between the three perissodactyl families, the shape of the patella strongly follows the phylogenetic affinities rather than variations in body mass. The patellar shape is more conservative than initially expected both within and between rhinoceroses, equids and tapirs. The development of a medial angle, engendering a strong mediolateral asymmetry of the patella, appears convergent in rhinoceroses and equids, while tapirs retain a symmetric bone close to the plesiomorphic condition of the order. This asymmetric patella is likely associated with the presence of a "knee locking" mechanism in both equids and rhinos. The emergence of this condition may be related to a shared locomotor habit (transverse gallop) in both groups. Our investigation underlines unexcepted evolutionary constraints on the shape of a sesamoid bone usually considered as mostly driven by functional factors.

Twelve toll-like receptor (TLR) genes in the family Equidae - comparative genomics, selection and evolution

Toll-like receptors (TLRs) represent an important part of the innate immune system. While human and murine TLRs have been intensively studied, little is known about TLRs in non-model species. The order Perissodactyla comprises a variety of free-living and domesticated species exposed to different pathogens in different habitats and is therefore suitable for analyzing the diversity and evolution of immunity-related genes. We analyzed TLR genes in the order Perissodactyla with a focus on the family Equidae. Twelve TLRs were identified by bioinformatic analyses of online genomic resources; their sequences were confirmed in equids by genomic DNA re-sequencing of a panel of nine species. The expression of TLR11 and TLR12 was confirmed in the domestic horse by cDNA sequencing. Phylogenetic reconstruction of the TLR gene family in Perissodactyla identified six sub-families. TLR4 clustered together with TLR5; the TLR1-6-10 subfamily showed a high degree of sequence identity. The average estimated evolutionary divergence of all twelve TLRs studied was 0.3% among the Equidae; the most divergent CDS were those of Equus caballus and Equus hemionus kulan (1.34%) in the TLR3, and Equus africanus somaliensis and Equus quagga antiquorum (2.1%) in the TLR1 protein. In each TLR gene, there were haplotypes shared between equid species, most extensively in TLR3 and TLR9 CDS, and TLR6 amino acid sequence. All twelve TLR genes were under strong negative overall selection. Signatures of diversifying selection in specific codon sites were detected in all TLRs except TLR8. Differences in the selection patterns between virus-sensing and non-viral TLRs were observed.

Mechanisms of Rapid Karyotype Evolution in Mammals

Chromosome reshuffling events are often a foundational mechanism by which speciation can occur, giving rise to highly derivative karyotypes even amongst closely related species. Yet, the features that distinguish lineages prone to such rapid chromosome evolution from those that maintain stable karyotypes across evolutionary time are still to be defined. In this review, we summarize lineages prone to rapid karyotypic evolution in the context of Simpson's rates of evolution-tachytelic, horotelic, and bradytelic-and outline the mechanisms proposed to contribute to chromosome rearrangements, their fixation, and their potential impact on speciation events. Furthermore, we discuss relevant genomic features that underpin chromosome variation, including patterns of fusions/fissions, centromere positioning, and epigenetic marks such as DNA methylation. Finally, in the era of telomere-to-telomere genomics, we discuss the value of gapless genome resources to the future of research focused on the plasticity of highly rearranged karyotypes.

The neonatal southern white rhinoceros ovary contains oogonia in germ cell nests

The northern white rhinoceros is functionally extinct with only two females left. Establishing methods to culture ovarian tissues, follicles, and oocytes to generate eggs will support conservation efforts using in vitro embryo production. To the best of our knowledge, this is the first description of the structure and molecular signature of any rhinoceros, more specifically, we describe the neonatal and adult southern white rhinoceros (Ceratotherium simum simum) ovary; the closest relation of the northern white rhinoceros. Interestingly, all ovaries contain follicles despite advanced age. Analysis of the neonate reveals a population of cells molecularly characterised as mitotically active, pluripotent with germ cell properties. These results indicate that unusually, the neonatal ovary still contains oogonia in germ cell nests at birth, providing an opportunity for fertility preservation. Therefore, utilising ovaries from stillborn and adult rhinoceros can provide cells for advanced assisted reproductive technologies and investigating the neonatal ovaries of other endangered species is crucial for conservation.

Anatomical, Histological and Histochemical Observations of the Eyelids and Orbital Glands in the Lowland Tapir (Tapirus terrestris Linnaeus, 1785) (Perissodactyla: Ceratomorpha)

The lowland tapir is one of four species belonging to the Tapiridae family of the Ceratomorpha suborder, similar to Rhinocerotidae. This study describes anatomy with morphometry, histology (hematoxylin and eosin, Masson-Goldner trichrome, Movat pentachrome, mucicarmine, picro-Mallory trichrome) and histochemistry (PAS, AB pH 1.0, AB pH 2.5; AB pH2.5/PAS and HDI) of the upper and lower eyelids, and superficial gland of the third eyelid with the third eyelid, deep gland of the third eyelid, and lacrimal gland. The aim of the work is to show the features of the above-mentioned structures typical only for Tapiridae, as well as to show the presence of similarities and differences between the families forming the order Perissodactyla. The eyelashes on the upper eyelid were long, while those of the lower eyelid were short and much less prominent. In the upper and lower eyelid sebaceous glands, a characteristic simple alveolar gland producing a mucus-like secretion and poorly developed tarsal glands were observed. The marginal zone of the posterior surface of the eyelids was covered by stratified columnar epithelium with 18-21 layers of nucleated cells, while the bulbar zone of these surfaces was covered by cubic multilayer epithelium with 6-11 non-keratinized layers of cells and with sparse goblet cells. In only lower eyelids, numerous lymphoid nodules, diffuse lymphocytes and high endothelial venules were observed. The superficial gland was an acinar complex which secreted mucous and contained plasma cells within the interlobular and interlobular connective tissue. The upper and lower branches of the third eyelid were the shape of a bent "caudal fin" and were composed of hyaline cartilage, and they contained conjunctiva associated lymphoid tissue (CALT). The deep gland was also an acinar complex producing a serous character and having numerous diffuse lymphocytes. The lacrimal gland was an acinar complex producing seromucous secretions and had numerous plasma cells located in the glandular interstitium. The results of our research indicate that the features of the anatomy of the eyelids and orbital region in the lowland tapir are also typical of the family Tapiridae, but also have features common to the families Equidae and Rhinocerotidae. We confirm the presence of poorly developed tarsal glands in both eyelids as well as presence of a palpebral part of the lacrimal gland in the upper eyelid, which is typical only to Tapirus terrestris.

Evolutionarily Young African Rhinoceros Gammaretroviruses

High-throughput sequences were generated from DNA and cDNA from four Southern white rhinoceros (Ceratotherium simum simum) located in the Taronga Western Plain Zoo in Australia. Virome analysis identified reads that were similar to Mus caroli endogenous gammaretrovirus (McERV). Previous analysis of perissodactyl genomes did not recover gammaretroviruses. Our analysis, including the screening of the updated white rhinoceros (Ceratotherium simum) and black rhinoceros (Diceros bicornis) draft genomes identified high-copy orthologous gammaretroviral ERVs. Screening of Asian rhinoceros, extinct rhinoceros, domestic horse, and tapir genomes did not identify related gammaretroviral sequences in these species. The newly identified proviral sequences were designated SimumERV and DicerosERV for the white and black rhinoceros retroviruses, respectively. Two long terminal repeat (LTR) variants (LTR-A and LTR-B) were identified in the black rhinoceros, with different copy numbers associated with each (n = 101 and 373, respectively). Only the LTR-A lineage (n = 467) was found in the white rhinoceros. The African and Asian rhinoceros lineages diverged approximately 16 million years ago. Divergence age estimation of the identified proviruses suggests that the exogenous retroviral ancestor of the African rhinoceros ERVs colonized their genomes within the last 8 million years, a result consistent with the absence of these gammaretroviruses from Asian rhinoceros and other perissodactyls. The black rhinoceros germ line was colonized by two lineages of closely related retroviruses and white rhinoceros by one. Phylogenetic analysis indicates a close evolutionary relationship with ERVs of rodents including sympatric African rats, suggesting a possible African origin of the identified rhinoceros gammaretroviruses. IMPORTANCE Rhinoceros genomes were thought to be devoid of gammaretroviruses, as has been determined for other perissodactyls (horses, tapirs, and rhinoceros). While this may be true of most rhinoceros, the African white and black rhinoceros genomes have been colonized by evolutionarily young gammaretroviruses (SimumERV and DicerosERV for the white and black rhinoceros, respectively). These high-copy endogenous retroviruses (ERVs) may have expanded in multiple waves. The closest relative of SimumERV and DicerosERV is found in rodents, including African endemic species. Restriction of the ERVs to African rhinoceros suggests an African origin for the rhinoceros gammaretroviruses.

Consideration of genetic variation and evolutionary history in future conservation of Indian one-horned rhinoceros (Rhinoceros unicornis)

Background

The extant members of the Asian rhinos have experienced severe population and range declines since Pleistocene through a combination of natural and anthropogenic factors. The one-horned rhino is the only Asian species recovered from such conditions but most of the extant populations are reaching carrying capacity. India currently harbours ~ 83% of the global wild one-horned rhino populations distributed across seven protected areas. Recent assessments recommend reintroduction-based conservation approaches for the species, and implementation of such efforts would greatly benefit from detailed genetic assessments and evolutionary history of these populations. Using mitochondrial data, we investigated the phylogeography, divergence and demographic history of one-horned rhinos across its Indian range.

Results

We report the first complete mitogenome from all the extant Indian wild one-horned rhino populations (n = 16 individuals). Further, we identified all polymorphic sites and assessed rhino phylogeography (2531 bp mtDNA, n = 111 individuals) across India. Results showed 30 haplotypes distributed as three distinct genetic clades (F_st value 0.68–1) corresponding to the states of Assam (n = 28 haplotypes), West Bengal and Uttar Pradesh (both monomorphic). The reintroduced population of Uttar Pradesh showed maternal signatures of Chitwan National Park, Nepal. Mitochondrial phylogenomics suggests one-horned rhino diverged from its recent common ancestors ~ 950 Kya and different populations (Assam, West Bengal and Uttar Pradesh/Nepal) coalesce at ~ 190–50 Kya, corroborating with the paleobiogeography history of the Indian subcontinent. Further, the demography analyses indicated historical decline in female effective population size ~ 300–200 Kya followed by increasing trends during ~ 110–60 Kya.

Conclusion

The phylogeography and phylogenomic outcomes suggest recognition of three ‘Evolutionary Significant Units (ESUs)’ in Indian rhino. With ongoing genetic isolation of the current populations, future management efforts should focus on identifying genetically variable founder animals and consider periodic supplementation events while planning future rhino reintroduction programs in India. Such well-informed, multidisciplinary approach will be the only way to ensure evolutionary, ecological and demographic stability of the species across its range.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02045-2.

Evolution of the Family Equidae, Subfamily Equinae, in North, Central and South America, Eurasia and Africa during the Plio-Pleistocene

Studies of horse evolution arose during the middle of the 19th century, and several hypotheses have been proposed for their taxonomy, paleobiogeography, paleoecology and evolution. The present contribution represents a collaboration of 19 multinational experts with the goal of providing an updated summary of Pliocene and Pleistocene North, Central and South American, Eurasian and African horses. At the present time, we recognize 114 valid species across these continents, plus 4 North African species in need of further investigation. Our biochronology and biogeography sections integrate Equinae taxonomic records with their chronologic and geographic ranges recognizing regional biochronologic frameworks. The paleoecology section provides insights into paleobotany and diet utilizing both the mesowear and light microscopic methods, along with calculation of body masses. We provide a temporal sequence of maps that render paleoclimatic conditions across these continents integrated with Equinae occurrences. These records reveal a succession of extinctions of primitive lineages and the rise and diversification of more modern taxa. Two recent morphological-based cladistic analyses are presented here as competing hypotheses, with reference to molecular-based phylogenies. Our contribution represents a state-of-the art understanding of Plio-Pleistocene Equus evolution, their biochronologic and biogeographic background and paleoecological and paleoclimatic contexts.

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.

Why wait to mark? Possible reasons behind latency from olfactory exploration to overmarking in four African equid species

Whereas most studies on overmarking in mammals analysed the rate of overmarking, that those investigate time between exploration of an olfactory stimulus and the response to it remain less common, with inconsistent results. We examined the latency in time between elimination by the sender and sniffing by the receiver, and from sniffing and overmarking, in four captive African equid species to explore differences among species, and among age and sex classes. We investigated these latency time periods in light of three potential hypotheses explaining overmarking behaviour in equids: social bonds, group cohesion, and intrasexual competition. Analysing 1684 events of sniffing and 719 of overmarking among 130 individuals, we found that (i) the time from elimination to overmarking was shorter among female friends and in parent-offspring dyads, proving support to the social bond hypothesis; (ii) intraspecific differences in time periods do not reflect the social organisation of species, thus not supporting the group cohesion hypothesis; (iii) males were more attracted to elimination of conspecifics than females, and female's eliminations were inspected longer, in line with the sexual competition hypothesis and/or reproductive behaviour. In addition, we found that the younger foals came to sniff eliminations faster than older ones, and in larger groups foals devoted longer time to sniffing the elimination before overmarking. We concluded that examination of the elimination could be driven by motivations other than the decision to overmark. Whereas overmarking serves to express bonds to a familiar individual, the latency of overmarking reflects more reproductive interests.

Mammalian gut metabolomes mirror microbiome composition and host phylogeny

In the past decade, studies on the mammalian gut microbiome have revealed that different animal species have distinct gut microbial compositions. The functional ramifications of this variation in microbial composition remain unclear: do these taxonomic differences indicate microbial adaptations to host-specific functionality, or are these diverse microbial communities essentially functionally redundant, as has been indicated by previous metagenomics studies? Here, we examine the metabolic content of mammalian gut microbiomes as a direct window into ecosystem function, using an untargeted metabolomics platform to analyze 101 fecal samples from a range of 25 exotic mammalian species in collaboration with a zoological center. We find that mammalian metabolomes are chemically diverse and strongly linked to microbiome composition, and that metabolome composition is further correlated to the phylogeny of the mammalian host. Specific metabolites enriched in different animal species included modified and degraded host and dietary compounds such as bile acids and triterpenoids, as well as fermentation products such as lactate and short-chain fatty acids. Our results suggest that differences in microbial taxonomic composition are indeed translated to host-specific metabolism, indicating that taxonomically distant microbiomes are more functionally diverse than redundant.

Sesamoid bones also show functional adaptation in their microanatomy-The example of the patella in Perissodactyla

The patella is the largest sesamoid bone of the skeleton. It is strongly involved in the knee, improving output force and velocity of the knee extensors, and thus plays a major role in locomotion and limb stability. However, the relationships between its structure and functional constraints, that would enable a better understanding of limb bone functional adaptations, are poorly known. This contribution proposes a comparative analysis, both qualitative and quantitative, of the microanatomy of the whole patella in perissodactyls, which show a wide range of morphologies, masses, and locomotor abilities, in order to investigate how the microanatomy of the patella adapts to evolutionary constraints. The inner structure of the patella consists of a spongiosa surrounded by a compact cortex. Contrary to our expectations, there is no increase in compactness with bone size, and thus body size and weight, but only an increase in the tightness of the spongiosa. No particular thickening of the cortex associated with muscle insertions is noticed but a strong thickening is observed anteriorly at about mid-length, where the strong intermediate patellar ligament inserts. The trabeculae are mainly oriented perpendicularly to the posterior articular surface, which highlights that the main stress is anteroposteriorly directed, maintaining the patella against the femoral trochlea. Conversely, anteriorly, trabeculae are rather circumferentially oriented, following the insertion of the patellar ligament and, possibly also, of the quadriceps tendon. A strong variation is observed among perissodactyl families but also intraspecifically, which is in accordance with previous studies suggesting a higher variability in sesamoid bones. Clear trends are nevertheless observed between the three families. Equids have a much thinner cortex than ceratomorphs. Rhinos and equids, both characterized by a development of the medial border, show an increase in trabecular density laterally suggesting stronger stresses laterally. The inner structure in tapirs is more homogeneous despite the absence of medial development of the medial border with no "compensation" of the inner structure, which suggests different stresses on their knees associated with a different morphology of their patellofemoral joint.

Long bone shape variation in the forelimb of Rhinocerotoidea: relation with size, body mass and body proportions

In quadrupeds, limb bones are strongly affected by functional constraints linked to weight support, but few studies have addressed the complementary effects of mass, size and body proportions on limb bone shape. During their history, Rhinocerotoidea have displayed a great diversity of body masses and relative size and proportions of limb bones, from small tapir-like forms to giant species. Here, we explore the evolutionary variation of shapes in forelimb bones and its relationship with body mass in Rhinocerotoidea. Our results indicate a general increase in robustness and greater development of muscular insertions in heavier species, counteracting the higher weight loadings induced by an increased body mass. The shape of the humerus changes allometrically and exhibits a strong phylogenetic signal. Shapes of the radius and ulna display a stronger link with body mass repartition than with the absolute mass itself. Congruent shape variation between the humerus and the proximal part of the ulna suggests that the elbow joint is comprised of two strongly covariant structures. In addition, our work confirms the uniqueness of giant Paraceratheriidae among Rhinocerotoidea, whose shape variation is related to both a high body mass and a cursorial forelimb construction.

Ancient and modern genomes unravel the evolutionary history of the rhinoceros family

Only five species of the once-diverse Rhinocerotidae remain, making the reconstruction of their evolutionary history a challenge to biologists since Darwin. We sequenced genomes from five rhinoceros species (three extinct and two living), which we compared to existing data from the remaining three living species and a range of outgroups. We identify an early divergence between extant African and Eurasian lineages, resolving a key debate regarding the phylogeny of extant rhinoceroses. This early Miocene (∼16 million years ago [mya]) split post-dates the land bridge formation between the Afro-Arabian and Eurasian landmasses. Our analyses also show that while rhinoceros genomes in general exhibit low levels of genome-wide diversity, heterozygosity is lowest and inbreeding is highest in the modern species. These results suggest that while low genetic diversity is a long-term feature of the family, it has been particularly exacerbated recently, likely reflecting recent anthropogenic-driven population declines.

Milk composition of white rhinoceros over lactation and comparison with other Perissodactyla

The proximate composition of milk from fifteen free-ranging white rhinoceroses at different lactation stages is reported with detailed analysis of fatty acid composition and minerals. Lactose is the main component at 7.93 ± 0.53%, followed by 0.93 ± 0.19% protein, 1.76 ± 0.39% fat, 0.40 ± 0.18% ash, 0.05 ± 0.01% glucose, and 0.04 ± 0.02% non protein nitrogen The interindividual variation of all the components is high, showing no trend of change over lactation. The K and P content decreased over lactation. The fatty acid composition of rhinoceros milk is characteristic with a high saturated fatty acids content of 62%-84%, of which the medium chain fatty acids form the major portion. The C10:0, C12:0, C16:0, and C18:1c9 are the major fatty acids in the milk fat. The results are compared with the two other Perissodactylae families, the Equidae and Tapiridae. Differences in gross composition are small, but the milk of the Rhinocerotidae have the lowest gross energy, while the milk fats consist of the highest amounts of saturated fatty acids, while the low levels of C16:1c9 and C18:1c9 indicate the lowest mammary Δ9-desaturase activity.

The origin of Rhinocerotoidea and phylogeny of Ceratomorpha (Mammalia, Perissodactyla)

Rhinoceroses have been considered to have originated from tapiroids in the middle Eocene; however, the transition remains controversial, and the first unequivocal rhinocerotoids appeared about 4 Ma later than the earliest tapiroids of the Early Eocene. Here we describe 5 genera and 6 new species of rhinoceroses recently discovered from the early Eocene to the early middle Eocene deposits of the Erlian Basin of Inner Mongolia, China. These new materials represent the earliest members of rhinocerotoids, forstercooperiids, and/or hyrachyids, and bridge the evolutionary gap between the early Eocene ceratomorphs and middle Eocene rhinocerotoids. The phylogenetic analyses using parsimony and Bayesian inference methods support their affinities with rhinocerotoids, and also illuminate the phylogenetic relationships and biogeography of Ceratomorpha, although some discrepancies are present between the two criteria. The nearly contemporary occurrence of various rhinocerotoids indicates that the divergence of different rhinocerotoid groups occurred no later than the late early Eocene, which is soon after the split between the rhinocerotoids and the tapiroids in the early early Eocene. However, the Bayesian tip-dating estimate suggests that the divergence of different ceratomorph groups occurred in the middle Paleocene.

Effects of body size on estimation of mammalian area requirements

Accurately quantifying species' area requirements is a prerequisite for effective area-based conservation. This typically involves collecting tracking data on species of interest and then conducting home-range analyses. Problematically, autocorrelation in tracking data can result in space needs being severely underestimated. Based on the previous work, we hypothesized the magnitude of underestimation varies with body mass, a relationship that could have serious conservation implications. To evaluate this hypothesis for terrestrial mammals, we estimated home-range areas with global positioning system (GPS) locations from 757 individuals across 61 globally distributed mammalian species with body masses ranging from 0.4 to 4000 kg. We then applied block cross-validation to quantify bias in empirical home-range estimates. Area requirements of mammals <10 kg were underestimated by a mean approximately15%, and species weighing approximately100 kg were underestimated by approximately50% on average. Thus, we found area estimation was subject to autocorrelation-induced bias that was worse for large species. Combined with the fact that extinction risk increases as body mass increases, the allometric scaling of bias we observed suggests the most threatened species are also likely to be those with the least accurate home-range estimates. As a correction, we tested whether data thinning or autocorrelation-informed home-range estimation minimized the scaling effect of autocorrelation on area estimates. Data thinning required an approximately93% data loss to achieve statistical independence with 95% confidence and was, therefore, not a viable solution. In contrast, autocorrelation-informed home-range estimation resulted in consistently accurate estimates irrespective of mass. When relating body mass to home range size, we detected that correcting for autocorrelation resulted in a scaling exponent significantly >1, meaning the scaling of the relationship changed substantially at the upper end of the mass spectrum.

Insertion of Telomeric Repeats in the Human and Horse Genomes: An Evolutionary Perspective

Interstitial telomeric sequences (ITSs) are short stretches of telomeric-like repeats (TTAGGG)n at nonterminal chromosomal sites. We previously demonstrated that, in the genomes of primates and rodents, ITSs were inserted during the repair of DNA double-strand breaks. These conclusions were derived from sequence comparisons of ITS-containing loci and ITS-less orthologous loci in different species. To our knowledge, insertion polymorphism of ITSs, i.e., the presence of an ITS-containing allele and an ITS-less allele in the same species, has not been described. In this work, we carried out a genome-wide analysis of 2504 human genomic sequences retrieved from the 1000 Genomes Project and a PCR-based analysis of 209 human DNA samples. In spite of the large number of individual genomes analyzed we did not find any evidence of insertion polymorphism in the human population. On the contrary, the analysis of ITS loci in the genome of a single horse individual, the reference genome, allowed us to identify five heterozygous ITS loci, suggesting that insertion polymorphism of ITSs is an important source of genetic variability in this species. Finally, following a comparative sequence analysis of horse ITSs and of their orthologous empty loci in other Perissodactyla, we propose models for the mechanism of ITS insertion during the evolution of this order.

Recent mitochondrial lineage extinction in the critically endangered Javan rhinoceros

The Javan rhinoceros (Rhinoceros sondaicus) is one of five extant rhinoceros species and among the rarest large mammals on Earth. Once widespread across Southeast Asia, it is now on the verge of extinction, with only one wild population remaining (estimated at ~60 individuals) on the island of Java, Indonesia. To assess the past genetic diversity of the female lineage of R. sondaicus, we generated mitochondrial genome data from eight museum specimens dating back to the 19th century, before the range of the Javan rhinoceros was dramatically reduced, for comparison against mitochondrial DNA (mtDNA) sequences of current R. sondaicus and other rhinoceros species. We succeeded in reconstructing five full and three partial ancient mitogenomes from the eight samples. We used BEAST to assess the phylogenetic relationship of the five extant rhinoceros species and the historical samples. The results show that the oldest and most diverse mtDNA lineages of R. sondaicus are found in historical samples, indicating a significant reduction of mtDNA diversity in modern Javan rhinos. We anticipate that the newly sequenced data will represent a useful resource for improving our understanding of evolutionary history of this species, should future studies be able to increase the available dataset. We hope this information may help in conservation efforts for this species.

Accelerated Evolution of Limb-Related Gene Hoxd11 in the Common Ancestor of Cetaceans and Ruminants (Cetruminantia)

Reduced numbers of carpal and tarsal bones (wrist and ankle joints) are extensively observed in the clade of Cetacea and Ruminantia (Cetruminantia). Homebox D11 (Hoxd11) is one of the important genes required for limb development in mammals. Mutations in Hoxd11 can lead to defects in particular bones of limbs, including carpus and tarsus. To test whether evolutionary changes in Hoxd11 underlie the loss of these bones in Cetruminantia, we sequenced and analyzed Hoxd11 coding sequences and compared them with other 5' HoxA and HoxD genes in a taxonomic coverage of Cetacea, Ruminantia and other mammalian relatives. Statistical tests on the Hoxd11 sequences found an accelerated evolution in the common ancestor of cetaceans and ruminants, which coincided with the reduction of carpal and tarsal bones in this clade. Five amino acid substitutions (G222S, G227A, G229S, A240T and G261V) and one amino acid deletion (G254Del) occurred in this lineage. In contrast, other 5' HoxA and HoxD genes do not show this same evolutionary pattern, but instead display a highly conserved pattern of evolution in this lineage. Accelerated evolution of Hoxd11, but not other 5' HoxA and HoxD genes, is probably related to the reduction of the carpal and tarsal bones in Cetruminantia. Moreover, we found two amino acid substitutions (G110S and D223N) in Hoxd11 that are unique to the lineage of Cetacea, which coincided with hindlimb loss in the common ancestor of cetaceans. Our results give molecular evidence of Hoxd11 adaptive evolution in cetaceans and ruminants, which could be correlated with limb morphological adaptation.

Influence of mass on tarsus shape variation: a morphometrical investigation among Rhinocerotidae (Mammalia: Perissodactyla)

Many tetrapod lineages show extreme increases in body mass in their evolutionary history, associated with important osteological changes. The ankle joint, essential for foot movement, is assumed to be particularly affected in this regard. We investigated the morphological adaptations of the astragalus and the calcaneus in Rhinocerotidae, and analysed them in light of a comparative analysis with other Perissodactyla. We performed 3D geometric morphometrics and correlated shape with centroid size of the bone and body mass of the species. Our results show that mass has an influence on bone shape in Rhinocerotidae and in Perissodactyla, but this is not as strong as expected. In heavy animals the astragalus has a flatter trochlea, orientated more proximally, associated with a more upright posture of the limb. The calcaneus is more robust, possibly to sustain the greater tension force exerted by the muscles during plantarflexion. Both bones show wider articular facets, providing greater cohesion and better dissipation of the loading forces. The body plan of the animals also has an influence. Short-legged Teleoceratina have a flatter astragalus than the other rhinocerotids. Paraceratherium has a thinner calcaneus than expected. This study clarifies adaptations to high body weight among Rhinocerotidae and calls for similar investigations in other groups with massive forms.

Trehalase Gene as a Molecular Signature of Dietary Diversification in Mammals

Diet is a key factor in determining and structuring animal diversity and adaptive radiations. The mammalian fossil record preserves phenotypic evidence of many dietary shifts, whereas genetic changes followed by dietary diversification in mammals remain largely unknown. To test whether living mammals preserve molecular evidence of dietary shifts, we examined the trehalase gene (Treh), which encodes an enzyme capable of digesting trehalose from insect blood, in bats and other mammals with diverse diets. Bats represent the largest dietary radiation among all mammalian orders, with independent origins of frugivory, nectarivory, carnivory, omnivory, and even sanguivory in an otherwise insectivorous clade. We found that Treh has been inactivated in unrelated bat lineages that independently radiated into noninsectivorous niches. Consistently, purifying selection has been markedly relaxed in noninsectivorous bats compared with their insectivorous relatives. Enzymatic assays of intestinal trehalase in bats suggest that trehalase activity tends to be lost or markedly reduced in noninsectivorous bats compared with their insectivorous relatives. Furthermore, our survey of Treh in 119 mammal species, which represent a deeper evolutionary timeframe, additionally identified a number of other independent losses of Treh in noninsectivorous species, recapitulating the evolutionary pattern that we found in bats. These results document a molecular record of dietary diversification in mammals, and suggest that such molecular signatures of dietary shifts would help us understand both historical and modern changes of animal diets.

Phylogenetics, patterns of genetic variation and population dynamics of Trypanosoma terrestris support both coevolution and ecological host-fitting as processes driving trypanosome evolution

BACKGROUND

A considerable amount of evidence has favored ecological host-fitting, rather than coevolution, as the main mechanism responsible for trypanosome divergence. Nevertheless, beyond the study of human pathogenic trypanosomes, the genetic basis of host specificity among trypanosomes isolated from forest-inhabiting hosts remains largely unknown.

METHODS

To test possible scenarios on ecological host-fitting and coevolution, we combined a host capture recapture strategy with parasite genetic data and studied the genetic variation, population dynamics and phylogenetic relationships of Trypanosoma terrestris, a recently described trypanosome species isolated from lowland tapirs in the Brazilian Pantanal and Atlantic Forest biomes.

RESULTS

We made inferences of T. terrestris population structure at three possible sources of genetic variation: geography, tapir hosts and 'putative' vectors. We found evidence of a bottleneck affecting the contemporary patterns of parasite genetic structure, resulting in little genetic diversity and no evidence of genetic structure among hosts or biomes. Despite this, a strongly divergent haplotype was recorded at a microgeographical scale in the landscape of Nhecolândia in the Pantanal. However, although tapirs are promoting the dispersion of the parasites through the landscape, neither geographical barriers nor tapir hosts were involved in the isolation of this haplotype. Taken together, these findings suggest that either host-switching promoted by putative vectors or declining tapir population densities are influencing the current parasite population dynamics and genetic structure. Similarly, phylogenetic analyses revealed that T. terrestris is strongly linked to the evolutionary history of its perissodactyl hosts, suggesting a coevolving scenario between Perissodactyla and their trypanosomes. Additionally, T. terrestris and T. grayi are closely related, further indicating that host-switching is a common feature promoting trypanosome evolution.

CONCLUSIONS

This study provides two lines of evidence, both micro- and macroevolutionary, suggesting that both host-switching by ecological fitting and coevolution are two important and non-mutually-exclusive processes driving the evolution of trypanosomes. In line with other parasite systems, our results support that even in the face of host specialization and coevolution, host-switching may be common and is an important determinant of parasite diversification.

Comparative forelimb myology and muscular architecture of a juvenile Malayan tapir (Tapirus indicus)

The absence of preserved soft tissues in the fossil record is frequently a hindrance for palaeontologists wishing to investigate morphological shifts in key skeletal systems, such as the limbs. Understanding the soft tissue composition of modern species can aid in understanding changes in musculoskeletal features through evolution, including those pertaining to locomotion. Establishing anatomical differences in soft tissues utilising an extant phylogenetic bracket can, in turn, assist in interpreting morphological changes in hard tissues and modelling musculoskeletal movements during evolutionary transitions (e.g. digit reduction in perissodactyls). Perissodactyls (horses, rhinoceroses, tapirs and their relatives) are known to have originated with a four-toed (tetradactyl) forelimb condition. Equids proceeded to reduce all but their central digit, resulting in monodactyly, whereas tapirs retained the ancestral tetradactyl state. The modern Malayan tapir (Tapirus indicus) has been shown to exhibit fully functional tetradactyly in its forelimb, more so than any other tapir, and represents an ideal case-study for muscular arrangement and architectural comparison with the highly derived monodactyl Equus. Here, we present the first quantification of muscular architecture of a tetradactyl perissodactyl (T. indicus), and compare it to measurements from modern monodactyl caballine horse (Equus ferus caballus). Each muscle of the tapir forelimb was dissected out from a cadaver and measured for architectural properties: muscle-tendon unit (MTU) length, MTU mass, muscle mass, pennation angle, and resting fibre length. Comparative parameters [physiological cross-sectional area (PCSA), muscle volume, and % muscle mass] were then calculated from the raw measurements. In the shoulder region, the infraspinatus of T. indicus exhibits dual origination sites on either side of the deflected scapular spine. Within ungulates, this condition has only been previously reported in suids. Differences in relative contribution to limb muscle mass between T. indicus and Equus highlight forelimb muscles that affect mobility in the lateral and medial digits (e.g. extensor digitorum lateralis). These muscles were likely reduced in equids during their evolutionary transition from tetradactyl forest-dwellers to monodactyl, open-habitat specialists. Patterns of PCSA across the forelimb were similar between T. indicus and Equus, with the notable exceptions of the biceps brachii and flexor carpi ulnaris, which were much larger in Equus. The differences observed in PCSA between the tapir and horse forelimb muscles highlight muscles that are essential for maintaining stability in the monodactyl limb while moving at high speeds. This quantitative dataset of muscle architecture in a functionally tetradactyl perissodactyl is a pivotal first step towards reconstructing the locomotor capabilities of extinct, four-toed ancestors of modern perissodactyls, and providing further insights into the equid locomotor transition.

Early Pleistocene enamel proteome from Dmanisi resolves Stephanorhinus phylogeny

Ancient DNA (aDNA) sequencing has enabled reconstruction of speciation, migration, and admixture events for extinct taxa¹. Outside the permafrost, however, irreversible aDNA post-mortem degradation² has so far limited aDNA recovery to the past ~0.5 million years (Ma)³. Contrarily, tandem mass spectrometry (MS) allowed sequencing ~1.5 million year (Ma) old collagen type I (COL1)⁴ and suggested the presence of protein residues in Cretaceous fossil remains⁵, although with limited phylogenetic use⁶. In the absence of molecular evidence, the speciation of several Early and Middle Pleistocene extinct species remain contentious. In this study, we address the phylogenetic relationships of the Eurasian Pleistocene Rhinocerotidae^7–9 using a ~1.77 Ma old dental enamel proteome of a Stephanorhinus specimen from the Dmanisi archaeological site in Georgia (South Caucasus)¹⁰. Molecular phylogenetic analyses place the Dmanisi Stephanorhinus as a sister group to the woolly (Coelodonta antiquitatis) and Merck’s rhinoceros (S. kirchbergensis) clade. We show that Coelodonta evolved from an early Stephanorhinus lineage and that the latter includes at least two distinct evolutionary lines. As such, the genus Stephanorhinus is currently paraphyletic and its systematic revision is therefore needed. We demonstrate that Early Pleistocene dental enamel proteome sequencing overcomes the limits of ancient collagen- and aDNA-based phylogenetic inference. It also provides additional information about the sex and taxonomic assignment of the specimens analysed. Dental enamel, the hardest tissue in vertebrates¹¹, is highly abundant in the fossil record. Our findings reveal that palaeoproteomic investigation of this material can push biomolecular investigation further back into the Early Pleistocene.

An interferon-gamma release assay for the diagnosis of the Mycobacterium bovis infection in white rhinoceros (Ceratotherium simum)

Mycobacterium bovis (M. bovis), the cause of bovine tuberculosis, is endemic in Kruger National Park (KNP), South Africa. The risk of spread of M. bovis infection currently prevents translocation of white rhinoceros (Ceratotherium simum) from this population. Therefore, accurate assays are necessary for screening this threatened species. Interferon gamma (IFN-γ) release assays (IGRA) are commonly used for tuberculosis diagnosis in humans and other wildlife species. Hence, the aim of this study was to develop an IGRA for M. bovis detection in white rhinoceros. Heparinized whole blood was collected from immobilized white rhinoceros in KNP (n = 131) and incubated overnight in QuantiFERON®-TB Gold (QFT) blood collection tubes, after which the plasma was harvested following centrifugation. Tissue samples for mycobacterial culture were available from a subset of 21 rhinoceros. The concentration of IFN-γ in plasma samples was measured using the Mabtech equine IFN-γ ELISA^PRO kit. An IGRA result was calculated as the difference in IFN-γ concentrations in the QFT Nil and TB antigen tubes. Using test results for the white rhinoceros with known infection status, a diagnostic cut-off value was calculated as 21 pg/ml. Additionally, cut-off values for IFN-γ concentrations for plasma from QFT Nil and QFT Mitogen tubes were calculated to increase confidence in IGRA result interpretation. The combination of the QFT stimulation platform and Mabtech equine IFN-γ ELISA is a promising diagnostic test to distinguish between of M. bovis-infected and -uninfected white rhinoceros.

Evolution, Systematics, and Biogeography of the Triatominae, Vectors of Chagas Disease

In this chapter, we review and update current knowledge about the evolution, systematics, and biogeography of the Triatominae (Hemiptera: Reduviidae)-true bugs that feed primarily on vertebrate blood. In the Americas, triatomines are the vectors of Trypanosoma cruzi, the etiological agent of Chagas disease. Despite declining incidence and prevalence, Chagas disease is still a major public health concern in Latin America. Triatomines occur also in the Old World, where vector-borne T. cruzi transmission has not been recorded. Triatomines evolved from predatory reduviid bugs, most likely in the New World, and diversified extensively across the Americas (including the Caribbean) and in parts of Asia and Oceania. Here, we first discuss our current understanding of how, how many times, and when the blood-feeding habit might have evolved among the Reduviidae. Then we present a summary of recent advances in the systematics of this diverse group of insects, with an emphasis on the contribution of molecular tools to the clarification of taxonomic controversies. Finally, and in the light of both up-to-date phylogenetic hypotheses and a thorough review of distribution records, we propose a global synthesis of the biogeography of the Triatominae. Over 130 triatomine species contribute to maintaining T. cruzi transmission among mammals (sometimes including humans) in almost every terrestrial ecoregion of the Americas. This means that Chagas disease will never be eradicated and underscores the fact that effective disease prevention will perforce require stronger, long-term vector control-surveillance systems.

A Primitive Trait in Two Breeds of Equus Caballus Revealed by Comparative Anatomy of the Distal Limb

Simple Summary

Understanding the complexities and evolutionary links between extinct and extant equids has been vital to genetic conservation and preservation of primitive traits. As domestication of the equid expanded, the loss of primitive traits that ensured survival in a wild environment has not been documented. In this study, the presence of functional interosseous muscle II and IV in the distal limb has been reported, and yet its existence could only be confirmed in relatives and two closely bred descendants of the extinct Tarpan. The morphology described was ligamentous in structure displaying clear longitudinal fibres with a skeletal origin and soft tissue insertion into the medial and lateral branches of the interosseous muscle III (suspensory ligament) dorsal to the sesamoids, similar in orientation to the flexor digitorum profundus ligamentum accessorium (inferior check ligament). Hence, providing a functional medial and lateral stability to the metacarpophalangeal joint (fetlock joint), which equates to one of the functions of the medial and lateral digits in the Mesohippus and Merychippus. The comparable anatomic links between species of the same family that experienced geographical isolation yet display structural conformity appears to be in response to a specific environment. Surmising this potential remnant of functional evolution is a primitive trait and not a breed anomaly.

Abstract

The 55-million-year history of equine phylogeny has been well-documented from the skeletal record; however, this is less true for the soft tissue structures that are now vestigial in modern horse. A recent study reported that two ligamentous structures resembling functional interosseous muscle II and IV were evident in Dutch Konik horses. The current study investigates this finding and compares it to members of the genus Equus to identify either a breed anomaly or functional primitive trait. Distal limbs (n = 574) were dissected from four species of Equus; E. caballus, E. asinus, E. przewalskii and E. quagga boehmi. E. caballus is represented by 18 breeds of horse, including the primitive Dutch Konik’. The interosseous muscle II and IV were evident in all four species, but only two breeds of E. caballus expressed this trait-the Dutch Konik and Bosnian Mountain Horse. These two breeds were the only close descendants of the extinct Equus ferus ferus (Tarpan) represented in this study. In conclusion, the interosseous muscle II and IV originated from the distal nodule of metacarpal II and IV, respectively, and inserted into the corresponding branches of interosseous muscle III proximal to the sesamoids. This suggests a functional role in medial and lateral joint stability and a primitive trait in modern equids.

A commercial ELISA for detection of interferon gamma in white rhinoceros

Bovine tuberculosis (bTB), caused by Mycobacterium bovis, is endemic in Kruger National Park, South Africa, home to the largest population of white rhinoceros (Ceratotherium simum) in the world. In 2016, the first cases of naturally occurring bTB were reported in white rhinoceros; however, there is a lack of understanding of infection and disease process in this species. Prevention and control of transmission depends on the availability of accurate tools to detect M. bovis infection. Interferon gamma (IFN-γ) assays are a reliable detection method for TB in other animal species, and studies have indicated that these tests can be used in white rhinoceros. We sought to screen and optimize a commercial IFN-γ enzyme-linked immunosorbent assay (ELISA) to detect endogenous white rhinoceros IFN-γ in mitogen-stimulated whole blood as a basis for developing a test for M. bovis infection. Optimizations included identifying ELISA antibodies and determining the effect of sample matrix, ELISA plate incubation temperature, ELISA linearity, assay reproducibility, and the assay's limit of quantification. The optimized assay employed an equine IFN-γ antibody pair that was used to create a commercial ELISA kit. This ELISA had a linear response to recombinant equine and endogenous rhinoceros IFN-γ (range: 7.8-125 pg/mL). When incubated at 37°C, the ELISA was highly reproducible, with an optimal recovery and a low limit of quantification, indicating that the Mabtech equine IFN-γ ELISA^PRO kit is a robust assay for measuring white rhinoceros IFN-γ.

Distribution, Diversity, and Evolution of Endogenous Retroviruses in Perissodactyl Genomes

The evolution of mammalian genomes has been shaped by interactions with endogenous retroviruses (ERVs). In this study, we investigated the distribution and diversity of ERVs in the mammalian order Perissodactyla, with a view to understanding their impact on the evolution of modern equids (family Equidae). We characterize the major ERV lineages in the horse genome in terms of their genomic distribution, ancestral genome organization, and time of activity. Our results show that subsequent to their ancestral divergence from rhinoceroses and tapirs, equids acquired four novel ERV lineages. We show that two of these ERV lineages proliferated extensively in the lineage leading to modern horses, and one contains loci that are actively transcribed in specific tissues. In addition, we show that the white rhinoceros has resisted germ line colonization by retroviruses for more than 54 million years-longer than any other extant mammalian species. The map of equine ERVs that we provide here will be of great utility to future studies aiming to investigate the potential functional roles of equine ERVs and their impact on equine evolution.IMPORTANCE ERVs in the host genome are highly informative about the long-term interactions of retroviruses and hosts. They are also interesting because they have influenced the evolution of mammalian genomes in various ways. In this study, we derive a calibrated timeline describing the process through which ERV diversity has been generated in the equine germ line. We determined the distribution and diversity of perissodactyl ERV lineages and inferred their retrotranspositional activity during evolution, thereby gaining insight into the long-term coevolutionary history of retroviruses and mammals. Our study provides a platform for future investigations to identify equine ERV loci involved in physiological processes and/or pathological conditions.

Evolution and extinction of the giant rhinoceros Elasmotherium sibiricum sheds light on late Quaternary megafaunal extinctions

Understanding extinction events requires an unbiased record of the chronology and ecology of victims and survivors. The rhinoceros Elasmotherium sibiricum, known as the 'Siberian unicorn', was believed to have gone extinct around 200,000 years ago-well before the late Quaternary megafaunal extinction event. However, no absolute dating, genetic analysis or quantitative ecological assessment of this species has been undertaken. Here, we show, by accelerator mass spectrometry radiocarbon dating of 23 individuals, including cross-validation by compound-specific analysis, that E. sibiricum survived in Eastern Europe and Central Asia until at least 39,000 years ago, corroborating a wave of megafaunal turnover before the Last Glacial Maximum in Eurasia, in addition to the better-known late-glacial event. Stable isotope data indicate a dry steppe niche for E. sibiricum and, together with morphology, a highly specialized diet that probably contributed to its extinction. We further demonstrate, with DNA sequencing data, a very deep phylogenetic split between the subfamilies Elasmotheriinae and Rhinocerotinae that includes all the living rhinoceroses, settling a debate based on fossil evidence and confirming that the two lineages had diverged by the Eocene. As the last surviving member of the Elasmotheriinae, the demise of the 'Siberian unicorn' marked the extinction of this subfamily.

Horse Y chromosome assembly displays unique evolutionary features and putative stallion fertility genes

Dynamic evolutionary processes and complex structure make the Y chromosome among the most diverse and least understood regions in mammalian genomes. Here, we present an annotated assembly of the male specific region of the horse Y chromosome (eMSY), representing the first comprehensive Y assembly in odd-toed ungulates. The eMSY comprises single-copy, equine specific multi-copy, PAR transposed, and novel ampliconic sequence classes. The eMSY gene density approaches that of autosomes with the highest number of retained X-Y gametologs recorded in eutherians, in addition to novel Y-born and transposed genes. Horse, donkey and mule testis RNAseq reveals several candidate genes for stallion fertility. A novel testis-expressed XY ampliconic sequence class, ETSTY7, is shared with the parasite Parascaris genome, providing evidence for eukaryotic horizontal transfer and inter-chromosomal mobility. Our study highlights the dynamic nature of the Y and provides a reference sequence for improved understanding of equine male development and fertility.

Evolutionary Associations of Endosymbiotic Ciliates Shed Light on the Timing of the Marsupial-Placental Split

Trichostome ciliates are among the most conspicuous protists in the gastrointestinal tract of a large variety of vertebrates. However, little is still known about phylogeny of the trichostome/vertebrate symbiotic systems, evolutionary correlations between trichostome extrinsic traits, and character-dependent diversification of trichostomes. These issues were investigated here, using the relaxed molecular clock technique along with stochastic mapping of character evolution, and binary-state speciation and extinction models. Clock analyses revealed that trichostomes colonized the vertebrate gastrointestinal tract ∼135 Ma, that is, near the paleontological minimum for the split of therian mammals into marsupials and placentals. According to stochastic mapping, the last common ancestor of trichostomes most likely invaded the hindgut of a mammal. Although multiple shifts to fish/amphibian or avian hosts and to the foregut compartments took place during the trichostome phylogeny, only transition to the foregut was recognized as a key innovation responsible for the explosive radiation of ophryoscolecid trichostomes after the Cretaceous/Tertiary boundary, when ungulates began their diversification. Since crown radiations of main trichostome lineages follow those of their mammalian hosts and are in agreement with their historic dispersal routes, the present time-calibrated phylogeny might help to elucidate controversies in the geological and molecular timing of the split between marsupials and placental mammals.

The Disappearing Lamellae: Implications of New Findings in the Family Equidae Suggest the Loss of Nuchal Ligament Lamellae on C6 and C7 Occurred After Domestication

Conflicting data predominantly describes the equine nuchal ligament lamellae (NLL) attachments as C2-C6 or C2-C7; however, preliminary evidence suggests C2-C5. This study aimed to identify morphological variations in the attachments of the NLL in four species of the family equidae and determine if the variations were species or breed linked. Cadaveric examination evaluated 98 equids; Equusasinus (n = 2), Equus ferus caballus (n = 93), Equus przewalskii (n = 2), and Equus quagga boehmi (n = 1). Twenty breeds of E. f. caballus were included according to breed type-modern horse (n = 81), ancient type (n = 10), primitive (n = 4). Of the 98 equids, observations were divided into three categories according to the NLL attachments: C2-C5, C2-C6, and C2-C7; in category C2-C7, E. f. caballus (n = 3/93), modern horse 0 of 81, ancient 1 of 8, primitive 2 of 4; in category C2-C6, modern horse 2 of 81; in category C2-C5, modern horse 79 of 81, ancient type 7 of 8, primitive 2 of 4. Equus asinus, E. przewalskii, and E. q. boehmi all expressed category C2-C7. In total, 88 of 98 expressed category C2-C5; 2 of 98 category C2-C6; and 8 of 98 category C2-C7. These findings suggest that category C2-C5 provides sufficient evidence to suggest a normal occurrence. The high incidence of an absent NLL on C6 and C7 may be considered a contributing factor in caudal cervical osteoarthritis. Furthermore, category C2-C7 still exists in an ancient and primitive breed closely related to the tarpan, implying the NLL may have been attached from C2 to C7 in prehistoric horse before domestication. This hypothesis is supported by the findings of E. asinus, E. przewalskii, and E. q. boehmi, displaying attachments from C2 to C7.

Contrasted patterns of evolution of the LINE-1 retrotransposon in perissodactyls: the history of a LINE-1 extinction

Background

LINE-1 (L1) is the dominant autonomously replicating non-LTR retrotransposon in mammals. Although our knowledge of L1 evolution across the tree of life has considerably improved in recent years, what we know of L1 evolution in mammals is biased and comes mostly from studies in primates (mostly human) and rodents (mostly mouse). It is unclear if patterns of evolution that are shared between those two groups apply to other mammalian orders. Here we performed a detailed study on the evolution of L1 in perissodactyls by making use of the complete genome of the domestic horse and of the white rhinoceros. This mammalian order offers an excellent model to study the extinction of L1 since the rhinoceros is one of the few mammalian species to have lost active L1.

Results

We found that multiple L1 lineages, carrying different 5’UTRs, have been simultaneously active during the evolution of perissodactyls. We also found that L1 has continuously amplified and diversified in horse. In rhinoceros, L1 was very prolific early on. Two successful families were simultaneously active until ~20my ago but became extinct suddenly at exactly the same time.

Conclusions

The general pattern of L1 evolution in perissodactyls is very similar to what was previously described in mouse and human, suggesting some commonalities in the way mammalian genomes interact with L1. We confirmed the extinction of L1 in rhinoceros and we discuss several possible mechanisms.

Electronic supplementary material

The online version of this article (10.1186/s13100-018-0117-4) contains supplementary material, which is available to authorized users.

Early fetal sexing in the rhinoceros by detection of male-specific genes in maternal serum

Genetic sexing of animals with long gestation time benefits the management of captive populations. Here, X and Y chromosome-specific primers, based on equine gene sequencing data, were developed and tested on captive rhinoceroses (10 males, 20 females) representing four species (Diceros bicornis, Certaotherium simum simum, Rhinoceros unicornis, and Dicerorhinus sumatrensis). The Y chromosome-specific primer set targeted SRY (Sex-determining region Y), and amplified a 177-bp product following PCR of DNA extracted from males, but not females, of all species. A primer set based on the equine AMEL (Amelogenin) gene resulted in a 232-bp product following PCR of all rhinoceros species. These gene-specific primer sets were then evaluated for their ability to determine gender in cell-free DNA from rhinoceros serum. Modifications to the original extraction and PCR protocols were required to obtain sufficient DNA quantities from serum, and both DNA yield and PCR amplification were substantially reduced or absent following multiple freeze-thaw cycles of serum. When fresh serum from 14 pregnant rhinoceroses (ultimately bearing seven male and seven female calves), representing four species at different stages of gestation (Days 61-490), were probed in a PCR-based assay, an accuracy of 71% was achieved for male-specific gene detection of SRY, which improved to 100% by including a reamplification step into the protocol. Such early sex determination should be a valuable tool for current management practices as well as future assisted reproduction of rhinoceroses.

How can we reliably identify a taxon based on humeral morphology? Comparative morphology of desmostylian humeri

Desmostylia is a clade of marine mammals belonging to either Tethytheria or Perissodactyla. Rich fossil records of Desmostylia were found in the Oligocene to Miocene strata of the Northern Pacific Rim, especially in the northwestern region, which includes the Japanese archipelago. Fossils in many shapes and forms, including whole or partial skeletons, skulls, teeth, and fragmentary bones have been discovered from this region. Despite the prevalent availability of fossil records, detailed taxonomic identification based on fragmentary postcranial materials has been difficult owing to to our limited knowledge of the postcranial diagnostic features of many desmostylian taxa. In this study, I propose the utilization of diagnostic characters found in the humerus to identify desmostylian genus. These characters can be used to identify isolated desmostylian humeri at the genus level, contributing to a better understanding of the stratigraphic and geographic distributions of each genus.