Free Access
Genet. Sel. Evol.
Volume 35, Number 6, November-December 2003
Page(s) 585 - 604
Genet. Sel. Evol. 35 (2003) 585-604
DOI: 10.1051/gse:2003041

A comparison of alternative methods to compute conditional genotype probabilities for genetic evaluation with finite locus models

Liviu R. Totira, Rohan L. Fernandoa, b, Jack C.M. Dekkersa, b, Soledad A. Fernándezc and Bernt Guldbrandtsend

a  Department of Animal Science, Iowa State University, Ames, IA 50011-3150, USA
b  Lawrence H. Baker Center for Bio-informatics and Biological Statistics, Iowa State University, Ames, IA 50011-3150, USA
c  Department of Statistics, The Ohio State University, Columbus, OH 43210, USA
d  Danish Institute of Animal Science, Foulum, Denmark

(Received 27 February 2002; accepted 5 May 2003)

An increased availability of genotypes at marker loci has prompted the development of models that include the effect of individual genes. Selection based on these models is known as marker-assisted selection (MAS). MAS is known to be efficient especially for traits that have low heritability and non-additive gene action. BLUP methodology under non-additive gene action is not feasible for large inbred or crossbred pedigrees. It is easy to incorporate non-additive gene action in a finite locus model. Under such a model, the unobservable genotypic values can be predicted using the conditional mean of the genotypic values given the data. To compute this conditional mean, conditional genotype probabilities must be computed. In this study these probabilities were computed using iterative peeling, and three Markov chain Monte Carlo (MCMC) methods - scalar Gibbs, blocking Gibbs, and a sampler that combines the Elston Stewart algorithm with iterative peeling (ESIP). The performance of these four methods was assessed using simulated data. For pedigrees with loops, iterative peeling fails to provide accurate genotype probability estimates for some pedigree members. Also, computing time is exponentially related to the number of loci in the model. For MCMC methods, a linear relationship can be maintained by sampling genotypes one locus at a time. Out of the three MCMC methods considered, ESIP, performed the best while scalar Gibbs performed the worst.

Key words: genotype probabilities / finite locus models / Markov chain Monte Carlo

Correspondence and reprints: Liviu R. Totir

© INRA, EDP Sciences 2003