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.

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 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.

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.

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.

Compliance with the U.S. Nuclear Regulatory Commission Revised Licensing Guidance for Radioactive Seed Localization.

The use of radioactive seed localization (RSL) as a reliable alternative to wire localization (WL) for guiding the surgical excision of non-palpable lesions has become more popular because of its demonstrated advantages for both the patient and surgeon. RSL is regulated under 10 CFR 35.1000 (Other Medical Uses) for which the U.S. Nuclear Regulatory Commission (U.S. NRC) issued its original licensing guidance entitled "I-125 and Pd-103 Low Dose Rate Brachytherapy Seeds Used for Localization of Non-palpable Lesions" in 2006. At that time, RSL was performed as an off-label use of the same radioactive I seeds used for brachytherapy, so the focus of this initial guidance was to establish the same requirements for RSL as for brachytherapy. Strict compliance with the licensing guidance was burdensome and made it difficult for some licensees to implement an RSL program. In response to a request from the user community and recommendations from an Advisory Committee on the Medical Use of Isotopes (ACMUI) subcommittee, the U.S. NRC formed a working group to consider revisions to its RSL guidance to make it more relevant to the conduct of the procedure. On October 7, 2016, the U.S. NRC issued a revised licensing guidance document titled "Low Activity Radioactive Seeds Used for Localization of Non-Palpable Lesions and Lymph Nodes." This manuscript highlights the changes in the latest U.S. NRC licensing guidance for RSL, including authorized user training and experience, written directives, surveys, instrumentation, and medical event criteria, and presents practical ways to satisfy the new requirements.

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.

Evolutionary Genomics and Conservation of the Endangered Przewalski's Horse.

Przewalski's horses (PHs, Equus ferus ssp. przewalskii) were discovered in the Asian steppes in the 1870s and represent the last remaining true wild horses. PHs became extinct in the wild in the 1960s but survived in captivity, thanks to major conservation efforts. The current population is still endangered, with just 2,109 individuals, one-quarter of which are in Chinese and Mongolian reintroduction reserves [1]. These horses descend from a founding population of 12 wild-caught PHs and possibly up to four domesticated individuals [2-4]. With a stocky build, an erect mane, and stripped and short legs, they are phenotypically and behaviorally distinct from domesticated horses (DHs, Equus caballus). Here, we sequenced the complete genomes of 11 PHs, representing all founding lineages, and five historical specimens dated to 1878-1929 CE, including the Holotype. These were compared to the hitherto-most-extensive genome dataset characterized for horses, comprising 21 new genomes. We found that loci showing the most genetic differentiation with DHs were enriched in genes involved in metabolism, cardiac disorders, muscle contraction, reproduction, behavior, and signaling pathways. We also show that DH and PH populations split ∼45,000 years ago and have remained connected by gene-flow thereafter. Finally, we monitor the genomic impact of ∼110 years of captivity, revealing reduced heterozygosity, increased inbreeding, and variable introgression of domestic alleles, ranging from non-detectable to as much as 31.1%. This, together with the identification of ancestry informative markers and corrections to the International Studbook, establishes a framework for evaluating the persistence of genetic variation in future reintroduced populations.

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.