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.

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.

Response to Bakker et al.

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

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.

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.

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.

Highly accurate long-read HiFi sequencing data for five complex genomes.

The PacBio® HiFi sequencing method yields highly accurate long-read sequencing datasets with read lengths averaging 10-25 kb and accuracies greater than 99.5%. These accurate long reads can be used to improve results for complex applications such as single nucleotide and structural variant detection, genome assembly, assembly of difficult polyploid or highly repetitive genomes, and assembly of metagenomes. Currently, there is a need for sample data sets to both evaluate the benefits of these long accurate reads as well as for development of bioinformatic tools including genome assemblers, variant callers, and haplotyping algorithms. We present deep coverage HiFi datasets for five complex samples including the two inbred model genomes Mus musculus and Zea mays, as well as two complex genomes, octoploid Fragaria × ananassa and the diploid anuran Rana muscosa. Additionally, we release sequence data from a mock metagenome community. The datasets reported here can be used without restriction to develop new algorithms and explore complex genome structure and evolution. Data were generated on the PacBio Sequel II System.

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.

A High-Quality, Long-Read De Novo Genome Assembly to Aid Conservation of Hawaii's Last Remaining Crow Species.

Abstract: Genome-level data can provide researchers with unprecedented precision to examine the causes and genetic consequences of population declines, which can inform conservation management. Here, we present a high-quality, long-read, de novo genome assembly for one of the world's most endangered bird species, the 'Alala (Corvus hawaiiensis; Hawaiian crow). As the only remaining native crow species in Hawai'i, the 'Alala survived solely in a captive-breeding program from 2002 until 2016, at which point a long-term reintroduction program was initiated. The high-quality genome assembly was generated to lay the foundation for both comparative genomics studies and the development of population-level genomic tools that will aid conservation and recovery efforts. We illustrate how the quality of this assembly places it amongst the very best avian genomes assembled to date, comparable to intensively studied model systems. We describe the genome architecture in terms of repetitive elements and runs of homozygosity, and we show that compared with more outbred species, the 'Alala genome is substantially more homozygous. We also provide annotations for a subset of immunity genes that are likely to be important in conservation management, and we discuss how this genome is currently being used as a roadmap for downstream conservation applications.

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.

Identification of California Condor Estrogen Receptors 1 and 2 and Their Activation by Endocrine Disrupting Chemicals.

Recently, California condors (Gymnogyps californianus) have been reintroduced to coastal regions of California where they feed on marine mammal carcasses. There is evidence that coastal-dwelling condors experience reproductive issues, such as eggshell thinning, likely resulting from exposure to endocrine-disrupting chemicals (EDCs). To address this problem, we have identified and cloned condor estrogen receptors (ESRs) 1 and 2 and characterized their activation by EDCs present in the coastal habitats where condors reside. Dichlorodiphenyltrichloroethane (DDT) and its metabolites all activated ESR1 and ESR2, although their relative potency differed between the receptors. Bisphenol A, dieldrin, trans-nonachlor, and polychlorinated biphenyl 52 (PCB52) moderately activated both ESRs, whereas PCB138 and PCB153 stimulated little to no activation. Overall, EDC activation of condor ESR2, which is the first ESR2 cloned from a raptor species, was greater than that of ESR1. Significant activation of both condor ESRs by EDCs occurred at high concentrations (≥1µM), which are within the range of plasma levels of certain EDCs (eg, dichlorodiphenyldichloroethylene [p'p-DDE]) in coastal-dwelling condors. Finally, phylogenetic analyses of ESRs of 41 avian species identified a single amino acid position in ESR2 under positive selection. Mutation of this amino acid affected receptor activation by EDCs, suggesting the identity of this amino acid may influence EDC sensitivity of avian species. Together, these findings broaden our understanding of EDC interactions with ESRs in avian species. For condors specifically, these data could be used to evaluate EDC exposure risk at future release sites to identify those least likely to compromise the continued recovery of this species.

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.

Characterization of Prdm9 in equids and sterility in mules.

Prdm9 (Meisetz) is the first speciation gene discovered in vertebrates conferring reproductive isolation. This locus encodes a meiosis-specific histone H3 methyltransferase that specifies meiotic recombination hotspots during gametogenesis. Allelic differences in Prdm9, characterized for a variable number of zinc finger (ZF) domains, have been associated with hybrid sterility in male house mice via spermatogenic failure at the pachytene stage. The mule, a classic example of hybrid sterility in mammals also exhibits a similar spermatogenesis breakdown, making Prdm9 an interesting candidate to evaluate in equine hybrids. In this study, we characterized the Prdm9 gene in all species of equids by analyzing sequence variation of the ZF domains and estimating positive selection. We also evaluated the role of Prdm9 in hybrid sterility by assessing allelic differences of ZF domains in equine hybrids. We found remarkable variation in the sequence and number of ZF domains among equid species, ranging from five domains in the Tibetan kiang and Asiatic wild ass, to 14 in the Grevy's zebra. Positive selection was detected in all species at amino acid sites known to be associated with DNA-binding specificity of ZF domains in mice and humans. Equine hybrids, in particular a quartet pedigree composed of a fertile mule showed a mosaic of sequences and number of ZF domains suggesting that Prdm9 variation does not seem by itself to contribute to equine hybrid sterility.

Conservation genomics of threatened animal species.

The genomics era has opened up exciting possibilities in the field of conservation biology by enabling genomic analyses of threatened species that previously were limited to model organisms. Next-generation sequencing (NGS) and the collection of genome-wide data allow for more robust studies of the demographic history of populations and adaptive variation associated with fitness and local adaptation. Genomic analyses can also advance management efforts for threatened wild and captive populations by identifying loci contributing to inbreeding depression and disease susceptibility, and predicting fitness consequences of introgression. However, the development of genomic tools in wild species still carries multiple challenges, particularly those associated with computational and sampling constraints. This review provides an overview of the most significant applications of NGS and the implications and limitations of genomic studies in conservation.

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.

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.

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.