Search

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.

Much remains unknown about the reproductive physiology of southern white rhinoceros (SWR) and the effect of ovarian stimulation prior to ovum pickup (OPU) have not been fully elucidated. Granulosa cells (GC) provide valuable insight into follicle growth and oocyte maturation status. The goals of this study were to evaluate transcriptomic changes in GC from three stages of follicle development and to identify biomarkers possibly associated with follicular growth and maturation as a result of ovarian stimulation. GC collected from SWRs following OPU were assigned stages based upon follicle size. Total RNA was isolated, and cDNA libraries were prepared and sequenced on a NovaSeq 6000. All bioinformatics analyses were performed utilizing the Galaxy web platform. Reads were aligned to CerSimCot1.0, and the manual curation was performed with EquCab3.0. Overall, 39,455 transcripts (21,612 genes) were identified across follicle stages, and manual curation yielded a 61% increase in gene identification from the original annotation. Granulosa cells from preovulatory follicles expressed the highest number of unique transcripts. The following seven biomarkers were determined based upon cluster analysis and patterns of expression: COL1A1, JMY, FBXW11, NRG1, TMPO, MACIR and COL4A1. These data can be used to potentially evaluate the effects of different ovarian stimulation protocols on follicle dynamics, improve OPU results, and support conservation efforts in this species., (© 2024. The Author(s).)

Efforts to implement effective assisted reproductive technologies (ARTs) for the conservation of the northern white rhinoceros (NWR; Ceratotherium simum cottoni) to prevent its forthcoming extinction, could be supported by research conducted on the closely related southern white rhinoceros (SWR; Ceratotherium simum simum). Within the follicle, extracellular vesicles (EVs) play a fundamental role in the bidirectional communication facilitating the crucial transport of regulatory molecules such as microRNAs (miRNAs) that control follicular growth and oocyte development. This study aimed to elucidate the dynamics of EV-miRNAs in stage-dependent follicular fluid (FF) during SWR ovarian antral follicle development. Three distinct follicular stages were identified based on diameter: Growing (G; 11-17 mm), Dominant (D; 18-29 mm), and Pre-ovulatory (P; 30-34 mm). Isolated EVs from the aspirated FF of segmented follicle stages were used to identify EV-miRNAs previously known via subsequent annotation to all equine (Equus caballus; eca), bovine (Bos taurus; bta), and human (Homo sapiens; hsa) miRNAs. A total of 417 miRNAs were detected, with 231 being mutually expressed across all three stages, including eca-miR-148a and bta-miR-451 as the top highly expressed miRNAs. Distinct expression dynamics in miRNA abundance were observed across the three follicular stages, including 31 differentially expressed miRNAs that target various pathways related to follicular growth and development, with 13 miRNAs commonly appearing amidst two different comparisons. In conclusion, this pioneering study provides a comprehensive understanding of the stage-specific expression dynamics of FF EV-miRNAs in the SWR. These findings provide insights that may lead to novel approaches in enhancing ARTs to catalyze rhinoceros conservation efforts., (© The Author(s) 2024. Published by Oxford University Press behalf of Society for the Study of Reproduction.)

Although black rhinoceros Diceros bicornis are mostly solitary in the wild, the Hiroshima City Asa Zoological Park (Asa Zoo) has kept a family group together during the daytime, with good reproductive performance over five decades. Management procedures at the zoo include temporary single housing of the mother before and after giving birth, which facilitates maintenance of a compatible family group. We recorded intra-group spatial relationships for 4 years and 4 months, during which time an adult female reared two consecutive calves. During daytime she remained in an enclosure with her new calf, one or two older offspring, and an adult male, the sire of all her offspring. Proximity (within two adult body-lengths) scores between the mother and her two calves were especially high during the first year after birth, and only slightly lower for her older offspring. The adult male had the lowest proximity scores. The spatial relationships were visualized by applying multidimensional scaling (MDS) to the proximity scores. Mother and calves were plotted close to each other, with older offspring slightly farther apart on the two-dimensional MDS representation; the adult male was more distant from the other group members. These findings indicate clear follower-type characteristics in the mother-calf pair and also older immature offspring, albeit to a lesser degree. Although black rhinoceros are generally solitary in the wild, our results duplicate observations of some wild black rhinoceros groups containing an adult female, her calf, and an older immature, with adult males being largely solitary., (© 2024 Wiley Periodicals LLC.)

Objectives: The black rhinoceros (Diceros bicornis) is an endangered mammal for which a captive breeding program is part of the conservation effort. Black rhinos in zoo's often suffer from chronic infections and heamochromatosis. Furthermore, breeding is hampered by low male fertility. To aid a research project studying these topics, we sequenced and assembled the genome of a captive male black rhino using ONT sequencing data only., Data Description: This work produced over 100 Gb whole genome sequencing reads from whole blood. These were assembled into a 2.47 Gb draft genome consisting of 834 contigs with an N50 of 29.53 Mb. The genome annotation was lifted over from an available genome annotation for black rhino, which resulted in the retrieval of over 99% of gene features. This new genome assembly will be a valuable resource in for conservation genetic research in this species., (© 2024. The Author(s).)

The woolly rhinoceros ( Coelodonta antiquitatis ) is an iconic species of the Eurasian Pleistocene megafauna, which was abundant in Eurasia in the Pleistocene until its demise beginning approximately 10 000 years ago. Despite the early recovery of several specimens from well-known European archaeological sites, including its type specimen (Blumenbach 1799), no genomes of European populations were available so far, and all available genomic data originated exclusively from Siberian populations. Using coprolites of cave hyenas ( Crocuta crocuta spelea ) recovered from Middle Palaeolithic layers of two caves in Germany (Bockstein-Loch and Hohlenstein-Stadel), we isolated and enriched predator and prey DNA to assemble the first European woolly rhinoceros mitogenomes, in addition to cave hyena mitogenomes. Both coprolite samples produced copious sequences assigned to C. crocuta (27% and 59% mitogenome coverage, respectively) and woolly rhinoceros ( Coelodonta antiquitatis ; 27% and 81% coverage, respectively). The sequences suggested considerable DNA degradation, which may limit the conclusions to be drawn; however, the mitogenomes of European woolly rhinoceros are genetically distinct from the Siberian woolly rhinoceros, and analyses of the more complete mitogenome suggest a split of the populations potentially coinciding with the earliest fossil records of woolly rhinoceros in Europe.

The genus Bartonella (Hyphomicrobiales: Bartonellaceae) encompasses facultative intracellular α-proteobacteria that parasite erythrocytes and endothelial cells from a wide range of vertebrate hosts and can cause disease in animals and humans. Considering the large diversity of vertebrate species that may act as reservoirs and arthropod species that may be associated with Bartonella transmission, the exposure of animals and humans to these microorganisms is likely underestimated. The present study aimed to investigate the occurrence of Bartonella sp. in wild tapirs (Tapirus terrestris; Perissodactyla: Tapiridae) from two biomes in Brazil: Pantanal and Cerrado. Ninety-nine GPS-monitored wild tapirs were sampled in Pantanal (n = 61/99) and Cerrado (n = 38/99). A qPCR (quantitative real-time polymerase chain reaction) assay targeting the nuoG gene was used for the screening for Bartonella spp. DNA. Positive samples were additionally subjected to conventional PCR assays targeting five molecular markers (ribC, gltA, rpoB, groEL, ITS). Eight (8/99; 08,08%) animals were positive in the qPCR assay for Bartonella spp.: 7 from Cerrado (7/8; 87.5%) and 1 from Pantanal (1/8; 12.5%). The 5 Bartonella ribC sequences obtained from tapirs' blood samples grouped together with Bartonella henselae obtained from cats, humans, wild felids and Ctenocephalides felis (Siphonaptera: Pulicidae) fleas. To the best of author's knowledge, this is the first report of Bartonella sp. in Tapirus terrestris. This finding contributes to the understanding of the occurrence of B henselae in wild mammals from Brazil as well as expands the knowledge regarding the potential vector-borne pathogens that may affect wild tapis from Cerrado and Pantanal biomes., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Background: Woolly rhinoceros (Coelodonta antiquitatis) is a typical indicator of cold-stage climate that was widely distributed in Northern Hemisphere during the Middle-Late Pleistocene. Although a plethora of fossils have been excavated from Northern China, their phylogenetic status, intraspecific diversity and phylogeographical structure are still vague., Results: In the present study, we generated four mitogenomes from Late Pleistocene woolly rhinoceros in Northern China and compared them with published data. Bayesian and network analyses indicate that the analyzed individuals contain at least four maternal haplogroups, and Chinese samples fall in three of them. One of our samples belongs to a previously unidentified early diverging clade (haplogroup D), which separated from other woolly rhinoceros around 0.57 Ma (95% CI: 0.76-0.41 Ma). The timing of this clade's origin coincides with the first occurrence of woolly rhinoceros, which are thought to have evolved in Europe. Our other three samples cluster in haplogroup C, previously only identified from one specimen from Wrangel Island (ND030) and initially considered to be an isolated clade. Herein, our findings suggest that ND030 is likely descended from a northward dispersal of the individuals carrying haplogroup C from Northern China. Additionally, Chinese woolly rhinoceros specimens exhibit higher nucleotide diversity than those from Siberia., Conclusion: Our findings highlight Northern China as a possible refugium and a key evolution center of the Pleistocene woolly rhinoceros., (© 2023. The Author(s).)

Background: The harsh conditions of high-altitude environments are known to drive the evolution of physiological and morphological traits in endothermic animals. These conditions are expected to result in the adaptive evolution of protein coding genes encoded in mitochondrial genomes that are vital for the oxidative phosphorylation pathway. In this study, we formally tested for signatures of adaptive evolution on mitochondrial protein coding genes in Tapirus pinchaque and other odd-toed ungulates inhabiting high-elevation environments., Results: The AT-rich mitochondrial genome of T. pinchaque is 16,750 bp long. A phylomitogenomic analysis supports the monophyly of the genus Tapirus and families in the Perissodactyla. The ratio of non-synonymous to synonymous substitutions demonstrated that all mitochondrial genes undergo purifying selection in T. pinchaque and other odd ungulates living at high elevations. Over this negative background selection, Branch Models suggested that cox3 and nad6 might be undergoing stronger purifying selection than other mitochondrial protein coding genes. Furthermore, Site Models suggested that one and four sites in nad2 and nad5, respectively, could be experiencing positive selection. However, these results were supported by Likelihood Ratio Tests but not Bayesian Empirical Bayes posterior probabilities. Additional analyses (in DataMonkey) indicated a relaxation of selection strength in nad6, evidence of episodic diversifying selection in cob, and revealed episodic positive/diversifying selection signatures for two sites in nad1, and one site each in nad2 and nad4., Conclusion: The mitochondrial genome of T. pinchaque is an important genomic resource for conservation of this species and this study contributes to the understanding of adaptive evolution of mitochondrial protein coding genes in odd-toed ungulates inhabiting high-altitude environments., (© 2023. BioMed Central Ltd., part of Springer Nature.)

The black rhinoceros (Diceros bicornis L.) is a critically endangered species historically distributed across sub-Saharan Africa. Hunting and habitat disturbance have diminished both its numbers and distribution since the 19th century, but a poaching crisis in the late 20th century drove them to the brink of extinction. Genetic and genomic assessments can greatly increase our knowledge of the species and inform management strategies. However, when a species has been severely reduced, with the extirpation and artificial admixture of several populations, it is extremely challenging to obtain an accurate understanding of historic population structure and evolutionary history from extant samples. Therefore, we generated and analyzed whole genomes from 63 black rhinoceros museum specimens collected between 1775 and 1981. Results showed that the black rhinoceros could be genetically structured into six major historic populations (Central Africa, East Africa, Northwestern Africa, Northeastern Africa, Ruvuma, and Southern Africa) within which were nested four further subpopulations (Maasailand, southwestern, eastern rift, and northern rift), largely mirroring geography, with a punctuated north-south cline. However, we detected varying degrees of admixture among groups and found that several geographical barriers, most prominently the Zambezi River, drove population discontinuities. Genomic diversity was high in the middle of the range and decayed toward the periphery. This comprehensive historic portrait also allowed us to ascertain the ancestry of 20 resequenced genomes from extant populations. Lastly, using insights gained from this unique temporal data set, we suggest management strategies, some of which require urgent implementation, for the conservation of the remaining black rhinoceros diversity., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

The study of endobiotic ciliates from the intestines of various zoo mammals is important for revealing the influence of various factors on endobiont communities. In this paper we describe the species diversity of endobiotic ciliates from the faeces of the eastern black rhinoceros Diceros bicornis michaeli (male and two females, mother and daughter) kept in the zoo. Seven species of the ciliates were found, among them Rhinozeta rhinozeta, R. triciliata and Prototapirella clypeata were observed in the rhinos only. The other two, Monoposthium sp. and Triplumaria sp., according to their morphology should have been identified as new ciliate species. Successful transmission of the endobionts from parents to the young rhino in the zoo was demonstrated. R. rhinozeta was investigated with the using of the methods of the light and immunofluorescence microscopy and molecular phylogeny. The composition and phylogeny of the family Cycloposthiidae were discussed., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier GmbH. All rights reserved.)

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.

A Gram-negative strain, anaerobic, non-motile, non-spore-forming, rod-shaped bacterial strain named as NGMCC 1.200684  T was isolated from the fresh feces of rhinoceros in Beijing Zoo. Based on 16S rRNA gene sequences, phylogenetic analysis indicated that strain NGMCC 1.200684  T belonged to the genus Bacteroides and was most strongly related to the type strain of Bacteroides uniformis ATCC 8492  T (96.88%). The G + C content of the genomic DNA was determined to be 46.62%. Between strains NGMCC 1.200684  T and B. uniformis ATCC 8492  T , the average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) were 93.89 and 67.60%, respectively. Strain NGMCC 1.200684  T can produce acid from fermentation of several substrates, including glucose, mannitol, lactose, saccharose, maltose, salicin, xylose, cellobiose, mannose, raffinose, sorbitol, trehalose, D-galactose, and maltotriose. The major cellular fatty acids (> 10%) were identified as anteiso-C 15:0 , iso-C 15:0 , iso-C 14:0 , and iso-C 17:0 3-OH. The polar lipid profiles of strain NGMCC 1.200684  T were determined to contain diphosphatidyl glycerol, phosphatidylglycerol, phosphatidylethanolamine, three unknown phospholipids, and two unknown amino-phospholipids. Based on phenotypic, phylogenetic, and chemotaxonomic characteristics, a novel species of the genus Bacteroides, Bacteroides rhinocerotis sp. nov. is proposed. The type strain is NGMCC 1.200684  T (= CGMCC 1.18013  T  = JCM 35702  T )., (© 2023. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Two strains of Chryseobacterium identified from different experiments are proposed to represent new species. Strain WLa1L2M3 T was isolated from the digestive tract of an Oryctes rhinoceros beetle larva. Strain 09-1422 T was isolated from a cage housing the stick insect Eurycantha calcarata . Sequence analysis of the 16S rRNA and rpoB genes found both strains to be similar but not identical to other Chryseobacterium species. Whole-genome sequencing suggested the isolates represent new species, with average nucleotide identity values ranging from 74.6 to 80.5 %. Genome-to-genome distance calculations produced values below 25.3 %, and digital DNA-DNA hybridization values were 13.7-29.9 %, all suggesting they are distinct species. The genomic DNA G+C content of WLa1L2M3 T is approximately 32.53 %, and of 09-1422 T is approximately 35.89 %. The predominant cellular fatty acids of strain WLa1L2M3 T are C 15 : 0 iso, summed feature 9 (C 16 : 0 10OH or C 17 : 1 iso ω6 c ), C 17 : 0 iso 3OH, summed feature 3 (C 16 : 1 ω7 c and/or C 16 : 1 ω6c), C 15 : 0 iso 3OH, C 15 : 0 anteiso and C 13 : 0 iso, and those of strain 09-1422 T are C 15 : 0 iso, summed feature 3 (C 16 : 1 ω7c and/or C 16 : 1 ω6 c ), C 17 : 0 iso 3OH, C 15 : 0 anteiso, C 15 : 0 iso 3OH, C 16 : 1 ω7 c , C 17 : 0 2OH and C 18 : 0 . In addition, physiological and biochemical tests revealed phenotypic differences from related Chryseobacterium type strains. These cumulative data indicate that the two strains represent novel species of the genus Chryseobacterium for which the names Chryseobacterium oryctis sp. nov. and Chryseobacterium kimseyorum sp. nov. are proposed with WLa1L2M3 T (=BCRC 81350 T =JCM 35215 T =CIP 112035 T ) and 09-1422 T (=UCDFST 09-1422 T =BCRC 81359 T =CIP 112165 T ), as type strains, respectively.

Context: With two northern white rhinos (NWR) remaining, the continued existence of this species relies on studying their relative, the southern white rhino (SWR)., Aims: (1) Characterise gene expression in granulosa cells (GC) from SWR cumulus oocyte complexes (COCs) prior to (Pre-) and after (Post-) in vitro maturation (IVM), comparing culture media and oocytes from donors treated with or without gonadotropin stimulation prior to ovum recovery; and (2) evaluate COC glucose consumption in spent media., Methods: COCs were retrieved from four SWRs. Granulosa cells were collected before and after IVM in SDZ or IZW medium. Total RNA was evaluated by qPCR., Key Results: Oocyte maturation was greater in SDZ than IZW media. Expression of genes associated with follicle development increased in Pre-IVM GC. Six genes were differentially expressed in Post-IVM GC from stimulated compared to unstimulated donors. COCs from stimulated animals consumed more glucose. Fifty seven percent of oocytes in SDZ medium consumed all available glucose., Conclusions: Gene expression changed upon in vitro maturation and gonadotropin stimulation. Higher glucose availability might be needed during IVM., Implications: This is the first study examining GC gene expression and COC metabolic requirements in rhinoceros, which are critical aspects to optimise IVM of rhinoceros oocytes.

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., (© 2022. The Author(s).)

Background: An optimal starting point for relating genome function to organismal biology is a high-quality nuclear genome assembly, and long-read sequencing is revolutionizing the production of this genomic resource in insects. Despite this, nuclear genome assemblies have been under-represented for agricultural insect pests, particularly from the order Coleoptera. Here we present a de novo genome assembly and structural annotation for the coconut rhinoceros beetle, Oryctes rhinoceros (Coleoptera: Scarabaeidae), based on Oxford Nanopore Technologies (ONT) long-read data generated from a wild-caught female, as well as the assembly process that also led to the recovery of the complete circular genome assemblies of the beetle's mitochondrial genome and that of the biocontrol agent, Oryctes rhinoceros nudivirus (OrNV). As an invasive pest of palm trees, O. rhinoceros is undergoing an expansion in its range across the Pacific Islands, requiring new approaches to management that may include strategies facilitated by genome assembly and annotation., Results: High-quality DNA isolated from an adult female was used to create four ONT libraries that were sequenced using four MinION flow cells, producing a total of 27.2 Gb of high-quality long-read sequences. We employed an iterative assembly process and polishing with one lane of high-accuracy Illumina reads, obtaining a final size of the assembly of 377.36 Mb that had high contiguity (fragment N50 length = 12 Mb) and accuracy, as evidenced by the exceptionally high completeness of the benchmarked set of conserved single-copy orthologous genes (BUSCO completeness = 99.1%). These quality metrics place our assembly ahead of the published Coleopteran genomes, including that of an insect model, the red flour beetle (Tribolium castaneum). The structural annotation of the nuclear genome assembly contained a highly-accurate set of 16,371 protein-coding genes, with only 2.8% missing BUSCOs, and the expected number of non-coding RNAs. The number and structure of paralogous genes in a gene family like Sigma GST is lower than in another scarab beetle (Onthophagus taurus), but higher than in the red flour beetle (Tribolium castaneum), which suggests expansion of this GST class in Scarabaeidae. The quality of our gene models was also confirmed with the correct placement of O. rhinoceros among other members of the rhinoceros beetles (subfamily Dynastinae) in a phylogeny based on the sequences of 95 protein-coding genes in 373 beetle species from all major lineages of Coleoptera. Finally, we provide a list of 30 candidate dsRNA targets whose orthologs have been experimentally validated as highly effective targets for RNAi-based control of several beetles., Conclusions: The genomic resources produced in this study form a foundation for further functional genetic research and management programs that may inform the control and surveillance of O. rhinoceros populations, and we demonstrate the efficacy of de novo genome assembly using long-read ONT data from a single field-caught insect., (© 2022. The Author(s).)

The exaggerated horns of beetles are attractive models for studying the origin of novel traits and morphological evolution. Closely related species often differ profoundly in the size, number, and shape of their horns, and in the body region from which they extend. In addition, beetle horns exhibit exquisite nutrition-dependent phenotypic plasticity, leading to disproportionate growth of the horns in the largest, best-condition individuals and much smaller - even stunted - horn sizes in poor-condition individuals. These exciting phenomena in beetle horns have recently been revealed at the molecular level with the advent of next-generation sequencing. This section reviews the latest research on a horned beetle, the Japanese rhinoceros beetle Trypoxylus dichotomus, whose genome was recently sequenced., (Copyright © 2022 Elsevier Inc. All rights reserved.)

The northern white rhinoceros (NWR) is probably the earth's most endangered mammal. To rescue the functionally extinct species, we aim to employ induced pluripotent stem cells (iPSCs) to generate gametes and subsequently embryos in vitro. To elucidate the regulation of pluripotency and differentiation of NWR PSCs, we generated iPSCs from a deceased NWR female using episomal reprogramming, and observed surprising similarities to human PSCs. NWR iPSCs exhibit a broad differentiation potency into the three germ layers and trophoblast, and acquire a naïve-like state of pluripotency, which is pivotal to differentiate PSCs into primordial germ cells (PGCs). Naïve culturing conditions induced a similar expression profile of pluripotency related genes in NWR iPSCs and human ESCs. Furthermore, naïve-like NWR iPSCs displayed increased expression of naïve and PGC marker genes, and a higher integration propensity into developing mouse embryos. As the conversion process was aided by ectopic BCL2 expression, and we observed integration of reprogramming factors, the NWR iPSCs presented here are unsuitable for gamete production. However, the gained insights into the developmental potential of both primed and naïve-like NWR iPSCs are fundamental for in future PGC-specification in order to rescue the species from extinction using cryopreserved somatic cells., (© 2022. The Author(s).)

Large vertebrates are extremely sensitive to anthropogenic pressure, and their populations are declining fast. The white rhinoceros (Ceratotherium simum) is a paradigmatic case: this African megaherbivore has suffered a remarkable decline in the last 150 years due to human activities. Its subspecies, the northern (NWR) and the southern white rhinoceros (SWR), however, underwent opposite fates: the NWR vanished quickly, while the SWR recovered after the severe decline. Such demographic events are predicted to have an erosive effect at the genomic level, linked to the extirpation of diversity, and increased genetic drift and inbreeding. However, there is little empirical data available to directly reconstruct the subtleties of such processes in light of distinct demographic histories. Therefore, we generated a whole-genome, temporal data set consisting of 52 resequenced white rhinoceros genomes, representing both subspecies at two time windows: before and during/after the bottleneck. Our data reveal previously unknown population structure within both subspecies, as well as quantifiable genomic erosion. Genome-wide heterozygosity decreased significantly by 10% in the NWR and 36% in the SWR, and inbreeding coefficients rose significantly by 11% and 39%, respectively. Despite the remarkable loss of genomic diversity and recent inbreeding it suffered, the only surviving subspecies, the SWR, does not show a significant accumulation of genetic load compared to its historical counterpart. Our data provide empirical support for predictions about the genomic consequences of shrinking populations, and our findings have the potential to inform the conservation efforts of the remaining white rhinoceroses., (© 2021 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

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., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Current management models for many endangered species focus primarily on demographic recovery, often ignoring their intrinsic ecological requirements. Across the protected area network of southern Africa, most southern white rhinoceros are managed in populations of less than 50 individuals, experiencing restricted dispersal opportunities, and limited breeding male numbers due to their exclusive home range requirements. In the absence of information on the breeding structure of these populations, poor management decisions may require females to either forego a breeding opportunity or select to inbreed with close relatives. Here, we use a combination of social pedigree data together with genetic analyses to reconstruct the parentage of all 28 offspring produced in a 5-year period in a managed free-ranging southern white rhinoceros population. During this period, all breeding females (founders and first-generation daughters) had access to both a founder male (father to most of the daughters) and two recently introduced inexperienced males. We report that while founder females were more likely to breed with the founder male, their daughters, in contrast, were more likely to breed with the introduced males, thus avoiding inbreeding. However, we also found evidence of father-daughter inbreeding in this population, and contend that in the absence of choice, rather than forego a breeding opportunity, female white rhinoceros will inbreed with their fathers. We argue that to effectively conserve the southern white rhinoceros, managers need to understand the breeding structure of these small populations, particularly in terms of parentage and kinship., (© The American Genetic Association. 2021. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

As one of the largest land mammals, the origin and evolution of the giant rhino Paraceratherium bugtiense in Pakistan have been unclear. We report a new species Paraceratherium linxiaense sp. nov. from northwestern China with an age of 26.5 Ma. Morphology and phylogeny reveal that P. linxiaense is the highly derived species of the genus Paraceratherium, and its clade with P. lepidum has a tight relationship to P. bugtiense. Based on the paleogeographical literature, P. bugtiense represents a range expansion of Paraceratherium from Central Asia via the Tibetan region. By the late Oligocene, P. lepidum and P. linxiaense were found in the north side of the Tibetan Plateau. The Tibetan region likely hosted some areas with low elevation, possibly under 2000 m during Oligocene, and the lineage of giant rhinos could have dispersed freely along the eastern coast of the Tethys Ocean and perhaps through some lowlands of this region.

Modes of mammalian reproduction are diverse and not always conserved among related species. Progesterone is universally required to supports pregnancy but sites of synthesis and metabolic pathways vary widely. The steroid metabolome of mid-to late gestation was characterized, focusing on 5α-reduced pregnanes in species representing the Perissodactyla, Cetartiodactyla and Carnivora using mass spectrometry. Metabolomes and steroidogenic enzyme ortholog sequences were used in heirarchial analyses. Steroid metabolite profiles were similar within orders, whales within cetartiodactyls for instance, but with notable exceptions such as rhinoceros clustering with goats, and tapirs with pigs. Steroidogenic enzyme sequence clustering reflected expected evolutionary relationships but once again with exceptions. Human sequences (expected outgroups) clustered with perissodactyl CYP11A1, CYP17A1 and SRD5A1 gene orthologues, forming outgroups only for HSD17B1 and SRD5A2. Spotted hyena CYP19A1 clustered within the Perissodactyla, between rhinoceros and equid orthologues, whereas CYP17A1 clustered within the Carnivora. This variability highlights the random adoption of divergent physiological strategies as pregnancy evolved among genetically similar species., (Copyright © 2021 Elsevier B.V. All rights reserved.)

The Greater one-horned (GoH) rhinoceros is one of the most charismatic endemic megaherbivores of the Indian subcontinent. Threatened by poaching, habitat loss and disease, the species is found only in small areas of its historical distribution. Increasing demands for rhino horns in chinese traditional medicine has put the existing population under continuing threat, and large profits and low conviction rates make poaching difficult to contain. DNA forensics such as the RhoDIS-Africa program has helped in combating illegal rhino horn trade, but the approach is yet to be optimised for Indian GoH rhinoceros. Here we followed the International Society for Forensic Genetics (ISFG) guidelines to establish a 14 dinucleotide microsatellite panel for Indian GoH rhinoceros DNA profiling. Selected from a large initial pool (n = 34), the microsatellite markers showed high polymorphism, stable peak characteristics, consistent allele calls and produced precise, reproducible genotypes from different types of rhino samples. The panel also showed low genotyping error and produced high statistical power during individual identification (P IDsibs value of 1.2*10 -4 ). As part of the official RhoDIS-India program, we used this panel to match poached rhino carcass with seized contraband as scientific evidence in court procedure. This program now moves to generate detailed allele-frequency maps of all GoH rhinoceros populations in India and Nepal for development of a genetic database and identification of poaching hotspots and trade routes across the subcontinent and beyond., (Copyright © 2021 Elsevier B.V. All rights reserved.)

Small populations are often exposed to high inbreeding and mutational load that can increase the risk of extinction. The Sumatran rhinoceros was widespread in Southeast Asia, but is now restricted to small and isolated populations on Sumatra and Borneo, and most likely extinct on the Malay Peninsula. Here, we analyse 5 historical and 16 modern genomes from these populations to investigate the genomic consequences of the recent decline, such as increased inbreeding and mutational load. We find that the Malay Peninsula population experienced increased inbreeding shortly before extirpation, which possibly was accompanied by purging. The populations on Sumatra and Borneo instead show low inbreeding, but high mutational load. The currently small population sizes may thus in the near future lead to inbreeding depression. Moreover, we find little evidence for differences in local adaptation among populations, suggesting that future inbreeding depression could potentially be mitigated by assisted gene flow among populations.

Objective: The Sumatran rhinoceros is critically endangered, with fewer than 100 individuals surviving across its current range. Accurate census estimates of the remaining populations are essential for development and implementation of conservation plans. In order to enable molecular censusing, we here develop microsatellite markers with amplicon sizes of short length, appropriate for non-invasive fecal sampling., Results: Due to limited sample quantity and potential lack of genome-wide diversity, Illumina sequence reads were generated from two Sumatran rhinoceros samples. Genomic screening identified reads with short tandem repeats and loci that were polymorphic within the dataset. Twenty-nine novel polymorphic microsatellite markers were characterized (A = 2.4; H O  = 0.30). These were sufficient to distinguish among individuals (P ID  < 0.0001), and to distinguish among siblings (P ID(sib)  < 0.0001). Among rhinos in Indonesia, almost all markers were established as polymorphic and effective for genotyping DNA from fecal samples. Notably, the markers amplified and displayed microsatellite polymorphisms using DNA extracted from 11 fecal samples collected non-invasively from wild Sumatran rhinoceros. These microsatellite markers provide an important resource for a census and genetic studies of wild Sumatran rhinos.

Extinction rates are rising, and current conservation technologies may not be adequate for reducing species losses. Future conservation efforts may be aided by the generation of induced pluripotent stem cells (iPSCs) from highly endangered species. Generation of a set of iPSCs from multiple members of a species can capture some of the dwindling genetic diversity of a disappearing species. We generated iPSCs from fibroblasts cryopreserved in the Frozen Zoo ® : nine genetically diverse individuals of the functionally extinct northern white rhinoceros ( Ceratotherium simum cottoni ) and two from the closely related southern white rhinoceros ( Ceratotherium simum simum ). We used a nonintegrating Sendai virus reprogramming method and developed analyses to confirm the cells' pluripotency and differentiation potential. This work is the first step of a long-term interdisciplinary plan to apply assisted reproduction techniques to the conservation of this highly endangered species. Advances in iPSC differentiation may enable generation of gametes in vitro from deceased and nonreproductive individuals that could be used to repopulate the species.

Local wild bovids have been determined to be important prey on the northeastern Tibetan Plateau (NETP), where hunting game was a major subsistence strategy until the late Neolithic, when farming lifestyles dominated in the neighboring Loess Plateau. However, the species affiliation and population ecology of these prehistoric wild bovids in the prehistoric NETP remain unknown. Ancient DNA (aDNA) analysis is highly informative in decoding this puzzle. Here, we applied aDNA analysis to fragmented bovid and rhinoceros specimens dating ∼5,200 y B.P. from the Neolithic site of Shannashuzha located in the marginal area of the NETP. Utilizing both whole genomes and mitochondrial DNA, our results demonstrate that the range of the present-day tropical gaur ( Bos gaurus ) extended as far north as the margins of the NETP during the late Neolithic from ∼29°N to ∼34°N. Furthermore, comparative analysis with zooarchaeological and paleoclimatic evidence indicated that a high summer temperature in the late Neolithic might have facilitated the northward expansion of tropical animals (at least gaur and Sumatran-like rhinoceros) to the NETP. This enriched the diversity of wildlife, thus providing abundant hunting resources for humans and facilitating the exploration of the Tibetan Plateau as one of the last habitats for hunting game in East Asia., Competing Interests: The authors declare no competing interest., (Copyright © 2020 the Author(s). Published by PNAS.)

Africa's black (Diceros bicornis) and white (Ceratotherium simum) rhinoceros are closely related sister-taxa that evolved highly divergent obligate browsing and grazing feeding strategies. Although their precursor species Diceros praecox and Ceratotherium mauritanicum appear in the fossil record ∼5.2 Ma, by 4 Ma both were still mixed feeders, and were even spatiotemporally sympatric at several Pliocene sites in what is today Africa's Rift Valley. Here, we ask whether or not D. praecox and C. mauritanicum were reproductively isolated when they came into Pliocene secondary contact. We sequenced and de novo assembled the first annotated black rhinoceros reference genome and compared it with available genomes of other black and white rhinoceros. We show that ancestral gene flow between D. praecox and C. mauritanicum ceased sometime between 3.3 and 4.1 Ma, despite conventional methods for the detection of gene flow from whole genome data returning false positive signatures of recent interspecific migration due to incomplete lineage sorting. We propose that ongoing Pliocene genetic exchange, for up to 2 My after initial divergence, could have potentially hindered the development of obligate feeding strategies until both species were fully reproductively isolated, but that the more severe and shifting paleoclimate of the early Pleistocene was likely the ultimate driver of ecological specialization in African rhinoceros., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Ancient DNA has significantly improved our understanding of the evolution and population history of extinct megafauna. However, few studies have used complete ancient genomes to examine species responses to climate change prior to extinction. The woolly rhinoceros (Coelodonta antiquitatis) was a cold-adapted megaherbivore widely distributed across northern Eurasia during the Late Pleistocene and became extinct approximately 14 thousand years before present (ka BP). While humans and climate change have been proposed as potential causes of extinction [1-3], knowledge is limited on how the woolly rhinoceros was impacted by human arrival and climatic fluctuations [2]. Here, we use one complete nuclear genome and 14 mitogenomes to investigate the demographic history of woolly rhinoceros leading up to its extinction. Unlike other northern megafauna, the effective population size of woolly rhinoceros likely increased at 29.7 ka BP and subsequently remained stable until close to the species' extinction. Analysis of the nuclear genome from a ∼18.5-ka-old specimen did not indicate any increased inbreeding or reduced genetic diversity, suggesting that the population size remained steady for more than 13 ka following the arrival of humans [4]. The population contraction leading to extinction of the woolly rhinoceros may have thus been sudden and mostly driven by rapid warming in the Bølling-Allerød interstadial. Furthermore, we identify woolly rhinoceros-specific adaptations to arctic climate, similar to those of the woolly mammoth. This study highlights how species respond differently to climatic fluctuations and further illustrates the potential of palaeogenomics to study the evolutionary history of extinct species., Competing Interests: Declaration of Interests The authors declare no competing interests., (Copyright © 2020 The Author(s). Published by Elsevier Inc. All rights reserved.)

A molecular sexing method by polymerase chain reaction (PCR) amplification of a portion of the sex-determining region Y (SRY) and the zinc finger (ZF) gene, as well as six equine Y-chromosome-specific microsatellite markers, were tested in the Malayan tapir (Tapirus indicus). While the microsatellite markers did not yield any male-specific amplicons for sex-typing, the SRY/ZF marker system produced reliable molecular sexing results by accurately sex-typing 31 reference Malayan tapirs, using whole blood, dried blood spot (DBS), or tissue samples as materials for DNA extraction. The marker system was also tested on 16 faecal samples, and the results were in general consistent with the pre-determined sexes of the animals, despite some amplification failures. A preliminary estimation of wild Malayan tapir population sex ratio was estimated from the Wildlife Genomic Resource Bank (WGRB) database of the Malaysian Department of Wildlife and National Parks (PERHILITAN), zoos, and the Sungai Dusun Wildlife Conservation Centre (WCC), as well as from the results of molecular sexing 12 samples of unknown sex. The overall sex ratio favoured females, but the deviation from parity was statistically not significant when tested using the binomial test (p > 0.05), which may be due to reduced statistical power caused by small sample sizes.

The sequencing of ancient DNA has enabled the reconstruction of speciation, migration and admixture events for extinct taxa 1 . However, the irreversible post-mortem degradation 2 of ancient DNA has so far limited its recovery-outside permafrost areas-to specimens that are not older than approximately 0.5 million years (Myr) 3 . By contrast, tandem mass spectrometry has enabled the sequencing of approximately 1.5-Myr-old collagen type I 4 , and suggested the presence of protein residues in fossils of the Cretaceous period 5 -although with limited phylogenetic use 6 . In the absence of molecular evidence, the speciation of several extinct species of the Early and Middle Pleistocene epoch remains contentious. Here we address the phylogenetic relationships of the Eurasian Rhinocerotidae of the Pleistocene epoch 7-9 , using the proteome of dental enamel from a Stephanorhinus tooth that is approximately 1.77-Myr old, recovered from the archaeological site of Dmanisi (South Caucasus, Georgia) 10 . Molecular phylogenetic analyses place this Stephanorhinus as a sister group to the clade formed by the woolly rhinoceros (Coelodonta antiquitatis) and Merck's rhinoceros (Stephanorhinus kirchbergensis). We show that Coelodonta evolved from an early Stephanorhinus lineage, and that this latter genus includes at least two distinct evolutionary lines. The genus Stephanorhinus is therefore currently paraphyletic, and its systematic revision is needed. We demonstrate that sequencing the proteome of Early Pleistocene dental enamel overcomes the limitations of phylogenetic inference based on ancient collagen or DNA. Our approach also provides additional information about the sex and taxonomic assignment of other specimens from Dmanisi. Our findings reveal that proteomic investigation of ancient dental enamel-which is the hardest tissue in vertebrates 11 , and is highly abundant in the fossil record-can push the reconstruction of molecular evolution further back into the Early Pleistocene epoch, beyond the currently known limits of ancient DNA preservation.

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.

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., (Copyright © 2018 American Society for Microbiology.)

The white rhinoceros ( Ceratotherium simum ) has a discontinuous African distribution, which is limited by the extent of sub-Saharan grasslands. The southern population (SWR) declined to its lowest number around the turn of the nineteenth century, but recovered to become the world's most numerous rhinoceros. In contrast, the northern population (NWR) was common during much of the twentieth century, declining rapidly since the 1970s, and now only two post-reproductive individuals remain. Despite this species's conservation status, it lacks a genetic assessment of its demographic history. We therefore sampled 232 individuals from extant and museum sources and analysed ten microsatellite loci and the mtDNA control region. Both marker types reliably partitioned the species into SWR and NWR, with moderate nuclear genetic diversity and only three mtDNA haplotypes for the species, including historical samples. We detected ancient interglacial demographic declines in both populations. Both populations may also have been affected by recent declines associated with the colonial expansion for the SWR, and with the much earlier Bantu migrations for the NWR. Finally, we detected post-divergence secondary contact between NWR and SWR, possibly occurring as recently as the last glacial maximum. These results suggest the species was subjected to regular periods of fragmentation and low genetic diversity, which may have been replenished upon secondary contact during glacial periods. The species's current situation thus reflects prehistoric declines that were exacerbated by anthropogenic pressure associated with the rise of late Holocene technological advancement in Africa. Importantly, secondary contact suggests a potentially positive outcome for a hybrid rescue conservation strategy, although further genome-wide data are desirable to corroborate these results., (© 2018 The Author(s).)

The Sumatran rhinoceros (Dicerorhinus sumatrensis), once widespread across Southeast Asia, now consists of as few as 30 individuals within Sumatra and Borneo. To aid in conservation planning, we sequenced 218 bp of control region mitochondrial (mt) DNA, identifying 17 distinct mitochondrial haplotypes across modern (N = 13) and museum (N = 26) samples. Museum specimens from Laos and Myanmar had divergent mtDNA, consistent with the placement of western mainland rhinos into the distinct subspecies D. s. lasiotis (presumed extinct). Haplotypes from Bornean rhinos were highly diverse, but dissimilar from those of other regions, supporting the distinctiveness of the subspecies D. s. harrissoni. Rhinos from Sumatra and Peninsular Malaysia shared mtDNA haplotypes, consistent with their traditional placement into a single subspecies D. s sumatrensis. Modern samples of D. s. sumatrensis were genotyped at 18 microsatellite loci. Rhinos within Sumatra formed 2 sub-populations, likely separated by the Barisan Mountains, though with only modest genetic differentiation between them. There are so few remaining Sumatran rhinoceros that separate management strategies for subspecies or subpopulations may not be viable, while each surviving rhino pedigree is likely to retain alleles found in no other individuals. Given the low population size and low reproductive potential of Sumatran rhinos, rapid genetic erosion is inevitable, though an under-appreciated concern is the potential for fixation of harmful genetic variants. Both concerns underscore 2 overriding priorities for the species: 1) translocation of wild rhinos to ex situ facilities, and 2) collection and storage of gametes and cell lines from every surviving captive and wild individual.

Rhinoceros (rhinos) have suffered a dramatic increase in poaching over the past decade due to the growing demand for rhino horn products in Asia. One way to reverse this trend is to enhance enforcement and intelligence gathering tools used for species identification of horns, in particular making them fast, inexpensive and accurate. Traditionally, species identification tests are based on DNA sequence data, which, depending on laboratory resources, can be either time or cost prohibitive. This study presents a rapid rhino species identification test, utilizing species-specific primers within the cytochrome b gene multiplexed in a single reaction, with a presumptive species identification based on the length of the resultant amplicon. This multiplex PCR assay can provide a presumptive species identification result in less than 24 hours. Sequence-based definitive testing can be conducted if/when required (e.g. court purposes). This work also presents an actual casework scenario in which the presumptive test was successfully utlitised, in concert with sequence-based definitive testing. The test was carried out on seized suspected rhino horns tested at the Institute of Ecology and Biological Resources, the CITES mandated laboratory in Vietnam, a country that is known to be a major source of demand for rhino horns. This test represents the basis for which future 'rapid species identification tests' can be trialed., Competing Interests: The authors have declared that no competing interests exist.

The critically endangered northern white rhinoceros is believed to be extinct in the wild, with the recent death of the last male leaving only two remaining individuals in captivity. Its extinction would appear inevitable, but the development of advanced cell and reproductive technologies such as cloning by nuclear transfer and the artificial production of gametes via stem cells differentiation offer a second chance for its survival. In this work, we analyzed genome-wide levels of genetic diversity, inbreeding, population history, and demography of the white rhinoceros sequenced from cryopreserved somatic cells, with the goal of informing how genetically valuable individuals could be used in future efforts toward the genetic rescue of the northern white rhinoceros. We present the first sequenced genomes of the northern white rhinoceros, which show relatively high levels of heterozygosity and an average genetic divergence of 0.1% compared with the southern subspecies. The two white rhinoceros subspecies appear to be closely related, with low genetic admixture and a divergent time <80,000 yr ago. Inbreeding, as measured by runs of homozygosity, appears slightly higher in the southern than the northern white rhinoceros. This work demonstrates the value of the northern white rhinoceros cryopreserved genetic material as a potential gene pool for saving this subspecies from extinction., (© 2018 Tunstall et al.; Published by Cold Spring Harbor Laboratory Press.)

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., (© 2018 Wiley Periodicals, Inc.)

The vertebrate extinction rate over the past century is approximately 22-100 times greater than background extinction rates [1], and large mammals are particularly at risk [2, 3]. Quaternary megafaunal extinctions have been attributed to climate change [4], overexploitation [5], or a combination of the two [6]. Rhinoceroses (Family: Rhinocerotidae) have a rich fossil history replete with iconic examples of climate-induced extinctions [7], but current pressures threaten to eliminate this group entirely. The Sumatran rhinoceros (Dicerorhinus sumatrensis) is among the most imperiled mammals on earth. The 2011 population was estimated at ≤216 wild individuals [8], and currently the species is extirpated, or nearly so, throughout the majority of its former range [8-12]. Understanding demographic history is important in placing current population status into a broader ecological and evolutionary context. Analysis of the Sumatran rhinoceros genome reveals extreme changes in effective population size throughout the Pleistocene. Population expansion during the early to middle Pleistocene was followed by decline. Ecological niche modeling indicated that changing climate most likely played a role in the decline of the Sumatran rhinoceros, as less suitable habitat on an emergent Sundaland corridor isolated Sumatran rhinoceros populations. By the end of the Pleistocene, the Sundaland corridor was submerged, and populations were fragmented and consequently reduced to low Holocene levels from which they would never recover. Past events denuded the Sumatran rhinoceros of genetic diversity through population decline, fragmentation, or some combination of the two and most likely made the species even more susceptible to later exploitation and habitat loss. VIDEO ABSTRACT., (Copyright © 2017 The Authors. Published by Elsevier Ltd.. All rights reserved.)

We presented the powerful techniques for species identification using the short amplicon of mitochondrial cytochrome b gene sequence. Two faecal samples and one single hair sample of the Asian tapir were tested using the new cytochrome b primers. The results showed a high sequence similarity with the mainland Asian tapir group. The comparative sequence analysis of the reserved wild mammals in Thailand and the other endangered mammal species from Southeast Asia comprehensibly verified the potential of our novel primers. The forward and reverse primers were 94.2 and 93.2%, respectively, by the average value of the sequence identity among 77 species sequences, and the overall mean distance was 35.9%. This development technique could provide rapid, simple, and reliable tools for species confirmation. Especially, it could recognize the problematic biological specimens contained less DNA material from illegal products and assist with wildlife crime investigation of threatened species and related forensic casework.

Rhinoceros (rhino) numbers have dwindled substantially over the past century. As a result, three of the five species are now considered to be critically endangered, one species is vulnerable and one species is near-threatened. Poaching has increased dramatically over the past decade due to a growing demand for rhino horn products, primarily in Asia. Improved wildlife forensic techniques, such as validated tests for species identification of seized horns, are critical to aid current enforcement and prosecution efforts and provide a deterrent to future rhino horn trafficking. Here, we present an internationally standardized species identification test based on a 230 base pair cytochrome-b region. This test improves on previous nested PCR protocols and can be used for the discrimination of samples with <20pg of template DNA, thus suitable for DNA extracted from horn products. The assay was designed to amplify water buffalo samples, a common 'rhino horn' substitute, but to exclude human DNA, a common contaminant. Phylogenetic analyses using this partial cytochrome-b region resolved the five extant rhino species. Testing successfully returned a sequence and correct identification for all of the known rhino horn samples and vouchered rhino samples from museum and zoo collections, and provided species level identification for 47 out of 52 unknown samples from seizures. Validation and standardization was carried out across five different laboratories, in four different countries, demonstrating it to be an effective and reproducible test, robust to inter laboratory variation in equipment and consumables (such as PCR reagents). This is one of the first species identification tests to be internationally standardized to produce data for evidential proceedings and the first published validated test for rhinos, one of the flagship species groups of the illegal wildlife trade and for which forensic tools are urgently required. This study serves as a model for how species identification tests should be standardized and disseminated for wildlife forensic testing., (Copyright © 2017 The Authors. Published by Elsevier B.V. All rights reserved.)

The Asian tapir (Tapirus indicus) has been classified as Endangered on the IUCN Red List of Threatened Species (2008). Genetic diversity data provide important information for the management of captive breeding and conservation of this species. We analyzed mitochondrial control region (CR) sequences from 37 captive Asian tapirs in Thailand. Multiple alignments of the full-length CR sequences sized 1268 bp comprised three domains as described in other mammal species. Analysis of 16 parsimony-informative variable sites revealed 11 haplotypes. Furthermore, the phylogenetic analysis using median-joining network clearly showed three clades correlated with our earlier cytochrome b gene study in this endangered species. The repetitive motif is located between first and second conserved sequence blocks, similar to the Brazilian tapir. The highest polymorphic site was located in the extended termination associated sequences domain. The results could be applied for future genetic management based in captivity and wild that shows stable populations.

The black rhinoceros is again on the verge of extinction due to unsustainable poaching in its native range. Despite a wide historic distribution, the black rhinoceros was traditionally thought of as depauperate in genetic variation, and with very little known about its evolutionary history. This knowledge gap has hampered conservation efforts because hunting has dramatically reduced the species' once continuous distribution, leaving five surviving gene pools of unknown genetic affinity. Here we examined the range-wide genetic structure of historic and modern populations using the largest and most geographically representative sample of black rhinoceroses ever assembled. Using both mitochondrial and nuclear datasets, we described a staggering loss of 69% of the species' mitochondrial genetic variation, including the most ancestral lineages that are now absent from modern populations. Genetically unique populations in countries such as Nigeria, Cameroon, Chad, Eritrea, Ethiopia, Somalia, Mozambique, Malawi and Angola no longer exist. We found that the historic range of the West African subspecies (D. b. longipes), declared extinct in 2011, extends into southern Kenya, where a handful of individuals survive in the Masai Mara. We also identify conservation units that will help maintain evolutionary potential. Our results suggest a complete re-evaluation of current conservation management paradigms for the black rhinoceros., Competing Interests: The authors declare no competing financial interests.

The objective of this study was to design an artificial insemination (AI) protocol using cryopreserved spermatozoa to obtain pregnancies in captive Indian rhinoceroses (Rhinoceros unicornis). Four methods developed varied by timing and approach, as follows; Method 1: females (n=2) were inseminated pre- and post-ovulation under general anesthesia, Method 2: females (n=2) were inseminated pre-ovulation without anesthetic via endoscopy, Method 3: females (n=1) were inseminated pre-ovulation without anesthetic via manual insertion of an insemination catheter, Method 4: females (n=2) were inseminated same as Method 3 with the addition of standing sedation. Semen deposition site varied as a result of changes in AI technology and experience. All females conceived following intrauterine AI using three methods. Four pregnancies (n=3 females) produced via Method 3 and 4 resulted in term births (n=2 male calves, n=2 female calves) at 481.8±12.8days post-AI. Unfortunately, two early pregnancy losses were documented in a fourth female conceiving via Method 2. Pregnancy rates were 0%, 22%, 17%, and 50% for Method 1-4, respectively. Method 3 and 4 rates improved to 29% and 67%, respectively when accounting for AI's conducted only on ovulatory estrous cycles. Spermatozoa (n=5 males) were cryopreserved 0.3-9.3 y prior to successful AI procedures. The lowest dose of frozen-thawed sperm resulting in conception was 500×10(6) motile sperm. Mean time from AI to ovulation in conceptive and non-conceptive cycles was 26±11.8h and 66±80.7h, respectively., (Copyright © 2016 Elsevier B.V. All rights reserved.)