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.
['Anderson|Eric C|EC|', 'Ng|Thomas C|TC|']
D000465:Algorithms D000818:Animals D001499:Bayes Theorem D019295:Computational Biology D003198:Computer Simulation D005819:Genetic Markers D005838:Genotype D016013:Likelihood Functions D008390:Markov Chains D010375:Pedigree D020641:Polymorphism, Single Nucleotide D012474:Salmon D035781:Siblings
Sum-Product algorithm Full-sibling reconstruction Multigeneration pedigree inference Relationship inference