Long-term Small Population Size, Deleterious Variation, and Altitude Adaptation in the Ethiopian Wolf, a Severely Endangered Canid

Genomes of critically endangered saola are shaped by population structure and purging

The saola is one of the most elusive large mammals, standing at the brink of extinction. We constructed a reference genome and resequenced 26 saola individuals, confirming the saola as a basal member of the Bovini. Despite its small geographic range, we found that the saola is partitioned into two populations with high genetic differentiation (FST = 0.49). We estimate that these populations diverged and started declining 5,000-20,000 years ago, possibly due to climate changes and exacerbated by increasing human activities. The saola has long tracts without genomic diversity; however, most of these tracts are not shared by the two populations. Saolas carry a high genetic load, yet their gradual decline resulted in the purging of the most deleterious genetic variation. Finally, we find that combining the two populations, e.g., in an eventual captive breeding program, would mitigate the genetic load and increase the odds of species survival.

Purging of Highly Deleterious Mutations Through an Extreme Bottleneck

Transitions to captivity often produce population bottlenecks. On the one hand, bottlenecks increase inbreeding and decrease effective population size, thus increasing extinction risk. On the other hand, elevated homozygosity associated with inbreeding may purge deleterious mutations. Previous empirical studies of purging in captive breeding programs have focused on phenotypic measurements. We test natural selection's ability to purge deleterious mutations following an extreme population bottleneck by analyzing patterns of genetic diversity in wild and captive-bred individuals of the Lord Howe Island stick insect, Dryococelus australis. Dryococelus australis has been bred in captivity for two decades, having passed through an extreme bottleneck-only two mating pairs with few new additions since then. The magnitude of the bottleneck together with high female fecundity but low offspring recruitment set up nearly ideal conditions for the purging of deleterious mutations. As expected, captive-bred individuals had a greater number of long runs of homozygosity compared with wild individuals, implying strong inbreeding in captivity which would facilitate purging in homozygous regions. Stop-codon mutations were preferentially depleted in captivity compared with other mutations in coding and noncoding regions. The more deleterious a mutation was predicted to be, the more likely it was found outside of runs of homozygosity, implying that inbreeding facilitates the expression and thus removal of deleterious mutations, even after such an extreme bottleneck and under the benign conditions of captivity. These data implicate inbreeding and recessive deleterious mutation load in fitness variation among captive and wild D. australis.

A Pipeline and Recommendations for Population and Individual Diagnostic SNP Selection in Non-Model Species

Despite substantial reductions in the cost of sequencing over the last decade, genetic panels remain relevant due to their cost-effectiveness and flexibility across a variety of sample types. In particular, single nucleotide polymorphism (SNP) panels are increasingly favoured for conservation applications. SNP panels are often used because of their adaptability, effectiveness with low-quality samples, and cost-efficiency for population monitoring and forensics. However, the selection of diagnostic SNPs for population assignment and individual identification can be challenging. The consequences of poor SNP selection are under-powered panels, inaccurate results, and monetary loss. Here, we develop a novel and user-friendly SNP selection pipeline (mPCRselect) that can be used to select SNPs for population assignment and/or individual identification. mPCRselect allows any researcher, who has sufficient SNP-level data, to design a successful and cost-effective SNP panel for a diploid species of conservation concern.

Whole-Genome Evaluation of Genetic Rescue: The Case of a Curiously Isolated and Endangered Butterfly

Genetic rescue, or the translocation of individuals among populations to augment gene flow, can help ameliorate inbreeding depression and loss of adaptive potential in small and isolated populations. Genetic rescue is currently under consideration for an endangered butterfly in Canada, the Half-moon Hairstreak (Satyrium semiluna). A small, unique population persists in Waterton Lakes National Park, Alberta, isolated from other populations by more than 400 km. However, whether genetic rescue would actually be helpful has not been evaluated. Here, we generate the first chromosome-level genome assembly and whole-genome resequence data for the species. We find that the Alberta population maintains extremely low genetic diversity and is genetically very divergent from the nearest populations in British Columbia and Montana. Runs of homozygosity suggest this is due to a long history of inbreeding, and coalescent analyses show that the population has been small and isolated, yet stable, for up to 40k years. When a population like this maintains its viability despite inbreeding and low genetic diversity, it has likely undergone purging of deleterious recessive alleles and could be threatened by the reintroduction of such alleles via genetic rescue. Ecological niche modelling indicates that the Alberta population also exhibits environmental associations that are atypical of the species. Together, these evolutionary and ecological divergences suggest that population crosses may result in outbreeding depression. We therefore infer that genetic rescue has a relatively unique potential to be harmful rather than helpful for this population at present. However, because of its reduced adaptive potential, the Alberta population may still benefit from future genetic rescue as climate and habitat conditions change. Proactive experimental population crosses should therefore be completed to assess reproductive compatibility and progeny fitness.

Genomic Insights Into Red Squirrels in Scotland Reveal Loss of Heterozygosity Associated With Extreme Founder Effects

Remnant populations of endangered species often have complex demographic histories associated with human impact. This can present challenges for conservation as populations modified by human activity may require bespoke management. The Eurasian red squirrel, Sciurus vulgaris (L., 1758), is endangered in the UK. Scotland represents a key stronghold, but Scottish populations have been subjected to intense anthropogenic influence, including widespread extirpations, reintroductions and competition from an invasive species. This study examined the genetic legacy of these events through low coverage whole-genome resequencing of 106 red squirrels. Previously undetected patterns of population structure and gene flow were uncovered. One offshore island, four mainland Scottish populations, and a key east-coast migration corridor were observed. An abrupt historical population bottleneck, related to extreme founder effects, has led to a severe and prolonged depression in genome-wide heterozygosity, which is amongst the lowest reported for any species. Current designated red squirrel conservation stronghold locations do not encompass all existing diversity. These findings highlight the genetic legacies of past anthropogenic influence on long-term diversity in endangered taxa. Continuing management interventions and regular genetic monitoring are recommended to safeguard and improve future diversity.

Genomic Evidence for the Purging of Deleterious Genetic Variation in the Endangered North Atlantic Right Whale

The reduced genetic diversity and frequent inbreeding associated with small population size may underpin the accumulation and expression of deleterious mutations (mutation load) in some declining populations. However, demographic perturbations and inbreeding coupled with purifying selection can also purge declining populations of deleterious mutations, leading to intriguing recoveries. To better understand the links between deleterious genetic variation and population status, we assess patterns of genetic diversity, inbreeding, and mutation load across the genomes of three species of Balaenidae whale with different demographic histories and recoveries following the end of commercial whaling in the 1980s. Unlike bowhead (BH) and Southern right whales (SRW), which show signs of recent recovery, reproductive rates of the endangered North Atlantic right whale (NARW) remain lower than expected. We show that the NARW is currently marked by low genetic diversity, historical inbreeding, and a high mutation load. Still, we reveal evidence that genetic purging has reduced the frequency of highly deleterious alleles in NARW, which could increase chances of future population recovery. We also identify a suite of mutations putatively linked to congenital defects that occur at high frequencies in nulliparous NARW females but are rare in NARW with high reproductive success. These same mutations are nearly absent in BH and SRW in this study, suggesting that the purging of key variants may shape the probability of population recovery. As anthropogenic disturbances continue to reduce the sizes of many populations in nature, resolving the links between population dynamics and mutation load could become increasingly important.

Evolutionary consequences of domestication on the selective effects of new amino acid changing mutations in canids

The domestication of wild canids led to dogs no longer living in the wild but instead residing alongside humans. Extreme changes in behavior and diet associated with domestication may have led to the relaxation of the selective pressure on traits that may be less important in the domesticated context. Thus, here we hypothesize that strongly deleterious mutations may have become less deleterious in domesticated populations. We test this hypothesis by estimating the distribution of fitness effects (DFE) for new amino acid changing mutations using whole-genome sequence data from 24 gray wolves and 61 breed dogs. We find that the DFE is strikingly similar across canids, with 26-28% of new amino acid changing mutations being neutral/nearly neutral (|s| < 1e-5), and 41-48% under strong purifying selection (|s| > 1e-2). Our results are robust to different model assumptions suggesting that the DFE is stable across short evolutionary timescales, even in the face of putative drastic changes in the selective pressure caused by artificial selection during domestication and breed formation. On par with previous works describing DFE evolution, our data indicate that the DFE of amino acid changing mutations depends more strongly on genome structure and organismal characteristics, and less so on shifting selective pressures or environmental factors. Given the constant DFE and previous data showing that genetic variants that differentiate wolf and dog populations are enriched in regulatory elements, we speculate that domestication may have had a larger impact on regulatory variation than on amino acid changing mutations.

Genomic signatures of inbreeding and mutation load in tree ferns

Ferns (Pteridophyta), as the second largest group of vascular plants, play important roles in ecosystem functioning. Homosporous ferns exhibit a remarkable range of mating systems, from extreme inbreeding to obligate outcrossing, which may have significant evolutionary and ecological implications. Despite their significance, the impact of genome-wide inbreeding on genetic diversity and mutation load within the fern lineage remain largely unexplored. In this study, we utilized whole-genome sequencing to investigate the genomic signatures of inbreeding and genetic load in three Alsophila tree fern species. Our analysis revealed extremely high inbreeding in A. spinulosa, in contrast to the predominantly outcrossing observed in A. costularis and A. latebrosa. This difference likely reflects divergent mating systems and demographic histories. Consistent with its extreme inbreeding propensity, A. spinulosa exhibits reduced genetic diversity and a pronounced decline in effective population size. Comparison of genetic load revealed an overall reduction in deleterious mutations in the highly inbred A. spinulosa, highlighting that long-term inbreeding may have contributed to the purging of strongly deleterious mutations, thereby prolonging the survival of A. spinulosa. Despite this, however, A. spinulosa carries a substantive realized genetic load that may potentially instigate future fitness decline. Our findings illuminate the complex evolutionary interplay between inbreeding and mutation load in homosporous ferns, yielding insights with important implications for the conservation and management of these species.

Genomic adaptation to small population size and saltwater consumption in the critically endangered Cat Ba langur

Many mammal species have declining populations, but the consequences of small population size on the genomic makeup of species remain largely unknown. We investigated the evolutionary history, genetic load and adaptive potential of the Cat Ba langur (Trachypithecus poliocephalus), a primate species endemic to Vietnam's famous Ha Long Bay and with less than 100 living individuals one of the most threatened primates in the world. Using high-coverage whole genome data of four wild individuals, we revealed the Cat Ba langur as sister species to its conspecifics of the northern limestone langur clade and found no evidence for extensive secondary gene flow after their initial separation. Compared to other primates and mammals, the Cat Ba langur showed low levels of genetic diversity, long runs of homozygosity, high levels of inbreeding and an excess of deleterious mutations in homozygous state. On the other hand, genetic diversity has been maintained in protein-coding genes and on the gene-rich human chromosome 19 ortholog, suggesting that the Cat Ba langur retained most of its adaptive potential. The Cat Ba langur also exhibits several unique non-synonymous variants that are related to calcium and sodium metabolism, which may have improved adaptation to high calcium intake and saltwater consumption.

Unraveling the genomic diversity and admixture history of captive tigers in the United States

Genomic studies of endangered species have primarily focused on describing diversity patterns and resolving phylogenetic relationships, with the overarching goal of informing conservation efforts. However, few studies have investigated genomic diversity housed in captive populations. For tigers (Panthera tigris), captive individuals vastly outnumber those in the wild, but their diversity remains largely unexplored. Privately owned captive tiger populations have remained an enigma in the conservation community, with some believing that these individuals are severely inbred, while others believe they may be a source of now-extinct diversity. Here, we present a large-scale genetic study of the private (non-zoo) captive tiger population in the United States, also known as "Generic" tigers. We find that the Generic tiger population has an admixture fingerprint comprising all six extant wild tiger subspecies. Of the 138 Generic individuals sequenced for the purpose of this study, no individual had ancestry from only one subspecies. We show that the Generic tiger population has a comparable amount of genetic diversity relative to most wild subspecies, few private variants, and fewer deleterious mutations. We observe inbreeding coefficients similar to wild populations, although there are some individuals within both the Generic and wild populations that are substantially inbred. Additionally, we develop a reference panel for tigers that can be used with imputation to accurately distinguish individuals and assign ancestry with ultralow coverage (0.25×) data. By providing a cost-effective alternative to whole-genome sequencing (WGS), the reference panel provides a resource to assist in tiger conservation efforts for both ex- and in situ populations.

Conservation genomics of wolves: The global impact of RK Wayne's research

RK Wayne has arguably been the most influential geneticist of canids, famously promoting the conservation of wolves in his homeland, the United States. His influence has been felt in other countries and regions outside the contiguous United States, where he inspired others, also including former graduate students and research fellows of his, to use modern molecular techniques to examine the evolutionary biology of canids to inform the conservation and management of wolves. In this review, we focus on the implications of Wayne's work on wolves outside the United States. He envisioned a clear future for wolf conservation research, involving the study of wolves' ecological and genetic diversity, and the description of ecotypes requiring conservation. He also documented widespread hybridization among canids and introgression of DNA from domestic dogs to wolves, a process that started dozens of thousands of years ago. His work therefore calls for innovative studies, such as examining the potential fitness benefits of introgression. Inspired by his results, for example, on the purging of deleterious alleles in small populations, wolf researchers should use novel molecular tools to challenge other conservation genetics paradigms. Overall, RK Wayne's work constitutes a call for answers, which as scientists or citizens concerned with conservation matters, we are obliged to address, as we contribute to monitoring and maintaining biodiversity during our period of dramatic transformations of the biosphere.

Genomic Reconstruction of the Successful Establishment of a Feralized Bovine Population on the Subantarctic Island of Amsterdam

The feral cattle of the subantarctic island of Amsterdam provide an outstanding case study of a large mammalian population that was established by a handful of founders and thrived within a few generations in a seemingly inhospitable environment. Here, we investigated the genetic history and composition of this population using genotyping and sequencing data. Our inference showed an intense but brief founding bottleneck around the late 19th century and revealed contributions from European taurine and Indian Ocean Zebu in the founder ancestry. Comparative analysis of whole-genome sequences further revealed a moderate reduction in genetic diversity despite high levels of inbreeding. The brief and intense bottleneck was associated with high levels of drift, a flattening of the site frequency spectrum and a slight relaxation of purifying selection on mildly deleterious variants. Unlike some populations that have experienced prolonged reductions in effective population size, we did not observe any significant purging of highly deleterious variants. Interestingly, the population's success in the harsh environment can be attributed to preadaptation from their European taurine ancestry, suggesting no strong bioclimatic challenge, and also contradicting evidence for insular dwarfism. Genome scan for footprints of selection uncovered a majority of candidate genes related to nervous system function, likely reflecting rapid feralization driven by behavioral changes and complex social restructuring. The Amsterdam Island cattle offers valuable insights into rapid population establishment, feralization, and genetic adaptation in challenging environments. It also sheds light on the unique genetic legacies of feral populations, raising ethical questions according to conservation efforts.

Chromosome-level assembly of the gray fox (Urocyon cinereoargenteus) confirms the basal loss of PRDM9 in Canidae

Reference genome assemblies have been created from multiple lineages within the Canidae family; however, despite its phylogenetic relevance as a basal genus within the clade, there is currently no reference genome for the gray fox (Urocyon cinereoargenteus). Here, we present a chromosome-level assembly for the gray fox (U. cinereoargenteus), which represents the most contiguous, non-domestic canid reference genome available to date, with 90% of the genome contained in just 34 scaffolds and a contig N50 and scaffold N50 of 59.4 and 72.9 Megabases, respectively. Repeat analyses identified an increased number of simple repeats relative to other canids. Based on mitochondrial DNA, our Vermont sample clusters with other gray fox samples from the northeastern United States and contains slightly lower levels of heterozygosity than gray foxes on the west coast of California. This new assembly lays the groundwork for future studies to describe past and present population dynamics, including the delineation of evolutionarily significant units of management relevance. Importantly, the phylogenetic position of Urocyon allows us to verify the loss of PRDM9 functionality in the basal canid lineage, confirming that pseudogenization occurred at least 10 million years ago.

Genomic evidence for evolutionary history and local adaptation of two endemic apricots: Prunus hongpingensis and P. zhengheensis

Apricot, belonging to the Armeniaca section of Rosaceae, is one of the economically important crop fruits that has been extensively cultivated. The natural wild apricots offer valuable genetic resources for crop improvement. However, some of them are endemic, with small populations, and are even at risk of extinction. In this study we unveil chromosome-level genome assemblies for two southern China endemic apricots, Prunus hongpingensis (PHP) and P. zhengheensis (PZH). We also characterize their evolutionary history and the genomic basis of their local adaptation using whole-genome resequencing data. Our findings reveal that PHP and PZH are closely related to Prunus armeniaca and form a distinct lineage. Both species experienced a decline in effective population size following the Last Glacial Maximum (LGM), which likely contributed to their current small population sizes. Despite the observed decrease in genetic diversity and heterozygosity, we do not observe an increased accumulation of deleterious mutations in these two endemic apricots. This is likely due to the combined effects of a low inbreeding coefficient and strong purifying selection. Furthermore, we identify a set of genes that have undergone positive selection and are associated with local environmental adaptation in PHP and PZH, respectively. These candidate genes can serve as valuable genetic resources for targeted breeding and improvement of cultivated apricots. Overall, our study not only enriches our comprehension of the evolutionary history of apricot species but also offers crucial insights for the conservation and future breeding of other endemic species amidst rapid climate changes.

Population genomics of the muskox' resilience in the near absence of genetic variation

Genomic studies of species threatened by extinction are providing crucial information about evolutionary mechanisms and genetic consequences of population declines and bottlenecks. However, to understand how species avoid the extinction vortex, insights can be drawn by studying species that thrive despite past declines. Here, we studied the population genomics of the muskox (Ovibos moschatus), an Ice Age relict that was at the brink of extinction for thousands of years at the end of the Pleistocene yet appears to be thriving today. We analysed 108 whole genomes, including present-day individuals representing the current native range of both muskox subspecies, the white-faced and the barren-ground muskox (O. moschatus wardi and O. moschatus moschatus) and a ~21,000-year-old ancient individual from Siberia. We found that the muskox' demographic history was profoundly shaped by past climate changes and post-glacial re-colonizations. In particular, the white-faced muskox has the lowest genome-wide heterozygosity recorded in an ungulate. Yet, there is no evidence of inbreeding depression in native muskox populations. We hypothesize that this can be explained by the effect of long-term gradual population declines that allowed for purging of strongly deleterious mutations. This study provides insights into how species with a history of population bottlenecks, small population sizes and low genetic diversity survive against all odds.

Evolution of the Mutation Spectrum Across a Mammalian Phylogeny

Although evolutionary biologists have long theorized that variation in DNA repair efficacy might explain some of the diversity of lifespan and cancer incidence across species, we have little data on the variability of normal germline mutagenesis outside of humans. Here, we shed light on the spectrum and etiology of mutagenesis across mammals by quantifying mutational sequence context biases using polymorphism data from thirteen species of mice, apes, bears, wolves, and cetaceans. After normalizing the mutation spectrum for reference genome accessibility and k-mer content, we use the Mantel test to deduce that mutation spectrum divergence is highly correlated with genetic divergence between species, whereas life history traits like reproductive age are weaker predictors of mutation spectrum divergence. Potential bioinformatic confounders are only weakly related to a small set of mutation spectrum features. We find that clock-like mutational signatures previously inferred from human cancers cannot explain the phylogenetic signal exhibited by the mammalian mutation spectrum, despite the ability of these signatures to fit each species' 3-mer spectrum with high cosine similarity. In contrast, parental aging signatures inferred from human de novo mutation data appear to explain much of the 1-mer spectrum's phylogenetic signal in combination with a novel mutational signature. We posit that future models purporting to explain the etiology of mammalian mutagenesis need to capture the fact that more closely related species have more similar mutation spectra; a model that fits each marginal spectrum with high cosine similarity is not guaranteed to capture this hierarchy of mutation spectrum variation among species.

"Evolution of the mutation spectrum across a mammalian phylogeny"

Little is known about how the spectrum and etiology of germline mutagenesis might vary among mammalian species. To shed light on this mystery, we quantify variation in mutational sequence context biases using polymorphism data from thirteen species of mice, apes, bears, wolves, and cetaceans. After normalizing the mutation spectrum for reference genome accessibility and k -mer content, we use the Mantel test to deduce that mutation spectrum divergence is highly correlated with genetic divergence between species, whereas life history traits like reproductive age are weaker predictors of mutation spectrum divergence. Potential bioinformatic confounders are only weakly related to a small set of mutation spectrum features. We find that clocklike mutational signatures previously inferred from human cancers cannot explain the phylogenetic signal exhibited by the mammalian mutation spectrum, despite the ability of these clocklike signatures to fit each species' 3-mer spectrum with high cosine similarity. In contrast, parental aging signatures inferred from human de novo mutation data appear to explain much of the mutation spectrum's phylogenetic signal when fit to non-context-dependent mutation spectrum data in combination with a novel mutational signature. We posit that future models purporting to explain the etiology of mammalian mutagenesis need to capture the fact that more closely related species have more similar mutation spectra; a model that fits each marginal spectrum with high cosine similarity is not guaranteed to capture this hierarchy of mutation spectrum variation among species.

The earliest Ethiopian wolf: implications for the species evolution and its future survival

In 2017, a hemimandible (MW5-B208), corresponding to the Ethiopian wolf (Canis simensis), was found in a stratigraphically-controlled and radio-isotopically-dated sequence of the Melka Wakena paleoanthropological site-complex, on the Southeastern Ethiopian Highlands, ~ 2300 m above sea level. The specimen is the first and unique Pleistocene fossil of this species. Our data provide an unambiguous minimum age of 1.6-1.4 Ma for the species' presence in Africa and constitutes the first empirical evidence that supports molecular interpretations. Currently, C. simensis is one of the most endangered carnivore species of Africa. Bioclimate niche modeling applied to the time frame indicated by the fossil suggests that the lineage of the Ethiopian wolf faced severe survival challenges in the past, with consecutive drastic geographic range contractions during warmer periods. These models help to describe future scenarios for the survival of the species. Projections ranging from most pessimistic to most optimistic future climatic scenarios indicate significant reduction of the already-deteriorating territories suitable for the Ethiopian Wolf, increasing the threat to the specie's future survival. Additionally, the recovery of the Melka Wakena fossil underscores the importance of work outside the East African Rift System in research of early human origins and associated biodiversity on the African continent.