We propose a birth-death-merge data-driven reversible jump (DDRJ) for multiple-QTL mapping where the phenotypic trait is modeled as a linear function of the additive and dominance effects of the unknown QTL genotypes. We compare the performance of the proposed methodology, usual reversible jump (RJ) and multiple-interval mapping (MIM), using simulated and real data sets. Compared with RJ, DDRJ shows a better performance to estimate the number of QTLs and their locations on the genome mainly when the QTLs effect is moderate, basically as a result of better mixing for transdimensional moves. The inclusion of a merge step of consecutive QTLs in DDRJ is efficient, under tested conditions, to avoid the split of true QTL's effects between false QTLs and, consequently, selection of the wrong model. DDRJ is also more precise to estimate the QTLs location than MIM in which the number of QTLs need to be specified in advance. As DDRJ is more efficient to identify and characterize QTLs with smaller effect, this method also appears to be useful and brings contributions to identifying single-nucleotide polymorphisms (SNPs) that usually have a small effect on phenotype.
Keywords: QTL mapping; birth–death–merge movements; data-driven reversible jump; mixing of MCMC; model selection.
Copyright © 2016 by the Genetics Society of America.