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

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

Evolution of the ancestral mammalian karyotype and syntenic regions.

Decrypting the rearrangements that drive mammalian chromosome evolution is critical to understanding the molecular bases of speciation, adaptation, and disease susceptibility. Using 8 scaffolded and 26 chromosome-scale genome assemblies representing 23/26 mammal orders, we computationally reconstructed ancestral karyotypes and syntenic relationships at 16 nodes along the mammalian phylogeny. Three different reference genomes (human, sloth, and cattle) representing phylogenetically distinct mammalian superorders were used to assess reference bias in the reconstructed ancestral karyotypes and to expand the number of clades with reconstructed genomes. The mammalian ancestor likely had 19 pairs of autosomes, with nine of the smallest chromosomes shared with the common ancestor of all amniotes (three still conserved in extant mammals), demonstrating a striking conservation of synteny for ∼320 My of vertebrate evolution. The numbers and types of chromosome rearrangements were classified for transitions between the ancestral mammalian karyotype, descendent ancestors, and extant species. For example, 94 inversions, 16 fissions, and 14 fusions that occurred over 53 My differentiated the therian from the descendent eutherian ancestor. The highest breakpoint rate was observed between the mammalian and therian ancestors (3.9 breakpoints/My). Reconstructed mammalian ancestor chromosomes were found to have distinct evolutionary histories reflected in their rates and types of rearrangements. The distributions of genes, repetitive elements, topologically associating domains, and actively transcribed regions in multispecies homologous synteny blocks and evolutionary breakpoint regions indicate that purifying selection acted over millions of years of vertebrate evolution to maintain syntenic relationships of developmentally important genes and regulatory landscapes of gene-dense chromosomes.

A New World Monkey Resembles Human in Bitter Taste Receptor Evolution and Function via a Single Parallel Amino Acid Substitution.

Bitter taste receptors serve as a vital component in the defense system against toxin intake by animals, and the family of genes encoding these receptors has been demonstrated, usually by family size variance, to correlate with dietary preference. However, few systematic studies of specific Tas2R to unveil their functional evolution have been conducted. Here, we surveyed Tas2R16 across all major clades of primates and reported a rare case of a convergent change to increase sensitivity to ß-glucopyranosides in human and a New World monkey, the white-faced saki. Combining analyses at multiple levels, we demonstrate that a parallel amino acid substitution (K172N) shared by these two species is responsible for this functional convergence of Tas2R16. Considering the specialized feeding preference of the white-faced saki, the K172N change likely played an important adaptive role in its early evolution to avoid potentially toxic cyanogenic glycosides, as suggested for the human TAS2R16 gene.

Evaluating recovery potential of the northern white rhinoceros from cryopreserved somatic cells.

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.

Facultative Parthenogenesis in California Condors.

Parthenogenesis is a relatively rare event in birds, documented in unfertilized eggs from columbid, galliform, and passerine females with no access to males. In the critically endangered California condor, parentage analysis conducted utilizing polymorphic microsatellite loci has identified two instances of parthenogenetic development from the eggs of two females in the captive breeding program, each continuously housed with a reproductively capable male with whom they had produced offspring. Paternal genetic contribution to the two chicks was excluded. Both parthenotes possessed the expected male ZZ sex chromosomes and were homozygous for all evaluated markers inherited from their dams. These findings represent the first molecular marker-based identification of facultative parthenogenesis in an avian species, notably of females in regular contact with fertile males, and add to the phylogenetic breadth of vertebrate taxa documented to have reproduced via asexual reproduction.

Identification of California condor (Gymnogyps californianus) estrogen receptor variants and their activation by xenoestrogens.

California condors released in costal sites are exposed to high levels of xenoestrogens, particularly p,p'-DDE, through scavenging of marine mammal carcasses. As a result, coastal condors carry a higher contaminant loads and experience eggshell thinning when compared to their inland counterparts. Given that condor estrogen receptors (Esrs) are activated by physiologically relevant levels of xenoestrogens, differences in vulnerability to endocrine disruption may exist depending on which Esr variant(s) an individual condor possesses. This work aims to characterize genetic polymorphisms in estrogen receptor genes (ESRs) in California condors; one identified for condor estrogen receptor 1 (ESR1) (N161S, E162D) and one in the ESR2 (T114M) gene. Each variant was confirmed in individual founder birds by direct PCR sequencing as well as in first generation offspring to understand the introduction of the alleles into the pedigree (6 birds for ESR1 and 5 birds for ESR2). Site-directed mutagenesis was performed on wild type receptors to produce each of the full-length ESR variants and activation of Esr1 and Esr2 variant and wild type receptors by xenoestrogens was compared. Maximal activation of the variant form of Esr1 was significantly higher (p < 0.05) in response to ethinyl estradiol (EE2), o,p'-DDE, p,p'-DDE, p,p'-DDT and p,p'-DDD compared to wild type Esr1. For Esr2 the wild type maximal activation was higher in response to o,p'-DDE, p,p'-DDE, o,p'-DDT, and p,p'-DDT. Although significant differences in activation of condor Esr variants by xenoestrogens occurred at high (micromolar) concentrations, they correspond to circulating concentrations previously reported in coastal birds. Release and relocation of California condors to the coast is a promising avenue for recovery, however, reproductive problems associated with xenoestrogen exposure pose a sub-lethal threat to long-term success. Based on above findings, future release decisions could be informed by ESR form(s) individual birds possess to reduce deleterious effects of xenoestrogen exposure and ultimately improve reproductive success in wild populations.

Chromosomal variation and perinatal mortality in San Diego zoo Soemmerring's gazelles.

Chromosomal translocations play a fundamental role in the evolution and speciation of antelopes (Antilopinae, Bovidae), with several species exhibiting polymorphism for centric fusions. For the past 35 years, the San Diego Zoo Global (SDZG) captive population of Soemmerring's gazelles has revealed complex karyotypes resulting from chromosomal translocations with diploid numbers ranging from 34 to 39. Poor reproductive performance of this species in captivity and elevated mortality the first month of life (perinatal) has been attributed to this chromosomal dynamism. We have extended the studies of karyotypic variation in the SDZG Soemmerring's gazelle population and analyzed the effect of chromosomal and genetic variation upon perinatal mortality. Karyotypes from 149 captive Soemmerring's gazelles were evaluated revealing two unreported autosomal combinations, now constituting a total of 15 distinct karyotypes for the 3 Robertsonian centric fusions originally described for this population. Among SDZG founders, distinct chromosomal variation and nuclear and mitochondrial genetic structure were detected corresponding to the institution of origin of the founders. Low levels of genetic distance and nucleotide diversity among individuals, in addition to high relatedness values, suggested that outbreeding is less of a concern than inbreeding for maintaining a sustainable captive population. Finally, analysis of karyotypes of offspring born into the SDZG Soemmerring's gazelle herds, in conjunction with the maternal karyotype showed association of chromosomal makeup with perinatal mortality. This supports the importance of continuing cytogenetic screening efforts, particularly to evaluate the presence of deleterious chromosomal rearrangements in stillborns.

Molecular phylogeny and chromosomal evolution of Alcelaphini (Antilopinae).

Robertsonian (Rb) translocations, in particular centric fusions, are thought to play a primary role in evolution and speciation of the Bovidae family. However, Rb fusions are often polymorphic within species, being suggested as phylogenetically uninformative characters. This work studies chromosome variation in 72 captive individuals of 6 species of Alcelaphini (Antilopinae): The hartebeest (genus Alcelaphus), hirola (Beatragus), black and blue wildebeests (Connochaetes), and the topi and bontebok (Damaliscus). We infer the phylogenic relationships among Alcelaphini species and determine patterns of chromosomal evolution using G-banded karyotypes and complete mitochondrial genome sequences. The molecular phylogeny showed an early divergence of Connochaetes, followed by the split of Alcelaphus plus Beatragus + Damaliscus as sister taxa. Mitochondrial and chromosomal phylogenies only differed in the position of the critically endangered Beatragus, likely due to homoplasic chromosome characters. Patterns of chromosome evolution, reconstructed using a probabilistic approach, suggest that chromosome changes leading to speciation in Alcelaphini do not exclusively involve consecutive reduction of diploid number through centric fusion but also the losses and reversions of Rb translocations in Beatragus and Damaliscus lineages. Our results provide evidence that complex scenarios of chromosomal rearrangements can be detected in relatively recent-diverged bovids, as in this group of antelopes.

Genetic load and viability of a future restored northern white rhino population.

As biodiversity loss outpaces recovery, conservationists are increasingly turning to novel tools for preventing extinction, including cloning and in vitro gametogenesis of biobanked cells. However, restoration of populations can be hindered by low genetic diversity and deleterious genetic load. The persistence of the northern white rhino (Ceratotherium simum cottoni) now depends on the cryopreserved cells of 12 individuals. These banked genomes have higher genetic diversity than southern white rhinos (C. s. simum), a sister subspecies that successfully recovered from a severe bottleneck, but the potential impact of genetic load is unknown. We estimated how demographic history has shaped genome-wide genetic load in nine northern and 13 southern white rhinos. The bottleneck left southern white rhinos with more fixed and homozygous deleterious alleles and longer runs of homozygosity, whereas northern white rhinos retained more deleterious alleles masked in heterozygosity. To gauge the impact of genetic load on the fitness of a northern white rhino population restored from biobanked cells, we simulated recovery using fitness of southern white rhinos as a benchmark for a viable population. Unlike traditional restoration, cell-derived founders can be reintroduced in subsequent generations to boost lost genetic diversity and relieve inbreeding. In simulations with repeated reintroduction of founders into a restored population, the fitness cost of genetic load remained lower than that borne by southern white rhinos. Without reintroductions, rapid growth of the restored population (>20-30% per generation) would be needed to maintain comparable fitness. Our results suggest that inbreeding depression from genetic load is not necessarily a barrier to recovery of the northern white rhino and demonstrate how restoration from biobanked cells relieves some constraints of conventional restoration from a limited founder pool. Established conservation methods that protect healthy populations will remain paramount, but emerging technologies hold promise to bolster these tools to combat the extinction crisis.

The contribution of historical processes to contemporary extinction risk in placental mammals.

Species persistence can be influenced by the amount, type, and distribution of diversity across the genome, suggesting a potential relationship between historical demography and resilience. In this study, we surveyed genetic variation across single genomes of 240 mammals that compose the Zoonomia alignment to evaluate how historical effective population size (Ne) affects heterozygosity and deleterious genetic load and how these factors may contribute to extinction risk. We find that species with smaller historical Ne carry a proportionally larger burden of deleterious alleles owing to long-term accumulation and fixation of genetic load and have a higher risk of extinction. This suggests that historical demography can inform contemporary resilience. Models that included genomic data were predictive of species' conservation status, suggesting that, in the absence of adequate census or ecological data, genomic information may provide an initial risk assessment.

Evolutionary constraint and innovation across hundreds of placental mammals.

Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.

Molecular phylogeny of extant equids and effects of ancestral polymorphism in resolving species-level phylogenies.

Short divergence times and processes such as incomplete lineage sorting and species hybridization are known to hinder the inference of species-level phylogenies due to the lack of sufficient informative genetic variation or the presence of shared but incongruent polymorphism among taxa. Extant equids (horses, zebras, and asses) are an example of a recently evolved group of mammals with an unresolved phylogeny, despite a large number of molecular studies. Previous surveys have proposed trees with rather poorly supported nodes, and the bias caused by genetic introgression or ancestral polymorphism has not been assessed. Here we studied the phylogenetic relationships of all extant species of Equidae by analyzing 22 partial mitochondrial and nuclear genes using maximum likelihood and Bayesian inferences that account for heterogeneous gene histories. We also examined genetic signatures of lineage sorting and/or genetic introgression in zebras by evaluating patterns of intraspecific genetic variation. Our study improved the resolution and support of the Equus phylogeny and in particular the controversial positions of the African wild ass (E. asinus) and mountain zebra (E. zebra): the African wild ass is placed as a sister species of the Asiatic asses and the mountain zebra as the sister taxon of Grevy's and Burchell's zebras. A shared polymorphism (indel) detected among zebra species in the Estrogen receptor 1 gene was likely due to incomplete lineage sorting and not genetic introgression as also indicated by other mitochondrial (Cytochrome b) and nuclear (Y chromosome and microsatellites) markers. Ancestral polymorphism in equids might have contributed to the long-standing lack of clarity in the phylogeny of this highly threatened group of mammals.

Response to Bakker et al.

Robinson and colleagues respond to the points raised about their paper by Bakker et al.

A Chromosome-Length Reference Genome for the Endangered Pacific Pocket Mouse Reveals Recent Inbreeding in a Historically Large Population.

High-quality reference genomes are fundamental tools for understanding population history, and can provide estimates of genetic and demographic parameters relevant to the conservation of biodiversity. The federally endangered Pacific pocket mouse (PPM), which persists in three small, isolated populations in southern California, is a promising model for studying how demographic history shapes genetic diversity, and how diversity in turn may influence extinction risk. To facilitate these studies in PPM, we combined PacBio HiFi long reads with Omni-C and Hi-C data to generate a de novo genome assembly, and annotated the genome using RNAseq. The assembly comprised 28 chromosome-length scaffolds (N50 = 72.6 MB) and the complete mitochondrial genome, and included a long heterochromatic region on chromosome 18 not represented in the previously available short-read assembly. Heterozygosity was highly variable across the genome of the reference individual, with 18% of windows falling in runs of homozygosity (ROH) >1 MB, and nearly 9% in tracts spanning >5 MB. Yet outside of ROH, heterozygosity was relatively high (0.0027), and historical Ne estimates were large. These patterns of genetic variation suggest recent inbreeding in a formerly large population. Currently the most contiguous assembly for a heteromyid rodent, this reference genome provides insight into the past and recent demographic history of the population, and will be a critical tool for management and future studies of outbreeding depression, inbreeding depression, and genetic load.

Species identification and chromosome variation of captive two-toed sloths.

Two-toed sloth species, Linnaeus's and Hoffmman's, are frequent residents of zoo collections in North America. However, species identification has always been problematic because of their large overlap in external morphology, which represents an obstacle to the captive breeding program. We describe here a PCR-based technique that allows species identification of two-toed sloths without requiring sequencing, by using a mitochondrial marker (COI gene) and restriction enzyme assay. We also report intra- and inter-specific patterns of chromosome variation in captive two-toed sloths. Molecularly, we identified 22 samples of Linnaeus's and Hoffmman's two-toed sloths corresponding to 14 and 8 individuals, respectively. One animal was identified as a hybrid using the nuclear gene Enam having alleles derived from both species. The chromosome number in Hoffman's two-toed sloths showed low variation ranging only between 50 and 51. In contrast, Linnaeus's two-toed sloths appeared to vary widely, with diploid numbers ranging from 53 to 67, suggesting distinct geographic groups. The species identification method presented here represents a low-cost easy-to-use tool that will help to improve management of the captive population of two-toed sloths.

Genome-wide diversity in the California condor tracks its prehistoric abundance and decline.

Due to their small population sizes, threatened and endangered species frequently suffer from a lack of genetic diversity, potentially leading to inbreeding depression and reduced adaptability.1 During the latter half of the twentieth century, North America's largest soaring bird,2 the California condor (Gymnogyps californianus; Critically Endangered3), briefly went extinct in the wild. Though condors once ranged throughout North America, by 1982 only 22 individuals remained. Following decades of captive breeding and release efforts, there are now >300 free-flying wild condors and ∼200 in captivity. The condor's recent near-extinction from lead poisoning, poaching, and loss of habitat is well documented,4 but much about its history remains obscure. To fill this gap and aid future management of the species, we produced a high-quality chromosome-length genome assembly for the California condor and analyzed its genome-wide diversity. For comparison, we also examined the genomes of two close relatives: the Andean condor (Vultur gryphus; Vulnerable3) and the turkey vulture (Cathartes aura; Least Concern3). The genomes of all three species show evidence of historic population declines. Interestingly, the California condor genome retains a high degree of variation, which our analyses reveal is a legacy of its historically high abundance. Correlations between genome-wide diversity and recombination rate further suggest a history of purifying selection against linked deleterious alleles, boding well for future restoration. We show how both long-term evolutionary forces and recent inbreeding have shaped the genome of the California condor, and provide crucial genomic resources to enable future research and conservation.

Genomic insights into the conservation status of the world's last remaining Sumatran rhinoceros populations.

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.

The genetic basis of phenotypic convergence in beach mice: similar pigment patterns but different genes.

Convergent evolution is a widespread phenomenon seen in diverse organisms inhabiting similar selective environments. However, it is unclear if similar phenotypes are produced by the same or different genes and mutations. Here we analyze the molecular mechanisms underlying convergent pigment pattern among subspecies of the beach mouse (Peromyscus polionotus) inhabiting the Gulf and Atlantic coasts of Florida. In these two geographic regions, separated by more than 300 km, "beach mice" have lighter colored coats than do their mainland counterparts, produced by natural selection for camouflage against the pale coastal sand dunes. We measured color pattern in eight beach mouse subspecies and showed that three of the Gulf Coast subspecies are more phenotypically similar to an Atlantic coast subspecies than to their Gulf Coast neighbors. However, light-colored beach mice do not form a monophyletic group. Previous results implicated a single derived amino acid change in the melanocortin-1 receptor (Mc1r) as a major contributor to pigment pattern in the Gulf Coast beach mice; despite phenotypic similarities, the derived Mc1r allele was not found in the Atlantic coast beach mouse populations. Here we show that Atlantic coast beach mice have high levels of Mc1r polymorphism but they lack unique alleles. Functional assays revealed that single amino acid mutations segregating in Atlantic coast beach mice do not cause any change in Mc1r activity compared with the activity of Mc1r from dark-colored mice. These joint results show that convergent pigment patterns in recently diverged beach mouse subspecies--whose developmental constraints are presumably similar--have evolved through a diversity of genetic mechanisms.

Adaptive variation in beach mice produced by two interacting pigmentation genes.

Little is known about the genetic basis of ecologically important morphological variation such as the diverse color patterns of mammals. Here we identify genetic changes contributing to an adaptive difference in color pattern between two subspecies of oldfield mice (Peromyscus polionotus). One mainland subspecies has a cryptic dark brown dorsal coat, while a younger beach-dwelling subspecies has a lighter coat produced by natural selection for camouflage on pale coastal sand dunes. Using genome-wide linkage mapping, we identified three chromosomal regions (two of major and one of minor effect) associated with differences in pigmentation traits. Two candidate genes, the melanocortin-1 receptor (Mc1r) and its antagonist, the Agouti signaling protein (Agouti), map to independent regions that together are responsible for most of the difference in pigmentation between subspecies. A derived mutation in the coding region of Mc1r, rather than change in its expression level, contributes to light pigmentation. Conversely, beach mice have a derived increase in Agouti mRNA expression but no changes in protein sequence. These two genes also interact epistatically: the phenotypic effects of Mc1r are visible only in genetic backgrounds containing the derived Agouti allele. These results demonstrate that cryptic coloration can be based largely on a few interacting genes of major effect.