Module 0

Evolution & Phylogenomics

Rhinocerotidae emerged ~50 Mya and radiated through the Oligocene and Miocene into some of the largest land mammals ever known. Five species survive today; at least a dozen more went extinct during the Pleistocene megafaunal collapse. This module traces the family’s deep history, the Paraceratherium (Indricotherium) record, the woolly rhino Coelodonta antiquitatis, and the Liu 2021 phylogenomic resolution of the five extant species.

1. Perissodactyla & the Early Radiation

Rhinos are Perissodactyla (odd-toed ungulates) along with horses and tapirs. The clade diverged from Cetartiodactyla (even-toed ungulates + cetaceans) around 55 Mya, followed by internal splits into Equidae, Tapiridae, and Rhinocerotidae by ~50 Mya. Early rhinos were small, tapir-like browsers (e.g., Hyrachyus, ~25 kg); by the Oligocene, three distinct rhino lineages coexisted — hyracodonts, amynodonts, and rhinocerotids.

Hyracodonts produced the most spectacular lineage: Paraceratherium(also Indricotherium), ~30–20 Mya, with an estimated shoulder height of 4.8 m and a body mass of 15–20 t — the largest land mammal that ever lived (Prothero 2013). Amynodonts were hippo-like aquatic rhinos that went extinct in the Miocene. Only Rhinocerotidae survives to the present.

2. Pleistocene & the Woolly Rhino

The Pleistocene (2.6–0.01 Mya) hosted a global rhino fauna: the woolly rhino (Coelodonta antiquitatis) grazed the Eurasian mammoth steppe from Spain to Korea; Elasmotherium sibiricum — the “Siberian unicorn,” with a dome-like basal horn core — survived until at least 39 ka BP (Kosintsev 2018); the Stephanorhinus genus occupied temperate Europe. All went extinct during the Late Pleistocene megafaunal wave, with climate and human hunting both implicated (Stuart & Lister 2014).

Ancient DNA now supplements morphological phylogenetics. Lord 2020 recovered a full woolly-rhino mitogenome from cave-bear gut contents; Liu 2021 resolved five extant species’ phylogeny from nuclear genomes.

3. The Five Living Species

Species	Range	Mass	IUCN	Est. wild pop.
White (Ceratotherium simum)	S + E Africa	1800–2500 kg	NT / CR (N. white)	~16 000 / 2
Black (Diceros bicornis)	S + E Africa	800–1400 kg	CR	~6 500
Indian (Rhinoceros unicornis)	N India, Nepal	1800–2500 kg	VU	~4 000
Javan (Rhinoceros sondaicus)	Ujung Kulon, Java	900–2300 kg	CR	~75
Sumatran (Dicerorhinus sumatrensis)	Sumatra, Borneo	600–950 kg	CR	<80

Sumatran rhinos retain the body hair and two-horn configuration of the extinct woolly rhino — they are the closest living analog. Liu 2021 dated the Asian/African split at ~16 Mya and confirmed that African rhinos (white + black) are sister to the Asian clade (Indian + Javan + Sumatran), with Sumatran the earliest-branching Asian.

Simulation: Phylogeny & Body-Mass Comparison

Plots rhinoceros divergence times (Liu 2021) and a body-mass comparison spanning extant species, the Pleistocene woolly rhino and Elasmotherium, and the 17 t Paraceratherium.

Python

script.py60 lines

import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

# Liu 2021 rhinoceros phylogenomic divergence times
species = {
    'White (C. simum)':       0,
    'Black (D. bicornis)':    4,
    'Indian (R. unicornis)':  19,
    'Javan (R. sondaicus)':   21,
    'Sumatran (D. sumatrensis)': 24,
    'Woolly (Coelodonta) extinct 0.01 Mya': 25,
    'Elasmotherium extinct 0.04 Mya':         28,
}
masses = {
    'White (C. simum)':       2300,
    'Black (D. bicornis)':    1100,
    'Indian (R. unicornis)':  2200,
    'Javan (R. sondaicus)':    900,
    'Sumatran (D. sumatrensis)': 800,
    'Woolly (Coelodonta) extinct 0.01 Mya': 3000,
    'Elasmotherium extinct 0.04 Mya':     4500,
    'Paraceratherium (Indricotherium)': 17000,
}
timeline = [(-50, 'Rhinocerotidae crown'),
            (-35, 'Hyracodonts peak'),
            (-28, 'Paraceratherium'),
            (-23, 'Modern genera diverge'),
            (-2,  'Woolly rhino radiation'),
            (-0.04,'Last Elasmotherium'),
            (-0.01,'Last Coelodonta')]

fig, (ax1, ax2) = plt.subplots(1, 2, figsize=(15, 6), facecolor='#0a0a1a')
for ax in (ax1, ax2):
    ax.set_facecolor('#111827')
    ax.tick_params(colors='#cbd5e1')
    for s in ax.spines.values(): s.set_color('#334155')
    ax.grid(True, color='#334155', alpha=0.3)

names = list(species.keys())
ages  = list(species.values())
ax1.barh(names, ages, color='#d6d3d1', edgecolor='#fcd34d', height=0.6)
ax1.set_xlabel('Divergence from white rhino (Mya)', color='#cbd5e1')
ax1.set_title('Rhinocerotidae phylogeny (Liu 2021)',
              color='#fcd34d', fontweight='bold')

kvs = sorted(masses.items(), key=lambda kv: kv[1])
ax2.barh([k for k,_ in kvs], [v for _,v in kvs],
         color=['#a8a29e']*(len(kvs)-1) + ['#dc2626'], edgecolor='#fcd34d')
ax2.axvline(2300, color='#fbbf24', ls='--', label='White rhino 2.3 t')
ax2.set_xlabel('Body mass (kg)', color='#cbd5e1')
ax2.set_title('Rhino body-mass comparison',
              color='#fcd34d', fontweight='bold')
ax2.legend(facecolor='#1e293b', edgecolor='#334155', labelcolor='#cbd5e1')

plt.tight_layout()
plt.savefig('output.png', dpi=120, bbox_inches='tight', facecolor='#0a0a1a')
print('White-black split: ~4 Mya (Liu 2021 genome)')
print('Paraceratherium (Indricotherium) largest land mammal ever: 15-20 t')

Click Run to execute the Python code

Code will be executed with Python 3 on the server

4. Genomic Signatures

The five extant genomes (Liu 2021) show that the northern white (C. s. cottoni) and southern white (C. s. simum) diverged ~1 Mya and carry 0.1% nucleotide diversity (mammalian norm ~0.3%). The Sumatran rhino has the lowest nuclear diversity of any non-domestic mammal measured — under 0.05% — reflecting prolonged isolation and a genetic bottleneck. Inbreeding coefficients for Javan and Sumatran rhinos approach 0.4, producing elevated load from recessive deleterious alleles.

Genomic conservation tools (Mays 2018) include identifying purging candidates, designing assortative matings, and assessing the feasibility of cryopreserved- gamete genetic rescue — the program that underwrites the IVF recovery attempt for the northern white rhino (M8).

Key References

• Prothero, D. R. (2013). Rhinoceros Giants: The Palaeobiology of Indricotheres. Indiana University Press.

• Liu, S. et al. (2021). “Ancient and modern genomes unravel the evolutionary history of the rhinoceros family.” Cell, 184, 4874–4885.

• Kosintsev, P. et al. (2019). “Evolution and extinction of the giant rhinoceros Elasmotherium sibiricum.” Nat. Ecol. Evol., 3, 31–38.

• Stuart, A. J. & Lister, A. M. (2014). “Chronology and causation of the Late Pleistocene megafaunal extinctions.” Quat. Sci. Rev., 96, 108–116.

• Lord, E. et al. (2020). “Pre-extinction demographic stability and genomic signatures of adaptation in the woolly rhinoceros.” Curr. Biol., 30, 3871–3879.

• Mays, H. L. et al. (2018). “Genomic analysis of demographic history and ecological niche modeling in the endangered Sumatran rhinoceros.” Curr. Biol., 28, 70–76.

← Course overview Module 1: Anatomy & Scaling →