Inbreeding load and purging: implications for the short-term survival and the conservation management of small populations

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

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

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.

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.

Questioning inbreeding: Could outbreeding affect productivity in the North African catfish in Thailand?

The North African catfish (Clarias gariepinus) is a significant species in aquaculture, which is crucial for ensuring food and nutrition security. Their high adaptability to diverse environments has led to an increase in the number of farms that are available for their production. However, long-term closed breeding adversely affects their reproductive performance, leading to a decrease in production efficiency. This is possibly caused by inbreeding depression. To investigate the root cause of this issue, the genetic diversity of captive North African catfish populations was assessed in this study. Microsatellite genotyping and mitochondrial DNA D-loop sequencing were applied to 136 catfish specimens, collected from three populations captured for breeding in Thailand. Interestingly, extremely low inbreeding coefficients were obtained within each population, and distinct genetic diversity was observed among the three populations, indicating that their genetic origins are markedly different. This suggests that outbreeding depression by genetic admixture among currently captured populations of different origins may account for the low productivity of the North African catfish in Thailand. Genetic improvement of the North African catfish populations is required by introducing new populations whose origins are clearly known. This strategy should be systematically integrated into breeding programs to establish an ideal founder stock for selective breeding.

Current State of Mugger Populations

The mugger (Crocodylus palustris) is a medium-sized crocodilian inhabiting South Asia. As a result of intensive hunting, its range declined drastically up till the 1970s. Currently, the world mugger population is fragmented and threatened mainly by habitat loss and the consequences of human-crocodile conflict, being classified as Vulnerable by the IUCN. The goal of this paper is to comprehensively determine the mugger's current range, and assess risks in notable habitats of the species across its range. To determine the range and notable habitats, extensive literature covering surveys, monitoring, population studies and reports of human-crocodile conflict was examined. Habitat suitability and risk assessment were performed by evaluating selected habitats using eight factors: the legal status of the area, elevation, surface water availability, water quality, salinity, availability of nesting and basking sites, interaction with humans and interspecific competition. Based on our findings, the chances of the mugger's survival varies greatly across its range and the threats they face are complex and often site-specific. Defining these threats is the first step for determining suitable risk mitigation efforts, some of which are explored in this review.

Genome-Wide Detection for Runs of Homozygosity in Baoshan Pigs Using Whole Genome Resequencing

Baoshan pigs (BS) are a local breed in Yunnan Province that may face inbreeding owing to its limited population size. To accurately evaluate the inbreeding level of the BS pig population, we used whole-genome resequencing to identify runs of homozygosity (ROH) regions in BS pigs, calculated the inbreeding coefficient based on pedigree and ROH, and screened candidate genes with important economic traits from ROH islands. A total of 22,633,391 SNPS were obtained from the whole genome of BS pigs, and 201 ROHs were detected from 532,450 SNPS after quality control. The number of medium-length ROH (1-5 Mb) was the highest (98.43%), the number of long ROH (>5 Mb) was the lowest (1.57%), and the inbreeding of BS pigs mainly occurred in distant generations. The inbreeding coefficient FROH, calculated based on ROH, was 0.018 ± 0.016, and the FPED, calculated based on the pedigree, was 0.027 ± 0.028, which were positively correlated. Forty ROH islands were identified, containing 507 genes and 891 QTLs. Several genes were associated with growth and development (IGFALS, PTN, DLX5, DKK1, WNT2), meat quality traits (MC3R, ACSM3, ECI1, CD36, ROCK1, CACNA2D1), and reproductive traits (NPW, TSHR, BMP7). This study provides a reference for the protection and utilization of BS pigs.

Demographic history shapes North American gray wolf genomic diversity and informs species' conservation

Effective population size estimates are critical information needed for evolutionary predictions and conservation decisions. This is particularly true for species with social factors that restrict access to breeding or experience repeated fluctuations in population size across generations. We investigated the genomic estimates of effective population size along with diversity, subdivision, and inbreeding from 162,109 minimally filtered and 81,595 statistically neutral and unlinked SNPs genotyped in 437 grey wolf samples from North America collected between 1986 and 2021. We found genetic structure across North America, represented by three distinct demographic histories of western, central, and eastern regions of the continent. Further, grey wolves in the northern Rocky Mountains have lower genomic diversity than wolves of the western Great Lakes and have declined over time. Effective population size estimates revealed the historical signatures of continental efforts of predator extermination, despite a quarter century of recovery efforts. We are the first to provide molecular estimates of effective population size across distinct grey wolf populations in North America, which ranged between Ne  ~ 275 and 3050 since early 1980s. We provide data that inform managers regarding the status and importance of effective population size estimates for grey wolf conservation, which are on average 5.2-9.3% of census estimates for this species. We show that while grey wolves fall above minimum effective population sizes needed to avoid extinction due to inbreeding depression in the short term, they are below sizes predicted to be necessary to avoid long-term risk of extinction.

Genomic insights into endangerment and conservation of the garlic-fruit tree (Malania oleifera), a plant species with extremely small populations

BACKGROUND

Advanced whole-genome sequencing techniques enable covering nearly all genome nucleotide variations and thus can provide deep insights into protecting endangered species. However, the use of genomic data to make conservation strategies is still rare, particularly for endangered plants. Here we performed comprehensive conservation genomic analysis for Malania oleifera, an endangered tree species with a high amount of nervonic acid. We used whole-genome resequencing data of 165 samples, covering 16 populations across the entire distribution range, to investigate the formation reasons of its extremely small population sizes and to evaluate the possible genomic offsets and changes of ecology niche suitability under future climate change.

RESULTS

Although M. oleifera maintains relatively high genetic diversity among endangered woody plants (θπ = 3.87 × 10-3), high levels of inbreeding have been observed, which have reduced genetic diversity in 3 populations (JM, NP, and BM2) and caused the accumulation of deleterious mutations. Repeated bottleneck events, recent inbreeding (∼490 years ago), and anthropogenic disturbance to wild habitats have aggravated the fragmentation of M. oleifera and made it endangered. Due to the significant effect of higher average annual temperature, populations distributed in low altitude exhibit a greater genomic offset. Furthermore, ecological niche modeling shows the suitable habitats for M. oleifera will decrease by 71.15% and 98.79% in 2100 under scenarios SSP126 and SSP585, respectively.

CONCLUSIONS

The basic realizations concerning the threats to M. oleifera provide scientific foundation for defining management and adaptive units, as well as prioritizing populations for genetic rescue. Meanwhile, we highlight the importance of integrating genomic offset and ecological niche modeling to make targeted conservation actions under future climate change. Overall, our study provides a paradigm for genomics-directed conservation.

Genomics reveals introgression and purging of deleterious mutations in the Arabian leopard (Panthera pardus nimr)

In endangered species, low-genetic variation and inbreeding result from recent population declines. Genetic screenings in endangered populations help to assess their vulnerability to extinction and to create informed management actions toward their conservation efforts. The leopard, Panthera pardus, is a highly generalist predator with currently eight different subspecies. Yet, genomic data are still lacking for the Critically Endangered Arabian leopard (P. p. nimr). Here, we sequenced the whole genome of two Arabian leopards and assembled the most complete genomic dataset for leopards to date. Our phylogenomic analyses show that leopards are divided into two deeply divergent clades: the African and the Asian. Conservation genomic analyses indicate a prolonged population decline, which has led to an increase in inbreeding and runs of homozygosity, with consequent purging of deleterious mutations in both Arabian individuals. Our study represents the first attempt to genetically inform captive breeding programmes for this Critically Endangered subspecies.

Management of diversity and inbreeding when importing new stock into an inbred population

This article relates to breeding programs that seek to manage genetic diversity. The method maximizes a multicomponent objective function, applicable across breeding scenarios. However, this paper focuses on breeding decisions following immigration of 10 unrelated individuals into a highly inbred simulated population (F ≈ 0.34). We use Optimal Contribution Selection to maximize retention of genetic diversity. However, some treatments add Coancestry Assortative Mating (CAM). This helps to avoid early dilution of immigrant genetic material, maximizing its ability to contribute to genetic diversity in the longer term. After 20 generations, this resulted in considerably increased genetic diversity, with mean coancestries 59% of what random pairing gave. To manage progeny inbreeding, common practice is to reject matings above an upper limit. As a suboptimal rules-based approach, this resulted in 26% decreased genetic diversity and 8% increased inbreeding in the long term, compared with random pairing. In contrast, including mean progeny inbreeding as a continuous variable in the overall objective function decreased final inbreeding by 37% compared with random pairing. Adding some emphasis on selection for a single trait resulted in a similar pattern of effects on coancestry and inbreeding, with 12% higher trait response under CAM. Results indicate the properties of alternative methods, but we encourage users to do their own investigations of particular scenarios, such as including inbreeding depression. Practical implementation of these methods is discussed: they have been widely adopted in domestic animal breeding and are highly flexible to accommodate a wide range of technical and logistical objectives and constraints.

Monitoring Genetic Erosion of Aromatic and Medicinal Plant Species in Alentejo (South Portugal)

The main goal of this work was to study the genetic erosion risk of plants with aromatic, medicinal and gastronomic applications in Portugal, particularly in the Alentejo region. The target species were coriander (Coriandrum sativum L.), hart's pennyroyal (Mentha cervina L.) and pennyroyal (Mentha pulegium L.). The methodology involved direct observations and surveys (2002/2003 and 2011). The GE formula applied in Hammer's studies was used to estimate genetic erosion. The main factors causing genetic erosion were the primary drivers of biodiversity loss: habitat loss, invasive species, and overexploitation influenced by human intervention such as the clearing of watercourses, vegetation control, grazing and desertification. The results indicate a reduction in individuals per species in Alentejo, with a net erosion loss of 11% for M. pulegium, 32% for M. cervina and 33% for C. sativum. The overall loss of accessions (genetic erosion risk) was higher in cultivated accessions (33%) than in wild accessions (11%), with an annual genetic erosion rate of 3.7% and 1.2%, respectively. The annual risk of genetic erosion for M. pulegium accessions collected in a natural habitat was 0.6%, which is much lower than the 3.7% for M. cervina. These results consolidate the importance of collecting and conserving genetic resources.

Whole-genome survey reveals extensive variation in genetic diversity and inbreeding levels among peregrine falcon subspecies

In efforts to prevent extinction, resource managers are often tasked with increasing genetic diversity in a population of concern to prevent inbreeding depression or improve adaptive potential in a changing environment. The assumption that all small populations require measures to increase their genetic diversity may be unwarranted, and limited resources for conservation may be better utilized elsewhere. We test this assumption in a case study focused on the peregrine falcon (Falco peregrinus), a cosmopolitan circumpolar species with 19 named subspecies. We used whole-genome resequencing to generate over two million single nucleotide polymorphisms (SNPs) from multiple individuals of all peregrine falcon subspecies. Our analyses revealed extensive variation among subspecies, with many island-restricted and nonmigratory populations possessing lower overall genomic diversity, elevated inbreeding coefficients (F _ROH)-among the highest reported, and extensive runs of homozygosity (ROH) compared to mainland and migratory populations. Similarly, the majority of subspecies that are either nonmigratory or restricted to islands show a much longer history of low effective population size (N _e). While mutational load analyses indicated an increased proportion of homozygous-derived deleterious variants (i.e., drift load) among nonmigrant and island populations compared to those that are migrant or reside on the mainland, no significant differences in the proportion of heterozygous deleterious variants (i.e., inbreeding load) was observed. Our results provide evidence that high levels of inbreeding may not be an existential threat for some populations or taxa. Additional factors such as the timing and severity of population declines are important to consider in management decisions about extinction potential.

Non-random mating within an Island rookery of Hawaiian hawksbill turtles: demographic discontinuity at a small coastline scale

Hawksbill sea turtles (Eretmochelys imbricata) from the Hawaiian archipelago form a small and genetically isolated population, consisting of only a few tens of individuals breeding annually. Most females nest on the island of Hawai'i, but little is known about the demographics of this rookery. This study used genetic relatedness, inferred from 135 microhaplotype markers, to determine breeding sex-ratios, estimate female nesting frequency and assess relationships between individuals nesting on different beaches. Samples were collected during the 2017 nesting season and final data included 13 nesting females and 1002 unhatched embryos, salvaged from 41 nests, of which 13 had no observed mother. Results show that most females used a single nesting beach laying 1-5 nests each. From female and offspring alleles, the paternal genotypes of 12 breeding males were reconstructed and many showed high relatedness to their mates. Pairwise relatedness of offspring revealed one instance of polygyny but otherwise suggested a 1 : 1 breeding-sex ratio. Relatedness analysis and spatial-autocorrelation of genotypes indicate that turtles from different nesting areas do not regularly interbreed, suggesting that strong natal homing tendencies in both sexes result in non-random mating across the study area. Complexes of nearby nesting beaches also showed unique patterns of inbreeding across loci, further indicating that Hawaiian hawksbill turtles have demographically discontinuous nesting populations separated by only tens of km.

Developing Methods for Maintaining Genetic Diversity in Novel Aquaculture Species: The Case of Seriola lalandi

Developing sound breeding programs for aquaculture species may be challenging when matings cannot be controlled due to communal spawning. We developed a genotyping-by-sequencing marker panel of 300 SNPs for parentage testing and sex determination by using data from an in-house reference genome as well as a 90 K SNP genotyping array based on different populations of yellowtail kingfish (Seriola lalandi). The minimum and maximum distance between adjacent marker pairs were 0.7 Mb and 13 Mb, respectively, with an average marker spacing of 2 Mb. Weak evidence of the linkage disequilibrium between adjacent marker pairs was found. The results showed high panel performance for parental assignment, with probability exclusion values equaling 1. The rate of false positives when using cross-population data was null. A skewed distribution of genetic contributions by dominant females was observed, thus increasing the risk of higher rates of inbreeding in subsequent captive generations when no parentage data are used. All these results are discussed in the context of breeding program design, using this marker panel to increase the sustainability of this aquaculture resource.

Chromosome-scale genomes reveal genomic consequences of inbreeding in the South China tiger: A comparative study with the Amur tiger

The South China tiger (Panthera tigris amoyensis, SCT) is the most critically endangered subspecies of tiger due to functional extinction in the wild. Inbreeding depression is observed among the captive population descended from six wild ancestors, resulting in high juvenile mortality and low reproduction. We assembled and characterized the first SCT genome and an improved Amur tiger (P. t. altaica, AT) genome named AmyTig1.0 and PanTig2.0. The two genomes are the most continuous and comprehensive among any tiger genomes yet reported at the chromosomal level. By using the two genomes and resequencing data of 15 SCT and 13 AT individuals, we investigated the genomic signature of inbreeding depression of the SCT. The results indicated that the effective population size of SCT experienced three phases of decline, ~5.0-1.0 thousand years ago, 100 years ago, and since captive breeding in 1963. We found 43 long runs of homozygosity fragments that were shared by all individuals in the SCT population and covered a total length of 20.63% in the SCT genome. We also detected a large proportion of identical-by-descent segments across the genome in the SCT population, especially on ChrB4. Deleterious nonsynonymous single nucleotide polymorphic sites and loss-of-function mutations were found across genomes with extensive potential influences, despite a proportion of these loads having been purged by inbreeding depression. Our research provides an invaluable resource for the formulation of genetic management policies for the South China tiger such as developing genome-based breeding and genetic rescue strategy.

A field theoretic approach to non-equilibrium population genetics in the strong selection regime

Natural populations are virtually never observed in equilibrium, yet equilibrium approximations comprise the majority of our understanding of population genetics. Using standard tools from statistical physics, a formalism is presented that re-expresses the stochastic equations describing allelic evolution as a partition functional over all possible allelic trajectories ('paths') governed by selection, mutation, and drift. A perturbative field theory is developed for strong additive selection, relevant to disease variation, that facilitates the straightforward computation of closed-form approximations for time-dependent moments of the allele frequency distribution across a wide range of non-equilibrium scenarios; examples are presented for constant population size, exponential growth, bottlenecks, and oscillatory size, all of which align well to simulations and break down just above the drift barrier. Equilibration times are computed and, even for static population size, generically extend beyond the order 1/s timescale associated with exponential frequency decay. Though the mutation load is largely robust to variable population size, perturbative drift-based corrections to the deterministic trajectory are readily computed. Under strong selection, the variance of a new mutation's frequency (related to homozygosity) is dominated by drift-driven dynamics and a transient increase in variance often occurs prior to equilibrating. The excess kurtosis over skew squared is roughly constant (i.e., independent of selection, provided 2Ns ≳ 5) for static population size, and thus potentially sensitive to deviation from equilibrium. These insights highlight the value of such closed-form approximations, naturally generated from Feynman diagrams in a perturbative field theory, which can simply and accurately capture the parameter dependences describing a variety of non-equilibrium population genetic phenomena of interest.

Transposable elements maintain genome-wide heterozygosity in inbred populations

Elevated levels of inbreeding increase the risk of inbreeding depression and extinction, yet many inbred species are widespread, suggesting that inbreeding has little impact on evolutionary potential. Here, we explore the potential for transposable elements (TEs) to maintain genetic variation in functional genomic regions under extreme inbreeding. Capitalizing on the mixed mating system of Arabidopsis lyrata, we assess genome-wide heterozygosity and signatures of selection at single nucleotide polymorphisms near transposable elements across an inbreeding gradient. Under intense inbreeding, we find systematically elevated heterozygosity downstream of several TE superfamilies, associated with signatures of balancing selection. In addition, we demonstrate increased heterozygosity in stress-responsive genes that consistently occur downstream of TEs. We finally reveal that TE superfamilies are associated with specific signatures of selection that are reproducible across independent evolutionary lineages of A. lyrata. Together, our study provides an important hypothesis for the success of self-fertilizing species.

Genomic trajectories of a near-extinction event in the Chatham Island black robin

BACKGROUND

Understanding the micro--evolutionary response of populations to demographic declines is a major goal in evolutionary and conservation biology. In small populations, genetic drift can lead to an accumulation of deleterious mutations, which will increase the risk of extinction. However, demographic recovery can still occur after extreme declines, suggesting that natural selection may purge deleterious mutations, even in extremely small populations. The Chatham Island black robin (Petroica traversi) is arguably the most inbred bird species in the world. It avoided imminent extinction in the early 1980s and after a remarkable recovery from a single pair, a second population was established and the two extant populations have evolved in complete isolation since then. Here, we analysed 52 modern and historical genomes to examine the genomic consequences of this extreme bottleneck and the subsequent translocation.

RESULTS

We found evidence for two-fold decline in heterozygosity and three- to four-fold increase in inbreeding in modern genomes. Moreover, there was partial support for temporal reduction in total load for detrimental variation. In contrast, compared to historical genomes, modern genomes showed a significantly higher realised load, reflecting the temporal increase in inbreeding. Furthermore, the translocation induced only small changes in the frequency of deleterious alleles, with the majority of detrimental variation being shared between the two populations.

CONCLUSION

Our results highlight the dynamics of mutational load in a species that recovered from the brink of extinction, and show rather limited temporal changes in mutational load. We hypothesise that ancestral purging may have been facilitated by population fragmentation and isolation on several islands for thousands of generations and may have already reduced much of the highly deleterious load well before human arrival and introduction of pests to the archipelago. The majority of fixed deleterious variation was shared between the modern populations, but translocation of individuals with low mutational load could possibly mitigate further fixation of high-frequency deleterious variation.

Reproductive biology and population structure of the endangered shrub Grevillea bedggoodiana (Proteaceae)

Narrowly endemic species are particularly vulnerable to catastrophic events. Compared to widespread species, they may also be less capable of adapting to shifts in environmental pressures as a result of specialisation on a narrow range of local condition and limited ability to disperse. However, life-history traits, such as preferential outcrossing and high fecundity can maintain genetic diversity and evolutionary potential, and boost species resilience. The endangered Grevillea bedggoodiana (Enfield Grevillea) is an understorey shrub restricted to an area of ca. 150 km2 in south-eastern Australia with a legacy of large-scale anthropogenic disturbance. Prior to this study little was known about its biology and population structure. Here, its breeding system was assessed through a controlled pollination experiment at one of its central populations, and eight populations were sampled for genetic analysis with microsatellite markers. The species was found to be preferentially outcrossing, with no evidence of pollination limitation. In most populations, allelic richness, observed heterozygosity and gene diversity were high (Ar: 3.8–6.3; Ho: 0.45–0.65, He: 0.60 − 0.75). However, the inbreeding coefficients were significant in at least four populations, ranging from Fi -0.061 to 0.259 despite high outcrossing rates. Estimated reproductive rates varied among sampled populations but were independent of gene diversity and inbreeding. Despite its small geographic range, the species’ populations showed moderate differentiation (AMOVA: FST = 0.123), which was largely attributable to isolation by distance. We interpret these results as suggesting that G. bedggoodiana is reproductively healthy and has maintained high levels of genetic diversity despite recent disturbance.

Eco-evolutionary extinction and recolonization dynamics reduce genetic load and increase time to extinction in highly inbred populations

Understanding how genetic and ecological effects can interact to shape genetic loads within and across local populations is key to understanding ongoing persistence of systems that should otherwise be susceptible to extinction through mutational meltdown. Classic theory predicts short persistence times for metapopulations comprising small local populations with low connectivity, due to accumulation of deleterious mutations. Yet, some such systems have persisted over evolutionary time, implying the existence of mechanisms that allow metapopulations to avoid mutational meltdown. We first hypothesize a mechanism by which the combination of stochasticity in the numbers and types of mutations arising locally (genetic stochasticity), resulting local extinction, and recolonization through evolving dispersal facilitates metapopulation persistence. We then test this mechanism using a spatially and genetically explicit individual-based model. We show that genetic stochasticity in highly structured metapopulations can result in local extinctions, which can favor increased dispersal, thus allowing recolonization of empty habitat patches. This causes fluctuations in metapopulation size and transient gene flow, which reduces genetic load and increases metapopulation persistence over evolutionary time. Our suggested mechanism and simulation results provide an explanation for the conundrum presented by the continued persistence of highly structured populations with inbreeding mating systems that occur in diverse taxa.

Are genetic variation and demographic performance linked?

Quantifying relationships between genetic variation and population viability is important from both basic biological and applied conservation perspectives, yet few populations have been monitored with both long-term demographic and population genetics approaches. To empirically test whether and how genetic variation and population dynamics are related, we present one such paired approach. First, we use eight years of historical demographic data from five populations of Boechera fecunda (Brassicaceae), a rare, self-compatible perennial plant endemic to Montana, USA, and use integral projection models to estimate the stochastic population growth rate (λ S) and extinction risk of each population. We then combine these demographic estimates with previously published metrics of genetic variation in the same populations to test whether genetic diversity within populations is linked to demographic performance. Our results show that in this predominantly inbred species, standing genetic variation and demography are weakly positively correlated. However, the inbreeding coefficient was not strongly correlated with demographic performance, suggesting that more inbred populations are not necessarily less viable or at higher extinction risk than less inbred populations. A contemporary re-census of these populations revealed that neither genetic nor demographic parameters were consistently strong predictors of current population density, although populations showing lower probabilities of extinction in demographic models had larger population sizes at present. In the absence of evidence for inbreeding depression decreasing population viability in this species, we recommend conservation of distinct, potentially locally adapted populations of B. fecunda rather than alternatives such as translocations or reintroductions.

Genomic analysis of the rare British Lop pig and identification of distinctive genomic markers

Concentration of production on a few commercial pig breeds has led to the marginalization of many native, numerically small breeds, increasing their risk of endangerment. In the UK, one such rare breed is the British Lop, a lop-eared breed, of similar origin to the Welsh breed. The objective of the present study was to address the genomic status of the British Lop and its relationship with other breeds and identify a small set of genomic markers that uniquely characterize and distinguish British Lop animals. Results have shown that the British Lop is a relatively distinct population with reduced genomic diversity and effective size consistent with its status as a rare breed. Furthermore, we demonstrated the genetic closeness of the British Lop to phenotypically similar breeds such as Landrace and Welsh as well Large White, Middle White and Pietrain. Finally, a set of 75 Single Nucleotide Polymorphisms distributed across multiple chromosomes were identified and validated as markers that can consistently distinguish British Lops from other closely related breeds. Results may inform breeding and management strategies aiming to enhance diversity as well as the development of a breed purity test.

Effects of genetic vs. environmental quality on condition-dependent morphological and life history traits in a neriid fly

Condition is assumed to reflect both genes and environment, enabling condition‐dependent signals to reveal genetic quality. However, because the phenotypic effects of variation in genetic quality could be masked by environmental heterogeneity, the contribution of genetic quality to phenotypic variation in fitness‐related traits and condition‐dependent signals remains unclear. We compared effects of ecologically relevant manipulations of environmental quality (nutrient dilution in the larval diet) and genetic quality (one generation of inbreeding) on male and female morphology, life history and reproductive performance in the neriid fly Telostylinus angusticollis. We found that larval diet quality had strong, positive effects on male and female body size, male secondary sexual traits, and aspects of male and female reproductive performance. By contrast, inbreeding had weak effects on most traits, and no trait showed clear and consistent effects of both environmental and genetic quality. Indeed, inbreeding effects on body size and male competitive performance were of opposite sign in rich vs. poor larval diet treatment groups. Our results suggest that environmental quality strongly affects condition, but the effects of genetic quality are subtle and environment‐dependent in this species. These findings raise questions about the genetic architecture of condition and the potential for condition‐dependent traits to function as signals of genetic quality.

Evolution of selfing syndrome and its influence on genetic diversity and inbreeding: A range-wide study in Oenothera primiveris

PREMISE

To avoid inbreeding depression, plants have evolved diverse breeding systems to favor outcrossing, such as self-incompatibility. However, changes in biotic and abiotic conditions can result in selective pressures that lead to a breakdown in self-incompatibility. The shift to increased selfing is commonly associated with reduced floral features, lower attractiveness to pollinators, and increased inbreeding. We tested the hypothesis that the loss of self-incompatibility, a shift to self-fertilization (autogamy), and concomitant evolution of the selfing syndrome (reduction in floral traits associated with cross-fertilization) will lead to increased inbreeding and population differentiation in Oenothera primiveris. Across its range, this species exhibits a shift in its breeding system and floral traits from a self-incompatible population with large flowers to self-compatible populations with smaller flowers.

METHODS

We conducted a breeding system assessment, evaluated floral traits in the field and under controlled conditions, and measured population genetic parameters using RADseq data.

RESULTS

Our results reveal a bimodal transition to the selfing syndrome from the west to the east of the range of O. primiveris. This shift includes variation in the breeding system and the mating system, a reduction in floral traits (flower diameter, herkogamy, and scent production), a shift to greater autogamy, reduced genetic diversity, and increased inbreeding.

CONCLUSIONS

The observed variation highlights the importance of range-wide studies to understand breeding system variation and the evolution of the selfing syndrome within populations and species.

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

Genetic variation, which is generated by mutation, recombination and gene flow, can reduce the mean fitness of a population, both now and in the future. This 'genetic load' has been estimated in a wide range of animal taxa using various approaches. Advances in genome sequencing and computational techniques now enable us to estimate the genetic load in populations and individuals without direct fitness estimates. Here, we review the classic and contemporary literature of genetic load. We describe approaches to quantify the genetic load in whole-genome sequence data based on evolutionary conservation and annotations. We show that splitting the load into its two components - the realized load (or expressed load) and the masked load (or inbreeding load) - can improve our understanding of the population genetics of deleterious mutations.

New developments in the field of genomic technologies and their relevance to conservation management

Recent technological advances in the field of genomics offer conservation managers and practitioners new tools to explore for conservation applications. Many of these tools are well developed and used by other life science fields, while others are still in development. Considering these technological possibilities, choosing the right tool(s) from the toolbox is crucial and can pose a challenging task. With this in mind, we strive to inspire, inform and illuminate managers and practitioners on how conservation efforts can benefit from the current genomic and biotechnological revolution. With inspirational case studies we show how new technologies can help resolve some of the main conservation challenges, while also informing how implementable the different technologies are. We here focus specifically on small population management, highlight the potential for genetic rescue, and discuss the opportunities in the field of gene editing to help with adaptation to changing environments. In addition, we delineate potential applications of gene drives for controlling invasive species. We illuminate that the genomic toolbox offers added benefit to conservation efforts, but also comes with limitations for the use of these novel emerging techniques.

Estimates of recent and historical effective population size in turbot, seabream, seabass and carp selective breeding programmes

Background

The high fecundity of fish species allows intense selection to be practised and therefore leads to fast genetic gains. Based on this, numerous selective breeding programmes have been started in Europe in the last decades, but in general, little is known about how the base populations of breeders have been built. Such knowledge is important because base populations can be created from very few individuals, which can lead to small effective population sizes and associated reductions in genetic variability. In this study, we used genomic information that was recently made available for turbot (Scophthalmus maximus), gilthead seabream (Sparus aurata), European seabass (Dicentrarchus labrax) and common carp (Cyprinus carpio) to obtain accurate estimates of the effective size for commercial populations.

Methods

Restriction-site associated DNA sequencing data were used to estimate current and historical effective population sizes. We used a novel method that considers the linkage disequilibrium spectrum for the whole range of genetic distances between all pairs of single nucleotide polymorphisms (SNPs), and thus accounts for potential fluctuations in population size over time.

Results

Our results show that the current effective population size for these populations is small (equal to or less than 50 fish), potentially putting the sustainability of the breeding programmes at risk. We have also detected important drops in effective population size about five to nine generations ago, most likely as a result of domestication and the start of selective breeding programmes for these species in Europe.

Conclusions

Our findings highlight the need to broaden the genetic composition of the base populations from which selection programmes start, and suggest that measures designed to increase effective population size within all farmed populations analysed here should be implemented in order to manage genetic variability and ensure the sustainability of the breeding programmes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00680-9.

Not out of the woods yet: Signatures of the prolonged negative genetic consequences of a population bottleneck in a rapidly re-expanding wader, the black-faced spoonbill Platalea minor

The long-term persistence of a population which has suffered a bottleneck partly depends on how historical demographic dynamics impacted its genetic diversity and the accumulation of deleterious mutations. Here we provide genomic evidence for the genetic effect of a recent population bottleneck in the endangered black-faced spoonbill (Platalea minor) after its rapid population recovery. Our data suggest that the bird's effective population size, N_e , had been relatively stable (7500-9000) since 22,000 years ago; however, a recent brief yet severe bottleneck (N_e  = 20) which we here estimated to occur around the 1940s wiped out >99% of its historical N_e in roughly three generations. Despite a >15-fold population recovery since 1988, we found that black-faced spoonbill population has higher levels of inbreeding (7.4 times more runs of homozygosity) than its sister species, the royal spoonbill (P. regia), which is not thought to have undergone a marked population contraction. Although the two spoonbills have similar levels of genome-wide genetic diversity, our results suggest that selection on more genes was relaxed in the black-faced spoonbill; moreover individual black-faced spoonbills carry more putatively deleterious mutations (Grantham's score > 50), and may therefore express more deleterious phenotypic effects than royal spoonbills. Here we demonstrate the value of using genomic indices to monitor levels of genetic erosion, inbreeding and mutation load in species with conservation concerns. To mitigate the prolonged negative genetic effect of a population bottleneck, we recommend that all possible measures should be employed to maintain population growth of a threatened species.

Genetic purging in captive endangered ungulates with extremely low effective population sizes

Inbreeding threatens the survival of small populations by producing inbreeding depression, but also exposes recessive deleterious effects in homozygosis allowing for genetic purging. Using inbreeding-purging theory, we analyze early survival in four pedigreed captive breeding programs of endangered ungulates where population growth was prioritized so that most adult females were allowed to contribute offspring according to their fitness. We find evidence that purging can substantially reduce inbreeding depression in Gazella cuvieri (with effective population size N_e = 14) and Nanger dama (N_e = 11). No purging is detected in Ammotragus lervia (N_e = 4), in agreement with the notion that drift overcomes purging under fast inbreeding, nor in G. dorcas (N_e = 39) where, due to the larger population size, purging is slower and detection is expected to require more generations. Thus, although smaller populations are always expected to show smaller fitness (as well as less adaptive potential) than larger ones due to higher homozygosis and deleterious fixation, our results show that a substantial fraction of their inbreeding load and inbreeding depression can be purged when breeding contributions are governed by natural selection. Since management strategies intended to maximize the ratio from the effective to the actual population size tend to reduce purging, the search for a compromise between these strategies and purging could be beneficial in the long term. This could be achieved either by allowing some level of random mating and some role of natural selection in determining breeding contributions, or by undertaking reintroductions into the wild at the earliest opportunity.

Reviewing the consequences of genetic purging on the success of rescue programs

Genetic rescue is increasingly considered a promising and underused conservation strategy to reduce inbreeding depression and restore genetic diversity in endangered populations, but the empirical evidence supporting its application is limited to a few generations. Here we discuss on the light of theory the role of inbreeding depression arising from partially recessive deleterious mutations and of genetic purging as main determinants of the medium to long-term success of rescue programs. This role depends on two main predictions: (1) The inbreeding load hidden in populations with a long stable demography increases with the effective population size; and (2) After a population shrinks, purging tends to remove its (partially) recessive deleterious alleles, a process that is slower but more efficient for large populations than for small ones. We also carry out computer simulations to investigate the impact of genetic purging on the medium to long term success of genetic rescue programs. For some scenarios, it is found that hybrid vigor followed by purging will lead to sustained successful rescue. However, there may be specific situations where the recipient population is so small that it cannot purge the inbreeding load introduced by migrants, which would lead to increased fitness inbreeding depression and extinction risk in the medium to long term. In such cases, the risk is expected to be higher if migrants came from a large non-purged population with high inbreeding load, particularly after the accumulation of the stochastic effects ascribed to repeated occasional migration events. Therefore, under the specific deleterious recessive mutation model considered, we conclude that additional caution should be taken in rescue programs. Unless the endangered population harbors some distinctive genetic singularity whose conservation is a main concern, restoration by continuous stable gene flow should be considered, whenever feasible, as it reduces the extinction risk compared to repeated occasional migration and can also allow recolonization events.

Genomic signatures of inbreeding in a critically endangered parrot, the kākāpō

Events of inbreeding are inevitable in critically endangered species. Reduced population sizes and unique life-history traits can increase the severity of inbreeding, leading to declines in fitness and increased risk of extinction. Here, we investigate levels of inbreeding in a critically endangered flightless parrot, the kākāpō (Strigops habroptilus), wherein a highly inbred island population and one individual from the mainland of New Zealand founded the entire extant population. Genotyping-by-sequencing (GBS), and a genotype calling approach using a chromosome-level genome assembly, identified a filtered set of 12,241 single-nucleotide polymorphisms (SNPs) among 161 kākāpō, which together encompass the total genetic potential of the extant population. Multiple molecular-based estimates of inbreeding were compared, including genome-wide estimates of heterozygosity (FH), the diagonal elements of a genomic-relatedness matrix (FGRM), and runs of homozygosity (RoH, FRoH). In addition, we compared levels of inbreeding in chicks from a recent breeding season to examine if inbreeding is associated with offspring survival. The density of SNPs generated with GBS was sufficient to identify chromosomes that were largely homozygous with RoH distributed in similar patterns to other inbred species. Measures of inbreeding were largely correlated and differed significantly between descendants of the two founding populations. However, neither inbreeding nor ancestry was found to be associated with reduced survivorship in chicks, owing to unexpected mortality in chicks exhibiting low levels of inbreeding. Our study highlights important considerations for estimating inbreeding in critically endangered species, such as the impacts of small population sizes and admixture between diverse lineages.

Long-term exhaustion of the inbreeding load in Drosophila melanogaster

Inbreeding depression, the decline in fitness of inbred individuals, is a ubiquitous phenomenon of great relevance in evolutionary biology and in the fields of animal and plant breeding and conservation. Inbreeding depression is due to the expression of recessive deleterious alleles that are concealed in heterozygous state in noninbred individuals, the so-called inbreeding load. Genetic purging reduces inbreeding depression by removing these alleles when expressed in homozygosis due to inbreeding. It is generally thought that fast inbreeding (such as that generated by full-sib mating lines) removes only highly deleterious recessive alleles, while slow inbreeding can also remove mildly deleterious ones. However, a question remains regarding which proportion of the inbreeding load can be removed by purging under slow inbreeding in moderately large populations. We report results of two long-term slow inbreeding Drosophila experiments (125-234 generations), each using a large population and a number of derived lines with effective sizes about 1000 and 50, respectively. The inbreeding load was virtually exhausted after more than one hundred generations in large populations and between a few tens and over one hundred generations in the lines. This result is not expected from genetic drift alone, and is in agreement with the theoretical purging predictions. Computer simulations suggest that these results are consistent with a model of relatively few deleterious mutations of large homozygous effects and partially recessive gene action.

Population genomics of the critically endangered kākāpō

The kākāpō is a flightless parrot endemic to New Zealand. Once common in the archipelago, only 201 individuals remain today, most of them descending from an isolated island population. We report the first genome-wide analyses of the species, including a high-quality genome assembly for kākāpō, one of the first chromosome-level reference genomes sequenced by the Vertebrate Genomes Project (VGP). We also sequenced and analyzed 35 modern genomes from the sole surviving island population and 14 genomes from the extinct mainland population. While theory suggests that such a small population is likely to have accumulated deleterious mutations through genetic drift, our analyses on the impact of the long-term small population size in kākāpō indicate that present-day island kākāpō have a reduced number of harmful mutations compared to mainland individuals. We hypothesize that this reduced mutational load is due to the island population having been subjected to a combination of genetic drift and purging of deleterious mutations, through increased inbreeding and purifying selection, since its isolation from the mainland ∼10,000 years ago. Our results provide evidence that small populations can survive even when isolated for hundreds of generations. This work provides key insights into kākāpō breeding and recovery and more generally into the application of genetic tools in conservation efforts for endangered species.

Landscape barriers to pollen and seed flow in the dioecious tropical tree Astronium fraxinifolium in Brazilian savannah

Gene flow studies provide information on gene exchange between populations, which is essential for developing genetic conservation strategies. Such analyses enable a better understanding of the life history and seed and pollen dispersal mechanisms of plant species. In this study, we investigate pollen and seed flow in a regenerant population of the pioneer species Astronium fraxinifolium in an area degraded during the construction of a hydroelectric dam. We mapped, sampled, sexed, and genotyped 386 individuals in the regenerant population (RP), as well as 128 adult trees located along two highways adjacent to the degraded area; one in Mato Grosso do Sul State (MS) and other in São Paulo State (SP). Parentage analyses was carried out for 370 individuals of the RP population, using as putative parents 348 individuals from RP and all 128 individuals sampled in MS and SP. Based on parentage analysis and eight microsatellite loci, our analyses revealed that for individuals of the RP with an identified father (pollen donor), 1.1% of the pollen was dispersed up to 532 m, while for those with an identified mother (seed donor), 0.5% of seeds were dispersed up to 4,782 m. However, a large proportion of pollen (76.5%) and seeds (57%) immigrated from trees outside the sampled populations. Pollen and seeds were dispersed through a pattern of isolation by distance. Genetic diversity was significantly similar between adults of both highway populations and individuals from RP, with significant levels of inbreeding detected only in RP. Our results demonstrate that the nearest trees contributed pollen and seeds for the recovery of the degraded area, indicating reproductive spatial isolation among the sampled populations due to the damming of the river. Such results help to understand the process of regeneration for A. fraxinifolium in regenerant populations to inform strategies for conservation and environmental recovery with this species.

Genomic data of different resolutions reveal consistent inbreeding estimates but contrasting homozygosity landscapes for the threatened Aotearoa New Zealand hihi

Inbreeding can lead to a loss of heterozygosity in a population and when combined with genetic drift may reduce the adaptive potential of a species. However, there is uncertainty about whether resequencing data can provide accurate and consistent inbreeding estimates. Here, we performed an in-depth inbreeding analysis for hihi (Notiomystis cincta), an endemic and nationally vulnerable passerine bird of Aotearoa New Zealand. We first focused on subsampling variants from a reference genome male, and found that low-density data sets tend to miss runs of homozygosity (ROH) in some places and overestimate ROH length in others, resulting in contrasting homozygosity landscapes. Low-coverage resequencing and 50 K SNP array densities can yield comparable inbreeding results to high-coverage resequencing approaches, but the results for all data sets are highly dependent on the software settings employed. Second, we extended our analysis to 10 hihi where low-coverage whole genome resequencing, RAD-seq and SNP array genotypes are available. We inferred ROH and individual inbreeding to evaluate the relative effects of sequencing depth versus SNP density on estimating inbreeding coefficients and found that high rates of missingness downwardly bias both the number and length of ROH. In summary, when using genomic data to evaluate inbreeding, studies must consider that ROH estimates are heavily dependent on analysis parameters, data set density and individual sequencing depth.

Genetic and environmental modulation of transposition shapes the evolutionary potential of Arabidopsis thaliana

BACKGROUND

How species can adapt to abrupt environmental changes, particularly in the absence of standing genetic variation, is poorly understood and a pressing question in the face of ongoing climate change. Here we leverage publicly available multi-omic and bio-climatic data for more than 1000 wild Arabidopsis thaliana accessions to determine the rate of transposable element (TE) mobilization and its potential to create adaptive variation in natural settings.

RESULTS

We demonstrate that TE insertions arise at almost the same rate as base substitutions. Mobilization activity of individual TE families varies greatly between accessions, in association with genetic and environmental factors as well as through complex gene-environment interactions. Although the distribution of TE insertions across the genome is ultimately shaped by purifying selection, reflecting their typically strong deleterious effects when located near or within genes, numerous recent TE-containing alleles show signatures of positive selection. Moreover, high rates of transposition appear positively selected at the edge of the species' ecological niche. Based on these findings, we predict through mathematical modeling higher transposition activity in Mediterranean regions within the next decades in response to global warming, which in turn should accelerate the creation of large-effect alleles.

CONCLUSIONS

Our study reveals that TE mobilization is a major generator of genetic variation in A. thaliana that is finely modulated by genetic and environmental factors. These findings and modeling indicate that TEs may be essential genomic players in the demise or rescue of native populations in times of climate crises.

Genomic insights into the conservation status of the world's last remaining Sumatran rhinoceros populations

Small populations are often exposed to high inbreeding and mutational load that can increase the risk of extinction. The Sumatran rhinoceros was widespread in Southeast Asia, but is now restricted to small and isolated populations on Sumatra and Borneo, and most likely extinct on the Malay Peninsula. Here, we analyse 5 historical and 16 modern genomes from these populations to investigate the genomic consequences of the recent decline, such as increased inbreeding and mutational load. We find that the Malay Peninsula population experienced increased inbreeding shortly before extirpation, which possibly was accompanied by purging. The populations on Sumatra and Borneo instead show low inbreeding, but high mutational load. The currently small population sizes may thus in the near future lead to inbreeding depression. Moreover, we find little evidence for differences in local adaptation among populations, suggesting that future inbreeding depression could potentially be mitigated by assisted gene flow among populations.

The inflated significance of neutral genetic diversity in conservation genetics

The current rate of species extinction is rapidly approaching unprecedented highs, and life on Earth presently faces a sixth mass extinction event driven by anthropogenic activity, climate change, and ecological collapse. The field of conservation genetics aims at preserving species by using their levels of genetic diversity, usually measured as neutral genome-wide diversity, as a barometer for evaluating population health and extinction risk. A fundamental assumption is that higher levels of genetic diversity lead to an increase in fitness and long-term survival of a species. Here, we argue against the perceived importance of neutral genetic diversity for the conservation of wild populations and species. We demonstrate that no simple general relationship exists between neutral genetic diversity and the risk of species extinction. Instead, a better understanding of the properties of functional genetic diversity, demographic history, and ecological relationships is necessary for developing and implementing effective conservation genetic strategies.

M. Groenen M, Madsen O. Genetic consequences of long-term small effective population size in the critically endangered pygmy hog

Increasing human disturbance and climate change have a major impact on habitat integrity and size, with far-reaching consequences for wild fauna and flora. Specifically, population decline and habitat fragmentation result in small, isolated populations. To what extend different endangered species can cope with small population size is still largely unknown. Studies on the genomic landscape of these species can shed light on past demographic dynamics and current genetic load, thereby also providing guidance for conservation programs. The pygmy hog (Porcula salvania) is the smallest and rarest wild pig in the world, with current estimation of only a few hundred living in the wild. Here, we analyzed whole-genome sequencing data of six pygmy hogs, three from the wild and three from a captive population, along with 30 pigs representing six other Suidae. First, we show that the pygmy hog had a very small population size with low genetic diversity over the course of the past ~1 million years. One indication of historical small effective population size is the absence of mitochondrial variation in the six sequenced individuals. Second, we evaluated the impact of historical demography. Runs of homozygosity (ROH) analysis suggests that the pygmy hog population has gone through past but not recent inbreeding. Also, the long-term, extremely small population size may have led to the accumulation of harmful mutations suggesting that the accumulation of deleterious mutations is exceeding purifying selection in this species. Thus, care has to be taken in the conservation program to avoid or minimize the potential for further inbreeding depression, and guard against environmental changes in the future.

Strongly deleterious mutations are a primary determinant of extinction risk due to inbreeding depression

Human-driven habitat fragmentation and loss have led to a proliferation of small and isolated plant and animal populations with high risk of extinction. One of the main threats to extinction in these populations is inbreeding depression, which is primarily caused by recessive deleterious mutations becoming homozygous due to inbreeding. The typical approach for managing these populations is to maintain high genetic diversity, increasingly by translocating individuals from large populations to initiate a "genetic rescue." However, the limitations of this approach have recently been highlighted by the demise of the gray wolf population on Isle Royale, which declined to the brink of extinction soon after the arrival of a migrant from the large mainland wolf population. Here, we use a novel population genetic simulation framework to investigate the role of genetic diversity, deleterious variation, and demographic history in mediating extinction risk due to inbreeding depression in small populations. We show that, under realistic models of dominance, large populations harbor high levels of recessive strongly deleterious variation due to these mutations being hidden from selection in the heterozygous state. As a result, when large populations contract, they experience a substantially elevated risk of extinction after these strongly deleterious mutations are exposed by inbreeding. Moreover, we demonstrate that, although genetic rescue is broadly effective as a means to reduce extinction risk, its effectiveness can be greatly increased by drawing migrants from small or moderate-sized source populations rather than large source populations due to smaller populations harboring lower levels of recessive strongly deleterious variation. Our findings challenge the traditional conservation paradigm that focuses on maximizing genetic diversity in small populations in favor of a view that emphasizes minimizing strongly deleterious variation. These insights have important implications for managing small and isolated populations in the increasingly fragmented landscape of the Anthropocene.

On the estimation of inbreeding depression using different measures of inbreeding from molecular markers

The inbreeding coefficient (F) of individuals can be estimated from molecular marker data, such as SNPs, using measures of homozygosity of individual markers or runs of homozygosity (ROH) across the genome. These different measures of F can then be used to estimate the rate of inbreeding depression (ID) for quantitative traits. Some recent simulation studies have investigated the accuracy of this estimation with contradictory results. Whereas some studies suggest that estimates of inbreeding from ROH account more accurately for ID, others suggest that inbreeding measures from SNP-by-SNP homozygosity giving a large weight to rare alleles are more accurate. Here, we try to give more light on this issue by carrying out a set of computer simulations considering a range of population genetic parameters and population sizes. Our results show that the previous studies are indeed not contradictory. In populations with low effective size, where relationships are more tight and selection is relatively less intense, F measures based on ROH provide very accurate estimates of ID whereas SNP-by-SNP-based F measures with high weight to rare alleles can show substantial upwardly biased estimates of ID. However, in populations of large effective size, with more intense selection and trait allele frequencies expected to be low if they are deleterious for fitness because of purifying selection, average estimates of ID from SNP-by-SNP-based F values become unbiased or slightly downwardly biased and those from ROH-based F values become slightly downwardly biased. The noise attached to all these estimates, nevertheless, can be very high in large-sized populations. We also investigate the relationship between the different F measures and the homozygous mutation load, which has been suggested as a proxy of inbreeding depression.

Multiple life-stage inbreeding depression impacts demography and extinction risk in an extinct-in-the-wild species

Inbreeding can depress individuals’ fitness traits and reduce population viability. However, studies that directly translate inbreeding depression on fitness traits into consequences for population viability, and further, into consequences for management choices, are lacking. Here, we estimated impacts of inbreeding depression (B, lethal equivalents) across life-history stages for an extinct-in-the-wild species, the sihek (Guam kingfisher, Todiramphus cinnamominus). We then projected population growth under different management alternatives with our B estimates incorporated, as well as without inbreeding depression (B = 0) or with a conventional default B. We found that inbreeding depression severely impacted multiple life-history stages, and directly translated into an effect on population viability under management alternatives. Simulations including our B estimates indicated rapid population decline, whereas projections without inbreeding depression or with default B suggested very gradual population decline. Further, our results demonstrate that incorporation of B across life-history stages can influence management decisions, as projections with our B estimates suggested a need to switch to increased breeding management to avoid species extinction and support wild releases. Our results demonstrate that magnitude of B across life-history stages can translate into demographic consequences, such that incorporation of multiple life-stage B into population models can be important for informed conservation management decision-making.

Inbreeding depression for kit survival at birth in a rabbit population under long-term selection

Background

Accumulation of detrimental mutations in small populations leads to inbreeding depression of fitness traits and a higher frequency of genetic defects, thus increasing risk of extinction. Our objective was to quantify the magnitude of inbreeding depression for survival at birth, in a closed rabbit population under long-term selection.

Methods

We used an information theory-based approach and multi-model inference to estimate inbreeding depression and its purging with respect to the trait ‘kit survival at birth’ over a 25-year period in a closed population of Pannon White rabbits, by analysing 22,718 kindling records. Generalised linear mixed models based on the logit link function were applied, which take polygenic random effects into account.

Results

Our results indicated that inbreeding depression occurred during the period 1992–1997, based on significant estimates of the z-standardised classical inbreeding coefficient z.F_L (CI_95% − 0.12 to − 0.03) and of the new inbreeding coefficient of the litter z.F_NEWL (CI_95% − 0.13 to − 0.04) as well as a 59.2% reduction in contributing founders. Inbreeding depression disappeared during the periods 1997–2007 and 2007–2017. For the period 1992–1997, the best model resulted in a significantly negative standardised estimate of the new inbreeding coefficient of the litter and a significantly positive standardised estimate of Kalinowski’s ancestral inbreeding coefficient of the litter (CI_95% 0.01 to 0.17), which indicated purging of detrimental load. Kindling season and parity had effects on survival at birth that differed across the three periods of 1992–1997, 1997–2007 and 2007–2017.

Conclusions

Our results support the existence of inbreeding depression and its purging with respect to kit survival at birth in this Pannon White rabbit population. However, we were unable to exclude possible confounding from the effects of parity and potentially other environmental factors during the study period, thus our results need to be extended and confirmed in other populations.

Accelerated inbreeding depression suggests synergistic epistasis for deleterious mutations in Drosophila melanogaster

Epistasis may have important consequences for a number of issues in quantitative genetics and evolutionary biology. In particular, synergistic epistasis for deleterious alleles is relevant to the mutation load paradox and the evolution of sex and recombination. Some studies have shown evidence of synergistic epistasis for spontaneous or induced deleterious mutations appearing in mutation-accumulation experiments. However, many newly arising mutations may not actually be segregating in natural populations because of the erasing action of natural selection. A demonstration of synergistic epistasis for naturally segregating alleles can be achieved by means of inbreeding depression studies, as deleterious recessive allelic effects are exposed in inbred lines. Nevertheless, evidence of epistasis from these studies is scarce and controversial. In this paper, we report the results of two independent inbreeding experiments carried out with two different populations of Drosophila melanogaster. The results show a consistent accelerated inbreeding depression for fitness, suggesting synergistic epistasis among deleterious alleles. We also performed computer simulations assuming different possible models of epistasis and mutational parameters for fitness, finding some of them to be compatible with the results observed. Our results suggest that synergistic epistasis for deleterious mutations not only occurs among newly arisen spontaneous or induced mutations, but also among segregating alleles in natural populations.

Origin and Evolution of Deleterious Mutations in Horses

Domestication has changed the natural evolutionary trajectory of horses by favoring the reproduction of a limited number of animals showing traits of interest. Reduced breeding stocks hampered the elimination of deleterious variants by means of negative selection, ultimately inflating mutational loads. However, ancient genomics revealed that mutational loads remained steady during most of the domestication history until a sudden burst took place some 250 years ago. To identify the factors underlying this trajectory, we gather an extensive dataset consisting of 175 modern and 153 ancient genomes previously published, and carry out the most comprehensive characterization of deleterious mutations in horses. We confirm that deleterious variants segregated at low frequencies during the last 3500 years, and only spread and incremented their occurrence in the homozygous state during modern times, owing to inbreeding. This independently happened in multiple breeds, following both the development of closed studs and purebred lines, and the deprecation of horsepower in the 20th century, which brought many draft breeds close to extinction. Our work illustrates the paradoxical effect of some conservation and improvement programs, which reduced the overall genomic fitness and viability.

Strategies of the invasive tropical fire ant (Solenopsis geminata) to minimize inbreeding costs

How invasive species overcome challenges associated with low genetic diversity is unclear. Invasive ant populations with low genetic diversity sometimes produce sterile diploid males, which do not contribute to colony labour or reproductive output. We investigated how inbreeding affects colony founding and potential strategies to overcome its effects in the invasive tropical fire ant, Solenopsis geminata. Our genetic analyses of field samples revealed that 13-100% of males per colony (n = 8 males per 10 colonies) were diploid, and that all newly mated queens (n = 40) were single-mated. Our laboratory experiment in which we assigned newly mated queens to nests consisting of 1, 2, 3, or 5 queens (n = 95 ± 9 replicates) revealed that pleometrosis (queens founding their nest together) and diploid male larvae execution can compensate for diploid male load. The proportion of diploid male producing (DMP) colonies was 22.4%, and DMP colonies produced fewer pupae and adult workers than non-DMP colonies. Pleometrosis significantly increased colony size. Queens executed their diploid male larvae in 43.5% of the DMP colonies, and we hypothesize that cannibalism benefits incipient colonies because queens can redirect nutrients to worker brood. Pleometrosis and cannibalism of diploid male larvae represent strategies through which invasive ants can successfully establish despite high inbreeding.