Bayesian mapping of genome-wide epistatic imprinted loci for quantitative traits

Genomic imprinting, an epigenetic phenomenon of parent-of-origin-specific gene expression, has been widely observed in plants, animals, and humans. To detect imprinting genes influencing quantitative traits, the least squares and maximum likelihood approaches for fitting a single quantitative trait locus (QTL) and Bayesian methods for simultaneously modeling multiple QTL have been adopted, respectively, in various studies. However, most of these studies have only estimated imprinting main effects and thus ignored imprinting epistatic effects. In the presence of extremely complex genomic imprinting architectures, we introduce a Bayesian model selection method to analyze the multiple interacting imprinted QTL (iQTL) model. This approach will greatly enhance the computational efficiency through setting the upper bound of the number of QTLs and performing selective sampling for QTL parameters. The imprinting types of detected main-effect QTLs can be estimated from the Bayes factor statistic formulated by the posterior probabilities for the genetic effects being compared. The performance of the proposed method is demonstrated by several simulation experiments. Moreover, this method is applied to dissect the imprinting genetic architecture for body weight in mouse and fruit weight in tomato. Matlab code for implementing this approach will be available from the authors upon request.