<title>Introduction</title> <p>Ghrelin is predominantly produced by the stomach and is the only known orexigenic hormone and endogenous ligand for the growth-hormone-secretagogue receptor (GHSR) ##REF##17167473##[1]##, ##REF##10604470##[2]##. Ghrelin is composed of 28 amino acids and is uniquely modified by the addition of an octanoyl group to the serine residue at position three. This acylation is necessary for ghrelin to bind to the GHSR and to cross the blood-brain barrier ##REF##15788704##[3]##.</p> <p>Ghrelin is implicated in growth hormone release ##REF##10604470##[2]##, energy balance ##REF##11483998##[4]##, food intake, and long-term regulation of body weight ##REF##11057670##[5]##, ##REF##11196643##[6]##. A major target for ghrelin is the arcuate nucleus of the hypothalamus, where ghrelin acts on appetite-stimulating groups of neurons, which express the potent orexigenic neuropeptides neuropeptide Y (NPY) and agouti-related protein (AGRP) ##REF##15058279##[7]##. These NPY/AGRP-containing neurons also express ghrelin receptors ##REF##10567856##[8]##.</p> <p>The <italic>GHSR</italic> gene, located on chromosome 3q26.31, encodes a protein that belongs to the family of G protein-coupled receptors ##REF##8688086##[9]##. Two distinct <italic>GHSR</italic> transcripts are known: GHSR type 1a acts as ghrelin receptor, whereas type 1b is a truncated form of the GHSR 1a and pharmacologically inactive ##REF##8688086##[9]##. GHSR 1a is expressed mainly in the hypothalamus and pituitary ##REF##15180951##[10]##, ##REF##12050285##[11]##, whereas the GHSR type 1b mRNA expression has also been found in many peripheral organs ##REF##15180951##[10]##, including immune cells ##REF##11549663##[12]##, ##REF##15232612##[13]##, indicating that ghrelin may have multiple functions in these tissues but with thus far unknown importance.</p> <p>Models of disrupted ghrelin signalling can show complex energy homeostasis-related phenotypes. <italic>GHSR</italic>-deficient mice show resistance to diet-induced obesity when supplied with a high-fat diet early in life, which implicates that functional ghrelin signalling is required for the full development of diet-induced obesity ##REF##17167473##[1]##, ##REF##16322794##[14]##.</p> <p>In humans, the <italic>GHSR</italic> gene lies within a quantitative trait locus strongly linked to multiple phenotypes related to obesity and the metabolic syndrome ##REF##15616037##[15]##. Baessler <italic>et al.</italic> have described a “susceptible” and “non-susceptible” haplotype for obesity of five <italic>GHSR</italic> single nucleotide polymorphisms (SNPs) ##REF##15616037##[15]##. Vartiainen <italic>et al.</italic>
##REF##15080774##[16]## showed an association between insulin values and a <italic>GHSR</italic> variant. A recent study reports a weak association between one <italic>GHSR</italic> variant and BMI ##UREF##0##[17]##. In recently conducted genome-wide association (GWA) studies, <italic>GHSR</italic> was not among the few genes associated with obesity ##REF##18454148##[18]##, ##REF##17434869##[19]## or type 2 diabetes ##REF##17463248##[20]##, ##REF##17668382##[21]##. Although GWA studies are very elegantly conducted to unravel the genetics of chronic diseases, candidate gene analyses also can often contribute to our knowledge about these disorders. Thus, we investigated SNPs in the <italic>GHSR</italic> gene and their association with obesity, glucose and insulin metabolism, and the conversion from impaired glucose tolerance (IGT) to type 2 diabetes in participants of the Finnish Diabetes Prevention Study (DPS). We demonstrated an association between one <italic>GHSR</italic> promoter SNP and weight loss during a 3-year follow-up and consequently showed that this SNP might have functional consequences on the expression of the <italic>GHSR</italic> gene.</p>
<title>Materials and Methods</title> <title>Subjects and study design</title> <p>The Finnish DPS is an intervention study carried out in five centres in Finland (clinical trial reg. no: NCT00518167). The main aim of the study was to assess the efficacy of an intensive individually designed diet and exercise program to prevent or delay the onset of type 2 diabetes in subjects with IGT. The study design, subjects, inclusion and exclusion criteria and intervention program have been described earlier in detail ##REF##10440120##[35]##, ##REF##11333990##[36]##. In brief, 522 overweight (BMI≥25 kg/m<sup>2</sup>) men and women aged 40 to 65 years with IGT participated in the study. DNA samples were available from 507 individuals (intervention group 259 and control group 248 subjects). IGT was defined as a 2-h plasma glucose concentration of 7.8 to 11.0 mmol/l (oral glucose tolerance test [OGTT] 75 g) with a fasting plasma glucose concentration less than 7.8 mmol/l ##UREF##1##[37]##. The persons enrolled in the study were randomly assigned to one of the two treatment modalities, the intervention group or the conventional care control group. The study protocol was approved by the Ethics Committee of the National Public Health Institute in Helsinki, Finland, and conducted in accordance with the guidelines proposed in the Declaration of Helsinki. All participants gave written informed consent. We certify that all applicable institutional and governmental regulations concerning the ethical use of human volunteers were followed during this research.</p> <title>Selection of <italic>GHSR</italic> SNPs and genotype analysis</title> <p>For selection of <italic>GHSR</italic> SNPs for genotype analysis the HapMap ##REF##16255080##[38]##, NCBI ##REF##17170002##[39]## and Ensembl ##REF##17148474##[40]## databases were used. Originally, six tag SNPs from two haploblocks covering 15 kb of the gene including three SNPs in haploblock 1 in the 5′ end of the gene (rs474225, rs490683 (alternative name in HapMap rs863441), rs9819506), and three SNPs in haploblock 2, one in the coding region in exon 1 leading to a synonymous amino acid substitution (rs495225), one in the intron (rs509035) and one in the 3′ end of the <italic>GHSR</italic> gene (rs565105) were selected (##FIG##0##Fig. 1A##). In addition, all SNPs in haploblock 1 in the 5′ regulatory region (spanning over 17 kb, chr 3: 173648897–17365897; NCBI accession number for <italic>GHSR</italic> NM_198407.1) available in the HapMap database and located in putative TF-binding sites were tested for their functionality in gelshift experiments. SNPs which showed different protein binding <italic>in vitro</italic> were genotyped additionally (rs6772676).</p> <p>Genotyping was performed with TaqMan Allelic Discrimination Assays (Applied Biosystems, Foster City, CA, USA). The PCR amplification was performed in a GeneAmp PCR system 2700 (Applied Biosystems, Foster City, CA, USA) in following conditions: 95°C for 10 min and 40 cycles of denaturation 92°C for 15 sec and annealing/extension 60°C for 1 min. Allele-specific fluorescence was detected on an ABI Prism 7000 sequence detector. A subset of randomly selected samples representing 6.3% of the study cohort was repeated.</p> <title>Promoter analysis</title> <p>
<italic>In silico</italic> promoter analysis of the sequences containing all SNPs in haploblock 1 for putative TF-binding sites was performed using Genomatix MatInspector Release professional 7.4.8.2 ##REF##15860560##[22]##. We determined all SNPs that lie in possible TF-binding sites and consequently all TF-binding sites that are disrupted by either one of each SNP allele. We used a threshold of 0.75 for core similarity and the optimized setting for matrix similarity.</p> <title>Protein extraction</title> <p>Hypothalami of three male Sprague-Dawley rats (supplier Harlan, Netherlands) were pooled and the frozen tissue was homogenized in 0.5 ml of low salt buffer (10 mM Hepes, pH 7.9, 1.5 mM MgCl<sub>2</sub>, 50 mM KCl, 0.5 mM dithiothreitol, proteinase inhibitors) using a rotor-stator homogenizer (ART Labortechnik, Müllheim, Germany). Cells were then resuspended in 2–3 packed volumes of the same low salt buffer with 0.5% Nonidet P-40 and homogenized by pipetting on ice. The nuclear pellet was subsequently resuspended in one packed cell volume of high salt buffer (10 mM Hepes, pH 7.9, 25% glycerol, 420 mM NaCl, 1.5 mM MgCl<sub>2</sub>, 0.2 mM EDTA, 0.5 mM dithiothreitol, proteinase inhibitors). Extraction was performed for 30 min on ice.</p> <title>Gelshift assay</title> <p>Gelshift assays were performed with 10 µg protein of the nuclear extracts. The proteins were incubated for 15 min in a total volume of 20 µl of binding buffer (150 mM KCl, 1 mM dithiothreitol, 25 ng/ml herring sperm DNA, 5% glycerol, 10 mM Hepes, pH 7.9). Constant amounts (1 ng) of <sup>32</sup>P-labeled double-stranded oligonucleotides (50,000 cpm) containing either one of the respective SNP alleles (supplementary ##SUPPL##0##Table S1##) were then added, and incubation was continued for 20 min at room temperature. Protein-DNA complexes were resolved by electrophoresis in 8% non-denaturing polyacrylamide gels (mono- to bisacrylamide ratio 19∶1) in 0.5×TBE (45 mM Tris, 45 mM boric acid, 1 mM EDTA, pH 8.3) for 90 min at 200 V and quantified on a FLA-3000 reader (Fuji, Tokyo, Japan) using ScienceLab99 software (Fuji).</p> <title>Statistical analysis</title> <p>The significances of differences in genotype and allele frequencies were analyzed using a chi-square test and a two-sided Fisher's exact test, respectively. Normal distributions were tested with the Kolmogorov-Smirnov test with Lilliefors correction. To normalize the skewed distributions, logarithmic transformations or reciprocal transformations were applied when needed. The differences in continuous variables between genotypes of single SNPs at baseline were evaluated with the general linear model (GLM) for the univariate analysis of variance (ANOVA). Biochemical measurements were adjusted for age, sex and BMI as main effects and weight-related variables for age and sex only. When variables could not be transformed to be normally distributed, Kruskal-Wallis test was used. GLM repeated measures ANOVA was used to analyze differences between genotypes of each SNP during the follow-up period concerning obesity-related and insulin- and glucose-metabolism related phenotypes. A final model was built manually with a backward selection procedure with following covariates included: age at baseline, sex, study group and BMI at baseline, as appropriate. Genotype*group interaction was tested and groups were analyzed separately if the interaction term was statistically significant. P-values are for trend between all three genotype groups and pairwise comparisons based on estimated marginal means were adjusted for multiple comparisons with Bonferroni correction. Relative changes in weight and measures of glucose and insulin metabolism from baseline to year 3 were analyzed with Kruskal-Wallis test, because data were not normally distributed even after applying several transformations. A p-value of 0.05 or lower was considered statistically significant. Correction for multiple hypothesis testing was performed with FDR using the Q-value 1.0 software ##REF##12883005##[41]##. π<sub>0</sub> was estimated with bootstrap method using a λ range from 0 to 0.9 by 0.05. The FDR for each p-value is reported as q-value. The q-value threshold is the proportion of significant features that turn out to be false leads. Data are given as means±SD, unless otherwise stated. All data were analyzed using SPSS for Windows 14.0 (SPSS Inc., Chicago, IL, USA). LD statistics were calculated and haplotype blocks were visualized by using the Haploview software ##REF##15297300##[42]##.</p>
<title>Results</title> <p>We genotyped 7 SNPs in the <italic>GHSR</italic> gene, four of which are located in the 5′-region, one in exon 1, one intronic and one SNP in the 3′-region (##FIG##0##Fig. 1A##). Rs490683, which was associated with weight loss and measures of glucose metabolism, is in strong linkage disequilibrium (LD) with rs509035 (##FIG##0##Fig. 1B##).</p> <title>Association with obesity-related phenotypes</title> <p>At baseline, no significant differences were observed between the genotypes of all seven SNPs for obesity-related traits (##TAB##0##Table 1## for rs490683, data not shown for other SNPs). However, using repeated measures ANOVA we observed that during the 3-year follow-up weight was significantly different between genotypes of rs9819506 (##FIG##1##Fig. 2A##) and tended to differ for BMI (p/q = 0.056/0.041).</p> <p>During the 3-year intervention. all subjects being analyzed together, individuals with <italic>rs490683-CC</italic> genotype lost significantly more weight than others (##FIG##1##Fig. 2B##). When the two intervention groups were analyzed separately, differences in weight loss according to the genotypes were only observed in the intervention group (##FIG##1##Fig. 2B##).</p> <p>No significant differences in weight loss during the intervention period were observed for other SNPs (data not shown).</p> <p>We carried out haplotype analysis with the 4 genotyped SNPs in haploblock 1 (rs6772676, rs474225, rs490683 and rs9819506, respectively; ##FIG##0##Fig. 1A##). The haplotype analysis did not offer more information than the single SNP analysis, thus these results were not included in the current manuscript. Concerning weight, the G-allele of rs9819506 seems to be responsible for higher weight in all subjects. The haplotype analysis regarding weight change did not reveal any significant results, largely due to the fact that not all relevant haplotypes could be tested due to their low frequency (data not shown).</p> <title>Association with phenotypes related to glucose metabolism</title> <p>At baseline, no significant associations with measures of glucose metabolism and the seven genotyped <italic>GHSR</italic> polymorphisms were observed (##TAB##0##Table 1## for rs490683, data not shown for other SNPs).</p> <p>Differences between genotypes were observed with the rs6772676 SNP for fasting (p/q = 0.038/0.035) and 2-h plasma glucose (p/q = 0.006/0.035), with heterozygotes having the highest levels in the whole study population when data were analyzed longitudinally. In similar analysis for rs490683, 2-h plasma glucose (##FIG##2##Fig. 3A##) and 2-h serum insulin levels (p/q = 0.050/0.041) differed significantly according to different genotypes with <italic>rs490683-CC</italic> subjects tending to have the lowest 2-h glucose and insulin levels throughout the follow-up. Three-year changes in fasting plasma glucose, 2-h plasma glucose and 2-h serum insulin in the entire study population differed between rs490683 genotypes (##FIG##2##Fig. 3B, C and E##, respectively). When the study groups were analyzed separately, significant differences were found between the rs490683 genotypes in fasting and 2-h plasma glucose and fasting and 2-h serum insulin (##FIG##2##Fig. 3B–E##) in the control group, but not in the intervention group. Levels of glucose and insulin decreased in persons with <italic>rs490683-CC</italic> genotype during the follow-up. Furthermore, differences between genotypes were also observed with the rs509035 SNP for 2-h plasma glucose and fasting insulin in the whole study population when data were analyzed longitudinally (p/q = 0.014/0.035 and p/q = 0.037/0.035, respectively). Changes in parameters of glucose metabolism differed according to rs509035 in the whole study population (fasting plasma glucose: p/q = 0.046/0.086, 2-h plasma glucose: p/q = 0.011/0.069, 2-h serum insulin: p/q = 0.037/0.086) and in the control group separately (2-h plasma glucose: p/q = 0.038/0.267, 2-h serum insulin: p/q = 0.014/0.152), with persons with <italic>rs509035-TT</italic> genotype showing a decrease in those parameters during the follow-up.</p> <p>None of the <italic>GHSR</italic> polymorphisms were statistically significantly associated with risk of type 2 diabetes (data not shown).</p> <title>Differential protein binding in the 5′-regulatory region</title> <p>We used the Genomatix MatInspector software ##REF##15860560##[22]## to identify putative transcription factor (TF)-binding sites in the 5′-region (haploblock 1) of the <italic>GHSR</italic> gene that may be implicated in the regulation of the gene expression. We screened 20 kb upstream of the translation initiation site and identified 8 SNPs that potentially disrupted one or more TF-binding sites in either one of the alleles or both. Subsequently, we tested oligonucleotides containing either of the alleles of each SNP (supplementary ##SUPPL##0##Table S1##) to identify possible differences in protein binding between the alleles of each SNP in gelshift assays using rat hypothalami as a TF protein source. The SNPs rs6772676 and rs490683 showed different protein binding when both alleles were compared (##FIG##3##Fig. 4##). Specifically, when the protein binding of the “wild-type” allele was normalized to 100%, <italic>rs6772676-T</italic> allele showed lower protein binding (only 28%) compared to the <italic>rs6772676-C</italic> allele. Furthermore, <italic>rs490683-G</italic> allele showed higher protein binding (745%) compared to that of the <italic>rs490683-C</italic> allele. The G-allele of rs490683 is present in the core sequence of the putative binding site of the TF nuclear factor 1 (NF-1), which is disrupted when G is replaced by C (##TAB##1##Table 2##).</p>
<title>Discussion</title> <p>We have shown that two SNPs (rs490683 and rs9819506) in the promoter region of the human <italic>GHSR</italic> gene are associated with the change in body weight and BMI. The <italic>GHSR</italic> SNP rs9819506, which has not been examined in previous studies, showed an association and individuals with <italic>rs9819506-AA</italic> genotype showed lowest values for body weight and a trend for lowest BMI during a 3-year follow-up. Interestingly, subjects with <italic>rs490683-CC</italic> genotype displayed highest weight loss in the whole study group as well as in the intervention group, when assessed separately. Previously, Baessler <italic>et al.</italic>
##REF##15616037##[15]## showed an association between <italic>GHSR</italic> haplotypes and obesity. The so-called “susceptible” haplotype tested in that study was found more often in obese individuals and included <italic>rs490683-C</italic> and <italic>rs509035-T</italic> alleles. These results might be considered contradictory to ours, where we see the <italic>rs490683-CC</italic> genotype as “beneficial” genotype. However, it is difficult to compare data from a longitudinal intervention study with individuals at high risk of type 2 diabetes describing the effect of the certain genotype on the weight development and a cross-sectional study carried out in a different population. A more recent study, however, could not point out a clear relationship between <italic>GHSR</italic> SNPs and obesity ##UREF##0##[17]##. Wang <italic>et al.</italic>
##REF##14715843##[23]## showed that obese children and adolescents had a higher <italic>rs495225-T</italic> allele frequency than underweight subjects, but this trend could not be confirmed in their further studies. Thus, they concluded that there is no evidence for involvement of <italic>GHSR</italic> SNPs in body weight regulation. In our study, SNPs rs495225 and rs565105 were not unequivocally associated with measures of obesity or type 2 diabetes.</p> <p>Vartiainen <italic>et al.</italic>
##REF##15080774##[16]## showed that subjects with <italic>rs495225-CC</italic> genotype had highest area under the insulin curve values and IGFBP-1 concentrations. In the present study, in the DPS population rs6772767, rs490683 and rs509035 were associated with several measures of glucose metabolism over four consecutive measurements. Beneficial changes in glucose metabolism could be observed specifically in individuals with <italic>rs490683-CC</italic> and <italic>rs509035-TT</italic> genotypes, especially in the control group. The control group in the DPS trial can be considered as a prospective cohort of people at high risk of diabetes. Since changes in measures of glucose metabolism during the study period were only observed in the control group, it indicates that the <italic>rs490683-CC</italic> genotype and the <italic>rs509035-TT</italic> genotype are beneficial or protective genotypes. A reduction in body weight was, however, evident when the whole study population was analysed, but this was mostly explained by the change in the intervention group. This finding is in some discrepancy with the association shown with measures of glucose metabolism. Thus, it might be that the different association with body weight and glucose metabolism are mediated via different mechanisms. In addition, we cannot exclude the impact of the moderate population size contributing to these varying associations. Unfortunately, we do not have the possibility to increase the size of the DPS population, but hopefully in the future it is possible to carry out joint analyses with the data from several other similar type of lifestyle intervention studies. However, we feel that the nature of DPS, as being a longitudinal study with an intensive lifestyle intervention in a carefully phenotyped high risk population, is unique and therefore our results offer valuable contribution and addition to the knowledge about weight and glucose control.</p> <p>Rs490683 and rs509035 are in strong LD with each other and this may further explain the similar results concerning glucose metabolism phenotypes. Rs509035 is located in the intron of the <italic>GHSR</italic> gene, thus it was not included in our analysis for putative TF binding sites. Nevertheless, it cannot be excluded that this SNP has functional consequences, but we did not test this. Thus, based on our results we may speculate that rs490683 could be the causal SNP accounting for the effects observed.</p> <p>As a limitation we need to point out that due to the distribution of these data, only non-parametric tests could be carried out where adjustments for covariates were not possible. To account for multiple testing and to control false positive rate, or type I error, we calculated the false discovery rate (FDR), expressed as q-value, for given p-values. Associations with a significant p-value, but a high q-value must be cautiously interpreted.</p> <p>To test the potential functional relevance of SNPs in the 5′-region of <italic>GHSR</italic> we first examined the regions where SNPs are possibly disrupting putative TF-binding sites with <italic>in silico</italic> promoter analysis. We then tested oligonucleotides with each allele of each SNP in gelshift assays using nuclear protein extracts from rat hypothalamus. We chose to use hypothalami of male rats as a comparable protein source, due to the fact that this tissue is unavailable from a human source. The DNA sequences (oligos) used for the gelshift experiments were designed to match the human sequence. We observed that nuclear proteins were binding to the sequence containing the <italic>rs490683-G</italic> allele with much higher affinity than to that of <italic>rs490683-C</italic> allele. At this position a putative NF-1 binding site exists and is disrupted by the SNP rs490683, changing the sequence from GCCA to CCCA. The NF-1 family of site-specific DNA-binding proteins functions both as cellular TFs in the regulation of gene expression and as replication factors for adenovirus DNA replication ##REF##10831836##[24]##. Binding specificity of NF-1 is equivalent to both the TGGCA-binding protein ##REF##6328417##[25]## and the CAAT-box TF ##REF##3024847##[26]##, ##REF##3398920##[27]##. NF-1 protein binds as a dimer to the dyad symmetric consensus sequence TTGGC(N<sub>5</sub>)<named-content content-type="gene">GCCAA</named-content> on duplex DNA, although NF-1 can also bind specifically to individual half sites (<named-content content-type="gene">TTGGC</named-content> or <named-content content-type="gene">GCCAA</named-content>) with somewhat reduced affinity ##REF##10831836##[24]##. NF-1 has been shown to activate transcription ##REF##10831836##[24]## and it can be speculated that in individuals with <italic>rs490683-GG</italic> genotype, where the NF-1 half site is intact, GHSR expression is enhanced, thus potentially leading to an increase in receptor signalling and ultimately to an increase in appetite ##REF##12876464##[28]##. In fact, our study shows that subjects with <italic>rs490683-GG</italic> genotype lost less weight during a lifestyle intervention than individuals with <italic>rs490683-CC</italic> genotype.</p> <p>It has been shown recently that GHSR signals with ∼50% activity even in the absence of an agonist ##REF##12907757##[29]##. That implies that control of the expression level of the receptor is directly correlated to signalling activity ##REF##15058279##[7]##. It has been shown that during prolonged fasting, GHSR expression in the hypothalamus is increased, which could contribute to the amplification of ghrelin action ##REF##12876464##[28]## and could be expected to result in a ghrelin-independent increase in receptor signalling and thereby an increase in appetite ##REF##15058279##[7]##. The <italic>GHSR</italic> polymorphism rs490683 was associated with weight change in our study population and the <italic>rs4906863-GG</italic> genotype, where the NF-1 site is intact, might lead to increased ghrelin receptor expression.</p> <p>The alleles of rs519384 did not show differential protein binding, indicating that the SNP is not in an important region for GHSR gene expression. Although rs9819506 was associated with weight and BMI, gelshift assays showed no difference in protein binding between the alleles, which may be due to species differences. Although we could observe a difference in protein binding between the alleles of rs6772676, it cannot explain the observed association with glucose levels.</p> <p>TFs and more importantly their DNA binding domains are highly conserved between humans and rodents. For example, the nucleic acid sequence of the rat NF-1 is 99% identical to the human sequence and on amino acid level there is a 91% overlap between the sequences. The variations in the rat NF-1 protein (compared to the human sequence) are not located within its DNA binding domain, thus they are unlikely to affect its binding to DNA. Due to the tissue unavailability from humans the rat protein extract provides a good model to test for differential protein binding to the alleles because of the high degree of conservation between the two species. We used hypothalami of male rats only, but obviously the autosomes of male and female rats are identical. Only genes located on the Y-chromosome will differ between the two sexes, whereas e.g. NF-1 is located on chromosome 17. Some TFs might show a difference in their expression levels. Thus, the gelshift bands might be stronger or weaker comparing males and females but the clear difference of protein binding to the different SNP alleles will still occur.</p> <p>To our knowledge, no GWA study for obesity has been carried out in Finns alone; however, samples from Finnish studies have been used in replication studies (e.g. ##REF##18454148##[18]##, ##REF##17434869##[19]##). The primary GWA studies investigating the genetics of type 2 diabetes are available. The studies regarding type 2 diabetes in Finnish subjects ##REF##17463248##[20]##, ##REF##17668382##[21]##, ##REF##18372903##[30]## have used the Illumina HumanHap300 Genotyping BeadChips containing over 317,000 SNPs. Two SNPs in the <italic>GHSR</italic> gene overlap with the SNPs investigated in our study, that is rs9819506 and rs565105. One more SNP, rs10513702, is in high LD (r<sup>2</sup> = 0.966) with rs6772676, which was also genotyped in our study. Other GWA studies investigating obesity have used variably the Affymetrix GeneChip Human Mapping 10K, 100K or 500K Array Set ##REF##18159244##[31]##, ##REF##17658951##[32]##, ##REF##18325910##[33]##, ##REF##17903300##[34]##. In the 500K Array Set are 6 SNPs of the <italic>GHSR</italic> gene included and one of those overlaps with the SNPs genotyped in our study (rs509035). The other <italic>GHSR</italic> SNPs are not in LD with SNPs which we genotyped. In GWA studies stringent quality control criteria are applied for the given setting. Thus it is possible that the few SNPs used in this study from the <italic>GHSR</italic> gene did not match those criteria (e.g. call rate, Hardy-Weinberg equilibrium for the given study population).</p> <p>In summary, results of our study suggest that polymorphisms in the <italic>GHSR</italic> gene are longitudinally associated with obesity and glucose metabolism in a prospective study setting. Genetic variation in the promoter may have an effect on appetite regulation by modulating gene expression of GHSR and subsequently ghrelin receptor signalling. This study is the first to show potential functional consequences of a <italic>GHSR</italic> SNP. Further studies are warranted to confirm our results and to find SNPs that are causal and have clinical relevance.</p>
<p>Conceived and designed the experiments: UM LP CC MU. Performed the experiments: UM TD. Analyzed the data: UM TD AMT. Contributed reagents/materials/analysis tools: LP MK JL JE CC JT MU. Wrote the paper: UM TD AMT. Read and commented on the paper: LP MK JL JE CC JT MU.</p> <title>Background</title> <p>Ghrelin may influence the development of obesity through its role in the control of energy balance, food intake, and regulation of body weight. The effects of ghrelin are mediated via the growth hormone secretagogue receptor (GHSR).</p> <title>Methodology/Principal Findings</title> <p>We genotyped 7 single nucleotide polymorphisms (SNPs) in the <italic>GHSR</italic> gene and assessed the association between those SNPs and obesity and type 2 diabetes-related phenotypes from 507 middle-aged overweight persons with impaired glucose tolerance participating in the Finnish Diabetes Prevention Study (DPS). Additionally, we performed <italic>in silico</italic> screening of the 5′-regulatory region of <italic>GHSR</italic> and evaluated SNPs disrupting putative transcription factor (TF) binding sites <italic>in vitro</italic> with gelshift assays to determine differences in protein binding between different alleles of SNPs. Rs9819506 in the promoter region of <italic>GHSR</italic> was associated with body weight (p = 0.036); persons with <italic>rs9819506-AA</italic> genotype having the lowest body weight. Individuals with <italic>rs490683-CC</italic> genotype displayed highest weight loss in the whole study population (p = 0.032). The false discovery rate for these results was &lt;10%. Rs490683 and rs509035 were associated with several measures of glucose and insulin metabolism during the follow-up. Rs490683 may be a functional SNP, since gelshift experiments showed differential protein binding between the alleles, with higher binding to the G-allele. <italic>Rs490683-C</italic> may disrupt a putative binding site for the TF nuclear factor 1 (NF-1), thus <italic>rs4906863-GG</italic> genotype where the NF-1 site is intact may lead to a higher <italic>GHSR</italic> gene expression.</p> <title>Conclusion/Significance</title> <p>Polymorphisms in the <italic>GHSR</italic> promoter may modify changes in body weight during long-term lifestyle intervention and affect ghrelin receptor signalling through modulation of <italic>GHSR</italic> gene expression.</p>
<title>Supporting Information</title>
<p>We thank Dr. Jari Nissinen and Dr. Maria Vlasova for providing rat hypothalami.</p> <p>The DPS group has been previously described in the earlier publications ##REF##10440120##[35]##, ##REF##11333990##[36]##. The Finnish Diabetes Prevention Group consists additionally of following people: Timo Valle, Helena Hämäläinen, Pirjo Ilanne-Parikka, Sirkka Keinänen-Kiukaanniemi.</p>
<fig id="pone-0002941-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002941.g001</object-id><label>Figure 1</label><caption><title>The human <italic>GHSR</italic> gene and location of SNPs and their LD measures.</title><p>A. Schematic representation of the human <italic>GHSR</italic> gene indicating the locations of analyzed single nucleotide polymorphisms (SNPs) that belong to two haploblocks. (black star), SNPs genotyped in the Finnish Diabetes Prevention Study population. # SNPs from HapMap database in transcription factor (TF) binding sites and tested in gelshift assays. B. Standardized pair-wise linkage disequilibria (LD) measures of D' and r<sup>2</sup> for seven <italic>GHSR</italic> single nucleotide polymorphisms (SNPs) of the Finnish Diabetes Prevention Study. Major/minor alleles and the minor allele frequency (in parentheses) are included below each SNP identification number. For each cell, the number on the top line is D' and the number on the bottom line is r<sup>2</sup>. Strength of LD is indicated in greyscale with black denoting strong and white weak LD.</p></caption></fig> <fig id="pone-0002941-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002941.g002</object-id><label>Figure 2</label><caption><title>Weight and weight loss according to rs9819506 and rs490683 genotypes in DPS.</title><p>A. Weight according to rs9819506 genotypes in the entire study population. Mean±SEM are estimated marginal means calculated from repeated measures ANOVA with sex and age as significant covariates. <italic>rs9819506-GG</italic> (n = 171, filled circles); <italic>rs9819506-GA</italic> (n = 200, open cirlces); <italic>rs9819506-AA</italic> (n = 59, filled triangles); pair-wise comparisons: <italic>GG</italic> vs. <italic>GA</italic> p = 1.000, <italic>GG</italic> vs. <italic>AA</italic> p = 0.030, <italic>GA</italic> vs. <italic>AA</italic> p = 0.116. B. Relative weight loss according to the different genotypes of rs490683 in all subjects and separately in intervention and control group. <italic>rs490683-CC</italic> (black bars); <italic>rs490683-CG</italic> (striped bars); <italic>rs490683-GG</italic> (white bars). The q-value denotes the false discovery rate (FDR) and is reported only for p&lt;0.1. n.s., non-significant.</p></caption></fig> <fig id="pone-0002941-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002941.g003</object-id><label>Figure 3</label><caption><title>Glucose metabolism measures according to rs490683 genotypes in DPS.</title><p>A. 2-h plasma glucose according to rs490683 genotypes in all subjects. Mean±SEM are estimated marginal means calculated from repeated measures ANOVA with group, age and BMI at baseline as significant covariates. <italic>rs490683-CC</italic> (n = 42, filled circles); <italic>rs490683-CG</italic> (n = 171, open circles); <italic>rs490683-GG</italic> (n = 210, filled triangles). Pair-wise comparisons all subjects: <italic>CC</italic> vs. <italic>CG</italic> p = 0.020, <italic>CC</italic> vs. <italic>GG</italic> p = 0.295, <italic>CG</italic> vs. <italic>GG</italic> p = 0.190. B–E. Relative changes fasting plasma glucose (B), 2-h plasma glucose (C), fasting serum insulin (D) and 2-h serum insulin (E) according to the different genotypes of rs490683 in all subjects and separately in intervention and control group. <italic>rs490683-CC</italic> (black bars); <italic>rs490683-CG</italic> (striped bars); <italic>rs490683-GG</italic> (white bars). The q-value denotes the false discovery rate (FDR) and is reported only for p&lt;0.1. n.s., non-significant.</p></caption></fig> <fig id="pone-0002941-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002941.g004</object-id><label>Figure 4</label><caption><title>Differential protein binding of rat hypothalamus nuclear protein extract.</title><p>Gelshift experiments were performed with <sup>32</sup>P-labeled oligonucleotides containing one or the other allele of respective single nucleotide polymorphism (SNP). Numbers below the gels indicate the means of three independent experiments with standard deviation (SD) in parenthesis. Protein binding of the “wild-type” allele was normalized to 100%. Representative gels are shown. NS denotes non-specific DNA-protein interactions, which are not due to specific binding between the labelled DNA and the proteins contained in the hypothalamus extract. <sup>*</sup> p&lt;0.05; <sup>**</sup> p&lt;0.001.</p></caption></fig>
<table-wrap id="pone-0002941-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002941.t001</object-id><label>Table 1</label><caption><title>Baseline characteristics of the DPS study subjects according to the rs490683 polymorphism of the <italic>GHSR</italic> gene.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td colspan="3" align="left" rowspan="1">Genotype</td><td align="left" rowspan="1" colspan="1">P-value<xref ref-type="table-fn" rid="nt105">*</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<italic>rs490683-CC</italic> (n)</td><td align="left" rowspan="1" colspan="1">
<italic>rs490683-CG</italic> (n)</td><td align="left" rowspan="1" colspan="1">
<italic>rs490683-GG</italic> (n)</td><td align="left" rowspan="1" colspan="1"/></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Sex (M/F)</td><td align="left" rowspan="1" colspan="1">13/33 (46)</td><td align="left" rowspan="1" colspan="1">68/147 (215)</td><td align="left" rowspan="1" colspan="1">85/161 (246)</td><td align="left" rowspan="1" colspan="1">0.635<xref ref-type="table-fn" rid="nt103">b</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Age (years)</td><td align="left" rowspan="1" colspan="1">54±7 (46)</td><td align="left" rowspan="1" colspan="1">56±7 (215)</td><td align="left" rowspan="1" colspan="1">55±7 (246)</td><td align="left" rowspan="1" colspan="1">0.246<xref ref-type="table-fn" rid="nt104">c</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Weight (kg)</td><td align="left" rowspan="1" colspan="1">88.5±13.8 (46)</td><td align="left" rowspan="1" colspan="1">84.4±12.9 (215)</td><td align="left" rowspan="1" colspan="1">87.4±15.2 (246)</td><td align="left" rowspan="1" colspan="1">0.080<xref ref-type="table-fn" rid="nt102">a</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1">BMI (kg/m<sup>2</sup>)</td><td align="left" rowspan="1" colspan="1">32.4±4.9 (46)</td><td align="left" rowspan="1" colspan="1">30.8±4.2 (215)</td><td align="left" rowspan="1" colspan="1">31.4±4.7 (246)</td><td align="left" rowspan="1" colspan="1">0.132<xref ref-type="table-fn" rid="nt102">a</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1">Waist circumference (cm)</td><td align="left" rowspan="1" colspan="1">101.5±10.2 (46)</td><td align="left" rowspan="1" colspan="1">100.7±10.8 (213)</td><td align="left" rowspan="1" colspan="1">101.6±11.3 (246)</td><td align="left" rowspan="1" colspan="1">0.062</td></tr><tr><td align="left" rowspan="1" colspan="1">Waist-to-hip ratio</td><td align="left" rowspan="1" colspan="1">0.91±0.07 (46)</td><td align="left" rowspan="1" colspan="1">0.92±0.07 (213)</td><td align="left" rowspan="1" colspan="1">0.92±0.08 (246)</td><td align="left" rowspan="1" colspan="1">0.415</td></tr><tr><td align="left" rowspan="1" colspan="1">Fasting plasma glucose (mmol/l)</td><td align="left" rowspan="1" colspan="1">6.3±0.7 (46)</td><td align="left" rowspan="1" colspan="1">6.1±0.7 (215)</td><td align="left" rowspan="1" colspan="1">6.1±0.8 (246)</td><td align="left" rowspan="1" colspan="1">0.653</td></tr><tr><td align="left" rowspan="1" colspan="1">2-h plasma glucose (mmol/l)</td><td align="left" rowspan="1" colspan="1">9.0±1.6 (46)</td><td align="left" rowspan="1" colspan="1">9.0±1.6 (215)</td><td align="left" rowspan="1" colspan="1">8.8±1.4 (246)</td><td align="left" rowspan="1" colspan="1">0.198</td></tr><tr><td align="left" rowspan="1" colspan="1">Fasting serum insulin (pmol/l)</td><td align="left" rowspan="1" colspan="1">16.5±8.7 (40)</td><td align="left" rowspan="1" colspan="1">14.5±6.5 (193)</td><td align="left" rowspan="1" colspan="1">14.7±7.9 (228)</td><td align="left" rowspan="1" colspan="1">0.574</td></tr><tr><td align="left" rowspan="1" colspan="1">2-h serum insulin (pmol/l)</td><td align="left" rowspan="1" colspan="1">105.4±61.2 (39)</td><td align="left" rowspan="1" colspan="1">92.9±57.3 (192)</td><td align="left" rowspan="1" colspan="1">95.9±71.7 (227)</td><td align="left" rowspan="1" colspan="1">0.728</td></tr><tr><td align="left" rowspan="1" colspan="1">HOMA-IR (mmol×mU×l<sup>−2</sup>)</td><td align="left" rowspan="1" colspan="1">4.61±2.43 (40)</td><td align="left" rowspan="1" colspan="1">4.01±2.01 (193)</td><td align="left" rowspan="1" colspan="1">4.09±2.47 (228)</td><td align="left" rowspan="1" colspan="1">0.550</td></tr><tr><td align="left" rowspan="1" colspan="1">HOMA-IS (U/mol)</td><td align="left" rowspan="1" colspan="1">126.8±78.2 (40)</td><td align="left" rowspan="1" colspan="1">119.6±73.2 (193)</td><td align="left" rowspan="1" colspan="1">115.7±58.7 (228)</td><td align="left" rowspan="1" colspan="1">0.567</td></tr><tr><td align="left" rowspan="1" colspan="1">Total serum cholesterol (mmol/l)</td><td align="left" rowspan="1" colspan="1">5.81±0.91 (46)</td><td align="left" rowspan="1" colspan="1">5.64±0.96 (215)</td><td align="left" rowspan="1" colspan="1">5.55±0.89 (245)</td><td align="left" rowspan="1" colspan="1">0.147</td></tr><tr><td align="left" rowspan="1" colspan="1">HDL cholesterol (mmol/l)</td><td align="left" rowspan="1" colspan="1">1.25±0.32 (46)</td><td align="left" rowspan="1" colspan="1">1.21±0.29 (215)</td><td align="left" rowspan="1" colspan="1">1.20±0.29 (245)</td><td align="left" rowspan="1" colspan="1">0.523</td></tr><tr><td align="left" rowspan="1" colspan="1">LDL cholesterol (mmol/l)</td><td align="left" rowspan="1" colspan="1">3.68±0.79 (45)</td><td align="left" rowspan="1" colspan="1">3.64±0.86 (215)</td><td align="left" rowspan="1" colspan="1">3.59±0.82 (244)</td><td align="left" rowspan="1" colspan="1">0.681</td></tr><tr><td align="left" rowspan="1" colspan="1">Triglycerides (mmol/l)</td><td align="left" rowspan="1" colspan="1">1.95±0.92 (46)</td><td align="left" rowspan="1" colspan="1">1.74±0.76 (215)</td><td align="left" rowspan="1" colspan="1">1.68±0.75 (245)</td><td align="left" rowspan="1" colspan="1">0.139</td></tr><tr><td align="left" rowspan="1" colspan="1">Systolic blood pressure (mmHg)</td><td align="left" rowspan="1" colspan="1">138±17 (45)</td><td align="left" rowspan="1" colspan="1">137.4±17.8 (213)</td><td align="left" rowspan="1" colspan="1">139±18 (244)</td><td align="left" rowspan="1" colspan="1">0.458</td></tr><tr><td align="left" rowspan="1" colspan="1">Diastolic blood pressure (mmHg)</td><td align="left" rowspan="1" colspan="1">87±8 (45)</td><td align="left" rowspan="1" colspan="1">85±10 (213)</td><td align="left" rowspan="1" colspan="1">86±10 (244)</td><td align="left" rowspan="1" colspan="1">0.592</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0002941-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002941.t002</object-id><label>Table 2</label><caption><title>Single nucleotide polymorphisms (SNPs) in the <italic>GHSR</italic> promoter and putative transcription factor (TF) binding sites.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">SNP</td><td align="left" rowspan="1" colspan="1">Allele</td><td align="left" rowspan="1" colspan="1">Putative TF binding sites</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">rs6772676<xref ref-type="table-fn" rid="nt106">1</xref>
<sup>,</sup>
<xref ref-type="table-fn" rid="nt107">2</xref>
</td><td align="left" rowspan="1" colspan="1">C</td><td align="left" rowspan="1" colspan="1">Activating transcription factor, <italic>ATF-1</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Pax-3 paired domain protein, expressed in embryogenesis, mutations correlate to Waardenburg Syndrome</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">T</td><td align="left" rowspan="1" colspan="1">PBX - HOXA9 binding site</td></tr><tr><td align="left" rowspan="1" colspan="1">rs9881097</td><td align="left" rowspan="1" colspan="1">T</td><td align="left" rowspan="1" colspan="1">Avian C-type LTR TATA box</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">C</td><td align="left" rowspan="1" colspan="1">Pu.1 (Pu120) Ets-like TF identified in lymphoid B-cells</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">PAX 2/5/8 binding site</td></tr><tr><td align="left" rowspan="1" colspan="1">rs474225<xref ref-type="table-fn" rid="nt106">1</xref>
</td><td align="left" rowspan="1" colspan="1">T</td><td align="left" rowspan="1" colspan="1">Fkh-domain factor FKHRL1 (FOXO)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Interferon regulatory factor 3 (IRF-3)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">C</td><td align="left" rowspan="1" colspan="1">-</td></tr><tr><td align="left" rowspan="1" colspan="1">rs490683<xref ref-type="table-fn" rid="nt106">1</xref>
<sup>,</sup>
<xref ref-type="table-fn" rid="nt107">2</xref>
</td><td align="left" rowspan="1" colspan="1">C</td><td align="left" rowspan="1" colspan="1">Orphan nuclear receptor rev-erb alpha (NR1D1), monomer binding site</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">G</td><td align="left" rowspan="1" colspan="1">Non-palindromic nuclear factor I binding sites</td></tr><tr><td align="left" rowspan="1" colspan="1">rs11918879</td><td align="left" rowspan="1" colspan="1">C</td><td align="left" rowspan="1" colspan="1">Brn-5, POU-VI protein class (also known as emb and CNS-1)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">T</td><td align="left" rowspan="1" colspan="1">Brn-5, POU-VI protein class (also known as emb and CNS-1)</td></tr><tr><td align="left" rowspan="1" colspan="1">rs477441</td><td align="left" rowspan="1" colspan="1">A</td><td align="left" rowspan="1" colspan="1">Hepatic nuclear factor 1</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">C</td><td align="left" rowspan="1" colspan="1">Ecotropic viral integration site 1 encoded factor, amino-terminal zinc finger domain</td></tr><tr><td align="left" rowspan="1" colspan="1">rs9819506<xref ref-type="table-fn" rid="nt106">1</xref>
</td><td align="left" rowspan="1" colspan="1">G</td><td align="left" rowspan="1" colspan="1">TEA domain-containing factors, transcriptional enhancer factors 1,3,4,5</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Octamer binding site (OCT1/OCT2)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Octamer-binding factor 1, POU-specific domain</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">A</td><td align="left" rowspan="1" colspan="1">-</td></tr><tr><td align="left" rowspan="1" colspan="1">rs519384</td><td align="left" rowspan="1" colspan="1">A</td><td align="left" rowspan="1" colspan="1">Liver enriched Cut - Homeodomain TF HNF6 (ONECUT)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">T</td><td align="left" rowspan="1" colspan="1">Liver enriched Cut - Homeodomain TF HNF6 (ONECUT)</td></tr></tbody></table></alternatives></table-wrap>
<supplementary-material content-type="local-data" id="pone.0002941.s001"><label>Table S1</label><caption><p>Oligonucleotides used in gelshift assays. List of oligonucleotides which were used to test the difference in protein binding among different GHSR SNP alleles.</p><p>(0.04 MB DOC)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><p>Values are means±SD; ANOVA for comparison among three genotype groups, adjusted for age, sex and BMI as main effects.</p></fn><fn id="nt102"><label>a</label><p>p-value adjusted for age and sex.</p></fn><fn id="nt103"><label>b</label><p>Chi-square Test.</p></fn><fn id="nt104"><label>c</label><p>Kruskal-Wallis Test.</p></fn><fn id="nt105"><label>*</label><p>FDR was &gt;60% for all p-values.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt106"><label>1</label><p>SNPs genotyped from DPS population.</p></fn><fn id="nt107"><label>2</label><p>different protein binding in gelshift assay.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>This work has been financially supported by grants from the Academy of Finland (number 209445 to MK, numbers 40758 and 117844 to MU; numbers 38387 and 46558 to JT), the EVO-fund of Kuopio University Hospital (numbers 5179, 5198, 5952 and 5926 to MU), the Ministry of Education of Finland, the Finnish Diabetes Foundation, Sigrid Juselius Foundation, the Austrian Academy of Sciences, the RTN NucSys of the European Union, the Finnish Graduate School on Applied Bioscience and the Yrjö Jahnsson Foundation.</p></fn></fn-group>
<graphic xlink:href="pone.0002941.g001"/> <graphic xlink:href="pone.0002941.g002"/> <graphic id="pone-0002941-t001-1" xlink:href="pone.0002941.t001"/> <graphic xlink:href="pone.0002941.g003"/> <graphic xlink:href="pone.0002941.g004"/> <graphic id="pone-0002941-t002-2" xlink:href="pone.0002941.t002"/>
<media xlink:href="pone.0002941.s001.doc"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["17"], "element-citation": ["\n"], "surname": ["Dossus", "McKay", "Canzian", "Wilkening", "Rinaldi"], "given-names": ["L", "JD", "F", "S", "S"], "year": ["2008"], "article-title": ["Polymorphisms of genes coding for ghrelin and its receptor in relation to anthropometry, circulating levels of IGF-I and IGFBP-3, and breast cancer risk: a case-control study nested within the European Prospective Investigation into Cancer and Nutrition (EPIC)."], "source": ["Carcinogenesis"], "comment": ["doi:10.1093/carcin/bgn083"]}, {"label": ["37"], "element-citation": ["\n"], "collab": ["WHO"], "year": ["1985"], "article-title": ["Diabetes mellitus. Report of a WHO Study Group."], "source": ["World Health Organ Tech Rep Ser"], "volume": ["727"], "fpage": ["i"], "lpage": ["xii"], "page-range": ["1\u2013113"]}]
PLoS One. 2008 Aug 13; 3(8):e2941