Bayesian Quantitative Trait Locus Mapping Using Inferred Haplotypes

Abstract

We describe a fast hierarchical Bayesian method for mapping quantitative trait loci by haplotype-based association, applicable when haplotypes are not observed directly but are inferred from multiple marker genotypes. The method avoids the use of a Monte Carlo Markov chain by employing priors for which the likelihood factorizes completely. It is parameterized by a single hyperparameter, the fraction of variance explained by the quantitative trait locus, compared to the frequentist fixed-effects model, which requires a parameter for the phenotypic effect of each combination of haplotypes; nevertheless it still provides estimates of haplotype effects. We use simulation to show that the method matches the power of the frequentist regression model and, when the haplotypes are inferred, exceeds it for small QTL effect sizes. The Bayesian estimates of the haplotype effects are more accurate than the frequentist estimates, for both known and inferred haplotypes, which indicates that this advantage is independent of the effect of uncertainty in haplotype inference and will hold in comparison with frequentist methods in general. We apply the method to data from a panel of recombinant inbred lines of Arabidopsis thaliana, descended from 19 inbred founders.

High-Resolution Genetic Mapping Using the Mouse Diversity Outbred Population

“…Larger sample sizes are clearly needed to support a regression model with that many degrees of freedom. Alternative implementations that retain some degree of parsimony and, at the same time, allow for nonadditive intralocus effects are the subject of current research (Durrant and Mott 2010;Lenarcic et al 2012). The additive genetic model provides a robust compromise for detecting genetic loci that have a marginal effect.…”

Measurement of serum PON-3 concentration: method evaluation, reference values, and influence of genotypes in a population-based study

“…The anti-PON3 antibody was generated by inoculating rabbits with the peptide CRVNASQEVEPVEPEN, which is specifi c to ma- with the EM algorithm ( 29 ) and the Markov Chain Monte Carlo method ( 30 ). Results are shown as means and SD in parentheses (parametric) or as medians and 95% CI (non-parametric).…”

Fine-mapping diabetes-related traits, including insulin resistance, in heterogeneous stock rats

“…One approach to obtaining stable estimates of allele effects is to constrain the values can take by penalization or formal hierarchical modeling. Until now, this has been achieved only for simple additive effect models that exclude covariate or random effects: specifically in inbred plants lines or the mouse Collaborative Cross using multiple imputation (17,39), or in the mouse Diversity Outbred population (akin to a rodent HS) using ridge regression (67).…”