A comparison of pedigree, genetic and genomic estimates of relatedness for informing pairing decisions in two critically endangered birds: Implications for conservation breeding programmes worldwide

Epigenetic Diversity and the Evolutionary Potential of Wild Populations

Fast-paced selective pressures imposed by climate change and anthropogenic activities call for adaptive evolutionary responses to emerge at ecological timescales. However, the evolution and heritability of genomic variation underlie mechanistic constraints, which dictate a slower pace of adaptation exclusively relying on standing genetic variation and novel mutations. Environmentally responsive epigenetic mechanisms can allow acclimatisation and adaptive phenotypes to arise faster than DNA sequence-based mechanisms alone. Nevertheless, the knowledge gap between identifying epigenetic marks and effectively deeming them functional is still wide in a natural context and often outside the scope of model organisms. With this Special Issue, we aimed to narrow this gap by presenting a compilation of original research articles, reviews and opinions on the topic of epigenetics in wild populations. We contextualised this collection within the overarching topic of conservation biology, as we firmly propose that epigenetic research can significantly enhance the effectiveness of conservation measures. Contributions highlighted the putative role of epigenetic variation in the acclimatisation and adaptive potential of species and populations directly and indirectly affected by climatic shifts and anthropogenic actions. They further exemplified how epigenetic variation can be used as biomarkers for monitoring variations in physiology, phenology and behaviour. Lastly, reviews and perspective articles illustrated the past and present of epigenetic research in wild populations while suggesting future research avenues.

Analysing the pedigree to identify undesirable losses of genetic diversity and to prioritize management decisions in captive breeding: a case study

When prevention of species extinction is the priority, captive breeding is a key component in conservation programmes, allowing the recording of pedigree information in studbooks. The genealogical information registered in Cuvier's gazelle studbook between 1975 and 2023 was analysed to (a) assess if the implemented mating policy was successful in preserving the genetic background of the founders (1 male:3 females) in the present population, and b) improve future management and breeding decisions. Although the maternal contribution of one founder female was lost and the mean inbreeding of the total live population was high (0.305 ± 0.095), the breeding policy applied produced better results than expected from a population starting from four founders. It was successful in keeping the individual increase in inbreeding low (0.047 ± 0.021), and, notably, the inbreeding tended to decrease during the last three decades of the breeding programme, ensuring the viability of this highly inbred population. Historical dissemination of individuals among the zoos of Europe and North America caused population structuring and genetic differentiation of the live North American population. However, it did not risk the viability of the captive population. The average relatedness coefficients allowed the identification of individuals with underrepresented genotypes, which is relevant to plan future mating guidelines to keep the founders' representation balanced in the next generations. This study highlights the importance of keeping long-term pedigree information to monitor changes in the genetic diversity of captive populations, which is crucial to implement optimal mating decisions and assuring their long-term viability within an ex situ conservation programme.

Benchmarking kinship estimation tools for ancient genomes using pedigree simulations

There is growing interest in uncovering genetic kinship patterns in past societies using low-coverage palaeogenomes. Here, we benchmark four tools for kinship estimation with such data: lcMLkin, NgsRelate, KIN, and READ, which differ in their input, IBD estimation methods, and statistical approaches. We used pedigree and ancient genome sequence simulations to evaluate these tools when only a limited number (1 to 50 K, with minor allele frequency ≥0.01) of shared SNPs are available. The performance of all four tools was comparable using ≥20 K SNPs. We found that first-degree related pairs can be accurately classified even with 1 K SNPs, with 85% F1 scores using READ and 96% using NgsRelate or lcMLkin. Distinguishing third-degree relatives from unrelated pairs or second-degree relatives was also possible with high accuracy (F1 > 90%) with 5 K SNPs using NgsRelate and lcMLkin, while READ and KIN showed lower success (69 and 79% respectively). Meanwhile, noise in population allele frequencies and inbreeding (first-cousin mating) led to deviations in kinship coefficients, with different sensitivities across tools. We conclude that using multiple tools in parallel might be an effective approach to achieve robust estimates on ultra-low-coverage genomes.

Estimates of heterozygosity from single nucleotide polymorphism markers are context-dependent and often wrong

Genetic diversity is frequently described using heterozygosity, particularly in a conservation context. Often, it is estimated using single nucleotide polymorphisms (SNPs); however, it has been shown that heterozygosity values calculated from SNPs can be biased by both study design and filtering parameters. Though solutions have been proposed to address these issues, our own work has found them to be inadequate in some circumstances. Here, we aimed to improve the reliability and comparability of heterozygosity estimates, specifically by investigating how sample size and missing data thresholds influenced the calculation of autosomal heterozygosity (heterozygosity calculated from across the genome, i.e. fixed and variable sites). We also explored how the standard practice of tri- and tetra-allelic site exclusion could bias heterozygosity estimates and influence eventual conclusions relating to genetic diversity. Across three distinct taxa (a frog, Litoria rubella; a tree, Eucalyptus microcarpa; and a grasshopper, Keyacris scurra), we found heterozygosity estimates to be meaningfully affected by sample size and missing data thresholds, partly due to the exclusion of tri- and tetra-allelic sites. These biases were inconsistent both between species and populations, with more diverse populations tending to have their estimates more severely affected, thus having potential to dramatically alter interpretations of genetic diversity. We propose a modified framework for calculating heterozygosity that reduces bias and improves the utility of heterozygosity as a measure of genetic diversity, whilst also highlighting the need for existing population genetic pipelines to be adjusted such that tri- and tetra-allelic sites be included in calculations.

Reduction of genetic diversity in 'Alalā (Hawaiian crow; Corvus hawaiiensis) between the late 1800s and the late 1900s

Genetic and genomic data are increasingly used to aid conservation management of endangered species by providing insights into evolutionary histories, factors associated with extinction risks, and potential for future adaptation. For the 'Alalā, or Hawaiian crow (Corvus hawaiiensis), genetic concerns include negative correlations between inbreeding and hatching success. However, it is unclear if low genetic diversity and inbreeding depression are consequences of a historical population bottleneck, or if 'Alalā had historically low genetic diversity that predated human influence, perhaps as a result of earlier declines or founding events. In this study, we applied a hybridization-based sequence capture to generate a genome-wide single nucleotide polymorphism (SNP) dataset for comparing historical specimens collected in the 1890s, when 'Alalā were more numerous, to samples taken between 1973 and 1998, when 'Alalā population densities were near the lowest documented levels in the wild, prior to all individuals being collected for captive rearing. We found low genome-wide diversity in both sample groups, however, the modern sample group (1973 to 1998 cohort) exhibited relatively fewer polymorphic alleles, a lower proportion of polymorphic loci, and lower observed heterozygosity, consistent with a population decline and potential bottleneck effects. These results combined with a current low population size highlight the importance of continued efforts by conservation managers to mitigate inbreeding and maintain founder representation to preserve what genetic diversity remains.

Pedigree analysis in the mhorr gazelle (Nanger dama mhorr): Genetic variability evolution of the captive population

Breeding programs have an essential role in the recovery of threatened populations through optimal genetic management and mating strategies. The dama gazelle (Nanger dama) is a North African ungulate listed as critically endangered. The mhorr subspecies is extinct in the wild and currently survives thanks to the creation in 1971 of an ex situ breeding program. The aim of the present study was to assess the evolution of genetic variability in this mhorr gazelle captive population, as well as the mating strategy used in two reference populations studied (Almeria and Europe). The entire pedigree, with 2739 animals, was analyzed to measure demographic characters, pedigree completeness level, probability of gene origin, level of relatedness and genetic structure of the population. The population size has been progressively increasing, with up to 264 individuals alive in Europe at the time of the study. The average number of equivalent complete generations was 5.55. The effective number of founders and ancestors was both 3, and the founder genome equivalent was 1.99. The genetic contributions of the four main ancestors were unbalanced. The average values of inbreeding and average relatedness for the whole pedigree were, respectively, 28.34% and 50.14%. The effective population size was 8.7 by individual increase in inbreeding and 9.8 by individual increase in coancestry. F-statistics evidenced a very small level of population subdivision (F ST = 0.033370). The mating strategy used, based on the minimum coancestry of the individuals, has minimized the losses of genetic variability and helped to balance the genetic contributions between ancestors. The strategy also avoided large subdivisions within the population and the appearance of new bottlenecks. This study shows how pedigree analysis can both be used to determine the genetic variability of the population and to assess the influence of the mating strategy used in the breeding program on such variability.

Female chimpanzees avoid inbreeding even in the presence of substantial bisexual philopatry

Inbreeding (reproduction between relatives) often decreases the fitness of offspring and is thus expected to lead to the evolution of inbreeding avoidance strategies. Chimpanzees (Pan troglodytes) are expected to avoid inbreeding as they are long-lived, invest heavily in offspring and may encounter adult, opposite sex kin frequently, especially in populations where both males and females commonly remain in the group in which they were born (bisexual philopatry). However, it is unclear whether substantial bisexual philopatry has been a feature of chimpanzees' evolutionary history or whether it is a result of recent anthropogenic interference, as the only groups for which it has been documented are significantly impacted by human encroachment and experience notable rates of potentially unsustainable inbreeding. Here we use 14 years of observational data and a large genomic dataset of 256 481 loci sequenced from 459 individuals to document dispersal and inbreeding dynamics in an eastern chimpanzee (P. t. schweinfurthii) community with low levels of anthropogenic disturbance. We document the first case of substantial bisexual philopatry in a relatively undisturbed chimpanzee community and show that, despite an increased inbreeding risk incurred by females who do not disperse before reaching reproductive age, natal females were still able to avoid producing inbred offspring.

Quantifying uncertainty in inferences of landscape genetic resistance due to choice of individual-based genetic distance metric

Estimates of gene flow resulting from landscape resistance inferences frequently inform conservation management decision-making processes. Therefore, results must be robust across approaches and reflect real-world gene flow instead of methodological artefacts. Here, we tested the impact of 32 individual-based genetic distance metrics on the robustness and accuracy of landscape resistance modelling results. We analysed three empirical microsatellite datasets and 36 simulated datasets that varied in landscape resistance and genetic spatial autocorrelation. We used ResistanceGA to generate optimised multi-feature resistance surfaces for each of these datasets using 32 different genetic distance metrics. Results of the empirical dataset demonstrated that the choice of genetic distance metric can have strong impacts on inferred optimised resistance surfaces. Simulations showed accurate parametrisation of resistance surfaces across most genetic distance metrics only when a small number of environmental features was impacting gene flow. Landscape scenarios with many features impacting gene flow led to a generally poor recovery of true resistance surfaces. Simulation results also emphasise that choosing a genetic distance metric should not be based on marginal R2 -based model fit. Until more robust methods are available, resistance surfaces can be optimised with different genetic distance metrics and the convergence of results needs to be assessed via pairwise matrix correlations. Based on the results presented here, high correlation coefficients across different genetic distance categories likely indicate accurate inference of true landscape resistance. Most importantly, empirical results should be interpreted with great caution, especially when they appear counter-intuitive in light of the ecology of a species.

Genomic approaches to mitigating genetic diversity loss in declining populations

The accelerating pace of global biodiversity loss is exacerbated by habitat fragmentation and subsequent inbreeding in small populations. To address this problem, conservation practitioners often turn to assisted breeding programmes with the aim of enhancing genetic diversity in declining populations. Although genomic information is infrequently included in these efforts, it has the potential to significantly enhance the success of such programmes. In this study, we showcase the value of genomic approaches for increasing genetic diversity in assisted breeding efforts, specifically focusing on a highly inbred population of Western burrowing owls. To maximize genetic diversity in the resulting offspring, we begin by creating an optimal pairing decision tree based on sex, kinship and patterns of homozygosity across the genome. To evaluate the effectiveness of our strategy, we compare genetic diversity, brood size and nestling success rates between optimized and non-optimized pairs. Additionally, we leverage recently discovered correlations between telomere length and fitness across species to investigate whether genomic optimization could have long-term fitness benefits. Our results indicate that pairing individuals with contrasting patterns of homozygosity across the genome is an effective way to increase genetic diversity in offspring. Although short-term field-based metrics of success did not differ significantly between optimized and non-optimized pairs, offspring from optimized pairs had significantly longer telomeres, suggesting that genetic optimization can help reduce the risk of inbreeding depression. These findings underscore the importance of genomic tools for informing efforts to preserve the adaptive potential of small, inbred populations at risk of further decline.

Species-wide genomics of kākāpō provides tools to accelerate recovery

The kākāpō is a critically endangered, intensively managed, long-lived nocturnal parrot endemic to Aotearoa New Zealand. We generated and analysed whole-genome sequence data for nearly all individuals living in early 2018 (169 individuals) to generate a high-quality species-wide genetic variant callset. We leverage extensive long-term metadata to quantify genome-wide diversity of the species over time and present new approaches using probabilistic programming, combined with a phenotype dataset spanning five decades, to disentangle phenotypic variance into environmental and genetic effects while quantifying uncertainty in small populations. We find associations for growth, disease susceptibility, clutch size and egg fertility within genic regions previously shown to influence these traits in other species. Finally, we generate breeding values to predict phenotype and illustrate that active management over the past 45 years has maintained both genome-wide diversity and diversity in breeding values and, hence, evolutionary potential. We provide new pathways for informing future conservation management decisions for kākāpō, including prioritizing individuals for translocation and monitoring individuals with poor growth or high disease risk. Overall, by explicitly addressing the challenge of the small sample size, we provide a template for the inclusion of genomic data that will be transformational for species recovery efforts around the globe.

Analyzing captive breeding outcomes to inform reintroduction practice: lessons from the pookila (Pseudomys novaehollandiae)

Captive breeding is often used to produce individuals for reintroduction programs in order to reestablish a species in an area where it has become locally extinct. To maximize the likelihood of establishing a self-sustaining population in the wild, an analysis of data from captive breeding programs is commonly undertaken to (1) increase the quantity of individuals and rate at which they can be released, and (2) maintain or improve the genetic and phenotypic quality of individuals. Here we demonstrate how the knowledge gained from these analyses can also be applied to decision-making during the design of subsequent reintroductions to further advance a reintroduction program toward success. We conducted an analysis of data from a captive breeding program for the threatened pookila (Pseudomys novaehollandiae, New Holland mouse) spanning 6 years. We found evidence for relationships between the reproductive output of pookila and behavioral, demographic, experiential, health, and physiological predictors. Based on a biological interpretation of these results, and with reference to a checklist of all known translocation tactics, we recommend 11 specific design elements to maximize the probability of pookila reproduction postrelease (thereby improving the likelihood of reintroduction success). These recommendations should be interpreted as hypotheses to be evaluated and refined in future reintroduction trials for the pookila. The uncertainty around the postrelease survival and reproduction of a species that is common in reintroduction practice warrants the creative use of existing data to inform adaptive management. Indeed, there is a wealth information in well-kept captive breeding records that is currently underused by reintroduction practitioners. The direct integration of knowledge derived from captive breeding (where available) with decision-making for reintroductions, as described here, will help navigate these uncertainties, which would benefit the conservation of both understudied and well-known species around the world.

Fish germ cell cryobanking and transplanting for conservation

The unprecedented loss of global biodiversity is linked to multiple anthropogenic stressors. New conservation technologies are urgently needed to mitigate this loss. The rights, knowledge and perspectives of Indigenous peoples in biodiversity conservation-including the development and application of new technologies-are increasingly recognised. Advances in germplasm cryopreservation and germ cell transplantation (termed 'broodstock surrogacy') techniques offer exciting tools to preserve biodiversity, but their application has been underappreciated. Here, we use teleost fishes as an exemplar group to outline (1) the power of these techniques to preserve genome-wide genetic diversity, (2) the need to apply a conservation genomic lens when selecting individuals for germplasm cryobanking and broodstock surrogacy and (3) the value of considering the cultural significance of these genomic resources. We conclude by discussing the opportunities and challenges of these techniques for conserving biodiversity in threatened teleost fish and beyond.

Genomics of founders for conservation breeding: the Jasper caribou case

Conservation breeding programs are increasingly used as recovery actions for wild animals; bringing founders into captivity to rear captive populations for future reintroduction into the wild. The International Union for the Conservation of Nature recommends that founders should come from genetically close populations and should have sufficient genetic diversity to avoid mating among relatives. Genomic data are highly informative for evaluating founders due to their high resolution and ability to capture adaptive divergence, yet, their application in that context remains limited. Woodland caribou are federally listed as a Species at Risk in Canada, with several populations facing extirpation, such as those in the Rocky Mountains of Alberta and British Columbia (BC). To prevent local extirpation, Jasper National Park (JNP) is proposing a conservation breeding program. We examined single nucleotide polymorphisms for 144 caribou from 11 populations encompassing a 200,0002 km area surrounding JNP to provide information useful for identifying appropriate founders for this program. We found that this area likely hosts a caribou metapopulation historically characterized by high levels of gene flow, which indicates that multiple sources of founders would be appropriate for initiating a breeding program. However, population structure and adaptive divergence analyses indicate that JNP caribou are closest to populations in the BC Columbia range, which also have suitable genetic diversity for conservation breeding. We suggest that collaboration among jurisdictions would be beneficial to implement the program to promote recovery of JNP caribou and possibly other caribou populations in the surrounding area, which is strategically at the periphery of the distribution of this endangered species.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10592-023-01540-3.

Hidden impacts of conservation management on fertility of the critically endangered kākāpō

Background

Animal conservation often requires intensive management actions to improve reproductive output, yet any adverse effects of these may not be immediately apparent, particularly in threatened species with small populations and long lifespans. Hand-rearing is an example of a conservation management strategy which, while boosting populations, can cause long-term demographic and behavioural problems. It is used in the recovery of the critically endangered kākāpō (Strigops habroptilus), a flightless parrot endemic to New Zealand, to improve the slow population growth that is due to infrequent breeding, low fertility and low hatching success.

Methods

We applied Bayesian mixed models to examine whether hand-rearing and other factors were associated with clutch fertility in kākāpō. We used projection predictive variable selection to compare the relative contributions to fertility from the parents' rearing environment, their age and previous copulation experience, the parental kinship, and the number of mates and copulations for each clutch. We also explored how the incidence of repeated copulations and multiple mates varied with kākāpō density.

Results

The rearing status of the clutch father and the number of mates and copulations of the clutch mother were the dominant factors in predicting fertility. Clutches were less likely to be fertile if the father was hand-reared compared to wild-reared, but there was no similar effect for mothers. Clutches produced by females copulating with different males were more likely to be fertile than those from repeated copulations with one male, which in turn had a higher probability of fertility than those from a single copulation. The likelihood of multiple copulations and mates increased with female:male adult sex ratio, perhaps as a result of mate guarding by females. Parental kinship, copulation experience and age all had negligible associations with clutch fertility.

Conclusions

These results provide a rare assessment of factors affecting fertility in a wild threatened bird species, with implications for conservation management. The increased fertility due to multiple mates and copulations, combined with the evidence for mate guarding and previous results of kākāpō sperm morphology, suggests that an evolutionary mechanism exists to optimise fertility through sperm competition in kākāpō. The high frequency of clutches produced from single copulations in the contemporary population may therefore represent an unnatural state, perhaps due to too few females. This suggests that opportunity for sperm competition should be maximised by increasing population densities, optimising sex ratios, and using artificial insemination. The lower fertility of hand-reared males may result from behavioural defects due to lack of exposure to conspecifics at critical development stages, as seen in other taxa. This potential negative impact of hand-rearing must be balanced against the short-term benefits it provides.

Sex ratio and relatedness in the Griffon vulture (Gyps fulvus) population of Serbia

Background

Once a widespread species across the region of Southeast Europe, the Griffon vulture is now confined to small and isolated populations across the Balkan Peninsula. The population from Serbia represents its biggest and most viable population that can serve as an important reservoir of genetic diversity from which the birds can be used for the region's reintroduction programmes. The available genetic data for this valuable population are scarce and as a protected species that belongs to the highly endangered vulture group, it needs to be well described so that it can be properly managed and used as a restocking population. Considering the serious recent bottleneck event that the Griffon vulture population from Serbia experienced we estimated the overall relatedness among the birds from this population. Sex ratio, another important parameter that shows the vitality and strength of the population was evaluated as well.

Methods

During the annual monitoring that was performed in the period from 2013-2021, we collected blood samples from individual birds that were marked in the nests. In total, 169 samples were collected and each was used for molecular sexing while 58 presumably unrelated birds from different nests were used for inbreeding and relatedness analyses. The relatedness was estimated using both biparentally (10 microsatellite loci) and uniparentally (Cytb and D-loop I of mitochondrial DNA) inherited markers.

Results

The level of inbreeding was relatively high and on average it was 8.3% while the mean number of relatives for each bird was close to three. The sex ratio was close to 1:1 and for the analysed period of 9 years, it didn't demonstrate a statistically significant deviation from the expected ratio of 1:1, suggesting that this is a stable and healthy population. Our data suggest that, even though a relatively high level of inbreeding can be detected among the individual birds, the Griffon vulture population from Serbia can be used as a source population for restocking and reintroduction programmes in the region. These data combined with previously observed genetic differentiation between the populations from the Iberian and Balkan Peninsulas suggest that the introduction of foreign birds should be avoided and that local birds should be used instead.

Leveraging an existing whole-genome resequencing population data set to characterize toll-like receptor gene diversity in a threatened bird

Species recovery programs are increasingly using genomic data to measure neutral genetic diversity and calculate metrics like relatedness. While these measures can inform conservation management, determining the mechanisms underlying inbreeding depression requires information about functional genes associated with adaptive or maladaptive traits. Toll-like receptors (TLRs) are one family of functional genes, which play a crucial role in recognition of pathogens and activation of the immune system. Previously, these genes have been analysed using species-specific primers and PCR. Here, we leverage an existing short-read reference genome, whole-genome resequencing population data set, and bioinformatic tools to characterize TLR gene diversity in captive and wild tchūriwat'/tūturuatu/shore plover (Thinornis novaeseelandiae), a threatened bird endemic to Aotearoa New Zealand. Our results show that TLR gene diversity in tchūriwat'/tūturuatu is low, and forms two distinct captive and wild genetic clusters. The bioinformatic approach presented here has broad applicability to other threatened species with existing genomic resources in Aotearoa New Zealand and beyond.

Comparing genome-based estimates of relatedness for use in pedigree-based conservation management

Researchers have long debated which estimator of relatedness best captures the degree of relationship between two individuals. In the genomics era, this debate continues, with relatedness estimates being sensitive to the methods used to generate markers, marker quality, and levels of diversity in sampled individuals. Here, we compare six commonly used genome-based relatedness estimators (kinship genetic distance (KGD), Wang Maximum Likelihood (TrioML), Queller and Goodnight (R_xy ), Kinship INference for Genome-wide association studies (KING-robust), and Pairwise Relatedness (R_AB ), allele-sharing co-ancestry (AS)) across five species bred in captivity-including three birds and two mammals-with varying degrees of reliable pedigree data, using reduced-representation and whole genome resequencing data. Genome-based relatedness estimates varied widely across estimators, sequencing methods, and species, yet the most consistent results for known first order relationships were found using R_xy , R_AB , and AS. However, AS was found to be less consistently correlated with known pedigree relatedness than either R_xy or R_AB . Our combined results indicate there is not a single genome-based estimator that is ideal across different species and data types. To determine the most appropriate genome-based relatedness estimator for each new dataset, we recommend assessing the relative: (1) correlation of candidate estimators with known relationships in the pedigree and (2) precision of candidate estimators with known first-order relationships. These recommendations are broadly applicable to conservation breeding programs, particularly where genome-based estimates of relatedness can complement and complete poorly pedigreed populations. Given a growing interest in the application of wild pedigrees, our results are also applicable to in-situ wildlife management.

Insights from the rescue and breeding management of Cuvier's gazelle (Gazella cuvieri) through whole-genome sequencing

Captive breeding programmes represent the most intensive type of ex situ population management for threatened species. One example is the Cuvier's gazelle programme that started in 1975 with only four founding individuals, and after more than four decades of management in captivity, a reintroduction effort was undertaken in Tunisia in 2016, to establish a population in an area historically included within its range. Here, we aim to determine the genetic consequences of this reintroduction event by assessing the genetic diversity of the founder stock as well as of their descendants. We present the first whole-genome sequencing dataset of 30 Cuvier's gazelles including captive-bred animals, animals born in Tunisia after a reintroduction and individuals from a genetically unrelated Moroccan population. Our analyses revealed no difference between the founder and the offspring cohorts in genome-wide heterozygosity and inbreeding levels, and in the amount and length of runs of homozygosity. The captive but unmanaged Moroccan gazelles have the lowest genetic diversity of all genomes analysed. Our findings demonstrate that the Cuvier's gazelle captive breeding programme can serve as source populations for future reintroductions of this species. We believe that this study can serve as a starting point for global applications of genomics to the conservation plan of this species.

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.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10592-022-01429-7.

Genetic approaches for increasing fitness in endangered species

The global rate of wildlife extinctions is accelerating, and the persistence of many species requires conservation breeding programs. A central paradigm of these programs is to preserve the genetic diversity of the founder populations. However, this may preserve original characteristics that make them vulnerable to extinction. We introduce targeted genetic intervention (TGI) as an alternative approach that promotes traits that enable species to persist in the face of threats by changing the incidence of alleles that impact on fitness. The TGI toolkit includes methods with established efficacy in model organisms and agriculture but are largely untried for conservation, such as synthetic biology and artificial selection. We explore TGI approaches as a species-restoration tool for intractable threats including infectious disease and climate change.

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.

Captive breeding and the conservation of the threatened houbara bustards

Translocation of captive-bred individuals to reinforce wild populations may be an important conservation approach for some species, but can be detrimental when employed to boost exploited wild populations, particularly where repeated long-term reinforcement aims to compensate for repeated unregulated offtake. We review evidence that captive breeding alters multiple physiological, life-history and temperamental traits through founder effects, genetic drift and unintended adaption to captivity; degrades learnt behaviours; and compromises biogeography, population structure and viability through introgression. We highlight these risks for the globally threatened African houbara Chlamydotis undulata and Asian houbara C. macqueenii, 2 bustard species hunted throughout much of their ranges and now subject to multiple large-scale captive-breeding programmes and translocations. In eastern Morocco, annual releases of captive-bred African houbara are 2‒3 times higher than original wild numbers, but no investigation of their potentially deleterious effects has, to our knowledge, been published, although most wild populations may now have been replaced by captive-bred domestic stock, which are reportedly not self-sustaining. Despite multiple decades of reinforcement, we are not aware of any analysis of the contribution of captive breeding to African houbara population dynamics, or of the genomic consequences. Asian houbara release programmes may also be promoting rather than preventing declines, and need to contextualise themselves through rigorous analyses of wild population numbers, demographic rates and threats, maintenance of phylogeographic concordance of released with supplemented populations, profiling of traits crucial to survival and the measurement and modelling of the impacts of reinforcement on physiological and behavioural fitness of wild populations.

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.

The relevance of pedigrees in the conservation genomics era

Over the past 50 years conservation genetics has developed a substantive toolbox to inform species management. One of the most long‐standing tools available to manage genetics—the pedigree—has been widely used to characterize diversity and maximize evolutionary potential in threatened populations. Now, with the ability to use high throughput sequencing to estimate relatedness, inbreeding, and genome‐wide functional diversity, some have asked whether it is warranted for conservation biologists to continue collecting and collating pedigrees for species management. In this perspective, we argue that pedigrees remain a relevant tool, and when combined with genomic data, create an invaluable resource for conservation genomic management. Genomic data can address pedigree pitfalls (e.g., founder relatedness, missing data, uncertainty), and in return robust pedigrees allow for more nuanced research design, including well‐informed sampling strategies and quantitative analyses (e.g., heritability, linkage) to better inform genomic inquiry. We further contend that building and maintaining pedigrees provides an opportunity to strengthen trusted relationships among conservation researchers, practitioners, Indigenous Peoples, and Local Communities.

Are immigrants outbred and unrelated? Testing standard assumptions in a wild metapopulation

Immigration into small recipient populations is expected to alleviate inbreeding and increase genetic variation, and hence facilitate population persistence through genetic and/or evolutionary rescue. Such expectations depend on three standard assumptions: that immigrants are outbred, unrelated to existing natives at arrival, and unrelated to each other. These assumptions are rarely explicitly verified, including in key field systems in evolutionary ecology. Yet, they could be violated due to non-random or repeated immigration from adjacent small populations. We combined molecular genetic marker data for 150-160 microsatellite loci with comprehensive pedigree data to test the three assumptions for a song sparrow (Melospiza melodia) population that is a model system for quantifying effects of inbreeding and immigration in the wild. Immigrants were less homozygous than existing natives on average, with mean homozygosity that closely resembled outbred natives. Immigrants can therefore be considered outbred on the focal population scale. Comparisons of homozygosity of real or hypothetical offspring of immigrant-native, native-native and immigrant-immigrant pairings implied that immigrants were typically unrelated to existing natives and to each other. Indeed, immigrants' offspring would be even less homozygous than outbred individuals on the focal population scale. The three standard assumptions of population genetic and evolutionary theory were consequently largely validated. Yet, our analyses revealed some deviations that should be accounted for in future analyses of heterosis and inbreeding depression, implying that the three assumptions should be verified in other systems to probe patterns of non-random or repeated dispersal and facilitate precise and unbiased estimation of key evolutionary parameters.

Exploiting genomic synteny in Felidae: cross-species genome alignments and SNV discovery can aid conservation management

Background

While recent advances in genomics has enabled vast improvements in the quantification of genome-wide diversity and the identification of adaptive and deleterious alleles in model species, wildlife and non-model species have largely not reaped the same benefits. This has been attributed to the resources and infrastructure required to develop essential genomic datasets such as reference genomes. In the absence of a high-quality reference genome, cross-species alignments can provide reliable, cost-effective methods for single nucleotide variant (SNV) discovery. Here, we demonstrated the utility of cross-species genome alignment methods in gaining insights into population structure and functional genomic features in cheetah (Acinonyx jubatas), snow leopard (Panthera uncia) and Sumatran tiger (Panthera tigris sumatrae), relative to the domestic cat (Felis catus).

Results

Alignment of big cats to the domestic cat reference assembly yielded nearly complete sequence coverage of the reference genome. From this, 38,839,061 variants in cheetah, 15,504,143 in snow leopard and 13,414,953 in Sumatran tiger were discovered and annotated. This method was able to delineate population structure but limited in its ability to adequately detect rare variants. Enrichment analysis of fixed and species-specific SNVs revealed insights into adaptive traits, evolutionary history and the pathogenesis of heritable diseases.

Conclusions

The high degree of synteny among felid genomes enabled the successful application of the domestic cat reference in high-quality SNV detection. The datasets presented here provide a useful resource for future studies into population dynamics, evolutionary history and genetic and disease management of big cats. This cross-species method of variant discovery provides genomic context for identifying annotated gene regions essential to understanding adaptive and deleterious variants that can improve conservation outcomes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07899-2.

Moving from trends to benchmarks by using regression tree analysis to find inbreeding thresholds in a critically endangered bird

Understanding how inbreeding affects endangered species in conservation breeding programs is essential for their recovery. The Hawaiian Crow ('Alalā) (Corvus hawaiiensis) is one of the world's most endangered birds. It went extinct in the wild in 2002, and, until recent release efforts starting in 2016, nearly all of the population remained under human care for conservation breeding. Using pedigree inbreeding coefficients (F), we evaluated the effects of inbreeding on Hawaiian Crow offspring survival and reproductive success. We used regression tree analysis to identify the level of inbreeding (i.e., inbreeding threshold) that explains a substantial decrease in 'Alalā offspring survival to recruitment. Similar to a previous study of inbreeding in 'Alalā, we found that inbreeding had a negative impact on offspring survival but that parental (vs. artificial) egg incubation improved offspring survival to recruitment. Furthermore, we found that inbreeding did not substantially affect offspring reproductive success, based on the assumption that offspring that survive to adulthood breed with distantly related mates. Our novel application of regression tree analysis showed that offspring with inbreeding levels exceeding F = 0.098 were 69% less likely to survive to recruitment than more outbred offspring, providing a specific threshold value for ongoing population management. Our results emphasize the importance of assessing inbreeding depression across all life history stages, confirm the importance of prioritizing parental over artificial egg incubation in avian conservation breeding programs, and demonstrate the utility of regression tree analysis as a tool for identifying inbreeding thresholds, if present, in any pedigree-managed population.

A novel framework for evaluating in situ breeding management strategies in endangered populations

Conservation breeding management aims to reduce inbreeding and maximize the retention of genetic diversity in endangered populations. However, breeding management of wild populations is still rare, and there is a need for approaches that provide data-driven evidence of the likelihood of success of alternative in situ strategies. Here, we provide an analytical framework that uses in silico simulations to evaluate, for real wild populations, (i) the degree of population-level inbreeding avoidance, (ii) the genetic quality of mating pairs, and (iii) the potential genetic benefits of implementing two breeding management strategies. The proposed strategies aim to improve the genetic quality of breeding pairs by splitting detrimental pairs and allowing the members to re-pair in different ways. We apply the framework to the wild population of the Critically Endangered helmeted honeyeater by combining genomic data and field observations to estimate the inbreeding (i.e., pair-kinship) and genetic quality (i.e., Mate Suitability Index) of all mating pairs for seven consecutive breeding seasons. We found no evidence of population-level inbreeding avoidance and that ~91.6% of breeding pairs were detrimental to the genetic health of the population. Furthermore, the framework revealed that neither proposed management strategy would significantly improve the genetic quality or reduce inbreeding of the mating pairs in this population. Our results demonstrate the usefulness of our analytical framework for testing the efficacy of different in situ breeding management strategies and for making evidence-based management decisions.

Meta-analytic evidence that animals rarely avoid inbreeding

Animals are usually expected to avoid mating with relatives (kin avoidance) as incestuous mating can lead to the expression of inbreeding depression. Yet, theoretical models predict that unbiased mating with regards to kinship should be common, and that under some conditions, the inclusive fitness benefits associated with inbreeding can even lead to a preference for mating with kin. This mismatch between empirical and theoretical expectations generates uncertainty as to the prevalence of inbreeding avoidance in animals. Here, we synthesized 677 effect sizes from 139 experimental studies of mate choice for kin versus non-kin in diploid animals, representing 40 years of research, using a meta-analytical approach. Our meta-analysis revealed little support for the widely held view that animals avoid mating with kin, despite clear evidence of publication bias. Instead, unbiased mating with regards to kinship appears widespread across animals and experimental conditions. The significance of a variety of moderators was explored using meta-regressions, revealing that the degree of relatedness and prior experience with kin explained some variation in the effect sizes. Yet, we found no difference in kin avoidance between males and females, choice and no-choice experiments, mated and virgin animals or between humans and animals. Our findings highlight the need to rethink the widely held view that inbreeding avoidance is a given in experimental studies.

Using historical genome-wide DNA to unravel the confused taxonomy in a songbird lineage that is extinct in the wild

Urgent conservation action for terminally endangered species is sometimes hampered by taxonomic uncertainty, especially in illegally traded animals that are often cross-bred in captivity. To overcome these problems, we used a genomic approach to analyze historical DNA from museum samples across the Asian Pied Starling (Gracupica contra) complex in tropical Asia, a popular victim of the ongoing songbird crisis whose distinct Javan population ("Javan Pied Starling") is extinct in the wild and subject to admixture in captivity. Comparing genomic profiles across the entire distribution, we detected three deeply diverged lineages at the species level characterized by a lack of genomic intermediacy near areas of contact. Our study demonstrates that the use of historical DNA can be instrumental in delimiting species in situations of taxonomic uncertainty, especially when modern admixture may obfuscate species boundaries. Results of our research will enable conservationists to commence a dedicated ex situ breeding program for the Javan Pied Starling, and serve as a blueprint for similar conservation problems involving terminally endangered species subject to allelic infiltration from close congeners.

Genetic Rescue and the Plight of Ponui Hybrids

Long-term sustainable and resilient populations is a key goal of conservation. How to best achieve this is controversial. There are, for instance, polarized views concerning the fitness and conservation value of hybrid populations founded through multi-origin translocations. A classic example concerns Apteryx (kiwi) in New Zealand. The A. mantelli of Ponui Island constitute a hybrid population where the birds are highly successful in their island habitat. A key dilemma for managers is understanding the reason for this success. Are the hybrid birds of Ponui Island of “no future conservation value” as recently asserted, or do they represent an outstanding example of genetic rescue and an important resource for future translocations? There has been a paradigm shift in scientific thinking concerning hybrids, but the ecological significance of admixed genomes remains difficult to assess. This limits what we can currently predict in conservation science. New understanding from genome science challenges the sufficiency of population genetic models to inform decision making and suggests instead that the contrasting outcomes of hybridization, “outbreeding depression” and “heterosis,” require understanding additional factors that modulate gene and protein expression and how these factors are influenced by the environment. We discuss these findings and the investigations that might help us to better understand the birds of Ponui, inform conservation management of kiwi and provide insight relevant for the future survival of Apteryx.

White-naped mangabeys' viable insurance population within European Zoo Network

The success and viability of an ex-situ conservation program lie in the establishment and potential maintenance of a demographically and genetically viable insurance population. Such population reserve may support reintroduction and reinforcement activities of wild populations. White-naped mangabeys are endangered restricted-range African primates which have experienced a dramatic population decrease in their natural habitats over the last few decades. Since 2001, some European zoos singularly monitor an ex-situ population aiming to seek the recovery of the current wild population. The aim of the present paper is to evaluate the genetic status and population demographics of European zoo-captive white-naped mangabeys based on pedigree data. The captive population is gradually growing and preserves specific reproductive and demographic parameters linked to the species. The intensive management program that is implemented has brought about the minimization of inbreeding and average relatedness levels, thus maintaining high levels of genetic diversity despite the existence of fragmented populations. This finding suggests white-naped mangabey ex-situ preservation actions may be a good example of multifaceted conservation throughout studbook management which could be used as a model for other ex-situ live-animal populations.

Polygenic basis for adaptive morphological variation in a threatened Aotearoa | New Zealand bird, the hihi (Notiomystis cincta)

To predict if a threatened species can adapt to changing selective pressures, it is crucial to understand the genetic basis of adaptive traits, especially in species historically affected by severe bottlenecks. We estimated the heritability of three hihi (Notiomystis cincta) morphological traits known to be under selection (nestling tarsus length, body mass and head-bill length) using 523 individuals and 39 699 single nucleotide polymorphisms (SNPs) from a 50 K Affymetrix SNP chip. We then examined the genetic architecture of the traits via chromosome partitioning analyses and genome-wide association scans (GWAS). Heritabilities estimated using pedigree relatedness or genomic relatedness were low. For tarsus length, the proportion of genetic variance explained by each chromosome was positively correlated with its size, and more than one chromosome explained significant variation for body mass and head-bill length. Finally, GWAS analyses suggested many loci of small effect contributing to trait variation for all three traits, although one locus (an SNP within an intron of the transcription factor HEY2) was tentatively associated with tarsus length. Our findings suggest a polygenic nature for the morphological traits, with many small effect size loci contributing to the majority of the variation, similar to results from many other wild populations. However, the small effective population size, polygenic architecture and already low heritabilities suggest that both the total response and rate of response to selection are likely to be limited in hihi.

An 85K SNP Array Uncovers Inbreeding and Cryptic Relatedness in an Antarctic Fur Seal Breeding Colony

High density single nucleotide polymorphism (SNP) arrays allow large numbers of individuals to be rapidly and cost-effectively genotyped at large numbers of genetic markers. However, despite being widely used in studies of humans and domesticated plants and animals, SNP arrays are lacking for most wild organisms. We developed a custom 85K Affymetrix Axiom array for an intensively studied pinniped, the Antarctic fur seal (Arctocephalus gazella). SNPs were discovered from a combination of genomic and transcriptomic resources and filtered according to strict criteria. Out of a total of 85,359 SNPs tiled on the array, 75,601 (88.6%) successfully converted and were polymorphic in 270 animals from a breeding colony at Bird Island in South Georgia. Evidence was found for inbreeding, with three genomic inbreeding coefficients being strongly intercorrelated and the proportion of the genome in runs of homozygosity being non-zero in all individuals. Furthermore, analysis of genomic relatedness coefficients identified previously unknown first-degree relatives and multiple second-degree relatives among a sample of ostensibly unrelated individuals. Such “cryptic relatedness” within fur seal breeding colonies may increase the likelihood of consanguineous matings and could therefore have implications for understanding fitness variation and mate choice. Finally, we demonstrate the cross-amplification potential of the array in three related pinniped species. Overall, our SNP array will facilitate future studies of Antarctic fur seals and has the potential to serve as a more general resource for the wider pinniped research community.