- Research article
- Open Access
- Open Peer Review
Association between a rare SNP in the second intron of human Agouti related protein gene and increased BMI
© Kalnina et al; licensee BioMed Central Ltd. 2009
- Received: 31 October 2008
- Accepted: 14 July 2009
- Published: 14 July 2009
The agouti related protein (AGRP) is an endogenous antagonist of the melanocortin 4 receptor and is one of the most potent orexigenic factors. The aim of the present study was to assess the genetic variability of AGRP gene and investigate whether the previously reported SNP rs5030980 and the rs11575892, a SNP that so far has not been studied with respect to obesity is associated with increased body mass index (BMI).
We determined the complete sequence of the AGRP gene and upstream promoter region in 95 patients with severe obesity (BMI > 35 kg/m2). Three polymorphisms were identified: silent mutation c.123G>A (rs34123523) in the second exon, non-synonymous mutation c.199G>A (rs5030980) and c.131-42C>T (rs11575892) located in the second intron. We further screened rs11575892 in a selected group of 1135 and rs5030980 in group of 789 participants from the Genome Database of Latvian Population and Latvian State Research Program Database.
The CT heterozygotes of rs11575892 had significantly higher mean BMI value (p = 0.027). After adjustment for age, gender and other significant non-genetic factors (presence of diseases), the BMI levels remained significantly higher in carriers of the rs11575892 T allele (p = 0.001). The adjusted mean BMI value of CC genotype was 27.92 ± 1.01 kg/m2 (mean, SE) as compared to 30.97 ± 1.03 kg/m2 for the CT genotype. No association was found between rs5030980 and BMI.
This study presents an association of rare allele of AGRP polymorphism in heterozygous state with increased BMI. The possible functional effects of this polymorphism are unclear but may relate to splicing defects.
- Body Mass Index
- Anorexia Nervosa
- Complete Linkage Disequilibrium
- Human Genome Variation Society
- Agouti Related Protein
The central melanocortin system and its components play an important role in regulation of food intake and energy balance and it is known as a major pathway of leptin signaling (reviewed in [1, 2]). Agouti related protein (AGRP) is predominantly located in ARC neurons that co-express neuropeptide Y (NPY), also a strong mediator of regulation of energy balance. AGRP is evolutionary conserved gene  that acts as endogenous melanocortin 4 receptor (MC4R) and melanocortin 3 receptor (MC3R) antagonist (inverse agonist) with strong orexigenic effect (for review see ). Both ICV injection [5, 6] and transgenic overexpression  of AGRP leads to increased food consumption and obesity in mice. Elevated plasma levels of AGRP have been reported in obese men  in comparison to non-obese men. In addition, AGRP levels are increased in humans as a result of fasting , indicating a possible role of AGRP in periphery.
The AGRP gene is located on chromosome 16q22, consists of four exons and expresses two alternatively spliced variants that differ in presence or absence of the 5' non-coding exon. The long transcript variant is predominantly expressed in the hypothalamus while the shorter transcript lacking the 5' untranslated exon is found in peripheral tissues such as testis, lung and kidney [10, 11]. A few single nucleotide polymorphisms (SNPs) have previously been identified in AGRP gene that are associated with energy homeostasis disorder phenotypes. Two SNPs in the promoter region -3019G>A and -38C>T have been found in individuals of African origin but not in Caucasians . Both these SNPs are in complete linkage disequilibrium (LD) and are associated with leanness, decreased risk for type two diabetes (T2DM) and reduced macronutrient intake [12–14]. Another SNP in the third exon of the AGRP gene that changes the alanine at position 67 to a threonine (rs5030980) has been identified, and is so far, found in Caucasians only. This SNP has been shown to be associated with anorexia nervosa . Homozygosity of Ala67 did not have any association to a metabolic phenotype at mean age of 25 years but it was associated with late onset obesity at mean age of 53 years . Patients homozygous for the Thr67 allele had significantly reduced body fatness . In the present study we demonstrate association of a C>T polymorphism at position c.131-42C>T (rs11575892) with increased BMI in cohort of patients from the Latvian Genome Database.
The study was based on data and samples from the Genome Database of Latvian Population (Biobank 1), the disease based biobank comprising 1173 subjects and 439 subjects from population based Latvian State Research Program Database (Biobank 2). Both biobanks were collected by using the same protocols, questionnaires and sample treatment techniques. Briefly; individuals were required to be over the age of 18; the health status was recorded, the diagnoses was based on approved clinical criteria according to the ICD-10 codes (International Classification of Diseases); anthropometric measurements (including weight and stature), ethnic, social, environmental information and familial health status was acquired based on self-reported questionnaire. Written informed consent was obtained from all participants. We selected 95 of the most obese (BMI > 35 kg/m2) subjects from the database for initial screening using sequencing. The study group for the subsequent genotyping of rs11575892, was selected from all available subjects with validated health records and questionnaires in biobank by exclusion of patients with diseases whose outcome or treatment may influence BMI, including all types of cancer and diseases of thyroid gland, but not excluding patients with cardiovascular diseases. Patients with metabolic diseases were included in the study as well. In total 1135 subjects were selected based on these criteria. This sample group included 73 individuals from the 95 that were sequenced. rs5030980 was genotyped subsequently in total of 789 samples comprising 350 samples from Biobank 1 and all 439 samples included in study from Biobank 2. The 348 samples from initial selection were not available for genotyping of rs5030980. Study protocol was approved by Central Medical Ethics Committee of Latvia.
The AGRP gene containing genomic DNA region including 5' upstream region and the entire coding region (positions from -1102 bp to +1027 respective to start codon) were amplified in two individual PCRs. Both strands of amplification products were directly sequenced using set of six primers (primer sequences and reaction conditions can be provided upon request from authors). In total 95 overweight individuals (BMI > 35 kg/m2) were sequenced with 100% success rate. All chromatograms were manually inspected using Contig Express software of Vector NTI Advance 9.0 package. Presence of polymorphisms was confirmed by opposite strand analysis. Restriction enzyme fragment polymorphism analysis of rs11575892 was used for genotyping the 698 samples from Biobank 1. Primer with a sequence mismatch in 3' part; 5'-CCCTCCCCTGGGAGGTGGG AG-3' (mismatched position underlined, wild type cuts only) was used to create mutation conditional AluI restriction sites in PCR reaction with reverse primer; 5'-GGTGAGGGAGTTTGGTGCTGG-3'. The amplified PCR product contained additional wt AluI site allowing to control the restriction reaction for false positive results. The AluI restriction products were visualized on PAAG gels. The samples with polymorphisms were verified by direct sequencing. Genotyping of the rs11575892 in all other samples and genotyping of rs5030980 was performed using pre-designed TaqMan SNP Genotyping Assays (Applied Biosystems, Foster City, California, USA) on ABI 7500 Real-Time PCR system (Applied Biosystems) according to the supplier's recommendations. In addition, the 5' upstream region of the AGRP sequence for all carriers of the rs11575892 from Biobank 1 was determined by sequencing. The designations of the SNP's are based on recommendation by the Human Genome Variation Society.
Statistical analyses were performed with the statistical package SPSS (Standard Version 12.0.0;SPSS, Chicago, IL, USA). In contrast to the natural BMI values, the logarithmically transformed BMI values displayed normal distribution and were further used for all quantitative analysis. A one-way ANOVA was used in univariate analysis to test for association of log-BMI values with other categorical variables. Linear regression analysis was used to study relationship between log-BMI and age. General linear model (GLM) approach was applied to study association between each SNP and log-BMI adjusting for main effects of age, gender and those factors that showed association with logBMI under separate ANOVA tests (presence of hypertension, angina pectoris, myocardial infarction, heart failure, dyslipidemia and T2DM). Adjusted mean logBMI values, corresponding CI intervals and significance of association were estimated according to the GLM model. Logistic regression analysis adjusting for age and gender was used for dichotomous trait analysis. Power calculations were performed using Quanto v.1.2.3 . Our sample size provided 80% power (at α = 0.05) to detect difference between rs11575892 genotypes in logBMI of >0.040 assuming minor allele frequency of 0.015 and to detect difference between rs5030980 genotypes in logBMI of >0.028 assuming minor allele frequency of 0.044. PLINK 1.00 software http://pngu.mgh.harvard.edu/purcell/plink/ was used to perform Hardy-Weinberg test, LD calculations, haplotype based association and permutation test for ANOVA and multivariate linear regression using label-swapping between SNP and logBMI values but leaving intact correlation with other covariates. 100000 permutations were performed for each analysis and we used corrected (EMP2) p-values. These values are corrected based on calculation of the proportion of permutations in which any of the test statistics exceeds the particular observed statistic and are more stringent then uncorrected p-values.
Direct sequencing of the coding and 5' untranslated regions of the human AGRP gene in 95 obese patients with BMI > 35 kg/m2 revealed three polymorphisms: silent mutation c.123G>A (rs34123523) in second exon (minor allele frequency MAF = 0.042), non-synonymous mutation c.199G>A (rs5030980) in third exon (MAF = 0.042) and c.131-42C>T (rs11575892) located in second intron (MAF = 0.016). As reported previously  rs34123523 and rs5030980 were in complete linkage disequilibrium (LD) with each other (D' = 1.00, r2 = 1.00) while none of these subjects with both SNPs had rs11575892. The allele frequency of rs5030980 was in agreement (MAF from 0.036 to 0.05) with previous reports [13, 20]. None of the SNPs currently present in the SNP databases or any new SNP were found in the 5' untranslated region of the AGRP gene. The SNPs reported for this region are however with low frequency (MAF = 0.005 for rs34731556 and rs34018897) or found in subjects of African origin (rs5030981), that may explain the low number of SNPs in our study.
Characteristics of the study population.
Mean ± SD
56.68 ± 12.52
28.07 ± 5.36
1.4407 ± 0.0804
Association of BMI with two SNPs
Mean logBMI ± SE
1.446 ± 0.005
1.491 ± 0.014
Mean BMIc (kg/m2)
Mean logBMI ± SE
1.458 ± 0.007
1.458 ± 0.011
1.471 ± 0.076
Mean BMIc (kg/m2)
Categorical analysis b
Genotype distribution (n)
OR (95% CI)
BMI <= median
BMI > median
These results indicate that association of rs11575892 with increased BMI is not due to possible confounding by specific diseases or other factors in study cohort. Nevertheless we cannot completely exclude the role of different environmental variables such as social background, diet and physical activity that could have impact on BMI, but such parameters were not accessible for present analysis. Due to low frequency of the rs11575892 polymorphism, we could only estimate the effect of this genetic variant in its heterozygous state, as we did not detect any homozygous individuals who may have pronounced effect of this SNP on BMI. The power analyses show that the present sample size had >80% power to detect changes in our logBMI of 0.045. Since rs11575892 has not previously been studied with respect to obesity and no association for this polymorphism has been found with anorexia by Vink et al, further studies with increased sample size or case-control studies with stratification for environmental factors such as social background, diet and physical activity would be valuable to verify the effects of rs1157892. There are reports that show effects of other AGRP polymorphisms on body weight are expressed in age dependant fashion . This study group consists mostly of elderly people (mean 59.0 years) and this could contribute to enhance the effect of this particular SNP on BMI.
In order to test previously reported associations of rs5030980 with decreased body fatness we subsequently performed genotyping in available samples (n = 789) that were used for genotyping of rs1157892. The estimated minor allele frequency of rs11575892 was 0.044 and we did not observe significant departure from Hardy-Weinberg equilibrium (p = 1.0). Pairwise LD analysis showed absence of LD between both SNPs (D' = 1.00, r2 = 0.001). No association of rs5030980 with BMI was found neither in GLM or logistic regression analysis (Table 2). We did not find any difference in mean values between GG and GA genotypes as opposite to previous report . We found only one AA genotype carrier in our study group with slightly increased BMI value (29.58 kg/m2). This is in disagreement with results of Marks et.al where carriers of AA displayed significantly lower BMI . Relatively low frequency of this SNP however permits to compare these effects with sufficient power. Haplotype analysis did not return significant difference in BMI values between three possible haplotypes (data not shown).
In summary, our results show novel association of a SNP in the second intron of AGRP with increased BMI suggesting that naturally occurring mutations in AGRP may be associated with enhanced effect of this peptide in the human population. This SNP could provide valuable information on the regulation of AGRP gene.
The studies were supported by the Latvian Council of Science Grant 05.0023.01, Latvian State Research Program 8.14, Swedish Research Council (VR, medicine), AFA insurance. We acknowledge Genome Database of Latvian Population, Latvian Biomedical Research and Study Center for providing data and DNA samples.
- Hillebrand JJ, Kas MJ, Adan RA: To eat or not to eat; regulation by the melanocortin system. Physiol Behav. 2006, 89 (1): 97-102. 10.1016/j.physbeh.2006.01.034.View ArticlePubMedGoogle Scholar
- Seeley RJ, Drazen DL, Clegg DJ: The critical role of the melanocortin system in the control of energy balance. Annu Rev Nutr. 2004, 24: 133-149. 10.1146/annurev.nutr.24.012003.132428.View ArticlePubMedGoogle Scholar
- Klovins J, Haitina T, Fridmanis D, Kilianova Z, Kapa I, Fredriksson R, Gallo-Payet N, Schioth HB: The melanocortin system in Fugu: determination of POMC/AGRP/MCR gene repertoire and synteny, as well as pharmacology and anatomical distribution of the MCRs. Mol Biol Evol. 2004, 21 (3): 563-579. 10.1093/molbev/msh050.View ArticlePubMedGoogle Scholar
- Stutz AM, Morrison CD, Argyropoulos G: The Agouti-related protein and its role in energy homeostasis. Peptides. 2005, 26 (10): 1771-1781. 10.1016/j.peptides.2004.12.024.View ArticlePubMedGoogle Scholar
- Rossi M, Kim MS, Morgan DG, Small CJ, Edwards CM, Sunter D, Abusnana S, Goldstone AP, Russell SH, Stanley SA, et al: A C-terminal fragment of Agouti-related protein increases feeding and antagonizes the effect of alpha-melanocyte stimulating hormone in vivo. Endocrinology. 1998, 139 (10): 4428-4431. 10.1210/en.139.10.4428.View ArticlePubMedGoogle Scholar
- Lu XY, Nicholson JR, Akil H, Watson SJ: Time course of short-term and long-term orexigenic effects of Agouti-related protein (86–132). Neuroreport. 2001, 12 (6): 1281-1284. 10.1097/00001756-200105080-00045.View ArticlePubMedGoogle Scholar
- Graham M, Shutter JR, Sarmiento U, Sarosi I, Stark KL: Overexpression of Agrt leads to obesity in transgenic mice. Nat Genet. 1997, 17 (3): 273-274. 10.1038/ng1197-273.View ArticlePubMedGoogle Scholar
- Katsuki A, Sumida Y, Gabazza EC, Murashima S, Tanaka T, Furuta M, Araki-Sasaki R, Hori Y, Nakatani K, Yano Y, et al: Plasma levels of agouti-related protein are increased in obese men. J Clin Endocrinol Metab. 2001, 86 (5): 1921-1924. 10.1210/jc.86.5.1921.View ArticlePubMedGoogle Scholar
- Hoggard N, Johnstone AM, Faber P, Gibney ER, Elia M, Lobley G, Rayner V, Horgan G, Hunter L, Bashir S, et al: Plasma concentrations of alpha-MSH, AgRP and leptin in lean and obese men and their relationship to differing states of energy balance perturbation. Clin Endocrinol (Oxf). 2004, 61 (1): 31-39. 10.1111/j.1365-2265.2004.02056.x.View ArticleGoogle Scholar
- Ollmann MM, Wilson BD, Yang YK, Kerns JA, Chen Y, Gantz I, Barsh GS: Antagonism of central melanocortin receptors in vitro and in vivo by agouti-related protein. Science. 1997, 278 (5335): 135-138. 10.1126/science.278.5335.135.View ArticlePubMedGoogle Scholar
- Shutter JR, Graham M, Kinsey AC, Scully S, Luthy R, Stark KL: Hypothalamic expression of ART, a novel gene related to agouti, is up-regulated in obese and diabetic mutant mice. Genes Dev. 1997, 11 (5): 593-602. 10.1101/gad.11.5.593.View ArticlePubMedGoogle Scholar
- Mayfield DK, Brown AM, Page GP, Garvey WT, Shriver MD, Argyropoulos G: A role for the Agouti-Related Protein promoter in obesity and type 2 diabetes. Biochem Biophys Res Commun. 2001, 287 (2): 568-573. 10.1006/bbrc.2001.5600.View ArticlePubMedGoogle Scholar
- Loos RJ, Rankinen T, Rice T, Rao DC, Leon AS, Skinner JS, Bouchard C, Argyropoulos G: Two ethnic-specific polymorphisms in the human Agouti-related protein gene are associated with macronutrient intake. Am J Clin Nutr. 2005, 82 (5): 1097-1101.PubMedGoogle Scholar
- Bonilla C, Panguluri RK, Taliaferro-Smith L, Argyropoulos G, Chen G, Adeyemo AA, Amoah A, Owusu S, Acheampong J, Agyenim-Boateng K, et al: Agouti-related protein promoter variant associated with leanness and decreased risk for diabetes in West Africans. Int J Obes (Lond). 2006, 30 (4): 715-721. 10.1038/sj.ijo.0803047.View ArticleGoogle Scholar
- Vink T, Hinney A, van Elburg AA, van Goozen SH, Sandkuijl LA, Sinke RJ, Herpertz-Dahlmann BM, Hebebrand J, Remschmidt H, van Engeland H, et al: Association between an agouti-related protein gene polymorphism and anorexia nervosa. Mol Psychiatry. 2001, 6 (3): 325-328. 10.1038/sj.mp.4000854.View ArticlePubMedGoogle Scholar
- Argyropoulos G, Rankinen T, Neufeld DR, Rice T, Province MA, Leon AS, Skinner JS, Wilmore JH, Rao DC, Bouchard C: A polymorphism in the human agouti-related protein is associated with late-onset obesity. J Clin Endocrinol Metab. 2002, 87 (9): 4198-4202. 10.1210/jc.2002-011834.View ArticlePubMedGoogle Scholar
- Marks DL, Boucher N, Lanouette CM, Perusse L, Brookhart G, Comuzzie AG, Chagnon YC, Cone RD: Ala67Thr polymorphism in the Agouti-related peptide gene is associated with inherited leanness in humans. Am J Med Genet A. 2004, 126 (3): 267-271. 10.1002/ajmg.a.20600.View ArticleGoogle Scholar
- Gauderman W, Morrison J: QUANTO 1.1: A computer program for power and sample size calculations for genetic-epidemiology studies. 2006, [http://hydra.usc.edu/gxe]Google Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, et al: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81 (3): 559-575. 10.1086/519795.View ArticlePubMedPubMed CentralGoogle Scholar
- van Rossum CT, Pijl H, Adan RA, Hoebee B, Seidell JC: Polymorphisms in the NPY and AGRP genes and body fatness in Dutch adults. Int J Obes (Lond). 2006, 30 (10): 1522-1528. 10.1038/sj.ijo.0803314.View ArticleGoogle Scholar
- Bai F, Rankinen T, Charbonneau C, Belsham DD, Rao DC, Bouchard C, Argyropoulos G: Functional dimorphism of two hAgRP promoter SNPs in linkage disequilibrium. J Med Genet. 2004, 41 (5): 350-353. 10.1136/jmg.2003.014092.View ArticlePubMedPubMed CentralGoogle Scholar
- Sozen MA, de Jonge LH, Greenway F, Ravussin E, Smith SR, Argyropoulos G: A rare mutation in AgRP, +79G>A, affects promoter activity. Eur J Clin Nutr. 2007, 61 (6): 809-812. 10.1038/sj.ejcn.1602585.View ArticlePubMedGoogle Scholar
- de Rijke CE, Jackson PJ, Garner KM, van Rozen RJ, Douglas NR, Kas MJ, Millhauser GL, Adan RA: Functional analysis of the Ala67Thr polymorphism in agouti related protein associated with anorexia nervosa and leanness. Biochem Pharmacol. 2005, 70 (2): 308-316. 10.1016/j.bcp.2005.04.033.View ArticlePubMedGoogle Scholar
- Creemers JW, Pritchard LE, Gyte A, Le Rouzic P, Meulemans S, Wardlaw SL, Zhu X, Steiner DF, Davies N, Armstrong D, et al: Agouti-related protein is posttranslationally cleaved by proprotein convertase 1 to generate agouti-related protein (AGRP)83-132: interaction between AGRP83-132 and melanocortin receptors cannot be influenced by syndecan-3. Endocrinology. 2006, 147 (4): 1621-1631. 10.1210/en.2005-1373.View ArticlePubMedGoogle Scholar
- Goto K, Inui A, Takimoto Y, Yuzuriha H, Asakawa A, Kawamura Y, Tsuji H, Takahara Y, Takeyama C, Katsuura G, et al: Acute intracerebroventricular administration of either carboxyl-terminal or amino-terminal fragments of agouti-related peptide produces a long-term decrease in energy expenditure in rats. Int J Mol Med. 2003, 12 (3): 379-383.PubMedGoogle Scholar
- Place RF, Li LC, Pookot D, Noonan EJ, Dahiya R: MicroRNA-373 induces expression of genes with complementary promoter sequences. Proc Natl Acad Sci USA. 2008, 105 (5): 1608-1613. 10.1073/pnas.0707594105.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/10/63/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.