Abstract
Most theoretical studies on epistatic QTL mapping have shown that this procedure is powerful, efficient to control the false positive rate (FPR), and precise to localize QTLs. The objective of this simulation-based study was to show that mapping epistatic QTLs is not an almost-perfect process. We simulated 50 samples of 400 F(2) plants/recombinant inbred lines, genotyped for 975 SNPs distributed in 10 chromosomes of 100 cM. The plants were phenotyped for grain yield, assuming 10 epistatic QTLs and 90 minor genes. Adopting basic procedures of r/qtl package, we maximized the power of detection for QTLs (56-74%, on average) but associated with a very high FPR (65%) and a low detection power for the epistatic pairs (7%). Increasing the average detection power for epistatic pairs (14%) highly increased the related FPR. Adopting a procedure to find the best balance between power and FPR, there was a significant decrease in the power of QTL detection (17-31%, on average), associated with a low average detection power for epistatic pairs (8%) and an average FPR of 31% for QTLs and 16% for epistatic pairs. The main reasons for these negative results are a simplified specification of the coefficients of epistatic effects, as theoretically proved, and the effects of minor genes since 2/3 of the FPR for QTLs were due to them. We hope that this study, including the partial derivation of the coefficients of epistatic effects, motivates investigations on how to increase the power of detection for epistatic pairs, effectively controlling the FPR.