[en] We propose a minimal protocol for exhaustive genome-wide association interaction analysis that involves screening for epistasis over large-scale genomic data combining strengths of different methods and statistical tools. The different steps of this protocol are illustrated on a real-life data application for Alzheimer's disease (AD) (2259 patients and 6017 controls from France). Particularly, in the exhaustive genome-wide epistasis screening we identified AD-associated interacting SNPs-pair from chromosome 6q11.1 (rs6455128, the KHDRBS2 gene) and 13q12.11 (rs7989332, the CRYL1 gene) (p = 0.006, corrected for multiple testing). A replication analysis in the independent AD cohort from Germany (555 patients and 824 controls) confirmed the discovered epistasis signal (p = 0.036). This signal was also supported by a meta-analysis approach in 5 independent AD cohorts that was applied in the context of epistasis for the first time. Transcriptome analysis revealed negative correlation between expression levels of KHDRBS2 and CRYL1 in both the temporal cortex (β = -0.19, p = 0.0006) and cerebellum (β = -0.23, p < 0.0001) brain regions. This is the first time a replicable epistasis associated with AD was identified using a hypothesis free screening approach.
Disciplines :
Genetics & genetic processes
Author, co-author :
Gusareva, Elena ; Université de Liège - ULiège > Dép. d'électric., électron. et informat. (Inst.Montefiore) > Bioinformatique
Carrasquillo, Minerva M.; Mayo Clinic Florida > Department of Neuroscience
Bellenguez, Céline; Unité d'Epidémiologie et de Santé Publique > INSERM U744
Cuyvers, Elise; Institute Born-Bunge, University of Antwerp > Department of Molecular Genetics, VIB
Colon, Samuel; Mayo Clinic Florida > Department of Neuroscience
Graff-Radford, Neill R.; Mayo Clinic Florida > Department of Neurology
Petersen, Ronald C.; Mayo Clinic Florida > Department of Neurology
Dickson, Dennis W.; Mayo Clinic Florida > Department of Neuroscience
Mahachie John, Jestinah ; Université de Liège - ULiège > Dép. d'électric., électron. et informat. (Inst.Montefiore) > Bioinformatique
Bessonov, Kyrylo ; Université de Liège - ULiège > Dép. d'électric., électron. et informat. (Inst.Montefiore) > Bioinformatique
Van Broeckhoven, Christine; Institute Born-Bunge, University of Antwerp > Department of Molecular Genetics VIB
Harold, Denise; Cardiff University School of Medicine, Institute of Psychological Medicine and Clinical Neurosciences > Medical Research Council Centre for Neuropsychiatric Genetics and Genomics
Williams, Julie; Cardiff University School of Medicine, Institute of Psychological Medicine and Clinical Neurosciences > Medical Research Council Centre for Neuropsychiatric Genetics and Genomics
Amouyel, Philippe; Unité d'Epidémiologie et de Santé Publique > INSERM U744
Sleegers, Kristel; Institute Born-Bunge, University of Antwerp > Department of Molecular Genetics VIB
Ertekin-Taner, Nilüfer; Mayo Clinic Florida > Department of Neurology
Lambert, Jean-Charles; Unité d'Epidémiologie et de Santé Publique > INSERM U744
Van Steen, Kristel ; Université de Liège - ULiège > Dép. d'électric., électron. et informat. (Inst.Montefiore) > Bioinformatique
Allen M., Zou F., Chai H.S., Younkin C.S., Miles R., Nair A.A., Crook J.E., Pankratz V.S., Carrasquillo M.M., Rowley C.N., Nguyen T., Ma L., Malphrus K.G., Bisceglio G., Ortolaza A.I., Palusak R., Middha S., Maharjan S., Georgescu C., Schultz D., Rakhshan F., Kolbert C.P., Jen J., Sando S.B., Aasly J.O., Barcikowska M., Uitti R.J., Wszolek Z.K., Ross O.A., Petersen R.C., Graff-Radford N.R., Dickson D.W., Younkin S.G., Ertekin-Taner N. Glutathione S-transferase omega genes in Alzheimer and Parkinson disease risk, age-at-diagnosis and brain gene expression: an association study with mechanistic implications. Mol. Neurodegener. 2012, 7:13.
Anastassiou D. Computational analysis of the synergy among multiple interacting genes. Mol. Syst. Biol. 2007, 3:83.
Berchtold D., Piccolis M., Chiaruttini N., Riezman I., Riezman H., Roux A., Walther T.C., Loewith R. Plasma membrane stress induces relocalization of Slm proteins and activation of TORC2 to promote sphingolipid synthesis. Nat. Cell Biol. 2012, 14:542-547.
Bettens K., Brouwers N., Van Miegroet H., Gil A., Engelborghs S., De Deyn P.P., Vandenberghe R., Van Broeckhoven C., Sleegers K. Follow-up study of susceptibility loci for Alzheimer's disease and onset age identified by genome-wide association. J.Alzheimer's Dis. 2010, 19:1169-1175.
Bettens K., Sleegers K., Van Broeckhoven C. Genetic insights in Alzheimer's disease. Lancet Neurol. 2013, 12:92-104.
Birnbaum S., Ludwig K.U., Reutter H., Herms S., Steffens M., Rubini M., Baluardo C., Ferrian M., Almeida de Assis N., Alblas M.A., Barth S., Freudenberg J., Lauster C., Schmidt G., Scheer M., Braumann B., Berge S.J., Reich R.H., Schiefke F., Hemprich A., Potzsch S., Steegers-Theunissen R.P., Potzsch B., Moebus S., Horsthemke B., Kramer F.J., Wienker T.F., Mossey P.A., Propping P., Cichon S., Hoffmann P., Knapp M., Nothen M.M., Mangold E. Key susceptibility locus for nonsyndromic cleft lip with or without cleft palate on chromosome 8q24. Nature Genet. 2009, 41:473-477.
Cattaert T., Calle M.L., Dudek S.M., Mahachie John J.M., Van Lishout F., Urrea V., Ritchie M.D., Van Steen K. Model-based multifactor dimensionality reduction for detecting epistasis in case-control data in the presence of noise. Ann. Hum. Genet. 2011, 75:78-89.
Dickson R.C. Thematic review series: sphingolipids. New insights into sphingolipid metabolism and function in budding yeast. J.Lipid Res. 2008, 49:909-921.
Harold D., Abraham R., Hollingworth P., Sims R., Gerrish A., Hamshere M.L., Pahwa J.S., Moskvina V., Dowzell K., Williams A., Jones N., Thomas C., Stretton A., Morgan A.R., Lovestone S., Powell J., Proitsi P., Lupton M.K., Brayne C., Rubinsztein D.C., Gill M., Lawlor B., Lynch A., Morgan K., Brown K.S., Passmore P.A., Craig D., McGuinness B., Todd S., Holmes C., Mann D., Smith A.D., Love S., Kehoe P.G., Hardy J., Mead S., Fox N., Rossor M., Collinge J., Maier W., Jessen F., Schurmann B., van den Bussche H., Heuser I., Kornhuber J., Wiltfang J., Dichgans M., Frolich L., Hampel H., Hull M., Rujescu D., Goate A.M., Kauwe J.S., Cruchaga C., Nowotny P., Morris J.C., Mayo K., Sleegers K., Bettens K., Engelborghs S., De Deyn P.P., Van Broeckhoven C., Livingston G., Bass N.J., Gurling H., McQuillin A., Gwilliam R., Deloukas P., Al-Chalabi A., Shaw C.E., Tsolaki M., Singleton A.B., Guerreiro R., Muhleisen T.W., Nothen M.M., Moebus S., Jockel K.H., Klopp N., Wichmann H.E., Carrasquillo M.M., Pankratz V.S., Younkin S.G., Holmans P.A., O'Donovan M., Owen M.J., Williams J. Genome-wide association study identifies variants at CLU and PICALM associated with Alzheimer's disease. Nat. Genet. 2009, 41:1088-1093.
Hillmer A.M., Brockschmidt F.F., Hanneken S., Eigelshoven S., Steffens M., Flaquer A., Herms S., Becker T., Kortum A.K., Nyholt D.R., Zhao Z.Z., Montgomery G.W., Martin N.G., Muhleisen T.W., Alblas M.A., Moebus S., Jockel K.H., Brocker-Preuss M., Erbel R., Reinartz R., Betz R.C., Cichon S., Propping P., Baur M.P., Wienker T.F., Kruse R., Nothen M.M. Susceptibility variants for male-pattern baldness on chromosome 20p11. Nature Genet. 2008, 40:1279-1281.
Lambert J.C., Heath S., Even G., Campion D., Sleegers K., Hiltunen M., Combarros O., Zelenika D., Bullido M.J., Tavernier B., Letenneur L., Bettens K., Berr C., Pasquier F., Fievet N., Barberger-Gateau P., Engelborghs S., De Deyn P., Mateo I., Franck A., Helisalmi S., Porcellini E., Hanon O., de Pancorbo M.M., Lendon C., Dufouil C., Jaillard C., Leveillard T., Alvarez V., Bosco P., Mancuso M., Panza F., Nacmias B., Bossu P., Piccardi P., Annoni G., Seripa D., Galimberti D., Hannequin D., Licastro F., Soininen H., Ritchie K., Blanche H., Dartigues J.F., Tzourio C., Gut I., Van Broeckhoven C., Alperovitch A., Lathrop M., Amouyel P. Genome-wide association study identifies variants at CLU and CR1 associated with Alzheimer's disease. Nat. Genet. 2009, 41:1094-1099. European Alzheimer's Disease Initiative, I.
Mahachie John J.M., Cattaert T., Lishout F.V., Gusareva E.S., Steen K.V. Lower-order effects adjustment in quantitative traits model-based multifactor dimensionality reduction. PLoS One 2012, 7:e29594.
McKhann G., Drachman D., Folstein M., Katzman R., Price D., Stadlan E.M. Clinical diagnosis of Alzheimer's disease: report of the NINCDS-ADRDA Work Group under the auspices of Department of Health and Human Services Task Force on Alzheimer's Disease. Neurology 1984, 34:939-944.
Miko I. Epistasis: gene interaction and phenotype effects. Nature Education 2008, 1:197.
Moore J.H. Aglobal view of epistasis. Nat. Genet. 2005, 37:13-14.
Moore J.H., Barney N., Tsai C.T., Chiang F.T., Gui J., White B.C. Symbolic modeling of epistasis. Hum. Hered. 2007, 63:120-133.
Phillips P.C. Epistasis-the essential role of gene interactions in the structure and evolution of genetic systems. Nat. Rev. Genet. 2008, 9:855-867.
Pierce A., Podlutskaya N., Halloran J.J., Hussong S.A., Lin P.Y., Burbank R., Strong R., Richardson A., Hart M.J., Galvan V. Overexpression of heat shock factor 1 phenocopies the effect of chronic inhibition of TOR by rapamycin and is sufficient to ameliorate Alzheimer's-like deficits in mice modeling the disease. J.Neurochem. 2012, 124:880-893.
Sennvik K., Fastbom J., Blomberg M., Wahlund L.O., Winblad B., Benedikz E. Levels of alpha- and beta-secretase cleaved amyloid precursor protein in the cerebrospinal fluid of Alzheimer's disease patients. Neurosci. Lett. 2000, 278:169-172.
Templeton A.R. Epistasis and Complex Traits 2000, Oxford University Press, New York.
Van Steen K. Travelling the world of gene-gene interactions. Brief. Bioinform. 2012, 13:1-19.
Van Steen K., Curran D., Kramer J., Molenberghs G., Van Vreckem A., Bottomley A., Sylvester R. Multicollinearity in prognostic factor analyses using the EORTC QLQ-C30: identification and impact on model selection. Stat. Med. 2002, 21:3865-3884.
Wan X., Yang C., Yang Q., Xue H., Fan X., Tang N.L., Yu W. BOOST: a fast approach to detecting gene-gene interactions in genome-wide case-control studies. Am. J. Hum. Genet. 2010, 87:325-340.
Watkinson J., Anastassiou D. Synergy disequilibrium plots: graphical visualization of pairwise synergies and redundancies of SNPs with respect to a phenotype. Bioinformatics 2009, 25:1445-1446.
Wichmann H.E., Gieger C., Illig T., Group M.K.S. KORA-gen-resource for population genetics, controls and a broad spectrum of disease phenotypes. Gesundheitswesen 2005, 67(suppl 1):S26-S30.
Wistow G. The human crystallin gene families. Hum. Genomics 2012, 6:26.
Zou F., Chai H.S., Younkin C.S., Allen M., Crook J., Pankratz V.S., Carrasquillo M.M., Rowley C.N., Nair A.A., Middha S., Maharjan S., Nguyen T., Ma L., Malphrus K.G., Palusak R., Lincoln S., Bisceglio G., Georgescu C., Kouri N., Kolbert C.P., Jen J., Haines J.L., Mayeux R., Pericak-Vance M.A., Farrer L.A., Schellenberg G.D., Petersen R.C., Graff-Radford N.R., Dickson D.W., Younkin S.G., Ertekin-Taner N. Brain expression genome-wide association study (eGWAS) identifies human disease-associated variants. PLoS Genet. 2012, 8:e1002707. Alzheimer's Disease Genetics, C.