You are here

Article Text

Article menu
XML
Genetics/Genomes/Proteomics/Metabolomics
Analysis of the interaction effect of 48 SNPs and obesity on type 2 diabetes in Chinese Hans
  1. Jia Liu1,
  2. Lu Wang1,
  3. Yun Qian1,
  4. Qian Shen1,
  5. Hai Chen1,
  6. Hongxia Ma2,
  7. Juncheng Dai2,
  8. Chong Shen2,
  9. Guangfu Jin2,
  10. Zhibin Hu2,
  11. Hongbing Shen2
  1. 1Department of Health Promotion & Chronic Non-Communicable Disease Control, Wuxi Center for Disease Control and Prevention, Wuxi, China
  2. 2Department of Epidemiology and Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, China
  1. Correspondence to Dr Yun Qian; wxqianyun111{at}sina.com

Abstract

Introduction Both environmental and genetic factors contribute to type 2 diabetes (T2D) risk. Dozens of T2D susceptibility loci have been identified by genome-wide association study. However, these loci account for only a small fraction of the familial T2D risk. We hypothesized that the gene-obesity interaction may contribute to the missing heritability.

Research design and method Forty-eight T2D-associated variants were genotyped using the TaqMan OpenArray Genotyping System and iPLEX Sequenom MassARRAY platform in two separate studies. Obesity was defined according to multiple indexes (body mass index (BMI), waist circumference and waist-hip ratio). Multiplicative interactions were tested using general logistic regression to assess the gene-obesity interaction effect on T2D risk among a total of 6206 Chinese Hans.

Results After adjusting for the main effects of genes and obesity, as well as covariates (age, sex, smoking and alcohol consumption status), robust multiplicative interaction effects were observed between rs10811661 in CDKN2A/CDKN2B and multiple obesity indices (p ranged from 0.001 to 0.043 for BMI, waist circumference and waist-hip ratio). Obese individuals with the TT genotype had a drastically higher risk of T2D than normal weight individuals without the risk allele (OR=17.58, p<0.001). There were no significant differences between subgroups in the stratification analysis. Plausible biological explanations were established using a public database. However, there were no significant interaction effects between the other 47 single nucleotide polymorphism (SNPs) and obesity.

Conclusion Our findings indicated that the CDKN2A/CDKN2B gene-obesity interaction significantly increases T2D risk in Chinese Hans. The interaction effect identified in our study may help to explain some of the missing heritability in the context of T2D susceptibility. In addition, the interaction effect may play a role in the precise prevention of T2D in Chinese individuals.

  • diabetes mellitus
  • type 2
  • genetics
  • obesity
http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

https://doi.org/10.1136/bmjdrc-2020-001638

Statistics from Altmetric.com

    Request Permissions

    If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

    Significance of this study

    What is already known about this subject?

    • The prevalence of type 2 diabetes (T2D) in China rapidly increased in recent years. Several risk factors have been found including obesity and genetic variants.

    What are the new findings?

    • In our current research, we focused on the interaction effect between SNPs and obesity for T2D. We found that CDKN2A/CDKN2B gene-obesity interaction significantly increases T2D risk in Chinese Hans.

    How might these results change the focus of research or clinical practice?

    • The interaction effect identified in our study may help us to explain some of the missing heritability in the context of T2D susceptibility. In addition, the interaction effect may play a role in the precise prevention of T2D in Chinese individuals.

    Introduction

    Type 2 diabetes (T2D) is a common disorder characterized by hyperglycemia, insulin resistance and deficient pancreatic beta-cell function, which accounts for approximately 90% of all diabetes cases.1 According to the estimation of the International Diabetes Federation, 463 million adults worldwide were affected by diabetes in 2019, and this number will increase to 700 million in 2045. In addition, 4.2 million deaths worldwide were attributed to diabetes in 2019. Notably, approximately 79% of patients with diabetes live in low-income and middle-income countries including China.2 The prevalence of T2D in China rapidly increased from less than 1% in 1980 to 11.6% in 2010.3 Diabetes has become an urgent public health problem worldwide, and it is essential to understand its pathogenesis to slow down the rapidly increasing rate.

    Obesity is an underlying cause of chronic noncommunicable disease, including T2D. Over the past few decades, rapid economic changes have drastically affected lifestyle and obesity rates, contributing to an epidemic of T2D.4 Approximately 90% of T2D cases are attributed to excess weight.5 Wilson et al6 found that obesity was a significant predictor of T2D in the prospective Framingham Offspring Study (OR=2.50, p=0.001). These findings together revealed the critical role of obesity in T2D.

    Genetic factors also contribute to the pathogenesis of T2D.1 T2D has an estimated rate of heritability of 30%–70%.7 Over the past several years, genome-wide association studies (GWAS) have provided new susceptibility loci for T2D.8–14 In our previous study, we successfully validated the association between single nucleotide polymorphism (SNPs) and T2D risk as well as lipid levels in Chinese Hans.15 16 However, these common variants account for only a small fraction of the familial risk of T2D, with a modest effect contributed by each locus (OR for each allele ranging from 1.1 to 1.5).17

    To explain the missing heritability, many hypotheses have been proposed, including gene-environment interactions.18 However, no systematic study focusing on multiple genes has been performed to evaluate the interaction between SNPs and obesity in context of T2D risk in Chinese Hans. In this current study, we attempted to evaluate the gene-environment interaction in 2925 T2D cases and 3281 matched controls in Chinese Hans.

    Materials and methods

    Study subjects

    This study was approved by the Institutional Review Board of Wuxi Center for Disease Control and Prevention. The details of the study subjects have been described previously.19 20 In brief, one part of samples (1200 T2D cases and 1200 orthoglycemic controls) were recruited from a community-based noncommunicable disease screening program conducted in Wuxi, Jiangsu, China from April to July in 2007 (Wuxi study). The other portion of samples (1725 cases and 2081 controls) were recruited from a population-based cohort study conducted in Changzhou and Nantong cities in Jiangsu Province during 2004 and 2008 (Changzhou-Nantong study). T2D cases were defined according to fasting plasma glucose (FPG)≥7.0 mmol/L or a history of T2D.21 Controls were selected from the subjects with FPG level <5.6 mmol/L and without a history of diabetes, coronary heart disease, hypertension, stroke or cancer. Controls were frequently matched to the cases by age, sex and residential area. Subjects were interviewed in person by trained interviewers to collect personal information and demographic data after signing informed consent forms. Personal information including name, ID number, address, career, history of diabetes, hypertension, stroke and cancer, smoking and drinking status were collected. Physical examinations including height, weight, waist circumference (WC), hip circumference and blood pressure were performed at the same time. Body mass index (BMI) was divided based on the Chinese criteria or WHO criteria:21 Chinese criteria: underweight: BMI<18.50 kg/m2; normal: 18.50≤BMI<24.00 kg/m2; overweight: 24.00≤BMI<28.00 and obesity: BMI≥28.00 kg/m2. Individuals were defined as smokers if they had smoked at an average of one cigarette or more per day and for at least 1 year in their lifetime. Similarly, participants were defined as drinkers if they drunk at least once per day and ≥4 days per week. Approximately 5 mL venous blood was collected from each subject after fasting for more than 8 hours. FPG and lipid levels (triglycerides (TG), total cholesterol (TC) and high-density lipoprotein cholesterol (HDL-C)) levels were measured in a standard method by Biochemistry Auto-analyzer (Olympus C2734-Au640).

    SNP selection

    The details of the selection process have been described previously.15 Briefly, 76 SNPs were reported in the GWAS catalog, and we included SNPs based on the following criteria: (1) minor allele frequency ≥0.05 in Chinese Han and (2) SNPs with strong linkage disequilibrium (r2>0.8); potential functional variants were retained as priorities. As a result, 48 SNPs from 34 genes were included in this study (online supplemental table 1).

    Supplemental material

    DNA isolation, genotyping and quality control (QC)

    Genomic DNA was isolated from the leukocyte pellets of venous blood by proteinase K digestion, followed by phenol-chloroform extraction and ethanol precipitation. All of the DNA samples were checked for quality by DNA electrophoresis. For the participants in Wuxi study, genotyping was performed using TaqMan OpenArray Genotyping System (Life Technologies, Carlsbad, California, USA). The overall call rates ranged from 98.5% to 99.3%. For Changzhou-Nantong study, SNPs were genotyped by iPLEX Sequenom MassARRAY platform (Sequenom, San Diego, California, USA). For quality control, there were two nontemplate controls in each plate. Based on the duplicate samples, the concordance rate was 100% for this two-genotyping platform.

    Statistical analysis

    The distributions of continuous variables (age, BMI, FPG level, WC, hip circumference and lipid levels) were described using the mean±SD or median value (P25, P75). Categorical variables (sex, smoking status and alcohol consumption status) were defined as counts (percentages). Multiplicative interactions were tested using a general logistic regression model by applying the equation:

    Embedded Image

    where Y is the logit of case-control status, G is the selected SNPs and E is the environment factor (obesity). β0 is the constant, and βG and βE represent the main effects of SNPs and obesity, respectively. βGE is the interaction term and Covari are the covariates for adjustment, including age, sex, smoking status and alcohol consumption status. The effect size directions of interaction should be consistent in Wuxi study and Changzhou-Nantong study for further analysis. To examine the differences between subgroups, the χ2-based Q-test was used to test the heterogeneity of effect sizes (ORs and 95% CIs) derived from the corresponding subgroups. All of the statistical analyses were performed with R software (V.3.5.1; The R Foundation for Statistical Computing, http://www.cran.r-project.org/) and Stata V.12.0 software (Stata, College Station, Texas, USA).

    Results

    The demographic and clinical characteristics of study population are summarized in table 1. No significant differences of sex, smoking status and alcohol consumption status were observed in either group or the combined population (p>0.05). The age of the T2D case group was approximately 1 year older than that of the control group, although the age group (5 years) was matched for the case and control groups. As expected, patients with T2D had significantly higher obesity indexes (BMI, WC, hip circumference), FPG, TGs and TC and lower levels of HDL-C than the control group in the combined population and the two separate groups.

    View this table:
    Table 1

    Demographic and clinical characteristics of study subjects

    The associations results between selected SNPs, obesity indices and T2D risk are shown in online supplemental table 2. First, we focused on the effect of the interaction between selected SNPs and BMI level on T2D risk. Of the 48 selected SNPs, multiplicative interaction analysis revealed that the interaction effect directions of seven SNPs were consistent in Wuxi study and Changzhou-Nantong study. As shown in table 2, of these seven SNPs, rs10811661 in CDKN2A/CDKN2B had a significant interaction effect with BMI level on T2D risk (p<0.05). Of interest, effects of the interactions between rs10811661 and BMI level were significant in both studies (Wuxi study: OR=1.22, p=0.033; Changzhou-Nantong study: OR=1.21, p=0.035). After false discovery rate (FDR) correction, rs10811661 still had a significant interaction effect with BMI level on T2D risk in the combined population (OR=1.22, p=0.021). Besides, in consideration of the small number of obesity, we combined the overweight and obesity group to conduct a sensitive analysis. The sensitive analysis showed that the interaction effect was robust (OR=1.27, p=0.002). However, there were no significant interaction effects between the other SNPs and BMI on T2D risk (online supplemental table 1).

    View this table:
    Table 2

    Interaction between selected SNPs and BMI levels on T2D risk

    Compared with normal weight individuals without risk allele T, obese individuals with the TT genotype had a substantially increased risk of T2D (OR=17.58, p<0.001) (table 3).

    View this table:
    Table 3

    Interaction between rs10811661 and BMI level on T2D risk

    We further analyzed the interaction effect between rs10811661 and other obesity indices (WC and waist-hip ratio (WHR)) on T2D risk. As shown in online supplemental tables 3 and 4, significant interaction effects were also observed between rs10811661 and other obesity indices (for WC, p=0.001; for WHR, p=0.043).

    In the stratification analysis, the interaction effect was evaluated in subgroups based on age, sex, smoking status and alcohol consumption status. As shown in figure 1, there were no significant differences between subgroups in regard to the interaction effect between rs10811661 and BMI level (all P for heterogeneity >0.05).

    Figure 1
    Figure 1

    Stratified analysis of the interaction between rs10811661 and BMI on type 2 diabetes risk. BMI, body mass index.

    To reveal the biological role of the identified SNP, we performed a bioinformatic analysis of the results above by using the public database GTEX (http://www.gtexportal.org/home/). According to GTEX database, rs10811661 in CDKN2A/CDKN2B significantly decreased the expression of CDKN2B in whole blood samples (β=−0.12, p=0.016) and pancreas (β=−0.17, p=0.035) (figure 2). However, there was no significant association between rs10811661 and CDKN2A (data not shown).

    Figure 2
    Figure 2

    The association of rs10811661 with host gene expression. Data were available in public database GTEX (A: whole blood; B: pancreas).

    Discussion

    In this study, we evaluated the effect of the interaction between 48 selected SNPs and obesity indexes on T2D risk in Han Chinese individuals. Robust interaction effects were observed between rs10811661 and BMI (Chinese criterion and WHO criterion), WC and WHR. To the best of our knowledge, this is the first study to report the effects of the interaction between rs10811661 in CDKN2A/CDKN2B and different obesity indexes in Han Chinese individuals. Our study provided genetic evidence at the population level that the gene-obesity interaction may contribute to T2D risk.

    Several studies have focused on the interaction effects of SNPs and obesity on T2D risk.22–24 However, these studies each mainly focused on only one gene (PPARG, FOXO1, TCF7L2), and the number of selected SNPs was limited. In addition, no research has been reported on the CDKN2A/CDKN2B-obesity interaction.

    The SNP rs10811661 is located upstream of CDKN2A/CDKN2B (cyclin-dependent kinase inhibitor 2A/2B). The effect of rs10811661 on T2D risk is well known. Li et al25 performed a meta-analysis including 17 studies with 29 990 cases and 40 977 controls to evaluate the association between rs10811661 and T2D risk. The results showed that the rs10811661 T allele was significantly associated with increased T2D risk in the additive model (OR=1.51, 95% CI 1.40 to 1.63, p<0.001). In our previous study, a similar result was identified (OR=1.32, 95% CI 1.17 to 1.49, p=1.22×10−5).15 Interestingly, in human populations, the rs10811661 risk allele T is also associated with reduced insulin secretory capacity.26–28

    According to the GTEX database, rs10811661 is an eQTL (expression quantitative trait loci) for CDKN2B. The risk allele T was identified to be associated with increased CDKN2B levels in whole blood and pancreas samples. Proteins encoded by the CDKN2B gene regulate pancreatic beta cell replication by inhibiting cyclin-dependent kinase 4/6 (CDK4/6).29 It has been shown that mice lacking CDK4 exhibit insulin-deficient diabetes due to a reduction in pancreatic beta cells.30 These results together suggest that rs10811661 may increase the risk of T2D by affecting pancreatic beta cells. Recently, Kong et al31 found that the risk allele of rs10811661 increased ANRIL expression in islet samples, which correlates with higher BMI and fat mass.32–35 Obesity is well studied and has been shown to contribute to T2D risk mainly by reducing insulin sensitivity and promoting insulin resistance.36 These mechanisms may together elucidate the increased risk of T2D in obese individuals with the rs10811661 TT genotype.

    BMI is usually used to assess general obesity. However, general obesity, as defined by BMI, is not a fitness indicator of body fat distribution.37 Recently, central obesity, also known as abdominal obesity, measured by WC and WHR have been widely studied.38 Cheng et al39 found that WHR was a better index than BMI for predicting the risk of T2D in Taiwanese population. Besides, Hu et al40 conducted the interASIA study and found that central obesity is more related to T2D than overall obesity in the Chinese population. However, in our study, the interaction effects between SNP and general obesity, central obesity were similar (ORs for BMI, WC and WHR were 1.22, 1.31 and 1.20, respectively). There was no significant difference between these obesity indices (p for heterogeneity was 0.84). Li et al24 conducted a study to explore the interaction effect of TCF7L2 and BMI/ WC on T2D. The effect was similar for BMI and WC (for BMI: OR=4.833; WC: OR=4.629). The present study is not sufficiently large to support the hypothesis that there may be difference of interaction effect between SNP-general obesity and SNP-central obesity. Further research, which focus on the difference, should be done.

    There were several strengths in our present study. First, this study was well designed and included two separate study populations to ensure the authenticity of the results. Second, multiple obesity indexes were evaluated for interaction effects to ensure the robustness of the study results. Third, a public database was used, and plausible interpretations of our population results were established. However, the biological mechanisms of the identified interaction effect remain unclear. Future functional studies will be necessary to validate our findings.

    In summary, our results suggested, for the first time, that rs10811661 in CDKN2A/CDKN2B had robust interaction effects with multiple obesity indexes. The identified gene-obesity interaction may contribute to some of the missing heritability of T2D risk.

    Acknowledgments

    The authors wish to thank all of the study participants, research staff and students who participated in this work. We thank the GTEX Projects for making the data public.

    References

    1. ↵
      1. Stumvoll M,
      2. Goldstein BJ,
      3. van Haeften TW
      . Type 2 diabetes: principles of pathogenesis and therapy. Lancet 2005;365:1333–46.doi:10.1016/S0140-6736(05)61032-Xpmid:http://www.ncbi.nlm.nih.gov/pubmed/15823385
      OpenUrlCrossRefPubMedWeb of Science
    2. ↵
      1. The International Diabetes Federation
      . Available: http://www.diabetesatlas.org/
    3. ↵
      1. Xu Y,
      2. Wang L,
      3. He J, et al
      . Prevalence and control of diabetes in Chinese adults. JAMA 2013;310:948–59.doi:10.1001/jama.2013.168118pmid:http://www.ncbi.nlm.nih.gov/pubmed/24002281
      OpenUrlCrossRefPubMedWeb of Science
    4. ↵
      1. Al-Daghri NM,
      2. Al-Attas OS,
      3. Alokail MS, et al
      . Diabetes mellitus type 2 and other chronic non-communicable diseases in the central region, Saudi Arabia (Riyadh cohort 2): a decade of an epidemic. BMC Med 2011;9:76. doi:10.1186/1741-7015-9-76pmid:http://www.ncbi.nlm.nih.gov/pubmed/21689399
      OpenUrlCrossRefPubMed
    5. ↵
      1. Hossain P,
      2. Kawar B,
      3. El Nahas M
      . Obesity and diabetes in the developing world--a growing challenge. N Engl J Med 2007;356:213–5.doi:10.1056/NEJMp068177pmid:http://www.ncbi.nlm.nih.gov/pubmed/17229948
      OpenUrlCrossRefPubMedWeb of Science
    6. ↵
      1. Wilson PWF,
      2. Meigs JB,
      3. Sullivan L, et al
      . Prediction of incident diabetes mellitus in middle-aged adults: the Framingham offspring study. Arch Intern Med 2007;167:1068–74.doi:10.1001/archinte.167.10.1068pmid:http://www.ncbi.nlm.nih.gov/pubmed/17533210
      OpenUrlCrossRefPubMedWeb of Science
    7. ↵
      1. Almgren P,
      2. Lehtovirta M,
      3. Isomaa B, et al
      . Heritability and familiality of type 2 diabetes and related quantitative traits in the Botnia study. Diabetologia 2011;54:2811–9.doi:10.1007/s00125-011-2267-5pmid:http://www.ncbi.nlm.nih.gov/pubmed/21826484
      OpenUrlCrossRefPubMedWeb of Science
    8. ↵
      1. Diabetes Genetics Initiative of Broad Institute of Harvard and MIT, Lund University, and Novartis Institutes of BioMedical Research,
      2. Saxena R,
      3. Voight BF, et al
      . Genome-Wide association analysis identifies loci for type 2 diabetes and triglyceride levels. Science 2007;316:1331–6.doi:10.1126/science.1142358pmid:http://www.ncbi.nlm.nih.gov/pubmed/17463246
      OpenUrlAbstract/FREE Full Text
    9. ↵
      1. Zeggini E,
      2. Weedon MN,
      3. Lindgren CM, et al
      . Replication of genome-wide association signals in UK samples reveals risk loci for type 2 diabetes. Science 2007;316:1336–41.doi:10.1126/science.1142364pmid:http://www.ncbi.nlm.nih.gov/pubmed/17463249
      OpenUrlAbstract/FREE Full Text
    10. ↵
      1. Scott LJ,
      2. Mohlke KL,
      3. Bonnycastle LL, et al
      . A genome-wide association study of type 2 diabetes in Finns detects multiple susceptibility variants. Science 2007;316:1341–5.doi:10.1126/science.1142382pmid:http://www.ncbi.nlm.nih.gov/pubmed/17463248
      OpenUrlAbstract/FREE Full Text
    11. ↵
      1. Sladek R,
      2. Rocheleau G,
      3. Rung J, et al
      . A genome-wide association study identifies novel risk loci for type 2 diabetes. Nature 2007;445:881–5.doi:10.1038/nature05616pmid:http://www.ncbi.nlm.nih.gov/pubmed/17293876
      OpenUrlCrossRefPubMedWeb of Science
    12. ↵
      1. Zeggini E,
      2. Scott LJ,
      3. Saxena R, et al
      . Meta-analysis of genome-wide association data and large-scale replication identifies additional susceptibility loci for type 2 diabetes. Nat Genet 2008;40:638–45.doi:10.1038/ng.120pmid:http://www.ncbi.nlm.nih.gov/pubmed/18372903
      OpenUrlCrossRefPubMedWeb of Science
    13. ↵
      1. Dupuis J,
      2. Langenberg C,
      3. Prokopenko I, et al
      . New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk. Nat Genet 2010;42:105–16.doi:10.1038/ng.520pmid:http://www.ncbi.nlm.nih.gov/pubmed/20081858
      OpenUrlCrossRefPubMedWeb of Science
    14. ↵
      1. Voight BF,
      2. Scott LJ,
      3. Steinthorsdottir V, et al
      . Twelve type 2 diabetes susceptibility loci identified through large-scale association analysis. Nat Genet 2010;42:579–89.doi:10.1038/ng.609pmid:http://www.ncbi.nlm.nih.gov/pubmed/20581827
      OpenUrlCrossRefPubMedWeb of Science
    15. ↵
      1. Qian Y,
      2. Lu F,
      3. Dong M, et al
      . Cumulative effect and predictive value of genetic variants associated with type 2 diabetes in Han Chinese: a case-control study. PLoS One 2015;10:e0116537. doi:10.1371/journal.pone.0116537pmid:http://www.ncbi.nlm.nih.gov/pubmed/25587982
      OpenUrlPubMed
    16. ↵
      1. Liu J,
      2. Wang L,
      3. Qian Y, et al
      . Association of 48 type 2 diabetes susceptibility loci with fasting plasma glucose and lipid levels in Chinese Hans. Diabetes Res Clin Pract 2018;139:114–21.doi:10.1016/j.diabres.2018.02.039pmid:http://www.ncbi.nlm.nih.gov/pubmed/29518490
      OpenUrlPubMed
    17. ↵
      1. Manolio TA,
      2. Collins FS,
      3. Cox NJ, et al
      . Finding the missing heritability of complex diseases. Nature 2009;461:747–53.doi:10.1038/nature08494pmid:http://www.ncbi.nlm.nih.gov/pubmed/19812666
      OpenUrlCrossRefPubMedWeb of Science
    18. ↵
      1. Maher B
      . Personal genomes: the case of the missing heritability. Nature 2008;456:18–21.doi:10.1038/456018apmid:http://www.ncbi.nlm.nih.gov/pubmed/18987709
      OpenUrlCrossRefPubMedWeb of Science
    19. ↵
      1. Lu F,
      2. Qian Y,
      3. Li H, et al
      . Genetic variants on chromosome 6p21.1 and 6p22.3 are associated with type 2 diabetes risk: a case-control study in Han Chinese. J Hum Genet 2012;57:320–5.doi:10.1038/jhg.2012.25pmid:http://www.ncbi.nlm.nih.gov/pubmed/22437209
      OpenUrlCrossRefPubMed
    20. ↵
      1. Qian Y,
      2. Lu F,
      3. Dong M, et al
      . Genetic variants of IDE-KIF11-HHEX at 10q23.33 associated with type 2 diabetes risk: a fine-mapping study in Chinese population. PLoS One 2012;7:e35060. doi:10.1371/journal.pone.0035060pmid:http://www.ncbi.nlm.nih.gov/pubmed/22506066
      OpenUrlPubMed
    21. ↵
      1. Expert Committee on the Diagnosis and Classification of Diabetes Mellitus
      . Report of the expert Committee on the diagnosis and classification of diabetes mellitus. Diabetes Care 2000;23 Suppl 1:S4–19.pmid:http://www.ncbi.nlm.nih.gov/pubmed/12017675
      OpenUrlPubMed
    22. ↵
      1. Lv X,
      2. Zhang L,
      3. Sun J, et al
      . Interaction between peroxisome proliferator-activated receptor gamma polymorphism and obesity on type 2 diabetes in a Chinese Han population. Diabetol Metab Syndr 2017;9:7.doi:10.1186/s13098-017-0205-5pmid:http://www.ncbi.nlm.nih.gov/pubmed/28123453
      OpenUrlPubMed
    23. ↵
      1. Gong L,
      2. Li R,
      3. Ren W, et al
      . The FoxO1 Gene-Obesity interaction increases the risk of type 2 diabetes mellitus in a Chinese Han population. J Korean Med Sci 2017;32:264–71.doi:10.3346/jkms.2017.32.2.264pmid:http://www.ncbi.nlm.nih.gov/pubmed/28049237
      OpenUrlPubMed
    24. ↵
      1. Li L,
      2. Wang J,
      3. Ping Z, et al
      . Interaction analysis of gene variants of TCF7L2 and body mass index and waist circumference on type 2 diabetes. Clin Nutr 2020;39:192–7.doi:10.1016/j.clnu.2019.01.014pmid:http://www.ncbi.nlm.nih.gov/pubmed/30718095
      OpenUrlPubMed
    25. ↵
      1. Li H,
      2. Tang X,
      3. Liu Q, et al
      . Association between type 2 diabetes and rs10811661 polymorphism upstream of CDKN2A/B: a meta-analysis. Acta Diabetol 2013;50:657–62.doi:10.1007/s00592-012-0400-7pmid:http://www.ncbi.nlm.nih.gov/pubmed/22623142
      OpenUrlPubMed
    26. ↵
      1. Grarup N,
      2. Rose CS,
      3. Andersson EA, et al
      . Studies of association of variants near the Hhex, CDKN2A/B, and IGF2BP2 genes with type 2 diabetes and impaired insulin release in 10,705 Danish subjects: validation and extension of genome-wide association studies. Diabetes 2007;56:3105–11.doi:10.2337/db07-0856pmid:http://www.ncbi.nlm.nih.gov/pubmed/17827400
      OpenUrlAbstract/FREE Full Text
    27. ↵
      1. Hribal ML,
      2. Presta I,
      3. Procopio T, et al
      . Glucose tolerance, insulin sensitivity and insulin release in European non-diabetic carriers of a polymorphism upstream of CDKN2A and CDKN2B. Diabetologia 2011;54:795–802.doi:10.1007/s00125-010-2038-8pmid:http://www.ncbi.nlm.nih.gov/pubmed/21234743
      OpenUrlPubMed
    28. ↵
      1. 't Hart LM,
      2. Simonis-Bik AM,
      3. Nijpels G, et al
      . Combined risk allele score of eight type 2 diabetes genes is associated with reduced first-phase glucose-stimulated insulin secretion during hyperglycemic clamps. Diabetes 2010;59:287–92.doi:10.2337/db09-0736pmid:http://www.ncbi.nlm.nih.gov/pubmed/19808892
      OpenUrlAbstract/FREE Full Text
    29. ↵
      1. Moritani M,
      2. Yamasaki S,
      3. Kagami M, et al
      . Hypoplasia of endocrine and exocrine pancreas in homozygous transgenic TGF-beta1. Mol Cell Endocrinol 2005;229:175–84.doi:10.1016/j.mce.2004.08.007pmid:http://www.ncbi.nlm.nih.gov/pubmed/15607541
      OpenUrlCrossRefPubMedWeb of Science
    30. ↵
      1. Rane SG,
      2. Dubus P,
      3. Mettus RV, et al
      . Loss of CDK4 expression causes insulin-deficient diabetes and CDK4 activation results in beta-islet cell hyperplasia. Nat Genet 1999;22:44–52.doi:10.1038/8751pmid:http://www.ncbi.nlm.nih.gov/pubmed/10319860
      OpenUrlCrossRefPubMedWeb of Science
    31. ↵
      1. Kong Y,
      2. Sharma RB,
      3. Ly S, et al
      . CDKN2A/B T2D Genome-Wide Association Study Risk SNPs Impact Locus Gene Expression and Proliferation in Human Islets. Diabetes 2018;67:872–84.doi:10.2337/db17-1055pmid:http://www.ncbi.nlm.nih.gov/pubmed/29432124
      OpenUrlAbstract/FREE Full Text
    32. ↵
      1. Lillycrop K,
      2. Murray R,
      3. Cheong C, et al
      . Anril promoter DNA methylation: a perinatal marker for later adiposity. EBioMedicine 2017;19:60–72.doi:10.1016/j.ebiom.2017.03.037pmid:http://www.ncbi.nlm.nih.gov/pubmed/28473239
      OpenUrlPubMed
    33. ↵
      1. Curtis EM,
      2. Murray R,
      3. Titcombe P, et al
      . Perinatal DNA methylation at CDKN2A is associated with offspring bone mass: findings from the Southampton women's survey. J Bone Miner Res 2017;32:2030–40.doi:10.1002/jbmr.3153pmid:http://www.ncbi.nlm.nih.gov/pubmed/28419547
      OpenUrlPubMed
    34. ↵
      1. Svensson P-A,
      2. Wahlstrand B,
      3. Olsson M, et al
      . Cdkn2B expression and subcutaneous adipose tissue expandability: possible influence of the 9p21 atherosclerosis locus. Biochem Biophys Res Commun 2014;446:1126–31.doi:10.1016/j.bbrc.2014.03.075pmid:http://www.ncbi.nlm.nih.gov/pubmed/24680834
      OpenUrlCrossRefPubMed
    35. ↵
      1. Horswell SD,
      2. Fryer LGD,
      3. Hutchison CE, et al
      . Cdkn2B expression in adipose tissue of familial combined hyperlipidemia patients. J Lipid Res 2013;54:3491–505.doi:10.1194/jlr.M041814pmid:http://www.ncbi.nlm.nih.gov/pubmed/24103848
      OpenUrlAbstract/FREE Full Text
    36. ↵
      1. Li P,
      2. Oh DY,
      3. Bandyopadhyay G, et al
      . Ltb4 promotes insulin resistance in obese mice by acting on macrophages, hepatocytes and myocytes. Nat Med 2015;21:239–47.doi:10.1038/nm.3800pmid:http://www.ncbi.nlm.nih.gov/pubmed/25706874
      OpenUrlCrossRefPubMed
    37. ↵
      1. Rankinen T,
      2. Kim SY,
      3. Pérusse L, et al
      . The prediction of abdominal visceral fat level from body composition and anthropometry: ROC analysis. Int J Obes Relat Metab Disord 1999;23:801–9.doi:10.1038/sj.ijo.0800929pmid:http://www.ncbi.nlm.nih.gov/pubmed/10490780
      OpenUrlCrossRefPubMedWeb of Science
    38. ↵
      1. Zhang X,
      2. Shu X-O,
      3. Gao Y-T, et al
      . General and abdominal adiposity and risk of stroke in Chinese women. Stroke 2009;40:1098–104.doi:10.1161/STROKEAHA.108.539692pmid:http://www.ncbi.nlm.nih.gov/pubmed/19211490
      OpenUrlAbstract/FREE Full Text
    39. ↵
      1. Cheng C-H,
      2. Ho C-C,
      3. Yang C-F, et al
      . Waist-to-hip ratio is a better anthropometric index than body mass index for predicting the risk of type 2 diabetes in Taiwanese population. Nutr Res 2010;30:585–93.doi:10.1016/j.nutres.2010.08.007pmid:http://www.ncbi.nlm.nih.gov/pubmed/20934599
      OpenUrlCrossRefPubMed
    40. ↵
      1. Hu D,
      2. Xie J,
      3. Fu P, et al
      . Central rather than overall obesity is related to diabetes in the Chinese population: the InterASIA study. Obesity 2007;15:2809–16.doi:10.1038/oby.2007.333pmid:http://www.ncbi.nlm.nih.gov/pubmed/18070772
      OpenUrlPubMed

    Supplementary materials

    • Supplementary Data

      This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

    Footnotes

    • JL and LW are joint first authors.

    • Contributors YQ designed and directed this study. JL and LW conducted the genotype assay and wrote this paper. QS, HM and CS directed the survey and quality control. HC and JD analyzed the data. GJ, ZH and HS directed the assay and quality control.

    • Funding This work was funded by the Wuxi Key Medical Disciplines grant (ZDXK010), the Wuxi Key Medical Talents grant (ZDRC004), Wuxi Young Medical Talents grant (QNRC050) for ‘Strengthening Health though Science and Education’, the important project of Wuxi Health Commission (Z201905), the Project of National Natural Science Foundation of China (81072379), the Jiangsu scientific research project of top health talents (LGY2018014), the Wuxi science and technology development fund project (WX18IIAN038) and the Wuxi Health Commission Youth Foundation (Q201802).

    • Competing interests None declared.

    • Patient consent for publication Not required.

    • Ethics approval This study was approved by the Institutional Review Board of Wuxi Center for Disease Control and Prevention (approval ID: 2011(001)).

    • Provenance and peer review Not commissioned; externally peer reviewed.

    • Data availability statement Data are available on reasonable request. Data is available upon reasonable request from corresponding author.

    • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.