Welcome to Acta Agronomica Sinica,

Acta Agron Sin ›› 2013, Vol. 39 ›› Issue (01): 1-11.doi: 10.3724/SP.J.1006.2013.00001

• CROP GENETICS & BREEDING · GERMPLASM RESOURCES · MOLECULAR GENETICS •     Next Articles

Properties of the Test Statistic under Null Hypothesis and the Calculation of LOD Threshold in Quantitative Trait Loci (QTL) Mapping

SUN Zi-Qi,LI Hui-Hui,ZHANG Lu-Yan,WANG Jian-Kang*   

  1. Institute of Crop Sciences / National Key Facility for Crop Gene Resources and Genetic Improvement, Chinese Academy of Agricultural Sciences / CIMMYT China Office, Beijing 100081, China
  • Received:2012-05-13 Revised:2012-09-05 Online:2013-01-12 Published:2012-11-14
  • Contact: 王建康, E-mail: wangjk@caas.net.cn, jkwang@cgiar.org, Tel: 010-82105846

Abstract:

Selecting an appropriate LOD threshold is of great interest in QTL mapping studies. Many approaches can be considered to calculate the critical value throughout a genome, such as simulation-based method, analytical approximation, and empirical method based on permutation test. Many tests are conducted in QTL mapping, which are not mutually independent because the linkage relationship of adjacent markers on chromosomes. In order to declare a significant QTL at a genome-wide significance level, it is necessary to understand the behavior of test statistic under null hypothesis in QTL mapping and to deal with the dependent multiple-test problem arising in the genome-wide test. Our objectives in this study were (1) to investigate the properties of LRT (likelihood ratio test) statistic of one-point scanning under null hypothesis in QTL mapping, (2) to determine the factors affecting the cumulative distribution of maximum LOD score, and (3) to identify the relationship between the effective number of independent tests and the length of chromosome by simulation method. Results indicated that the LRT test statistic in one-dimensional scanning of additive-dominant QTL and two-dimensional scanning of epistatic QTL followed chi-square distributions, and the degree of freedom (df) was equal to the number of genetic parameters to be estimated. For example, degree of freedom in recombinant inbred lines (RIL) population was equal to 1 in one dimensional or two dimensional scanning. Degree of freedom in F2 populations was equal to 2 in one-dimensional scanning and 4 in two-dimensional scanning. Number of chromosome, population size and phenotyping error variance did not have any effect on the distribution of LRT under null hypothesis, and therefore will not affect the selection of LOD threshold. On the contrary, population type, genome size and marker density had significant impacts. For BC1, RIL and F2 populations, the threshold was the smallest in BC1 population and the highest in F2 population. Higher marker density  and longer chromosome resulted in higher LOD threshold. It was identified that the effective number of independent tests (Meff) was proportional to the length of chromosome in one-dimensional scanning of additive-dominant QTL. In two-dimensional scanning of epistatic QTL, it was identified that Meffwas in a squared relationship to the length of chromosome. With the help of Bonferroni correction, we could acquire the relationship between point-wise and genome-wide significance levels. Therefore, it is convenient to calculate the threshold LOD in QTL mapping, given the genome-wide significance level, the population type, marker density and genome size.

Key words: QTL mapping, Likelihood ratio test, LOD score, Null hypothesis, Significance level, Number of independent tests

[1]Lander E S, Botstein D. Mapping mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics, 1989, 121: 185–199



[2]van Ooijen J W. LOD significance thresholds for QTL analysis in experimental populations of diploid species. Heredity, 1999, 83: 613–624



[3]Dupuis J, Siegmund D. Statistical methods for mapping quantitative trait loci from a dense set of markers. Genetics, 1999, 151: 373–386



[4]Feingold E, Brown P O, Siegmund D. Gaussian models for genetic linkage analysis using complete high-resolution maps of identity by descent. Am J Hum Genet, 1993, 53: 234–251



[5]Rebaï A, Goffinet B, Mangin B. Approximate thresh olds of interval mapping tests for QTL detection. Genetics, 1994, 138: 235–240



[6]Dupuis J. Statistical Problems Associated with Mapping Complex and Quantitative Traits from Genomic Mismatch Scanning Data. PhD Dissertation of Stanford University, 1994



[7]Kruglyak L, Lander E S. A nonparametric approach for mapping quantitative trait loci. Genetics, 1995, 139: 1421–1428



[8]Lander E, Kruglyak L. Genetic dissection of complex traits: guidelines for interpreting and reporting linkage results. Nature Genet, 1995, 11: 241–247



[9]Piepho H P. A quick method for computing approximate thresholds for quantitative trait loci detection. Genetics, 2001, 157: 425–432



[10]Zou F, Fine J P, Hu J H, Lin D Y. An efficient resampling method for assessing genome-wide statistical significance in mapping quantitative trait loci. Genetics, 2004, 168: 2307–2316



[11]Chang M N, Wu R, Wu S S, Casella G. Score statistics for mapping quantitative trait loci. Stat Appl Genet Mol Biol, 2009, 8(1): Article 16



[12]Churchill G A, Doerge R W. Empirical threshold values for quantitative trait mapping. Genetics, 1994, 138: 963–971



[13]Cheverud J M. A simple correction for multiple comparisons in interval mapping genome scans. Heredity, 2001, 87: 52–58



[14]Nyholt D R. A simple correction for multiple testing for single-nucleotide polymorphisms in linkage disequilibrium with each other. Am J Hum Genet, 2004, 74: 765–769



[15]Li J, Ji L. Adjusting multiple testing in multilocus analyses using the eigenvalues of a correlation matrix. Heredity, 2005, 95: 221–227



[16]Gao X, Starmer J, Martin E R. A multiple testing correction method for genetic association studies using correlated single nucleotide polymorphisms. Genetic Epidemiology, 2008, 32: 361–369



[17]Li H, Ribaut J M, Li Z, Wang J. Inclusive composite interval mapping (ICIM) for digenic epistasis of quantitative traits in biparental populations. Theor Appl Genet, 2008, 116: 243–260



[18]Zhang L, Li H, Li Z, Wang J. Interactions between markers can be caused by the dominance effect of QTL. Genetics, 2008, 180: 1177–1190



[19]Zhang L, Li H, Wang J. Statistical power of inclusive composite interval mapping in detecting digenic epistasis showing common F2 segregation ratios. J Integr Plant Biol, 2012, 54: 270–279



[20]Doerge R W, Rebaï A. Significance thresholds for QTL interval mapping tests. Heredity, 1996, 76: 459–464



[21]Wang J, Chapman S C, Bonnett D G, Rebetzke G J. Simultaneous selection of major and minor genes: use of QTL to increase selection efficiency of coleoptiles length of wheat (Triticum aestivum L.). Theor Appl Genet, 2009, 119: 65–74



[22]Ribaut J M, Hoisington D A, Deutsch J A, Jiang C, González-de-León D. Identification of quantitative trait loci under drought conditions in tropical maize: 1. Flowering parameters and the anthesis-silking interval. Theor Appl Genet, 1996, 92: 905–914



[23]Zhang L, Wang S, Li H, Deng Q, Zheng A, Li S, Li P, Li Z, Wang J. Effects of missing marker and segregation distortion on QTL mapping in F2 populations. Theor Appl Genet, 2010, 121: 1071–1082



[24]Li H, Ye G, Wang J. A modified algorithm for the improvement of composite interval mapping. Genetics, 2007, 175: 361–374



[25]Li H-H(李慧慧), Zhang L-Y(张鲁燕), Wang J-K(王建康). Analytical answers to frequently asked questions in quantitative trait locus mapping. Acta Agron Sin (作物学报), 2010, 36(6): 918–931 (in Chinese with English abstract)
[1] Zheng Yu-Zhen, Qi Fei-Yan, Sun Zi-Qi, Liu Hua, Qin Li, Shi Lei, Wang Juan, Wang Meng-Meng, Han Suo-Yi, Xu Jing, Miao Li-Juan, Huang Bing-Yan, Dong Wen-Zhao, Zheng Zheng, Zhang Xin-You. QTL mapping of total very long-chain fatty acids and seven fatty acid components in peanut seeds [J]. Acta Agronomica Sinica, 2026, 52(6): 1646-1657.
[2] Liu Chang-You, Wang Shen, Shi Hui-Ying, Shen Ying-Chao, Sun Lei, Wang Yan, Zhang Zhi-Xiao, Su Qiu-Zhu, Tian Jing, Fan Bao-Jie. QTL mapping for bruchid resistance in an adzuki bean distant hybridization population using rice bean genetic resources [J]. Acta Agronomica Sinica, 2026, 52(3): 936-944.
[3] YANG Hai-Yang, WU Lin-Xuan, LI Bo-Wen, SHI Han-Feng, YUAN Xi-Long, LIU Jin-Zhao, CAI Hai-Rong, CHEN Shi-Yi, GUO Tao, WANG Hui. OsWRI3, identified based on QTL mapping, regulates seed shattering in rice [J]. Acta Agronomica Sinica, 2025, 51(7): 1712-1724.
[4] HU Meng, SHA Dan, ZHANG Sheng-Rui, GU Yong-Zhe, ZHANG Shi-Bi, LI Jing, SUN Jun-Ming, QIU Li-Juan, LI Bin. QTL mapping and candidate gene screening for branch number in soybean [J]. Acta Agronomica Sinica, 2025, 51(7): 1747-1756.
[5] SHAO Shun-Wei, CHEN Zhuo, LAN Zhen-Dong, CAI Xing-Kui, ZOU Hua-Fen, LI Chen-Xi, TANG Jing-Hua, ZHU Xi, ZHANG Yu, DONG Jian-Ke, JIN Hui, SONG Bo-Tao. QTL mapping of tuber eye depth based on BSA-seq technique [J]. Acta Agronomica Sinica, 2025, 51(7): 1725-1735.
[6] ZHANG Jin-Ze, ZHOU Qing-Guo, XIAO Li-Jing, JIN Hai-Run, OU-YANG Qing-Jing, LONG Xu, YAN Zhong-Bin, TIAN En-Tang. QTL mapping and candidate gene analysis of glucosinolate content in various tissues of Brassica juncea [J]. Acta Agronomica Sinica, 2025, 51(5): 1166-1177.
[7] YONG Rui, HU Wen-Jing, WU Di, WANG Zun-Jie, LI Dong-Sheng, ZHAO Die, YOU Jun-Chao, XIAO Yong-Gui, WANG Chun-Ping. Identification and validation of quantitative trait loci for grain number per spike showing pleiotropic effect on thousand grain weight in bread wheat (Triticum aestivum L.) [J]. Acta Agronomica Sinica, 2025, 51(2): 312-323.
[8] GUO Shu-Hui, PAN Zhuan-Xia, ZHAO Zhan-Sheng, YANG Liu-Liu, HUANG-FU Zhang-Long, GUO Bao-Sheng, HU Xiao-Li, LU Ya-Dan, DING Xiao, WU Cui-Cui, LAN Gang, LYU Bei-Bei, TAN Feng-Ping, LI Peng-Bo. Genetic analysis of a major fiber length locus on chromosome D11 of upland cotton [J]. Acta Agronomica Sinica, 2025, 51(2): 383-394.
[9] WANG Zhe, HU Yan-Ling, GONG Fang-Yi, YI Rui, ZHAO Shu-Hong, LIU Rui-Qin, LIU Yu-Hang, ZHANG Tian, ZHANG Ya-Zhou, ZHENG You-Liang, LIU Deng-Cai, HUANG Lin, WU Bi-Hua. QTL mapping of grain protein content in the introgression line BAd7-209 derived from wild emmer [J]. Acta Agronomica Sinica, 2025, 51(12): 3238-3250.
[10] XU Xiao-Wei, FENG Jing, WANG Feng-Tao, TONG Zhao-Yang, ZHANG Jian-Zhou, LI Chun-Ying, LIN Rui-Ming. QTL mapping of adult plant resistance to stripe rust in the Chinese wheat landrace Canlaomai [J]. Acta Agronomica Sinica, 2025, 51(11): 2933-2943.
[11] YANG Jing-Fa, YU Xin-Lian, YAO You-Hua, YAO Xiao-Hua, WANG Lei, WU Kun-Lun, LI Xin. QTL mapping of tiller angle in qingke (Hordeum vulgare L.) [J]. Acta Agronomica Sinica, 2025, 51(1): 260-272.
[12] HAN Li, TANG Sheng-Sheng, LI Jia, HU Hai-Bin, LIU Long-Long, WU Bin. Construction of SNP high-density genetic map and localization of QTL for β-glucan content in oats [J]. Acta Agronomica Sinica, 2024, 50(7): 1710-1718.
[13] BI Jun-Ge, ZENG Zhan-Kui, LI Qiong, HONG Zhuang-Zhuang, YAN Qun-Xiang, ZHAO Yue, WANG Chun-Ping. QTL mapping and KASP marker development of grain quality-relating traits in two wheat RIL populations [J]. Acta Agronomica Sinica, 2024, 50(7): 1669-1683.
[14] QIN Na, YE Zhen-Yan, ZHU Can-Can, FU Sen-Jie, DAI Shu-Tao, SONG Ying-Hui, JING Ya, WANG Chun-Yi, LI Jun-Xia. QTL mapping for flavonoid content and seed color in foxtail millet [J]. Acta Agronomica Sinica, 2024, 50(7): 1719-1727.
[15] ZHENG Xue-Qing, WANG Xing-Rong, ZHANG Yan-Jun, GONG Dian-Ming, QIU Fa-Zhan. Mapping of QTL for ear-related traits and prediction of key candidate genes in maize [J]. Acta Agronomica Sinica, 2024, 50(6): 1435-1450.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!