Significant Linkage of BMI to Chromosome 10p in the U.K. Population and Evaluation of GAD2 as a Positional Candidate
http://www.100md.com
糖尿病学杂志 2006年第6期
1 Oxford Centre for Diabetes, Endocrinology and Metabolism, University of Oxford, Oxford, U.K
2 Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, U.K
3 School of Clinical Medical Sciences, University of Newcastle, Newcastle, U.K
4 Department of Diabetes and Metabolic Medicine, Bart’s and The London Queen Mary’s School of Medicine and Dentistry, London, U.K
5 Departments of Medicine and Clinical Biochemistry, Addenbrooke’s Hospital, Cambridge, U.K
6 Centre for Molecular Genetics, Peninsula Medical School, Exeter, U.K
LD, linkage disequilibrium; LOD, logarithm of odds; QTDT, quantitative transmission-disequilibrium test; QTL, quantitative trait locus; SNP, single nucleotide polymorphism
ABSTRACT
Obesity is a major health problem, and many family-based studies have suggested that it has a strong genetic basis. We performed a genome-wide quantitative trait linkage scan for loci influencing BMI in 573 pedigrees from the U.K. We identified genome-wide significant linkage (logarithm of odds = 3.74, between D10S208 and D10S196, genome-wide P = 0.0186) on chromosome 10p. The size of our study population and the statistical significance of our findings provide substantial contributions to the body of evidence for a locus on chromosome 10p. We examined eight single nucleotide polymorphisms (SNPs) in GAD2, which maps to this linkage region, tagging the majority of variation in the gene, and observed marginally significant (0.01 < P < 0.05) associations between four common variants and BMI. However, these SNPs did not account for our evidence of linkage to BMI, and they did not replicate (in direction of effect) the previous associations. We therefore conclude that these SNPs are not the etiological variants underlying this locus. We cannot rule out the possibility that other untagged variations in GAD2 may, in part, be involved, but it is most likely that alternative gene(s) within the broad gene-rich region of linkage on 10p are responsible for variation in body mass and susceptibility to obesity.
Human obesity is a disease of increasing concern, with obese individuals at increased risk of developing many chronic disorders. By 2004, 49 obesity-related Mendelian syndromes had been localized to genomic regions and 204 quantitative trait loci (QTLs) to obesity-related phenotypes identified from 50 genome scans in humans (1). In addition, nearly 100 candidate genes have been implicated in obesity etiology (1). The close epidemiological and pathological relationship between obesity and type 2 diabetes suggests that their etiology may, in part, be shared. We have analyzed a 573 full sibpair pedigree resource previously ascertained and analyzed for susceptibility to type 2 diabetes (2) for loci influencing BMI.
In our dataset of 573 pedigrees, the phenotypes of age and BMI (both taken at time of ascertainment) were available for nine parents (70eC89 years of age) and 1,215 offspring (37eC88 years of age) with type 2 diabetes. Mean (±SD) BMI was 27.8 ± 4.3 kg/m2 in male subjects (n = 657) and 29.8 ± 5.6 kg/m2 in female subjects (n = 567). Measures of BMI were log transformed to reduced skewness and kurtosis. To account for possible effects of age, sex, and selection (due to ascertainment on the basis of type 2 diabetes status), we adjusted and standardized our sample data against U.K. population data obtained through the Health Survey for England 1998 (3). Linkage analysis was performed using the combined "squared sums-squared differences" Haseman-Elston regression approach, implemented in MERLIN-REGRESS (4). This method properly accommodates the selected nature of our data by preventing inflation of the type I error rate and achieving maximal power during linkage analysis; therefore, it has advantages in this situation over variance componentseCbased methods (4). We estimated sibling genotype identity-by-descent coefficients using genotypes generated from 418 autosomal microsatellite markers [mean spacing 9.26 cM(H)] (2). We specified the sex- and age-adjusted population mean (zero), variance (one), and heritability (60%) (5) in our analyses. Power calculations indicate a power of 70% to detect linkage to a locus accounting for 30% of the total phenotypic variance with a logarithm of odds (LOD) score of 2 and 45% power to detect a linkage with an LOD score of 3 (6).
Our genome-wide scan for BMI susceptibility genes identified two linkage peaks with an LOD 1.18: one on chromosome 1p36.13 with an LOD score of 1.38 between markers D1S2697 and D1S199 and a second, much larger peak on chromosome 10p11.22-23 with an LOD of 3.74 between markers D10S208 and D10S196 (Fig. 1). We determined by simulation that the signal on chromosome 10 was significant at the genome-wide level (P = 0.0186) and the point-wise level (P = 0.0001) using an approach that has been previously described (6). Our original genome-wide linkage analysis of type 2 diabetes susceptibility in this dataset (2) identified only marginally significant linkage to this region on 10p, with an LOD of 0.68 (asymptotic point-wise P = 0.039). Conversely, we observed no linkage of BMI (with asymptotic point-wise P < 0.05) to the region on 10q23 in which we had previously observed maximal evidence for linkage to type 2 diabetes. We can therefore be confident that the linkage evidence on 10p reflects a locus influencing BMI susceptibility directly rather than the surrogate effects of a type 2 diabetes susceptibility locus.
Our study adds substantially, in terms of significance of results and sample size, to the growing body of evidence for linkage to this region. Nine previous studies have reported evidence for linkage to BMI or related traits on the short arm or pericentric region of chromosome 10 (7,8,9,10,11,12,13,14,15). The majority of these studies, our own included, demonstrate maximum evidence for linkage in the region 55eC80 cM from the p-terminal of chromosome 10 (Fig. 1B), strongly implicating a locus in this region that influences susceptibility to body mass and obesity.
Several excellent candidate genes fall within this region on chromosome 10, including GAD2, UCN3 (urocortin III), and PPYR1/NPY4R (pancreatic polypeptide receptor 1). GAD2 encodes GAD65, which catalyzes the formation of the neurotransmitter aminobutyric acid (GABA), is synthesized by GABAergic neurones within nuclei of the hypothalamus that are known to be involved in feeding behavior (16), and has been the focus of three recent studies.
The fine-mapping linkage study of adult obesity in the French population by Boutin et al. (13) observed that the maximally linked microsatellite marker in their study, D10S197, was located within intron 7 of GAD2. Boutin et al. identified three single nucleotide polymorphisms (SNPs), a promoter variant (rs2236418, eC243 A>G), and two intronic variants (rs992990, +61450 C>A and rs928197, +83897 T>A) that are associated with obesity. In each case, the rare allele was associated with an obese phenotype. They observed marginally significant association between the variant allele of rs992990 and evidence for linkage (P = 0.02), together with functional evidence implicating the other two SNPs. The second study, also in a French population, observed that the same (G) allele of the promoter SNP (rs2236418) was associated with severe childhood obesity (17). Most recently, Swarbrick et al. (18) investigated the importance of GAD2 in obesity susceptibility in two large case/control studies (in American and Canadian populations) and in a smaller study of German pedigrees and found no statistical evidence supporting the candidacy of GAD2.
Given this prior, albeit conflicting, evidence of the importance of GAD2 in body mass determination, we examined the candidacy of GAD2 in our population. We first considered the three SNPs (rs2236418, rs992990, and rs928197) that showed some evidence of association in the previous studies. The 573 sibships were genotyped for these three SNPs using a competitive allele-specific PCR method (KASPar; Kbioscience, Hertfordshire, U.K.). Genotyping call rates exceeded 95% at each SNP, and no discordant genotypes were observed in 92 duplicate samples.
We found no evidence for population stratification with the three SNPs analyzed (all P 0.5). This fits with previous studies of these case samples when STRUCTURE has been applied to the genome-wide microsatellite data (N. Martin, L. Cardon, personal communication) and a more recent genome-wide analysis using the quantitative transmission-disequilibrium test (QTDT) (19,20) (all P 0.6, corrected for multiple testing). We therefore determined the total evidence for association of each SNP with adjusted and standardized BMI using QTDT, as well as with QPDTPHASE (21), which makes less restrictive assumptions about the data distribution (22). We observed associations of marginal significance with two of these SNPs, rs2236418 (P = 0.025) and rs928197 (P = 0.011) (Table 1). However, in both cases, the minor allele was associated with lower mean BMI, in contrast to the findings of Boutin et al. To capture additional genetic diversity within GAD2, we used HAPLOVIEW (23) and the GAD2 resequencing data from the Environmental Genome Project (www.egp.gs.washington.edu) to select an additional five SNPs for genotyping. The resulting eight-SNP panel tagged all of the haplotypes in the three largest linkage disequilibrium (LD) blocks (blocks 1, 2, and 5, spanning 73% of the total gene sequence) (Fig. 1A). The gene also contained two small blocks (blocks 3 and 4) and the intervening regions of low LD that would have required substantial additional genotyping to tag. These eight SNPs tagged the proximal promoter, the 3' untranslated region, and all exons but one, capturing 93% of all common (>10%) variation with an r2 > 0.5 and 35% of that with an r2 > 0.8 across the gene as a whole. The arrangement of the SNPs, the LD blocks that they tag, and their r2 and D' values are shown in Fig. 2. We tested the five additional SNPs individually for total association with BMI using QTDT and QPDTPHASE, having first confirmed no evidence for population stratification at these markers (all P 0.6). We also tested haplotypes of block-tagging SNPs (Table 2) with QPDTPHASE, pooling haplotypes with frequencies <1%. We found SNP rs2839669 (P = 0.036), located upstream of the associated promoter SNP rs2236418, and rs7071922 (P = 0.031), located in block 2, to be marginally associated, with the minor allele again associated with lower BMI (Table 1). Haplotypes of the block-tagging SNPs were also marginally associated with BMI (Table 2). No global association was seen with microsatellite D10S197 (genotyped in the original genome scan), located proximal to rs7908975 in block 2 (P > 0.1). No SNP showed evidence of dominance (P > 0.1).
We determined the extent to which these associated SNPs accounted for the linkage signal. The evidence for linkage to a quantitative trait is reduced when an association with an SNP is incorporated into the test for linkage. When the associated SNP is in perfect LD with the QTL, or is the QTL itself, the evidence for linkage can be eliminated completely (24). We modeled each associated SNP as a covariate, with genotype "11" coded as eC1, "12" as 0, and "22" as +1, consistent with a standard biometrical model and the gene-dosing approach used by QTDT (20). We also modeled the three sets of block-tagging SNPs as a set of covariates representing each individual SNP. We regressed each individual’s GAD2 covariate(s) from BMI using coefficients calculated with QTDT. To allow a meaningful interpretation of these GAD2 covariate analyses, we obtained the appropriate baseline multipoint LOD scores for each SNP in question, that is, the evidence for linkage alone, without accounting for the association, restricted to those family members with a successful GAD2 genotype and covariate.
For these analyses, we used the latest version of MERLIN-REGRESS, which accommodates the effects of LD between SNPs (known to inflate linkage evidence if not properly accommodated) by grouping them into clusters for the purposes of calculating allele-sharing probabilities (25). However, this method assumes no recombination within clusters and no LD between clusters. To accommodate the high multiallelic D' (i.e., nonzero LD) between the haplotype blocks in GAD2, we used two alternative clustering approaches. In the first, we treated the eight SNPs as three separate clusters (representing the three blocks they were chosen to tag), with the intervening microsatellite (D10S197, falling within block 2) considered as part of the second cluster. In the second, we grouped all SNPs (together with D10S197) into one single cluster.
We found that the baseline LOD scores differed according to the clustering approach used to model LD (Table 3, columns 2 and 4). However, the effects of the GAD2 covariates themselves did not differ appreciably between clustering methods (Table 3, columns 3 and 5). All four associated SNPs, and the haplotypes containing them, resulted in very minor LOD score reductions (LOD < 0.6; Table 3). These results suggest that not one of the SNPs tested is the QTL responsible for the linkage to BMI.
In conclusion, we have detected genome-wide significant linkage to BMI on the pericentric region of chromosome 10 in U.K. pedigrees with type 2 diabetes. This is the largest dataset to observe linkage to this region, so our study substantially adds to the body of evidence implicating loci in this region in BMI variation. We have examined the candidacy of GAD2 as the gene underlying this linkage, following two previous studies reporting association. We observed marginally significant associations to four SNPs in GAD2. Of the three SNPs in common with the Boutin study, we only observed associations with rs2236418 (P = 0.025) and rs928197 (P = 0.011), but in the opposite direction. Differences in the underlying LD structure of chromosome 10 in U.K. and French populations are highly unlikely explanations for this reversal of association. It is, in principle at least, conceivable that differences in genetic background, patterns of interacting environmental exposures, and clinical features between the two samples could explain a reversal of association direction. However, such differences would have to be extreme, which seems implausible given that these are otherwise quite similar European samples. We note that the Boutin study examined extreme obesity in a case/control setting, whereas we have quantitatively examined a continuous distribution of BMI measurements, only a very small proportion of which could be described as extremely obese. While such differences in study design might in theory explain a failure to reproduce the finding of an effect, it is difficult to see how these differences could explain a change in the direction of an effect. We conclude that our findings do not constitute replication (defined as same allele/haplotype, same phenotype) of the associations seen in previous studies.
Given the number of SNPs that we have tested, and the low correlation between them, our association findings may reflect chance. The alternative is that the associations we observe are true and reflect modest LD with other etiological variants. While it is possible that these variants lie in GAD2, as we have not exhaustively surveyed all genetic variation therein, it is more likely that they are located elsewhere in the 30-Mb region encompassing the consensus of linkage peaks.
ACKNOWLEDGMENTS
The Warren 2 family collection was funded by Diabetes UK through the Warren 2 Bequest and the genotyping through project grant RD02/0002385. S.W. is funded by a Wellcome Trust Career Development Fellowship.
We thank Drs. Stacey Cherny and Javier Gayan (Wellcome Trust Centre for Human Genetics) for useful discussions. Prof. Paola Primatesta (University College London) kindly made available population data from the Health Survey for England.
FOOTNOTES
DOI: 10.2337/db05-1674
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
REFERENCES
Perusse L, Rankinen T, Zuberi A, Chagnon YC, Weisnagel SJ, Argyropoulos G, Walts B, Snyder EE, Bouchard C: The human obesity gene map: the 2004 update. Obes Res 13:381eC490, 2005
Wiltshire S, Hattersley AT, Hitman GA, Walker M, Levy JC, Sampson M, O’Rahilly S, Frayling TM, Bell JI, Lathrop GM, Bennett A, Dhillon R, Fletcher C, Groves CJ, Jones E, Prestwich P, Simecek N, Rao PV, Wishart M, Bottazzo GF, Foxon R, Howell S, Smedley D, Cardon LR, Menzel S, McCarthy MI: A genomewide scan for loci predisposing to type 2 diabetes in a U.K. population (the Diabetes UK Warren 2 Repository): analysis of 573 pedigrees provides independent replication of a susceptibility locus on chromosome 1q. Am J Hum Genet 69:553eC569, 2001
Erens B, Primatesta P: Health Survey for England 1998. London, The Stationary Office, 2001
Sham PC, Purcell S, Cherny SS, Abecasis GR: Powerful regression-based quantitative trait linkage analysis of general pedigrees. Am J Hum Genet 71:238eC253, 2002
Schousboe K, Visscher PM, Erbas B, Kyvik KO, Hopper JL, Henriksen JE, Heitmann BL, Sorensen TI: Twin study of genetic and environmental influences on adult body size, shape, and composition. Int J Obes Relat Metab Disord 28:39eC48, 2004
Wiltshire S, Frayling TM, Groves CJ, Levy JC, Hitman GA, Sampson M, Walker M, Menzel S, Hattersley AT, Cardon LR, McCarthy MI: Evidence from a large U.K. family collection that genes influencing age of onset of type 2 diabetes map to chromosome 12p and to the MODY3/NIDDM2 locus on 12q24. Diabetes 53:855eC860, 2004
Hager J, Dina C, Francke S, Dubois S, Houari M, Vatin V, Vaillant E, Lorentz N, Basdevant A, Clement K, Guy-Grand B, Froguel P: A genome-wide scan for human obesity genes reveals a major susceptibility locus on chromosome 10. Nat Genet 20:304eC308, 1998
Hinney A, Ziegler A, Oeffner F, Wedewardt C, Vogel M, Wulftange H, Geller F, Stubing K, Siegfried W, Goldschmidt HP, Remschmidt H, Hebebrand J: Independent confirmation of a major locus for obesity on chromosome 10. J Clin Endocrinol Metab 85:2962eC2965, 2000 [erratum in J Clin Endocrinol Metab 85:3972, 2000]
Chagnon YC, Rice T, Perusse L, Borecki IB, Ho-Kim MA, Lacaille M, Pare C, Bouchard L, Gagnon J, Leon AS, Skinner JS, Wilmore JH, Rao DC, Bouchard C: HERITAGE Family Study: genomic scan for genes affecting body composition before and after training in Caucasians from HERITAGE. J Appl Physiol 90:1777eC1787, 2001
Hsueh WC, Mitchell BD, Schneider JL, St Jean PL, Pollin TI, Ehm MG, Wagner MJ, Burns DK, Sakul H, Bell CJ, Shuldiner AR: Genome-wide scan of obesity in the Old Order Amish. J Clin Endocrinol Metab 86:1199eC1205, 2001
Lindsay RS, Kobes S, Knowler WC, Bennett PH, Hanson RL: Genome-wide linkage analysis assessing parent-of-origin effects in the inheritance of type 2 diabetes and BMI in Pima Indians. Diabetes 50:2850eC2857, 2001
Price RA, Li WD, Bernstein A, Crystal A, Golding EM, Weisberg SJ, Zuckerman WA: A locus affecting obesity in human chromosome region 10p12. Diabetologia 44:363eC366, 2001
Boutin P, Dina C, Vasseur F, Dubois S, Corset L, Seron K, Bekris L, Cabellon J, Neve B, Vasseur-Delannoy V, Chikri M, Charles MA, Clement K, Lernmark A, Froguel P: GAD2 on chromosome 10p12 is a candidate gene for human obesity. PLoS Biol 1:E68, 2003
Dong C, Wang S, Li WD, Li D, Zhao H, Price RA: Interacting genetic loci on chromosomes 20 and 10 influence extreme human obesity. Am J Hum Genet 72:115eC124, 2003
Gorlova OY, Amos CI, Wang NW, Shete S, Turner ST, Boerwinkle E: Genetic linkage and imprinting effects on body mass index in children and young adults. Eur J Hum Genet 11:425eC432, 2003
Kalra SP, Dube MG, Pu S, Xu B, Horvath TL, Kalra PS.: Interacting appetite-regulating pathways in the hypothalamic regulation of body weight. Endocr Rev 20:68eC100, 1999
Meyre D, Boutin P, Tounian A, Deweirder M, Aout M, Jouret B, Heude B, Weill J, Tauber M, Tounian P, Froguel P: Is glutamate decarboxylase 2 (GAD2) a genetic link between low birth weight and subsequent development of obesity in children J Clin Endocrinol Metab 90:2384eC2390, 2005
Swarbrick MM, Waldenmaier B, Pennacchio LA, Lind DL, Cavazos MM, Geller F, Merriman R, Ustaszewska A, Malloy M, Scherag A, Hsueh WC, Rief W, Mauvais-Jarvis F, Pullinger CR, Kane JP, Dent R, McPherson R, Kwok PY, Hinney A, Hebebrand J, Vaisse C: Lack of support for the association between GAD2 polymorphisms and sever human obesity. PLOS Biol 3:e315, 2005
Abecasis GR, Cardon LR, Cookson WO: A general test of association for quantitative traits in nuclear families. Am J Hum Genet 66:279eC292, 2000
Abecasis GR, Cookson WO, Cardon LR: Pedigree tests of transmission disequilibrium. Eur J Hum Genet 8:545eC551, 2000
Dudbridge F: Pedigree disequilibrium tests for multilocus haplotypes. Genet Epidemiol 25:115eC121, 2003
Monks SA, Kaplan NL: Removing the sampling restrictions from family-based tests of association for a quantitative trait locus. Am J Hum Genet 66:576eC592, 2000
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualisation of LD and haplotype maps. Bioinformatics 21:263eC265, 2005
Fulker DW, Cherny SS, Sham PC, Hewitt JK: Combined linkage and association sib-pair analysis for quantitative traits. Am J Hum Genet 64:259eC267, 1999
Abecasis GR, Wigginton JE: Handling marker-marker linkage disequilibrium: pedigree analysis with clustered markers. Am J Hum Genet 77:754eC767, 2005(Christopher J. Groves, El)
2 Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, U.K
3 School of Clinical Medical Sciences, University of Newcastle, Newcastle, U.K
4 Department of Diabetes and Metabolic Medicine, Bart’s and The London Queen Mary’s School of Medicine and Dentistry, London, U.K
5 Departments of Medicine and Clinical Biochemistry, Addenbrooke’s Hospital, Cambridge, U.K
6 Centre for Molecular Genetics, Peninsula Medical School, Exeter, U.K
LD, linkage disequilibrium; LOD, logarithm of odds; QTDT, quantitative transmission-disequilibrium test; QTL, quantitative trait locus; SNP, single nucleotide polymorphism
ABSTRACT
Obesity is a major health problem, and many family-based studies have suggested that it has a strong genetic basis. We performed a genome-wide quantitative trait linkage scan for loci influencing BMI in 573 pedigrees from the U.K. We identified genome-wide significant linkage (logarithm of odds = 3.74, between D10S208 and D10S196, genome-wide P = 0.0186) on chromosome 10p. The size of our study population and the statistical significance of our findings provide substantial contributions to the body of evidence for a locus on chromosome 10p. We examined eight single nucleotide polymorphisms (SNPs) in GAD2, which maps to this linkage region, tagging the majority of variation in the gene, and observed marginally significant (0.01 < P < 0.05) associations between four common variants and BMI. However, these SNPs did not account for our evidence of linkage to BMI, and they did not replicate (in direction of effect) the previous associations. We therefore conclude that these SNPs are not the etiological variants underlying this locus. We cannot rule out the possibility that other untagged variations in GAD2 may, in part, be involved, but it is most likely that alternative gene(s) within the broad gene-rich region of linkage on 10p are responsible for variation in body mass and susceptibility to obesity.
Human obesity is a disease of increasing concern, with obese individuals at increased risk of developing many chronic disorders. By 2004, 49 obesity-related Mendelian syndromes had been localized to genomic regions and 204 quantitative trait loci (QTLs) to obesity-related phenotypes identified from 50 genome scans in humans (1). In addition, nearly 100 candidate genes have been implicated in obesity etiology (1). The close epidemiological and pathological relationship between obesity and type 2 diabetes suggests that their etiology may, in part, be shared. We have analyzed a 573 full sibpair pedigree resource previously ascertained and analyzed for susceptibility to type 2 diabetes (2) for loci influencing BMI.
In our dataset of 573 pedigrees, the phenotypes of age and BMI (both taken at time of ascertainment) were available for nine parents (70eC89 years of age) and 1,215 offspring (37eC88 years of age) with type 2 diabetes. Mean (±SD) BMI was 27.8 ± 4.3 kg/m2 in male subjects (n = 657) and 29.8 ± 5.6 kg/m2 in female subjects (n = 567). Measures of BMI were log transformed to reduced skewness and kurtosis. To account for possible effects of age, sex, and selection (due to ascertainment on the basis of type 2 diabetes status), we adjusted and standardized our sample data against U.K. population data obtained through the Health Survey for England 1998 (3). Linkage analysis was performed using the combined "squared sums-squared differences" Haseman-Elston regression approach, implemented in MERLIN-REGRESS (4). This method properly accommodates the selected nature of our data by preventing inflation of the type I error rate and achieving maximal power during linkage analysis; therefore, it has advantages in this situation over variance componentseCbased methods (4). We estimated sibling genotype identity-by-descent coefficients using genotypes generated from 418 autosomal microsatellite markers [mean spacing 9.26 cM(H)] (2). We specified the sex- and age-adjusted population mean (zero), variance (one), and heritability (60%) (5) in our analyses. Power calculations indicate a power of 70% to detect linkage to a locus accounting for 30% of the total phenotypic variance with a logarithm of odds (LOD) score of 2 and 45% power to detect a linkage with an LOD score of 3 (6).
Our genome-wide scan for BMI susceptibility genes identified two linkage peaks with an LOD 1.18: one on chromosome 1p36.13 with an LOD score of 1.38 between markers D1S2697 and D1S199 and a second, much larger peak on chromosome 10p11.22-23 with an LOD of 3.74 between markers D10S208 and D10S196 (Fig. 1). We determined by simulation that the signal on chromosome 10 was significant at the genome-wide level (P = 0.0186) and the point-wise level (P = 0.0001) using an approach that has been previously described (6). Our original genome-wide linkage analysis of type 2 diabetes susceptibility in this dataset (2) identified only marginally significant linkage to this region on 10p, with an LOD of 0.68 (asymptotic point-wise P = 0.039). Conversely, we observed no linkage of BMI (with asymptotic point-wise P < 0.05) to the region on 10q23 in which we had previously observed maximal evidence for linkage to type 2 diabetes. We can therefore be confident that the linkage evidence on 10p reflects a locus influencing BMI susceptibility directly rather than the surrogate effects of a type 2 diabetes susceptibility locus.
Our study adds substantially, in terms of significance of results and sample size, to the growing body of evidence for linkage to this region. Nine previous studies have reported evidence for linkage to BMI or related traits on the short arm or pericentric region of chromosome 10 (7,8,9,10,11,12,13,14,15). The majority of these studies, our own included, demonstrate maximum evidence for linkage in the region 55eC80 cM from the p-terminal of chromosome 10 (Fig. 1B), strongly implicating a locus in this region that influences susceptibility to body mass and obesity.
Several excellent candidate genes fall within this region on chromosome 10, including GAD2, UCN3 (urocortin III), and PPYR1/NPY4R (pancreatic polypeptide receptor 1). GAD2 encodes GAD65, which catalyzes the formation of the neurotransmitter aminobutyric acid (GABA), is synthesized by GABAergic neurones within nuclei of the hypothalamus that are known to be involved in feeding behavior (16), and has been the focus of three recent studies.
The fine-mapping linkage study of adult obesity in the French population by Boutin et al. (13) observed that the maximally linked microsatellite marker in their study, D10S197, was located within intron 7 of GAD2. Boutin et al. identified three single nucleotide polymorphisms (SNPs), a promoter variant (rs2236418, eC243 A>G), and two intronic variants (rs992990, +61450 C>A and rs928197, +83897 T>A) that are associated with obesity. In each case, the rare allele was associated with an obese phenotype. They observed marginally significant association between the variant allele of rs992990 and evidence for linkage (P = 0.02), together with functional evidence implicating the other two SNPs. The second study, also in a French population, observed that the same (G) allele of the promoter SNP (rs2236418) was associated with severe childhood obesity (17). Most recently, Swarbrick et al. (18) investigated the importance of GAD2 in obesity susceptibility in two large case/control studies (in American and Canadian populations) and in a smaller study of German pedigrees and found no statistical evidence supporting the candidacy of GAD2.
Given this prior, albeit conflicting, evidence of the importance of GAD2 in body mass determination, we examined the candidacy of GAD2 in our population. We first considered the three SNPs (rs2236418, rs992990, and rs928197) that showed some evidence of association in the previous studies. The 573 sibships were genotyped for these three SNPs using a competitive allele-specific PCR method (KASPar; Kbioscience, Hertfordshire, U.K.). Genotyping call rates exceeded 95% at each SNP, and no discordant genotypes were observed in 92 duplicate samples.
We found no evidence for population stratification with the three SNPs analyzed (all P 0.5). This fits with previous studies of these case samples when STRUCTURE has been applied to the genome-wide microsatellite data (N. Martin, L. Cardon, personal communication) and a more recent genome-wide analysis using the quantitative transmission-disequilibrium test (QTDT) (19,20) (all P 0.6, corrected for multiple testing). We therefore determined the total evidence for association of each SNP with adjusted and standardized BMI using QTDT, as well as with QPDTPHASE (21), which makes less restrictive assumptions about the data distribution (22). We observed associations of marginal significance with two of these SNPs, rs2236418 (P = 0.025) and rs928197 (P = 0.011) (Table 1). However, in both cases, the minor allele was associated with lower mean BMI, in contrast to the findings of Boutin et al. To capture additional genetic diversity within GAD2, we used HAPLOVIEW (23) and the GAD2 resequencing data from the Environmental Genome Project (www.egp.gs.washington.edu) to select an additional five SNPs for genotyping. The resulting eight-SNP panel tagged all of the haplotypes in the three largest linkage disequilibrium (LD) blocks (blocks 1, 2, and 5, spanning 73% of the total gene sequence) (Fig. 1A). The gene also contained two small blocks (blocks 3 and 4) and the intervening regions of low LD that would have required substantial additional genotyping to tag. These eight SNPs tagged the proximal promoter, the 3' untranslated region, and all exons but one, capturing 93% of all common (>10%) variation with an r2 > 0.5 and 35% of that with an r2 > 0.8 across the gene as a whole. The arrangement of the SNPs, the LD blocks that they tag, and their r2 and D' values are shown in Fig. 2. We tested the five additional SNPs individually for total association with BMI using QTDT and QPDTPHASE, having first confirmed no evidence for population stratification at these markers (all P 0.6). We also tested haplotypes of block-tagging SNPs (Table 2) with QPDTPHASE, pooling haplotypes with frequencies <1%. We found SNP rs2839669 (P = 0.036), located upstream of the associated promoter SNP rs2236418, and rs7071922 (P = 0.031), located in block 2, to be marginally associated, with the minor allele again associated with lower BMI (Table 1). Haplotypes of the block-tagging SNPs were also marginally associated with BMI (Table 2). No global association was seen with microsatellite D10S197 (genotyped in the original genome scan), located proximal to rs7908975 in block 2 (P > 0.1). No SNP showed evidence of dominance (P > 0.1).
We determined the extent to which these associated SNPs accounted for the linkage signal. The evidence for linkage to a quantitative trait is reduced when an association with an SNP is incorporated into the test for linkage. When the associated SNP is in perfect LD with the QTL, or is the QTL itself, the evidence for linkage can be eliminated completely (24). We modeled each associated SNP as a covariate, with genotype "11" coded as eC1, "12" as 0, and "22" as +1, consistent with a standard biometrical model and the gene-dosing approach used by QTDT (20). We also modeled the three sets of block-tagging SNPs as a set of covariates representing each individual SNP. We regressed each individual’s GAD2 covariate(s) from BMI using coefficients calculated with QTDT. To allow a meaningful interpretation of these GAD2 covariate analyses, we obtained the appropriate baseline multipoint LOD scores for each SNP in question, that is, the evidence for linkage alone, without accounting for the association, restricted to those family members with a successful GAD2 genotype and covariate.
For these analyses, we used the latest version of MERLIN-REGRESS, which accommodates the effects of LD between SNPs (known to inflate linkage evidence if not properly accommodated) by grouping them into clusters for the purposes of calculating allele-sharing probabilities (25). However, this method assumes no recombination within clusters and no LD between clusters. To accommodate the high multiallelic D' (i.e., nonzero LD) between the haplotype blocks in GAD2, we used two alternative clustering approaches. In the first, we treated the eight SNPs as three separate clusters (representing the three blocks they were chosen to tag), with the intervening microsatellite (D10S197, falling within block 2) considered as part of the second cluster. In the second, we grouped all SNPs (together with D10S197) into one single cluster.
We found that the baseline LOD scores differed according to the clustering approach used to model LD (Table 3, columns 2 and 4). However, the effects of the GAD2 covariates themselves did not differ appreciably between clustering methods (Table 3, columns 3 and 5). All four associated SNPs, and the haplotypes containing them, resulted in very minor LOD score reductions (LOD < 0.6; Table 3). These results suggest that not one of the SNPs tested is the QTL responsible for the linkage to BMI.
In conclusion, we have detected genome-wide significant linkage to BMI on the pericentric region of chromosome 10 in U.K. pedigrees with type 2 diabetes. This is the largest dataset to observe linkage to this region, so our study substantially adds to the body of evidence implicating loci in this region in BMI variation. We have examined the candidacy of GAD2 as the gene underlying this linkage, following two previous studies reporting association. We observed marginally significant associations to four SNPs in GAD2. Of the three SNPs in common with the Boutin study, we only observed associations with rs2236418 (P = 0.025) and rs928197 (P = 0.011), but in the opposite direction. Differences in the underlying LD structure of chromosome 10 in U.K. and French populations are highly unlikely explanations for this reversal of association. It is, in principle at least, conceivable that differences in genetic background, patterns of interacting environmental exposures, and clinical features between the two samples could explain a reversal of association direction. However, such differences would have to be extreme, which seems implausible given that these are otherwise quite similar European samples. We note that the Boutin study examined extreme obesity in a case/control setting, whereas we have quantitatively examined a continuous distribution of BMI measurements, only a very small proportion of which could be described as extremely obese. While such differences in study design might in theory explain a failure to reproduce the finding of an effect, it is difficult to see how these differences could explain a change in the direction of an effect. We conclude that our findings do not constitute replication (defined as same allele/haplotype, same phenotype) of the associations seen in previous studies.
Given the number of SNPs that we have tested, and the low correlation between them, our association findings may reflect chance. The alternative is that the associations we observe are true and reflect modest LD with other etiological variants. While it is possible that these variants lie in GAD2, as we have not exhaustively surveyed all genetic variation therein, it is more likely that they are located elsewhere in the 30-Mb region encompassing the consensus of linkage peaks.
ACKNOWLEDGMENTS
The Warren 2 family collection was funded by Diabetes UK through the Warren 2 Bequest and the genotyping through project grant RD02/0002385. S.W. is funded by a Wellcome Trust Career Development Fellowship.
We thank Drs. Stacey Cherny and Javier Gayan (Wellcome Trust Centre for Human Genetics) for useful discussions. Prof. Paola Primatesta (University College London) kindly made available population data from the Health Survey for England.
FOOTNOTES
DOI: 10.2337/db05-1674
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
REFERENCES
Perusse L, Rankinen T, Zuberi A, Chagnon YC, Weisnagel SJ, Argyropoulos G, Walts B, Snyder EE, Bouchard C: The human obesity gene map: the 2004 update. Obes Res 13:381eC490, 2005
Wiltshire S, Hattersley AT, Hitman GA, Walker M, Levy JC, Sampson M, O’Rahilly S, Frayling TM, Bell JI, Lathrop GM, Bennett A, Dhillon R, Fletcher C, Groves CJ, Jones E, Prestwich P, Simecek N, Rao PV, Wishart M, Bottazzo GF, Foxon R, Howell S, Smedley D, Cardon LR, Menzel S, McCarthy MI: A genomewide scan for loci predisposing to type 2 diabetes in a U.K. population (the Diabetes UK Warren 2 Repository): analysis of 573 pedigrees provides independent replication of a susceptibility locus on chromosome 1q. Am J Hum Genet 69:553eC569, 2001
Erens B, Primatesta P: Health Survey for England 1998. London, The Stationary Office, 2001
Sham PC, Purcell S, Cherny SS, Abecasis GR: Powerful regression-based quantitative trait linkage analysis of general pedigrees. Am J Hum Genet 71:238eC253, 2002
Schousboe K, Visscher PM, Erbas B, Kyvik KO, Hopper JL, Henriksen JE, Heitmann BL, Sorensen TI: Twin study of genetic and environmental influences on adult body size, shape, and composition. Int J Obes Relat Metab Disord 28:39eC48, 2004
Wiltshire S, Frayling TM, Groves CJ, Levy JC, Hitman GA, Sampson M, Walker M, Menzel S, Hattersley AT, Cardon LR, McCarthy MI: Evidence from a large U.K. family collection that genes influencing age of onset of type 2 diabetes map to chromosome 12p and to the MODY3/NIDDM2 locus on 12q24. Diabetes 53:855eC860, 2004
Hager J, Dina C, Francke S, Dubois S, Houari M, Vatin V, Vaillant E, Lorentz N, Basdevant A, Clement K, Guy-Grand B, Froguel P: A genome-wide scan for human obesity genes reveals a major susceptibility locus on chromosome 10. Nat Genet 20:304eC308, 1998
Hinney A, Ziegler A, Oeffner F, Wedewardt C, Vogel M, Wulftange H, Geller F, Stubing K, Siegfried W, Goldschmidt HP, Remschmidt H, Hebebrand J: Independent confirmation of a major locus for obesity on chromosome 10. J Clin Endocrinol Metab 85:2962eC2965, 2000 [erratum in J Clin Endocrinol Metab 85:3972, 2000]
Chagnon YC, Rice T, Perusse L, Borecki IB, Ho-Kim MA, Lacaille M, Pare C, Bouchard L, Gagnon J, Leon AS, Skinner JS, Wilmore JH, Rao DC, Bouchard C: HERITAGE Family Study: genomic scan for genes affecting body composition before and after training in Caucasians from HERITAGE. J Appl Physiol 90:1777eC1787, 2001
Hsueh WC, Mitchell BD, Schneider JL, St Jean PL, Pollin TI, Ehm MG, Wagner MJ, Burns DK, Sakul H, Bell CJ, Shuldiner AR: Genome-wide scan of obesity in the Old Order Amish. J Clin Endocrinol Metab 86:1199eC1205, 2001
Lindsay RS, Kobes S, Knowler WC, Bennett PH, Hanson RL: Genome-wide linkage analysis assessing parent-of-origin effects in the inheritance of type 2 diabetes and BMI in Pima Indians. Diabetes 50:2850eC2857, 2001
Price RA, Li WD, Bernstein A, Crystal A, Golding EM, Weisberg SJ, Zuckerman WA: A locus affecting obesity in human chromosome region 10p12. Diabetologia 44:363eC366, 2001
Boutin P, Dina C, Vasseur F, Dubois S, Corset L, Seron K, Bekris L, Cabellon J, Neve B, Vasseur-Delannoy V, Chikri M, Charles MA, Clement K, Lernmark A, Froguel P: GAD2 on chromosome 10p12 is a candidate gene for human obesity. PLoS Biol 1:E68, 2003
Dong C, Wang S, Li WD, Li D, Zhao H, Price RA: Interacting genetic loci on chromosomes 20 and 10 influence extreme human obesity. Am J Hum Genet 72:115eC124, 2003
Gorlova OY, Amos CI, Wang NW, Shete S, Turner ST, Boerwinkle E: Genetic linkage and imprinting effects on body mass index in children and young adults. Eur J Hum Genet 11:425eC432, 2003
Kalra SP, Dube MG, Pu S, Xu B, Horvath TL, Kalra PS.: Interacting appetite-regulating pathways in the hypothalamic regulation of body weight. Endocr Rev 20:68eC100, 1999
Meyre D, Boutin P, Tounian A, Deweirder M, Aout M, Jouret B, Heude B, Weill J, Tauber M, Tounian P, Froguel P: Is glutamate decarboxylase 2 (GAD2) a genetic link between low birth weight and subsequent development of obesity in children J Clin Endocrinol Metab 90:2384eC2390, 2005
Swarbrick MM, Waldenmaier B, Pennacchio LA, Lind DL, Cavazos MM, Geller F, Merriman R, Ustaszewska A, Malloy M, Scherag A, Hsueh WC, Rief W, Mauvais-Jarvis F, Pullinger CR, Kane JP, Dent R, McPherson R, Kwok PY, Hinney A, Hebebrand J, Vaisse C: Lack of support for the association between GAD2 polymorphisms and sever human obesity. PLOS Biol 3:e315, 2005
Abecasis GR, Cardon LR, Cookson WO: A general test of association for quantitative traits in nuclear families. Am J Hum Genet 66:279eC292, 2000
Abecasis GR, Cookson WO, Cardon LR: Pedigree tests of transmission disequilibrium. Eur J Hum Genet 8:545eC551, 2000
Dudbridge F: Pedigree disequilibrium tests for multilocus haplotypes. Genet Epidemiol 25:115eC121, 2003
Monks SA, Kaplan NL: Removing the sampling restrictions from family-based tests of association for a quantitative trait locus. Am J Hum Genet 66:576eC592, 2000
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualisation of LD and haplotype maps. Bioinformatics 21:263eC265, 2005
Fulker DW, Cherny SS, Sham PC, Hewitt JK: Combined linkage and association sib-pair analysis for quantitative traits. Am J Hum Genet 64:259eC267, 1999
Abecasis GR, Wigginton JE: Handling marker-marker linkage disequilibrium: pedigree analysis with clustered markers. Am J Hum Genet 77:754eC767, 2005(Christopher J. Groves, El)