Genetic diversity, phylogeny and conservation of the Javan rhinoceros (Rhinoceros sondaicus)

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

Background

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

Results

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

Conclusion

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

Supplementary Information

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

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.

Museum specimens tell the history of rhinoceroses

Rhinoceroses are facing an uncertain future due to human-induced population-size reductions. In this issue of Cell, Liu et al. analyze rhinoceros genomes to study their past and find support for an early divergence between African and Asian species, signatures of past hybridization, and evidence of loss of genetic diversity.

Comparative Phylogeography of Forest-Dependent Mammals Reveals Paleo-Forest Corridors throughout Sundaland

The evolutionary history of the colugo, a gliding arboreal mammal distributed throughout Sundaland, was influenced by the location of and connections between forest habitats. By comparing colugo phylogenetic patterns, species ecology, sample distributions, and times of divergence to those of other Sundaic taxa with different life-history traits and dispersal capabilities, we inferred the probable distribution of paleo-forest corridors and their influence on observed biogeographic patterns. We identified a consistent pattern of early diversification between east and west Bornean lineages in colugos, lesser mouse deer, and Sunda pangolins, but not in greater mouse deer. This deep east-west split within Borneo has not been commonly described in mammals. Colugos on West Borneo diverged from colugos in Peninsular Malaysia and Sumatra in the late Pliocene, however most other mammalian populations distributed across these same geographic regions diverged from a common ancestor more recently in the Pleistocene. Low genetic divergence between colugos on large landmasses and their neighboring satellite islands indicated that past forest distributions were recently much larger than present refugial distributions. Our analysis of colugo evolutionary history reconstructs Borneo as the most likely ancestral area of origin for Sunda colugos, and suggests that forests present during the middle Pliocene within the Sunda Shelf were more evergreen and contiguous, while forests were more fragmented, transient, seasonal, or with lower density canopies in the Pleistocene.

Recent mitochondrial lineage extinction in the critically endangered Javan rhinoceros

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

Interspecific variation in the limb long bones among modern rhinoceroses-extent and drivers

Among amniotes, numerous lineages are subject to an evolutionary trend toward body mass and size increases. Large terrestrial species may face important constraints linked to weight bearing, and the limb segments are particularly affected by such constraints due to their role in body support and locomotion. Such groups showing important limb modifications related to high body mass have been called "graviportal." Often considered graviportal, rhinoceroses are among the heaviest terrestrial mammals and are thus of particular interest to understand the limb modifications related to body mass and size increase. Here, we present a morphofunctional study of the shape variation of the limb long bones among the five living rhinos to understand how the shape may vary between these species in relation with body size, body mass and phylogeny. We used three dimensional geometric morphometrics and comparative analyses to quantify the shape variation. Our results indicate that the five species display important morphological differences depending on the considered bones. The humerus and the femur exhibit noticeable interspecific differences between African and Asiatic rhinos, associated with a significant effect of body mass. The radius and ulna are more strongly correlated with body mass. While the tibia exhibits shape variation both linked with phylogeny and body mass, the fibula displays the greatest intraspecific variation. We highlight three distinct morphotypes of bone shape, which appear in accordance with the phylogeny. The influence of body mass also appears unequally expressed on the different bones. Body mass increase among the five extant species is marked by an increase of the general robustness, more pronounced attachments for muscles and a development of medial parts of the bones. Our study underlines that the morphological features linked to body mass increase are not similar between rhinos and other heavy mammals such as elephants and hippos, suggesting that the weight bearing constraint can lead to different morphological responses.

A tale of two seas: contrasting patterns of population structure in the small-spotted catshark across Europe

Elasmobranchs represent important components of marine ecosystems, but they can be vulnerable to overexploitation. This has driven investigations into the population genetic structure of large-bodied pelagic sharks, but relatively little is known of population structure in smaller demersal taxa, which are perhaps more representative of the biodiversity of the group. This study explores spatial population genetic structure of the small-spotted catshark (Scyliorhinus canicula), across European seas. The results show significant genetic differences among most of the Mediterranean sample collections, but no significant structure among Atlantic shelf areas. The data suggest the Mediterranean populations are likely to have persisted in a stable and structured environment during Pleistocene sea-level changes. Conversely, the Northeast Atlantic populations would have experienced major changes in habitat availability during glacial cycles, driving patterns of population reduction and expansion. The data also provide evidence of male-biased dispersal and female philopatry over large spatial scales, implying complex sex-determined differences in the behaviour of elasmobranchs. On the basis of this evidence, we suggest that patterns of connectivity are determined by trends of past habitat stability that provides opportunity for local adaptation in species exhibiting philopatric behaviour, implying that resilience of populations to fisheries and other stressors may differ across the range of species.

Genetics and the last stand of the Sumatran rhinoceros Dicerorhinus sumatrensis

The Sumatran rhinoceros Dicerorhinus sumatrensis is on the brink of extinction. Although habitat loss and poaching were the reasons of the decline, today's reproductive isolation is the main threat to the survival of the species. Genetic studies have played an important role in identifying conservation priorities, including for rhinoceroses. However, for a species such as the Sumatran rhinoceros, where time is of the essence in preventing extinction, to what extent should genetic and geographical distances be taken into account in deciding the most urgently needed conservation interventions? We propose that the populations of Sumatra and Borneo be considered as a single management unit.

Extremely low genetic diversity indicating the endangered status of Ranodon sibiricus (Amphibia: Caudata) and implications for phylogeography

BACKGROUND

The Siberian salamander (Ranodon sibiricus), distributed in geographically isolated areas of Central Asia, is an ideal alpine species for studies of conservation and phylogeography. However, there are few data regarding the genetic diversity in R. sibiricus populations.

METHODOLOGY/PRINCIPAL FINDINGS

We used two genetic markers (mtDNA and microsatellites) to survey all six populations of R. sibiricus in China. Both of the markers revealed extreme genetic uniformity among these populations. There were only three haplotypes in the mtDNA, and the overall nucleotide diversity in the mtDNA was 0.00064, ranging from 0.00000 to 0.00091 for the six populations. Although we recovered 70 sequences containing microsatellite repeats, there were only two loci that displayed polymorphism. We used the approximate Bayesian computation (ABC) method to study the demographic history of the populations. This analysis suggested that the extant populations diverged from the ancestral population approximately 120 years ago and that the historical population size was much larger than the present population size; i.e., R. sibiricus has experienced dramatic population declines.

CONCLUSION/SIGNIFICANCE

Our findings suggest that the genetic diversity in the R. sibiricus populations is the lowest among all investigated amphibians. We conclude that the isolation of R. sibiricus populations occurred recently and was a result of recent human activity and/or climatic changes. The Pleistocene glaciation oscillations may have facilitated intraspecies genetic homogeneity rather than enhanced divergence. A low genomic evolutionary rate and elevated inbreeding frequency may have also contributed to the low genetic variation observed in this species. Our findings indicate the urgency of implementing a protection plan for this endangered species.