Inbreeding increases genome homozygosity within populations, which can exacerbate inbreeding depression by exposing homozygous deleterious alleles that are responsible for declines in fitness traits. In small populations, genetic purging that occurs under pressure of natural selection acts as an opposing force, contributing to a reduction of deleterious alleles. Both inbreeding and genetic purging are paramount in the field of conservation genomics. The Amur tiger (Panthera tigris altaica) lives in small populations in the forests of Northeast Asia and is among the most endangered animals on the planet. Using genome-wide assessment and comparison, we reveal substantially higher and more extensive inbreeding in wild Amur tigers (FROH = 0.50) than in captive individuals (FROH = 0.24). However, a relatively reduced number of loss-of-function mutations in wild Amur tigers is observed when compared with captive individuals, indicating a genetic purging of relatively large-effect inbreeding load. The higher ratio of homozygous mutation load and number of fixed damaging alleles in the wild population indicate a less-efficient genetic purging, with purifying selection also contributing to this process. These findings provide valuable insights for future conservation of Amur tigers.

Inbreeding depression has been documented in various fitness traits in a wide range of species and taxa, however, the mutational basis is not yet well understood. We investigate how putatively deleterious variation influences fitness and is shaped by individual ancestry by re-sequencing complete genomes of 37 individuals in a natural arctic fox (Vulpes lagopus) population subjected to both inbreeding depression and genetic rescue. We find that individuals with high proportion of homozygous loss of function genotypes (LoFs), which are predicted to exert a strong effect on fitness, generally have lower lifetime reproductive success and live shorter lives compared with individuals with lower proportion of LoFs. We also find that juvenile survival is negatively associated with the proportion of homozygous missense genotypes and positively associated with genome wide heterozygosity. Our results demonstrate that homozygosity of strongly and moderately deleterious mutations can be an important cause of trait specific inbreeding depression in wild populations, and mark an important step towards making more informed decisions using applied conservation genetics.

Inbreeding depression has been documented in various fitness traits in a wide range of species and taxa, however, the mutational basis is not yet well understood. We investigate how putatively deleterious variation influences fitness and is shaped by individual ancestry by re-sequencing complete genomes of 37 individuals in a natural arctic fox (Vulpes lagopus) population subjected to both inbreeding depression and genetic rescue. We find that individuals with high proportion of homozygous loss of function genotypes (LoFs), which are predicted to exert a strong effect on fitness, generally have lower lifetime reproductive success and live shorter lives compared with individuals with lower proportion of LoFs. We also find that juvenile survival is negatively associated with the proportion of homozygous missense genotypes and positively associated with genome wide heterozygosity. Our results demonstrate that homozygosity of strongly and moderately deleterious mutations can be an important cause of trait specific inbreeding depression in wild populations, and mark an important step towards making more informed decisions using applied conservation genetics.

Critically endangered species are usually restricted to small and isolated populations. High inbreeding without gene flow among populations further aggravates their threatened condition and reduces the likelihood of their long-term survival. Chinese alligator (Alligator sinensis) is one of the most endangered crocodiles in the world and has experienced a continuous decline over the past c. 1 million years. In order to identify the genetic status of the remaining populations and aid conservation efforts, we assembled the first high-quality chromosome-level genome of Chinese alligator and explored the genomic characteristics of three extant breeding populations. Our analyses revealed the existence of at least three genetically distinct populations, comprising two breeding populations in China (Changxing and Xuancheng) and one breeding population in an American wildlife refuge. The American population does not belong to the last two populations of its native range (Xuancheng and Changxing), thus representing genetic diversity extinct in the wild and provides future opportunities for genetic rescue. Moreover, the effective population size of these three populations has been continuously declining over the past 20 ka. Consistent with this decline, the species shows extremely low genetic diversity, a large proportion of long runs of homozygous fragments, and mutational load across the genome. Finally, to provide genomic insights for future breeding management and conservation, we assessed the feasibility of mixing extant populations based on the likelihood of introducing new deleterious alleles and signatures of local adaptation. Overall, this study provides a valuable genomic resource and important genomic insights into the ecology, evolution, and conservation of critically endangered alligators. 

Understanding how climatic and environmental changes, as well as human activities, induce changes in the distribution and population size of avian species refines our ability to predict future impacts on threatened species. Using multilocus genetic data, we show that the population of a threatened New Zealand endemic open-habitat specialist, the Black-fronted Tern Chlidonias albostriatus – in contrast to forest specialists – expanded during the last glacial period. The population has decreased subsequently despite the availability of extensive open habitat after human arrival to New Zealand. We conclude that population changes for open habitat specialists such as Black-fronted Terns in pre-human New Zealand were habitat-dependent, similar to Northern Hemisphere cold-adapted species, whereas post-human settlement populations were constrained by predators independent of habitat availability, similar to other island endemic species. 

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.

Background: Many wild species have suffered drastic population size declines over the past centuries, which have led to 'genomic erosion' processes characterized by reduced genetic diversity, increased inbreeding, and accumulation of harmful mutations. Yet, genomic erosion estimates of modern-day populations often lack concordance with dwindling population sizes and conservation status of threatened species. One way to directly quantify the genomic consequences of population declines is to compare genome-wide data from pre-decline museum samples and modern samples. However, doing so requires computational data processing and analysis tools specifically adapted to comparative analyses of degraded, ancient or historical, DNA data with modern DNA data as well as personnel trained to perform such analyses.

Results: Here, we present a highly flexible, scalable, and modular pipeline to compare patterns of genomic erosion using samples from disparate time periods. The GenErode pipeline uses state-of-the-art bioinformatics tools to simultaneously process whole-genome re-sequencing data from ancient/historical and modern samples, and to produce comparable estimates of several genomic erosion indices. No programming knowledge is required to run the pipeline and all bioinformatic steps are well-documented, making the pipeline accessible to users with different backgrounds. GenErode is written in Snakemake and Python3 and uses Conda and Singularity containers to achieve reproducibility on high-performance compute clusters. The source code is freely available on GitHub (https://github.com/NBISweden/GenErode).

Conclusions: GenErode is a user-friendly and reproducible pipeline that enables the standardization of genomic erosion indices from temporally sampled whole genome re-sequencing data.

Accelerating climate change is causing severe habitat fragmentation in the Arctic, threatening the persistence of many cold-adapted species. The Scandinavian arctic fox (Vulpes lagopus) is highly fragmented, with a once continuous, circumpolar distribution, it struggled to recover from a demographic bottleneck in the late 19th century. The future persistence of the entire Scandinavian population is highly dependent on the northernmost Fennoscandian subpopulations (Scandinavia and the Kola Peninsula), to provide a link to the viable Siberian population. By analyzing 43 arctic fox genomes, we quantified genomic variation and inbreeding in these populations. Signatures of genome erosion increased from Siberia to northern Sweden indicating a stepping-stone model of connectivity. In northern Fennoscandia, runs of homozygosity (ROH) were on average ~1.47-fold longer than ROH found in Siberia, stretching almost entire scaffolds. Moreover, consistent with recent inbreeding, northern Fennoscandia harbored more homozygous deleterious mutations, whereas Siberia had more in heterozygous state. This study underlines the value of documenting genome erosion following population fragmentation to identify areas requiring conservation priority. With the increasing fragmentation and isolation of Arctic habitats due to global warming, understanding the genomic and demographic consequences is vital for maintaining evolutionary potential and preventing local extinctions.