COANCESTRY: a program for simulating, estimating and analysing relatedness and inbreeding coefficients

The software package COANCESTRY implements seven relatedness estimators and three inbreeding estimators to estimate relatedness and inbreeding coefficients from multilocus genotype data. Two likelihood estimators that allow for inbred individuals and account for genotyping errors are for the first time included in this user-friendly program for PCs running Windows operating system. A simulation module is built in the program to simulate multilocus genotype data of individuals with a predefined relationship, and to compare the estimators and the simulated relatedness values to facilitate the selection of the best estimator in a particular situation. Bootstrapping and permutations are used to obtain the 95% confidence intervals of each relatedness or inbreeding estimate, and to test the difference in averages between groups.     

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.     

Estimation of pairwise relatedness between individuals and characterization of isolation-by-distance processes using dominant genetic markers

A new estimator of the pairwise relatedness coefficient between individuals adapted to dominant genetic markers is developed. This estimator does not assume genotypes to be in Hardy-Weinberg proportions but requires a knowledge of the departure from these proportions (i.e. the inbreeding coefficient). Simulations show that the estimator provides accurate estimates, except for some particular types of individual pairs such as full-sibs, and performs better than a previously developed estimator. When comparing marker-based relatedness estimates with pedigree expectations, a new approach to account for the change of the reference population is developed and shown to perform satisfactorily. Simulations also illustrate that this new relatedness estimator can be used to characterize isolation by distance within populations, leading to essentially unbiased estimates of the neighbourhood size. In this context, the estimator appears fairly robust to moderate errors made on the assumed inbreeding coefficient. The analysis of real data sets suggests that dominant markers (random amplified polymorphic DNA, amplified fragment length polymorphism) may be as valuable as co-dominant markers (microsatellites) in studying microgeographic isolation-by-distance processes. It is argued that the estimators developed should find major applications, notably for conservation biology.     

Ex situ population management in the absence of pedigree information

For captive breeding to play a significant role in conservation, ex situ populations must be scientifically managed to meet objective goals for retaining representative genetic variation. Imperfect genealogical information requires fundamental assumptions to be made that may bias downstream measures of genetic importance, upon which management decisions are based. The impacts of such assumptions are most pronounced within breeding programmes characterized by a high proportion of individuals of unknown ancestry, as exemplified by the large captive population of the St Vincent parrot (Amazona guildingii). The degree to which microsatellite-based estimates of relatedness may improve upon the assumptions of conventional pedigree-based management was investigated using genotypic data collected at eight microsatellite loci and two marker-based relatedness estimators. The measure, rxyLR, was found to explain the highest amount of variation in true relatedness. Integration of pairwise estimates of founder relatedness with studbook data transformed current understanding of the relatedness structure of the A. guildingii population from two subgroups characterized by a high and low degree of relatedness, respectively, to a situation where all 72 individuals are prioritized for breeding according to their estimated mean kinships. Furthermore, the discovery of opposing, directional bias exhibited by rxyLR and rxyQG in assigning dyads to a given relationship category suggests that an approach that utilizes a combination of pairwise relatedness estimators may provide the most genetic information for balancing the dual considerations of maximizing gene diversity and minimizing inbreeding in developing breeding recommendations.     

Sibship within samples of brown trout (<Emphasis Type="Italic">Salmo trutta</Emphasis>) and implications for supportive breeding

We analysed family relationships among brown trout from two small tributary populations that have been suggested as a source of individuals for supportive breeding, using variation at eight microsatellite loci. As a control, we analysed a sample of supposedly unrelated individuals representing a large anadromous population, and we simulated unrelated individuals based on the allelic distributions in all three samples. Two different approaches were used: (1) pairwise estimates of relatedness between individuals and (2) a method for partitioning individuals into half-sib and full-sib families. The anadromous population did not show evidence of a significant number of closely related individuals. In both tributary populations, however, the distributions of pairwise relatedness estimates suggested the presence of several related individuals, and sibship reconstruction suggested fewer families consisting of more individuals than were observed for the simulated individuals. The expected increase of inbreeding coefficient in the two samples due to family structure was 0.026 and 0.030 respectively. Moreover, tests for recent bottlenecks yielded significant outcomes in both populations suggesting a history of low effective population sizes. Depending on the effective population size of captive spawners and past effective population sizes in the populations it could be beneficial to conduct sib-avoidance matings, though this cannot eliminate inbreeding but only delay it. Alternatively, individuals from different populations could be crossed. Sibship reconstruction provided the clearest evidence for family structure, but pairwise relatedness is the best measure for designing mating schemes, as it allows for mating as unrelated individuals as possible rather than just avoiding mating between sibs.     

Patterns of relatedness and parentage in an asocial, polyandrous striped hyena population

We investigated patterns of relatedness and reproduction in a population of striped hyenas in which individuals are behaviourally solitary but form polyandrous spatial groups consisting of one adult female and multiple adult males. Group-mate males were often close relatives, but were unrelated or distantly related in some cases, indicating that male coalitions are not strictly a result of philopatry or dispersal with cohorts of relatives. Most male-female pairs within spatial groups were unrelated or only distantly related. Considering patterns of relatedness between groups, relatedness was significantly higher among adult males living in non-neighbouring ranges than among neighbouring males. Mean relatedness among male-female dyads was highest for group-mates, but relatedness among non-neighbouring males and females was also significantly higher than among dyads of opposite-sex neighbours. Female-female relatedness also increased significantly with increasing geographic separation. These unusual and unexpected patterns may reflect selection to settle in a nonadjacent manner to reduce inbreeding and/or competition among relatives for resources (both sexes), or mates (males). Finally, resident males fathered the majority of the resident female's cubs, but extra-group paternity was likely in 31% of the cases examined, and multiple paternity was likely in half of the sampled litters.     

No evidence of inbreeding avoidance despite demonstrated survival costs in a polygynous rodent

Individuals are generally predicted to avoid inbreeding because of detrimental fitness effects. However, several recent studies have shown that limited inbreeding is tolerated by some vertebrate species. Here, we examine the costs and benefits of inbreeding in a largely polygynous rodent, the yellow-bellied marmot (Marmota flaviventris). We use a pedigree constructed from 8 years of genetic data to determine the relatedness of all marmots in our study population and examine offspring survival, annual male reproductive success, relatedness between breeding pairs and the effects of group composition on likelihood of male reproduction to assess inbreeding in this species. We found decreased survival in inbred offspring, but equal net reproductive success among males that inbred and those that avoided it. Relatedness between breeding pairs was greater than that expected by chance, indicating that marmots do not appear to avoid breeding with relatives. Further, male marmots do not avoid inbreeding: males mate with equal frequency in groups composed of both related and unrelated females and in groups composed of only female relatives. Our results demonstrate that inbreeding can be tolerated in a polygynous species if the reproductive costs of inbreeding are low and individuals that mate indiscriminately do not suffer decreased reproductive success.     

Estimating pairwise relatedness from dominant genetic markers

Knowledge of the genetic relatedness between a pair of individuals is important in many research areas of quantitative genetics, conservation genetics, evolution and ecology. Many estimators have been developed to estimate such pairwise relatedness (r) using codominant markers, such as microsatellites and enzymes. In contrast, only two estimators are proposed to use dominant markers, such as random amplified polymorphic DNAs (RAPDs) and amplified fragment length polymorphisms (AFLPs), in relatedness inference. They are both biased estimators, and their statistical properties and robustness to the sampling errors in allele frequency have not been investigated. In this short paper, I propose two new pairwise relatedness estimators for dominant markers, and compare them in precision, accuracy and robustness to sampling with the two previous estimators using simulations. It was found that the new estimator based on the least squares approach is unbiased when allele frequencies are known or estimated from a sample without correcting for sampling effects. It has, however, a low precision and as a result, an intermediate overall performance among the four estimators in terms of the mean squared deviation (MSD) of estimates from actual values of r. The new estimator based on a similarity index is slightly biased but has generally the lowest MSD among the four estimators compared, regardless of the number of loci, type of actual relationships, allele frequencies known or estimated from samples. Simulations also show that the confidence intervals estimated by bootstrapping are appropriate for different estimators provided that the number of loci used in the estimation is not small.     

Inference of relationships in population data using identity-by-descent and identity-by-state

It is an assumption of large, population-based datasets that samples are annotated accurately whether they correspond to known relationships or unrelated individuals. These annotations are key for a broad range of genetics applications. While many methods are available to assess relatedness that involve estimates of identity-by-descent (IBD) and/or identity-by-state (IBS) allele-sharing proportions, we developed a novel approach that estimates IBD0, 1, and 2 based on observed IBS within windows. When combined with genome-wide IBS information, it provides an intuitive and practical graphical approach with the capacity to analyze datasets with thousands of samples without prior information about relatedness between individuals or haplotypes. We applied the method to a commonly used Human Variation Panel consisting of 400 nominally unrelated individuals. Surprisingly, we identified identical, parent-child, and full-sibling relationships and reconstructed pedigrees. In two instances non-sibling pairs of individuals in these pedigrees had unexpected IBD2 levels, as well as multiple regions of homozygosity, implying inbreeding. This combined method allowed us to distinguish related individuals from those having atypical heterozygosity rates and determine which individuals were outliers with respect to their designated population. Additionally, it becomes increasingly difficult to identify distant relatedness using genome-wide IBS methods alone. However, our IBD method further identified distant relatedness between individuals within populations, supported by the presence of megabase-scale regions lacking IBS0 across individual chromosomes. We benchmarked our approach against the hidden Markov model of a leading software package (PLINK), showing improved calling of distantly related individuals, and we validated it using a known pedigree from a clinical study. The application of this approach could improve genome-wide association, linkage, heterozygosity, and other population genomics studies that rely on SNP genotype data.     

Comparison of marker-based pairwise relatedness estimators on a pedigreed plant population

Several estimators have been proposed that use molecular marker data to infer the degree of relatedness for pairs of individuals. The objective of this study was to evaluate the performance of seven estimators when applied to marker data of a set of 33 key individuals from a large complex apple pedigree. The evaluation considered different scenarios of allele frequencies and different numbers of marker loci. The method of moments estimators were Similarity, Queller-Goodknight, Lynch-Ritland and Wang. The maximum likelihood estimators were Thompson, Anderson-Weir and Jacquard. The pedigree-based coancestry coefficients were taken as the point of reference in calculating correlations and root mean square error (RMSE). The marker data comprised 86 multi-allelic SSR markers on 17 linkage groups, covering 11 Morgans. Additionally, we simulated 10 datasets conditional on the real pedigree to support the results on the real dataset. None of the estimators outperformed the others. Knowledge of allele frequencies appeared to be the most influential, i.e., the highest correlations and lowest RMSE were found when frequencies from the founder population were available. When equal allele frequencies were used, all estimators resulted in very similar, but on average lower, correlations. The use of allele frequencies estimated from the set of 33 individuals gave, on average, the poorest results. The maximum likelihood estimators and the Lynch-Ritland estimator were the most sensitive to allele frequencies. The results from the simulation study fully supported the trends in results of the real dataset. This study indicated that high correlations (up to 0.90) and small RMSE (below 0.03), may be obtained when population allelic frequencies are available. In this scenario, the performances of the various estimators were similar, but seemed to favor the maximum likelihood estimators. In the absence of reliable allele frequencies the method of moments estimators were shown to be more robust. The number of marker loci influenced the average performance of the estimators; however, the ranking was not affected. Correlations up to 0.80 were obtained when two markers per chromosome and appropriate allele frequencies were available. Adding more markers to the current dataset may lead to marginal improvements.     

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.     

DOES GENETIC RELATEDNESS OF MATES INFLUENCE COMPETITIVE FERTILIZATION SUCCESS IN GUPPIES?

A growing number of studies highlight the nontransitive properties of ejaculates when they are in competition to fertilize a female's eggs. Increasingly, these studies suggest that postcopulatory processes act as a filter against sperm from closely related males or those with similar genotypes, limiting the deleterious effects of inbreeding on offspring fitness. We investigated the potential for such postcopulatory mechanisms of inbreeding avoidance in the guppy (Poecilia reticulata), a promiscuous livebearing fish. We used artificial insemination as a method of delivering to a female the combined ejaculates from a first cousin (relatedness coefficient r = 0.125) and an unrelated male. This method of sperm delivery controls behavioral processes of pre- and postcopulatory female choice, which can bias paternity toward unrelated males. Our genetic analysis revealed no effect of parental relatedness on paternity outcomes. The observed mean paternity share for related males (0.47) and associated variance did not differ significantly from an expected binomial distribution that assumes no biased use of sperm with respect to relatedness (0.5). Although our data provide no evidence for postcopulatory mechanisms of inbreeding avoidance, the ability of female guppies to influence ejaculate transfer and retention offers an alternative and easily testable mechanism of inbreeding avoidance in this species.     

Benefiting friends or dominants: prosocial choices mainly depend on rank position in long-tailed macaques (<Emphasis Type="Italic">Macaca fascicularis</Emphasis>)

Long-term observational studies in a number of animal species suggest that exchange patterns of social acts depend on long-term emotional bonds. Therefore, it is expected that the frequency of prosocial behavior will depend on the strength of such a bond. In this study we tested whether variation in relationship quality among unrelated individuals, i.e., “friends” and “nonfriends,” is predictive of the prosocial behavior of long-tailed macaques in two experiments. First, we related relationship quality to prosociality in a dyadic prosociality test, and second, we gave subjects the choice to give to either a friend or a nonfriend in a triadic choice test. We show that prosocial behavior of long-tailed macaques in the dyadic test is not related to relationship quality. When given the choice to give to either a friend or a nonfriend in the triadic test, there is a minor indication that long-tailed macaques show a preference to give to their friends, yet this indication is neither significant nor consistent. In contrast, subordinate long-tailed macaques make a more “competitive” choice and avoid giving to the individual closest in rank. Therefore, in the short-term situation of experimental tests, prosocial behavior of long-tailed macaques seems unaffected by the relationship quality of the dyad/triad tested, and the relative dominance position of these dyads/triads seems to have a much stronger effect on their prosocial behavior.     

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.     

Relatedness measured by oligonucleotide probe DNA fingerprints and an estimate of the mating system of Sea Lavender (Limonium carolinianum)

Using DNA fingerprint markers within species and populations of wild plants requires information on the relationship between fingerprint similarity and relatedness. We identified a hypervariable marker based on oliog(GATA)₄-hybridization of DpnII-cut genomic DNA from Sea Lavender (Limonium carolinianum). Banding patterns were somatically stable and highly variable among unrelated individuals. Band molecular-weight sizing errors (as a percent of band molecular weight) were estimated at 0.44%±0.003 within gels and 0.76%±0.964 between gels. Band sizing errors defined a 99% confidence bin of ±0.95% (1.90% total) of molecular weight. Band-sharing estimates were based on this bin size and on variance estimates that compensate for non-independent comparisons. Band-sharing among nine unrelated individuals () was 0.198±0.O11. Experimental pollinations designed to produce selfed, fulland half-sib progeny groups led to five selfed progeny groups and no outcrossed progeny (mean band-sharing, ovS=0.468±0.074). A linear regression between band-sharing (S) and relatedness (r) assuming 17% inbreeding was r=0.006+0.914*S (R²=0.973) and established the maximum amount of inbreeding. ovS(0.392±0.022) estimated from wild pollinated seeds from four maternal families was intermediate to unrelated individuals and experimental selfed progeny, giving evidence for mixed mating in wild plants. More extensive plant pedigrees with known levels of inbreeding will be needed to measure variation in the relationship between S and r among populations and families.     

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 .     

TECHNICAL ADVANCES: A maximum-likelihood relatedness estimator allowing for negative relatedness values

Previously reported maximum-likelihood pairwise relatedness (r) estimator of Thompson and Milligan (M) was extended to allow for negative r estimates under the regression interpretation of r. This was achieved by establishing the equivalency of the likelihoods used in the kinship program and the likelihoods of Thompson. The new maximum-likelihood (ML) estimator was evaluated by Monte Carlo simulations. It was found that the new ML estimator became unbiased significantly faster compared to the original M estimator when the amount of genotype information was increased. The effects of allele frequency estimation errors on the new and existing relatedness estimators were also considered.     

Sex-biased natal dispersal and inbreeding avoidance in American black bears as revealed by spatial genetic analyses

We tested the hypothesis that sex-biased natal dispersal reduces close inbreeding in American black bears, a solitary species that exhibits nearly complete male dispersal and female philopatry. Using microsatellite DNA and spatial data from reproductively mature bears (>or= 4 years old), we examined the spatial genetic structure of two distinct populations in New Mexico from 1993 to 2000. As predicted, relatedness (r) and the frequency of close relationships (parent-offspring or full siblings) decreased with distance among female dyads, but little change was observed among male or opposite-sex dyads. Neighbouring females were more closely related than neighbouring males. The potential for inbreeding was low. Most opposite-sex pairs that lived sufficiently close to facilitate mating were unrelated, and few were close relatives. We found no evidence that bears actively avoided inbreeding in their selection of mates from this nearby pool, as mean r and relationship frequencies did not differ between potential and actual mating pairs (determined by parentage analysis). These basic patterns were apparent in both study areas despite a nearly two-fold difference in density. However, the sex bias in dispersal was less pronounced in the lower-density area, based on proportions of bears with male and female relatives residing nearby. This result suggests that male bears may respond to reduced competition by decreasing their rate or distance of dispersal. Evidence supports the hypothesis that inbreeding avoidance is achieved by means of male-biased dispersal but also indicates that competition (for mates or resources) modifies dispersal patterns.     

On the use of large marker panels to estimate inbreeding and relatedness: empirical and simulation studies of a pedigreed zebra finch population typed at 771 SNPs

In recent years there has been a dramatic increase in the availability of high density genetic marker data for both model and non‐model organisms. A potential application of these data is to infer relatedness in the absence of a complete pedigree. Using a marker panel of 771 SNPs genotyped in three generations of an extensive zebra finch pedigree, correlations between pedigree relatedness and seven marker‐based estimates of relatedness were examined, as was the relationship between heterozygosity and inbreeding. Although marker‐based and pedigree relatedness were highly correlated, the variance in estimated relatedness was high. Further, the correlation between heterozygosity and inbreeding was weak, even though mean inbreeding coefficient is typical of that seen in wild vertebrate pedigrees; the weak relationship was in part due to the small variance in inbreeding in the pedigree. Our data suggest that using marker information to reconstruct the pedigree, and then calculating relatedness from the pedigree, is likely to give more accurate relatedness estimates than using marker‐based estimators directly.     

Friends and family: A software program for identification of unrelated individuals from molecular marker data

The identification of related and unrelated individuals from molecular marker data is often difficult, particularly when no pedigree information is available and the data set is large. High levels of relatedness or inbreeding can influence genotype frequencies and thus genetic marker evaluation, as well as the accurate inference of hidden genetic structure. Identification of related and unrelated individuals is also important in breeding programmes, to inform decisions about breeding pairs and translocations. We present Friends and Family, a Windows executable program with a graphical user interface that identifies unrelated individuals from a pairwise relatedness matrix or table generated in programs such as coancestry and genalex . Friends and Family outputs a list of samples that are all unrelated to each other, based on a user‐defined relatedness cut‐off value. This unrelated data set can be used in downstream analyses, such as marker evaluation or inference of genetic structure. The results can be compared to that of the full data set to determine the effect related individuals have on the analyses. We demonstrate one of the applications of the program: how the removal of related individuals altered the Hardy–Weinberg equilibrium test outcome for microsatellite markers in an empirical data set. Friends and Family can be obtained from https://github.com/DeondeJager/Friends-and-Family .