Methods of parentage analysis in natural populations

The recent proliferation of hypervariable molecular markers has ushered in a surge of techniques for the analysis of parentage in natural and experimental populations. Consequently, the potential for meaningful studies of paternity and maternity is at an all-time high. However, the details and implementation of the multifarious techniques often differ in subtle ways that can influence the results of parentage analyses. Now is a good time to reflect on the available techniques and to consider their strengths and weaknesses. Here, we review the leading techniques in parentage analysis, with a particular emphasis on those that have been implemented in readily useable software packages. Our survey leads to some important insights with respect to the utility of the different approaches. This review should serve as a useful guide to anyone who wishes to embark on the study of parentage.     

On minimizing assignment errors and the trade‐off between false positives and negatives in parentage analysis

Genetic parentage analyses provide a practical means with which to identify parent–offspring relationships in the wild. In Harrison et al.'s study (2013a), we compare three methods of parentage analysis and showed that the number and diversity of microsatellite loci were the most important factors defining the accuracy of assignments. Our simulations revealed that an exclusion‐Bayes theorem method was more susceptible to false‐positive and false‐negative assignments than other methods tested. Here, we analyse and discuss the trade‐off between type I and type II errors in parentage analyses. We show that controlling for false‐positive assignments, without reporting type II errors, can be misleading. Our findings illustrate the need to estimate and report both the rate of false‐positive and false‐negative assignments in parentage analyses.     

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.     

Estimation of the number of SNP genetic markers required for parentage verification

Formulae were developed to compute exclusion probabilities for parentage confirmation for any number of diallelic markers under the assumption that the minor allele frequency (MAF) varied among markers, but has a uniform distribution. Three scenarios were analysed: a progeny with (1) a single putative parent; (2) two putative parents; and (3) one actual parent and one putative parent. Exclusion probabilities were computed for minimum values for the MAFs of 0.1, 0.2 and 0.3, and required either one or at least two conflicts for exclusion. The numbers of markers required to obtain 99% exclusion probabilities based on a single conflict for the three minimum MAFs were 54, 45 and 39 for scenario 1; 17, 16 and 15 for scenario 2; and 28, 25 and 24 for scenario 3. The requirement of at least two conflicts for exclusion increased the number of markers required by approximately 45% for all three scenarios and all three minimum MAFs. The results obtained by the analytical formulae were very close to results obtained by simulation and to values in the literature for specific marker sets.     

Relative accuracy of three common methods of parentage analysis in natural populations

Parentage studies and family reconstructions have become increasingly popular for investigating a range of evolutionary, ecological and behavioural processes in natural populations. However, a number of different assignment methods have emerged in common use and the accuracy of each may differ in relation to the number of loci examined, allelic diversity, incomplete sampling of all candidate parents and the presence of genotyping errors. Here, we examine how these factors affect the accuracy of three popular parentage inference methods (colony , famoz and an exclusion‐Bayes’ theorem approach by Christie (Molecular Ecology Resources, 2010a, 10, 115) to resolve true parent–offspring pairs using simulated data. Our findings demonstrate that accuracy increases with the number and diversity of loci. These were clearly the most important factors in obtaining accurate assignments explaining 75–90% of variance in overall accuracy across 60 simulated scenarios. Furthermore, the proportion of candidate parents sampled had a small but significant impact on the susceptibility of each method to either false‐positive or false‐negative assignments. Within the range of values simulated, colony outperformed FaMoz, which outperformed the exclusion‐Bayes’ theorem method. However, with 20 or more highly polymorphic loci, all methods could be applied with confidence. Our results show that for parentage inference in natural populations, careful consideration of the number and quality of markers will increase the accuracy of assignments and mitigate the effects of incomplete sampling of parental populations.     

Unexpected collective larval dispersal but little support for sweepstakes reproductive success in the highly dispersive brooding mollusc Crepidula fornicata

In many marine invertebrates, long‐distance dispersal is achieved during an extended pelagic larval phase. Although such dispersal should result in high gene flow over broad spatial scales, fine‐scale genetic structure has often been reported, a pattern attributed to interfamilial variance in reproductive success and limited homogenization during dispersal. To examine this hypothesis, the genetic diversity of dispersing larvae must be compared with the postdispersal stages, that is benthic recruits and adults. Such data remain, however, scarce due to the difficulty to sample and analyse larvae of minute size. Here, we carried out such an investigation using the marine gastropod Crepidula fornicata. Field sampling of three to four larval pools was conducted over the reproductive season and repeated over 3 years. The genetic composition of larval pools, obtained with 16 microsatellite loci, was compared with that of recruits and adults sampled from the same site and years. In contrast to samples of juveniles and adults, large genetic temporal variations between larval pools produced at different times of the same reproductive season were observed. In addition, full‐ and half‐sibs were detected in early larvae and postdispersal juveniles, pointing to correlated dispersal paths between several pairs of individuals. Inbred larvae were also identified. Such collective larval dispersal was unexpected given the long larval duration of the study species. Our results suggest that each larval pool is produced by a small effective number of reproducers but that, over a reproductive season, the whole larval pool is produced by large numbers of reproducers across space and time.     

A practical guide to methods of parentage analysis

The use of molecular techniques for parentage analysis has been a booming science for over a decade. The most important technological breakthrough was the introduction of microsatellite markers to molecular ecology, an advance that was accompanied by a proliferation and refinement of statistical techniques for the analysis of parentage data. Over the last several years, we have seen steady progress in a number of areas related to parentage analysis, and the prospects for successful studies continue to improve. Here, we provide an updated guide for scientists interested in embarking on parentage analysis in natural or artificial populations of organisms, with a particular focus on computer software packages that implement various methods of analysis. Our survey of the literature shows that there are a few established methods that perform extremely well in the analysis of most types of parentage studies. However, particular experimental designs or study systems can benefit from some of the less well-known computer packages available. Overall, we find that parentage analysis is feasible and satisfying in most systems, and we try to provide a simple roadmap to help other scientists navigate the confusing topography of statistical techniques.     

Isolation and characterization of polymorphic microsatellite loci for parentage analysis in a rhacophorid tree frog (Chirixalus eiffingeri) with unusual parental care

Chirixalus eiffingeri is an arboreal breeding rhacophorid frog with unique parental care behaviours. We developed 11 polymorphic microsatellites as genetic makers for parentage analysis to resolve the ecology of parental care. The numbers of alleles per locus ranged from two to 17. The observed and expected heterozygosities averaged 0.433 and 0.656, respectively. Total exclusionary probability of these loci is 0.984 when no parental genotypes are known, and is 0.999 when one of the parental genotypes is known. The results indicate that the 11 markers should provide sufficient resolution for inferring genetic parentage in C. eiffingeri.     

Supportive breeding boosts natural population abundance with minimal negative impacts on fitness of a wild population of Chinook salmon

While supportive breeding programmes strive to minimize negative genetic impacts to populations, case studies have found evidence for reduced fitness of artificially produced individuals when they reproduce in the wild. Pedigrees of two complete generations were tracked with molecular markers to investigate differences in reproductive success (RS) of wild and hatchery‐reared Chinook salmon spawning in the natural environment to address questions regarding the demographic and genetic impacts of supplementation to a natural population. Results show a demographic boost to the population from supplementation. On average, fish taken into the hatchery produced 4.7 times more adult offspring, and 1.3 times more adult grand‐offspring than naturally reproducing fish. Of the wild and hatchery fish that successfully reproduced, we found no significant differences in RS between any comparisons, but hatchery‐reared males typically had lower RS values than wild males. Mean relative reproductive success (RRS) for hatchery F₁ females and males was 1.11 (P = 0.84) and 0.89 (P = 0.56), respectively. RRS of hatchery‐reared fish (H) that mated in the wild with either hatchery or wild‐origin (W) fish was generally equivalent to W × W matings. Mean RRS of H × W and H × H matings was 1.07 (P = 0.92) and 0.94 (P = 0.95), respectively. We conclude that fish chosen for hatchery rearing did not have a detectable negative impact on the fitness of wild fish by mating with them for a single generation. Results suggest that supplementation following similar management practices (e.g. 100% local, wild‐origin brood stock) can successfully boost population size with minimal impacts on the fitness of salmon in the wild.     

Large‐scale parentage analysis reveals reproductive patterns and heritability of spawn timing in a hatchery population of steelhead (Oncorhynchus mykiss)

Understanding life history traits is an important first step in formulating effective conservation and management strategies. The use of artificial propagation and supplementation as such a strategy can have numerous effects on the supplemented natural populations and minimizing life history divergence is crucial in minimizing these effects. Here, we use single nucleotide polymorphism (SNP) genotypes for large‐scale parentage analysis and pedigree reconstruction in a hatchery population of steelhead, the anadromous form of rainbow trout. Nearly complete sampling of the broodstock for several consecutive years in two hatchery programmes allowed inference about multiple aspects of life history. Reconstruction of cohort age distribution revealed a strong component of fish that spawn at 2 years of age, in contrast to programme goals and distinct from naturally spawning steelhead in the region, which raises a significant conservation concern. The first estimates of variance in family size for steelhead in this region can be used to calculate effective population size and probabilities of inbreeding, and estimation of iteroparity rate indicates that it is reduced by hatchery production. Finally, correlations between family members in the day of spawning revealed for the first time a strongly heritable component to this important life history trait in steelhead and demonstrated the potential for selection to alter life history traits rapidly in response to changes in environmental conditions. Taken together, these results demonstrate the extraordinary promise of SNP‐based pedigree reconstruction for providing biological inference in high‐fecundity organisms that is not easily achievable with traditional physical tags.     

Assessment of microsatellite and SNP markers for parentage assignment in ex situ African Penguin (Spheniscus demersus) populations

Captive management of ex situ populations of endangered species is traditionally based on pedigree information derived from studbook data. However, molecular methods could provide a powerful set of complementary tools to verify studbook records and also contribute to improving the understanding of the genetic status of captive populations. Here, we compare the utility of single nucleotide polymorphisms (SNPs) and microsatellites (MS) and two analytical methods for assigning parentage in ten families of captive African penguins held in South African facilities. We found that SNPs performed better than microsatellites under both analytical frameworks, but a combination of all markers was most informative. A subset of combined SNP (n = 14) and MS loci (n = 10) provided robust assessments of parentage. Captive or supportive breeding programs will play an important role in future African penguin conservation efforts as a source of individuals for reintroduction. Cooperation among these captive facilities is essential to facilitate this process and improve management. This study provided us with a useful set of SNP and MS markers for parentage and relatedness testing among these captive populations. Further assessment of the utility of these markers over multiple (>3) generations and the incorporation of a larger variety of relationships among individuals (e.g., half‐siblings or cousins) is strongly suggested.     

An evaluation of allowing for mismatches as a way to manage genotyping errors in parentage assignment by exclusion

In parentage assignment by exclusion, using multiple and very polymorphic loci, genotyping errors are a major cause of non‐assignment. Using stochastic simulations, we tested the possibility to allow for mismatches at one or more allele as a way to recover assignment power. This was very efficient provided the set of loci used had a high assignment power (> 99%) and the error rate was not too high (below 3–4%). In these cases, most of the theoretical assignment power could be recovered. We also showed the efficiency of the method in a practical experiment with rainbow trout.     

Pronounced reproductive skew in a natural population of green swordtails, Xiphophorus helleri

For many species in nature, a sire's progeny may be distributed among a few or many dams. This poses logistical challenges--typically much greater across males than across females--for assessing means and variances in mating success (number of mates) and reproductive success (number of progeny). Here we overcome these difficulties by exhaustively analyzing a population of green swordtail fish (Xiphophorus helleri) for genetic paternity (and maternity) using a suite of highly polymorphic microsatellite loci. Genetic analyses of 1476 progeny from 69 pregnant females and 158 candidate sires revealed pronounced skews in male reproductive success both within and among broods. These skews were statistically significant, greater than in females, and correlated in males but not in females with mating success. We also compare the standardized variances in swordtail reproductive success to the few such available estimates for other taxa, notably several mammal species with varied mating systems and degrees of sexual dimorphism. The comparison showed that the opportunity for selection on male X. helleri is among the highest yet reported in fishes, and it is intermediate compared to estimates available for mammals. This study is one of a few exhaustive genetic assessments of joint-sex parentage in a natural fish population, and results are relevant to the operation of sexual selection in this sexually dimorphic, high-fecundity species.     

Effects of genotyping errors on parentage exclusion analysis

Genetic markers are widely used to determine the parentage of individuals in studies of mating systems, reproductive success, dispersals, quantitative genetic parameters and in the management of conservation populations. These markers are, however, imperfect for parentage analyses because of the presence of genotyping errors and undetectable alleles, which may cause incompatible genotypes (mismatches) between parents and offspring and thus result in false exclusions of true parentage. Highly polymorphic markers widely used in parentage analyses, such as microsatellites, are especially prone to genotyping errors. In this investigation, I derived the probabilities of excluding a random (related) individual from parentage and the probabilities of Mendelian-inconsistent errors (mismatches) and Mendelian-consistent errors (which do not cause mismatches) in parent-offspring dyads, when a marker having null alleles, allelic dropouts and false alleles is used in a parentage analysis. These probabilities are useful in evaluating the impact of various types of genotyping errors on the information content of a set of markers in and thus the power of a parentage analysis, in determining the threshold number of genetic mismatches that is appropriate for a parentage exclusion analysis and in estimating the rates of genotyping errors and frequencies of null alleles from observed mismatches between known parent-offspring dyads. These applications are demonstrated by numerical examples using both hypothetical and empirical data sets and discussed in the context of practical parentage exclusion analyses.     

Unbiased estimation of relative reproductive success of different groups: evaluation and correction of bias caused by parentage assignment errors

Parentage assignment is widely applied to studies on mating systems, population dynamics and natural selection. However, little is known about the consequence of assignment errors, especially when some parents are not sampled. We investigated the effects of two types of error in parentage assignment, failing to assign a true parent (type A) and assigning an untrue parent (type B), on an estimate of the relative reproductive success (RRS) of two groups of parents. Employing a mathematical approach, we found that (i) when all parents are sampled, minimizing either type A or type B error insures the minimum bias on RRS, and (ii) when a large number of parents is not sampled, type B error substantially biases the estimated RRS towards one. Interestingly, however, (iii) when all parents were sampled and both error rates were moderately high, type A error biased the estimated RRS even more than type B error. We propose new methods to obtain an unbiased estimate of RRS and the number of offspring whose parents are not sampled (zW(z)), by correcting the error effects. Applying them to genotypic data from steelhead trout (Oncorhynchus mykiss), we illustrated how to estimate and control the assignment errors. In the data, we observed up to a 30% assignment error and a strong trade-off between the two types of error, depending on the stringency of the assignment decision criterion. We show that our methods can efficiently estimate an unbiased RRS and zW(z) regardless of assignment method, and how to maximize the statistical power to detect a difference in reproductive success between groups.     

Exclusion probabilities of 22 bovine microsatellite markers in fluorescent multiplexes for semiautomated parentage testing

Six multiplexes developed for semiautomated fluorescence genotyping were evaluated for parentage testing. These multiplexes contained primer pairs for the amplification of 22 microsatellites on 17 bovine autosomes. Exclusion probabilities were determined from genotypes of 1022 Holstein cattle and 311 beef cattle belonging to five breeds. Two cases were considered: case 1, genotypes are known for an alleged parent and an offspring but genotypes of a confirmed parent are unknown; and case 2, genotypes are known for an alleged parent, a confirmed parent and an offspring. If the alleged parent is not the true parent, then the 22 markers will exclude the alleged parent with a probability of >0·9986 for case 1 and with a probability of >0·99999 for case 2. On the basis of these exclusion probabilities, the probability that an alleged parent will be falsely included as a parent is in the range of 1/716 to 1/2845 for case 1 and 1/1·2 million to 1/159753 for case 2. In addition to these results, a rapid and efficient non-organic method for extraction of DNA from semen is described.     

Male parentage in dependent-lineage populations of the harvester ant Pogonomyrmex barbatus

We investigated the extent to which workers reproduce in a dependent-lineage population of the monogynous harvester ant Pogonomyrmex barbatus. Dependent-lineage populations contain two interbreeding, yet genetically distinct mitochondrial lineages, each associated with specific alleles at nuclear loci. Workers develop from matings between lineages, and queens develop from matings within lineages, so queens must mate with males of both lineages to produce daughter queens and workers. Males develop from unfertilized eggs and are haploid. Worker production of males could lead to male-mediated gene flow between the lineages if worker-produced males were reproductively capable. This could result in the loss of the dependent-lineage system, because its persistence depends on the maintenance of allelic differences between the lineages. To investigate the extent of worker reproduction in P. barbatus, we genotyped 19-20 males and workers from seven colonies, at seven microsatellite loci, and 1239 additional males at two microsatellite loci. Our methods were powerful enough to detect worker reproduction if workers produced more than 0.39% of males in the population. We detected no worker-produced males; all males appeared to be produced by queens. Thus, worker reproduction is sufficiently infrequent to have little impact on the dependent-lineage system. These results are consistent with predictions based on inclusive fitness theory because the effective queen mating frequency calculated from worker genotypes was 4.26, which is sufficiently high for workers to police those that attempt to reproduce.     

Microsatellite determination of male reproductive success in a natural population of the territorial ornate dragon lizard, Ctenophorus ornatus

It is now evident that the genetic mating system can be very different to the observed mating system. However, it is less well known what makes particular individuals more (or less) successful than expected from the observed system. In this study the observed territorial structure of a field population of the agamid lizard, Ctenophorus ornatus, was compared with the mating system as evidenced by microsatellite parentage assignment. This study also investigated whether any male traits predicted reproductive success. Sixty-five per cent of clutches were sired at least partially by a male other than the main territory-holding male and 35% of clutches were sired by a male with no overlap of the female's territory. Multiple paternity was moderately frequent at 25% of clutches. Male chest patch size predicted territory size and the number of females in the territory, but did not predict reproductive success. Instead, male head depth and body size were independently related to the number of offspring sired. As male head depth also predicted the number of females in a territory, these males are likely to be gaining increased reproductive success as a consequence of the higher number of females in their territories. Larger body size males, however, did not have a greater number of females in their territory and instead had more extra-territorial copulations. Whether these extra-territorial copulations are due to female choice or success in male competition is unknown.     

A genetic analysis of breeding success in the cooperative meerkat (Suricata suricatta)

Measurement of reproductive skew in social groups is fundamentalto understanding the evolution and maintenance of sociality,as it determines the immediate fitness benefits to helpers ofstaying and helping in a group. However, there is a lack ofstudies in natural populations that provide reliable measuresof reproductive skew and the correlates of reproductive success,particularly in vertebrates. We present results of a study thatuses a combination of field and genetic (microsatellite) dataon a cooperatively breeding mongoose, the meerkat (Suricatasuricatta). We sampled 458 individuals from 16 groups at twosites and analyzed parentage of pups in 110 litters with upto 12 microsatellites. We show that there is strong reproductiveskew in favor of dominants, but that the extent of skew differsbetween the sexes and between different sites. Our data suggestthat the reproductive skew arises from incest avoidance andreproductive suppression of the subordinates by the dominants.