Objective Obesity has become a leading preventable cause of morbidity and mortality in many parts of the world. It is thought to originate from multiple genetic and environmental determinants. The aim of the current study was to introduce haplotype-based multi-locus stepwise regression (MSR) as a method to investigate combinations of unlinked single nucleotide polymorphisms (SNPs) for obesity phenotypes. Methods In 2,122 healthy randomly selected men and women of the EPIC-Potsdam cohort, the association between 41 SNPs from 18 obesity-candidate genes and either body mass index (BMI, mean = 25.9 kg/m2, SD = 4.1) or waist circumference (WC, mean = 85.2 cm, SD = 12.6) was assessed. Single SNP analyses were done by using linear regression adjusted for age, sex, and other covariates. Subsequently, MSR was applied to search for the ‘best’ SNP combinations. Combinations were selected according to specific AICc and p-value criteria. Model uncertainty was accounted for by a permutation test. Results The strongest single SNP effects on BMI were found for TBC1D1 rs637797 (β = ？0.33, SE = 0.13), FTO rs9939609 (β = 0.28, SE = 0.13), MC4R rs17700144 (β = 0.41, SE = 0.15), and MC4R rs10871777 (β = 0.34, SE = 0.14). All these SNPs showed similar effects on waist circumference. The two ‘best’ six-SNP combinations for BMI (global p-value = 3.45？10–6 and 6.82？10–6) showed effects ranging from ？1.70 (SE = 0.34) to 0.74 kg/m2 (SE = 0.21) per allele combination. We selected two six-SNP combinations on waist circumference (global p-value = 7.80？10–6 and 9.76？10–6) with an allele combination effect of ？2.96 cm (SE = 0.76) at maximum. Additional adjustment for BMI revealed 15 three-SNP combinations (global p-values ranged from 3.09？10–4 to 1.02？10–2). However, after carrying out the permutation test all SNP combinations lost significance indicating that the statistical associations might have occurred by chance. Conclusion MSR provides a tool to search for risk-related SNP combinations of common traits or diseases. However, the search process does not always find meaningful SNP combinations in a dataset.
References
[1]
World Health Organization (WHO) (2009) Obesity and overweight. Fact Sheet Number 311. Available: http://www.who.int/mediacentre/factsheet？s/fs311/en/index.html. Accessed 2012 May 20.
[2]
Xia Q, Grant SF (2013) The genetics of human obesity. Ann N Y Acad Sci. Epub ahead of print.
[3]
Rankinen T, Zuberi A, Chagnon YC, Weisnagel SJ, Argyropoulos G, et al. (2006) The human obesity gene map: the 2005 update. Obesity (Silver Spring) 14: 529–644.
[4]
Li S, Loos RJ (2008) Progress in the genetics of common obesity: size matters. Curr Opin Lipidol 19: 113–121.
[5]
Speliotes EK, Willer CJ, Berndt SI, Monda KL, Thorleifsson G, et al. (2010) Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index. Nat Genet 42: 937–948.
[6]
Dina C (2008) New insights into the genetics of body weight. Current Opinion in Clinical Nutrition and Metabolic Care 11: 378–384.
[7]
Schaid DJ, Rowland CM, Tines DE, Jacobson RM, Poland GA (2002) Score tests for association between traits and haplotypes when linkage phase is ambiguous. Am J Hum Genet 70: 425–434.
[8]
Zaykin DV, Westfall PH, Young SS, Karnoub MA, Wagner MJ, et al. (2002) Testing association of statistically inferred haplotypes with discrete and continuous traits in samples of unrelated individuals. Hum Hered 53: 79–91.
[9]
Draper N, Smith H (1998) Applied Regression Analysis: Wiley-Interscience.
[10]
Knüppel S, Esparza-Gordillo J, Marenholz I, Holzhütter H, Bauerfeind A, et al. (2012) Multi-locus stepwise regression: a haplotype-based algorithm for finding genetic associations applied to atopic dermatitis. BMC Med Genet 13: 8.
[11]
Abo R, Knight S, Wong J, Cox A, Camp NJ (2008) hapConstructor: automatic construction and testing of haplotypes in a Monte Carlo framework. Bioinformatics 24: 2105–2107.
[12]
Boeing H, Korfmann A, Bergmann MM (1999) Recruitment procedures of EPIC-Germany. European Investigation into Cancer and Nutrition. Ann Nutr Metab 43: 205–215.
[13]
Riboli E, Hunt KJ, Slimani N, Ferrari P, Norat T, et al. (2002) European Prospective Investigation into Cancer and Nutrition (EPIC): study populations and data collection. Public Health Nutr 5: 1113–1124.
[14]
Kroke A, Bergmann MM, Lotze G, Jeckel A, Klipstein-Grobusch K, et al. (1999) Measures of quality control in the German component of the EPIC study. European Prospective Investigation into Cancer and Nutrition. Ann Nutr Metab 43: 216–224.
[15]
Herbert A, Gerry NP, McQueen MB, Heid IM, Pfeufer A, et al. (2006) A common genetic variant is associated with adult and childhood obesity. Science 312: 279–283.
[16]
Frayling TM, Timpson NJ, Weedon MN, Zeggini E, Freathy RM, et al. (2007) A common variant in the FTO gene is associated with body mass index and predisposes to childhood and adult obesity. Science 316: 889–894.
[17]
Willer CJ, Speliotes EK, Loos RJ, Li S, Lindgren CM, et al. (2009) Six new loci associated with body mass index highlight a neuronal influence on body weight regulation. Nat Genet 41: 25–34.
[18]
Meyre D, Delplanque J, Chevre JC, Lecoeur C, Lobbens S, et al. (2009) Genome-wide association study for early-onset and morbid adult obesity identifies three new risk loci in European populations. Nat Genet 41: 157–159.
[19]
Eberle D, Clement K, Meyre D, Sahbatou M, Vaxillaire M, et al. (2004) SREBF-1 gene polymorphisms are associated with obesity and type 2 diabetes in French obese and diabetic cohorts. Diabetes 53: 2153–2157.
[20]
Brouillette C, Bosse Y, Perusse L, Gaudet D, Vohl MC (2004) Effect of liver fatty acid binding protein (FABP) T94A missense mutation on plasma lipoprotein responsiveness to treatment with fenofibrate. J Hum Genet 49: 424–432.
[21]
Lindner I, Helwig U, Rubin D, Fischer A, Marten B, et al. (2007) Prostaglandin E synthase 2 (PTGES2) Arg298His polymorphism and parameters of the metabolic syndrome. Mol Nutr Food Res 51: 1447–1451.
[22]
Hani EH, Clement K, Velho G, Vionnet N, Hager J, et al. (1997) Genetic studies of the sulfonylurea receptor gene locus in NIDDM and in morbid obesity among French Caucasians. Diabetes 46: 688–694.
[23]
Masuzaki H, Paterson J, Shinyama H, Morton NM, Mullins JJ, et al. (2001) A transgenic model of visceral obesity and the metabolic syndrome. Science 294: 2166–2170.
[24]
Narisawa S, Huang L, Iwasaki A, Hasegawa H, Alpers DH, et al. (2003) Accelerated fat absorption in intestinal alkaline phosphatase knockout mice. Mol Cell Biol 23: 7525–7530.
[25]
Sun G, Gagnon J, Chagnon YC, Perusse L, Despres JP, et al. (1999) Association and linkage between an insulin-like growth factor-1 gene polymorphism and fat free mass in the HERITAGE Family Study. Int J Obes Relat Metab Disord 23: 929–935.
[26]
Chadt A, Leicht K, Deshmukh A, Jiang LQ, Scherneck S, et al. (2008) Tbc1d1 mutation in lean mouse strain confers leanness and protects from diet-induced obesity. Nat Genet 40: 1354–1359.
[27]
Scherag A, Dina C, Hinney A, Vatin V, Scherag S, et al. (2010) Two new Loci for body-weight regulation identified in a joint analysis of genome-wide association studies for early-onset extreme obesity in French and german study groups. PLoS Genet 6: e1000916.
[28]
Grant SF, Bradfield JP, Zhang H, Wang K, Kim CE, et al. (2009) Investigation of the locus near MC4R with childhood obesity in Americans of European and African ancestry. Obesity (Silver Spring) 17: 1461–1465.
[29]
Rebhan M, Chalifa-Caspi V, Prilusky J (1997) GeneCards: integrating information about genes, proteins and diseases. Trends Genet 13: 163. Available: http://www.genecards.org. Accessed 2013 May 20.
[30]
Hampe J, Wollstein A, Lu T, Frevel HJ, Will M, et al. (2001) An integrated system for high throughput TaqMan based SNP genotyping. Bioinformatics 17: 654–655.
[31]
Excoffier L, Slatkin M (1995) Maximum-likelihood estimation of molecular haplotype frequencies in a diploid population. Mol Biol Evol 12: 921–927.
[32]
Sugiura N (1978) Further analysts of the data by akaike’ s information criterion and the finite corrections. Communications in Statistics, Theory and Methods 7: 13–26.
[33]
Efron B, Tibshirani RJ (1993) An introduction to the bootstrap: Chapman and Hall.
[34]
R Development Core Team (2011) R: A Language and Environment for Statistical Computing. Vienna, Austria.
[35]
Fisher E, Schulze MB, Stefan N, Haring HU, Doring F, et al. (2009) Association of the FTO rs9939609 single nucleotide polymorphism with C-reactive protein levels. Obesity (Silver Spring) 17: 330–334.
[36]
Curtin K, Wolff RK, Herrick JS, Abo R, Slattery ML (2010) Exploring multilocus associations of inflammation genes and colorectal cancer risk using hapConstructor. BMC Med Genet 11: 170.
[37]
Zuk O, Hechter E, Sunyaev SR, Lander ES (2012) The mystery of missing heritability: Genetic interactions create phantom heritability. Proc Natl Acad Sci USA 109: 1193–1198.
