Bayesian association mapping for quantitative traits in a mixture of two populations

We introduce a novel Bayesian approach to estimate and account for population structure simultaneously with association mapping of multiple quantitative trait loci. The method is designed for an analysis of unrelated individuals from a mixture of two populations (no admixture), where the individual population memberships are unknown. In our approach, the population structure is estimated and accounted for by using data on additional "grouping" markers which are assumed to be in Hardy-Weinberg equilibrium within the populations but have different allele frequencies between the populations. We use Bayesian hierarchical modeling and Markov chain Monte Carlo estimation, where we allow both population stratification and genetic heterogeneity. In our model the number of quantitative trait loci and their positions are treated as random variables, and we obtain their posterior distributions. Here we select the candidate and the grouping markers based on results from a preliminary SOLAR analysis.