Genetic load: genomic estimates and applications in non-model animals

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.

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.

Rare "Jackpot" Individuals Drive Rapid Adaptation in Threespine Stickleback

Recombination has long been considered the primary mechanism to bring beneficial alleles together, which can increase the speed of adaptation from standing genetic variation. Recombination is fundamental to the transporter hypothesis proposed to explain precise parallel adaptation in Threespine Stickleback. We studied an instance of freshwater adaptation in the Threespine Stickleback system using whole genomic data from an evolutionary time series to observe the genomic dynamics underlying rapid parallel adaptation. Our experiment showed that rapid adaptation to a freshwater environment depended on a few individuals with large haploblocks of freshwater-adaptive alleles (jackpot carriers) present among the anadromous (i.e., sea-run) founders at low frequencies. Biological kinship analyses indicate that mating among jackpot carriers and between jackpot carriers and non-jackpot individuals led to a rapid increase in freshwater-adaptive alleles within the first few generations. This process allowed the population to overcome a substantial bottleneck likely caused by the low fitness of first-generation stickleback with a few freshwater-adaptive alleles born in the lake. Additionally, we found evidence that the genetic load that emerged from population growth after the bottleneck may have been reduced through an increase in homozygosity by inbreeding, ultimately purging deleterious alleles. Recombination likely played a limited role in this case of very rapid adaptation.

Incorporating genetic load contributes to predicting Arabidopsis thaliana's response to climate change

Understanding how species respond to climate change can facilitate species conservation and crop breeding. Current prediction frameworks about population vulnerability focused on predicting range shifts or local adaptation but ignored genetic load, which is also crucial for adaptation. By analyzing 1115 globally distributed Arabidopsis thaliana natural accessions, we find that effective population size (Ne) is the major contributor of genetic load variation, both along genome and among populations, and can explain 74-94% genetic load variation in natural populations. Intriguingly, Ne affects genetic load by changing both effectiveness of purifying selection and GC biased gene conversion strength. In particular, by incorporating genetic load, genetic offset and species distribution models (SDM), we predict that, the populations at species' range edge are generally at higher risk. The populations at the eastern range perform poorer in all aspects, southern range have higher genetic offset and lower SDM suitability, while northern range have higher genetic load. Among the diverse natural populations, the Yangtze River basin population is the most vulnerable population under future climate change. Overall, here we deciphered the driving forces of genetic load in A. thaliana, and incorporated SDM, local adaptation and genetic load to predict the fate of populations under future climate change.

Time-lagged genomic erosion and future environmental risks in a bird on the brink of extinction

Global biodiversity is rapidly declining due to habitat degradation and genomic erosion, highlighting the urgent need to monitor endangered species and their genetic health. Temporal genomics and ecological modelling offer finer resolution than single-time-point measurements, providing a comprehensive view of species' recent and future trajectories. We investigated genomic erosion and environmental suitability in the critically endangered regent honeyeater (Anthochaera phrygia) by sequencing whole genomes of historical and modern specimens and building multi-temporal species distribution models (SDMs) across the last century. The species has declined from hundreds of thousands of individuals to fewer than 300 over the past 100 years. SDMs correctly predicted known patterns of local extinction in southeast Australia. Our demographic reconstructions revealed a gradual population decline from 2000 to 2500 years ago, sharply accelerating in the last 500 years due to climate variability and habitat loss. Despite this substantial demographic collapse, the regent honeyeater has lost only 9% of its genetic diversity, with no evidence of inbreeding or connectivity loss. Also, it exhibits higher diversity than many other threatened bird species. Forward-in-time genomic simulations indicate that this time lag between population decline and genetic diversity loss conceals the risk of ongoing genomic erosion into a future of rapidly degrading environmental suitability. Our work underscores the need for targeted conservation efforts and continuous genetic monitoring to prevent species extinction.

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.

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.

Mind the lag: understanding genetic extinction debt for conservation

The delay between disturbance events and genetic responses within populations is a common but surprisingly overlooked phenomenon in ecology and evolutionary and conservation genetics. If not accounted for when interpreting genetic data, this time lag problem can lead to erroneous conservation assessments. We (i) identify life-history traits related to longevity and reproductive strategies as the main determinants of time lags, (ii) evaluate potential confounding factors affecting genetic parameters during time lags, and (iii) propose approaches that allow controlling for time lags. Considering the current unprecedented rate of loss of genetic diversity and adaptive potential, we expect our novel interpretive and methodological framework for time lags to stimulate further research and discussion on the most appropriate approaches to analyse genetic diversity for conservation.

Population genomics reveals strong impacts of genetic drift without purging and guides conservation of bull and giant kelp

Kelp forests are declining in many parts of the northeast Pacific.1,2,3,4 In small populations, genetic drift can reduce adaptive variation and increase fixation of recessive deleterious alleles,5,6,7 but natural selection may purge harmful variants.8,9,10 To understand evolutionary dynamics and inform restoration strategies, we investigated genetic structure and the outcomes of genetic drift and purging by sequencing the genomes of 429 bull kelp (Nereocystis luetkeana) and 211 giant kelp (Macrocystis sp.) from the coastlines of British Columbia and Washington. We identified 6 to 7 geographically and genetically distinct clusters in each species. Low effective population size was associated with low genetic diversity and high inbreeding coefficients (including increased selfing rates), with extreme variation in these genetic health indices among bull kelp populations but more moderate variation in giant kelp. We found no evidence that natural selection is purging putative recessive deleterious alleles in either species. Instead, genetic drift has fixed many such alleles in small populations of bull kelp, leading us to predict (1) reduced within-population inbreeding depression in small populations, which may be associated with an observed shift toward increased selfing rate, and (2) hybrid vigor in crosses between small populations. Our genomic findings imply several strategies for optimal sourcing and crossing of populations for restoration and aquaculture, but these require experimental validation. Overall, our work reveals strong genetic structure and suggests that conservation strategies should consider the multiple health risks faced by small populations whose evolutionary dynamics are dominated by genetic drift.

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.

The conservation value of small population remnants: Variability in inbreeding depression and heterosis of a perennial herb, the narrow-leaved purple coneflower (Echinacea angustifolia)

Anthropogenically fragmented populations may have reduced fitness due to loss of genetic diversity and inbreeding. The extent of such fitness losses due to fragmentation and potential gains from conservation actions are infrequently assessed together empirically. Controlled crosses within and among populations can identify whether populations are at risk of inbreeding depression and whether inter-population crossing alleviates fitness loss. Because fitness depends on the environment and life stage, studies quantifying cumulative fitness over a large portion of the lifecycle in conditions that mimic natural environments are most informative. To assess the fitness consequences of habitat fragmentation, we leveraged controlled within-family, within-population, and between-population crosses to quantify inbreeding depression and heterosis in seven populations of Echinacea angustifolia within a 6,400-ha area. We then assessed cumulative offspring fitness after 14 yr of growth in a natural experimental plot (N = 1,136). The mean fitness of progeny from within-population crosses varied considerably, indicating genetic differentiation among source populations, even though these sites are all less than 9 km apart. The fitness consequences of within-family and between-population crosses varied in magnitude and direction. Only one of the seven populations showed inbreeding depression of high effect, while four populations showed substantial heterosis. Outbreeding depression was rare and slight. Our findings indicate that local crossings between isolated populations yield unpredictable fitness consequences ranging from slight decreases to substantial increases. Interestingly, inbreeding depression and heterosis did not relate closely to population size, suggesting that all fragmented populations could contribute to conservation goals as either pollen recipients or donors.

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.

What Can Genome Sequence Data Reveal About Population Viability?

Biologists have long sought to understand the impacts of deleterious genetic variation on fitness and population viability. However, our understanding of these effects in the wild is incomplete, in part due to the rarity of sufficient genetic and demographic data needed to measure their impact. The genomics revolution is promising a potential solution by predicting the effects of deleterious genetic variants (genetic load) bioinformatically from genome sequences alone bypassing the need for costly demographic data. After a historical perspective on the theoretical and empirical basis of our understanding of the dynamics and fitness effects of deleterious genetic variation, we evaluate the potential for these new genomic measures of genetic load to predict population viability. We argue that current genomic analyses alone cannot reliably predict the effects of deleterious genetic variation on population growth, because these depend on demographic, ecological and genetic parameters that need more than just genome sequence data to be measured. Thus, while purely genomic analyses of genetic load promise to improve our understanding of the composition of the genetic load, they are currently of little use for evaluating population viability. Demographic data and ecological context remain crucial to our understanding of the consequences of deleterious genetic variation for population fitness. However, when combined with such demographic and ecological data, genomic information can offer important insights into genetic variation and inbreeding that are crucial for conservation decision making.

Genomic and fitness consequences of a near-extinction event in the northern elephant seal

Understanding the genetic and fitness consequences of anthropogenic bottlenecks is crucial for biodiversity conservation. However, studies of bottlenecked populations combining genomic approaches with fitness data are rare. Theory predicts that severe bottlenecks deplete genetic diversity, exacerbate inbreeding depression and decrease population viability. However, actual outcomes are complex and depend on how a species' unique demography affects its genetic load. We used population genetic and veterinary pathology data, demographic modelling, whole-genome resequencing and forward genetic simulations to investigate the genomic and fitness consequences of a near-extinction event in the northern elephant seal. We found no evidence of inbreeding depression within the contemporary population for key fitness components, including body mass, blubber thickness and susceptibility to parasites and disease. However, we detected a genomic signature of a recent extreme bottleneck (effective population size = 6; 95% confidence interval = 5.0-7.5) that will have purged much of the genetic load, potentially leading to the lack of observed inbreeding depression in our study. Our results further suggest that deleterious genetic variation strongly impacted the post-bottleneck population dynamics of the northern elephant seal. Our study provides comprehensive empirical insights into the intricate dynamics underlying species-specific responses to anthropogenic bottlenecks.

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.

Fitness consequences of structural variation inferred from a House Finch pangenome

Genomic structural variants (SVs) play a crucial role in adaptive evolution, yet their average fitness effects and characterization with pangenome tools are understudied in wild animal populations. We constructed a pangenome for House Finches (Haemorhous mexicanus), a model for studies of host-pathogen coevolution, using long-read sequence data on 16 individuals (32 de novo-assembled haplotypes) and one outgroup. We identified 887,118 SVs larger than 50 base pairs, mostly (60%) involving repetitive elements, with reduced SV diversity in the eastern US as a result of its introduction by humans. The distribution of fitness effects of genome-wide SVs was estimated using maximum likelihood approaches and revealed that SVs in both coding and noncoding regions were on average more deleterious than smaller indels or single nucleotide polymorphisms. The reference-free pangenome facilitated identification of a > 10-My-old, 11-megabase-long pericentric inversion on chromosome 1. We found that the genotype frequencies of the inversion, estimated from 135 birds widely sampled temporally and geographically, increased steadily over the 25 y since House Finches were first exposed to the bacterial pathogen Mycoplasma gallisepticum and showed signatures of balancing selection, capturing genes related to immunity and telomerase activity. We also observed shorter telomeres in populations with a greater number of years exposure to Mycoplasma. Our study illustrates the utility of long-read sequencing and pangenome methods for understanding wild animal populations, estimating fitness effects of genome-wide SVs, and advancing our understanding of adaptive evolution through structural variation.

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.

Extensive Population Structure Highlights an Apparent Paradox of Stasis in the Impala (Aepyceros melampus)

Impalas are unusual among bovids because they have remained morphologically similar over millions of years-a phenomenon referred to as evolutionary stasis. Here, we sequenced 119 whole genomes from the two extant subspecies of impala, the common (Aepyceros melampus melampus) and black-faced (A. m. petersi) impala. We investigated the evolutionary forces working within the species to explore how they might be associated with its evolutionary stasis as a taxon. Despite being one of the most abundant bovid species, we found low genetic diversity overall, and a phylogeographic signal of spatial expansion from southern to eastern Africa. Contrary to expectations under a scenario of evolutionary stasis, we found pronounced genetic structure between and within the two subspecies with indications of ancient, but not recent, gene flow. Black-faced impala and eastern African common impala populations had more runs of homozygosity than common impala in southern Africa, and, using a proxy for genetic load, we found that natural selection is working less efficiently in these populations compared to the southern African populations. Together with the fossil record, our results are consistent with a fixed-optimum model of evolutionary stasis, in which impalas in the southern African core of the range are able to stay near their evolutionary fitness optimum as a generalist ecotone species, whereas eastern African impalas may struggle to do so due to the effects of genetic drift and reduced adaptation to the local habitat, leading to recurrent local extinction in eastern Africa and re-colonisation from the South.

Sonification of Genomic Data to Represent Genetic Load in Zoo Populations

Maintaining a diverse gene pool is important in the captive management of zoo populations, especially in endangered species such as the pink pigeon (Nesoenas mayeri). However, due to the limited number of breeding individuals and relaxed natural selection, the loss of variation and accumulation of harmful variants is inevitable. Inbreeding results in a loss of fitness (i.e., inbreeding depression), principally because related parents are more likely to transmit a copy of the same recessive deleterious genetic variant to their offspring. Genomics-informed captive breeding can manage harmful variants by artificial selection, reducing the genetic load by avoiding the inheritance of two copies of the same harmful variant. To explain this concept in an interactive way to zoo visitors, we developed a sonification game to represent the fitness impacts of harmful variants by detuning notes in a familiar musical melody (i.e., Beethoven's Für Elise). Conceptually, zoo visitors play a game aiming to create the most optimal pink pigeon offspring in terms of inbreeding depression. They select virtual crosses between pink pigeon individuals and listen for the detuning of the melody, which represents the realised load of the resultant offspring. Here we present the sonification algorithm and the results of an online survey to see whether participants could identify the most and least optimal offspring from three potential pink pigeon offspring. Of our 98 respondents, 85 (86.7%) correctly identified the least optimal offspring, 73 (74.5%) correctly identified the most optimal, and 62 (63.3%) identified both the most and least optimal offspring using only the sonification.

Lineage Differentiation and Genomic Vulnerability in a Relict Tree From Subtropical Forests

The subtropical forests of East Asia are renowned for their high plant diversity, particularly the abundance of ancient relict species. However, both the evolutionary history of these relict species and their capacity for resilience in the face of impending climatic changes remain unclear. Using whole-genome resequencing data, we investigated the lineage differentiation and demographic history of the relict and endangered tree, Bretschneidera sinensis (Akaniaceae). We employed a combination of population genomic and landscape genomic approaches to evaluate variation in mutation load and genomic offset, aiming to predict how different populations may respond to climate change. Our analysis revealed a profound genomic divergence between the East and West lineages, likely as the result of recurrent bottlenecks due to climatic fluctuations during the glacial period. Furthermore, we identified several genes potentially linked to growth characteristics and hypoxia response that had been subjected to positive selection during the lineage differentiation. Our assessment of genomic vulnerability uncovered a significantly higher mutation load and genomic offset in the edge populations of B. sinensis compared to their core counterparts. This implies that the edge populations are likely to experience the most significant impact from the predicted climate conditions. Overall, our research sheds light on the historical lineage differentiation and contemporary genomic vulnerability of B. sinensis. Broadening our understanding of the speciation history and future resilience of relict and endangered species such as B. sinensis, is crucial in developing effective conservation strategies in anticipation of future climatic changes.

Integrating evolutionary genomics of forest trees to inform future tree breeding amid rapid climate change

Global climate change is leading to rapid and drastic shifts in environmental conditions, posing threats to biodiversity and nearly all life forms worldwide. Forest trees serve as foundational components of terrestrial ecosystems and play a crucial and leading role in combating and mitigating the adverse effects of extreme climate events, despite their own vulnerability to these threats. Therefore, understanding and monitoring how natural forests respond to rapid climate change is a key priority for biodiversity conservation. Recent progress in evolutionary genomics, driven primarily by cutting-edge multi-omics technologies, offers powerful new tools to address several key issues. These include precise delineation of species and evolutionary units, inference of past evolutionary histories and demographic fluctuations, identification of environmentally adaptive variants, and measurement of genetic load levels. As the urgency to deal with more extreme environmental stresses grows, understanding the genomics of evolutionary history, local adaptation, future responses to climate change, and conservation and restoration of natural forest trees will be critical for research at the nexus of global change, population genomics, and conservation biology. In this review, we explore the application of evolutionary genomics to assess the effects of global climate change using multi-omics approaches and discuss the outlook for breeding of climate-adapted trees.

The relationship between neutral genetic diversity and performance in wild arthropod populations

Larger effective populations (Ne) are characterized by higher genetic diversity, which is expected to predict population performance (average individual performance that influences fitness). Empirical studies of the relationship between neutral diversity and performance mostly represent species with small Ne, while there is limited data from the species-rich and ecologically important arthropods that are assumed to have large Ne but are threatened by massive declines. We performed a systematic literature search and used meta-analytical models to test the prediction of a positive association between neutral genetic diversity and performance in wild arthropods. From 14 relevant studies of 286 populations, we detected a weak (r = 0.15) but nonsignificant positive association both in the full data set (121 effect sizes) and a reduced data set accounting for dependency (14 effect sizes). Theory predicts that traits closely associated with fitness show a relatively stronger correlation with neutral diversity; this relationship was upheld for longevity and marginally for reproduction. Our analyses point to major knowledge gaps in our understanding of relationships between neutral diversity and performance. Future studies using genome-wide data sets across populations could guide more powerful designs to evaluate relationships between adaptive, deleterious and neutral diversity and performance.

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.

Leveraging genomic load estimates to optimize captive breeding programmes

Rapid biodiversity loss threatens many species with extinction. Captive populations of species of conservation concern (such as those housed in zoos and dedicated breeding centres) act as an insurance should wild populations go extinct or need supplemental individuals to boost populations. Limited resources mean that captive populations are almost always small and started from few founding individuals. As a result, captive populations require careful management to minimize negative genetic impacts, with decisions about which individuals to breed together often guided by the principle of minimizing relatedness. Typically this strategy aims to retain 90% of genetic diversity over 200 years (Soulé et al., Zoo Biology, 1986, 5, 101), but it has a weakness in that it does not directly manage for genetic load. In this issue of Molecular Ecology Resources, Speak et al. (Molecular Ecology Resources, 2024, e13967) present a novel proof-of-concept study for taking this next step and incorporating estimates of individual genetic load into the planning of captive breeding, using an approach that is likely to be widely applicable to many captive populations.

Relaxed selection in evolution of genes regulating limb development gives clue to variation in forelimb morphology of cetaceans and other mammals

Cetaceans have evolved unique limb structures, such as flippers, due to genetic changes during their transition to aquatic life. However, the full understanding of the genetic and evolutionary mechanisms behind these changes is still developing. By examining 25 limb-related protein-coding genes across various mammalian species, we compared genetic changes between aquatic mammals, like whales, and other mammals with unique limb structures such as bats, rodents and elephants. Our findings revealed significant modifications in limb-related genes, including variations in the Hox, GDF5 and Evx genes. Notably, a relaxed selection in several key genes was observed, suggesting a lifting of developmental constraints, which might have facilitated the emergence of morphological innovations in cetacean limb morphology. We also uncovered non-synonymous changes, insertions and deletions in these genes, particularly in the polyalanine tract of HOXD13, which are distinctive to cetaceans or convergent with other aquatic mammals. These genetic variations correlated with the diverse and specialized limb structures observed in cetaceans, indicating a complex interplay of relaxed selection and specific mutations in mammalian limb evolution.

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.

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.

Large-scale genome sequencing of giant pandas improves the understanding of population structure and future conservation initiatives

The extinction risk of the giant panda has been demoted from "endangered" to "vulnerable" on the International Union for Conservation of Nature Red List, but its habitat is more fragmented than ever before, resulting in 33 isolated giant panda populations according to the fourth national survey released by the Chinese government. Further comprehensive investigations of the genetic background and in-depth assessments of the conservation status of wild populations are still necessary and urgently needed. Here, we sequenced the genomes of 612 giant pandas with an average depth of ~26× and generated a high-resolution map of genomic variation with more than 20 million variants covering wild individuals from six mountain ranges and captive representatives in China. We identified distinct genetic clusters within the Minshan population by performing a fine-grained genetic structure. The estimation of inbreeding and genetic load associated with historical population dynamics suggested that future conservation efforts should pay special attention to the Qinling and Liangshan populations. Releasing captive individuals with a genetic background similar to the recipient population appears to be an advantageous genetic rescue strategy for recovering the wild giant panda populations, as this approach introduces fewer deleterious mutations into the wild population than mating with differentiated lineages. These findings emphasize the superiority of large-scale population genomics to provide precise guidelines for future conservation of the giant panda.

Contrasting demographic history and mutational load in three threatened whitebark pines (Pinus subsect. Gerardianae): implications for conservation

Demographic history and mutational load are of paramount importance for the adaptation of the endangered species. However, the effects of population evolutionary history and genetic load on the adaptive potential in endangered conifers remain unclear. Here, using population transcriptome sequencing, whole chloroplast genomes and mitochondrial DNA markers, combined with niche analysis, we determined the demographic history and mutational load for three threatened whitebark pines having different endangered statuses, Pinus bungeana, P. gerardiana and P. squamata. Demographic inference indicated that severe bottlenecks occurred in all three pines at different times, coinciding with periods of major climate and geological changes; in contrast, while P. bungeana experienced a recent population expansion, P. gerardiana and P. squamata maintained small population sizes after bottlenecking. Abundant homozygous-derived variants accumulated in the three pines, particularly in P. squamata, while the species with most heterozygous variants was P. gerardiana. Abundant moderately and few highly deleterious variants accumulated in the pine species that have experienced the most severe demographic bottlenecks (P. gerardiana and P. squamata), most likely because of purging effects. Finally, niche modeling showed that the distribution of P. bungeana might experience a significant expansion in the future, and the species' identified genetic clusters are also supported by differences in the ecological niche. The integration of genomic, demographic and niche data has allowed us to prove that the three threatened pines have contrasting patterns of demographic history and mutational load, which may have important implications in their adaptive potential and thus are also key for informing conservation planning.

Study of the cyto- and genotoxic activity of water from the Kapshagai reservoir (Kazakhstan) on laboratory mice

Chemical compounds in the environment, which exhibit toxic and genotoxic activity, increase the mutational pressure on biota. This study aimed to investigate the genotoxic, mutagenic, and toxic effects of water from the Ile River and the Kapshagai Reservoir, both sites of active economic activities. Cytogenetic analysis of bone marrow from mice exposed to water samples from the Ile River and the Kapshagai Reservoir revealed a statistically significant increase in aberrant (p<0.05) and polyploid cells (p<0.01), as well as a decrease in the mitotic index (p<0.001), compared to the negative control. The water samples caused statistically significant increases in single- and double-strand DNA breaks in cells across various organs in the experimental mice compared to unexposed animals (p<0.001). These observations suggest the existence of chemical compounds within the water samples from the Kapshagai Reservoir and the Ile River, which exhibit genotoxic, mutagenic, and toxic properties.

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.

Temporal dynamics of woolly mammoth genome erosion prior to extinction

A number of species have recently recovered from near-extinction. Although these species have avoided the immediate extinction threat, their long-term viability remains precarious due to the potential genetic consequences of population declines, which are poorly understood on a timescale beyond a few generations. Woolly mammoths (Mammuthus primigenius) became isolated on Wrangel Island around 10,000 years ago and persisted for over 200 generations before becoming extinct around 4,000 years ago. To study the evolutionary processes leading up to the mammoths' extinction, we analyzed 21 Siberian woolly mammoth genomes. Our results show that the population recovered quickly from a severe bottleneck and remained demographically stable during the ensuing six millennia. We find that mildly deleterious mutations gradually accumulated, whereas highly deleterious mutations were purged, suggesting ongoing inbreeding depression that lasted for hundreds of generations. The time-lag between demographic and genetic recovery has wide-ranging implications for conservation management of recently bottlenecked populations.

Chromosome-Level Reference Genome of the Ponza Grayling (Hipparchia sbordonii), an Italian Endemic and Endangered Butterfly

Islands are crucial evolutionary hotspots, providing unique opportunities for differentiation of novel biodiversity and long-term segregation of endemic species. Islands are also fragile ecosystems, where biodiversity is more exposed to environmental and anthropogenic pressures than on continents. The Ponza grayling, Hipparchia sbordonii, is an endemic butterfly species that is currently found only in two tiny islands of the Pontine archipelago, off the coast of Italy, occupying an area smaller than 10 km2. It has been classified as Endangered (IUCN) because of the extremely limited area of occurrence, population fragmentation, and the recent demographic decline. Thanks to a combination of different assemblers of long and short genomic reads, bulk transcriptome RNAseq, and synteny analysis with phylogenetically close butterflies, we produced a highly contiguous, chromosome-scale annotated reference genome for the Ponza grayling, including 28 autosomes and the Z sexual chromosomes. The final assembly spanned 388.61 Gb with a contig N50 of 14.5 Mb and a BUSCO completeness score of 98.5%. Synteny analysis using four other butterfly species revealed high collinearity with Hipparchia semele and highlighted 10 intrachromosomal inversions longer than 10 kb, of which two appeared on the lineage leading to H. sbordonii. Our results show that a chromosome-scale reference genome is attainable also when chromatin conformation data may be impractical or present specific technical challenges. The high-quality genomic resource for H. sbordonii opens up new opportunities for the accurate assessment of genetic diversity and genetic load and for the investigations of the genomic novelties characterizing the evolutionary path of this endemic island species.

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.

Efficient purging of deleterious mutations contributes to the survival of a rare conifer

Cupressaceae is a conifer family rich in plants of horticultural importance, including Cupressus, Chamaecyparis, Juniperus, and Thuja, yet genomic surveys are lacking for this family. Cupressus gigantea, one of the many rare conifers that are threatened by climate change and anthropogenic habitat fragmentation, plays an ever-increasing role in ecotourism in Tibet. To infer how past climate change has shaped the population evolution of this species, we generated a de novo chromosome-scale genome (10.92 Gb) and compared the species' population history and genetic load with that of a widespread close relative, C. duclouxiana. Our demographic analyses, based on 83 re-sequenced individuals from multiple populations of the two species, revealed a sharp decline of population sizes during the first part of the Quaternary. However, populations of C. duclouxiana then started to recover, while C. gigantea populations continued to decrease until recently. The total genomic diversity of C. gigantea is smaller than that of C. duclouxiana, but contrary to expectations, C. gigantea has fewer highly and mildly deleterious mutations than C. duclouxiana, and simulations and statistical tests support purifying selection during prolonged inbreeding as the explanation. Our results highlight the evolutionary consequences of decreased population size on the genetic burden of a long-lived endangered conifer with large genome size and suggest that genetic purging deserves more attention in conservation management.

Global and Local Ancestry and its Importance: A Review

The fastest way to significantly change the composition of a population is through admixture, an evolutionary mechanism. In animal breeding history, genetic admixture has provided both short-term and long-term advantages by utilizing the phenomenon of complementarity and heterosis in several traits and genetic diversity, respectively. The traditional method of admixture analysis by pedigree records has now been replaced greatly by genome-wide marker data that enables more precise estimations. Among these markers, SNPs have been the popular choice since they are cost-effective, not so laborious, and automation of genotyping is easy. Certain markers can suggest the possibility of a population's origin from a sample of DNA where the source individual is unknown or unwilling to disclose their lineage, which are called Ancestry-Informative Markers (AIMs). Revealing admixture level at the locus-specific level is termed as local ancestry and can be exploited to identify signs of recent selective response and can account for genetic drift. Considering the importance of genetic admixture and local ancestry, in this mini-review, both concepts are illustrated, encompassing basics, their estimation/identification methods, tools/software used and their applications.

Population Genomic History of the Endangered Anatolian and Cyprian Mouflons in Relation to Worldwide Wild, Feral, and Domestic Sheep Lineages

Once widespread in their homelands, the Anatolian mouflon (Ovis gmelini anatolica) and the Cyprian mouflon (Ovis gmelini ophion) were driven to near extinction during the 20th century and are currently listed as endangered populations by the International Union for Conservation of Nature. While the exact origins of these lineages remain unclear, they have been suggested to be close relatives of domestic sheep or remnants of proto-domestic sheep. Here, we study whole genome sequences of n = 5 Anatolian mouflons and n = 10 Cyprian mouflons in terms of population history and diversity, comparing them with eight other extant sheep lineages. We find reciprocal genetic affinity between Anatolian and Cyprian mouflons and domestic sheep, higher than all other studied wild sheep genomes, including the Iranian mouflon (O. gmelini). Studying diversity indices, we detect a considerable load of short runs of homozygosity blocks (<2 Mb) in both Anatolian and Cyprian mouflons, reflecting small effective population size (Ne). Meanwhile, Ne and mutation load estimates are lower in Cyprian compared with Anatolian mouflons, suggesting the purging of recessive deleterious variants in Cyprian sheep under a small long-term Ne, possibly attributable to founder effects, island isolation, introgression from domestic lineages, or differences in their bottleneck dynamics. Expanding our analyses to worldwide wild and feral Ovis genomes, we observe varying viability metrics among different lineages and a limited consistency between viability metrics and International Union for Conservation of Nature conservation status. Factors such as recent inbreeding, introgression, and unique population dynamics may have contributed to the observed disparities.

Can genetic rescue help save Arabia's last big cat?

Genetic diversity underpins evolutionary potential that is essential for the long-term viability of wildlife populations. Captive populations harbor genetic diversity potentially lost in the wild, which could be valuable for release programs and genetic rescue. The Critically Endangered Arabian leopard (Panthera pardus nimr) has disappeared from most of its former range across the Arabian Peninsula, with fewer than 120 individuals left in the wild, and an additional 64 leopards in captivity. We (i) examine genetic diversity in the wild and captive populations to identify global patterns of genetic diversity and structure; (ii) estimate the size of the remaining leopard population across the Dhofar mountains of Oman using spatially explicit capture-recapture models on DNA and camera trap data, and (iii) explore the impact of genetic rescue using three complementary computer modeling approaches. We estimated a population size of 51 (95% CI 32-79) in the Dhofar mountains and found that 8 out of 25 microsatellite alleles present in eight loci in captive leopards were undetected in the wild. This includes two alleles present only in captive founders known to have been wild-sourced from Yemen, which suggests that this captive population represents an important source for genetic rescue. We then assessed the benefits of reintroducing novel genetic diversity into the wild population as well as the risks of elevating the genetic load through the release of captive-bred individuals. Simulations indicate that genetic rescue can improve the long-term viability of the wild population by reducing its genetic load and realized load. The model also suggests that the genetic load has been partly purged in the captive population, potentially making it a valuable source population for genetic rescue. However, the greater loss of its genetic diversity could exacerbate genomic erosion of the wild population during a rescue program, and these risks and benefits should be carefully evaluated. An important next step in the recovery of the Arabian leopard is to empirically validate these conclusions, implement and monitor a genomics-informed management plan, and optimize a strategy for genetic rescue as a tool to recover Arabia's last big cat.

Pervasive relaxed selection in termite genomes

Genetic changes that enabled the evolution of eusociality have long captivated biologists. More recently, attention has focussed on the consequences of eusociality on genome evolution. Studies have reported higher molecular evolutionary rates in eusocial hymenopteran insects compared with their solitary relatives. To investigate the genomic consequences of eusociality in termites, we analysed nine genomes, including newly sequenced genomes from three non-eusocial cockroaches. Using a phylogenomic approach, we found that termite genomes have experienced lower rates of synonymous substitutions than those of cockroaches, possibly as a result of longer generation times. We identified higher rates of non-synonymous substitutions in termite genomes than in cockroach genomes, and identified pervasive relaxed selection in the former (24-31% of the genes analysed) compared with the latter (2-4%). We infer that this is due to reductions in effective population size, rather than gene-specific effects (e.g. indirect selection of caste-biased genes). We found no obvious signature of increased genetic load in termites, and postulate efficient purging of deleterious alleles at the colony level. Additionally, we identified genomic adaptations that may underpin caste differentiation, such as genes involved in post-translational modifications. Our results provide insights into the evolution of termites and the genomic consequences of eusociality more broadly.

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 transcriptome assembly of the European freshwater mussel Unio elongatulus C. Pfeiffer, 1825

Freshwater mussels of the order Unionida are a global conservation concern. Species of this group are strictly freshwater, sessile, slow-growing animals and, extremely sensitive to environmental changes. Human-mediated changes in freshwater habitats are imposing enormous pressure on the survival of freshwater mussels. Although a few flagship species are protected in Europe, other highly imperilled species receive much less attention. Moreover, knowledge about biology, ecology, and evolution and proper conservation assessments of many European species are still sparse. This knowledge gap is further aggravated by the lack of genomic resources available, which are key tools for conservation. Here we present the transcriptome assembly of Unio elongatulus C. Pfeiffer, 1825, one of the least studied European freshwater mussels. Using the individual sequencing outputs from eight physiologically representative mussel tissues, we provide an annotated panel of tissue-specific Relative Gene Expression profiles. These resources are pivotal to studying the species' biological and ecological features, as well as helping to understand its vulnerability to current and future threats.

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

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

Forces driving transposable element load variation during Arabidopsis range expansion

Genetic load refers to the accumulated and potentially life-threatening deleterious mutations in populations. Understanding the mechanisms underlying genetic load variation of transposable element (TE) insertion, a major large-effect mutation, during range expansion is an intriguing question in biology. Here, we used 1,115 global natural accessions of Arabidopsis (Arabidopsis thaliana) to study the driving forces of TE load variation during its range expansion. TE load increased with range expansion, especially in the recently established Yangtze River basin population. Effective population size, which explains 62.0% of the variance in TE load, high transposition rate, and selective sweeps contributed to TE accumulation in the expanded populations. We genetically mapped and identified multiple candidate causal genes and TEs, and revealed the genetic architecture of TE load variation. Overall, this study reveals the variation in TE genetic load during Arabidopsis expansion and highlights the causes of TE load variation from the perspectives of both population genetics and quantitative genetics.

High genetic load without purging in caribou, a diverse species at risk

High intra-specific genetic diversity is associated with adaptive potential, which is key for resilience to global change. However, high variation may also support deleterious alleles through genetic load, thereby increasing the risk of inbreeding depression if population sizes decrease. Purging of deleterious variation has been demonstrated in some threatened species. However, less is known about the costs of declines and inbreeding in species with large population sizes and high genetic diversity even though this encompasses many species globally that are expected to undergo population declines. Caribou is a species of ecological and cultural significance in North America with a wide distribution supporting extensive phenotypic variation but with some populations undergoing significant declines resulting in their at-risk status in Canada. We assessed intra-specific genetic variation, adaptive divergence, inbreeding, and genetic load across populations with different demographic histories using an annotated chromosome-scale reference genome and 66 whole-genome sequences. We found high genetic diversity and nine phylogenomic lineages across the continent with adaptive diversification of genes, but also high genetic load among lineages. We found highly divergent levels of inbreeding across individuals, including the loss of alleles by drift but not increased purging in inbred individuals, which had more homozygous deleterious alleles. We also found comparable frequencies of homozygous deleterious alleles between lineages regardless of nucleotide diversity. Thus, further inbreeding may need to be mitigated through conservation efforts. Our results highlight the "double-edged sword" of genetic diversity that may be representative of other species atrisk affected by anthropogenic activities.

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.

How Can Genomics Help or Hinder Wildlife Conservation?

Genomic data are becoming increasingly affordable and easy to collect, and new tools for their analysis are appearing rapidly. Conservation biologists are interested in using this information to assist in management and planning but are typically limited financially and by the lack of genomic resources available for non-model taxa. It is therefore important to be aware of the pitfalls as well as the benefits of applying genomic approaches. Here, we highlight recent methods aimed at standardizing population assessments of genetic variation, inbreeding, and forms of genetic load and methods that help identify past and ongoing patterns of genetic interchange between populations, including those subjected to recent disturbance. We emphasize challenges in applying some of these methods and the need for adequate bioinformatic support. We also consider the promises and challenges of applying genomic approaches to understand adaptive changes in natural populations to predict their future adaptive capacity.

Moderating the neutralist-selectionist debate: exactly which propositions are we debating, and which arguments are valid?

Half a century after its foundation, the neutral theory of molecular evolution continues to attract controversy. The debate has been hampered by the coexistence of different interpretations of the core proposition of the neutral theory, the 'neutral mutation-random drift' hypothesis. In this review, we trace the origins of these ambiguities and suggest potential solutions. We highlight the difference between the original, the revised and the nearly neutral hypothesis, and re-emphasise that none of them equates to the null hypothesis of strict neutrality. We distinguish the neutral hypothesis of protein evolution, the main focus of the ongoing debate, from the neutral hypotheses of genomic and functional DNA evolution, which for many species are generally accepted. We advocate a further distinction between a narrow and an extended neutral hypothesis (of which the latter posits that random non-conservative amino acid substitutions can cause non-ecological phenotypic divergence), and we discuss the implications for evolutionary biology beyond the domain of molecular evolution. We furthermore point out that the debate has widened from its initial focus on point mutations, and also concerns the fitness effects of large-scale mutations, which can alter the dosage of genes and regulatory sequences. We evaluate the validity of neutralist and selectionist arguments and find that the tested predictions, apart from being sensitive to violation of underlying assumptions, are often derived from the null hypothesis of strict neutrality, or equally consistent with the opposing selectionist hypothesis, except when assuming molecular panselectionism. Our review aims to facilitate a constructive neutralist-selectionist debate, and thereby to contribute to answering a key question of evolutionary biology: what proportions of amino acid and nucleotide substitutions and polymorphisms are adaptive?