Genomic signatures of extensive inbreeding in Isle Royale wolves, a population on the threshold of extinction

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.

Museomics of an extinct European flat oyster population

Understanding the factors that predispose species and populations to decline and extinction is a major challenge of biodiversity research. In the present study, we investigated the historical population genomics of an extinct population of the European oyster (Ostrea edulis L.) from the Wadden Sea collected between 1868 and 1888, and compared it to French and English populations sampled at the same time. Our museomic results indicate that the now-extinct population was genetically isolated from the French and English populations and showed signs of local adaptation in the form of Fst outlier loci between the Wadden Sea and the other two populations. Thus the Wadden Sea oysters may have been predisposed for extinction because they were not naturally replenished from other populations. A comparison of population-wide genomic diversity may hint towards a sudden population contraction of the Wadden Sea population, possibly being the result of stronger - or earlier - population decline in this population than in the others. In summary, our historical population genomic exploration hints at some potential causes of population decline in flat oysters from the Wadden Sea, which might have led to their extinction.

Constraints to gene flow increase the risk of genome erosion in the Ngorongoro Crater lion population

Small, isolated populations are at greater risk of genome erosion than larger populations. Successful conservation efforts may lead to demographic recovery and mitigate the negative genetic effects of bottlenecks. However, constrained gene flow can hamper genomic recovery. Here, we use population genomic analyses and forward simulations to assess the genomic impacts of near extinction in the isolated Ngorongoro Crater lion (Panthera leo) sub-population. We show that 200 years of quasi-isolation and the recent epizootic in 1962 resulted in a two-fold increase in inbreeding and an excess in the frequency of highly deleterious mutations relative to other populations of the Greater Serengeti. There was little evidence for purging of genetic load. Furthermore, forward simulations indicate that higher gene flow from outside of the Crater is needed to prevent future genomic erosion in the population, with a minimum of one to five effective male migrants per decade required to reduce the risk of long-term inbreeding depression and reduction in genetic diversity. Our results suggest that in spite of a rapid post-epizootic demographic recovery since the 1970s, continued isolation of the population driven by habitat fragmentation and potentially male territoriality, exacerbate the effects of genome erosion.

Allele co-segregation and haplotype diversity of MHC IIβ genes in the small-spotted catshark Scyliorhinus canicula

The major histocompatibility complex (MHC) constitutes a functionally relevant multigene family playing an essential role in the adaptive immune responses of jawed vertebrates, being directly involved in pathogen recognition. MHC diversity, driven by pathogen-mediated selection, is vital for species survival and is characterized by high genetic diversity in many taxa, namely at the sequence, allelic and haplotype levels. Chondrichthyans, the most basal jawed vertebrates with an adaptive immune system, exhibit a high diversity of MHC gene lineages conservatively organized in a compact region of the genome. Such genomic architecture suggests linkage among MHC genes, where alleles from different genes possibly co-segregate together. Such condition may have major implications on immune response, individual fitness and evolution. In this study, we examine MHC IIβ haplotype diversity in a model shark species, the small spotted catshark, Scyliorhinus canicula. Making use of pedigree data, we reconstructed MHC IIβ haplotypes to understand allele transmission from parent to offspring. Results indicate allele co-segregation consistent with tight linkage among MHC IIβ genes, suggesting the presence of functional stable haplotypes inherited from parents to offspring. The reconstructed haplotypes suggested extensive haplotype diversity characterized by variable allele numbers and allelic lineage composition, as well as marked allelic divergence, consistent with previous population-level data on this species. These findings underscore the complexity of MHC genetics (and of MHC evolution) in chondrichthyans. Accurate reconstruction of MHC haplotypes and assessment of its functional significance are crucial for better understanding adaptive immune responses and MHC evolutionary dynamics in chondrichthyans.

Translocations spur population growth but fail to prevent genetic erosion in imperiled Florida Scrub-Jays

Land and natural resource use in addition to climate change can restrict populations to degraded and fragmented habitats, catalyzing extinction through the reinforced interplay of small population size and genetic decay. Translocating individuals is a powerful strategy for overcoming direct threats from human development and reconnecting isolated populations but is not without risks.1 Habitat Management Plan analyses under section 7 of the U.S. Endangered Species Act determined that multiple subpopulations of Federally Threatened Florida Scrub-Jays (Aphelocoma coerulescens, hereafter FSJ) belonging to a metapopulation on Florida's west coast were declining demographic sinks, occupying areas where agriculture and fire suppression had degraded and fragmented the habitat.2 In order to increase the viability of the overall metapopulation, 51 FSJs from five of these small subpopulations in areas to be mined were translocated throughout 2003-2010 into a larger site of more contiguous, recently restored habitat at the core of the metapopulation, which contained a small resident population.3 Prior to translocations and for nearly two decades afterward, this core population, referred to as the M4 core region (CR) population, was extensively monitored, yielding a nearly complete pedigree. We used this pedigree, along with temporal genomic analyses and simulations, to show that translocations coupled with habitat restoration generated rapid population growth, but high reproductive skew increased inbreeding and led to genetic erosion. This mechanistic understanding of mixed conservation outcomes highlights the importance of monitoring and the potential need for genetic rescue to offset consequences of reproductive skew following translocations, regardless of demographic recovery.

Genomes of Galápagos Mockingbirds Reveal the Impact of Island Size and Past Demography on Inbreeding and Genetic Load in Contemporary Populations

Restricted range size brings about noteworthy genetic consequences that may affect the viability of a population and eventually its extinction. Particularly, the question if an increase in inbreeding can avert the accumulation of genetic load via purging is hotly debated in the conservation genetic field. Insular populations with limited range sizes represent an ideal setup for relating range size to these genetic factors. Leveraging a set of eight differently sized populations of Galápagos mockingbirds (Mimus), we investigated how island size shaped effective population size (Ne), inbreeding and genetic load. We assembled a genome of M. melanotis and genotyped three individuals per population by whole-genome resequencing. Demographic inference showed that the Ne of most populations remained high after the colonisation of the archipelago 1-2 Mya. Ne decline in M. parvulus happened only 10-20 Kya, whereas the critically endangered M. trifasciatus showed a longer history of reduced Ne. Despite these historical fluctuations, the current island size determines Ne in a linear fashion. In contrast, significant inbreeding coefficients, derived from runs of homozygosity, were identified only in the four smallest populations. The index of additive genetic load suggested purging in M. parvulus, where the smallest populations showed the lowest load. By contrast, M. trifasciatus carried the highest genetic load, possibly due to a recent rapid bottleneck. Overall, our study demonstrates a complex effect of demography on inbreeding and genetic load, providing implications in conservation genetics in general and in a conservation project of M. trifasciatus in particular.

To genetic rescue or not?

Inbreeding depression and genetic rescue are central themes in conservation biology. Translocation is a tool to assist genetic rescue but is connected to risks. A new study by Quinn et al. used genomic data to evaluate translocations as a potential action in montane red fox, bringing important implications also for other threatened species.

Contrasting Genomic Diversity and Inbreeding Levels Among Two Closely Related Falcon Species With Overlapping Geographic Distributions

Genomic resources are valuable to examine historical demographic patterns and their effects to better inform management and conservation of threatened species. We evaluated population trends and genome-wide variation in the near-threatened Orange-breasted Falcon (Falco deiroleucus) and its more common sister species, the Bat Falcon (F. rufigularis), to explore how the two species differ in genomic diversity as influenced by their contrasting long-term demographic histories. We generated and aligned whole genome resequencing data for 12 Orange-breasted Falcons and 9 Bat Falcons to an annotated Gyrfalcon (F. rusticolus) reference genome that retained approximately 22.4 million biallelic autosomal SNPs (chromosomes 1-22). Our analyses indicated much lower genomic diversity in Orange-breasted Falcons compared to Bat Falcons. All sampled Orange-breasted Falcons were significantly more inbred than the sampled Bat Falcons, with values similar to those observed in island-mainland species comparisons. The distribution of runs of homozygosity showed variation suggesting long-term low population size and the possibility of bottlenecks in Orange-breasted Falcons contrasting with consistently larger populations in Bat Falcons. Analysis of genetic load suggests that Orange-breasted Falcons are less likely to experience inbreeding depression than Bat Falcons due to reduced inbreeding load but are at elevated risk from fixation of deleterious gene variants and perhaps a reduced adaptive potential. These genomic analyses highlight differences in the historical demography of two closely related species that have influenced their current genomic diversity and should result in differing strategies for their continued conservation.

Comparative Population Genomics Reveal the Determinants of Genome Erosion in Two Sympatric Neotropical Falcons

Studying genetic diversity in endangered species has become an important component of conservation science over the past decades. Thanks to recent developments in sequencing technologies and bioinformatics, genetic parameters of conservation relevance such as neutral and functional genome-wide variation are now routinely estimated. Since inbreeding and deleterious mutations represent significant threats to small and declining populations, assessing the dynamics of these parameters has received particular attention in many recent conservation genomics studies. In this issue of Molecular Ecology, Martin et al. analyse the genomes of two Neotropical falcon species to assess the impact of their contrasting population histories on genome-wide diversity. They show that the Orange-breasted falcon which has had a low long-term population size and has experienced recent population bottlenecks is more inbred but has relatively fewer deleterious variations compared to its sister taxon, the Bat falcon, which is characterised by a larger long-term population size. This study not only provides insights into the role of past demography on the dynamics of deleterious variation in two species with contrasting population histories but also highlights the increasing importance of comparative approaches in population and conservation genomics.

Genomic Insights Into Inbreeding and Adaptive Divergence of Trout Populations to Inform Genetic Rescue

Genetic rescue, specifically translocation to facilitate gene flow among populations and reduce the effects of inbreeding, is an increasingly used approach in conservation. However, this approach comes with trade-offs, wherein gene flow may reduce fitness when populations have adaptive differentiation (i.e., outbreeding depression). A better understanding of the interaction between isolation, inbreeding, and adaptive divergence in key traits, such as life history traits, will help to inform genetic rescue efforts. Stream-dwelling salmonids, such as the westslope cutthroat trout (Oncorhynchus lewisi; WCT), are well-suited for examining these trade-offs because they are increasingly isolated by habitat degradation, exhibit substantial variation in life history traits among populations, and include many species of conservation concern. However, few genomic studies have examined the potential trade-offs in inbreeding versus outbreeding depression in salmonids. We used > 150,000 SNPs to examine genomic variation and inbreeding coefficients in 565 individuals across 25 WCT populations that differed in their isolation status and demographic histories. Analyses of runs of homozygosity revealed that several isolated WCT populations had "flatlined" having extremely low genetic variation and high inbreeding coefficients. Additionally, we conducted genome scans to identify potential outlier loci that could explain life history differences among 10 isolated populations. Genome scans identified one candidate genomic region that influenced maximum length and age-1 to age-2 growth. However, the limited number of candidate loci suggests that the life history traits examined may be driven by many genes of small effect or phenotypic plasticity. Although adaptive differentiation should be considered, the high inbreeding coefficients in several populations suggest that genetic rescue may benefit the most genetically depauperate WCT populations.

Genomic consequences of intensive inbreeding in miniature inbred pigs

BACKGROUND

Inbreeding, a central theme in evolutionary and conservation biology, is a crucial practice in breeding to stabilize and enhance the specific traits or to establish inbred lines. It also carries the risk of inbreeding depression, reduced fitness, and increased potential for extinction. Nevertheless, inbreeding has been extensively studied in small and endangered populations but its effects in large domesticated animals are poorly understood. Here, we aim to investigate the genomic consequences of inbreeding in the Banna miniature inbred pig (BN), a breed that has been inbred for over 40 years.

RESULTS

We have sequenced 41 genomes of BN and Diannan miniature pig (DN) at high-coverage (> 31×) and combined them with published whole-genomes of swine to comprehensively investigate the genetic consequences of inbreeding. We find that BN is genetically closely related to DN, which is consistent with breeding records. All families of BN have undergone an extreme bottleneck due to intensive inbreeding, resulting in higher genomic inbreeding coefficients, reduced genetic diversity, and a lower effective population size (Ne) compare to non-inbred pigs. Furthermore, BN and DN exhibit an increased genetic load relative to Asian wild boars. Prolonged inbreeding and bottlenecks have led to some purging of deleterious mutations in BN compared to DN, and a conversion from masked load to realized load.

CONCLUSIONS

We present a comprehensive analysis to understand and assess the consequences of inbreeding in miniature inbred pigs from a perspective of population genomics. Utilizing genomic measurements proves effective in estimating the consequences of inbreeding, especially when a detailed and accurate historical record of pedigree are lacking. Our results provide valuable resources and a detailed perspective on the genomic impacts of inbreeding, potentially guiding efforts in breeding, breed improvement, and conservation.

Museomics and morphological analyses of historical and contemporary peninsular Italian wolf (Canis lupus italicus) samples

After centuries of decline and protracted bottlenecks, the peninsular Italian wolf population has naturally recovered. However, an exhaustive comprehension of the effects of such a conservation success is still limited by the reduced availability of historical data. Therefore, in this study, we morphologically and genetically analyzed historical and contemporary wolf samples, also exploiting the optimization of an innovative bone DNA extraction method, to describe the morphological variability of the subspecies and its genetic diversity during the last 30 years. We obtained high amplification and genotyping success rates for tissue, blood and also petrous bone DNA samples. Multivariate, clustering and variability analyses confirmed that the Apennine wolf population is genetically and morphologically well-distinguishable from both European wolves and dogs, with no natural immigration from other populations, while its genetic variability has remained low across the last three decades, without significant changes between historical and contemporary specimens. This study highlights the scientific value of well-maintained museum collections, demonstrates that petrous bones represent reliable DNA sources, and emphasizes the need to genetically long-term monitor the dynamics of peculiar wolf populations to ensure appropriate conservation management actions.

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.

Beyond population size: Whole-genome data reveal bottleneck legacies in the peninsular Italian wolf

Preserving genetic diversity and adaptive potential while avoiding inbreeding depression is crucial for the long-term conservation of natural populations. Despite demographic increases, traces of past bottleneck events at the genomic level should be carefully considered for population management. From this perspective, the peninsular Italian wolf is a paradigmatic case. After being on the brink of extinction in the late 1960s, peninsular Italian wolves rebounded and recolonized most of the peninsula aided by conservation measures, including habitat and legal protection. Notwithstanding their demographic recovery, a comprehensive understanding of the genomic consequences of the historical bottleneck in Italian wolves is still lacking. To fill this gap, we sequenced whole genomes of 13 individuals sampled in the core historical range of the species in Central Italy to conduct population genomic analyses, including a comparison with wolves from two highly-inbred wolf populations (i.e. Scandinavia and Isle Royale). We found that peninsular Italian wolves, despite their recent recovery, still exhibit relatively low genetic diversity, a small effective population size, signatures of inbreeding, and a non-negligible genetic load. Our findings indicate that the peninsular Italian wolf population is still susceptible to bottleneck legacies, which could lead to local inbreeding depression in case of population reduction or fragmentations. This study emphasizes the importance of considering key genetic parameters to design appropriate long-term conservation management plans.

Understanding the Past to Preserve the Future: Genomic Insights Into the Conservation Management of a Critically Endangered Waterbird

To ensure the success of genetic rescue, we must minimise the potential negative effects of outbreeding depression that may arise from selecting source populations. The difficulty in assessing the likelihood of outbreeding depression has hindered its consideration in endangered species conservation. However, genomic research offers feasible indications. Here, we conduct conservation genomic analyses on the East Asian (EA) population and the relict Western/Central Asian (WCA) population of the critically endangered Siberian crane (Leucogeranus leucogeranus). We aim to assess genetic rescue's potential advantages and disadvantages between the two populations. Our analysis shows genomic evidence of limited genetic differentiation between them. The persistent decline in population size due to historical climatic oscillations leads to a decrease in genetic diversity and an increase in inbreeding. The WCA population has excessive deleterious homozygous mutations than the EA population, suggesting suffering from inbreeding depression resulting from less effective purifying selection. Forward simulations support the increase in genetic load due to elevated levels of inbreeding compromises fitness during population collapse. We strongly recommend an urgent genetic rescue for the WCA population through population supplementation from the EA population. However, the continuous monitoring of fitness outcomes is required through captive breeding. This work provides useful insights into the genetic management of a critically endangered species and emphasises the importance of evaluating the likelihood of outbreeding depression through genomic approaches.

Sexual selection matters in genetic rescue, but productivity benefits fade over time: a multi-generation experiment to inform conservation

Globally, many species are threatened by population decline because of anthropogenic changes leading to population fragmentation, genetic isolation and inbreeding depression. Genetic rescue, the controlled introduction of genetic variation, is a method used to relieve such effects in small populations. However, without understanding how the characteristics of rescuers impact rescue attempts interventions run the risk of being sub-optimal, or even counterproductive. We use the red flour beetle (Tribolium castaneum) to test the impact of rescuer sex, and sexual selection background, on population productivity. We record the impact of genetic rescue on population productivity in 24 and 36 replicated populations for ten generations following intervention. We find little or no impact of rescuer sex on the efficacy of rescue but show that a background of elevated sexual selection makes individuals more effective rescuers. In both experiments, rescue effects diminish 6-10 generations after the rescue. Our results confirm that the efficacy of genetic rescue can be influenced by characteristics of the rescuers and that the level of sexual selection in the rescuing population is an important factor. We show that any increase in fitness associated with rescue may last for a limited number of generations, suggesting implications for conservation policy and practice.

Ecological and anthropogenic effects on the genomic diversity of lemurs in Madagascar

Ecological variation and anthropogenic landscape modification have had key roles in the diversification and extinction of mammals in Madagascar. Lemurs represent a radiation with more than 100 species, constituting roughly one-fifth of the primate order. Almost all species of lemurs are threatened with extinction, but little is known about their genetic diversity and demographic history. Here, we analyse high-coverage genome-wide resequencing data from 162 unique individuals comprising 50 species of Lemuriformes, including multiple individuals from most species. Genomic diversity varies widely across the infraorder and yet is broadly consistent among individuals within species. We show widespread introgression in multiple genera and generally high levels of genomic diversity likely resulting from allele sharing that occurred during periods of connectivity and fragmentation during climatic shifts. We find distinct patterns of demographic history in lemurs across the ecogeographic regions of Madagascar within the last million years. Within the past 2,000 years, lemurs underwent major declines in effective population size that corresponded to the timing of human population expansion in Madagascar. In multiple regions of the island, we identified chronological trajectories of inbreeding that are consistent across genera and species, suggesting localized effects of human activity. Our results show how the extraordinary diversity of these long-neglected, endangered primates has been influenced by ecological and anthropogenic factors.

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.

Recently Delisted Songbird Harbors Extensive Genomic Evidence of Inbreeding, Potentially Complicating Future Recovery

The Kirtland's warbler (Setophaga kirtlandii) is a rare migratory passerine species and habitat specialist of the North American Jack Pine Forests. Their near extinction in the 1970s classified them as endangered and protected under the Endangered Species Act of 1973. After decades of intense conservation management, their population size recovered, and they were delisted from federal protection in 2019. We explore the genomic consequences of this harsh bottleneck and recovery by comparing the genomic architecture of two closely related species whose population sizes have remained large and stable, Hooded Warblers (Setophaga citrina) and American Redstarts (Setophaga ruticilla). We used whole-genome sequencing to characterize the distribution of runs of homozygosity and deleterious genetic variation. We find evidence that Kirtland's warblers exhibit genetic patterns consistent with recent inbreeding. Our results also show that Kirtland's warblers carry an excess proportion of deleterious variation, which could complicate management for this conservation-reliant species. This analysis provides a genetically informed perspective that should be thoroughly considered when delisting other species from federal protections. Through the increasing accessibility of genome sequencing technology, it will be more feasible to monitor the genetic landscape of recovering populations to ensure their long-term survival independent of conservation intervention.

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.

Estimation of spatial demographic maps from polymorphism data using a neural network

A fundamental goal in population genetics is to understand how variation is arrayed over natural landscapes. From first principles we know that common features such as heterogeneous population densities and barriers to dispersal should shape genetic variation over space, however there are few tools currently available that can deal with these ubiquitous complexities. Geographically referenced single nucleotide polymorphism (SNP) data are increasingly accessible, presenting an opportunity to study genetic variation across geographic space in myriad species. We present a new inference method that uses geo-referenced SNPs and a deep neural network to estimate spatially heterogeneous maps of population density and dispersal rate. Our neural network trains on simulated input and output pairings, where the input consists of genotypes and sampling locations generated from a continuous space population genetic simulator, and the output is a map of the true demographic parameters. We benchmark our tool against existing methods and discuss qualitative differences between the different approaches; in particular, our program is unique because it infers the magnitude of both dispersal and density as well as their variation over the landscape, and it does so using SNP data. Similar methods are constrained to estimating relative migration rates, or require identity-by-descent blocks as input. We applied our tool to empirical data from North American grey wolves, for which it estimated mostly reasonable demographic parameters, but was affected by incomplete spatial sampling. Genetic based methods like ours complement other, direct methods for estimating past and present demography, and we believe will serve as valuable tools for applications in conservation, ecology and evolutionary biology. An open source software package implementing our method is available from https://github.com/kr-colab/mapNN.

Detectability of runs of homozygosity is influenced by analysis parameters and population-specific demographic history

Wild populations are increasingly threatened by human-mediated climate change and land use changes. As populations decline, the probability of inbreeding increases, along with the potential for negative effects on individual fitness. Detecting and characterizing runs of homozygosity (ROHs) is a popular strategy for assessing the extent of individual inbreeding present in a population and can also shed light on the genetic mechanisms contributing to inbreeding depression. Here, we analyze simulated and empirical datasets to demonstrate the downstream effects of program selection and long-term demographic history on ROH inference, leading to context-dependent biases in the results. Through a sensitivity analysis we evaluate how various parameter values impact ROH-calling results, highlighting its utility as a tool for parameter exploration. Our results indicate that ROH inferences are sensitive to factors such as sequencing depth and ROH length distribution, with bias direction and magnitude varying with demographic history and the programs used. Estimation biases are particularly pronounced at lower sequencing depths, potentially leading to either underestimation or overestimation of inbreeding. These results are particularly important for the management of endangered species, as underestimating inbreeding signals in the genome can substantially undermine conservation initiatives. We also found that small true ROHs can be incorrectly lumped together and called as longer ROHs, leading to erroneous inference of recent inbreeding. To address these challenges, we suggest using a combination of ROH detection tools and ROH length-specific inferences, along with sensitivity analysis, to generate robust and context-appropriate population inferences regarding inbreeding history. We outline these recommendations for ROH estimation at multiple levels of sequencing effort, which are typical of conservation genomics studies.

Whole Genomes Inform Genetic Rescue Strategy for Montane Red Foxes in North America

A few iconic examples have proven the value of facilitated gene flow for counteracting inbreeding depression and staving off extinction; yet, the practice is often not implemented for fear of causing outbreeding depression. Using genomic sequencing, climatic niche modeling, and demographic reconstruction, we sought to assess the risks and benefits of using translocations as a tool for recovery of endangered montane red fox (Vulpes vulpes) populations in the western United States. We demonstrated elevated inbreeding and homozygosity of deleterious alleles across all populations, but especially those isolated in the Cascade and Sierra Nevada ranges. Consequently, translocations would be expected to increase population growth by masking deleterious recessive alleles. Demographic reconstructions further indicated shallow divergences of less than a few thousand years among montane populations, suggesting low risk of outbreeding depression. These genomic-guided findings set the stage for future management, the documentation of which will provide a roadmap for recovery of other data-deficient taxa.

Constraining models of dominance for nonsynonymous mutations in the human genome

Dominance is a fundamental parameter in genetics, determining the dynamics of natural selection on deleterious and beneficial mutations, the patterns of genetic variation in natural populations, and the severity of inbreeding depression in a population. Despite this importance, dominance parameters remain poorly known, particularly in humans or other non-model organisms. A key reason for this lack of information about dominance is that it is extremely challenging to disentangle the selection coefficient (s) of a mutation from its dominance coefficient (h). Here, we explore dominance and selection parameters in humans by fitting models to the site frequency spectrum (SFS) for nonsynonymous mutations. When assuming a single dominance coefficient for all nonsynonymous mutations, we find that numerous h values can fit the data, so long as h is greater than ~0.15. Moreover, we also observe that theoretically-predicted models with a negative relationship between h and s can also fit the data well, including models with h = 0.05 for strongly deleterious mutations. Finally, we use our estimated dominance and selection parameters to inform simulations revisiting the question of whether the out-of-Africa bottleneck has led to differences in genetic load between African and non-African human populations. These simulations suggest that the relative burden of genetic load in non-African populations depends on the dominance model assumed, with slight increases for more weakly recessive models and slight decreases shown for more strongly recessive models. Moreover, these results also demonstrate that models of partially recessive nonsynonymous mutations can explain the observed severity of inbreeding depression in humans, bridging the gap between molecular population genetics and direct measures of fitness in humans. Our work represents a comprehensive assessment of dominance and deleterious variation in humans, with implications for parameterizing models of deleterious variation in humans and other mammalian species.

Genomic divergence and mutation load in the Begonia masoniana complex from limestone karsts

Understanding genome-wide diversity, inbreeding, and the burden of accumulated deleterious mutations in small and isolated populations is essential for predicting and enhancing population persistence and resilience. However, these effects are rarely studied in limestone karst plants. Here, we re-sequenced the nuclear genomes of 62 individuals of the Begonia masoniana complex (B. liuyanii, B. longgangensis, B. masoniana and B. variegata) and investigated genomic divergence and genetic load for these four species. Our analyses revealed four distinct clusters corresponding to each species within the complex. Notably, there was only limited admixture between B. liuyanii and B. longgangensis occurring in overlapping geographic regions. All species experienced historical bottlenecks during the Pleistocene, which were likely caused by glacial climate fluctuations. We detected an asymmetric historical gene flow between group pairs within this timeframe, highlighting a distinctive pattern of interspecific divergence attributable to karst geographic isolation. We found that isolated populations of B. masoniana have limited gene flow, the smallest recent population size, the highest inbreeding coefficients, and the greatest accumulation of recessive deleterious mutations. These findings underscore the urgency to prioritize conservation efforts for these isolated population. This study is among the first to disentangle the genetic differentiation and specific demographic history of karst Begonia plants at the whole-genome level, shedding light on the potential risks associated with the accumulation of deleterious mutations over generations of inbreeding. Moreover, our findings may facilitate conservation planning by providing critical baseline genetic data and a better understanding of the historical events that have shaped current population structure of rare and endangered karst plants.

Ecological genetics of isolated loach populations indicate compromised adaptive potential

Many endangered species live in fragmented and isolated populations with low genetic variability, signs of inbreeding, and small effective population sizes - all features elevating their extinction risk. The flat-headed loach (Oreonectes platycephalus), a small noemacheilid fish, is widely across southern China, but only in the headwaters of hillstreams; as a result, they are spatially isolated from conspecific populations. We surveyed single nucleotide polymorphisms in 16 Hong Kong populations of O. platycephalus to determine whether loach populations from different streams were genetically isolated from each other, showed low levels of genetic diversity, signs of inbreeding, and had small contemporary effective population sizes. Estimates of average observed heterozygosity (HO = 0.0473), average weighted nucleotide diversity (πw = 0.0546) and contemporary effective population sizes (Ne = 10.2 ~ 129.8) were very low, and several populations showed clear signs of inbreeding as judged from relatedness estimates. The degree of genetic differentiation among populations was very high (average FST = 0.668), even over short geographic distances (<1.5 km), with clear patterns of isolation by distance. These results suggest that Hong Kong populations of O. platycephalus have experienced strong genetic drift and loss of genetic variability because sea-level rise after the last glaciation reduced connectedness among paleodrainages, isolating populations in headwaters. All this, together with the fact that the levels of genetic diversity and contemporary effective population sizes within O. platycephalus populations are lower than most other freshwater fishes, suggests that they face high local extinction risk and have limited capacity for future adaptation.

Estimation of spatial demographic maps from polymorphism data using a neural network

A fundamental goal in population genetics is to understand how variation is arrayed over natural landscapes. From first principles we know that common features such as heterogeneous population densities and barriers to dispersal should shape genetic variation over space, however there are few tools currently available that can deal with these ubiquitous complexities. Geographically referenced single nucleotide polymorphism (SNP) data are increasingly accessible, presenting an opportunity to study genetic variation across geographic space in myriad species. We present a new inference method that uses geo-referenced SNPs and a deep neural network to estimate spatially heterogeneous maps of population density and dispersal rate. Our neural network trains on simulated input and output pairings, where the input consists of genotypes and sampling locations generated from a continuous space population genetic simulator, and the output is a map of the true demographic parameters. We benchmark our tool against existing methods and discuss qualitative differences between the different approaches; in particular, our program is unique because it infers the magnitude of both dispersal and density as well as their variation over the landscape, and it does so using SNP data. Similar methods are constrained to estimating relative migration rates, or require identity by descent blocks as input. We applied our tool to empirical data from North American grey wolves, for which it estimated mostly reasonable demographic parameters, but was affected by incomplete spatial sampling. Genetic based methods like ours complement other, direct methods for estimating past and present demography, and we believe will serve as valuable tools for applications in conservation, ecology, and evolutionary biology. An open source software package implementing our method is available from https://github.com/kr-colab/mapNN .

Patterns of reproduction and autozygosity distinguish the breeding from nonbreeding gray wolves of Yellowstone National Park

For species of management concern, accurate estimates of inbreeding and associated consequences on reproduction are crucial for predicting their future viability. However, few studies have partitioned this aspect of genetic viability with respect to reproduction in a group-living social mammal. We investigated the contributions of foundation stock lineages, putative fitness consequences of inbreeding, and genetic diversity of the breeding versus nonreproductive segment of the Yellowstone National Park gray wolf population. Our dataset spans 25 years and seven generations since reintroduction, encompassing 152 nuclear families and 329 litters. We found more than 87% of the pedigree foundation genomes persisted and report influxes of allelic diversity from two translocated wolves from a divergent source in Montana. As expected for group-living species, mean kinship significantly increased over time but with minimal loss of observed heterozygosity. Strikingly, the reproductive portion of the population carried a significantly lower genome-wide inbreeding coefficients, autozygosity, and more rapid decay for linkage disequilibrium relative to the nonbreeding population. Breeding wolves had significantly longer lifespans and lower inbreeding coefficients than nonbreeding wolves. Our model revealed that the number of litters was negatively significantly associated with heterozygosity (R = -0.11). Our findings highlight genetic contributions to fitness, and the importance of the reproductively active individuals in a population to counteract loss of genetic variation in a wild, free-ranging social carnivore. It is crucial for managers to mitigate factors that significantly reduce effective population size and genetic connectivity, which supports the dispersion of genetic variation that aids in rapid evolutionary responses to environmental challenges.

The difficulty of detecting inbreeding depression and its effect on conservation decisions

Statistical inferences about inbreeding depression are often derived from analyses with low power and a high risk of failing to detect inbreeding depression. That risk is widely appreciated by scientists familiar with the relevant statistical and genetical theory, but may be overlooked and underappreciated by decision-makers. Consequently, there is value in demonstrating this risk using a real example. We use data from the wolf population on Isle Royale to demonstrate the difficulty of making reliable statistical inferences about inbreeding depression. This wolf population is known-by other methods-to have gone effectively extinct due to deleterious genetic processes associated with inbreeding. Beyond that demonstration, we use two case-studies-wolves on Isle Royale and vaquita (porpoises) from the Gulf of California, Mexico-to show how statistical inferences about inbreeding depression can affect conservation decisions. According to most decision theory, decisions depend importantly on: 1) probabilities that certain states exist (e.g. inbreeding depression is present) and 2) the utility assigned to various outcomes (e.g. the value of acting to mitigate inbreeding when it is present). The probabilities are provided by statistical inference; whereas utilities are almost entirely determined by normative values and judgements. Our analysis suggests that decisions to mitigate inbreeding depression are often driven more by utilities (normative values) than probabilities (statistical inferences). As such, advocates for mitigating inbreeding depression will benefit from better communicating to decision-makers the value of populations persisting and the extent to which decisions should depend on normative values.

Whole-genome analysis reveals the diversification of Galapagos rail (Aves: Rallidae) and confirms the success of goat eradication programs

Similar to other insular birds around the world, the Galapagos rail (Laterallus spilonota Gould, 1841) exhibits reduced flight capacity following its colonization of the archipelago ~1.2 mya. Despite their short evolutionary history, rails have colonized seven different islands spanning the entire width of the archipelago. Galapagos rails were once common on islands with sufficiently high altitudes to support shrubs in humid habitats. After humans introduced goats, this habitat was severely reduced due to overgrazing. Habitat loss devastated some rail populations, with less than 50 individuals surviving, rendering the genetic diversity of Galapagos rail a pressing conservation concern. Additionally, one enigma is the reappearance of rails on the island of Pinta after they were considered extirpated. Our approach was to investigate the evolutionary history and geographic distribution of Galapagos rails as well as examine the genome-wide effects of historical population bottlenecks using 39 whole genomes across different island populations. We recovered an early divergence of rail ancestors leading to the isolated populations on Pinta and a second clade comprising the rest of the islands, historically forming a single landmass. Subsequently, the separation of the landmass ~900 kya may have led to the isolation of the Isabela population with more panmictic populations found on Santa Cruz and Santiago islands. We found that rails genomes contain long runs of homozygosity (>2 Mb) that could be related to the introduction of goats. Finally, our findings show that the modern eradication of goats was critical to avoiding episodes of inbreeding in most populations.

Large variance in inbreeding within the Iberian wolf population

The gray wolf (Canis lupus) population on the Iberian Peninsula was the largest in western and central Europe during most of the 20th century, with its size apparently never under a few hundred individuals. After partial legal protection in the 1970s in Spain, the northwest Iberian population increased to about 300 to 350 packs and then stabilized. In contrast to many current European wolf populations, which have been connected through gene flow, the Iberian wolf population has been isolated for decades. Here, we measured changes in genomic diversity and inbreeding through the last decades in a geographic context. We find that the level of genomic diversity in Iberian wolves is low compared with other Eurasian wolf populations. Despite population expansion in the last 50 years, some modern wolves had very high inbreeding, especially in the recently recolonized and historical edge areas. These individuals contrast with others with low inbreeding within the same population. The high variance in inbreeding despite population expansion seems associated with small-scale fragmentation of the range that is revealed by the genetic similarity between modern and historical samples from close localities despite being separated by decades, remaining differentiated from other individuals that are just over 100 km away, a small distance for a species with great dispersal capacity inhabiting a continuous range. This illustrates that, despite its demographically stable condition, the population would probably benefit from favoring connectivity within the population as well as genetic exchange with other European wolf populations to avoid excessive fragmentation and local inbreeding depression.

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.

Genomic Consequences of Isolation and Inbreeding in an Island Dingo Population

Dingoes come from an ancient canid lineage that originated in East Asia around 8,000 to 11,000 years BP. As Australia's largest terrestrial predator, dingoes play an important ecological role. A small, protected population exists on a world heritage listed offshore island, K'gari (formerly Fraser Island). Concern regarding the persistence of dingoes on K'gari has risen due to their low genetic diversity and elevated inbreeding levels. However, whole-genome sequence data is lacking from this population. Here, we include five new whole-genome sequences of K'gari dingoes. We analyze a total of 18 whole-genome sequences of dingoes sampled from mainland Australia and K'gari to assess the genomic consequences of their demographic histories. Long (>1 Mb) runs of homozygosity (ROHs)-indicators of inbreeding-are elevated in all sampled dingoes. However, K'gari dingoes showed significantly higher levels of very long ROH (>5 Mb), providing genomic evidence for small population size, isolation, inbreeding, and a strong founder effect. Our results suggest that, despite current levels of inbreeding, the K'gari population is purging strongly deleterious mutations, which, in the absence of further reductions in population size, may facilitate the persistence of small populations despite low genetic diversity and isolation. However, there may be little to no purging of mildly deleterious alleles, which may have important long-term consequences, and should be considered by conservation and management programs.

Imputation of ancient canid genomes reveals inbreeding history over the past 10,000 years

The multi-millenia long history between dogs and humans has placed them at the forefront of archeological and genomic research. Despite ongoing efforts including the analysis of ancient dog and wolf genomes, many questions remain regarding their geographic and temporal origins, and the microevolutionary processes that led to the diversity of breeds today. Although ancient genomes provide valuable information, their use is hindered by low depth of coverage and post-mortem damage, which inhibits confident genotype calling. In the present study, we assess how genotype imputation of ancient dog and wolf genomes, utilising a large reference panel, can improve the resolution provided by ancient datasets. Imputation accuracy was evaluated by down-sampling high coverage dog and wolf genomes to 0.05-2x coverage and comparing concordance between imputed and high coverage genotypes. We measured the impact of imputation on principal component analyses and runs of homozygosity. Our findings show high (R2>0.9) imputation accuracy for dogs with coverage as low as 0.5x and for wolves as low as 1.0x. We then imputed a dataset of 90 ancient dog and wolf genomes, to assess changes in inbreeding during the last 10,000 years of dog evolution. Ancient dog and wolf populations generally exhibited lower inbreeding levels than present-day individuals. Interestingly, regions with low ROH density maintained across ancient and present-day samples were significantly associated with genes related to olfaction and immune response. Our study indicates that imputing ancient canine genomes is a viable strategy that allows for the use of analytical methods previously limited to high-quality genetic data.

The first complete genome of the extinct European wild ass (Equus hemionus hydruntinus)

We present palaeogenomes of three morphologically unidentified Anatolian equids dating to the first millennium BCE, sequenced to a coverage of 0.6-6.4×. Mitochondrial DNA haplotypes of the Anatolian individuals clustered with those of Equus hydruntinus (or Equus hemionus hydruntinus), the extinct European wild ass, secular name 'hydruntine'. Further, the Anatolian wild ass whole genome profiles fell outside the genomic diversity of other extant and past Asiatic wild ass (E. hemionus) lineages. These observations suggest that the three Anatolian wild asses represent hydruntines, making them the latest recorded survivors of this lineage, about a millennium later than the latest observations in the zooarchaeological record. Our mitogenomic and genomic analyses indicate that E. h. hydruntinus was a clade belonging to ancient and present-day E. hemionus lineages that radiated possibly between 0.6 and 0.8 Mya. We also find evidence consistent with recent gene flow between hydruntines and Middle Eastern wild asses. Analyses of genome-wide heterozygosity and runs of homozygosity suggest that the Anatolian wild ass population may have lost genetic diversity by the mid-first millennium BCE, a possible sign of its eventual demise.

The role of Allee effects for Gaussian and Lévy dispersals in an environmental niche

In the study of biological populations, the Allee effect detects a critical density below which the population is severely endangered and at risk of extinction. This effect supersedes the classical logistic model, in which low densities are favorable due to lack of competition, and includes situations related to deficit of genetic pools, inbreeding depression, mate limitations, unavailability of collaborative strategies due to lack of conspecifics, etc. The goal of this paper is to provide a detailed mathematical analysis of the Allee effect. After recalling the ordinary differential equation related to the Allee effect, we will consider the situation of a diffusive population. The dispersal of this population is quite general and can include the classical Brownian motion, as well as a Lévy flight pattern, and also a "mixed" situation in which some individuals perform classical random walks and others adopt Lévy flights (which is also a case observed in nature). We study the existence and nonexistence of stationary solutions, which are an indication of the survival chance of a population at the equilibrium. We also analyze the associated evolution problem, in view of monotonicity in time of the total population, energy consideration, and long-time asymptotics. Furthermore, we also consider the case of an "inverse" Allee effect, in which low density populations may access additional benefits.

The Maintenance of Deleterious Variation in Wild Chinese Rhesus Macaques

Understanding how deleterious variation is shaped and maintained in natural populations is important in conservation and evolutionary biology, as decreased fitness caused by these deleterious mutations can potentially lead to an increase in extinction risk. It is known that demographic processes can influence these patterns. For example, population bottlenecks and inbreeding increase the probability of inheriting identical-by-descent haplotypes from a recent common ancestor, creating long tracts of homozygous genotypes called runs of homozygosity (ROH), which have been associated with an accumulation of mildly deleterious homozygotes. Counterintuitively, positive selection can also maintain deleterious variants in a population through genetic hitchhiking. Here, we analyze the whole genomes of 79 wild Chinese rhesus macaques across five subspecies and characterize patterns of deleterious variation with respect to ROH and signals of recent positive selection. We show that the fraction of homozygotes occurring in long ROH is significantly higher for deleterious homozygotes than tolerated ones, whereas this trend is not observed for short and medium ROH. This confirms that inbreeding, by generating these long tracts of homozygosity, is the main driver of the high burden of homozygous deleterious alleles in wild macaque populations. Furthermore, we show evidence that homozygous LOF variants are being purged. Next, we identify seven deleterious variants at high frequency in regions putatively under selection near genes involved with olfaction and other processes. Our results shed light on how evolutionary processes can shape the distribution of deleterious variation in wild nonhuman primates.

The Maintenance of Deleterious Variation in Wild Chinese Rhesus Macaques

Understanding how deleterious variation is shaped and maintained in natural populations is important in conservation and evolutionary biology, as decreased fitness caused by these deleterious mutations can potentially lead to an increase in extinction risk. It is known that demographic processes can influence these patterns. For example, population bottlenecks and inbreeding increase the probability of inheriting identical-by-descent haplotypes from a recent common ancestor, creating long tracts of homozygous genotypes called runs of homozygosity (ROH), which have been associated with an accumulation of mildly deleterious homozygotes. Counter intuitively, positive selection can also maintain deleterious variants in a population through genetic hitchhiking. Here we analyze the whole genomes of 79 wild Chinese rhesus macaques across five subspecies and characterize patterns of deleterious variation with respect to ROH and signals of recent positive selection. We show that the fraction of homozygotes occurring in long ROH is significantly higher for deleterious homozygotes than tolerated ones, whereas this trend is not observed for short and medium ROH. This confirms that inbreeding, by generating these long tracts of homozygosity, is the main driver of the high burden of homozygous deleterious alleles in wild macaque populations. Furthermore, we show evidence that homozygous LOF variants are being purged. Next, we identify 7 deleterious variants at high frequency in regions putatively under selection near genes involved with olfaction and other processes. Our results shed light on how evolutionary processes can shape the distribution of deleterious variation in wild non-human primates.

The origin, connectivity, and individual specialization of island wolves after deer extirpation

Wolves are assumed to be ungulate obligates, however, a recently described pack on Pleasant Island, Alaska USA, is persisting on sea otters and other marine resources without ungulate prey, violating this long-held assumption. We address questions about these wolves regarding their origin and fate, degree of isolation, risk of inbreeding depression, and diet specialization by individual and sex. We applied DNA metabarcoding and genotyping by amplicon sequencing using 957 scats collected from 2016 to 2022, and reduced representation sequencing of tissue samples to establish a detailed understanding of Pleasant Island wolf ecology and compare them with adjacent mainland wolves. Dietary overlap was higher among individual wolves on Pleasant Island (Pianka's index mean 0.95 ± 0.03) compared to mainland wolves (0.70 ± 0.21). The individual diets of island wolves were dominated by sea otter, ranging from 40.6% to 63.2% weighted percent of occurrence (wPOO) (mean 55.5 ± 8.7). In contrast, individual mainland wolves primarily fed on ungulates (42.2 ± 21.3) or voles during a population outbreak (31.2 ± 23.2). We traced the origin of the Pleasant Island pack to a mainland pair that colonized around 2013 and produced several litters. After this breeding pair was killed, their female offspring and an immigrant male became the new breeders in 2019. We detected 20 individuals of which 8 (40%) were trapped and killed while two died of natural causes during the 6-year study. Except for the new breeding male, the pedigree analysis and genotype results showed no additional movement to or from the island, indicating limited dispersal but no evidence of inbreeding. Our findings suggest wolves exhibit more flexible foraging behavior than previously believed, and hunting strategies can substantially differ between individuals within or between packs. Nevertheless, anthropogenic and natural mortality combined with limited connectivity to the mainland may inhibit the continued persistence of Pleasant Island wolves.

Relatedness-based mate choice and female philopatry: inbreeding trends of wolf packs in a human-dominated landscape

Inbreeding can reduce offspring fitness and has substantial implications for the genetic diversity and long-term viability of populations. In social cooperative canids, inbreeding is conditioned by the geographic proximity between opposite-sex kin outside natal groups and the presence of related individuals in neighbouring groups. Consequently, challenges in moving into other regions where the species is present can also affect inbreeding rates. These can be particularly problematic in areas of high human density, where movement can be restricted, even for highly vagile species. In this study, we investigate the socio-ecological dynamics of Iberian wolf packs in the human-dominated landscape of Alto Minho, in northwest Portugal, where wolves exhibit a high prevalence of short-distance dispersal and limited gene flow with neighbouring regions. We hypothesise that mating occurs regardless of relatedness, resulting in recurrent inbreeding due to high kin encounter rates. Using data from a 10-year non-invasive genetic monitoring programme and a combination of relatedness estimates and genealogical reconstructions, we describe genetic diversity, mate choice, and dispersal strategies among Alto Minho packs. In contrast with expectations, our findings reveal relatedness-based mate choice, low kin encounter rates, and a reduced number of inbreeding events. We observed a high prevalence of philopatry, particularly among female breeders, with the most common breeding strategy involving the pairing of a philopatric female with an unrelated immigrant male. Overall, wolves were not inbred, and temporal changes in genetic diversity were not significant. Our findings are discussed, considering the demographic trend of wolves in Alto Minho and its human-dominated landscape.

Genomic consequences of isolation and inbreeding in an island dingo population

Dingoes come from an ancient canid lineage that originated in East Asia around 8000-11,000 years BP. As Australia's largest terrestrial predator, dingoes play an important ecological role. A small, protected population exists on a world heritage listed offshore island, K'gari (formerly Fraser Island). Concern regarding the persistence of dingoes on K'gari has risen due to their low genetic diversity and elevated inbreeding levels. However, whole-genome sequencing data is lacking from this population. Here, we include five new whole-genome sequences of K'gari dingoes. We analyze a total of 18 whole genome sequences of dingoes sampled from mainland Australia and K'gari to assess the genomic consequences of their demographic histories. Long (>1 Mb) runs of homozygosity (ROH) - indicators of inbreeding - are elevated in all sampled dingoes. However, K'gari dingoes showed significantly higher levels of very long ROH (>5 Mb), providing genomic evidence for small population size, isolation, inbreeding, and a strong founder effect. Our results suggest that, despite current levels of inbreeding, the K'gari population is purging strongly deleterious mutations, which, in the absence of further reductions in population size, may facilitate the persistence of small populations despite low genetic diversity and isolation. However, there may be little to no purging of mildly deleterious alleles, which may have important long-term consequences, and should be considered by conservation and management programs.

A High-Quality Blue Whale Genome, Segmental Duplications, and Historical Demography

The blue whale, Balaenoptera musculus, is the largest animal known to have ever existed, making it an important case study in longevity and resistance to cancer. To further this and other blue whale-related research, we report a reference-quality, long-read-based genome assembly of this fascinating species. We assembled the genome from PacBio long reads and utilized Illumina/10×, optical maps, and Hi-C data for scaffolding, polishing, and manual curation. We also provided long read RNA-seq data to facilitate the annotation of the assembly by NCBI and Ensembl. Additionally, we annotated both haplotypes using TOGA and measured the genome size by flow cytometry. We then compared the blue whale genome with other cetaceans and artiodactyls, including vaquita (Phocoena sinus), the world's smallest cetacean, to investigate blue whale's unique biological traits. We found a dramatic amplification of several genes in the blue whale genome resulting from a recent burst in segmental duplications, though the possible connection between this amplification and giant body size requires further study. We also discovered sites in the insulin-like growth factor-1 gene correlated with body size in cetaceans. Finally, using our assembly to examine the heterozygosity and historical demography of Pacific and Atlantic blue whale populations, we found that the genomes of both populations are highly heterozygous and that their genetic isolation dates to the last interglacial period. Taken together, these results indicate how a high-quality, annotated blue whale genome will serve as an important resource for biology, evolution, and conservation research.

Patterns of Wolf Dispersal Respond to Harvest Density across an Island Complex

Wolves are highly mobile predators and can disperse across a variety of habitats and over long distances. However, less is known about dispersal capabilities across water and among islands. The biogeography of island systems fosters spatially structured local populations, and their degree of connectivity may influence the dynamics and long-term viability of the regional population. We sought to quantify wolf dispersal rate, distance, and dispersal sex bias throughout Prince of Wales Island, a 6670 km2 island in southeast Alaska, and the surrounding islands that constitute the wildlife management unit (9025 km2). We also investigated patterns of dispersal in relation to hunting and trapping intensity and wolf population density. We used DNA data collected during 2012-2021 long-term monitoring efforts and genotyped 811 wolves, 144 of which (18%) were dispersers. Annual dispersal rates were 9-23% and had a weakly positive relationship with wolf density. Wolves dispersed 41.9 km on average (SD = 23.7 km), and males and females did not disperse at different rates. Of the dispersing wolves, 107 died, and the majority (n = 81) died before they were able to settle. The leading manner of death was trapping (97% of mortalities), and wolves tended to disperse from areas with low harvest density to areas where harvest density was relatively higher. Dispersal occurred both to and from small islands and the larger Prince of Wales Island, indicating bidirectional as opposed to asymmetrical movement, and the genetic overlap of wolf groups demonstrates connectivity throughout this naturally patchy system. Island ecosystems have different predator-prey dynamics and recolonization processes than large, intact systems due to their isolation and restricted sizes; thus, a better understanding of the degree of population connectivity including dispersal patterns among islands in the Prince of Wales archipelago could help inform the management and research strategies of these wolves.

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.

Fitness consequences and ancestry loss in the Apennine brown bear after a simulated genetic rescue intervention

Reduction in population size, with its predicted effects on population fitness, is the most alarming anthropogenic impact on endangered species. By introducing compatible individuals, genetic rescue (GR) is a promising but debated approach for reducing the genetic load unmasked by inbreeding and for restoring the fitness of declining populations. Although GR can improve genetic diversity and fitness, it can also produce loss of ancestry, hampering local adaptation, or replace with introduced variants the unique genetic pools evolved in endemic groups. We used forward genetic simulations based on empirical genomic data to assess fitness benefits and loss of ancestry risks of GR in the Apennine brown bear (Ursus arctos marsicanus). There are approximately 50 individuals of this isolated subspecies, and they have lower genetic diversity and higher inbreeding than other European brown bears, and GR has been suggested to reduce extinction risks. We compared 10 GR scenarios in which the number and genetic characteristics of migrants varied with a non-GR scenario of simple demographic increase due to nongenetic factors. The introduction of 5 individuals of higher fitness or lower levels of deleterious mutations than the target Apennine brown bear from a larger European brown bear population produced a rapid 10-20% increase in fitness in the subspecies and up to 22.4% loss of ancestry over 30 generations. Without a contemporary demographic increase, fitness started to decline again after a few generations. Doubling the population size without GR gradually increased fitness to a comparable level, but without losing ancestry, thus resulting in the best strategy for the Apennine brown bear conservation. Our results highlight the importance for management of endangered species of realistic forward simulations grounded in empirical whole-genome data.

Conservation priority and run of homozygosity pattern assessment of global chicken genetic resources

The conservation of genetic resources is becoming increasingly important for the sustainable development of the poultry industry. In the present study, we systematically analyzed the population structure, conservation priority, runs of homozygosity (ROH) of chicken breeds globally, and proposed rational conservation strategies. We used a 600K Affymetrix Axiom HD genotyping SNP array dataset of 2,429 chickens from 134 populations. The chickens were divided into 5 groups based on their country of origin and sampling location: Asian chickens (AS-LOC), African chickens (AF), European local chickens (EU-LOC), Asian breeds sampled in Germany (AS-DE), and European breeds sampled in Germany (EU-DE). The results indicated that the population structure was consistent with the actual geographical distribution of the populations. AS-LOC had the highest positive contribution to the total gene (HT, 1.00%,) and allelic diversity (AT, 0.0014%), the lowest inbreeding degree and the fastest linkage disequilibrium (LD) decay rate; the lowest contribution are derived by European ex situ chicken breeds (EU-DE:HT = -0.072%, AT = -0.0014%), which showed the highest inbreeding and slowest LD decay. Breeds farmed in ex situ (AS-DE, EU-DE) conditions exhibited reduced genetic diversity and increased inbreeding due to small population size. Given limited funds, it is a better choice for government to conserve the breeds with the highest contribution to genetic diversity in each group. Therefore, we evaluated the contribution of each breed to genetic and allelic diversity in 5 groups. Among each group, KUR(AF), BANG(AS-LOC), ALxx(EU-LOC), BHwsch(AS-DE), and ARw(EU-DE) had the highest contribution to gene diversity in the order of the above grouping. Similarly, according to the allelic diversity standard (in the same order), ZIMxx, PIxx, ALxx, SHsch, and ARsch had the highest contribution. After analyzing ROH, we found a total of 144,708 fragments and 27 islands. The gene and genome regions identified by the ROH islands and QTLs indicate that chicken breeds have potential for adaptation to different production systems. Based on these findings, it is recommended to prioritize the conservation of breeds with the highest genetic diversity in each group, while paying more attention to the conservation of Asian and African breeds. Furthermore, providing a valuable reference for the conservation and utilization of chicken.

Functional genomic diversity is correlated with neutral genomic diversity in populations of an endangered rattlesnake

Theory predicts that genetic erosion in small, isolated populations of endangered species can be assessed using estimates of neutral genetic variation, yet this widely used approach has recently been questioned in the genomics era. Here, we leverage a chromosome-level genome assembly of an endangered rattlesnake (Sistrurus catenatus) combined with whole genome resequencing data (N = 110 individuals) to evaluate the relationship between levels of genome-wide neutral and functional diversity over historical and future timescales. As predicted, we found positive correlations between genome-wide estimates of neutral genetic diversity (π) and inferred levels of adaptive variation and an estimate of inbreeding mutation load, and a negative relationship between neutral diversity and an estimate of drift mutation load. However, these correlations were half as strong for projected future levels of neutral diversity based on contemporary effective population sizes. Broadly, our results confirm that estimates of neutral genetic diversity provide an accurate measure of genetic erosion in populations of a threatened vertebrate. They also provide nuance to the neutral-functional diversity controversy by suggesting that while these correlations exist, anthropogenetic impacts may have weakened these associations in the recent past and into the future.

Adaptation to the High-Arctic island environment despite long-term reduced genetic variation in Svalbard reindeer

Typically much smaller in number than their mainland counterparts, island populations are ideal systems to investigate genetic threats to small populations. The Svalbard reindeer (Rangifer tarandus platyrhynchus) is an endemic subspecies that colonized the Svalbard archipelago ca. 6,000-8,000 years ago and now shows numerous physiological and morphological adaptations to its arctic habitat. Here, we report a de-novo chromosome-level assembly for Svalbard reindeer and analyze 133 reindeer genomes spanning Svalbard and most of the species' Holarctic range, to examine the genomic consequences of long-term isolation and small population size in this insular subspecies. Empirical data, demographic reconstructions, and forward simulations show that long-term isolation and high inbreeding levels may have facilitated the reduction of highly deleterious-and to a lesser extent, moderately deleterious-variation. Our study indicates that long-term reduced genetic diversity did not preclude local adaptation to the High Arctic, suggesting that even severely bottlenecked populations can retain evolutionary potential.

Range-wide and temporal genomic analyses reveal the consequences of near-extinction in Swedish moose

Ungulate species have experienced severe declines over the past centuries through overharvesting and habitat loss. Even if many game species have recovered thanks to strict hunting regulation, the genome-wide impacts of overharvesting are still unclear. Here, we examine the temporal and geographical differences in genome-wide diversity in moose (Alces alces) over its whole range in Sweden by sequencing 87 modern and historical genomes. We found limited impact of the 1900s near-extinction event but local variation in inbreeding and load in modern populations, as well as suggestion of a risk of future reduction in genetic diversity and gene flow. Furthermore, we found candidate genes for local adaptation, and rapid temporal allele frequency shifts involving coding genes since the 1980s, possibly due to selective harvesting. Our results highlight that genomic changes potentially impacting fitness can occur over short time scales and underline the need to track both deleterious and selectively advantageous genomic variation.