Identification of kin structure among Guam rail founders: a comparison of pedigrees and DNA profiles

Kin structure among founders can have a significant effect on subsequent population structure. Here we use the correlation between DNA profile similarity and relatedness calculated from pedigrees to test hypotheses regarding kin structure among founders to the captive Guam rail (Rallus owstoni) population. Five different pedigrees were generated under the following hypotheses: (i) founders are unrelated; (ii) founders are unrelated except for same-nest chicks; (iii) founders from the same major site are siblings; (iv) founders from the same local site are siblings; and (v) founders are related as defined by a UPGMA cluster analysis of DNA similarity data. Relatedness values from pedigrees 1, 2 and 5 had the highest correlation with DNA similarity but the correlation between relatedness and similarity were not significantly different among pedigrees. Pedigree 5 resulted in the highest correlation overall when using only relatedness values that changed as a result of different founder hypotheses. Thus, founders were assigned relatedness based on pedigree 5 because it had the highest correlations with DNA similarity, was the most conservative approach, and incorporated all field data. The analyses indicated that estimating relatedness using DNA profiles remains problematic, therefore we compared mean kinship, a measure of genetic importance, with mean DNA profile similarity to determine if genetic importance among individuals could be determined via use of DNA profiles alone. The significant correlation suggests this method may provide more information about population structure than was previously thought. Thus, DNA profiles can provide a reasonable explanation for founder relatedness and mean DNA profile similarity may be helpful in determining relative genetic importance of individuals when detailed pedigrees are absent.     

Inbreeding and inbreeding depression in endangered red wolves (Canis rufus)

In natural populations, the expression and severity of inbreeding depression can vary widely across taxa. Describing processes that influence the extent of inbreeding and inbreeding depression aid in our understanding of the evolutionary history of mating systems such as cooperative breeding and nonrandom mate selection. Such findings also help shape wildlife conservation theory because inbreeding depression reduces the viability of small populations. We evaluated the extent of inbreeding and inbreeding depression in a small, re‐introduced population of red wolves (Canis rufus) in North Carolina. Since red wolves were first re‐introduced in 1987, pedigree inbreeding coefficients (f) increased considerably and almost every wild born wolf was inbred (average f = 0.154 and max f = 0.383). The large inbreeding coefficients were due to both background relatedness associated with few founders and numerous close relative matings. Inbreeding depression was most evident for adult body size and generally absent for direct fitness measures such as reproductive success and survival; no lethal equivalents (LE = 0.00) were detected in juvenile survival. The lack of strong inbreeding depression in direct measures of fitness could be due to a founder effect or because there were no outbred individuals for comparison. Our results highlight the variable expression of inbreeding depression across traits and the need to measure a number of different traits when evaluating inbreeding depression in a wild population.     

High-Resolution Genetic Mapping of Complex Traits from a Combined Analysis of F2 and Advanced Intercross Mice

Genetic influences on anxiety disorders are well documented; however, the specific genes underlying these disorders remain largely unknown. To identify quantitative trait loci (QTL) for conditioned fear and open field behavior, we used an F₂ intercross (n = 490) and a 34th-generation advanced intercross line (AIL) (n = 687) from the LG/J and SM/J inbred mouse strains. The F₂ provided strong support for several QTL, but within wide chromosomal regions. The AIL yielded much narrower QTL, but the results were less statistically significant, despite the larger number of mice. Simultaneous analysis of the F₂ and AIL provided strong support for QTL and within much narrower regions. We used a linear mixed-model approach, implemented in the program QTLRel, to correct for possible confounding due to familial relatedness. Because we recorded the full pedigree, we were able to empirically compare two ways of accounting for relatedness: using the pedigree to estimate kinship coefficients and using genetic marker estimates of “realized relatedness.” QTL mapping using the marker-based estimates yielded more support for QTL, but only when we excluded the chromosome being scanned from the marker-based relatedness estimates. We used a forward model selection procedure to assess evidence for multiple QTL on the same chromosome. Overall, we identified 12 significant loci for behaviors in the open field and 12 significant loci for conditioned fear behaviors. Our approach implements multiple advances to integrated analysis of F₂ and AILs that provide both power and precision, while maintaining the advantages of using only two inbred strains to map QTL.     

Using microsatellite diversity in wild Anegada iguanas (<Emphasis Type="Italic">Cyclura pinguis</Emphasis>) to establish relatedness in a captive breeding group of this critically endangered species

Awareness of the genealogical relationships between founder animals in captive breeding programs is essential for the selection of mating pairs that maintain genetic diversity. If captive founder relationships are unknown they can be inferred using genetic data from wild populations. Here, we report the results of such an analysis for six Cyclura pinguis (Sauria: Iguanidae) acquired as adults in 1999 by the San Diego Zoo Institute for Conservation Research to begin a captive breeding program for this critically endangered species. The six founder animals were reportedly hatched in captivity from eggs collected on Anegada in 1985. No records exist, however, as to where on Anegada the eggs were collected or from how many nests they originated. To assist determination of genealogical relationships, we genotyped the six captive founders, their first six offspring, and 33 wild adult iguanas from Anegada at 23 informative microsatellite loci. With these data, we estimated allele frequencies among the wild samples and then estimated the relatedness of the captive population. Using likelihood inference, we determined that three closely related pairs exist among the six captive founders and that each pair is not closely related to the other two. In addition, we were able to assign parentage for all six of the founders’ offspring tested, one of which had been previously misdiagnosed. Using the assigned parentage and inferred relatedness of the six founders, we calculated mean kinship for each of the six founders and their five living offspring. Finally, based on the allelic diversity of the wild iguanas sampled, we conclude that the C. pinguis population on Anegada is not excessively inbred; however, further investigation is warranted.     

Molecular and pedigree measures of relatedness provide similar estimates of inbreeding depression in a bottlenecked population

Individual‐based estimates of the degree of inbreeding or parental relatedness from pedigrees provide a critical starting point for studies of inbreeding depression, but in practice wild pedigrees are difficult to obtain. Because inbreeding increases the proportion of genomewide loci that are identical by descent, inbreeding variation within populations has the potential to generate observable correlations between heterozygosity measured using molecular markers and a variety of fitness related traits. Termed heterozygosity‐fitness correlations (HFCs), these correlations have been observed in a wide variety of taxa. The difficulty of obtaining wild pedigree data, however, means that empirical investigations of how pedigree inbreeding influences HFCs are rare. Here, we assess evidence for inbreeding depression in three life‐history traits (hatching and fledging success and juvenile survival) in an isolated population of Stewart Island robins using both pedigree‐ and molecular‐derived measures of relatedness. We found results from the two measures were highly correlated and supported evidence for significant but weak inbreeding depression. However, standardized effect sizes for inbreeding depression based on the pedigree‐based kin coefficients (k) were greater and had smaller standard errors than those based on molecular genetic measures of relatedness (RI), particularly for hatching and fledging success. Nevertheless, the results presented here support the use of molecular‐based measures of relatedness in bottlenecked populations when information regarding inbreeding depression is desired but pedigree data on relatedness are unavailable.     

Estimating average within-group relatedness from DNA fingerprints

We present a detailed analysis of a regression and an equivalent band-sharing technique for estimating mean within-group relatedness from DNA fingerprint data. First, we show that the relatedness estimate obtained from this technique tends to underestimate the pedigree relatedness in random-mating populations, but that the degree of underestimation will be small unless DNA fingerprint similarity among unrelated individuals is unusually high. We derive the exact asymptotic downward bias as a function of allele frequency and probabilities of identity by descent. We suggest some methods for removing the small downward bias in the regression estimate. Secondly, we show that the technique yields accurate estimates of parent-offspring, full, and half-sibling relatedness in natural populations of the red-winged blackbird Agelaius phoeniceus and indigo bunting Passerina cyanea.     

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.     

Partial genotyping at polymorphic markers can improve heritability estimates in sibling groups

Accurate estimates of heritability () are necessary to assess adaptive responses of populations and evolution of fitness‐related traits in changing environments. For plants, estimates generally rely on maternal progeny designs, assuming that offspring are either half‐sibs or unrelated. However, plant mating systems often depart from half‐sib assumptions, this can bias estimates. Here, we investigate how to accurately estimate in nonmodel species through the analysis of sibling designs with a moderate genotyping effort. We performed simulations to investigate how microsatellite marker information available for only a subset of offspring can improve estimates based on maternal progeny designs in the presence of nonrandom mating, inbreeding in the parental population or maternal effects. We compared the basic family method, considering or not adjustments based on average relatedness coefficients, and methods based on the animal model. The animal model was used with average relatedness information, or with hybrid relatedness information: associating one‐generation pedigree and family assumptions, or associating one‐generation pedigree and average relatedness coefficients. Our results highlighted that methods using marker‐based relatedness coefficients performed as well as pedigree‐based methods in the presence of nonrandom mating (i.e. unequal male reproductive contributions, selfing), offering promising prospects to investigate in situ heritabilities in natural populations. In the presence of maternal effects, only the use of pairwise relatednesses through pedigree information improved the accuracy of estimates. In that case, the amount of father‐related offspring in the sibling design is the most critical. Overall, we showed that the method using both one‐generation pedigree and average relatedness coefficients was the most robust to various ecological scenarios.     

Quality and quantity of genetic relatedness data affect the analysis of social structure

Kinship plays a fundamental role in the evolution of social systems and is considered a key driver of group living. To understand the role of kinship in the formation and maintenance of social bonds, accurate measures of genetic relatedness are critical. Genotype‐by‐sequencing technologies are rapidly advancing the accuracy and precision of genetic relatedness estimates for wild populations. The ability to assign kinship from genetic data varies depending on a species’ or population's mating system and pattern of dispersal, and empirical data from longitudinal studies are crucial to validate these methods. We use data from a long‐term behavioural study of a polygynandrous, bisexually philopatric marine mammal to measure accuracy and precision of parentage and genetic relatedness estimation against a known partial pedigree. We show that with moderate but obtainable sample sizes of approximately 4,235 SNPs and 272 individuals, highly accurate parentage assignments and genetic relatedness coefficients can be obtained. Additionally, we subsample our data to quantify how data availability affects relatedness estimation and kinship assignment. Lastly, we conduct a social network analysis to investigate the extent to which accuracy and precision of relatedness estimation improve statistical power to detect an effect of relatedness on social structure. Our results provide practical guidance for minimum sample sizes and sequencing depth for future studies, as well as thresholds for post hoc interpretation of previous analyses.     

A comparison of pedigree‐ and DNA‐based measures for identifying inbreeding depression in the critically endangered Attwater's Prairie‐chicken

The primary goal of captive breeding programmes for endangered species is to prevent extinction, a component of which includes the preservation of genetic diversity and avoidance of inbreeding. This is typically accomplished by minimizing mean kinship in the population, thereby maintaining equal representation of the genetic founders used to initiate the captive population. If errors in the pedigree do exist, such an approach becomes less effective for minimizing inbreeding depression. In this study, both pedigree‐ and DNA‐based methods were used to assess whether inbreeding depression existed in the captive population of the critically endangered Attwater's Prairie‐chicken (Tympanuchus cupido attwateri), a subspecies of prairie grouse that has experienced a significant decline in abundance and concurrent reduction in neutral genetic diversity. When examining the captive population for signs of inbreeding, variation in pedigree‐based inbreeding coefficients (f_pedigree) was less than that obtained from DNA‐based methods (f_DNA). Mortality of chicks and adults in captivity were also positively correlated with parental relatedness (r_DNA) and f_DNA, respectively, while no correlation was observed with pedigree‐based measures when controlling for additional variables such as age, breeding facility, gender and captive/release status. Further, individual homozygosity by loci (HL) and parental r_DNA values were positively correlated with adult mortality in captivity and the occurrence of a lethal congenital defect in chicks, respectively, suggesting that inbreeding may be a contributing factor increasing the frequency of this condition among Attwater's Prairie‐chickens. This study highlights the importance of using DNA‐based methods to better inform management decisions when pedigrees are incomplete or errors may exist due to uncertainty in pairings.     

Computer software for performing likelihood tests of pedigree relationship using genetic markers

Molecular techniques are making ever more genetic markers available for use in parentage assignment, and measures of relatedness. We present a program, Kinship, designed to use likelihood techniques to test for any non-inbred pedigree relationship between pairs of individuals, using single-locus codominant genetic markers. Kinship calculates the likelihood that each pair of individuals in a data set are related by a given pedigree hypothesis, and likelihood ratios for any pair of hypotheses. The program also uses a simulation routine to attach statistical significance to its results.     

Random amplified polymorphic DNA and pedigree relationships in spring barley

Summary We investigated random amplified polymorphic DNA (RAPD) in 27 inbred barley lines with varying amounts of common ancestry and in 20 doubled-haploid (DH) lines from a biparental cross. Of 33 arbitrary 10 base primers that were tested, 19 distinguished a total of 31 polymorphisms. All polymorphisms were scored as dominant genetic markers except for 1, where Southern analysis indicated the presence of two codominant amplification products. The inheritance of 19 RAPD polymorphisms and one morphological trait was studied in the DH lines. There was no evidence for segregation distortion, but a group of four tightly linked loci was detected. The frequencies of RAPD polymorphism in pairs of inbred lines were used to compute values of genetic distance (d), which were compared to kinship coefficients (r) between the same pairs of lines. A linear relationship between r and d was evident, but low values of r gave poor predictions of d. Cluster analysis showed that groups of inbred lines based on r were similar to those based on d with some notable exceptions. RAPD markers can be used to gain information about genetic similarities or differences that are not evident from pedigree information.     

SSRs,SNPs and DArTs comparison on estimation of relatedness and genetic parameters’ precision from a small half-sib sample population of <Emphasis Type="Italic">Eucalyptus grandis</Emphasis>

Simple sequence repeats (SSR) are the most widely used molecular markers for relatedness inference due to their multi-allelic nature and high informativeness. However, there is a growing trend toward using high-throughput and inter-specific transferable single-nucleotide polymorphisms (SNP) and Diversity Arrays Technology (DArT) in forest genetics owing to their wide genome coverage. We compared the efficiency of 15 SSRs, 181 SNPs and 2816 DArTs to estimate the relatedness coefficients, and their effects on genetic parameters’ precision, in a relatively small data set of an open-pollinated progeny trial of Eucalyptus grandis (Hill ex Maiden) with limited relationship from the pedigree. Both simulations and real data of Eucalyptus grandis were used to study the statistical performance of three relatedness estimators based on co-dominant markers. Relatedness estimates in pairs of individuals belonging to the same family (related) were higher for DArTs than for SNPs and SSRs. DArTs performed better compared to SSRs and SNPs in estimated relatedness coefficients in pairs of individuals belonging to different families (unrelated) and showed higher ability to discriminate unrelated from related individuals. The likelihood-based estimator exhibited the lowest root mean squared error (RMSE); however, the differences in RMSE among the three estimators studied were small. For the growth traits, heritability estimates based on SNPs yielded, on average, smaller standard errors compared to those based on SSRs and DArTs. Estimated relatedness in the realized relationship matrix and heritabilities can be accurately inferred from co-dominant or sufficiently dense dominant markers in a relatively small E. grandis data set with shallow pedigree.     

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.     

INBREEDING INFLUENCES WITHIN‐BROOD HETEROZYGOSITY‐FITNESS CORRELATIONS (HFCS) IN AN ISOLATED PASSERINE POPULATION

Molecular estimates of inbreeding may be made using genetic markers such as microsatellites, however the interpretation of resulting heterozygosity‐fitness correlations (HFCs) with respect to inbreeding depression is not straightforward. We investigated the relationship between pedigree‐determined inbreeding coefficients (f) and HFCs in a closely monitored, reintroduced population of Stewart Island robins (Petroica australis rakiura) on Ulva Island, New Zealand. Using a full sibling design, we focused on differences in juvenile survival associated specifically with individual sibling variation in standardized multilocus heterozygosity (SH) when expected f was identical. We found that within broods, siblings with higher SH at microsatellite loci experienced a higher probability of juvenile survival. This effect, however, was detected primarily within broods that experienced inbreeding or when inbreeding had occurred in their pedigree histories (i.e., at the parents’ level). Thus we show, for the first time in a wild population, that the strength of an HFC is partially dependent on the presence of inbreeding events in the recent pedigree history. Our results illustrate the importance of realized effects of inbreeding on genetic variation and fitness and the value of full‐sibling designs for the study of HFCs in the context of small, inbred populations.     

Quantifying inbreeding avoidance through extra-pair reproduction

Extra-pair reproduction is widely hypothesized to allow females to avoid inbreeding with related socially paired males. Consequently, numerous field studies have tested the key predictions that extra-pair offspring are less inbred than females’ alternative within-pair offspring, and that the probability of extra-pair reproduction increases with a female's relatedness to her socially paired male. However, such studies rarely measure inbreeding or relatedness sufficiently precisely to detect subtle effects, or consider biases stemming from failure to observe inbred offspring that die during early development. Analyses of multigenerational song sparrow (Melospiza melodia) pedigree data showed that most females had opportunity to increase or decrease the coefficient of inbreeding of their offspring through extra-pair reproduction with neighboring males. In practice, observed extra-pair offspring had lower inbreeding coefficients than females’ within-pair offspring on average, while the probability of extra-pair reproduction increased substantially with the coefficient of kinship between a female and her socially paired male. However, simulations showed that such effects could simply reflect bias stemming from inbreeding depression in early offspring survival. The null hypothesis that extra-pair reproduction is random with respect to kinship therefore cannot be definitively rejected in song sparrows, and existing general evidence that females avoid inbreeding through extra-pair reproduction requires reevaluation given such biases.     

How many SNPs are enough?

Since the days of allozyme analysis, we have been enamored with the idea that if we just had enough polymorphic mendelian loci, we could gauge the inbreeding level of individuals by measuring heterozygosity and simultaneously measure the degree of genetic relatedness between pairs of individuals. Given Mendel’s Laws, we have always known that we would need numerous independently segregating loci to achieve any reasonable degree of accuracy. Santure et al. (2010, this issue) use a 771 marker SNP panel to assess heterozygosity levels and to assess pairwise relatedness, and compare both with theoretical expectations obtained from a carefully recorded pedigree of a zebra finch breeding colony, as a function of increasing numbers of SNP markers. They also compare the SNP results with those from a 20‐locus microsatellite panel, showing that adding SNPs to a fairly large microsatellite panel improves accuracy, but given an existing panel of 125 SNPs, little is to be gained by adding microsatellites. They show that the accuracy available for estimating individual levels of inbreeding is somewhat limited. They also show that the average pairwise relatedness measures bracket pedigree relationship very nicely, but the variances for individual pairs remain substantial, even with a very large panel.     

Measuring genetic relatedness in natural populations: Methodology

The estimation of relatedness within social groups, such as the colonies of a population of social insects, is an important field for evaluating hypotheses concerning the evolution and maintenance of social behaviour. The methodology of this estimation from genetic data in the absence of pedigree information has been poorly understood; we develop this methodology for b, the regression coefficient of relatedness, and discuss its applications. Both b and G (the pedigree coefficient of relatedness) are potentially asymmetric coefficients, whereas φ, r, and F_ST are necessarily symmetric. We develop an estimator for b suitable for small samples, and also one for standard deviation, and examine the properties of both using sampling simulations. The b estimator returns values slightly below E(b), and the standard deviation estimator yields conservative confidence intervals. A comparative study of b and F_ST shows that, given the same set of data, b is estimated with greater reliability than is F_ST. As is the case for F_ST, b can be used to examine population structure at various levels, and b possesses the advantage of an estimator for its standard error, which can also be used to test for heterogeneity among the loci surveyed. The actual numbers of identical genes held in common by interacting individuals, and not simply their proportions, need to be considered in using coefficients of relatedness in inclusive fitness calculations. This necessity is handled by the weighted coefficients of relatedness, G′ and b′, which have been referred to in the literature as r (as have most relatedness measures).     

Genetic assessment of a conservation breeding program of the houbara bustard (Chlamydotis undulata undulata) in Morocco,based on pedigree and molecular analyses

Protection and restoration of species in the wild may require conservation breeding programs under genetic management to minimize deleterious effects of genetic changes that occur in captivity, while preserving populations' genetic diversity and evolutionary resilience. Here, through interannual pedigree analyses, we first assessed the efficiency of a 21-year genetic management, including minimization of mean kinship, inbreeding avoidance, and regular addition of founders, of a conservation breeding program targeting on Houbara bustard (Chlamydotis undulata undulata) in Morocco. Secondly, we compared pedigree analyses, the classical way of assessing and managing genetic diversity in captivity, to molecular analyses based on seven microsatellites. Pedigree-based results indicated an efficient maintenance of the genetic diversity (99% of the initial genetic diversity retained) while molecular-based results indicated an increase in allelic richness and an increase in unbiased expected heterozygosity across time. The pedigree-based average inbreeding coefficient F remained low (between 0.0004 and 0.003 in 2017) while the proportion of highly inbred individuals (F > .1) decreased over time and reached 0.2% in 2017. Furthermore, pedigree-based F and molecular-based individual multilocus heterozygosity were weakly negatively correlated, (Pearson's r = −.061 when considering all genotyped individuals), suggesting that they cannot be considered as alternatives, but rather as complementary sources of information. These findings suggest that a strict genetic monitoring and management, based on both pedigree and molecular tools can help mitigate genetic changes and allow to preserve genetic diversity and evolutionary resilience in conservation breeding programs.     

Measures of Relatedness

Measures of genetic relatedness are essential to models of evolution by kin selection and determinations of inclusive fitness. Under a kin selection paradigm, individuals are expected to distribute actions influencing the fitness of relatives based on the relatedness of these relatives. In addition, it is necessary to have an accurate measure of relatedness to estimate heritability (h²) of phenotypic characters and to predict the efficacy of selection. Relatedness is often defined as the genotypic correlation between individuals. Assessed on the basis of common ancestry, relatedness can only be determined sensu strictu from pedigree analysis. Recent methodological and statistical advances allow the estimation of relatedness from allele frequency data. Many coefficients of relatedness can be found in the literature; I review and evaluate these, with emphasis on situations for which each is appropriate.