The cryptic genetic structure of the North American captive gorilla population

Western lowland gorillas (Gorilla gorilla gorilla) were imported from across their geographical range to North American zoos from the late 1800s through 1974. The majority of these gorillas were imported with little or no information regarding their original provenance and no information on their genetic relatedness. Here, we analyze 32 microsatellite loci in 144 individuals using a Bayesian clustering method to delineate clusters of individuals among a sample of founders of the captive North American zoo gorilla collection. We infer that the majority of North American zoo founders sampled are distributed into two distinct clusters, and that some individuals are of admixed ancestry. This new information regarding the existence of ancestral genetic population structure in the North American zoo population lays the groundwork for enhanced efforts to conserve the evolutionary units of the western lowland gorilla gene pool. Our data also show that the genetic diversity estimates in the founder population were comparable to those in wild gorilla populations (Mondika and Cross River), and that pairwise relatedness among the founders is no different from that expected for a random mating population. However, the relatively high level of relatedness (R = 0.54) we discovered in a pair of known breeding pairs reveals the need for incorporating genetic relatedness estimates in the captive management of western lowland gorillas.     

Performance of marker-based relatedness estimators in natural populations of outbred vertebrates

Knowledge of relatedness between pairs of individuals plays an important role in many research areas including evolutionary biology, quantitative genetics, and conservation. Pairwise relatedness estimation methods based on genetic data from highly variable molecular markers are now used extensively as a substitute for pedigrees. Although the sampling variance of the estimators has been intensively studied for the most common simple genetic relationships, such as unrelated, half- and full-sib, or parent-offspring, little attention has been paid to the average performance of the estimators, by which we mean the performance across all pairs of individuals in a sample. Here we apply two measures to quantify the average performance: first, misclassification rates between pairs of genetic relationships and, second, the proportion of variance explained in the pairwise relatedness estimates by the true population relatedness composition (i.e., the frequencies of different relationships in the population). Using simulated data derived from exceptionally good quality marker and pedigree data from five long-term projects of natural populations, we demonstrate that the average performance depends mainly on the population relatedness composition and may be improved by the marker data quality only within the limits of the population relatedness composition. Our five examples of vertebrate breeding systems suggest that due to the remarkably low variance in relatedness across the population, marker-based estimates may often have low power to address research questions of interest.     

Application of a noninvasive, PCR-based test for sex identification in an endangered parrot, Amazona guildingii

The St. Vincent Amazon Parrot Consortium was established during the 1980s in an effort to cooperatively manage the captive population of the St. Vincent parrot (Amazona guildingii) and support conservation of the wild population. Ex situ management of A. guildingii has been hindered by the sexual monomorphism of the group, mediated in the past through the use of time-consuming, expert-driven, and sometimes hazardous sexing procedures. In this article, we apply a noninvasive, molecular sexing technique, using the polymerase chain reaction (PCR) and DNA from a single feather tip to the captive populations of A. guildingii residing on St. Vincent (66 individuals) and Barbados (13 individuals). This study allowed for the rapid assessment of gender, while posing no threat to individual health, and will facilitate the efforts of the consortium breeding programs in the United States, Europe and on the islands of St. Vincent and Barbados. Zoo Biol 20:41-45, 2001. Copyright 2001 Wiley-Liss, Inc.     

A genetic evaluation of relatedness for broodstock management of captive,endangered Snake River sockeye salmon, <Emphasis Type="Italic">Oncorhynchus nerka</Emphasis>

The use of captive broodstocks is becoming more frequently employed as the number of species facing endangerment or extinction throughout the world increases. Efforts to rebuild the endangered Snake River sockeye salmon, Oncorhynchus nerka, population have been ongoing for over a decade, but the use of microsatellite data to develop inbreeding avoidance matrices is a more recent component to the program. This study used known genealogical relationships among sockeye salmon offspring to test four different pairwise relatedness estimators and a maximum-likelihood (M-L) relatedness estimator. The goal of this study was to develop a breeding strategy with these estimators that would minimize the loss of genetic diversity, minimize inbreeding, and determine how returning anadromous adults are incorporated into the broodstock along with full-term hatchery adults. Results of this study indicated that both the M _xy and R _QG estimators had the lowest Type II error rates and the M-L and R _R estimators had the lowest Type I error rates. An approach that utilizes a combination of estimators may provide the most valuable information for managers. We recommend that the M-L and R _R methods be used to rank the genetic importance of returning adults and the M _xy or R _QG estimators be used to determine which fish to pair for spawning. This approach provides for the best genetic management of this captive, endangered population and should be generally applicable to the genetic management of other endangered stocks with no pedigree.     

Maximum-likelihood estimation of relatedness

Relatedness between individuals is central to many studies in genetics and population biology. A variety of estimators have been developed to enable molecular marker data to quantify relatedness. Despite this, no effort has been given to characterize the traditional maximum-likelihood estimator in relation to the remainder. This article quantifies its statistical performance under a range of biologically relevant sampling conditions. Under the same range of conditions, the statistical performance of five other commonly used estimators of relatedness is quantified. Comparison among these estimators indicates that the traditional maximum-likelihood estimator exhibits a lower standard error under essentially all conditions. Only for very large amounts of genetic information do most of the other estimators approach the likelihood estimator. However, the likelihood estimator is more biased than any of the others, especially when the amount of genetic information is low or the actual relationship being estimated is near the boundary of the parameter space. Even under these conditions, the amount of bias can be greatly reduced, potentially to biologically irrelevant levels, with suitable genetic sampling. Additionally, the likelihood estimator generally exhibits the lowest root mean-square error, an indication that the bias in fact is quite small. Alternative estimators restricted to yield only biologically interpretable estimates exhibit lower standard errors and greater bias than do unrestricted ones, but generally do not improve over the maximum-likelihood estimator and in some cases exhibit even greater bias. Although some nonlikelihood estimators exhibit better performance with respect to specific metrics under some conditions, none approach the high level of performance exhibited by the likelihood estimator across all conditions and all metrics of performance.     

swinger: a user‐friendly computer program to establish captive breeding groups that minimize relatedness without pedigree information

Captive breeding programmes are often a necessity for the continued persistence of a population or species. They typically have the goal of maintaining genetic diversity and minimizing inbreeding. However, most captive breeding programmes have been based on the assumption that the founding breeders are unrelated and outbred, even though in situ anthropogenic impacts often mean these founders may have high relatedness and substantial inbreeding. In addition, polygamous group‐breeding species in captivity often have uncertain pedigrees, making it difficult to select the group composition for subsequent breeding. Molecular‐based estimates of relatedness and inbreeding may instead be used to select breeding groups (≥two individuals) that minimize relatedness and filter out inbred individuals. swinger constructs breeding groups based on molecular estimates of relatedness and inbreeding. The number of possible combinations of breeding groups quickly becomes intractable by hand. swinger was designed to overcome this major issue in ex situ conservation biology. The user can specify parameters within swinger to reach breeding solutions that suit the mating system of the target species and available resources. We provide evidence of the efficiency of the software with an empirical example and using simulations. The only data required are a typical molecular marker data set, such as a microsatellite or SNP data set, from which estimates of inbreeding and pairwise relatedness may be obtained. Such molecular data sets are becoming easier to gather from non‐model organisms with next‐generation sequencing technology. swinger is an open‐source software with a user‐friendly interface and is available at http://www.molecularecology.flinders.edu.au/molecular-ecology-lab/software/swinger/swinger/ and https://github.com/Yuma248/Swinger .     

The Genetic Integrity of the Ex Situ Population of the European Wildcat (Felis silvestris silvestris) Is Seriously Threatened by Introgression from Domestic Cats (Felis silvestris catus)

Studies on the genetic diversity and relatedness of zoo populations are crucial for implementing successful breeding programmes. The European wildcat, Felis s. silvestris, is subject to intensive conservation measures, including captive breeding and reintroduction. We here present the first systematic genetic analysis of the captive population of Felis s. silvestris in comparison with a natural wild population. We used microsatellites and mtDNA sequencing to assess genetic diversity, structure and integrity of the ex situ population. Our results show that the ex situ population of the European wildcat is highly structured and that it has a higher genetic diversity than the studied wild population. Some genetic clusters matched the breeding lines of certain zoos or groups of zoos that often exchanged individuals. Two mitochondrial haplotype groups were detected in the in situ populations, one of which was closely related to the most common haplotype found in domestic cats, suggesting past introgression in the wild. Although native haplotypes were also found in the captive population, the majority (68%) of captive individuals shared a common mtDNA haplotype with the domestic cat (Felis s. catus). Only six captive individuals (7.7%) were assigned as wildcats in the STRUCTURE analysis (at K = 2), two of which had domestic cat mtDNA haplotypes and only two captive individuals were assigned as purebred wildcats by NewHybrids. These results suggest that the high genetic diversity of the captive population has been caused by admixture with domestic cats. Therefore, the captive population cannot be recommended for further breeding and reintroduction.     

Monitoring and managing genetic variation in group breeding populations without individual pedigrees

The genetic management of captive populations to conserve genetic variation is currently based on analyses of individual pedigrees to infer inbreeding and kinship coefficients and values of individuals as breeders. Such analyses require that individual pedigrees are known and individual pairing (mating) can be controlled. Many species in captivity, however, breed in groups due to various reasons, such as space constraints and fertility considerations for species living naturally in social groups, and thus have no pedigrees available for the traditional genetic analyses and management. In the absence of individual pedigree, such group breeding populations can still be genetically monitored, evaluated and managed by suitable population genetics models using population level information (such as census data). This article presents a simple genetic model of group breeding populations to demonstrate how to estimate the genetic variation maintained within and among populations and to optimise management based on these estimates. A numerical example is provided to illustrate the use of the proposed model. Some issues relevant to group breeding, such as the development and robustness evaluation of the population genetics model appropriate for a particular species under specific management and recording systems and the genetic monitoring with markers, are also briefly discussed.     

Population genetics and relatedness in a critically endangered island raptor, Ridgway’s Hawk Buteo ridgwayi

Many island avian populations are of conservation interest because they have a higher risk of extinction than mainland populations. Susceptibility of island birds to extinction is primarily related to human induced change through habitat loss, persecution, and introduction of exotic species, in combination with genetic factors. We used microsatellite profiles from 11 loci to assess genetic diversity and relatedness in the critically endangered hawk Buteo ridgwayi endemic to the island of Hispaniola in the Caribbean. Using samples collected between 2005 and 2009, our results revealed a relatively high level of heterozygosity, evidence of a recent genetic bottleneck, and the occurrence of inbreeding within the population. Pair relatedness analysis found 4 of 7 sampled breeding pairs to be related similar to that of first cousin or greater. Pedigree estimates indicated that up to 18 % of potential pairings would be between individuals with relatedness values similar to that of half-sibling. We discuss our findings in the context of conservation genetic management suggesting both carefully managed translocations and the initiation of a captive population as a safeguard of the remaining genetic diversity.     

Isolation and characterization of 23 polymorphic microsatellite loci for a West Indian iguana (Cyclura pinguis) from the British Virgin Islands

Twenty-three polymorphic microsatellite markers were identified and characterized for Cyclura pinguis, a critically endangered species of lizard (Sauria: Iguanidae) native to Anegada Island in the British Virgin Islands. We examined variation at these loci for 39 C. pinguis, finding up to five alleles per locus and an average expected heterozygosity of 0.55. Allele frequency estimates for these microsatellite loci will be used to characterize genetic diversity of captive and wild C. pinguis populations and to estimate relatedness among adult iguanas at the San Diego Zoo that form the nucleus of a captive breeding programme for this critically endangered species.     

Conservation Genetics of Threatened Hippocampus guttulatus in Vulnerable Habitats in NW Spain: Temporal and Spatial Stability of Wild Populations with Flexible Polygamous Mating System in Captivity

This study was focused on conservation genetics of threatened Hippocampus guttulatus on the Atlantic coast of NW Iberian Peninsula. Information about spatial structure and temporal stability of wild populations was obtained based on microsatellite markers, and used for monitoring a captive breeding program firstly initiated in this zone at the facilities of the Institute of Marine Research (Vigo, Spain). No significant major genetic structure was observed regarding the biogeographical barrier of Cape Finisterre. However, two management units under continuous gene flow are proposed based on the allelic differentiation between South-Atlantic and Cantabrian subpopulations, with small to moderate contemporary effective size based on single-sample methods. Temporal stability was observed in South-Atlantic population samples of H. guttulatus for the six-year period studied, suggesting large enough effective population size to buffer the effects of genetic drift within the time frame of three generations. Genetic analysis of wild breeders and offspring in captivity since 2009 allowed us to monitor the breeding program founded in 2006 in NW Spain for this species. Similar genetic diversity in the renewed and founder broodstock, regarding the wild population of origin, supports suitable renewal and rearing processes to maintain genetic variation in captivity. Genetic parentage proved single-brood monogamy in the wild and in captivity, but flexible short- and long-term mating system under captive conditions, from strict monogamy to polygamy within and/or among breeding seasons. Family analysis showed high reproductive success in captivity under genetic management assisted by molecular relatedness estimates to avoid inbreeding. This study provides genetic information about H. guttulatus in the wild and captivity within an uncovered geographical range for this data deficient species, to be taken into account for management and conservation purposes.     

Genetic evidence of hybridization between the critically endangered Cuban crocodile and the American crocodile: implications for population history and in situ/ex situ conservation

Inter-specific hybridization may be especially detrimental when one species is extremely rare and the other is abundant owing to the potential for genetic swamping. The Cuban crocodile (Crocodylus rhombifer) is a critically endangered island endemic largely restricted to Zapata Swamp, where it is sympatric with the widespread American crocodile (C. acutus). An on-island, C. rhombifer captive breeding program is underway with the goals of maintaining taxonomic integrity and providing a source of individuals for reintroduction, but its conservation value is limited by lack of genetic information. Here we collected mtDNA haplotypic and nuclear genotypic data from wild and captive C. rhombifer and C. acutus in Cuba to: (1) investigate the degree of inter-specific hybridization in natural (in situ) and captive (ex situ) populations; (2) quantify the extent, distribution and in situ representation of genetic variation ex situ; and (3) reconstruct founder relatedness to inform management. We found high levels of hybridization in the wild (49.1%) and captivity (16.1%), and additional evidence for a cryptic lineage of C. acutus in the Antilles. We detected marginally higher observed heterozygosity and allelic diversity ex situ relative to the wild population, with captive C. rhombifer exhibiting over twice the frequency of private alleles. Although mean relatedness was high in captivity, we identified 37 genetically important individuals that possessed individual mean kinship (MK) values lower than the population MK. Overall, these results will guide long-term conservation management of Cuban crocodiles for maintaining the genetic integrity and viability of this species of high global conservation value.     

Potential genetic consequences of a recent bottleneck in the Amur tiger of

The Amur tiger, Panthera tigris altaica, is a highly endangered felid whose range and population size has been severely reduced in recent times. At present, the wild population is estimated at 490 individuals, having rebounded from the 20–30 tigers remaining following a severe bottleneck in the 1940's. The current study presents preliminary data on the patterns and levels of genetic variation in the mitochondrial DNA control region using DNA extracted from non-invasively sampled faecal material, collected throughout the entire range of P. t. altaica in the Russian Far East. Analysis of 82 scat samples representing at least 27 individuals revealed extremely low levels of CR haplotype diversity, characterized by a single widespread haplotype (96.4%) and two rare variants, each differing by a single step within the hypervariable I (2.4%) and central conserved regions (1.2%), respectively. A comparison with previous data on cytochrome bvariation in 14 captive individuals revealed a potentially greater amount of genetic variation represented in captivity relative to that found in the wild population. The extremely low levels of mitochondrial DNA variation in the wild population is discussed in light of the demographic processes that might have shaped these patterns as well as the potential bias introduced through analysis of fecal samples. These results highlight the continuing need to assess levels of genetic variation even in recovering populations that are increasing in number and underscore the important role that captive breeding programs may play in preserving remnant genetic diversity of endangered species.     

Analysis of population genetic structure with RAPD markers

Recent advances in the application of the polymerase chain reaction make it possible to score individuals at a large number of loci. The RAPD (random amplified polymorphic DNA) method is one such technique that has attracted widespread interest. The analysis of population structure with RAPD data is hampered by the lack of complete genotypic information resulting from dominance, since this enhances the sampling variance associated with single loci as well as induces bias in parameter estimation. We present estimators for several population-genetic parameters (gene and genotype frequencies, within- and between-population heterozygosities, degree of inbreeding and population subdivision, and degree of individual relatedness) along with expressions for their sampling variances. Although completely unbiased estimators do not appear to be possible with RAPDs, several steps are suggested that will insure that the bias in parameter estimates is negligible. To achieve the same degree of statistical power, on the order of 2 to 10 times more individuals need to be sampled per locus when dominant markers are relied upon, as compared to codominant (RFLP, isozyme) markers. Moreover, to avoid bias in parameter estimation, the marker alleles for most of these loci should be in relatively low frequency. Due to the need for pruning loci with low-frequency null alleles, more loci also need to be sampled with RAPDs than with more conventional markers, and some problems of bias cannot be completely eliminated.     

Estimating genetic correlations in natural populations

Information on the genetic correlation between traits provides fundamental insight into the constraints on the evolutionary process. Estimates of such correlations are conventionally obtained by raising individuals of known relatedness in artificial environments. However, many species are not readily amenable to controlled breeding programmes, and considerable uncertainty exists over the extent to which estimates derived under benign laboratory conditions reflect the properties of populations in natural settings. Here, non-invasive methods that allow the estimation of genetic correlations from phenotypic measurements derived from individuals of unknown relatedness are introduced. Like the conventional approach, these methods demand large sample sizes in order to yield reasonably precise estimates, and special precautions need to be taken to eliminate bias from shared environmental effects. Provided the sample consists of at least 20% or so relatives, informative estimates of the genetic correlation are obtainable with sample sizes of several hundred individuals, particularly if supplemental information on relatedness is available from polymorphic molecular markers.     

A comparison of microsatellite-based pairwise relatedness estimators

Studies of inbreeding depression or kin selection require knowledge of relatedness between individuals. If pedigree information is lacking, one has to rely on genotypic information to infer relatedness. In this study we investigated the performance (absolute and relative) of 10 marker-based relatedness estimators using allele frequencies at microsatellite loci obtained from natural populations of two bird species and one mammal species. Using Monte Carlo simulations we show that many factors affect the performance of estimators and that different sets of loci promote the use of different estimators: in general, there is no single best-performing estimator. The use of locus-specific weights turns out to greatly improve the performance of estimators when marker loci are used that differ strongly in allele frequency distribution. Microsatellite-based estimates are expected to explain between 25 and 79% of variation in true relatedness depending on the microsatellite dataset and on the population composition (i.e. the frequency distribution of relationship in the population). We recommend performing Monte Carlo simulations to decide which estimator to use in studies of pairwise relatedness.     

Integrating microsatellite and pedigree analyses to facilitate the captive management of the endangered Mississippi sandhill crane (Grus canadensis pulla)

The minimization of kinship in captive populations is usually achieved through the use of pedigree information. However, pedigree knowledge alone is not sufficient if pedigree information is missing, questionable, or when the founders of the captive population are related to one another. If this is the case, higher levels of inbreeding and lower levels of genetic diversity may be present in a captive population than those calculated by pedigree analyses alone. In this study, the genetic status of the critically endangered Mississippi sandhill crane (MSC) (Grus canadensis pulla) was analyzed using studbook data from the U.S. Fish and Wildlife Service managed captive breeding program as well as microsatellite DNA data. These analyses provided information on shared founder genotypes, allowing for refined analysis of genetic variation in the population, and the development of a new DNA-based studbook pedigree that will assist in the genetic management of the MSC population.     

Genetic quality of the Miyaluo captive forest musk deer (Moschus berezovskii) population as assessed by microsatellite loci

The Miyaluo captive forest musk deer population (Sichuan Province, China) is one of the largest captive breeding populations in the world. In order to evaluate the genetic quality and provide available genetic management strategy, seven polymorphism microsatellite loci were applied to assess the genetic variation of the Miyaluo forest musk deer. The results indicated that a total of 168 alleles were detected from these seven microsatellite loci in 361 individuals, and the number of the alleles per locus ranged from 12 to 41 with a mean of 24. The average observed heterozygosity, expected heterozygosity, and PIC were 0.782, 0.854, and 0.837, respectively. Considering the results of the loci Hardy–Weinberg equilibrium test, the comparison of the common allele frequency as well as the private allele between the adults and juveniles, we concluded that the heterozygosity and the genetic diversity of the Miyaluo captive breeding population are increasing due to the input of new individuals from other populations. However, the frequency of some alleles declined sharply, and some were even lost indicating that there is a risk for diversity loss. Thus, we proposed an improved management and breeding strategy for the captive breeding population of the forest musk deer.     

Estimation of relatedness among non-pedigreed Yakutian cryo-bank bulls using molecular data: implications for conservation and breed management

Background

Yakutian cattle, the last remaining native cattle breed in Siberia, are well adapted to the extreme sub-arctic conditions. Nowadays only ca. 1200 purebred animals are left in Yakutia. The semen of six Yakutian bulls was stored in a cryo-bank without any pedigree documentation because of the traditional free herding style of the population.

Methods

To clarify the genetic relatedness between these bulls and to provide recommendations to use their semen in future conservation and breed management programs, we have analysed 30 autosomal microsatellites and mitochondrial DNA sequences in 60 individuals including the six for which semen has been stored. Four relatedness estimators were calculated. In addition, we assessed the value of the cryo-bank bulls for the preservation of genetic variation of the contemporary Yakutian cattle by calculating allelic and gene diversity estimates and mean molecular coancestries.

Results

On the basis of microsatellite variability, including the Yakutian cryo-bank bulls increases the allelic variation in the contemporary population by 3% and in the male subpopulation by 13%. In terms of the mean molecular coancestries, they are less related to the contemporary cow population than the breeding bulls and therefore could be used to reduce inbreeding in the living population. Although 30 loci are insufficient to resolve definitely their relatedness categories, the data suggest four pairs of cryo-bank bulls as possible half-sibs.

Conclusions

Our results show that even relatively limited cryo-bank storage of semen can carry allelic variation through a bottleneck. We propose a breeding scheme based on the rotation of breeding females and the division of cryo-bank bulls into three groups. Thus, if molecular data (e.g. autosomal microsatellite genotypes) for the contemporary population are available and based on relatively small-scale laboratory analyses, it is possible to avoid serious mistakes in their use for breeding applications. The approach suggested here based on the use of Yakutian cryo-bank semen can be easily extended to cryo-bank materials of other animals in future breeding programs.     

Marker-based estimation of the coefficient of coancestry in hybrid breeding programmes

Molecular markers allow to estimate the pairwise relatedness between the members of a breeding pool when their selection history is no longer available or has become too complex for a classical pedigree analysis. The field of population genetics has several estimation procedures at its disposal, but when the genotyped individuals are highly selected inbred lines, their application is not warranted as the theoretical assumptions on which these estimators were built, usually linkage equilibrium between marker loci or even Hardy–Weinberg equilibrium, are not met. An alternative approach requires the availability of a genotyped reference set of inbred lines, which allows to correct the observed marker similarities for their inherent upward bias when used as a coancestry measure. However, this approach does not guarantee that the resulting coancestry matrix is at least positive semi-definite (psd), a necessary condition for its use as a covariance matrix. In this paper we present the weighted alikeness in state (WAIS) estimator. This marker-based coancestry estimator is compared to several other commonly applied relatedness estimators under realistic hybrid breeding conditions in a number of simulations. We also fit a linear mixed model to phenotypical data from a commercial maize breeding programme and compare the likelihood of the different variance structures. WAIS is shown to be psd which makes it suitable for modelling the covariance between genetic components in linear mixed models involved in breeding value estimation or association studies. Results indicate that it generally produces a low root mean squared error under different breeding circumstances and provides a fit to the data that is comparable to that of several other marker-based alternatives. Recommendations for each of the examined coancestry measures are provided.