The growth hormone secretagogue receptor (GHSR) (ghrelin receptor) plays an important role in the regulation of food intake and energy homeostasis. The GHSR gene lies on human chromosome 3q26 within a quantitative trait locus strongly linked to multiple phenotypes related to obesity and the metabolic syndrome. Because the biological function and location of the GHSR gene make it an excellent candidate gene, we tested the relation between common single nucleotide polymorphisms (SNPs) in the GHSR gene and human obesity. We performed a comprehensive analysis of SNPs, linkage disequilibrium (LD), and haplotype structure across the entire GHSR gene region (99.3 kb) in 178 pedigrees with multiple obese members (DNA of 1,095 Caucasians) and in an independent sample of the general population (MONICA Augsburg left ventricular hypertrophy substudy; DNA of 1,418 Caucasians). The LD analysis revealed a disequilibrium block consisting of five SNPs, consistent in both study cohorts. We found linkage among all five SNPs, their haplotypes, and BMI. Further, we found suggestive evidence for transmission disequilibrium for the minor SNP alleles (P < 0.05) and the two most common haplotypes with the obesity affection status (“susceptible” P = 0.025, “nonsusceptible” P = 0.045) in the family cohort using the family-based association test program. Replication of these findings in the general population resulted in stronger evidence for an association of the SNPs (best P = 0.00001) and haplotypes with the disease (“susceptible” P = 0.002, “nonsusceptible” P = 0.002). To our knowledge, these data are the first to demonstrate linkage and association of SNPs and haplotypes within the GHSR gene region and human obesity. This linkage, together with significant transmission disequilibrium in families and replication of this association in an independent population, provides evidence that common SNPs and haplotypes within the GHSR region are involved in the pathogenesis of human obesity.
Obesity is a common multifactorial disorder of considerable heterogeneity and, as a pivotal component of the metabolic syndrome, a major risk factor for type 2 diabetes, hypertension, and coronary heart disease as well as premature cardiovascular morbidity and death (1). Together with its associated pathologic features, obesity is among the major causes of illness and death worldwide, as its prevalence continues to rise dramatically (2).
The etiology of obesity is complex, determined by the interplay of genetic and environmental factors. Epidemiological studies have demonstrated a substantial heritable component to the risk for obesity; specifically, 50–70% of the variation in BMI may be attributable to genetic factors (3).
To unravel the genetic etiology of obesity, we previously performed a genomewide linkage scan on a large cohort of Caucasian families from which we localized on chromosome 3q26-q29 a major quantitative trait locus (QTL) strongly linked to six phenotypes of obesity and the metabolic syndrome (4). This QTL has been replicated in several studies and represents one of the most stable findings in complex human genetics (5–9).
A comprehensive review of the available genomic information in the QTL region reveals a positional candidate gene of ∼4.3 kb in length encoding the growth hormone secretagogue, or ghrelin, receptor (GHSR). GHSR is known to be involved in growth hormone secretion (10,11). Its major physiological role, however, appears to be in regulating food intake and energy homeostasis by partaking in neuronal mechanisms involving neuropeptide Y and agouti-related protein (12–15). The endogenous GHSR ligand, ghrelin, plays a key role as the major orexigenic hormone. It is secreted in the gastrointestinal tract and is carried to the hypothalamic areas that govern food intake, thereby counterbalancing the effects of a multitude of anorectic hormones, such as leptin, insulin, and PYY3–36 (16).
The importance of ghrelin in the central regulation of feeding has been demonstrated in animals and humans (17,18). Ghrelin administration increases appetite and food intake in normal subjects and patients with decreased appetite, such as those suffering from cancer cachexia (17). It reduces insulin secretion and enhances energy intake by ∼30% (19). Moreover, given that plasma ghrelin levels have been shown to be lower in obese subjects (20,21), recent evidence suggests that obesity is associated with an impairment of the entire ghrelin system (22).
Because its biological function and location make the GHSR gene a strong candidate gene for obesity, we carried out a family-based linkage disequilibrium (LD) study in 178 pedigrees as well as in an independent case-control study from the general population. We systematically explored the LD and haplotype structure of the genomic region encompassing the GHSR gene and comprehensively assessed the role of common sequence variants and haplotypes in obesity. We report evidence for linkage and association of five single nucleotide polymorphisms (SNPs) and the two most common five-marker haplotypes with obesity in our family cohort. In addition, we describe the association of the same SNPs and haplotypes with obesity and, more striking, with the quantitative phenotype BMI in the general population. Thus, the replication of our findings, together with the location and biological function of the GHSR gene, indicate that this gene region is involved in the pathogenesis of the complex disease of obesity.
RESEARCH DESIGN AND METHODS
The family subjects of the study, a large group of Caucasian families consisting of obese and nonobese members residing in the midwestern U.S., were ascertained through the TOPS (Take Off Pounds Sensibly) Club membership as part of the Metabolic Risk and Complications of Obesity Genes project at the Medical College of Wisconsin. The ascertainment strategies and exclusion criteria have been previously published (4). Informed consent was obtained from all participants. Data for the analyses presented here were based on results obtained from 1,302 phenotyped individuals distributed among 178 families (average 7.3 members per family; 441 founders; both parents available for genotyping in 148 families; average generations 2.3). These families have been identified as primary contributors to the QTL on chromosome 3. In all, 1,095 DNA samples were available for genotyping in the present study, including samples of 307 (28%) men and 788 (72%) women. Research protocols were approved by the institutional review board of the Medical College of Wisconsin.
For the general population arm of the study, we used data from subjects in the Monitoring Trends and Determinants in Cardiovascular Disease (MONICA) Augsburg left ventricular hypertrophy (LVH) substudy, as part of the Third MONICA Augsburg survey, which now is continued in the framework of KORA (Cooperative Health Research in the Augsburg Area). The study population of the LVH substudy was sampled from the general population of the city of Augsburg, Germany, in 1994/1995, which originated from a sex- and age-stratified cluster sample of all German residents of the city of Augsburg. The Augsburg project was part of the international collaborative World Health Organization MONICA study (23). The study design, sampling frame, and data collection methods have been described in detail elsewhere (23,24). All the participants gave written informed consent. The LVH substudy represents individuals aged 25–74 years, with ∼300 subjects for each 10-year increment (n = 1,674) (24). Of these, 1,418 DNA samples were available for genotyping in the present study, including 724 men (51%) and 694 women (49%). The study was approved by the local ethics committee.
BMI was calculated as weight (kg) divided by height (m) squared. In both cohorts, obesity was defined by a BMI >32 kg/m2. Subjects were classified as “unaffected” if they presented with a BMI <28 kg/m2. These cutoff values were chosen to ensure clear phenotypes and avoid misclassification regarding affection status. The obesity affection status of subjects with a BMI of 28–32 kg/m2 was treated as “unknown.”
SNPs and genotyping methods
SNPs.
To obtain complete coverage of the GHSR gene, 10 SNPs covering the GHSR gene and its flanking regions were selected from the SNP public databases (dbSNP; available at http://www.ncbi.nlm.nih.gov/SNP) (Fig. 1). We preferred validated SNPs with a minor allele frequency of >5%. Priority was given to SNPs submitted multiple times then to SNPs discovered by The SNP Consortium (25,26). Regarding the intergenic regions, we prioritized SNPs located in highly conserved noncoding regions. Of the 10 selected SNPs, 1 was located in exon 1, 1 was in the intron, 3 were within 41.5 kb past the 3′ end of the gene, and 5 covered a region of 53.5 kb past the 5′ end of the gene. The coding SNP (rs572169) led to a synonymous amino acid substitution. The eight SNPs located beyond the boundaries of the gene were picked to determine the extent of LD and explore the impact of sequence variations in noncoding and intergenic regions on the disease. In total, a 99.3-kb region was covered with SNPs, with an average resolution of one SNP per 10 kb.
Genotyping.
In families, SNP genotyping was performed using the Biplex Invader assay (Third Wave Technologies), as previously described (27). In the general population, SNPs were genotyped using the 5′-exonuclease activity (TaqMan) assay on a HT7900 (Applied Biosystems, Darmstadt, Germany). H2O controls, CEPH DNA samples of known genotype, as well as blind duplicate samples were included in each round of amplification to check for consistency and to ensure intra- and interplate genotype quality control.
In neither cohort were genotyping discrepancies detected between the repeated samples. The overall misgenotyping rate of <0.5% was due to insufficient PCR amplification.
Statistical analysis.
For each of the 10 SNPs, we tested whether the observed allele frequencies departed from the Hardy-Weinberg proportion. No deviations from the expected genotype proportions were detected for any of the SNPs used in the analyses. We also assessed LD between all pairs of SNPs, applying the standard definition of r2 (28,29).
Families.
All genotype data were analyzed using PEDCHECK to check for Mendelian inconsistencies and genotype incompatibilities (30). To test for linkage between these SNPs and the quantitative trait BMI, we used the quantitative transmission disequilibrium test (QTDT) program, which is based on the standard variance component methods and identity by descent among relatives (31). To analyze transmission disequilibrium between the discrete trait obesity and the SNPs, we selected trios with both available parents and one randomly chosen “obesity-affected” offspring and computed the conventional transmission disequilibrium test (TDT) statistic (32). In addition, we applied the software program family-based association test (FBAT) version 1.5 to handle the different types of family structures and use all the information in each family (33). Here, we tested family-based association under the null hypothesis of linkage but no association using the −e flag. This option computes the test statistic through use of the empirical variance (34), which is needed because the markers are in an area of known linkage and the sample contains multiple nuclear families in some pedigrees and multiple affected individuals in these nuclear families.
Structures of haplotypes were analyzed from parental genotypes based on LD pattern using the expectation-maximization algorithm (35). A haplotype block was defined as a region in which all pairwise r2 values were >0.45 (29). Haplotypes with a frequency >5% were tested for association with obesity using the haplotype FBAT program, which can be used in candidate gene studies with tightly linked markers (36). P values were corrected for multiple testing by a method that specifically considers SNPs being in LD with each other and is based on the spectral decomposition (available at http://genepi.qimr.edu.au/general/daleN/SNPSpD/) of matrices of pairwise LD between SNPs (37). Moreover, we tested whether the linkage at this region could be explained by the haplotypes associated with BMI using the QTDT software by modeling linkage and association simultaneously as suggested by Fulker et al. (38). For this purpose, individual haplotypes were reconstructed using HAPLORE (39).
General population.
Haplotype frequencies were estimated using the expectation-maximization algorithm (35). By performing an automated, randomized selection of control subjects, all available obese subjects (n = 130 case subjects) were carefully matched by age (±5 years) and sex, whereby one case subject was matched with up to three control subjects (n = 364 control subjects). Frequencies of genotypes in the total study sample as well as in case and control subjects were compared by the Armitage test for trend, and odds ratios with their 95% CIs were reported (40). In addition to single-locus analysis, statistically inferred five-marker haplotypes were tested for association with the discrete trait obesity as well as with the continuous trait BMI using the haplotype trend regression method (41). In this analysis, all individuals were considered, including those with a BMI of 28–32 kg/m2. Permutation tests (50,000 permutations) were used to test for empirical significance. We also divided the MONICA LVH population into quartiles of BMI distribution, calculated the haplotype frequencies, and tested for significance using Cochran’s test for trend.
RESULTS
Phenotypic characteristics.
The phenotypic characteristics of TOPS families and the MONICA LVH population are presented in Table 1. In TOPS families, the mean BMI and prevalence of obesity in men and women were markedly different (P < 0.001). In the MONICA LVH population, the prevalence of obesity was significantly higher in women than in men (P < 0.001).
Linkage disequilibrium evaluation and haplotype structure in families and the general population.
Figure 1 depicts the gene structure and all SNPs used in this study, including their position and general characteristics based on the July 2003 release of the Golden Path Genome Browser (available at http://www.genome.ucsc.edu). The pairwise LD block structure defined by the 10 SNPs covering a 99.3-kb region in TOPS families and the MONICA LVH population is shown in Fig. 2. A region of strong LD (r2 >0.75) was detected between five SNPs (rs509035, rs572169, rs519384, rs512692, and rs863441), and the LD pattern was comparable in both study populations. The pairwise LD between these five SNPs spanned a 11.63-kb region and encompassed most of the intron, exon 1, and 5′ adjacent region of the GHSR gene.
According to this high-LD block, five-marker haplotypes were constructed. Only 3 of 25 possible haplotypes were estimated to have frequencies >0.01 in obese and/or nonobese subjects. Less frequently occurring haplotypes were not shown, owing to concern over the accuracy of low-frequency alleles in the expectation-maximization algorithm. The two most frequently occurring haplotypes, haplotypes 1 and 2, comprised 94% of total chromosomes in subjects of both study cohorts (Tables 2 and 3). For further linkage and association analyses, we focused on those markers that contributed to this haplotype block. We tested the hypothesis of a relation between these five SNPs and/or five-marker haplotypes and obesity phenotypes. To ensure that the SNPs, which were not included in the high-LD block, were not associated with obesity, we performed association analysis for these SNPs. None of the SNPs showed any evidence for association with the obesity affection status or with BMI (data not shown).
Family data: genetic linkage and transmission disequilibrium of SNPs and haplotypes in the GHSR region.
To test for linkage in families, we used the variance component methodology. Evidence for linkage with the quantitative phenotype BMI was detected for all five SNPs forming the high-LD block depicted in Fig. 2 (P < 0.05; data not shown). Thus, our SNP genotype data confirmed linkage to the previously shown QTL on chromosome 3q (4).
To test for transmission disequilibrium in families, we applied both the conventional TDT statistic (in the 148 trios with one randomly selected affected offspring) and the FBAT statistic (considering all family members) for each of the SNPs contributing to the haplotype. The TDT analysis revealed increased transmission for the minor alleles of the five SNPs to obese offspring (Table 2). A slightly stronger pattern of association of the single SNPs with the obesity-affection status was observed when the FBAT statistic was used (P < 0.05 for all five SNPs) (Table 2).
In addition, we observed transmission disequilibrium for the two most frequent haplotypes, one consisting of the five major alleles (haplotype 1) and the other consisting of the five minor alleles (haplotype 2) (Table 2). Corresponding to the “susceptible” haplotype, haplotype 2 had a greater number of transmissions to affected offspring (P = 0.025). In contrast, the transmission rate of haplotype 1 was significantly reduced in these offspring, suggesting that this haplotype is “nonsusceptible” or resistant to obesity (P = 0.045).
After reconstructing the individual haplotypes, we found suggestive evidence for linkage with the quantitative trait BMI (P = 0.06). Modeling linkage and association simultaneously resulted in no residual evidence for linkage at this haplotype marker (P = 0.57). This indicated that the evidence of linkage at this site was accounted for by association; that is, the haplotype marker contained the disease mutation itself or was in strong LD with it.
General population data: association of SNPs and haplotypes with obesity and BMI.
Aiming to replicate the results obtained in the family cohort, we then performed an association analysis in an independent sample of the general population (MONICA Augsburg LVH substudy). Results of the association of the five SNPs are summarized in Table 4 for the entire study sample as well as for matched case and control subjects. Odds ratios were calculated for the comparison of allele frequencies and the “homozygous trait” and “allele positivity” comparisons. Overall, the five SNPs consistently showed nominally significant association with obesity in all three comparisons, in both the entire study sample and the matched case and control subjects (entire study sample, best P = 0.0000; matched sample, best P = 0.0007 for rs863441). When the result was corrected for multiple testing for SNPs in LD, most P values remained significant. In the entire study sample, the increased risk presented by the presence of the minor allele of these SNPs ranged between 41 (P = 0.014) and 56% (P = 0.001).
To further test for association, we next undertook full haplotype analysis. In agreement with our family data, haplotype 1 was more frequently present in nonobese individuals (“nonsusceptible,” P = 0.002 after 50,000 permutations) and haplotype 2 was more frequently found in obese individuals (“susceptible,” P = 0.001 after 50,000 permutations) (Table 3). Moreover, we observed a significant relation between the number of copies of haplotypes 1 and 2 with the qualitative trait obesity in matched case and control subjects and that in the entire study sample. Individuals homozygous for the “susceptible” haplotype 2 or lacking the “nonsusceptible” haplotype 1 presented more often with obesity than individuals with one or no copy of the respective haplotype (haplotype 2: P = 0.003 in the matched sample, P = 0.006 in the entire sample; haplotype 1: P = 0.003 in the matched sample, P = 0.003 in the entire sample after 50,000 permutations) (Fig. 3). When tested for association with BMI, individuals carrying two copies of the “susceptible” haplotype 2 or no copy of the “nonsusceptible” haplotype 1 analogously presented with higher BMI than individuals with one or no copy of the respective haplotype (haplotype 2: P = 0.009 in the matched sample, P = 0.007 in the entire sample; haplotype 1: P = 0.005 in the matched sample, P = 0.006 in the entire sample after 50,000 permutations) (Fig. 3). We also divided the MONICA LVH population into quartiles of BMI distribution and observed a significantly increasing frequency of the “susceptible” haplotype 2, from 23% in the lowest quartile to 30% in the highest quartile (P <0.004 for trend; data not shown).
DISCUSSION
This study investigated the relation between common sequence variants and haplotypes covering the GHSR gene region with obesity phenotypes in families and in an independent sample of the general population. We were led to pursue the GHSR because 1) the major physiological role of the GHSR appears to be in the central regulation of food intake and body weight; 2) the GHSR gene is located within the QTL on chromosome 3q26-q29, which has been linked to six phenotypes of the metabolic syndrome, including BMI, in the same family cohort; and 3) only one association study investigating the relation between GHSR gene variants and obesity has thus far been performed (42). Because of its location and biological function, the GHSR gene represents an excellent candidate gene that may contribute to the pathophysiology and risk of obesity.
The present study offers the first comprehensive analysis of LD, genetic variants, and haplotype structure across the entire GHSR gene region in two independent cohorts: families and the general population. The initial LD analysis in the 99.3-kb region revealed an LD block consisting of five SNPs in the GHSR gene region, which compared very well in both study cohorts. Subsequently, we focused on these five SNPs and the five-SNP haplotypes. We report linkage between all five SNPs and BMI and, furthermore, provide weak yet suggestive evidence for transmission disequilibrium for the minor alleles of the SNPs as well as for the two most common five-SNP haplotypes with the obesity affection status. The replication of these findings in an independent sample of the general population further supports an association of GHSR gene variants with human obesity. Moreover, we report that our haplotypes or one in LD with them could account in part for the observed linkage signal. Thus, our data implicate common haplotypes in this gene region in the pathogenesis of human obesity.
We initially determined the extent of the high-LD region by covering the entire gene region, including the surrounding genomic regions close to the neighboring genes, with SNPs. The identified high-LD region encompasses part of the intron, exon 1, and the 5′ adjacent region extending 8.8 kb past the 5′ end of the gene, but not encompassing the flanking genes. Therefore, it is unlikely that the association between variants of the GHSR gene and obesity is seen because of LD with the proper causal mutation in one of the neighboring genes. In addition, we tested the SNPs that were not included in the high-LD block for association (data not shown) and found that none of the SNPs showed evidence for association with the obesity affection status or BMI. This observation supports the hypothesis that genetic variations within the LD block encompassing the GHSR gene, and not within neighboring genes, are related to our obesity phenotypes.
We focused only on common sequence variants, as it is more likely that these variants play a role in the general population. Furthermore, we included SNPs located in noncoding and intergenic regions rather than exclusively focusing on the coding region. This strategy was driven by the hypothesis that variations underlying complex diseases would not be limited to the structure of the encoded protein; that is, they could be due to variants leading to altered gene expression (43). Gene regulation is the result of the combinatorial action of multiple transcription factors binding at multiple sites in and near a gene and therefore can be affected by multiple SNPs. In fact, it has been recently demonstrated that gene regulatory elements reside in noncoding and intergenic regions (44,45). These enhancers are able to modulate gene expression over long distances, turning intergenic regions into reservoirs for sequence elements containing important functions (46). Little is known about the impact of sequence variations in these regions. In our study, SNPs located in the intergenic region past the 5′ end of the GHSR gene showed stronger association than the SNPs located in the coding or intronic region of the gene. These data suggest that the promoter, regulatory elements or transcriptional initiation could be involved. The exclusive screening of the coding region could be one possible explanation for the lack of association in the study by Wang et al. (42), who were not able to find a relation between two coding SNPs of the GHSR gene (one of them, rs 572169, was investigated in the present study) and obesity in children and adolescents. The extremely young age of Wang et al.’s population and the missing haplotype analysis offer further explanations for those researchers’ negative results.
The fact that we showed a consistent association between the five SNPs and the two most frequent five-marker haplotypes with the obesity affection status and quantitative BMI in families, the entire MONICA LVH substudy sample, and matched case and control subjects argues against any reasonable likelihood that the findings are artifacts of population stratification or multiple testing. There were no differences with respect to the SNP or haplotype allele frequencies between the two study cohorts and almost all significant P values remained after correction for multiple testing using a method for SNPs in LD or running 50,000 permutations. Because our region showed high LD, we did not use the conventional Bonferroni correction, which assumes independence of markers and therefore would make the correction overly conservative (37).
Although the associated minor allele haplotype (haplotype 2) seems to confer susceptibility to obesity, the major allele haplotype (haplotype 1) acts in a reverse fashion by lowering the risk of obesity. Either effect is present independent of carrying one or two copies. However, the effect is strongest in those presenting with two copies of the respective haplotype and decreases with the number of copies.
It is interesting that the ghrelin receptor, encoded by the GHSR gene, along with its endogenous ligand ghrelin provides the only hormonal, appetite-stimulatory input that counterbalances a large number of inhibitory signals that are mediated by leptin, insulin, and PYY3–36 (14,47,48). GHSR is expressed in neuropeptide Y− and agouti-related protein−containing neurons in the hypothalamus that respond to ghrelin by increasing their firing rate (16). Recently, it was shown that during fasting, GHSR expression is increased eightfold, which would be expected to result in an increase in receptor signaling and thereby an increase in appetite (49). Accordingly, genetic variations in the ghrelin receptor gene, and thereby altered expression of the receptor protein, should, in turn, result in altered signaling and consequently altered regulation of appetite. Thus, increased ghrelin receptor expression should be expected to be associated with obesity. Furthermore, it was recently shown that the ghrelin receptor exhibits a high constitutive activity signal of ∼50% efficacy between meals and thus provides a set point for food intake between meals (15). An increase in this constitutive activity based on genetic variation, such as the “susceptible” haplotype 2, could result in decreased sensitivity to the multiple inhibitory signals and consequently promote snack-eating behavior between meals. On the other hand, drugs blocking this constitutive activity of the ghrelin receptor might reduce the craving for desserts and intermeal snacks by increasing sensitivity to inhibitory signals (16).
Thus, it seems convincing that genetic variations in the ghrelin receptor gene may change either ghrelin receptor expression or receptor properties and thereby have an effect on appetite regulation by altered signaling, altered response to ghrelin, or an impaired capability to counterbalance inhibitory signals. A greater susceptibility to obesity could be the consequence. Additional functional studies are needed to clarify whether individuals carrying the “susceptible” haplotype present with altered receptor activity by, for example, measuring the inositol (1,4,5)-triphosphate turnover or determining the activation of cAMP responsive element−mediated gene transcription (15,50). Along a similar line, it should be investigated whether those individuals carrying the “nonsusceptible” haplotype exhibit more favorable receptor properties.
In summary, our work offers a comprehensive analysis of the LD structure, common genetic variants, and haplotypes within the GHSR gene region. To our knowledge, these data are the first to demonstrate linkage and association of genetic variants within the GHSR gene region and human obesity. The findings of linkage together with both transmission disequilibrium in families and the replication of this association in an independent population provide evidence that variants within the GHSR gene region might influence susceptibility to obesity and be involved in the pathogenesis of this complex disease. As we focus on common SNPs and haplotypes, the conclusions should be of great importance for a significant proportion of the population. Moreover, they provide background for the development of efficient antiobesity drugs, especially for individuals with the “susceptible” haplotype.
. | Men . | Women . |
---|---|---|
TOPS | ||
n | 307 | 788 |
Age (years) | 50.9 ± 17.6 | 46.0 ± 14.4 |
Obesity affection status (%) | 24.7 | 50.1 |
BMI (kg/m2) | ||
Total sample | 29.1 ± 5.6 | 33.1 ± 8.4 |
Affected | 37.0 ± 4.5 (75) | 39.8 ± 6.2 (393) |
Unaffected | 24.5 ± 2.4 (141) | 24.3 ± 2.4 (244) |
Body weight (kg) | 93.2 ± 21.1 | 89.7 ± 24.2 |
Waist-to-hip ratio | 0.95 ± 0.09 | 0.85 ± 0.10 |
Hypertension (%) | 41.8 | 36.0 |
Type 2 diabetes (%) | 7.9 | 7.8 |
Current smoker (%) | 37.2 | 31.0 |
MONICA | ||
n | 724 | 694 |
Age (years) | 52.5 ± 13.9 | 51.6 ± 13.6 |
Obesity affection status (%) | 9.9 | 16.1 |
BMI (kg/m2) | ||
Total sample | 27.0 ± 3.5 | 26.4 ± 4.7 |
Affected | 34.6 ± 2.5 (50) | 35.0 ± 2.7 (89) |
Unaffected | 24.9 ± 2.1 (457) | 23.8 ± 2.5 (463) |
Body weight (kg) | 81.9 ± 11.5 | 68.8 ± 11.9 |
Waist-to-hip ratio | 0.92 ± 0.06 | 0.80 ± 0.06 |
Hypertension (%) | 41.9 | 33.4 |
Type 2 diabetes (%) | 5.2 | 3.0 |
Current smoker (%) | 30.2 | 23.6 |
. | Men . | Women . |
---|---|---|
TOPS | ||
n | 307 | 788 |
Age (years) | 50.9 ± 17.6 | 46.0 ± 14.4 |
Obesity affection status (%) | 24.7 | 50.1 |
BMI (kg/m2) | ||
Total sample | 29.1 ± 5.6 | 33.1 ± 8.4 |
Affected | 37.0 ± 4.5 (75) | 39.8 ± 6.2 (393) |
Unaffected | 24.5 ± 2.4 (141) | 24.3 ± 2.4 (244) |
Body weight (kg) | 93.2 ± 21.1 | 89.7 ± 24.2 |
Waist-to-hip ratio | 0.95 ± 0.09 | 0.85 ± 0.10 |
Hypertension (%) | 41.8 | 36.0 |
Type 2 diabetes (%) | 7.9 | 7.8 |
Current smoker (%) | 37.2 | 31.0 |
MONICA | ||
n | 724 | 694 |
Age (years) | 52.5 ± 13.9 | 51.6 ± 13.6 |
Obesity affection status (%) | 9.9 | 16.1 |
BMI (kg/m2) | ||
Total sample | 27.0 ± 3.5 | 26.4 ± 4.7 |
Affected | 34.6 ± 2.5 (50) | 35.0 ± 2.7 (89) |
Unaffected | 24.9 ± 2.1 (457) | 23.8 ± 2.5 (463) |
Body weight (kg) | 81.9 ± 11.5 | 68.8 ± 11.9 |
Waist-to-hip ratio | 0.92 ± 0.06 | 0.80 ± 0.06 |
Hypertension (%) | 41.9 | 33.4 |
Type 2 diabetes (%) | 5.2 | 3.0 |
Current smoker (%) | 30.2 | 23.6 |
Data are means ± SD, with (n) where appropriate, unless otherwise indicated. The obesity affection status is defined as a BMI >32 kg/m2. Hypertension is defined as systolic blood pressure >140 mmHg, diastolic blood pressure >90 mmHg, or a history of hypertension. Because of nonindependent observations in TOPS families, characteristics are descriptive. Subjects classified as “unknown” are not shown.
SNP/haplotype . | Frequency . | Allele . | Trio TDT (n probands with allele) . | . | . | FBAT P . | ||
---|---|---|---|---|---|---|---|---|
. | . | . | Transmitted . | Not transmitted . | P . | . | ||
rs509035 | 0.71 | G | 24 | 37 | 0.096 (0.137) | 0.021 (0.030) | ||
0.29 | A | 37 | 24 | |||||
rs572169 | 0.71 | G | 23 | 38 | 0.055 (0.079) | 0.032 (0.046) | ||
0.29 | A | 38 | 23 | |||||
rs519384 | 0.73 | T | 23 | 37 | 0.070 (0.100) | 0.010 (0.014) | ||
0.27 | A | 37 | 23 | |||||
rs512692 | 0.74 | A | 21 | 38 | 0.027 (0.039) | 0.033 (0.047) | ||
0.26 | T | 38 | 21 | |||||
rs863441 | 0.73 | G | 22 | 37 | 0.051 (0.073) | 0.012 (0.017) | ||
0.27 | C | 37 | 22 | — | — | |||
Haplotype 1 | 0.69 | G-G-T-A-G | — | — | — | 0.025 (NA) | ||
Haplotype 2 | 0.25 | A-A-A-T-C | — | — | — | 0.045 (NA) | ||
Haplotype 3 | 0.03 | A-A-T-A-G | — | — | — | 0.905 (NA) |
SNP/haplotype . | Frequency . | Allele . | Trio TDT (n probands with allele) . | . | . | FBAT P . | ||
---|---|---|---|---|---|---|---|---|
. | . | . | Transmitted . | Not transmitted . | P . | . | ||
rs509035 | 0.71 | G | 24 | 37 | 0.096 (0.137) | 0.021 (0.030) | ||
0.29 | A | 37 | 24 | |||||
rs572169 | 0.71 | G | 23 | 38 | 0.055 (0.079) | 0.032 (0.046) | ||
0.29 | A | 38 | 23 | |||||
rs519384 | 0.73 | T | 23 | 37 | 0.070 (0.100) | 0.010 (0.014) | ||
0.27 | A | 37 | 23 | |||||
rs512692 | 0.74 | A | 21 | 38 | 0.027 (0.039) | 0.033 (0.047) | ||
0.26 | T | 38 | 21 | |||||
rs863441 | 0.73 | G | 22 | 37 | 0.051 (0.073) | 0.012 (0.017) | ||
0.27 | C | 37 | 22 | — | — | |||
Haplotype 1 | 0.69 | G-G-T-A-G | — | — | — | 0.025 (NA) | ||
Haplotype 2 | 0.25 | A-A-A-T-C | — | — | — | 0.045 (NA) | ||
Haplotype 3 | 0.03 | A-A-T-A-G | — | — | — | 0.905 (NA) |
SNPs are shown as rs numbers from the dbSNP database. Haplotypes are derived from the five SNPs contributing to the high-LD block (r2 >0.75). P values were corrected for multiple testing (shown in parentheses). P values for trio TDT are based on the TDT statistic; P values for FBAT are based on the FBAT (−e option) statistic. NA, correction not applied.
Haplotypes . | . | . | . | . | . | Frequency . | . | Asymptotic P . | Empirical P . | ||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
ID . | rs509035 . | rs572169 . | rs519384 . | rs512692 . | rs863441 . | Obese . | Non-obese . | . | . | ||||||
1 | G | G | T | A | G | 0.61 | 0.69 | 0.002 | 0.002 | ||||||
2 | A | A | A | T | C | 0.33 | 0.25 | 0.002 | 0.001 | ||||||
3 | A | A | T | A | G | 0.04 | 0.04 | 0.474 | 0.452 |
Haplotypes . | . | . | . | . | . | Frequency . | . | Asymptotic P . | Empirical P . | ||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
ID . | rs509035 . | rs572169 . | rs519384 . | rs512692 . | rs863441 . | Obese . | Non-obese . | . | . | ||||||
1 | G | G | T | A | G | 0.61 | 0.69 | 0.002 | 0.002 | ||||||
2 | A | A | A | T | C | 0.33 | 0.25 | 0.002 | 0.001 | ||||||
3 | A | A | T | A | G | 0.04 | 0.04 | 0.474 | 0.452 |
Haplotypes are derived from the five SNPs contributing to the high-LD block (r2 > 0.75). Asymptotic P values are based on haplotype trend regression analysis; empirical P values are based on 50,000 permutations. The permutation test was implemented in the haplotype trend regression analysis program.
SNP . | Case subjects . | . | . | . | Control subjects . | . | . | . | Allele 2 vs. 1 . | . | . | Genotype 22 vs. 11 . | . | . | Genotypes 22 + 12 vs. 11 . | . | . | ||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
. | 11 . | 12 . | 22 . | MAF . | 11 . | 12 . | 22 . | MAF . | OR (95% CI) . | P . | P(corr) . | OR (95% CI) . | P . | P(corr) . | OR (95% CI) . | P . | P(corr) . | ||||||||||||
rs509035 | |||||||||||||||||||||||||||||
Matched | 47 | 59 | 17 | 0.38 | 173 | 166 | 21 | 0.29 | 1.50 (1.10–2.03) | 0.009 | 0.013 | 2.98 (1.46–6.10) | 0.002 | 0.003 | 1.50 (0.99–2.27) | 0.058 | 0.085 | ||||||||||||
Total study | 49 | 65 | 17 | 0.38 | 424 | 375 | 69 | 0.30 | 1.45 (1.11–1.90) | 0.007 | 0.010 | 2.13 (1.16–3.91) | 0.013 | 0.019 | 1.60 (1.10–2.33) | 0.015 | 0.022 | ||||||||||||
rs572169 | |||||||||||||||||||||||||||||
Matched | 50 | 62 | 17 | 0.37 | 182 | 157 | 25 | 0.28 | 1.49 (1.11–2.01) | 0.009 | 0.013 | 2.48 (1.24–4.94) | 0.009 | 0.013 | 1.58 (1.05–2.38) | 0.028 | 0.041 | ||||||||||||
Total study | 52 | 67 | 17 | 0.37 | 451 | 371 | 80 | 0.29 | 1.42 (1.09–1.85) | 0.010 | 0.015 | 1.84 (1.02–3.35) | 0.042 | 0.062 | 1.62 (1.12–2.34) | 0.011 | 0.016 | ||||||||||||
rs519384 | |||||||||||||||||||||||||||||
Matched | 56 | 60 | 15 | 0.34 | 207 | 155 | 16 | 0.25 | 1.59 (1.18–2.16) | 0.003 | 0.004 | 3.47 (1.62–7.44) | 0.0009 | 0.001 | 1.62 (1.09–2.42) | 0.018 | 0.026 | ||||||||||||
Total study | 59 | 65 | 15 | 0.34 | 501 | 361 | 58 | 0.26 | 1.48 (1.13–1.94) | 0.004 | 0.006 | 2.20 (1.17–4.12) | 0.012 | 0.018 | 1.62 (1.13–2.33) | 0.008 | 0.012 | ||||||||||||
rs512692 | |||||||||||||||||||||||||||||
Matched | 57 | 55 | 15 | 0.33 | 200 | 154 | 15 | 0.25 | 1.51 (1.11–2.06) | 0.008 | 0.012 | 3.51 (1.62–7.61) | 0.0009 | 0.001 | 1.45 (0.97–2.18) | 0.070 | 0.103 | ||||||||||||
Total study | 60 | 60 | 15 | 0.33 | 482 | 356 | 56 | 0.26 | 1.41 (1.07–1.86) | 0.014 | 0.021 | 2.15 (1.15–4.04) | 0.015 | 0.022 | 1.46 (1.02–2.11) | 0.040 | 0.059 | ||||||||||||
rs863441 | |||||||||||||||||||||||||||||
Matched | 56 | 56 | 18 | 0.35 | 204 | 149 | 13 | 0.24 | 1.74 (1.28–2.37) | 0.0003 | 0.0004 | 5.04 (2.33–10.9) | 0.00001 | 0.00002 | 1.66 (1.11–2.49) | 0.013 | 0.019 | ||||||||||||
Total study | 59 | 61 | 18 | 0.35 | 493 | 355 | 55 | 0.26 | 1.56 (1.20–2.04) | 0.001 | 0.002 | 2.74 (1.51–5.00) | 0.0007 | 0.001 | 1.61 (1.12–2.31) | 0.009 | 0.013 |
SNP . | Case subjects . | . | . | . | Control subjects . | . | . | . | Allele 2 vs. 1 . | . | . | Genotype 22 vs. 11 . | . | . | Genotypes 22 + 12 vs. 11 . | . | . | ||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
. | 11 . | 12 . | 22 . | MAF . | 11 . | 12 . | 22 . | MAF . | OR (95% CI) . | P . | P(corr) . | OR (95% CI) . | P . | P(corr) . | OR (95% CI) . | P . | P(corr) . | ||||||||||||
rs509035 | |||||||||||||||||||||||||||||
Matched | 47 | 59 | 17 | 0.38 | 173 | 166 | 21 | 0.29 | 1.50 (1.10–2.03) | 0.009 | 0.013 | 2.98 (1.46–6.10) | 0.002 | 0.003 | 1.50 (0.99–2.27) | 0.058 | 0.085 | ||||||||||||
Total study | 49 | 65 | 17 | 0.38 | 424 | 375 | 69 | 0.30 | 1.45 (1.11–1.90) | 0.007 | 0.010 | 2.13 (1.16–3.91) | 0.013 | 0.019 | 1.60 (1.10–2.33) | 0.015 | 0.022 | ||||||||||||
rs572169 | |||||||||||||||||||||||||||||
Matched | 50 | 62 | 17 | 0.37 | 182 | 157 | 25 | 0.28 | 1.49 (1.11–2.01) | 0.009 | 0.013 | 2.48 (1.24–4.94) | 0.009 | 0.013 | 1.58 (1.05–2.38) | 0.028 | 0.041 | ||||||||||||
Total study | 52 | 67 | 17 | 0.37 | 451 | 371 | 80 | 0.29 | 1.42 (1.09–1.85) | 0.010 | 0.015 | 1.84 (1.02–3.35) | 0.042 | 0.062 | 1.62 (1.12–2.34) | 0.011 | 0.016 | ||||||||||||
rs519384 | |||||||||||||||||||||||||||||
Matched | 56 | 60 | 15 | 0.34 | 207 | 155 | 16 | 0.25 | 1.59 (1.18–2.16) | 0.003 | 0.004 | 3.47 (1.62–7.44) | 0.0009 | 0.001 | 1.62 (1.09–2.42) | 0.018 | 0.026 | ||||||||||||
Total study | 59 | 65 | 15 | 0.34 | 501 | 361 | 58 | 0.26 | 1.48 (1.13–1.94) | 0.004 | 0.006 | 2.20 (1.17–4.12) | 0.012 | 0.018 | 1.62 (1.13–2.33) | 0.008 | 0.012 | ||||||||||||
rs512692 | |||||||||||||||||||||||||||||
Matched | 57 | 55 | 15 | 0.33 | 200 | 154 | 15 | 0.25 | 1.51 (1.11–2.06) | 0.008 | 0.012 | 3.51 (1.62–7.61) | 0.0009 | 0.001 | 1.45 (0.97–2.18) | 0.070 | 0.103 | ||||||||||||
Total study | 60 | 60 | 15 | 0.33 | 482 | 356 | 56 | 0.26 | 1.41 (1.07–1.86) | 0.014 | 0.021 | 2.15 (1.15–4.04) | 0.015 | 0.022 | 1.46 (1.02–2.11) | 0.040 | 0.059 | ||||||||||||
rs863441 | |||||||||||||||||||||||||||||
Matched | 56 | 56 | 18 | 0.35 | 204 | 149 | 13 | 0.24 | 1.74 (1.28–2.37) | 0.0003 | 0.0004 | 5.04 (2.33–10.9) | 0.00001 | 0.00002 | 1.66 (1.11–2.49) | 0.013 | 0.019 | ||||||||||||
Total study | 59 | 61 | 18 | 0.35 | 493 | 355 | 55 | 0.26 | 1.56 (1.20–2.04) | 0.001 | 0.002 | 2.74 (1.51–5.00) | 0.0007 | 0.001 | 1.61 (1.12–2.31) | 0.009 | 0.013 |
SNPs shown as rs numbers from the dbSNP database. Odds ratios (ORs) with 95% CIs are based on the Armitage test for trend. MAF, minor allele frequency; P(corr), P value corrected for multiple testing.
dbSNP, single nucleotide polymorphism public database; FBAT, family-based association test; GHSR, growth hormone secretagogue receptor; LD, linkage disequilibrium; LVH, left ventricular hypertrophy; MONICA, Monitoring Trends and Determinants in Cardiovascular Disease; QTDT, quantitative transmission disequilibrium test; QTL, quantitative trait locus; SNP, single nucleotide polymorphism; TDT, transmission disequilibrium test; TOPS, Take Off Pounds Sensibly.
Article Information
This study was supported by National Institutes of Health Grants R01-DK-54026 and M01-RR-00058 and grants from the Deutsche Forschungsgemeinschaft (BA2227/1-1) and the Deutsche Akademie fuer Naturforscher Leopoldina (BMBF-LPD9901/8-72, to M.F.).
We would like to thank all TOPS members and all volunteers of the MONICA Augsburg study for their participation. We also would like to acknowledge Jacqueline Marks, Roland James, and John Blimke for assistance in recruiting and phenotyping the patients, establishing the databases, and preparing the DNA samples. Moreover, we thank Jill Stein for helping to establish the invader assays as well as Gretta Borchardt, Regina Cole, Martina Köhler, and Josef Simon for excellent technical assistance and managing the lab supply.