Inbreeding depression across the lifespan in a wild mammal population

Little evidence of inbreeding depression for birth mass, survival and growth in Antarctic fur seal pups

Inbreeding depression, the loss of offspring fitness due to consanguineous mating, is generally detrimental for individual performance and population viability. We investigated inbreeding effects in a declining population of Antarctic fur seals (Arctocephalus gazella) at Bird Island, South Georgia. Here, localised warming has reduced the availability of the seal's staple diet, Antarctic krill, leading to a temporal increase in the strength of selection against inbred offspring, which are increasingly failing to recruit into the adult breeding population. However, it remains unclear whether selection operates before or after nutritional independence at weaning. We therefore used microsatellite data from 885 pups and their mothers, and SNP array data from 98 mother-offspring pairs, to quantify the effects of individual and maternal inbreeding on three important neonatal fitness traits: birth mass, survival and growth. We did not find any clear or consistent effects of offspring or maternal inbreeding on any of these traits. This suggests that selection filters inbred individuals out of the population as juveniles during the time window between weaning and recruitment. Our study brings into focus a poorly understood life-history stage and emphasises the importance of understanding the ecology and threats facing juvenile pinnipeds.

Detecting inbreeding depression in structured populations

Measuring inbreeding and its consequences on fitness is central for many areas in biology including human genetics and the conservation of endangered species. However, there is no consensus on the best method, neither for quantification of inbreeding itself nor for the model to estimate its effect on specific traits. We simulated traits based on simulated genomes from a large pedigree and empirical whole-genome sequences of human data from populations with various sizes and structures (from the 1,000 Genomes project). We compare the ability of various inbreeding coefficients ([Formula: see text]) to quantify the strength of inbreeding depression: allele-sharing, two versions of the correlation of uniting gametes which differ in the weight they attribute to each locus and two identical-by-descent segments-based estimators. We also compare two models: the standard linear model and a linear mixed model (LMM) including a genetic relatedness matrix (GRM) as random effect to account for the nonindependence of observations. We find LMMs give better results in scenarios with population or family structure. Within the LMM, we compare three different GRMs and show that in homogeneous populations, there is little difference among the different [Formula: see text] and GRM for inbreeding depression quantification. However, as soon as a strong population or family structure is present, the strength of inbreeding depression can be most efficiently estimated only if i) the phenotypes are regressed on [Formula: see text] based on a weighted version of the correlation of uniting gametes, giving more weight to common alleles and ii) with the GRM obtained from an allele-sharing relatedness estimator.

Genetic architecture of inbreeding depression may explain its persistence in a population of wild red deer

Inbreeding depression is of major concern in declining populations, but relatively little is known about its genetic architecture in wild populations, such as the degree to which it is composed of large or small effect loci and their distribution throughout the genome. Here, we combine fitness and genomic data from a wild population of red deer to investigate the genomic distribution of inbreeding effects. Based on the runs of homozygosity (ROH)-based inbreeding coefficient, FROH, we use chromosome-specific inbreeding coefficients (FROHChr) to explore whether the effect of inbreeding varies between chromosomes. Under the assumption that within an individual the probability of being identical-by-descent is equal across all chromosomes, we used a multi-membership model to estimate the deviation of FROHChr from the average inbreeding effect. This novel approach ensures effect sizes are not overestimated whilst maximising the power of our available dataset of >3000 individuals genotyped on >35,000 autosomal SNPs. We find that most chromosomes confer a minor reduction in fitness-related traits, which when these effects are summed, results in the observed inbreeding depression in birth weight, survival and lifetime breeding success. However, no chromosomes had a significant detrimental effect compared to the overall effect of inbreeding, indicating no major effect loci. We conclude that in this population, inbreeding depression is likely the result of multiple mildly or moderately deleterious mutations spread across all chromosomes, which are difficult to detect with statistical confidence. Such mutations will be inefficiently purged, which may explain the persistence of inbreeding depression in this population.

Early-acting inbreeding depression can evolve as an inbreeding avoidance mechanism

Despite the potential for mechanical, developmental and/or chemical mechanisms to prevent self-fertilization, incidental self-fertilization is inevitable in many predominantly outcrossing species. In such cases, inbreeding can compromise individual fitness. Unquestionably, much of this inbreeding depression is maladaptive. However, we show that when reproductive compensation allows for the replacement of inviable embryos lost early in development, selection can favour deleterious recessive variants that induce 'self-sacrificial' death of inbred embryos. Our theoretical results provide numerous testable predictions which could challenge the assumption that inbreeding depression is always maladaptive. Our work is applicable any species that cannot fully avoid inbreeding, exhibits substantial inbreeding depression, and has the potential to compensate embryos lost early in development. In addition to its general applicability, our theory suggests that self-sacrificial variants might be responsible for the remarkably low realized selfing rates of gymnosperms with high primary selfing rates, as gymnosperms exhibit strong inbreeding depression, have effective reproductive compensation mechanisms, and cannot evolve chemical self-incompatibility.

Development and evaluation of a novel single nucleotide polymorphism panel for North American bison

Genome-wide single nucleotide polymorphism (SNP) genotyping platforms have become increasingly popular in characterizing livestock and wildlife populations, replacing traditional methods such as microsatellite fragment analysis. Herein, we report the development and evaluation of a novel bison SNP panel for population management and conservation. Initially, 2474 autosomal SNPs were selected from existing bison whole-genome sequences and variable sites among bison on the GGSP bovine 50K Chip, based on minor allele frequency, data completeness, and chromosome location. Additionally, 20 mitochondrial SNPs were chosen to identify known mitochondrial haplotypes in bison according to previous research. The SNPs were further evaluated using genotyping-by-sequencing with 190 bison, representing the historical lineages that survived the major population crash of the late 1800s. Variants with high potential for genotyping error were filtered out, and the remaining SNPs were placed on a custom Illumina™ array. The final panel consisting of 798 autosomal and 13 mitochondrial SNPs was used to establish baseline genetic parameters, compare populations, and assign mitochondrial haplotypes in 995 bison across ten populations. These SNPs were also found to be highly informative for individual animal identification and parentage assignment. This SNP panel provides a powerful new method to establish a baseline for estimating genetic health of bison populations and a new tool for bison managers to make informed management decisions based on genetic information specific to their populations.

Landscape genomics analysis provides insights into future climate change-driven risk in rhesus macaque

Climate change significantly affects the suitability of wildlife habitats. Thus, understanding how animals adapt ecologically and genetically to climate change is important for targeted species protection. Rhesus macaques (Macaca mulatta) are widely distributed and multi-climatically adapted primates. This study explored how rhesus macaques adapt to climate change by integrating ecological and genetic methods and applying species distribution models (SDMs) and a gradient forest (GF) model. The findings suggested that temperature seasonality primarily affects habitat suitability and indicated that climate change will have a dramatic impact on macaque populations in the future. We also applied genotype-environment association (GEA) analyses and selection signature analyses to identify genes associated with climate change and provide possible explanations for the adaptation of rhesus macaques to climatic environments. The population genomics analyses suggested that the Taihang population has the highest genomic vulnerability with inbreeding and low heterozygosity. Combined with the higher ecological vulnerability, additional conservation strategies are required for this population under higher risk of climate change. Our work measured the impact of climate change and enabled the identification of populations that exhibit high vulnerability to severe climate change. Such information is useful for selecting populations of rhesus macaques as subject of long-term monitoring or evolutionary rescue under future climate change.

Multi-locus homozygosity promotes actuarial senescence in a wild mammal

Genome-wide homozygosity, caused for example by inbreeding, is expected to have deleterious effects on survival and/or reproduction. Evolutionary theory predicts that any fitness costs are likely to be detected in late life because natural selection will filter out negative impacts on younger individuals with greater reproductive value. Here we infer associations between multi-locus homozygosity (MLH), sex, disease and age-dependent mortality risks using Bayesian analysis of the life histories of wild European badgers Meles meles in a population naturally infected with Mycobacterium bovis (the causative agent of bovine tuberculosis [bTB]). We find important effects of MLH on all parameters of the Gompertz-Makeham mortality hazard function, but particularly in later life. Our findings confirm the predicted association between genomic homozygosity and actuarial senescence. Increased homozygosity is particularly associated with an earlier onset, and greater rates of actuarial senescence, regardless of sex. The association between homozygosity and actuarial senescence is further amplified among badgers putatively infected with bTB. These results recommend further investigation into the ecological and behavioural processes that result in genome-wide homozygosity, and focused work on whether homozygosity is harmful or beneficial during early life-stages.

No evidence of inbreeding depression despite a historical severe bottleneck in the endangered Bermuda petrel (Pterodroma cahow)

The Bermuda petrel Pterodroma cahow is an island endemic seabird that belongs to the Procellariiformes, one of the most endangered orders of birds. Historical records suggest a significant population size decline following human settlement in Bermuda, bringing the species to near extinction. Since the 1950s, the population has been recovering aided by the implementation of an ongoing conservation plan. However, it still faces several threats, and negative genetic effects resulting from that drastic decline are to be expected, including inbreeding and genetic drift. We studied genetic diversity and levels of inbreeding, and their effects on individual fitness and mating choice. We also tested for a genetic signature of the recent demographic bottleneck. For this, we analyzed variation in thousands of nuclear single-nucleotide polymorphisms derived from double digest restriction site-associated DNA sequencing and 1 mitochondrial gene (cytochrome oxidase I). The results revealed that the Bermuda petrel suffered a recent genetic bottleneck and shows low mitochondrial diversity compared with other petrel species. Conversely, nuclear diversity was similar to that of other endangered petrels. Inbreeding levels were not high overall, although some individuals were highly inbred. However, we found no evidence that individual inbreeding or relatedness between mates affected hatching success, or that mate choice is influenced by kinship in this very small population.

Demographic history and genomic consequences of 10,000 generations of isolation in a wild mammal

Increased human activities caused the isolation of populations in many species-often associated with genetic depletion and negative fitness effects. The effects of isolation are predicted by theory, but long-term data from natural populations are scarce. We show, with full genome sequences, that common voles (Microtus arvalis) in the Orkney archipelago have remained genetically isolated from conspecifics in continental Europe since their introduction by humans over 5,000 years ago. Modern Orkney vole populations are genetically highly differentiated from continental conspecifics as a result of genetic drift processes. Colonization likely started on the biggest Orkney island and vole populations on smaller islands were gradually split off, without signs of secondary admixture. Despite having large modern population sizes, Orkney voles are genetically depauperate and successive introductions to smaller islands resulted in further reduction of genetic diversity. We detected high levels of fixation of predicted deleterious variation compared with continental populations, particularly on smaller islands, yet the fitness effects realized in nature are unknown. Simulations showed that predominantly mildly deleterious mutations were fixed in populations, while highly deleterious mutations were purged early in the history of the Orkney population. Relaxation of selection overall due to benign environmental conditions on the islands and the effects of soft selection may have contributed to the repeated, successful establishment of Orkney voles despite potential fitness loss. Furthermore, the specific life history of these small mammals, resulting in relatively large population sizes, has probably been important for their long-term persistence in full isolation.

A polygenic basis for birth weight in a wild population of red deer (Cervus elaphus)

The genetic architecture of traits under selection has important consequences for the response to selection and potentially for population viability. Early QTL mapping studies in wild populations have reported loci with large effect on trait variation. However, these results are contradicted by more recent genome-wide association analyses, which strongly support the idea that most quantitative traits have a polygenic basis. This study aims to re-evaluate the genetic architecture of a key morphological trait, birth weight, in a wild population of red deer (Cervus elaphus), using genomic approaches. A previous study using 93 microsatellite and allozyme markers and linkage mapping on a kindred of 364 deer detected a pronounced QTL on chromosome 21 explaining 29% of the variance in birth weight, suggesting that this trait is partly controlled by genes with large effects. Here, we used data for more than 2,300 calves genotyped at >39,000 SNP markers and two approaches to characterise the genetic architecture of birth weight. First, we performed a genome-wide association (GWA) analysis, using a genomic relatedness matrix to account for population structure. We found no SNPs significantly associated with birth weight. Second, we used genomic prediction to estimate the proportion of variance explained by each SNP and chromosome. This analysis confirmed that most genetic variance in birth weight was explained by loci with very small effect sizes. Third, we found that the proportion of variance explained by each chromosome was slightly positively correlated with its size. These three findings highlight a highly polygenic architecture for birth weight, which contradicts the previous QTL study. These results are probably explained by the differences in how associations are modelled between QTL mapping and GWA. Our study suggests that models of polygenic adaptation are the most appropriate to study the evolutionary trajectory of this trait.

Selection, recombination and population history effects on runs of homozygosity (ROH) in wild red deer (Cervus elaphus)

The distribution of runs of homozygosity (ROH) may be shaped by a number of interacting processes such as selection, recombination and population history, but little is known about the importance of these mechanisms in shaping ROH in wild populations. We combined an empirical dataset of >3000 red deer genotyped at >35,000 genome-wide autosomal SNPs and evolutionary simulations to investigate the influence of each of these factors on ROH. We assessed ROH in a focal and comparison population to investigate the effect of population history. We investigated the role of recombination using both a physical map and a genetic linkage map to search for ROH. We found differences in ROH distribution between both populations and map types indicating that population history and local recombination rate have an effect on ROH. Finally, we ran forward genetic simulations with varying population histories, recombination rates and levels of selection, allowing us to further interpret our empirical data. These simulations showed that population history has a greater effect on ROH distribution than either recombination or selection. We further show that selection can cause genomic regions where ROH is common, only when the effective population size (Ne) is large or selection is particularly strong. In populations having undergone a population bottleneck, genetic drift can outweigh the effect of selection. Overall, we conclude that in this population, genetic drift resulting from a historical population bottleneck is most likely to have resulted in the observed ROH distribution, with selection possibly playing a minor role.

Effect of MHC and inbreeding on disassortative reproduction: A data revisit, extension and inclusion of fertilization in sand lizards

The harmful effects of close inbreeding have been recognized for centuries and, with the rise of Mendelian genetics, was realized to be an effect of homozygosis. This historical background led to great interest in ways to quantify inbreeding, its depression effects on the phenotype and flow‐on effects on mate choice and other aspects of behavioral ecology. The mechanisms and cues used to avoid inbreeding are varied and include major histocompatibility complex (MHC) molecules and the peptides they transport as predictors of the degree of genetic relatedness. Here, we revisit and complement data from a Swedish population of sand lizards (Lacerta agilis) showing signs of inbreeding depression to assess the effects of genetic relatedness on pair formation in the wild. Parental pairs were less similar at the MHC than expected under random mating but mated at random with respect to microsatellite relatedness. MHC clustered in groups of RFLP bands but no partner preference was observed with respect to partner MHC cluster genotype. Male MHC band patterns were unrelated to their fertilization success in clutches selected for analysis on the basis of showing mixed paternity. Thus, our data suggest that MHC plays a role in pre‐copulatory, but not post‐copulatory partner association, suggesting that MHC is not the driver of fertilization bias and gamete recognition in sand lizards.

Inbreeding depression explains killer whale population dynamics

Understanding the factors that cause endangered populations to either grow or decline is crucial for preserving biodiversity. Conservation efforts often address extrinsic threats, such as environmental degradation and overexploitation, that can limit the recovery of endangered populations. Genetic factors such as inbreeding depression can also affect population dynamics but these effects are rarely measured in the wild and thus often neglected in conservation efforts. Here we show that inbreeding depression strongly influences the population dynamics of an endangered killer whale population, despite genomic signatures of purging of deleterious alleles via natural selection. We find that the 'Southern Residents', which are currently endangered despite nearly 50 years of conservation efforts, exhibit strong inbreeding depression for survival. Our population models suggest that this inbreeding depression limits population growth and predict further decline if the population remains genetically isolated and typical environmental conditions continue. The Southern Residents also had more inferred homozygous deleterious alleles than three other, growing, populations, further suggesting that inbreeding depression affects population fitness. These results demonstrate that inbreeding depression can substantially limit the recovery of endangered populations. Conservation actions focused only on extrinsic threats may therefore fail to account for key intrinsic genetic factors that also limit population growth.

Spatial genetic structure of European wild boar, with inferences on late-Pleistocene and Holocene demographic history

European wildlife has been subjected to intensifying levels of anthropogenic impact throughout the Holocene, yet the main genetic partitioning of many species is thought to still reflect the late-Pleistocene glacial refugia. We analyzed 26,342 nuclear SNPs of 464 wild boar (Sus scrofa) across the European continent to infer demographic history and reassess the genetic consequences of natural and anthropogenic forces. We found that population fragmentation, inbreeding and recent hybridization with domestic pigs have caused the spatial genetic structure to be heterogeneous at the local scale. Underlying local anthropogenic signatures, we found a deep genetic structure in the form of an arch-shaped cline extending from the Dinaric Alps, via Southeastern Europe and the Baltic states, to Western Europe and, finally, to the genetically diverged Iberian peninsula. These findings indicate that, despite considerable anthropogenic influence, the deeper, natural continental structure is still intact. Regarding the glacial refugia, our findings show a weaker signal than generally assumed, but are nevertheless suggestive of two main recolonization routes, with important roles for Southern France and the Balkans. Our results highlight the importance of applying genomic resources and framing genetic results within a species' demographic history and geographic distribution for a better understanding of the complex mixture of underlying processes.

Genomic signatures of inbreeding depression for a threatened Aotearoa New Zealand passerine

For small and isolated populations, the increased chance of mating between related individuals can result in a substantial reduction in individual and population fitness. Despite the increasing availability of genomic data to measure inbreeding accurately across the genome, inbreeding depression studies for threatened species are still scarce due to the difficulty of measuring fitness in the wild. Here, we investigate inbreeding and inbreeding depression for the extensively monitored Tiritiri Mātangi island population of a threatened Aotearoa New Zealand passerine, the hihi (Notiomystis cincta). First, using a custom 45 k single nucleotide polymorphism (SNP) array, we explore genomic inbreeding patterns by inferring homozygous segments across the genome. Although all individuals have similar levels of ancient inbreeding, highly inbred individuals are affected by recent inbreeding, which can probably be explained by bottleneck effects such as habitat loss after European arrival and their translocation to the island in the 1990s. Second, we investigate genomic inbreeding effects on fitness, measured as lifetime reproductive success, and its three components, juvenile survival, adult annual survival and annual reproductive success, in 363 hihi. We find that global inbreeding significantly affects juvenile survival but none of the remaining fitness traits. Finally, we employ a genome-wide association approach to test the locus-specific effects of inbreeding on fitness, and identify 13 SNPs significantly associated with lifetime reproductive success. Our findings suggest that inbreeding depression does impact hihi, but at different genomic scales for different traits, and that purging has therefore failed to remove all variants with deleterious effects from this population of conservation concern.

Effect of reduced genomic representation on using runs of homozygosity for inbreeding characterization

Genomic measures of inbreeding based on identical-by-descent (IBD) segments are increasingly used to measure inbreeding and mostly estimated on SNP arrays and whole-genome sequencing (WGS) data. However, some softwares recurrently used for their estimation assume that genomic positions which have not been genotyped are nonvariant. This might be true for WGS data, but not for reduced genomic representations and can lead to spurious IBD segments estimation. In this project, we simulated the outputs of WGS, two SNP arrays of different sizes and RAD-sequencing for three populations with different sizes and histories. We compare the results of IBD segments estimation with two softwares: runs of homozygosity (ROHs) estimated with PLINK and homozygous-by-descent (HBD) segments estimated with RZooRoH. We demonstrate that to obtain meaningful estimates of inbreeding, RZooRoH requires a SNPs density 11 times smaller compared to PLINK: ranks of inbreeding coefficients were conserved among individuals above 22 SNPs/Mb for PLINK and 2 SNPs/Mb for RZooRoH. We also show that in populations with simple demographic histories, distribution of ROHs and HBD segments are correctly estimated with both SNP arrays and WGS. PLINK correctly estimated distribution of ROHs with SNP densities above 22 SNPs/Mb, while RZooRoH correctly estimated distribution of HBD segments with SNPs densities above 11 SNPs/Mb. However, in a population with a more complex demographic history, RZooRoH resulted in better distribution of IBD segments estimation compared to PLINK even with WGS data. Consequently, we advise researchers to use either methods relying on excess homozygosity averaged across SNPs or model-based HBD segments calling methods for inbreeding estimations.

Transitioning from microsatellites to SNP-based microhaplotypes in genetic monitoring programmes: Lessons from paired data spanning 20 years

Many long-term genetic monitoring programmes began before next-generation sequencing became widely available. Older programmes can now transition to new marker systems usually consisting of 1000s of SNP loci, but there are still important questions about comparability, precision, and accuracy of key metrics estimated using SNPs. Ideally, transitioned programmes should capitalize on new information without sacrificing continuity of inference across the time series. We combined existing microsatellite-based genetic monitoring information with SNP-based microhaplotypes obtained from archived samples of Rio Grande silvery minnow (Hybognathus amarus) across a 20-year time series to evaluate point estimates and trajectories of key genetic metrics. Demographic and genetic monitoring bracketed multiple collapses of the wild population and included cases where captive-born repatriates comprised the majority of spawners in the wild. Even with smaller sample sizes, microhaplotypes yielded comparable and in some cases more precise estimates of variance genetic effective population size, multilocus heterozygosity and inbreeding compared to microsatellites because many more microhaplotype loci were available. Microhaplotypes also recorded shifts in allele frequencies associated with population bottlenecks. Trends in microhaplotype-based inbreeding metrics were associated with the fraction of hatchery-reared repatriates to the wild and should be incorporated into future genomic monitoring. Although differences in accuracy and precision of some metrics were observed between marker types, biological inferences and management recommendations were consistent.

A comparison of marker-based estimators of inbreeding and inbreeding depression

BACKGROUND

The availability of genome-wide marker data allows estimation of inbreeding coefficients (F, the probability of identity-by-descent, IBD) and, in turn, estimation of the rate of inbreeding depression (ΔID). We investigated, by computer simulations, the accuracy of the most popular estimators of inbreeding based on molecular markers when computing F and ΔID in populations under random mating, equalization of parental contributions, and artificially selected populations. We assessed estimators described by Li and Horvitz (F_LH1 and F_LH2), VanRaden (F_VR1 and F_VR2), Yang and colleagues (F_YA1 and F_YA2), marker homozygosity (F_HOM), runs of homozygosity (F_ROH) and estimates based on pedigree (F_PED) in comparison with estimates obtained from IBD measures (F_IBD).

RESULTS

If the allele frequencies of a base population taken as a reference for the computation of inbreeding are known, all estimators based on marker allele frequencies are highly correlated with F_IBD and provide accurate estimates of the mean ΔID. If base population allele frequencies are unknown and current frequencies are used in the estimations, the largest correlation with F_IBD is generally obtained by F_LH1 and the best estimator of ΔID is F_YA2. The estimators F_VR2 and F_LH2 have the poorest performance in most scenarios. The assumption that base population allele frequencies are equal to 0.5 results in very biased estimates of the average inbreeding coefficient but they are highly correlated with F_IBD and give relatively good estimates of ΔID. Estimates obtained directly from marker homozygosity (F_HOM) substantially overestimated ΔID. Estimates based on runs of homozygosity (F_ROH) provide accurate estimates of inbreeding and ΔID. Finally, estimates based on pedigree (F_PED) show a lower correlation with F_IBD than molecular estimators but provide rather accurate estimates of ΔID. An analysis of data from a pig population supports the main findings of the simulations.

CONCLUSIONS

When base population allele frequencies are known, all marker-allele frequency-based estimators of inbreeding coefficients generally show a high correlation with F_IBD and provide good estimates of ΔID. When base population allele frequencies are unknown, F_LH1 is the marker frequency-based estimator that is most correlated with F_IBD, and F_YA2 provides the most accurate estimates of ΔID. Estimates from F_ROH are also very precise in most scenarios. The estimators F_VR2 and F_LH2 have the poorest performances.

The Unusual Value of Long-Term Studies of Individuals: The Example of the Isle of Rum Red Deer Project

Long-term studies of individuals enable incisive investigations of questions across ecology and evolution. Here, we illustrate this claim by reference to our long-term study of red deer on the Isle of Rum, Scotland. This project has established many of the characteristics of social organization, selection, and population ecology typical of large, polygynous, seasonally breeding mammals, with wider implications for our understanding of sexual selection and the evolution of sex differences, as well as for their population dynamics and population management. As molecular genetic techniques have developed, the project has pivoted to investigate evolutionary genetic questions, also breaking new ground in this field. With ongoing advances in genomics and statistical approaches and the development of increasingly sophisticated ways to assay new phenotypic traits, the questions that long-term studies such as the red deer study can answer become both broader and ever more sophisticated. They also offer powerful means of understanding the effects of ongoing climate change on wild populations.

Contributions of distemper control and habitat expansion to the Amur leopard viability

The Amur leopard (Panthera pardus orientalis) is a critically endangered top predator that struggles on the brink of extinction due to threats such as canine distemper virus (CDV), habitat loss, and inbreeding depression. Here we develop a viability analysis metamodel that combines a traditional individual-based demographic model with an epidemiological model to assess the benefits of alternative population management actions in response to multiple distinct threats. Our results showed an extinction risk of 10.3%-99.9% if no management actions were taken over 100 years under different levels of inbreeding depression. Reducing the risk of CDV infection in Amur leopards through the low-coverage vaccination of leopards and the management of sympatric domestic dogs could effectively improve the survival probability of the leopard population, and with habitat expansion added to these management measures, the population expanded further. Our findings highlight that protecting the Amur leopard necessitates a multifaceted synergistic effort, and controlling multiple threats together may significantly escalate overall viability of a species, especially for small-isolated threatened population. More broadly, our modeling framework could offer critical perspectives and scientific support for conservation planning, as well as specific adaptive management actions for endangered species around the world.

In the absence of management strategies, canine distemper virus threatens the future existence of the endangered Amur leopard.

Using PVA and captive breeding to balance trade-offs in the rescue of the island dibbler onto a new island ark

In the face of the current global extinction crisis, it is critical we give conservation management strategies the best chance of success. Australia is not exempt from global trends with currently the world’s greatest mammal extinction rate (~ 1 per 8 years). Many more are threatened including the dibbler (Parantechinus apicalis) whose remnant range has been restricted to Western Australia at just one mainland site and two small offshore islands—Whitlock Island (5 ha) and Boullanger Island (35 ha). Here, we used 14 microsatellite markers to quantify genetic variation in the remaining island populations from 2013 to 2018 and incorporated these data into population viability analysis (PVA) models, used to assess factors important to dibbler survival and to provide guidance for translocations. Remnant population genetic diversity was low (< 0.3), and populations were highly divergent from each other (pairwise F_STs 0.29–0.52). Comparison of empirical data to an earlier study is consistent with recent declines in genetic diversity and models projected increasing extinction risk and declining genetic variation in the next century. Optimal translocation scenarios recommend 80 founders for new dibbler populations—provided by captive breeding—and determined the proportion of founders from parental populations to maximise genetic diversity and minimise harvesting impact. The goal of our approach is long-term survival of genetically diverse, self-sustaining populations and our methods are transferable. We consider mixing island with mainland dibblers to reinforce genetic variation.

Decomposing phenotypic skew and its effects on the predicted response to strong selection

The major frameworks for predicting evolutionary change assume that a phenotype's underlying genetic and environmental components are normally distributed. However, the predictions of these frameworks may no longer hold if distributions are skewed. Despite this, phenotypic skew has never been decomposed, meaning the fundamental assumptions of quantitative genetics remain untested. Here we demonstrate that the substantial phenotypic skew in the body size of juvenile blue tits (Cyanistes caeruleus) is driven by environmental factors. Although skew had little impact on our predictions of selection response in this case, our results highlight the impact of skew on the estimation of inheritance and selection. Specifically, the nonlinear parent-offspring regressions induced by skew, alongside selective disappearance, can strongly bias estimates of heritability. The ubiquity of skew and strong directional selection on juvenile body size imply that heritability is commonly overestimated, which may in part explain the discrepancy between predicted and observed trait evolution.

Genetic variance in fitness indicates rapid contemporary adaptive evolution in wild animals

The rate of adaptive evolution, the contribution of selection to genetic changes that increase mean fitness, is determined by the additive genetic variance in individual relative fitness. To date, there are few robust estimates of this parameter for natural populations, and it is therefore unclear whether adaptive evolution can play a meaningful role in short-term population dynamics. We developed and applied quantitative genetic methods to long-term datasets from 19 wild bird and mammal populations and found that, while estimates vary between populations, additive genetic variance in relative fitness is often substantial and, on average, twice that of previous estimates. We show that these rates of contemporary adaptive evolution can affect population dynamics and hence that natural selection has the potential to partly mitigate effects of current environmental change.

Dynamics of reduced genetic diversity in increasingly fragmented populations of Florida Scrub‐Jays , Aphelocoma coerulescens

Understanding the genomic consequences of population decline is important for predicting species' vulnerability to intensifying global change. Empirical information about genomic changes in populations in the early stages of decline, especially for those still experiencing immigration, remains scarce. We used 7834 autosomal SNPs and demographic data for 288 Florida scrub jays (Aphelocoma coerulescens; FSJ) sampled in 2000 and 2008 to compare levels of genetic diversity, inbreeding, relatedness, and lengths of runs of homozygosity (ROH) between two subpopulations within dispersal distance of one another but have experienced contrasting demographic trajectories. At Archbold Biological Station (ABS), the FSJ population has been stable because of consistent habitat protection and management, while at nearby Placid Lakes Estates (PLE), the population declined precipitously due to suburban development. By the onset of our sampling in 2000, birds in PLE were already less heterozygous, more inbred, and on average more related than birds in ABS. No significant changes occurred in heterozygosity or inbreeding across the 8‐year sampling interval, but average relatedness among individuals decreased in PLE, thus by 2008 average relatedness did not differ between sites. PLE harbored a similar proportion of short ROH but a greater proportion of long ROH than ABS, suggesting one continuous population of shared demographic history in the past, which is now experiencing more recent inbreeding. These results broadly uphold the predictions of simple population genetic models based on inferred effective population sizes and rates of immigration. Our study highlights how, in just a few generations, formerly continuous populations can diverge in heterozygosity and levels of inbreeding with severe local population decline despite ongoing gene flow.

Inbreeding is associated with shorter early-life telomere length in a wild passerine

Inbreeding can have negative effects on survival and reproduction, which may be of conservation concern in small and isolated populations. However, the physiological mechanisms underlying inbreeding depression are not well-known. The length of telomeres, the DNA sequences protecting chromosome ends, has been associated with health or fitness in several species. We investigated effects of inbreeding on early-life telomere length in two small island populations of wild house sparrows (Passer domesticus) known to be affected by inbreeding depression. Using genomic measures of inbreeding we found that inbred nestling house sparrows (n = 371) have significantly shorter telomeres. Using pedigree-based estimates of inbreeding we found a tendency for inbred nestling house sparrows to have shorter telomeres (n = 1195). This negative effect of inbreeding on telomere length may have been complemented by a heterosis effect resulting in longer telomeres in individuals that were less inbred than the population average. Furthermore, we found some evidence of stronger effects of inbreeding on telomere length in males than females. Thus, telomere length may reveal subtle costs of inbreeding in the wild and demonstrate a route by which inbreeding negatively impacts the physiological state of an organism already at early life-history stages.

Genetic management on the brink of extinction: sequencing microsatellites does not improve estimates of inbreeding in wild and captive Vancouver Island marmots (Marmota vancouverensis)

Captive breeding is often a last resort management option in the conservation of endangered species which can in turn lead to increased risk of inbreeding depression and loss of genetic diversity. Thus, recording breeding events via studbook for the purpose of estimating relatedness, and facilitating mating pair selection to minimize inbreeding, is common practice. However, as founder relatedness is often unknown, loss of genetic variation and inbreeding cannot be entirely avoided. Molecular genotyping is slowly being adopted in captive breeding programs, however achieving sufficient resolution can be challenging in small, low diversity, populations. Here, we evaluate the success of the Vancouver Island marmot (Marmota vancouverensis; VIM; among the worlds most endangered mammals) captive breeding program in preventing inbreeding and maintaining genetic diversity. We explored the use of high-throughput amplicon sequencing of microsatellite regions to assay greater genetic variation in both captive and wild populations than traditional length-based fragment analysis. Contrary to other studies, this method did not considerably increase diversity estimates, suggesting: (1) that the technique does not universally improve resolution, and (2) VIM have exceedingly low diversity. Studbook estimates of pairwise relatedness and inbreeding in the current population were weakly, but positively, correlated to molecular estimates. Thus, current studbooks are moderately effective at predicting genetic similarity when founder relatedness is known. Finally, we found that captive and wild populations did not differ in allelic frequencies, and conservation efforts to maintain diversity have been successful with no significant decrease in diversity over the last three generations.

Inbreeding Coefficient and Distance in MHC Genes of Parents as Predictors of Reproductive Success in Domestic Cat

Simple Summary

Inbreeding and low diversity in MHC (major histocompatibility complex) genes can have a significant impact on the survival and quality of offspring in mammals. At the same time, a decrease in genetic diversity can be disastrous for animals at individual and species level. For felines, studies of the effects of inbreeding and low variety in MHC genes are conducted on populations with a low number of animals, where there is a high probability of a shortage of available partners, and, accordingly, their choice. The use of model species, especially domestic cats, allows us to identify the main consequences of inbreeding and the lack of a choice of partners for future offspring. The survival of offspring in a domestic cat is primarily affected by the degree of similarity/difference in the genes of the parents’ MHC. Parents with the maximum distance in MHC genes have a larger proportion of surviving kittens, and this effect is most pronounced immediately after birth. In parents with the minimum distance in MHC genes, a significant percentage of kittens are either stillborn or die on the first day after birth. However, inbreeding and the similarity of parents in MHC genes in domestic cats did not affect the body mass of kittens.

Abstract

Inbreeding and low diversity in MHC genes are considered to have a negative effect on reproductive success in animals. This study presents an analysis of the number and body mass of offspring in domestic cat, depending on the inbreeding coefficient and the degree of similarity in MHC genes of class I and II in parents. Inbred partners had a lower number of live kittens at birth than outbred ones. At the same time, the inbreeding coefficient did not affect the litter size and the number of offspring who survived until the period of transition to solid food. The most significant predictor for the number of surviving offspring was the degree of parental similarity in MHC genes: the parents with the maximum distance in MHC genes had more survived kittens. Moreover, this effect was most pronounced immediately after birth. A significant percentage of kittens from parents with a minimum distance in MHC genes were either stillborn or died on the first day after birth. By the age of transition to solid food, this effect is no longer so pronounced. Furthermore, neither the inbreeding coefficient nor the distance in MHC genes of parents had any effect on the body mass of kittens.

Who are you? A framework to identify and report genetic sample mix-ups

Sample mix-ups occur when samples have accidentally been duplicated, mislabelled or swapped. When samples are subsequently genotyped or sequenced, this can lead to individual IDs being incorrectly linked to genetic data, resulting in incorrect or biased research results, or reduced power to detect true biological patterns. We surveyed the community and found that almost 80% of responding researchers have encountered sample mix-ups. However, many recent studies in the field of molecular ecology do not appear to systematically report individual assignment checks as part of their publications. Although checks may be done, lack of consistent reporting means that it is difficult to assess whether sample mix-ups have occurred or been detected. Here, we present an easy-to-follow sample verification framework that can utilise existing metadata, including species, population structure, sex and pedigree information. We demonstrate its application to a dataset representing individuals of a threatened Aotearoa New Zealand bird species, the hihi, genotyped on a 50K SNP array. We detected numerous incorrect genotype-ID associations when comparing observed and genetic sex or comparing to relationships in a verified microsatellite pedigree. The framework proposed here helped to confirm 488 individuals (39%), correct another 20 bird-genotype links, and detect hundreds of incorrect sample IDs, emphasizing the value of routinely checking genetic and genomic datasets for their accuracy. We therefore promote the implementation and reporting of this simple yet effective sample verification framework as a standardized quality control step for studies in the field of molecular ecology.

The genome sequence of the red deer, Cervus elaphus Linnaeus 1758

We present a genome assembly from an individual female Cervus elaphus (the red deer; Chordata; Mammalia; Artiodactyla; Cervidae). The genome sequence is 2,887 megabases in span. The majority of the assembly is scaffolded into 34 chromosomal pseudomolecules, with the X sex chromosome assembled.

Genomic analysis reveals a polygenic architecture of antler morphology in wild red deer (Cervus elaphus)

Sexually selected traits show large variation and rapid evolution across the animal kingdom, yet genetic variation often persists within populations despite apparent directional selection. A key step in solving this long-standing paradox is to determine the genetic architecture of sexually selected traits to understand evolutionary drivers and constraints at the genomic level. Antlers are a form of sexual weaponry in male red deer (Cervus elaphus). On the island of Rum, Scotland, males with larger antlers have increased breeding success, yet there has been no evidence of any response to selection at the genetic level. To try and understand the mechanisms underlying this observation, we investigate the genetic architecture of ten antler traits and their principal components using genomic data from >38,000 SNPs. We estimate the heritabilities and genetic correlations of the antler traits using a genomic relatedness approach. We then use genome-wide association and haplotype-based regional heritability to identify regions of the genome underlying antler morphology, and an empirical Bayes approach to estimate the underlying distributions of allele effect sizes. We show that antler morphology is highly repeatable over an individual's lifetime, heritable and has a polygenic architecture and that almost all antler traits are positively genetically correlated with some loci identified as having pleiotropic effects. Our findings suggest that a large mutational target and genetic covariances among antler traits, in part maintained by pleiotropy, are likely to contribute to the maintenance of genetic variation in antler morphology in this population.

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.

Admixture mapping reveals loci for carcass mass in red deer x sika hybrids in Kintyre, Scotland

We deployed admixture mapping on a sample of 386 deer from a hybrid swarm between native red deer (Cervus elaphus) and introduced Japanese sika (Cervus nippon) sampled in Kintyre, Scotland to search for quantitative trait loci (QTLs) underpinning phenotypic differences between the species. These two species are highly diverged genetically [Fst between pure species, based on 50K single nucleotide polymorphism (SNPs) = 0.532] and phenotypically: pure red have on average twice the carcass mass of pure sika in our sample (38.7 kg vs 19.1 kg). After controlling for sex, age, and population genetic structure, we found 10 autosomal genomic locations with QTL for carcass mass. Effect sizes ranged from 0.191 to 1.839 kg and as expected, in all cases the allele derived from sika conferred lower carcass mass. The sika population was fixed for all small carcass mass alleles, whereas the red deer population was typically polymorphic. GO term analysis of genes lying in the QTL regions are associated with oxygen transport. Although body mass is a likely target of selection, none of the SNPs marking QTL are introgressing faster or slower than expected in either direction.

Description of breed ancestry and genetic health traits in arctic sled dog breeds

Background

This study describes the presence and frequency of health traits among three populations of dogs traditionally used for sledding and explores their ancestry and breed composition as provided by the commercially available Embark dog DNA test. The three populations include the purebred Siberian Husky and the admixed populations of Alaskan sled dogs and Polar Huskies. While the Siberian Husky represents a well-established breed with extensive historical and health data, the Alaskan sled dog is less studied but has been the subject of nutritional, physiological, and genetic studies related to ancestry and performance. In contrast, the Polar Husky is a relatively obscure and rare group of dogs used for arctic exploration with very little-known information. The three populations were compared using Embark results, providing new insight into the health traits circulating within the populations and the potential ancestral linkage of the health traits between the sledding populations. Embark results are based upon 228,588 single-nucleotide polymorphisms (SNPs) spanning the canine genome, characterized using a custom-designed Illumina beadchip array.

Results

Specifically, breed composition was summarized for the two admixed populations with most of the dogs being predominantly categorized as Alaskan husky- type dog or “Supermutt”. Mitochondrial and Y chromosome haplogroups and haplotypes were found with Alaskan sled dogs carrying most of the haplogroups and types found in Siberian and Polar Huskies. Genomic principal component analysis reflected population structure corresponding to breed and substructure within the Alaskan sled dogs related to sprint or distance competition. Genetic markers associated with Alanine Aminotransferase activity, Alaskan Husky Encephalopathy, dilated cardiomyopathy, Collie eye anomaly, degenerative myelopathy, ichthyosis, and factor VII deficiency were identified in the populations of sledding breeds.

Conclusion

These results provide a preliminary description of genetic characteristics found in sledding breeds, improving the understanding and care of working sled dogs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40575-021-00108-z.

Social evolution in mammals

[Figure: see text].

Exercise-Associated Sudden Death in Finnish Standardbred and Coldblooded Trotters - A Case Series With Pedigree Analysis

Exercise-associated sudden deaths (EASDs) are deaths occurring unexpectedly during or immediately after exercise. Sudden cardiac death (SCD) is one cause of EASD. Cardiac arrhythmias caused by genetic variants have been linked to SCD in humans. We hypothesize that genetic variants may be associated with SCD in animals, including horses. Genetic variants are transmitted to offspring and their frequency might increase within a family. Therefore, the frequency of such variants might increase with the inbreeding factor. Higher inbreeding could have a negative impact on racing performance. Pedigree data and career earnings from racehorses diagnosed with SCD between 2002 and 2017 were compared using non-parametric tests with 1) control horses that died due to catastrophic musculoskeletal injuries and 2) horses that raced during the same period without reported problems. Diagnosis of SCD was based on necropsy reports, including macroscopic and microscopic examinations. Death was registered in the study period for 61 horses. Eleven of these horses were excluded due to missing autopsy reports. In 25 cases, the diagnosis remained unknown and death was possibly caused by cardiac arrhythmia, in two cases cardiac disease was identified, in seven cases a rupture of a major vessel had occurred. In addition, 16 horses died or were euthanized due to severe musculoskeletal injuries. No significant differences in inbreeding coefficients or in career earnings were found between the groups or between horses with EASD compared with other horses racing during the same period. The study provides no evidence for increased inbreeding factor in Finnish racehorses with SCD.

Genomic partitioning of inbreeding depression in humans

Across species, offspring of related individuals often exhibit significant reduction in fitness-related traits, known as inbreeding depression (ID), yet the genetic and molecular basis for ID remains elusive. Here, we develop a method to quantify enrichment of ID within specific genomic annotations and apply it to human data. We analyzed the phenomes and genomes of ∼350,000 unrelated participants of the UK Biobank and found, on average of over 11 traits, significant enrichment of ID within genomic regions with high recombination rates (>21-fold; p < 10^-5), with conserved function across species (>19-fold; p < 10^-4), and within regulatory elements such as DNase I hypersensitive sites (∼5-fold; p = 8.9 × 10^-7). We also quantified enrichment of ID within trait-associated regions and found suggestive evidence that genomic regions contributing to additive genetic variance in the population are enriched for ID signal. We find strong correlations between functional enrichment of SNP-based heritability and that of ID (r = 0.8, standard error: 0.1). These findings provide empirical evidence that ID is most likely due to many partially recessive deleterious alleles in low linkage disequilibrium regions of the genome. Our study suggests that functional characterization of ID may further elucidate the genetic architectures and biological mechanisms underlying complex traits and diseases.

Inbreeding reduces fitness of seed beetles under thermal stress

Human‐induced environmental change can influence populations both at the global level through climatic warming and at the local level through habitat fragmentation. As populations become more isolated, they can suffer from high levels of inbreeding, which contributes to a reduction in fitness, termed inbreeding depression. However, it is still unclear if this increase in homozygosity also results in a corresponding increase in sensitivity to stressful conditions, which could intensify the already detrimental effects of environmental warming. Here, in a fully factorial design, we assessed the life‐long impact of increased inbreeding load and elevated temperature on key life history traits in the seed beetle, Callosobruchus maculatus. We found that beetles raised at higher temperatures had far reduced fitness and survival than beetles from control temperatures. Importantly, these negative effects were exacerbated in inbred beetles as a result of increased inbreeding load, with further detrimental effects manifesting on individual eclosion probability and lifetime reproductive success. These results reveal the harmful impact that increasing temperature and likelihood of habitat fragmentation due to anthropogenetic changes in environmental conditions could have on populations of organisms worldwide.

The role of maternally transferred antibodies in maternal performance in red deer

Maternal effects are ubiquitous. Yet, the pathways through which maternal effects occur in wild mammals remain largely unknown. We hypothesise that maternal immune transfer is a key mechanism by which mothers can affect their offspring fitness, and that individual variation in maternally derived antibodies mainly depends on a mother's characteristics and the environmental conditions she experiences. To test this, we assayed six colostrum-derived antibodies in the plasma of 1447 neonates in a wild red deer population. Neonatal antibody levels were mainly affected by maternal genes, environmental variation and costs of prior reproductive investment. We found consistent heterogeneity in maternal performance across traits, with mothers producing the heaviest calves also having calves with more antibodies. Unexpectedly, antibody levels were not associated with calf survival. We provide a unique example of how evolutionary theory on maternal effects can be used to gain insight into the causes of maternal effects in wild populations.

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.

The long-standing significance of genetic diversity in conservation

Since allozymes were first used to assess genetic diversity in the 1960s and 1970s, biologists have attempted to characterize gene pools and conserve the diversity observed in domestic crops, livestock, zoos and (more recently) natural populations. Recently, some authors have claimed that the importance of genetic diversity in conservation biology has been greatly overstated. Here, we argue that a voluminous literature indicates otherwise. We address four main points made by detractors of genetic diversity's role in conservation by using published literature to firmly establish that genetic diversity is intimately tied to evolutionary fitness, and that the associated demographic consequences are of paramount importance to many conservation efforts. We think that responsible management in the Anthropocene should, whenever possible, include the conservation of ecosystems, communities, populations and individuals, and their underlying genetic diversity.

Mutation load decreases with haplotype age in wild Soay sheep

Runs of homozygosity (ROH) are pervasive in diploid genomes and expose the effects of deleterious recessive mutations, but how exactly these regions contribute to variation in fitness remains unclear. Here, we combined empirical analyses and simulations to explore the deleterious effects of ROH with varying genetic map lengths in wild Soay sheep. Using a long-term dataset of 4879 individuals genotyped at 417K SNPs, we found that inbreeding depression increases with ROH length. A 1% genomic increase in long ROH (>12.5 cM) reduced the odds of first-year survival by 12.4% compared to only 7.7% for medium ROH (1.56-12.5 cM), whereas short ROH (<1.56 cM) had no effect on survival. We show by forward genetic simulations that this is predicted: compared to shorter ROH, long ROH will have higher densities of deleterious alleles, with larger average effects on fitness and lower population frequencies. Taken together, our results are consistent with the idea that the mutation load decreases in older haplotypes underlying shorter ROH, where purifying selection has had more time to purge deleterious mutations. Finally, our study demonstrates that strong inbreeding depression can persist despite ongoing purging in a historically small population.

Genomic and fitness consequences of inbreeding in an endangered carnivore

Reduced fitness through genetic drift and inbreeding is a major threat to small and isolated populations. Although previous studies have generally used genetically verified pedigrees to document effects of inbreeding and gene flow, these often fail to capture the whole inbreeding history of the species. By assembling a draft arctic fox (Vulpes lagopus) genome and resequencing complete genomes of 23 additional foxes born before and after a well-documented immigration event in Scandinavia, we here look into the genomic consequences of inbreeding and genetic rescue. We found a difference in genome-wide diversity, with 18% higher heterozygosity and 81% lower F_ROH in immigrant F1 compared to native individuals. However, more distant descendants of immigrants (F2, F3) did not show the same pattern. We also found that foxes with lower inbreeding had higher probability to survive their first year of life. Our results demonstrate the important link between genetic variation and fitness as well as the transient nature of genetic rescue. Moreover, our results have implications in conservation biology as they demonstrate that inbreeding depression can effectively be detected in the wild by a genomic approach.

Genetic architecture and lifetime dynamics of inbreeding depression in a wild mammal

Inbreeding depression is ubiquitous, but we still know little about its genetic architecture and precise effects in wild populations. Here, we combine long-term life-history data with 417 K imputed SNP genotypes for 5952 wild Soay sheep to explore inbreeding depression on a key fitness component, annual survival. Inbreeding manifests in long runs of homozygosity (ROH), which make up nearly half of the genome in the most inbred individuals. The ROH landscape varies widely across the genome, with islands where up to 87% and deserts where only 4% of individuals have ROH. The fitness consequences of inbreeding are severe; a 10% increase in individual inbreeding FROH is associated with a 60% reduction in the odds of survival in lambs, though inbreeding depression decreases with age. Finally, a genome-wide association scan on ROH shows that many loci with small effects and five loci with larger effects contribute to inbreeding depression in survival.

Demographic History, Not Mating System, Explains Signatures of Inbreeding and Inbreeding Depression in a Large Outbred Population

AbstractInbreeding depression is often found in small, inbred populations, but whether it can be detected in and have evolutionary consequences for large, wide-ranging populations is poorly known. Here, we investigate the possibility of inbreeding in a large population to determine whether mild levels of inbreeding can still have genetic and phenotypic consequences and how genomically widespread these effects can be. We apply genome-wide methods to investigate whether individual and parental heterozygosity is related to morphological, growth, or life-history traits in a pelagic seabird, Leach's storm-petrel (Oceanodroma leucorhoa). Examining 560 individuals as part of a multiyear study, we found a substantial effect of maternal heterozygosity on chick traits: chicks from less heterozygous (relatively inbred) mothers were significantly smaller than chicks from more heterozygous (noninbred) mothers. We show that these heterozygosity-fitness correlations were due to general genome-wide effects and demonstrate a correlation between heterozygosity and inbreeding, suggesting inbreeding depression. We used population genetic models to further show that the variance in inbreeding was probably due to past demographic events rather than the current mating system and ongoing mate choice. Our findings demonstrate that inbreeding depression can be observed in large populations and illustrate how the integration of genomic techniques and fieldwork can elucidate its underlying causes.