Using microsatellite diversity in wild Anegada iguanas (Cyclura pinguis) to establish relatedness in a captive breeding group of this critically endangered species

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

Conservation management strategies for many highly threatened species include conservation breeding to prevent extinction and enhance recovery. Pairing decisions for these conservation breeding programmes can be informed by pedigree data to minimize relatedness between individuals in an effort to avoid inbreeding, maximize diversity and maintain evolutionary potential. However, conservation breeding programmes struggle to use this approach when pedigrees are shallow or incomplete. While genetic data (i.e., microsatellites) can be used to estimate relatedness to inform pairing decisions, emerging evidence indicates this approach may lack precision in genetically depauperate species, and more effective estimates will likely be obtained from genomic data (i.e., thousands of genome-wide single nucleotide polymorphisms, or SNPs). Here, we compare relatedness estimates and subsequent pairing decisions using pedigrees, microsatellites and SNPs from whole-genome resequencing approaches in two critically endangered birds endemic to New Zealand: kakī/black stilt (Himantopus novaezelandiae) and kākāriki karaka/orange-fronted parakeet (Cyanoramphus malherbi). Our findings indicate that SNPs provide more precise estimates of relatedness than microsatellites when assessing empirical parent-offspring and full sibling relationships. Further, our results show that relatedness estimates and subsequent pairing recommendations using PMx are most similar between pedigree- and SNP-based approaches. These combined results indicate that in lieu of robust pedigrees, SNPs are an effective tool for informing pairing decisions, which has important implications for many poorly pedigreed conservation breeding programmes worldwide.

Heterozygosity-Fitness Correlations Reveal Inbreeding Depression in Neonatal Body Size in a Critically Endangered Rock Iguana

Inbreeding depression, though challenging to identify in nature, may play an important role in regulating the dynamics of small and isolated populations. Conversely, greater expression of genetic load can enhance opportunities for natural selection. Conditional expression concentrates these opportunities for selection and may lead to failure of detection. This study investigates the possibility for age-dependent expression of inbreeding depression in a critically endangered population of rock iguanas, Cyclura nubila caymanensis. We employ heterozygote-fitness correlations to examine the contributions of individual genetic factors to body size, a fitness-related trait. Nonsignificant reductions in homozygosity (up to 7%) were detected between neonates and individuals surviving past their first year, which may reflect natural absorption of inbreeding effects by this small, fecund population. The majority of variation in neonate body size was attributed to maternal or environmental effects (i.e., clutch identity and incubation length); however, heterozygosity across 22 microsatellite loci also contributed significantly and positively to model predictions. Conversely, effects of heterozygosity on fitness were not detectable when adults were examined, suggesting that inbreeding depression in body size may be age dependent in this taxon. Overall, these findings emphasize the importance of taking holistic, cross-generational approaches to genetic monitoring of endangered populations.

Initial founders of captive populations are genetically representative of natural populations in critically endangered dusky gopher frogs, Lithobates sevosus

The rapid rate of decline in amphibian populations has urged many researchers and conservationists to establish captive, or ex situ, populations. Such populations are guarded against effects of habitat loss and degradation, and if actively managed, can serve as a reservoir for rare alleles that might be lost in the wild. Without proper management, ex situ population sizes can dwindle and will no longer perform this function. The dusky gopher frog, Lithobates sevosus, is a critically endangered species, imperiled by habitat loss and population isolation. To assist in recovery of the species and prevent further genetic erosion, a captive breeding program was initiated. We investigated how well natural genetic variation was captured within the ex situ population and determined relatedness within each ex situ population. We genotyped individuals from two natural populations and two founding, captive populations to compare metrics of genetic variation and relatedness. The data show the initial founder populations are genetically representative of the natural populations, although variation is low in each, and that relatedness values are similar. Therefore, founding captive populations were successful at capturing genetic variation in the wild. Future research should continue to compare genetic variation of captive and natural populations to monitor efficacy of their management programs. Zoo Biol. 35:378-384, 2016. © 2016 Wiley Periodicals, Inc.

Applying SNP-Derived Molecular Coancestry Estimates to Captive Breeding Programs

Captive breeding programs for wildlife species typically rely on pedigrees to inform genetic management. Although pedigree-based breeding strategies are quite effective at retaining long-term genetic variation, management of zoo-based breeding programs continues to be hampered when pedigrees are poorly known. The objective of this study was to evaluate 2 options for generating single nucleotide polymorphism (SNP) data to resolve unknown relationships within captive breeding programs. We generated SNP data for a zoo-based population of addax (Addax nasomasculatus) using both the Illumina BovineHD BeadChip and double digest restriction site-associated DNA (ddRAD) sequencing. Our results demonstrated that estimates of allele sharing (AS) between pairs of individuals exhibited low variances. Average AS variances were highest when using 50 loci (SNPchipall = 0.00159; ddRADall = 0.0249), but fell below 0.0003 for the SNP chip dataset when sampling ≥250 loci and below 0.0025 for the ddRAD dataset when sampling ≥500 loci. Furthermore, the correlation between the SNPchipall and ddRADall AS datasets was 0.88 (95%CI = 0.84-0.91) when subsampling 500 loci. Collectively, our results indicated that both SNP genotyping methods produced sufficient data for accurately estimating relationships, even within an extremely bottlenecked population. Our results also suggested that analytic assumptions historically integrated into the addax pedigree are not adversely impacting long-term pedigree-based management; kinships calculated from the analytic pedigree were significantly correlated (P << 0.001) with AS estimates. Overall, our conclusions are intended to serve as both a proof of concept and a model for applying molecular data to the genetic management of captive breeding programs.

Genetic diversity and structure in the Endangered Allen Cays Rock Iguana, Cyclura cychlura inornata

The Endangered Allen Cays Rock Iguana (Cyclura cychlura inornata) is endemic to the Allen Cays, a tiny cluster of islands in the Bahamas. Naturally occurring populations exist on only two cays (<4 ha each). However, populations of unknown origin were recently discovered on four additional cays. To investigate patterns of genetic variation among these populations, we analyzed nuclear and mitochondrial markers for 268 individuals. Analysis of three mitochondrial gene regions (2,328 bp) and data for eight nuclear microsatellite loci indicated low genetic diversity overall. Estimates of effective population sizes based on multilocus genotypes were also extremely low. Despite low diversity, significant population structuring and variation in genetic diversity measures were detected among cays. Genetic data confirm the source population for an experimentally translocated population while raising concerns regarding other, unauthorized, translocations. Reduced heterozygosity is consistent with a documented historical population decline due to overharvest. This study provides the first range-wide genetic analysis of this subspecies. We suggest strategies to maximize genetic diversity during ongoing recovery including additional translocations to establish assurance populations and additional protective measures for the two remaining natural populations.

Testing the regional genetic representativeness of captive koala populations in South-East Queensland

Context Captive breeding for release back to the wild is an important component of ex situ conservation but requires genetic diversity that is representative of the wild population and has the ultimate goal of producing ecologically sustainable and resilient populations. However, defining and testing for representativeness of captive populations is difficult. Koalas (Phascolarctos cinereus) are bred for educational and tourism purposes in zoos and wildlife parks in South-East Queensland, but there are drastic declines evident in some wild koala populations in this region. Aim We compared genetic diversity at microsatellite loci and mitochondrial DNA in two captive koala populations with that of the local, wild koalas of South-East Queensland, determining the degree to which genetic diversity of neutral loci had been preserved and was represented in the captive populations. Key results The expected heterozygosity and the allelic richness was significantly greater in one captive colony than one wild South-East Queensland population. There was low but significant differentiation of the captive from wild populations using FST, with greater differentiation described by Jost’s Dest. In contrast, a newly introduced Kullback–Leibler divergence measure, which assesses similarity of allele frequencies, showed no significant divergence of colony and wild populations. The captive koalas lacked many of the mitochondrial haplotypes identified from South-East Queensland koalas and possessed seven other haplotypes. Conclusions Captive colonies of koalas have maintained levels of overall neutral genetic diversity similar to wild populations at microsatellite loci and low but significant differentiation likely resulted from drift and founder effects in small captive colonies or declining wild populations. Mitochondrial DNA suggests that captive founders were from a wider geographic source or that haplotypes have been lost locally. Implications Overall, tested captive koalas maintain sufficient microsatellite diversity to act as an in situ reservoir for neutral genetic diversity of regional populations.

Estimating relatedness and inbreeding using molecular markers and pedigrees: the effect of demographic history

Estimates of inbreeding and relatedness are commonly calculated using molecular markers, although the accuracy of such estimates has been questioned. As a further complication, in many situations, such estimates are required in populations with reduced genetic diversity, which is likely to affect their accuracy. We investigated the correlation between microsatellite- and pedigree-based coefficients of inbreeding and relatedness in laboratory populations of Drosophila melanogaster that had passed through bottlenecks to manipulate their genetic diversity. We also used simulations to predict expected correlations between marker- and pedigree-based estimates and to investigate the influence of linkage between loci and null alleles. Our empirical data showed lower correlations between marker- and pedigree-based estimates in our control (nonbottleneck) population than were predicted by our simulations or those found in similar studies. Correlations were weaker in bottleneck populations, confirming that extreme reductions in diversity can compromise the ability of molecular estimates to detect recent inbreeding events. However, this result was highly dependent on the strength of the bottleneck and we did not observe or predict any reduction in correlations in our population that went through a relatively severe bottleneck of N = 10 for one generation. Our results are therefore encouraging, as molecular estimates appeared robust to quite severe reductions in genetic diversity. It should also be remembered that pedigree-based estimates may not capture realized identity-by-decent and that marker-based estimates may actually be more useful in certain situations.

Kinship-based management strategies for captive breeding programs when pedigrees are unknown or uncertain

Zoo-based captive breeding programs typically rely on accurate pedigrees to maintain long-term population genetic diversity and prevent close inbreeding. For many mixed-sex captive populations, it is difficult to assign parentage of offspring with certainty without conducting DNA-based parentage analyses. Using the demographic parameters of a North American captive population of Arabian oryx (Oryx leucoryx), 2 kinship-based breeding-pair selection strategies were modeled for their performance in handling pedigrees with varying degrees of parentage uncertainty. We also compared these strategies with 2 nonkinship-based methods. Pedigrees simulated under different management strategies were compared for their long-term ability to maintain gene diversity (GD) and avoid inbreeding. For the Arabian oryx, results indicate that recording multiple possible parents instead of removing the unknown genomic portion of the pedigree can more efficiently utilize all animals available for breeding without compromising GD and inbreeding avoidance. Both kinship-based breeding-pair selection strategies significantly outperformed the nonkinship-based strategies.