A comprehensive phylogeny of extant horses, rhinos and tapirs (Perissodactyla) through data combination

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.

Non-invasive assessment of fecal glucocorticoid, progesterone, and androgen metabolites and microbiome in free-ranging southern white rhinoceros (Ceratotherium simum simum) in South Africa

Increased poaching in northern South Africa has necessitated relocation of large numbers of southern white rhinoceros (Ceratotherium simum simum) to the Eastern Cape Province. The climate and grassland ecology of this province differ from that of northern South Africa which may impact the health of this species. This assessment of fecal steroid levels and microbiome in 10 free-ranging southern white rhinoceros in the Eastern Cape will provide insights into white rhinoceros physiology in this biome. Fecal steroid metabolites were analyzed using enzyme immunoassay (EIA) and ultra-performance convergence chromatography tandem mass spectrometry (UPC²-MS/MS). Fecal microbial composition was assessed via next generation sequencing. EIAs with antibodies raised against progesterone (P4; mouse monoclonal - CL425 clone), testosterone (T; rabbit polyclonal), corticosterone (B; sheep polyclonal) were utilized. Pregnant females had large quantities of fecal progesterone metabolites (FPMs) detected by CL425 EIA. Pregnant females also had native P4 and 11α-hydroxydihydroprogesterone (11αOHDHP4; 4-pregnen-11α-ol-3,20-dione) detected by UPC²-MS/MS but these concentrations were 1000-fold less than the concentrations of FPMs detected by the CL425 EIA. By contrast, non-pregnant females had FPM concentrations detected by CL425 EIA which were similar to native P4 and 11αOHDHP4 concentrations detected by UPC²-MS/MS. Mean fecal androgen metabolite (FAM) concentrations detected by the T EIA were similar between males and females. 11-ketoandrostenedione (11KA4) detected by UPC²-MS/MS was higher in females than males. However, there was no difference between males and females in the concentration of fecal glucocorticoid metabolites (FGMs) detected by the B EIA. Bacteroidia, followed by Clostridia, was the most abundant classes of fecal microbes. The unfiltered microbiome of females was more diverse than that of males. The core fecal microbiome of young rhinoceros had a higher observed species richness (Shannon diversity index, and Simpson diversity index) than that of old rhinoceros. In the alpha male, immobilization was associated with an increase in FGMs detected by 11-deoxycortisol (S) detected by UPC²-MS/MS coupled with decreased abundance of Spirochaetia. We detected substantially different FAM and FPM concentrations from those previously reported for both captive and wild white rhinoceros. Comparison of our UPC²-MS/MS and EIA results underscores the fact that most EIAs are highly cross reactive for many steroid metabolites. Our data also demonstrates a distinct effect of stress not only on FGMs but also on the fecal microbiome. This is the first non-invasive assessment of fecal steroid metabolites by UPC²-MS/MS and the fecal microbiome in wild white rhinoceros.

Experimental cross-species infection of donkeys with equine hepacivirus and analysis of host immune signatures

BACKGROUND

The Equine Hepacivirus (EqHV) is an equine-specific and liver-tropic virus belonging to the diverse genus of Hepaciviruses. It was recently found in a large donkey (Equus asinus) cohort with a similar seroprevalence (30%), but lower rate of RNA-positive animals (0.3%) compared to horses. These rare infection events indicate either a lack of adaptation to the new host or a predominantly acute course of infection.

METHODS

In order to analyze the susceptibility and the course of EqHV infection in donkeys, we inoculated two adult female donkeys and one control horse intravenously with purified EqHV from a naturally infected horse. Liver biopsies were taken before and after inoculation to study changes in the transcriptome.

RESULTS

Infection kinetics were similar between the equids. All animals were EqHV PCR-positive from day three. EqHV RNA-levels declined when the animals seroconverted and both donkeys cleared the virus from the blood by week 12. Infection did not have an impact on the clinical findings and no significant histopathological differences were seen. Blood biochemistry revealed a mild increase in GLDH at the time of seroconversion in horses, which was less pronounced in donkeys. Transcriptomic analysis revealed a distinct set of differentially expressed genes, including viral host factors and immune genes.

CONCLUSION

To summarize, our findings indicate that donkeys are a natural host of EqHV, due to the almost identical infection kinetics. The different immune responses do however suggest different mechanisms in reacting to hepaciviral infections.

Comparison of three hematocrit measurement methods in the southern white rhinoceros (Ceratotherium simum simum)

Background

Hematocrit (HCT) determination is an integral part of health and disease assessments in captive and wild white rhinoceroses. Several affordable automated hematology analyzers have been developed for in‐clinic and field use and have the advantage of being able to measure a large number of additional measurands. However, the accuracy of these analyzers for rhinoceros HCT measurements has not yet been investigated.

Objectives

We aimed to compare the HCT results generated by the EPOC portable analyzer system and the Abaxis VetScan HM5 with the gold standard of a manual packed cell volume (PCV) measured using the microhematocrit method.

Methods

Hematocrits were measured with the EPOC and the Abaxis VetScan HM5 (bovine setting) and compared with the PCVs of 69 white rhinoceros whole blood samples. Results were compared using Bland–Altman difference plots and Passing‐Bablok regression analysis. A total allowable analytical error of 10% was set as the performance goal.

Results

A significant positive bias, with a mean of 7.7% for the EPOC and 17.9% for the Abaxis, was found compared with the manual PCV method.

Conclusions

The allowable error goal of 10% was not exceeded with the EPOC analyzer. Although not analytically equivalent to the gold standard, the EPOC results could therefore be used as approximations in critical situations where manual measurements cannot be performed. The Abaxis exceeded this allowable error and overestimated HCTs in rhinoceroses. Therefore, method‐specific reference intervals should be used.

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.

On-site genetic analysis for species identification using lab-on-a-chip

This paper presents a microfluidic device capable of performing genetic analysis on dung samples to identify White Rhinoceros (Ceratotherium simum). The development of a microfluidic device, which can be used in the field, offers a portable and cost-effective solution for DNA analysis and species identification to aid conservation efforts. Optimization of the DNA extraction processes produced equivalent yields compared to conventional kit-based methods within just 5 minutes. The use of a color-changing loop-mediated isothermal amplification reaction for simultaneous detection of the cytochrome B sequence of C. simum enabled positive results to be obtained within as little as 30 minutes. Field testing was performed at Knowsley Safari to demonstrate real-world applicability of the microfluidic device for testing of biological samples.

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.

Embryos and embryonic stem cells from the white rhinoceros

The northern white rhinoceros (NWR, Ceratotherium simum cottoni) is the most endangered mammal in the world with only two females surviving. Here we adapt existing assisted reproduction techniques (ART) to fertilize Southern White Rhinoceros (SWR) oocytes with NWR spermatozoa. We show that rhinoceros oocytes can be repeatedly recovered from live SWR females by transrectal ovum pick-up, matured, fertilized by intracytoplasmic sperm injection and developed to the blastocyst stage in vitro. Next, we generate hybrid rhinoceros embryos in vitro using gametes of NWR and SWR. We also establish embryonic stem cell lines from the SWR blastocysts. Blastocysts are cryopreserved for later embryo transfer. Our results indicate that ART could be a viable strategy to rescue genes from the iconic, almost extinct, northern white rhinoceros and may also have broader impact if applied with similar success to other endangered large mammalian species.

The Brain of the Black (Diceros bicornis) and White (Ceratotherium simum) African Rhinoceroses: Morphology and Volumetrics from Magnetic Resonance Imaging

The morphology and volumetrics of the understudied brains of two iconic large terrestrial African mammals: the black (Diceros bicornis) and white (Ceratotherium simum) rhinoceroses are described. The black rhinoceros is typically solitary whereas the white rhinoceros is social, and both are members of the Perissodactyl order. Here, we provide descriptions of the surface of the brain of each rhinoceros. For both species, we use magnetic resonance images (MRI) to develop a description of the internal anatomy of the rhinoceros brain and to calculate the volume of the amygdala, cerebellum, corpus callosum, hippocampus, and ventricular system as well as to determine the gyrencephalic index. The morphology of both black and white rhinoceros brains is very similar to each other, although certain minor differences, seemingly related to diet, were noted, and both brains evince the general anatomy of the mammalian brain. The rhinoceros brains display no obvious neuroanatomical specializations in comparison to other mammals previously studied. In addition, the volumetric analyses indicate that the size of the various regions of the rhinoceros brain measured, as well as the extent of gyrification, are what would be predicted for a mammal with their brain mass when compared allometrically to previously published data. We conclude that the brains of the black and white rhinoceros exhibit a typically mammalian organization at a superficial level, but histological studies may reveal specializations of interest in relation to rhinoceros behavior.

Arterial Patterns of the Face in Perissodactyla

Considerable consistency in the arterial pattern of the head has been observed in species of Artiodactyla, but few studies have examined the order Perissodactyla. Here, we describe arteries supplying the intermandibular, mental, masseteric, buccal, labial, and nasal regions in eight perissodactylans, including representing of all families comprising this order. Observations were made on a total of 45 preparations of head arteries, obtained by injection of arteries with acetone-dissolved stained vinyl superchloride or stained latex LBS3060. In the Equidae species alone it was found that the facial artery descends from the linguofacial trunk. In tapirs and rhinos the facial artery branches off directly from the main arteries of the head. In tapirs alone it was found that the inferior alveolar artery gives off the buccal and sublingual arteries, and then extends into the mental artery. In the rhino a specific feature of the arterial pattern of the head was the exit of the occipital artery from the superficial temporal artery. In all equines studied, the transverse facial artery gave off a larger blood vessel to the masseter muscle and ran along the facial crest, while in tapirs and rhinos the transverse facial artery fanned out branches in the masseteric fossa. The variations observed can be considered in future studies on the origin of Perissodactyla. In this context, we note that the most similar patterns of exit and course of the facial, mental, transverse facial and infraorbital arteries exist in tapirs and rhinos (Ceratomorpha suborder), at least among the perissodactylans studied here. Anat Rec, 300:1529-1534, 2017. © 2017 Wiley Periodicals, Inc.

Middle Pleistocene protein sequences from the rhinoceros genus Stephanorhinus and the phylogeny of extant and extinct Middle/Late Pleistocene Rhinocerotidae

Background

Ancient protein sequences are increasingly used to elucidate the phylogenetic relationships between extinct and extant mammalian taxa. Here, we apply these recent developments to Middle Pleistocene bone specimens of the rhinoceros genus Stephanorhinus. No biomolecular sequence data is currently available for this genus, leaving phylogenetic hypotheses on its evolutionary relationships to extant and extinct rhinoceroses untested. Furthermore, recent phylogenies based on Rhinocerotidae (partial or complete) mitochondrial DNA sequences differ in the placement of the Sumatran rhinoceros (Dicerorhinus sumatrensis). Therefore, studies utilising ancient protein sequences from Middle Pleistocene contexts have the potential to provide further insights into the phylogenetic relationships between extant and extinct species, including Stephanorhinus and Dicerorhinus.

Methods

ZooMS screening (zooarchaeology by mass spectrometry) was performed on several Late and Middle Pleistocene specimens from the genus Stephanorhinus, subsequently followed by liquid chromatography-tandem mass spectrometry (LC-MS/MS) to obtain ancient protein sequences from a Middle Pleistocene Stephanorhinus specimen. We performed parallel analysis on a Late Pleistocene woolly rhinoceros specimen and extant species of rhinoceroses, resulting in the availability of protein sequence data for five extant species and two extinct genera. Phylogenetic analysis additionally included all extant Perissodactyla genera (Equus, Tapirus), and was conducted using Bayesian (MrBayes) and maximum-likelihood (RAxML) methods.

Results

Various ancient proteins were identified in both the Middle and Late Pleistocene rhinoceros samples. Protein degradation and proteome complexity are consistent with an endogenous origin of the identified proteins. Phylogenetic analysis of informative proteins resolved the Perissodactyla phylogeny in agreement with previous studies in regards to the placement of the families Equidae, Tapiridae, and Rhinocerotidae. Stephanorhinus is shown to be most closely related to the genera Coelodonta and Dicerorhinus. The protein sequence data further places the Sumatran rhino in a clade together with the genus Rhinoceros, opposed to forming a clade with the black and white rhinoceros species.

Discussion

The first biomolecular dataset available for Stephanorhinus places this genus together with the extinct genus Coelodonta and the extant genus Dicerorhinus. This is in agreement with morphological studies, although we are unable to resolve the order of divergence between these genera based on the protein sequences available. Our data supports the placement of the genus Dicerorhinus in a clade together with extant Rhinoceros species. Finally, the availability of protein sequence data for both extinct European rhinoceros genera allows future investigations into their geographic distribution and extinction chronologies.

Rewinding the process of mammalian extinction

With only three living individuals left on this planet, the northern white rhinoceros (Ceratotherium simum cottoni) could be considered doomed for extinction. It might still be possible, however, to rescue the (sub)species by combining novel stem cell and assisted reproductive technologies. To discuss the various practical options available to us, we convened a multidisciplinary meeting under the name "Conservation by Cellular Technologies." The outcome of this meeting and the proposed road map that, if successfully implemented, would ultimately lead to a self-sustaining population of an extremely endangered species are outlined here. The ideas discussed here, while centered on the northern white rhinoceros, are equally applicable, after proper adjustments, to other mammals on the brink of extinction. Through implementation of these ideas we hope to establish the foundation for reversal of some of the effects of what has been termed the sixth mass extinction event in the history of Earth, and the first anthropogenic one. Zoo Biol. 35:280-292, 2016. © 2016 The Authors. Zoo Biology published by Wiley Periodicals, Inc.

Natural killer cell receptor genes in the family Equidae: not only Ly49

Natural killer (NK) cells have important functions in immunity. NK recognition in mammals can be mediated through killer cell immunoglobulin-like receptors (KIR) and/or killer cell lectin-like Ly49 receptors. Genes encoding highly variable NK cell receptors (NKR) represent rapidly evolving genomic regions. No single conservative model of NKR genes was observed in mammals. Single-copy low polymorphic NKR genes present in one mammalian species may expand into highly polymorphic multigene families in other species. In contrast to other non-rodent mammals, multiple Ly49-like genes appear to exist in the horse, while no functional KIR genes were observed in this species. In this study, Ly49 and KIR were sought and their evolution was characterized in the entire family Equidae. Genomic sequences retrieved showed the presence of at least five highly conserved polymorphic Ly49 genes in horses, asses and zebras. These findings confirmed that the expansion of Ly49 occurred in the entire family. Several KIR-like sequences were also identified in the genome of Equids. Besides a previously identified non-functional KIR-Immunoglobulin-like transcript fusion gene (KIR-ILTA) and two putative pseudogenes, a KIR3DL-like sequence was analyzed. In contrast to previous observations made in the horse, the KIR3DL sequence, genomic organization and mRNA expression suggest that all Equids might produce a functional KIR receptor protein molecule with a single non-mutated immune tyrosine-based inhibition motif (ITIM) domain. No evidence for positive selection in the KIR3DL gene was found. Phylogenetic analysis including rhinoceros and tapir genomic DNA and deduced amino acid KIR-related sequences showed differences between families and even between species within the order Perissodactyla. The results suggest that the order Perissodactyla and its family Equidae with expanded Ly49 genes and with a potentially functional KIR gene may represent an interesting model for evolutionary biology of NKR genes.

Mitochondrial phylogenomics of modern and ancient equids

The genus Equus is richly represented in the fossil record, yet our understanding of taxonomic relationships within this genus remains limited. To estimate the phylogenetic relationships among modern horses, zebras, asses and donkeys, we generated the first data set including complete mitochondrial sequences from all seven extant lineages within the genus Equus. Bayesian and Maximum Likelihood phylogenetic inference confirms that zebras are monophyletic within the genus, and the Plains and Grevy’s zebras form a well-supported monophyletic group. Using ancient DNA techniques, we further characterize the complete mitochondrial genomes of three extinct equid lineages (the New World stilt-legged horses, NWSLH; the subgenus Sussemionus; and the Quagga, Equus quagga quagga). Comparisons with extant taxa confirm the NWSLH as being part of the caballines, and the Quagga and Plains zebras as being conspecific. However, the evolutionary relationships among the non-caballine lineages, including the now-extinct subgenus Sussemionus, remain unresolved, most likely due to extremely rapid radiation within this group. The closest living outgroups (rhinos and tapirs) were found to be too phylogenetically distant to calibrate reliable molecular clocks. Additional mitochondrial genome sequence data, including radiocarbon dated ancient equids, will be required before revisiting the exact timing of the lineage radiation leading up to modern equids, which for now were found to have possibly shared a common ancestor as far as up to 4 Million years ago (Mya).

A high density SNP array for the domestic horse and extant Perissodactyla: utility for association mapping, genetic diversity, and phylogeny studies

An equine SNP genotyping array was developed and evaluated on a panel of samples representing 14 domestic horse breeds and 18 evolutionarily related species. More than 54,000 polymorphic SNPs provided an average inter-SNP spacing of ∼43 kb. The mean minor allele frequency across domestic horse breeds was 0.23, and the number of polymorphic SNPs within breeds ranged from 43,287 to 52,085. Genome-wide linkage disequilibrium (LD) in most breeds declined rapidly over the first 50-100 kb and reached background levels within 1-2 Mb. The extent of LD and the level of inbreeding were highest in the Thoroughbred and lowest in the Mongolian and Quarter Horse. Multidimensional scaling (MDS) analyses demonstrated the tight grouping of individuals within most breeds, close proximity of related breeds, and less tight grouping in admixed breeds. The close relationship between the Przewalski's Horse and the domestic horse was demonstrated by pair-wise genetic distance and MDS. Genotyping of other Perissodactyla (zebras, asses, tapirs, and rhinoceros) was variably successful, with call rates and the number of polymorphic loci varying across taxa. Parsimony analysis placed the modern horse as sister taxa to Equus przewalski. The utility of the SNP array in genome-wide association was confirmed by mapping the known recessive chestnut coat color locus (MC1R) and defining a conserved haplotype of ∼750 kb across all breeds. These results demonstrate the high quality of this SNP genotyping resource, its usefulness in diverse genome analyses of the horse, and potential use in related species.