{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,2,27]],"date-time":"2026-02-27T04:25:29Z","timestamp":1772166329601,"version":"3.50.1"},"reference-count":16,"publisher":"Springer Science and Business Media LLC","issue":"1","license":[{"start":{"date-parts":[[2021,3,1]],"date-time":"2021-03-01T00:00:00Z","timestamp":1614556800000},"content-version":"tdm","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by\/4.0\/"},{"start":{"date-parts":[[2021,3,1]],"date-time":"2021-03-01T00:00:00Z","timestamp":1614556800000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/501100004085","name":"Ministry of Education, Science and Technology","doi-asserted-by":"publisher","award":["2019R1F1A1062005)"],"award-info":[{"award-number":["2019R1F1A1062005)"]}],"id":[{"id":"10.13039\/501100004085","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100004085","name":"Ministry of Education, Science and Technology","doi-asserted-by":"publisher","award":["2017M3A9C4065964"],"award-info":[{"award-number":["2017M3A9C4065964"]}],"id":[{"id":"10.13039\/501100004085","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":["link.springer.com"],"crossmark-restriction":false},"short-container-title":["BioData Mining"],"published-print":{"date-parts":[[2021,12]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:sec>\n                    <jats:title>Background<\/jats:title>\n                    <jats:p>For gene-gene interaction analysis, the multifactor dimensionality reduction (MDR) method has been widely employed to reduce multi-levels of gene-gene interactions into high- or low-risk groups using a binary attribute. For the survival phenotype, the Cox-MDR method has been proposed using a martingale residual of a Cox model since Surv-MDR was first proposed using a log-rank test statistic. Recently, the KM-MDR method was proposed using the Kaplan-Meier median survival time as a classifier. All three methods used the cross-validation procedure to identify single nucleotide polymorphism (SNP) using SNP interactions among all possible SNP pairs. Furthermore, these methods require the permutation test to verify the significance of the selected SNP pairs. However, the unified model-based multifactor dimensionality reduction method (UM-MDR) overcomes this shortcoming of MDR by unifying the significance testing with the MDR algorithm within the framework of the regression model. Neither cross-validation nor permutation testing is required to identify SNP by SNP interactions in the UM-MDR method. The UM-MDR method comprises two steps: in the first step, multi-level genotypes are classified into high- or low-risk groups, and an indicator variable for the high-risk group is defined. In the second step, the significance of the indicator variable of the high-risk group is tested in the regression model included with other adjusting covariates. The Cox-UMMDR method was recently proposed by combining Cox-MDR with UM-MDR to identify gene-gene interactions associated with the survival phenotype. In this study, we propose two simple methods either by combining KM-MDR with UM-MDR, called KM-UMMDR or by modifying Cox-UMMDR by adjusting for the covariate effect in step 1, rather than in step 2, a process called Cox2-UMMDR. The KM-UMMDR method allows the covariate effect to be adjusted for in the regression model of step 2, although KM-MDR cannot adjust for the covariate effect in the classification procedure of step 1. In contrast, Cox2-UMMDR differs from Cox-UMMDR in the sense that the martingale residuals are obtained from a Cox model by adjusting for the covariate effect in step 1 of Cox2-UMMDR whereas Cox-UMMDR adjusts for the covariate effect in the regression model in step 2. We performed simulation studies to compare the power of several methods such as KM-UMMDR, Cox-UMMDR, Cox2-UMMDR, Cox-MDR, and KM-MDR by considering the effect of covariates and the marginal effect of SNPs. We also analyzed a real example of Korean leukemia patient data for illustration and a short discussion is provided.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Results<\/jats:title>\n                    <jats:p>In the simulation study, two different scenarios are considered: the first scenario compares the power of the cases with and without the covariate effect. The second scenario is to compare the power of cases with the main effect of SNPs versus without the main effect of SNPs. From the simulation results, Cox-UMMDR performs the best across all scenarios among KM-UMMDR, Cox2-UMMDR, Cox-MDR and KM-MDR. As expected, both Cox-UMMDR and Cox-MDR perform better than KM-UMMDR and KM-MDR when a covariate effect exists because the former adjusts for the covariate effect but the latter cannot. However, Cox2-UMMDR behaves similarly to KM-UMMDR and KM-MDR even though there is a covariate effect. This implies that the covariate effect would be more efficiently adjusted for in the regression model of the second step rather than under the classification procedure of the first step. When there is a main effect of any SNP, Cox-UMMDR, Cox2-UMMDR and KM-UMMDR perform better than Cox-MDR and KM-MDR if the main effects of SNPs are properly adjusted for in the regression model. From the simulation results of two different scenarios, Cox-UMMDR seems to be the most robust when there is either any covariate effect adjusting for or any SNP that has a main effect on the survival phenotype. In addition, the power of all methods decreased as the censoring fraction increased from 0.1 to 0.3, as heritability increased. The power of all methods seems to be greater under MAF\u2009=\u20090.2 than under MAF\u2009=\u20090.4. For illustration, both KM-UMMDR and Cox2-UMMDR were applied to identify SNP by SNP interactions with the survival phenotype to a real dataset of Korean leukemia patients.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Conclusion<\/jats:title>\n                    <jats:p>Both KM-UMMDR and Cox2-UMMDR were easily implemented by combining KM-MDR and Cox-MDR with UM-MDR, respectively, to detect significant gene-gene interactions associated with survival time without cross-validation and permutation testing. The simulation results demonstrate the utility of KM-UMMDR, Cox2-UMMDR and Cox-UMMDR, which outperforms Cox-MDR and KM-MDR when some SNPs with only marginal effects might mask the detection of causal epistasis. In addition, Cox-UMMDR, Cox2-UMMDR and Cox-MDR performed better than KM-UMMDR and KM-MDR when there were potentially confounding covariate effects.<\/jats:p>\n                  <\/jats:sec>","DOI":"10.1186\/s13040-021-00248-9","type":"journal-article","created":{"date-parts":[[2021,3,1]],"date-time":"2021-03-01T15:02:58Z","timestamp":1614610978000},"update-policy":"https:\/\/doi.org\/10.1007\/springer_crossmark_policy","source":"Crossref","is-referenced-by-count":3,"title":["A comparative study on the unified model based multifactor dimensionality reduction methods for identifying gene-gene interactions associated with the survival phenotype"],"prefix":"10.1186","volume":"14","author":[{"given":"Jung Wun","family":"Lee","sequence":"first","affiliation":[]},{"ORCID":"https:\/\/orcid.org\/0000-0001-7941-8933","authenticated-orcid":false,"given":"Seungyeoun","family":"Lee","sequence":"additional","affiliation":[]}],"member":"297","published-online":{"date-parts":[[2021,3,1]]},"reference":[{"issue":"5720","key":"248_CR1","doi-asserted-by":"publisher","first-page":"385","DOI":"10.1126\/science.1109557","volume":"308","author":"RJ Klein","year":"2005","unstructured":"Klein RJ, Zeiss C, Chew EY, Tsai JY, Sackler RS, Haynes C, Henning AK, SanGiovanni JP, Mane SM, Mayne ST, Bracken MB, Ferris FL, Ott J, Narnstable C, Hoh J. Complement factor H polymorphism in age-related macular degeneration. Science. 2005;308(5720):385\u20139.","journal-title":"Science"},{"key":"248_CR2","doi-asserted-by":"publisher","first-page":"88","DOI":"10.1080\/07853890252953473","volume":"34","author":"JH Moore","year":"2002","unstructured":"Moore JH, Williams SW. New strategies for identifying gene-gene interactions in hypertension. Ann Med. 2002;34:88\u201395.","journal-title":"Ann Med"},{"key":"248_CR3","doi-asserted-by":"publisher","first-page":"166","DOI":"10.1056\/NEJMra0905980","volume":"363","author":"TA Manolio","year":"2010","unstructured":"Manolio TA. Genome-wide association studies and assessment of the risk of disease. N Engl J Med. 2010;363:166\u201376.","journal-title":"N Engl J Med"},{"issue":"6","key":"248_CR4","doi-asserted-by":"publisher","first-page":"446","DOI":"10.1038\/nrg2809","volume":"11","author":"EE Eichler","year":"2010","unstructured":"Eichler EE, Flint J, Gibson G, Kong A, Leal SM, Moore JH, Nadeau JH. Missing heritability and strategies for finding the underlying causes of complex disease. Nat Rev Genet. 2010;11(6):446\u201350.","journal-title":"Nat Rev Genet"},{"key":"248_CR5","doi-asserted-by":"publisher","first-page":"138","DOI":"10.1086\/321276","volume":"69","author":"MD Ritchie","year":"2001","unstructured":"Ritchie MD, Hahn LW, Roodi N, Bailey LR, Dupont WD, Parl FF, Moore JH. Multifactor-dimensionality reduction reveals high-order interactions among estrogen-metabolism genes in sporadic breast cancer. Am J Hum Genet. 2001;69:138\u201347.","journal-title":"Am J Hum Genet"},{"key":"248_CR6","doi-asserted-by":"publisher","first-page":"101","DOI":"10.1007\/s00439-010-0905-5","volume":"129","author":"J Gui","year":"2011","unstructured":"Gui J, Moore JH, Kelsey KT, Marsit CJ, Karagas MR, Andrew AS. A novel survival multifactor dimensionality reduction method for detecting gene-gene interactions with application to bladder cancer prognosis. Hum Genet. 2011;129:101\u201310.","journal-title":"Hum Genet"},{"key":"248_CR7","doi-asserted-by":"publisher","first-page":"i582","DOI":"10.1093\/bioinformatics\/bts415","volume":"28","author":"SY Lee","year":"2012","unstructured":"Lee SY, Kwon MS, Oh JM, Park T. Gene-gene interaction analysis for the survival phenotype based on the COX model. Bioinformatics. 2012;28:i582\u20138.","journal-title":"Bioinformatics."},{"issue":"5","key":"248_CR8","first-page":"1","volume":"25","author":"JS Oh","year":"2014","unstructured":"Oh JS, Lee SY. An extension of multifactor dimensionality reduction method for detecting gene-gene interactions with the survival time. J. Korean Data Information Sci Soc. 2014;25(5):1\u201311.","journal-title":"Korean Data Information Sci Soc"},{"key":"248_CR9","first-page":"1","volume":"Article ID 5282","author":"M Park","year":"2020","unstructured":"Park M, Lee JW, Park T, Lee SY. Gene-gene interaction analysis for the survival phenotype based on the Kaplan-Meier median estimate. BioMed Res Int. 2020;Article ID 5282345:1\u201310.","journal-title":"BioMed Res Int"},{"key":"248_CR10","doi-asserted-by":"publisher","first-page":"i605","DOI":"10.1093\/bioinformatics\/btw424","volume":"32","author":"W Yu","year":"2016","unstructured":"Yu W, Lee SY, Park T. A unified model based multifactor dimensionality reduction framework for detecting gene-gene interactions. Bioinformatics. 2016;32:i605\u201310.","journal-title":"Bioinformatics."},{"issue":"27","key":"248_CR11","first-page":"1","volume":"11","author":"SY Lee","year":"2018","unstructured":"Lee SY, Son DH, Kim YK, Yu W, Park T. Unified cox model based multifactor dimensionality reduction method for gene-gene interaction analysis of the survival phenotype. BioData Mining. 2018;11(27):1\u201313.","journal-title":"BioData Mining"},{"key":"248_CR12","doi-asserted-by":"publisher","first-page":"306","DOI":"10.1002\/gepi.20211","volume":"31","author":"DR Velez","year":"2007","unstructured":"Velez DR, White BC, Motsinger AA, Bush WS, Ritchie MD, Williams SM, et al. A balanced accuracy metric for epistasis modeling in imbalanced datasets using multifactor dimensionality reduction. Genet Epidemiol. 2007;31:306\u201315.","journal-title":"Genet Epidemiol"},{"key":"248_CR13","first-page":"3","volume":"11","author":"T Kido","year":"2018","unstructured":"Kido T, Sikora-Wohlfeld W, Kawashima M, Kikuchi S, Kamatani N, Patwardhan A, Chen R, Sirota M, Kodama K, Hadley D, Butte AJ. Are minor alleles more likely to be risk alleles? BMC Med Genet. 2018;11:3.","journal-title":"BMC Med Genet"},{"issue":"11","key":"248_CR14","doi-asserted-by":"publisher","first-page":"1394","DOI":"10.1001\/jamaneurol.2014.1491","volume":"71","author":"AC Naj","year":"2014","unstructured":"Naj AC, Jun G, Reitz C, Kunkle BW, Pericak-Vance MA, et al. Age-at-onset in late onset Alzheimer disease is modified by multiple genetic loci. JAMA Neurol. 2014;71(11):1394\u2013404.","journal-title":"JAMA Neurol"},{"issue":"3","key":"248_CR15","doi-asserted-by":"publisher","first-page":"515","DOI":"10.1093\/biomet\/81.3.515","volume":"81","author":"P Grambsch","year":"1994","unstructured":"Grambsch P, Therneau TM. Proportional hazards tests and diagnostics based on weighted residuals. Biometrika. 1994;81(3):515\u201326.","journal-title":"Biometrika"},{"key":"248_CR16","doi-asserted-by":"publisher","first-page":"1547","DOI":"10.1002\/ijc.25175","volume":"127","author":"M Timofeeva","year":"2010","unstructured":"Timofeeva M, Kropp S, Sauter W, Beckmann L, Rosenberger A, Illig T, Wichmann HE, et al. Genetic polymorphisms of MPO, GSTT1, GSTM1, GSTP1, EPHX1 and NQO1 as risk factors of early-onset lung cancer. Int J Cancer. 2010;127:1547\u201361.","journal-title":"Int J Cancer"}],"container-title":["BioData Mining"],"original-title":[],"language":"en","link":[{"URL":"http:\/\/link.springer.com\/content\/pdf\/10.1186\/s13040-021-00248-9.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"text-mining"},{"URL":"http:\/\/link.springer.com\/article\/10.1186\/s13040-021-00248-9\/fulltext.html","content-type":"text\/html","content-version":"vor","intended-application":"text-mining"},{"URL":"http:\/\/link.springer.com\/content\/pdf\/10.1186\/s13040-021-00248-9.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2021,3,1]],"date-time":"2021-03-01T15:04:33Z","timestamp":1614611073000},"score":1,"resource":{"primary":{"URL":"https:\/\/biodatamining.biomedcentral.com\/articles\/10.1186\/s13040-021-00248-9"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2021,3,1]]},"references-count":16,"journal-issue":{"issue":"1","published-print":{"date-parts":[[2021,12]]}},"alternative-id":["248"],"URL":"https:\/\/doi.org\/10.1186\/s13040-021-00248-9","relation":{"has-preprint":[{"id-type":"doi","id":"10.21203\/rs.3.rs-90733\/v2","asserted-by":"object"},{"id-type":"doi","id":"10.21203\/rs.3.rs-90733\/v1","asserted-by":"object"}]},"ISSN":["1756-0381"],"issn-type":[{"value":"1756-0381","type":"electronic"}],"subject":[],"published":{"date-parts":[[2021,3,1]]},"assertion":[{"value":"7 October 2020","order":1,"name":"received","label":"Received","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"11 February 2021","order":2,"name":"accepted","label":"Accepted","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"1 March 2021","order":3,"name":"first_online","label":"First Online","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"Not applicable.","order":1,"name":"Ethics","group":{"name":"EthicsHeading","label":"Ethics approval and consent to participate"}},{"value":"Not applicable.","order":2,"name":"Ethics","group":{"name":"EthicsHeading","label":"Consent for publication"}},{"value":"The authors declare that they have no competing\u00a0 interests.","order":3,"name":"Ethics","group":{"name":"EthicsHeading","label":"Competing\u00a0 interests"}}],"article-number":"17"}}