Pedigree Inference Research Articles

Improved computations for relationship inference using low-coverage sequencing data

Pedigree inference, for example determining whether two persons are second cousins or unrelated, can be done by comparing their genotypes at a selection of genetic markers. When the data for one or more of the persons is from low-coverage next generation sequencing (lcNGS), currently available computational methods either ignore genetic linkage or do not take advantage of the probabilistic nature of lcNGS data, relying instead on first estimating the genotype. We provide a method and software (see familias.name/lcNGS) bridging the above gap. Simulations indicate how our results are considerably more accurate compared to some previously available alternatives. Our method, utilizing a version of the Lander-Green algorithm, uses a group of symmetries to speed up calculations. This group may be of further interest in other calculations involving linked loci.

Parent-offspring inference in inbred populations.

Genealogical relationships are fundamental components of genetic studies. However, it is often challenging to infer correct and complete pedigrees even when genome-wide information is available. For example, inbreeding can obscure genetic differences between individuals, making it difficult to even distinguish first-degree relatives such as parent-offspring from full siblings. Similarly, genotyping errors can interfere with the detection of genetic similarity between parents and their offspring. Inbreeding is common in natural, domesticated, and experimental populations and genotyping of these populations often has more errors than in human data sets, so efficient methods for building pedigrees under these conditions are necessary. Here, we present a new method for parent-offspring inference in inbred pedigrees called specific parent-offspring relationship estimation (spore). spore is vastly superior to existing pedigree-inference methods at detecting parent-offspring relationships, in particular when inbreeding is high or in the presence of genotyping errors, or both. spore therefore fills an important void in the arsenal of pedigree inference tools.

Bonsai: An efficient method for inferring large human pedigrees from genotype data

SummaryPedigree inference from genotype data is a challenging problem, particularly when pedigrees are sparsely sampled and individuals may be distantly related to their closest genotyped relatives. We present a method that infers small pedigrees of close relatives and then assembles them into larger pedigrees. To assemble large pedigrees, we introduce several formulas and tools including a likelihood for the degree separating two small pedigrees, a generalization of the fast DRUID point estimate of the degree separating two pedigrees, a method for detecting individuals who share background identity-by-descent (IBD) that does not reflect recent common ancestry, and a method for identifying the ancestral branches through which distant relatives are connected. Our method also takes several approaches that help to improve the accuracy and efficiency of pedigree inference. In particular, we incorporate age information directly into the likelihood rather than using ages only for consistency checks and we employ a heuristic branch-and-bound-like approach to more efficiently explore the space of possible pedigrees. Together, these approaches make it possible to construct large pedigrees that are challenging or intractable for current inference methods.

Isolation and characterization of 15 new microsatellite markers for the globally endangered Lear's macaw Anodorhynchus leari.

A set of 16 microsatellite markers was characterized for Lear's macaw (Anodorhynchus leari) using DNA samples from captive individuals. Extending this molecular toolkit, including the use of samples from wild individuals, is expected to provide the required power of resolution for pedigree inference of both wild and captive individuals, and could support research on the genetic structure of wild populations.We characterize a set of 15 microsatellite markers optimized for the Lear's macaw, developed from a microsatellite-enriched library in a three-step procedure. Primer pairs were initially designed for 62 microsatellites with > 7 tandem repetitions. After amplification of DNA of five wild individuals from different localities, 22 loci seemed to be polymorphic and were further tested on 12 wild nestling samples. Fifteen unlinked loci showed unambiguous peaks and low to moderate polymorphism levels. The combination of the four most polymorphic markers allowed individual identification even of putative sibs.These markers complement previously described microsatellites developed for A. leari and constitute a fundamental genetic toolkit for the investigation of the genetics of both wild and captive populations, thus assisting integrated management plans for the conservation of this globally endangered species.

Composite likelihood method for inferring local pedigrees.

Pedigrees contain information about the genealogical relationships among individuals and are of fundamental importance in many areas of genetic studies. However, pedigrees are often unknown and must be inferred from genetic data. Despite the importance of pedigree inference, existing methods are limited to inferring only close relationships or analyzing a small number of individuals or loci. We present a simulated annealing method for estimating pedigrees in large samples of otherwise seemingly unrelated individuals using genome-wide SNP data. The method supports complex pedigree structures such as polygamous families, multi-generational families, and pedigrees in which many of the member individuals are missing. Computational speed is greatly enhanced by the use of a composite likelihood function which approximates the full likelihood. We validate our method on simulated data and show that it can infer distant relatives more accurately than existing methods. Furthermore, we illustrate the utility of the method on a sample of Greenlandic Inuit.

Impacts of ecology and behavior on Antarctic fur seal remating and relatedness

Antarctic fur seals (Arctocephalus gazella) are polygynous and both sexes are typically faithful to a breeding site. These characteristics could promote remating among individuals over time, leading to increased relatedness levels and negatively affecting genetic diversity. To examine this issue, the reproductive output of 55 females was monitored annually for 12years and their pups were sampled (n=280) and genotyped using 17 microsatellite markers. A full likelihood pedigree inference method was used to confirm maternities inferred in the field and estimate the number of full sibling pups born across years. Relatedness coefficients were estimated for pairs of individuals in the pedigree and compared to simulated values for each relationship category. There were nine cases where a female mated with the same male twice and one case where a female mated with the same male three times over the study period. The observed relatedness coefficients estimated among the sampled pups matched the simulated distribution for half-siblings. In addition, no first order relatives were found among the fur seal mothers studied, nor did observed relatedness coefficient distributions differ significantly from simulated values. Together, these results suggest a low remating rate and a negligible effect of remating on pair-wise relatedness. Territorial male replacement over time as well as female small-scale movements, driven by suitable pupping habitat, likely contribute to the low remating frequency observed in Antarctic fur seals.

Bayesian pedigree inference with small numbers of single nucleotide polymorphisms via a factor-graph representation

We develop a computational framework for addressing pedigree inference problems using small numbers (80–400) of single nucleotide polymorphisms (SNPs). Our approach relaxes the assumptions, which are commonly made, that sampling is complete with respect to the pedigree and that there is no genotyping error. It relies on representing the inferred pedigree as a factor graph and invoking the Sum-Product algorithm to compute and store quantities that allow the joint probability of the data to be rapidly computed under a large class of rearrangements of the pedigree structure. This allows efficient MCMC sampling over the space of pedigrees, and, hence, Bayesian inference of pedigree structure. In this paper we restrict ourselves to inference of pedigrees without loops using SNPs assumed to be unlinked. We present the methodology in general for multigenerational inference, and we illustrate the method by applying it to the inference of full sibling groups in a large sample (n=1157) of Chinook salmon typed at 95 SNPs. The results show that our method provides a better point estimate and estimate of uncertainty than the currently best-available maximum-likelihood sibling reconstruction method. Extensions of this work to more complex scenarios are briefly discussed.

Hsphase: an R package for pedigree reconstruction, detection of recombination events, phasing and imputation of half-sib family groups.

BackgroundIdentification of recombination events and which chromosomal segments contributed to an individual is useful for a number of applications in genomic analyses including haplotyping, imputation, signatures of selection, and improved estimates of relationship and probability of identity by descent. Genotypic data on half-sib family groups are widely available in livestock genomics. This structure makes it possible to identify recombination events accurately even with only a few individuals and it lends itself well to a range of applications such as parentage assignment and pedigree verification.ResultsHere we present hsphase, an R package that exploits the genetic structure found in half-sib livestock data to identify and count recombination events, impute and phase un-genotyped sires and phase its offspring. The package also allows reconstruction of family groups (pedigree inference), identification of pedigree errors and parentage assignment. Additional functions in the package allow identification of genomic mapping errors, imputation of paternal high density genotypes from low density genotypes, evaluation of phasing results either from hsphase or from other phasing programs. Various diagnostic plotting functions permit rapid visual inspection of results and evaluation of datasets.ConclusionThe hsphase package provides a suite of functions for analysis and visualization of genomic structures in half-sib family groups implemented in the widely used R programming environment. Low level functions were implemented in C++ and parallelized to improve performance. hsphase was primarily designed for use with high density SNP array data but it is fast enough to run directly on sequence data once they become more widely available. The package is available (GPL 3) from the Comprehensive R Archive Network (CRAN) or from http://www-personal.une.edu.au/~cgondro2/hsphase.htm.

High male reproductive success in a low-density Antarctic fur seal (Arctocephalus gazella) breeding colony

Understanding how population density influences mating systems may lead to important insights into the plasticity of breeding behavior, but few natural systems allow for such studies. Antarctic fur seals (Arctocephalus gazella) provide an interesting model system because they breed in colonies of varying densities. Previous studies have largely focused on a high-density site at Bird Island, South Georgia. Here, 13 highly polymorphic microsatellite loci were used to conduct a genetic analysis of a low-density breeding colony of this species at Livingston Island, approximately 1,600 km south of South Georgia. The majority of adults seen ashore (n = 54) were sampled together with every pup born (n = 97) over four consecutive years. Paternities were confidently assigned for 34 out of the 97 pups. Two out of 23 sampled males accounted for the paternity of 28 % of all pups sampled during the study and 82 % of the pups with an assigned father. Moreover, a full likelihood pedigree inference method assigned a further eight paternities to a single unsampled male seal that is inferred to have held a territory during the season before the study began. The most successful males in our study easily surpassed the previous record for the total number of pups sired per male seal for the species. Furthermore, we identified two triads of full siblings implying that their parents remated in three consecutive years. These findings suggest that territorial male fur seals may achieve greater success in monopolizing access to breeding females when population density is relatively low.

On the choice of prior density for the Bayesian analysis of pedigree structure

This article is concerned with the choice of structural prior density for use in a fully Bayesian approach to pedigree inference. It is found that the choice of prior has considerable influence on the accuracy of the estimation. To guide this choice, a scale invariance property is introduced. Under a structural prior with this property, the marginal prior distribution of the local properties of a pedigree node (number of parents, offspring, etc.) does not depend on the number of nodes in the pedigree. Such priors are found to arise naturally by an application of the Minimum Description Length (MDL) principle, under which construction of a prior becomes equivalent to the problem of determining the length of a code required to encode a pedigree, using the principles of information theory. The approach is demonstrated using simulated and actual data, and is compared to two well-known applications, CERVUS and COLONY.

Molecular marker‐based pedigrees for animal conservation biologists

AbstractPedigrees, depicting the genealogical relationships between individuals in a population, are of fundamental importance to several research areas including conservation biology. For example, they are useful for estimating inbreeding, heritability, selection, studying kin selection and for measuring gene flow between populations. Pedigrees constructed from direct observations of reproduction are usually unavailable for wild populations. Therefore, pedigrees for these populations are usually estimated using molecular marker data. Despite their obvious importance, and the fact that pedigrees are conceptually well understood, the methods, and limitations of marker‐based pedigree inference are often less well understood. Here we introduce animal conservation biologists to molecular marker‐based pedigrees. We briefly describe the history of pedigree inference research, before explaining the underlying theory and basic mechanics of pedigree construction using standard methods. We explain the assumptions and limitations that accompany many of these methods, before going on to explain methods that relax several of these assumptions. Finally, we look to future and discuss some recent exciting advances such as the use of single‐nucleotide polymorphisms, inference of multigenerational pedigrees and incorporation of non‐genetic data such as field observations into the calculations. We also provide some guidelines on efficient marker selection in order to maximize accuracy and power. Throughout we use examples from the field of animal conservation and refer readers to appropriate software where possible. It is our hope that this review will help animal conservation biologists to understand, choose, and use the methods and tools of this fast‐moving field.

Founder contribution and pedigree inference in a captive breeding colony of lion‐tailed macaques, using mitochondrial DNA and DNA fingerprint analyses

AbstractCaptive colony genetic management is dependent on knowledge of the colony pedigree. In most zoos, records for colonies held prior to the past decade are incomplete or nonexistent. We have used a nested set of molecular techniques to infer genetic relationships in a colony of 15 lion‐tailed macaques after twenty years of unmonitored reproduction. The troop was first divided into lines of matrilineal descent using mitochondrial DNA RFLP patterns. Paternity and maternity were then assigned on the basis of DNA fingerprint analysis. As mtDNA is maternally inherited, it is of no use in paternity assessment, but enabled subdivision of the colony into maternal lines. DNA fingerprinting would have been of little use by itself because of high levels of relatedness in several generations of offspring within the colony. The combination of the two techniques, however, enabled strong inference of both paternity and maternity in all cases. These data were instrumental in the inclusion of this troop in the AAZPA Species Survival Plan.