PCSK1 rs6232 Is Associated with Childhood and Adult Class III Obesity in the Mexican Population

Abstract

Background

Common variants rs6232 and rs6235 in the PCSK1 gene have been associated with obesity in European populations. We aimed to evaluate the contribution of these variants to obesity and related traits in Mexican children and adults.

Methodology/Principal Findings

Rs6232 and rs6235 were genotyped in 2382 individuals, 1206 children and 1176 adults. Minor allele frequencies were 0.78% for rs6232 and 19.99% for rs6235. Rs6232 was significantly associated with childhood obesity and adult class III obesity (OR = 3.01 95%CI 1.64–5.53; P = 4×10−4 in the combined analysis). In addition, this SNP was significantly associated with lower fasting glucose levels (P = 0.01) and with increased insulin levels and HOMA-B (P = 0.05 and 0.01, respectively) only in non-obese children. In contrast, rs6235 showed no significant association with obesity or with glucose homeostasis parameters in any group.

Conclusion/Significance

Although rs6232 is rare in the Mexican population, it should be considered as an important risk factor for extreme forms of obesity.

About the Authors:

Marisela Villalobos-Comparán

Affiliations Departamento de Biología, Facultad de Química, Universidad Nacional Autónoma de México (UNAM), Mexico City, Mexico, Unidad de Biología Molecular y Medicina Genómica, Instituto Nacional de Ciencias Médicas y Nutrición “Salvador Zubirán” (INCMNSZ), Mexico City, Mexico

Hugo Villamil-Ramírez

Affiliations Departamento de Biología, Facultad de Química, Universidad Nacional Autónoma de México (UNAM), Mexico City, Mexico, Unidad de Biología Molecular y Medicina Genómica, Instituto Nacional de Ciencias Médicas y Nutrición “Salvador Zubirán” (INCMNSZ), Mexico City, Mexico

Teresa Villarreal-Molina

Affiliation: Instituto Nacional de Medicina Genómica (INMEGEN), Mexico City, Mexico

Elena Larrieta-Carrasco

Affiliation: Unidad de Biología Molecular y Medicina Genómica, Instituto Nacional de Ciencias Médicas y Nutrición “Salvador Zubirán” (INCMNSZ), Mexico City, Mexico

Paola León-Mimila

Affiliations Departamento de Biología, Facultad de Química, Universidad Nacional Autónoma de México (UNAM), Mexico City, Mexico, Unidad de Biología Molecular y Medicina Genómica, Instituto Nacional de Ciencias Médicas y Nutrición “Salvador Zubirán” (INCMNSZ), Mexico City, Mexico

Sandra Romero-Hidalgo

Affiliation: Instituto Nacional de Medicina Genómica (INMEGEN), Mexico City, Mexico

Leonor Jacobo-Albavera

Affiliations Departamento de Biología, Facultad de Química, Universidad Nacional Autónoma de México (UNAM), Mexico City, Mexico, Unidad de Biología Molecular y Medicina Genómica, Instituto Nacional de Ciencias Médicas y Nutrición “Salvador Zubirán” (INCMNSZ), Mexico City, Mexico

Adriana E. Liceaga-Fuentes

Affiliation: Clínica de Obesidad, Hospital Rubén Leñero, Mexico City, Mexico

Francisco J. Campos-Pérez

Affiliation: Clínica de Obesidad, Hospital Rubén Leñero, Mexico City, Mexico

Blanca E. López-Contreras

Affiliations Unidad de Biología Molecular y Medicina Genómica, Instituto Nacional de Ciencias Médicas y Nutrición “Salvador Zubirán” (INCMNSZ), Mexico City, Mexico, Instituto Nacional de Medicina Genómica (INMEGEN), Mexico City, Mexico

Teresa Tusié-Luna

Affiliations Unidad de Biología Molecular y Medicina Genómica, Instituto Nacional de Ciencias Médicas y Nutrición “Salvador Zubirán” (INCMNSZ), Mexico City, Mexico, Instituto de investigaciones Biomédicas, Universidad Nacional Autónoma de México (UNAM), Mexico City, Mexico

Blanca E. del Río-Navarro

Affiliation: Departamento de Alergia e Inmunología Clínica, Hospital Infantil de México Federico Gómez, Mexico City, Mexico

Carlos A. Aguilar-Salinas

Affiliation: Departmento de Endocrinología y Metabolismo, INCMNSZ, Mexico City, Mexico

Samuel Canizales-Quinteros

* E-mail: [email protected]

Affiliations Departamento de Biología, Facultad de Química, Universidad Nacional Autónoma de México (UNAM), Mexico City, Mexico, Unidad de Biología Molecular y Medicina Genómica, Instituto Nacional de Ciencias Médicas y Nutrición “Salvador Zubirán” (INCMNSZ), Mexico City, Mexico

Introduction

The prevalence of obesity has increased worldwide, including Mexico, where more than 70% of the adults and 26% of the children are overweight or obese [1]. Even though obesity has a strong genetic contribution, the identification of genes related to obesity risk has proven difficult [2]–[4]. The discovery of genes causing monogenic forms of obesity such as the prohormone convertase subtilisin/kexin type 1 gene (PCSK1) has greatly improved our understanding of the pathophysiology of obesity [2]–[4]. This gene encodes an enzyme expressed in neuroendocrine cells that converts inactive prohormones into functional key hormones that regulate central and/or peripheral energy metabolism. Although loss-of-function mutations in this gene causing childhood obesity and impaired glucose tolerance are rare [5]–[7], two common nonsynonymous variants (rs6232 and rs6235) were recently found to be strongly associated with childhood and adulthood obesity in European population [8]. Functional in vitro analysis of these variants revealed a significant impairment of the catalytic activity of the enzyme for rs6232 (N221D), but no enzyme activity alteration for rs6235 (S690T) [8].

Several studies have sought to replicate the association of these variants (mainly rs6235) with obesity and obesity-associated traits in Asian and European adult populations, with inconsistent results [9]–[15]. These population-based studies included only a reduced number of obese class III individuals and did not include children as the initial report of Benzinou et al. [8]. This may partially explain such inconsistencies, as genetic influences on BMI may be stronger precisely in children and individuals with class III obesity [2], [16]. Thus, the aim of the present study was to analyze the association of rs6232 and rs6235 with obesity and related traits in a case-control analysis of Mexican-Mestizo adults and children.

Materials and Methods

Subjects

The study included 1206 non-related Mexican-Mestizo children aged 5 to 12 years (596 boys and 610 girls), recruited at a summer camp for children of employees of the Mexican Health Ministry (Convivencia Infantil 2008–2009, Secretaría de Salud) and the Hospital Infantil de México. Body mass index (BMI) was calculated as weight in kilograms divided by height in meter squared. BMI z-scores and percentiles were calculated using age and sex specific BMI reference data, as recommended by the Centers for Disease Control and Prevention [17]. The population was stratified based on percentile BMI, 802 children were non-obese (BMI percentile<95th) and 404 were obese (BMI percentile≥95th).

The study also included 1176 unrelated Mexican Mestizos aged 18–82 years (806 nonpregnant women and 370 men): 788 individuals described by Villalobos-Comparán et al. [18] and 388 additional subjects recruited from Obesity Clinics at the INCMNSZ and Ruben Leñero Hospital. Two hundred and fifty four (21.59%) of these individuals had been diagnosed with type 2 diabetes (T2D) according to WHO (World Health Organization) criteria [19]. Individuals were grouped according to BMI: 562 were non-obese (BMI <30 kg/m2), 380 had class I/II obesity (30≤ BMI <40 kg/m2), and 234 had class III obesity (BMI ≥40 kg/m2). Biochemical parameters were measured in blood samples obtained after 12-h fast as previously described [20]. Homeostasis model assessment of beta-cell function (HOMA-B) and insulin sensitivity (HOMA-S) as measures of beta-cell function and insulin sensitivity were estimated using a computer model [21]. The characteristics of the children and adult populations are shown in Table 1.

[Figure omitted. See PDF.]

Table 1. Characteristics of the cases and controls in the children and adult populations.

https://doi.org/10.1371/journal.pone.0039037.t001

As part of the eligibility criteria, subjects with thyroid gland disease or showing weight instability three months prior to the study were excluded, as well as elderly subjects with signs of dementia by Mini-Mental state examination [22]. Only individuals born in Mexico whose parents and grandparents identified themselves as Mexican Mestizos were included. This project was approved by the Institutional Committee of Biomedical Research in Humans of the INCMNSZ. All adult participants and parents of the children provided written informed consent prior to their inclusion in the study.

Single Nucleotide Polymorphism Genotyping

Rs6232 and rs6235 were genotyped using TaqMan assays (ABI Prism 7900HT Sequence Detection System; Applied Biosystems, Foster City, CA). Genotyping call rate exceeded 97% per SNP and no discordant genotypes were observed in 25 duplicate samples. In addition, because the Mexican-Mestizo population is admixed, we analyzed 10 ancestry informative markers to rule out population stratification [18]. Genotyping was performed by KBiosciences (Hertfordshire, UK, http://www.kbioscience.co.uk/) using a KASPar assay system. Genotyping call rates of each ancestry informative marker exceeded 95%, and no discordant genotypes were observed in 54 duplicate samples. Deviation from Hardy–Weinberg equilibrium was not observed for rs6232 and rs6235 in any group (P>0.66 and P>0.21, respectively).

Statistical Analysis

Logistic regression was used to test for associations between the rs6232 and rs6235 SNPs and obesity. Children and adult combined odds ratios were estimated using the Mantel-Haenszel method. The AdmixMap program was used to test the possible effect of population stratification on associations of rs6232 with obesity only in the adult population [23]. Generalized linear regression was applied to test for associations of rs6232 and rs6235 with quantitative traits only in individuals without T2D. Because fasting serum insulin levels and HOMA indices were not normally distributed, they were log transformed for analysis. Interactions between the SNPs and age or gender were tested by including a two-way interaction term (SNP*age or SNP*gender) in the model. All analyses were adjusted for age and gender, and with other covariates as appropriate. The reported P-values are nominal and two-sided. Association analyses were performed with SPSS V15.0, statistical package; Chicago, IL.

Pairwise linkage disequilibrium (LD) between both SNPs was estimated using Haploview V3.2. (http://www.broad.mit.edu/mpg/haploview). Power calculations were performed using QUANTO software (http://hydra.usc.edu/gxe/).

Results

The study included a total of 2382 individuals, 1176 adults and 1206 children. Minor allele frequencies (MAFs) for rs6232 and rs6235 SNPs were 0.78%, and 19.99%, respectively. Both variants are in weak linkage disequilibrium (r2 = 0.30). Because risk genotypes frequencies showed no significant differences in lean and overweight individuals (children or adults; P>0.40), both groups were considered together as non-obese subjects for the analyses. SNP rs6232 was significantly associated with obesity in children (OR = 3.78, 95%CI 1.42–9.88; P = 7×10−3) and was significantly associated with class III obesity in adults (OR = 2.61, 95%CI 1.10–6.19; P = 0.02,) showing a trend of association after adjusting for admixture (P = 0.07). The odds ratio estimated in the combined analysis was 3.01 (95%CI 1.64–5.53; P = 4×10−4, Table 2). In contrast, the rs6235 variant was not associated with obesity in children or adults (OR = 1.16, 95%CI 0.96–1.42; P = 0.20, in the combined analysis, Table 2). No significant genotype-age or genotype-gender interactions on the risk of obesity were observed for any SNP (P>0.05 in all tested models).

[Figure omitted. See PDF.]

Table 2. Association of the rs6232 and rs6235 with obesity in children and adult populations.

https://doi.org/10.1371/journal.pone.0039037.t002

We also explored the effect of the rs6232 and rs6235 on glucose homeostasis in non-obese and obese children and in non-diabetic adults. The rs6232 risk allele was significantly associated with lower fasting glucose levels (P = 0.01) and with increased insulin levels and HOMA-B (P = 0.05 and 0.01, respectively) only in non-obese children (Table 3). No significant associations between rs6235 and glucose homeostasis parameters were found in children or adults (Table 4).

[Figure omitted. See PDF.]

Table 3. Glucose homeostasis parameters in nondiabetic children and adult populations according to PCSK1 rs6232.

https://doi.org/10.1371/journal.pone.0039037.t003

[Figure omitted. See PDF.]

Table 4. Glucose homeostasis parameters in nondiabetic children and adult populations according to PCSK1 rs6235.

https://doi.org/10.1371/journal.pone.0039037.t004

Discussion

Association with Obesity

Obesity is a complex disorder involving both genetic and environmental factors [2]. Although obesity is highly prevalent in Mexico, both in children and adults, studies on the genetic component of this disease in the Mexican population are scarce [18], [20], [24]. Recently, N221D (rs6232) and S690T (rs6235) PCSK1 nonsynonymous polymorphisms were found to contribute to the etiology of polygenic obesity in European populations [8]. In the present case-control study, the PCSK1 rs6235 was common but not associated with obesity in Mexican children or adults. This result is consistent with several recent studies in European and Chinese populations [9]–[14]. In contrast, rs6232 was infrequent in the Mexican population (0.78% as compared to 4–8% in Europeans), but was significantly associated with obesity in both children and adults. In fact, the association with obesity for G allele carriers was significant only in extreme phenotypes (childhood and adult class III obesity), and was higher than the risk previously reported (OR = 3.01 vs. OR = 1.34, respectively) [8]. Thus, it is likely to provide only a weak population-attributable risk for common obesity and not to be a major contributor to obesity in the general population of Mexico. However, the presence of this functional variant should be considered as a serious risk factor for extreme forms of obesity in the Mexican population, as has been recently reported for heterozygous PCSK1 mutations in Europeans [25].

A limitation of the study was low statistical power. For the combined analysis (children and adults), the present study had only 13.6% and 78.3% statistical power to detect previously reported associations of rs6232 and rs6235 with obesity, respectively [8]. Because statistical power to detect obesity subclasses was lower, we cannot rule out the possibility of rs6235 associations with extreme forms of obesity in this population.

Associations with Glucose Homeostasis

Both rs6235 and rs6232 have been recently associated with glucometabolic traits. Rs6235 was found to be associated with decreasing fasting glucose levels and increased HOMA-B [15], [26], while rs6232 was associated with decreased circulating post-prandial glucose and elevated glucagon levels in a Danish population-based study [27], and discordantly associated with decreased fasting insulin levels and reduced insulin sensitivity in German adults with increased risk of T2D [26]. In the present study rs6235 was not associated with glucometabolic traits in any group; however the rs6232 G-allele was significantly associated with decreased insulin sensitivity only in non-obese children. The inconsistencies of the effect of PCSK1 variants on glucose metabolism may result from different study designs and allele frequencies, particularly in the case of rs6232 which is very infrequent in the Mexican population. Additional studies are required to confirm and further characterize the effect of these polymorphisms on glucose homeostasis.

In conclusion, although the rs6232 polymorphism is rare in the Mexican population, it is significantly associated with extreme forms of obesity. However, it is likely not to be a major contributor to obesity in the general Mexican population. Further studies with larger samples sizes are needed to confirm this association.

Acknowledgments

We thank Salvador Ramírez-Jimenez and Luz E. Guillén-Pineda for their technical assistance.

Author Contributions

Conceived and designed the experiments: MVC SCQ. Performed the experiments: MVC HVR ELC PLM. Analyzed the data: MVC SRH TVM SCQ. Contributed reagents/materials/analysis tools: LJA AELF FJCP BELC TTL BERN CAAS SCQ. Wrote the paper: MVC TVM SRH CAAS SCQ.

Citation: Villalobos-Comparán M, Villamil-Ramírez H, Villarreal-Molina T, Larrieta-Carrasco E, León-Mimila P, Romero-Hidalgo S, et al. (2012) PCSK1 rs6232 Is Associated with Childhood and Adult Class III Obesity in the Mexican Population. PLoS ONE7(6): e39037. https://doi.org/10.1371/journal.pone.0039037

References

1. Secretaria de Salud, Instituto Nacional de Salud Pública (2006) Encuesta Nacional de Salud y Nutrición. Salud Secretaria deInstituto Nacional de Salud Pública2006Encuesta Nacional de Salud y Nutrición(ENSANUT). (ENSANUT).

2. Bell CG, Walley AJ, Froguel P (2005) The genetics of human obesity. Nat Rev Genet 6: 221–234.CG BellAJ WalleyP. Froguel2005The genetics of human obesity.Nat Rev Genet6221234

3. Farooqi S, O’Rahilly S (2006) Genetics of obesity in humans. Endocr Rev 27: 710–718.S. FarooqiS. O’Rahilly2006Genetics of obesity in humans.Endocr Rev27710718

4. Ramachandrappa S, Farooqi IS (2011) Genetic approaches to understanding human obesity. J Clin Invest 121: 2080–2086.S. RamachandrappaIS Farooqi2011Genetic approaches to understanding human obesity.J Clin Invest12120802086

5. Jackson RS, Creemers JW, Ohagi S, Raffin-Sanson ML, Sanders L, et al. (1997) Obesity and impaired prohormone processing associated with mutations in the human prohormone convertase 1 gene. Nat Genet 16: 303–306.RS JacksonJW CreemersS. OhagiML Raffin-SansonL. Sanders1997Obesity and impaired prohormone processing associated with mutations in the human prohormone convertase 1 gene.Nat Genet16303306

6. Jackson RS, Creemers JW, Farooqi IS, Raffin-Sanson ML, Varro A, et al. (2003) Small-intestinal dysfunction accompanies the complex endocrinopathy of human proprotein convertase 1 deficiency. J Clin Invest 112: 1550–1560.RS JacksonJW CreemersIS FarooqiML Raffin-SansonA. Varro2003Small-intestinal dysfunction accompanies the complex endocrinopathy of human proprotein convertase 1 deficiency.J Clin Invest11215501560

7. Farooqi IS, Volders K, Stanhope R, Heuschkel R, White A, et al. (2007) Hyperphagia and early-onset obesity due to a novel homozygous missense mutation in prohormone convertase 1/3. J Clin Endocrinol Metab 92: 3369–3373.IS FarooqiK. VoldersR. StanhopeR. HeuschkelA. White2007Hyperphagia and early-onset obesity due to a novel homozygous missense mutation in prohormone convertase 1/3.J Clin Endocrinol Metab9233693373

8. Benzinou M, Creemers JW, Choquet H, Lobbens S, Dina C, et al. (2008) Common nonsynonymous variants in PCSK1 confer risk of obesity. Nat Gen 40: 943–945.M. BenzinouJW CreemersH. ChoquetS. LobbensC. Dina2008Common nonsynonymous variants in PCSK1 confer risk of obesity.Nat Gen40943945

9. Kilpeläinen TO, Bingham SA, Khaw KT, Wareham NJ, Loos RJ (2009) Association of variants in the PCSK1 gene with obesity in the EPIC-Norfolk study. Hum Mol Genet 18: 3496–3501.TO KilpeläinenSA BinghamKT KhawNJ WarehamRJ Loos2009Association of variants in the PCSK1 gene with obesity in the EPIC-Norfolk study.Hum Mol Genet1834963501

10. Renström F, Payne F, Nordström A, Brito EC, Rolandsson O, et al. (2009) Replication and extension of genome-wide association study results for obesity in 4923 adults from northern Sweden. Hum Mol Genet 18: 1489–1496.F. RenströmF. PayneA. NordströmEC BritoO. Rolandsson2009Replication and extension of genome-wide association study results for obesity in 4923 adults from northern Sweden.Hum Mol Genet1814891496

11. Sandholt CH, Sparsø T, Grarup N, Albrechtsen A, Almind K, et al. (2010) Combined Analysis of 20 Common obesity susceptibility variants. Diabetes 59: 1667–1673.CH SandholtT. SparsøN. GrarupA. AlbrechtsenK. Almind2010Combined Analysis of 20 Common obesity susceptibility variants.Diabetes5916671673

12. Qi Q, Li H, Loos RJ, Liu C, Hu FB, et al. (2010) Association of PCSK1 rs6234 with obesity and related traits in a Chinese Han population. PLoS One 5: e10590.Q. QiH. LiRJ LoosC. LiuFB Hu2010Association of PCSK1 rs6234 with obesity and related traits in a Chinese Han population.PLoS One5e10590

13. Chang YC, Chiu YF, Shih KC, Lin MW, Sheu WH, et al. (2010) Common PCSK1 haplotypes are associated with obesity in the Chinese population. Obesity 18: 1404–1409.YC ChangYF ChiuKC ShihMW LinWH Sheu2010Common PCSK1 haplotypes are associated with obesity in the Chinese population.Obesity1814041409

14. Wen W, Cho YS, Zheng W, Dorajoo R, Kato N, et al. (2012) Meta-analysis identifies common variants associated with body mass index in east Asians. Nat Genet 44: 307–311.W. WenYS ChoW. ZhengR. DorajooN. Kato2012Meta-analysis identifies common variants associated with body mass index in east Asians.Nat Genet44307311

15. Strawbridge RJ, Dupuis J, Prokopenko I, Barker A, Ahlqvist E, et al. (2011) Genome-wide association identifies nine common variants associated with fasting proinsulin levels and provides new insights into the pathophysiology of type 2 diabetes. Diabetes 60: 2624–2634.RJ StrawbridgeJ. DupuisI. ProkopenkoA. BarkerE. Ahlqvist2011Genome-wide association identifies nine common variants associated with fasting proinsulin levels and provides new insights into the pathophysiology of type 2 diabetes.Diabetes6026242634

16. Wardle J, Carnell S, Haworth CM, Plomin R (2008) Evidence for a strong genetic influence on childhood adiposity despite the force of the obesogenic environment. Am J Clin Nutr 87: 398–404.J. WardleS. CarnellCM HaworthR. Plomin2008Evidence for a strong genetic influence on childhood adiposity despite the force of the obesogenic environment.Am J Clin Nutr87398404

17. Kuczmarski RJ, Ogden CL, Guo SS, Grummer-Strawn LM, Flegal KM, et al. (2002) CDC growth charts for the United States: methods and development. Vital Health Stat 246: 1–190.RJ KuczmarskiCL OgdenSS GuoLM Grummer-StrawnKM Flegal2002CDC growth charts for the United States: methods and development.Vital Health Stat2461190

18. Villalobos-Comparán M, Flores-Dorantes MT, Villarreal-Molina MT, Rodríguez-Cruz M, García-Ulloa AC, et al. (2008) The FTO gene is associated with adulthood obesity in the Mexican population. Obesity 16: 2296–2301.M. Villalobos-ComparánMT Flores-DorantesMT Villarreal-MolinaM. Rodríguez-CruzAC García-Ulloa2008The FTO gene is associated with adulthood obesity in the Mexican population.Obesity1622962301

19. The Expert Committee on the Diagnosis and Classification of Diabetes Mellitus. Report of the Expert Committee on the Diagnosis and Classification of Diabetes Mellitus (2003) Diabetes Care 26: s5–20.The Expert Committee on the Diagnosis and Classification of Diabetes Mellitus. Report of the Expert Committee on the Diagnosis and Classification of Diabetes Mellitus2003Diabetes Care26s520

20. Villarreal-Molina MT, Aguilar-Salinas CA, Rodríguez-Cruz M, Riaño D, Villalobos-Comparán M, et al. (2007) The ATP-binding cassette transporter A1 R230C variant affects HDL cholesterol levels and BMI in the Mexican population: association with obesity and obesity-related comorbidities. Diabetes 56: 1881–1887.MT Villarreal-MolinaCA Aguilar-SalinasM. Rodríguez-CruzD. RiañoM. Villalobos-Comparán2007The ATP-binding cassette transporter A1 R230C variant affects HDL cholesterol levels and BMI in the Mexican population: association with obesity and obesity-related comorbidities.Diabetes5618811887

21. Levy JC, Matthews DR, Hermans MP (1998) Correct homeostasis model assessment (HOMA) evaluation uses the computer program. Diabetes Care 21: 2191–2192.JC LevyDR MatthewsMP Hermans1998Correct homeostasis model assessment (HOMA) evaluation uses the computer program.Diabetes Care2121912192

22. Folstein MF, Folstein SE, McHugh PR (1975) “Mini-mental state”: a practical method for grading the cognitive state of patients for the clinician. J Psychiatr Res 12: 189–198.MF FolsteinSE FolsteinPR McHugh1975“Mini-mental state”: a practical method for grading the cognitive state of patients for the clinician.J Psychiatr Res12189198

23. Hoggart CJ, Shriver MD, Kittles RA, Clayton DG, McKeigue PM (2004) Design and analysis of admixture mapping studies. Am J Hum Genet 74: 965–978.CJ HoggartMD ShriverRA KittlesDG ClaytonPM McKeigue2004Design and analysis of admixture mapping studies.Am J Hum Genet74965978

24. Canizales-Quinteros S, Aguilar-Salinas CA, Ortiz-López MG, Rodríguez-Cruz M, Villarreal-Molina MT, et al. (2007) Association of PPARG2 Pro12Ala variant with larger body mass index in Mestizo and Amerindian populations of Mexico. Hum Biol 79: 111–119.S. Canizales-QuinterosCA Aguilar-SalinasMG Ortiz-LópezM. Rodríguez-CruzMT Villarreal-Molina2007Association of PPARG2 Pro12Ala variant with larger body mass index in Mestizo and Amerindian populations of Mexico.Hum Biol79111119

25. Creemers JW, Choquet H, Stijnen P, Vatin V, Pigeyre M, et al. (2012) Heterozygous Mutations Causing Partial Prohormone Convertase 1 Deficiency Contribute to Human Obesity. Diabetes 61: 383–390.JW CreemersH. ChoquetP. StijnenV. VatinM. Pigeyre2012Heterozygous Mutations Causing Partial Prohormone Convertase 1 Deficiency Contribute to Human Obesity.Diabetes61383390

26. Gjesing AP, Vestmar MA, Jørgensen T, Heni M, Holst JJ, et al. (2011) The effect of PCKS1 variants on waist, waist-hip ratio and glucose metabolism is modified by sex and glucose tolerance status. PLoS One 6: e23907.AP GjesingMA VestmarT. JørgensenM. HeniJJ Holst2011The effect of PCKS1 variants on waist, waist-hip ratio and glucose metabolism is modified by sex and glucose tolerance status.PLoS One6e23907

27. Heni M, Haupt A, Schäfer SA, Ketterer C, Thamer C, et al. (2010) Association of obesity risk SNPs in PCSK1 with insulin sensitivity and proinsulin conversion. BMC Med Genet 11: 86.M. HeniA. HauptSA SchäferC. KettererC. Thamer2010Association of obesity risk SNPs in PCSK1 with insulin sensitivity and proinsulin conversion.BMC Med Genet1186