Diabetes mellitus (DM) is a complex metabolic disease that is caused by a complex interplay between genetic and environmental factors. This research aimed to investigate the association of genetic polymorphisms in PDX1 and MC4R with T2DM risk.
The genotypes of 10 selected SNPs in PDX1 and MC4R were identified using the Agena MassARRAY platform. We utilized odds ratio (OR) and 95% confidence intervals (CIs) to assess the correlation between genetic polymorphisms and T2DM risk.
We found that PDX1-rs9581943 decreased susceptibility to T2DM among in a Chinese Han population (OR = 0.76, p = 0.045). We also found that selected genetic polymorphisms in PDX1 and MC4R could modify the risk of T2DM, which might also be influenced by age, sex, BMI, smoking status, and drinking status (p < 0.05).
We concluded that PDX1 and MC4R genetic variants were significantly associated with T2DM risk in a Chinese Han population. These single polymorphic markers may be considered to be new targets in the assessment and prevention of T2DM among Chinese Han people.
Diabetes mellitus (DM) is a metabolic disease characterized by the presence of chronic hyperglycemia, which results from either weakened insulin secretion or insulin action or both . The global prevalence of diabetes reached 9.3% (463 million) in 2019, and it is expected to increase to 10.9% (700 million) by 2045 . China has the highest number of adults with diabetes, approximatedly116 million, ranking first in diabetes prevalence worldwide . Type 2 diabetes mellitus (T2DM) accounts for nearly 90% of the total diabetes patients. There are multiple reasons for the incidence of T2DM including aging, sedentary lifestyles and genetic factors . It has been reported that subjects withT2DM-affected siblings have a two- to three fold increased risk of developing T2DM compared with the general population . Having one parent with diabetes increases the risk of developing T2DM by 30–40%, and having two parents with diabetes increases the risk to 70% . Furthermore, some research reported that genetic polymorphisms in candidate genes could influence the formation and course of T2DM [6, 7].
Pancreatic and duodenal homeobox-1 (PDX1) modulates pancreas development and β-cell function. The PDX1 gene encodes a protein of 283 amino acids in humans. It also regulates many genes, such as those encoding insulin and glucokinase (GK), involved in maintaining the function of β-cells. In adults, PDX1 is highly expressed in β-cells, where it is required for efficient insulin gene transcription . Indeed, PDX1 has been proposed to be an oncogene, since its overexpression increased pancreatic cancer cell proliferation, invasion, and growth in humans . Gurevich et al. also illustrated that PDX1 was upregulated in neuroendocrine tumors, including pancreatic ductal and acinar cell tumors and gastric signet ring cell carcinomas . It has previously been noted that PDX1 deficiency inhibits the development of pancreatic buds, leading to extreme hyperglycemia . These findings demonstrated that PDX1 plays a pivotal role in the development of pancreas-related disease. However, no literature supports the effect of PDX1 polymorphisms on T2DM.
Melanocortin receptor 4 (MC4R) belongs to class A of G protein-coupled receptors and is a member of the melanocortin receptor family . MC4R can control energy homeostasis, sympathetic nervous system activity, and blood pressure in rodents and humans . For instance, MC4R knockdown mice were severely obese and the loss of one MC4R allele resulted in an intermediate obesity phenotype . Greenfield et al. demonstrated reduction in blood pressure and circulating catecholamine levels in humans with MC4R deficiency . In addition, previous research has established that MC4R deletion or mutation results in obesity, hyperphagia, and insulin resistance . These observations highlight a potential role for MC4R in obesity-related diseases. In addition, obesity is believed to be an independent risk factor for T2DM . Based on the above information, we hypothesized that MC4R may be involved in the occurrence of T2DM.
Therefore, we mainly examined the role of PDX1 and MC4R genetic polymorphisms in T2DM development in a Chinese population. We identified four polymorphisms in PDX1 (rs11619319, rs2293941, rs9581943 and rs7981781) and six polymorphisms in MC4R (rs6567160, rs663129, rs17782313, rs12969709, rs11663816, and rs12970134) to investigate the correlations between genetic polymorphisms and T2DM susceptibility. The current study will provide new targets for the early assessment and prevention of T2DM.
A total of 500 T2DM patients and 501 healthy controls were enrolled from the First Affiliated Hospital of Xi’an Jiaotong University in the present study. All patients were diagnosed with T2DM based on fasting plasma glucose ≥ 7.0 mmol/L or postprandial plasma glucose ≥ 11.1 mmol/L or HbA1c ≥ 6.5% . Patients with type 1 diabetes mellitus; gestational diabetes; acute or chronic diseases of the liver, kidney, or heart; other endocrine diseases; inflammatory diseases; or malignant tumors were excluded. The inclusion criteria for controls were no history of diabetes, metabolic disorders or severe diseases. The demographic and clinical characteristics of all subjects, including age, sex, smoking status, drinking status, complications, and body mass index (BMI), were collected from medical records and questionnaires.
This research received approval from the Ethics Committee of the First Affiliated Hospital of Xi’an Jiaotong University, and conformed to the Declaration of Helsinki. Informed consent was acquired from each participant at recruitment after fully describing our research to them.
We selected four SNPs in PDX1 and six SNPs in MC4R and all SNPs had minor allele frequencies (MAFs) ≥ 5% in the 1000 Genomes Chinese Han Beijing population. Peripheral blood samples (5 mL) were collected from each subject, and genomic DNA was extracted using the GoldMag whole-blood DNA purification kit (GoldMag Co.Ltd., Xi’an, China) following the manufacturer’s protocol. Genotyping of PDX1 and MC4R polymorphisms was performed by the Agena MassARRAY platform (Agena Bioscience, San Diego, CA, USA). Moreover, Agena Typer 4.0 software was used to analyze and manage data.
Gene expression analysis
We performed PDX1 and MC4R mRNA expression analysis with blood samples from 100 unrelated Chinese Han individuals. Total RNA was isolated from peripheral blood using a Qiagen kit (Qiagen) according to the manufacturer’s instructions. RNA was reverse transcribed to synthesize first-strand cDNA using the PrimeScript-1st strand cDNA Synthesis Kit (Takara Bio, Shiga, Japan), as described by the manufacturers. The mRNA expression of the PDX1 and MC4R genes and the internal control GAPDH were assessed using quantitative real-time PCR (ABI PRISM 7500 Real-Time PCR System; Applied Biosystems). The relative mRNA expression was calculated by the 2−Δ(ΔCt) comparative method and normalized to GAPDH expression.
The primer sequences for the mRNA expression of PDX1, MC4R and GAPDH are shown in Additional file 1: Table S1. Amplification was performed in a reaction mixture containing 10 pM each primer, 10 μl SYBR Green/High ROX (Amplicon), 7 μl nuclease-free water, and 2 μl cDNA solution. Experiments were performed in triplicate.
Statistical differences in demographic characteristics of the participants were assessed using the χ2 test and Student′s t-test. Hardy–Weinberg equilibrium (HWE) of each SNP among controls was evaluated using the χ2 test. The association of the selected SNPs with T2DM susceptibility was examined by odds ratio (ORs) and 95% confidence intervals (CIs) by logistic regression analysis in multiple inheritance models and different subgroups (age, sex, smoking, drinking and BMI). The potential functions of the selected SNPs were forecasted using HaploReg v4.1 (https://pubs.broadinstitute.org/mammals/haploreg/haploreg.php). Haploview software and PLINK software were used for Haploview analysis and linkage disequilibrium [19, 20]. The mRNA expression was analyzed using Student’s t-test in the case and control groups. The effects of the polymorphisms on mRNA expression were examined by one-way analysis of variance (ANOVA). A p value < 0.05 was considered statistically significant.
Characteristics of the study population
As presented in Table 1, there were 500 T2DM patients (358 men and 142 women) and 501 healthy controls (358 men and 143 women) in this study. The average ages were 59.87 ± 12.87 years for cases and 59.85 ± 9.34 years for controls. There were no significant differences in age (p = 0.973) or sex (p = 0.960) between the case and control groups. In addition, significant differences were observed in total cholesterol (p < 0.001), low-density lipoprotein cholesterol (LDL-C, p = 0.012), high-density lipoprotein cholesterol (HDL-C, p = 0.024), fasting blood glucose (p < 0.001) and urea (p < 0.001) between the two groups.
T2DM risk assessment
Four candidate SNPs in PDX1 (rs11619319, rs2293941, rs9581943, and rs7981781) and six SNPs in MC4R (rs6567160, rs663129, rs17782313, rs12969709, rs11663816, and rs12970134) were successfully genotyped, as shown in Additional file 1: Table S2. Deviation from HWE was assessed in controls and all candidate SNPs reached the expected p values (p > 0.05). There were no significant associations between allele frequencies of any SNP and susceptibility to T2DM (p > 0.05).
Additionally, we investigated the correlation of PDX1 and MC4R polymorphisms with T2DM risk in multiple inheritance models by logistic regression analyses (Table 2). The results revealed that the AG genotype of PDX1-rs9581943 decreased susceptibility to T2DM in the study subjects (OR = 0.76, 95% CI = 0.58–0.99, p = 0.045).
Stratification analysis was carried out by age, sex, smoking, drinking and BMI. The results of stratification by age and sex are shown in Table 3. We found that PDX1-rs9581943 significantly decreased the risk of T2DM among patients aged ≤ 60 years in the codominant (OR = 0.66, 95% CI = 0.45–0.98, p = 0.039) and dominant models (OR = 0.69, 95% CI = 0.48–1.00, p = 0.049). Rs6567160, rs663129, rs17782313, rs12969709 and rs11663816 in MC4R reduced the susceptibility to T2DM among individuals aged ≤ 60 years under the codominant (rs6567160: OR = 0.33, 95% CI = 0.13–0.81, p = 0.015; rs663129: OR = 0.33, 95% CI = 0.13–0.82, p = 0.017; rs17782313: OR = 0.34, 95% CI = 0.14–0.83, p = 0.018; rs12969709: OR = 0.27, 95% CI = 0.10–0.75, p = 0.012; rs11663816: OR = 0.31, 95% CI = 0.11–0.88, p = 0.027) and recessive (rs6567160: OR = 0.33, 95% CI = 0.14–0.81, p = 0.016; rs663129: OR = 0.33, 95% CI = 0.14–0.81, p = 0.016; rs17782313: OR = 0.33, 95% CI = 0.14–0.81, p = 0.016; rs12969709: OR = 0.27, 95% CI = 0.10–0.75, p = 0.012; and rs11663816: OR = 0.32, 95% CI = 0.11–0.91, p = 0.032) models. After stratifying by sex, rs9581943 (OR = 0.73, 95% CI = 0.5–1.00, p = 0.049) and rs7981781 (OR = 0.70, 95% CI = 0.56–0.97, p = 0.033) were found to be associated with a decreased risk of T2DM in males under the codominant model.
In addition, as shown in Table 4, PDX1-rs7981781 reduced the susceptibility to T2DM among smokers under the codominant (OR = 0.50, 95% CI = 0.29–0.89, p = 0.018) and dominant (OR = 0.55, 95% CI = 0.32–0.95, p = 0.030) models. However, MC4R-rs6567160 could increase the occurrence of T2DM among nonsmokers under the codominant (OR = 1.60, 95% CI = 1.04–2.45, p = 0.032) and dominant (OR = 1.56, 95% CI = 1.04–2.34, p = 0.031) models. MC4R-rs663129 induced a significantly higher susceptibility to T2DM among individuals who were nonsmokers in the codominant (OR = 1.64, 95% CI = 1.07–2.52, p = 0.023), dominant (OR = 1.60, 95% CI = 1.07–2.40, p = 0.023) and additive (OR = 1.40, 95% CI = 1.00–1.95, p = 0.049) models. Moreover, rs17782313 in MC4R was related to a higher risk of T2DM among nonsmokers under the allelic (OR = 1.43, 95% CI = 1.00–1.95, p = 0.036), codominant (OR = 1.72, 95% CI = 1.12–2.64, p = 0.014), dominant (OR = 1.66, 95% CI = 1.11–2.50, p = 0.014) and additive (OR = 1.44, 95% CI = 1.03–2.01, p = 0.034) models.
Additionally, PDX1-rs11619319 (OR = 0.51, 95% CI = 0.27–20.97, p = 0.039) and rs2293941 (OR = 0.51, 95% CI = 0.27–0.97, p = 0.040) were predominantly related to a reduced risk of T2DM among drinkers under the codominant model. Rs7981781 was correlated with a lower risk of T2DM among drinkers under the codominant (OR = 0.47, 95% CI = 0.25–0.88, p = 0.019) and dominant (OR = 0.49, 95% CI = 0.27–0.90, p = 0.022) models.
When stratified by BMI (Table 5), PDX1-rs7981781 was correlated with a lower risk of T2DM among subjects with BMI > 24 kg/m2 under the codominant model (OR = 0.64, 95% CI = 0.41–1.00, p = 0.049).
We next conducted linkage disequilibrium (LD) analysis for the polymorphisms in MC4R1 and PDX1. Our results indicated two blocks (block1: rs11619319 and rs2293941; block2: rs9581943 and rs7981781) in PDX1 (Fig. 1) and two blocks (block1: rs6567160, rs663129, and rs17782313; block2: rs11663816 and rs12970134) in MC4R (Fig. 2). Besides, Table 6 shows that there was no association between haplotype frequency and T2DM risk (p > 0.05).
The relative mRNA expression of PDX1 and MC4R
The MC4R mRNA expression levels in T2DM case subjects decreased compared with those in their nondiabetic counterparts (p = 0.040, Fig. 3a). In addition, although no significant differences were observed in the expression levels of PDX1 mRNA between the two groups, we did observe a decreased pattern of PDX1 expression in individual samples between the cases and controls (p = 0.054, Fig. 3b).
The association of relative mRNA expression and PDX1 and MC4R polymorphisms
The PDX1 and MC4R polymorphisms were not associated with the relative PDX1 and MC4R mRNA expression in the T2DM patients and controls (Figs. 4, 5).
This research focused on the association of PDX1 and MC4R polymorphisms with susceptibility to T2DM in Chinese Han people. We found that PDX1-rs9581943 was correlated with a decreased risk of T2DM among the study subjects. In addition, the effects of PDX1 and MC4R polymorphisms on T2DM susceptibility were dependent on age, sex, smoking status, drinking status and BMI. These findings suggest that genetic polymorphisms in PDX1 and MC4R may play a crucial role in the development of T2DM.
In humans, the PDX1 gene is located on chromosome 13q12.1. It is a key transcription factor involved in pancreatic development, islet hormone and insulin expression. Data from several studies suggested that deletion and mutation in PDX1 caused overt diabetes and maturity-onset diabetes of the young [21, 22]. Additionally, Steinthorsdottir et al. found that rare frameshift variants in PDX1 were associated with a higher risk of T2DM in Icelanders . Recently, a homozygous mutation in PDX1 was detected in a 65-day-old Iranian patient with neonatal diabetes . However, there are few studies on rs11619319, rs2293941, rs9581943, and rs7981781. In the present study, we found that only rs9581943 decreased the incidence of T2DM among the study subjects. Moreover, we found that the relative mRNA expression of the PDX1 gene was lower in T2DM patients than in controls, but the difference was insignificant. Interestingly, stratified analysis results revealed that rs9581943, rs11619319, rs2293941, and rs7981781were associated with susceptibility to T2DM in different subgroups. Manning et al.  illustrated that rs2293941 was associated with fasting glucose levels in individuals of European ancestry. However, this correlation was not observed among participants in the Chinese Han population in the present study (not shown). The inconsistencies in these reports may result from subjects of different ethnicities and different environments. Taken together, these results demonstrated that the PDX1 polymorphism is important in the development and risk assessment of T2DM.
MC4R is a G-protein-coupled receptor that is highly expressed in the hypothalamus, where it regulates appetite, energy expenditure and body weight . It is located on chromosome 18q21 in humans. Disruption of the MC4R gene leads to the obesity phenotype, which is related to T2DM . Vaisse et al. claimed that rare heterozygous MC4R variants have been identified in obese children and adults in many populations . Obesity is an important risk factor for the progression of T2DM .
Herein, we explored whether MC4R polymorphisms could contribute to T2DM risk in a Chinese Han population. In this study, we found that the mRNA level of MC4R was decreased in T2DM patients compared to healthy controls. However, the overall analysis revealed that the association between MC4R polymorphisms and T2DM risk was insignificant. Subsequently, we examined the correlation of MC4R polymorphisms and T2DM risk by stratification analysis. We found that rs17782313 in MC4R obviously reduced the susceptibility toT2DM among individuals younger than 60 years old. It has previously been demonstrated that the MC4R-rs17782313 polymorphism is strongly related to obesity in adults and children of European descent . Moreover, Hardy et al. also demonstrated that rs17782313 was associated with weight and BMI. The association of this polymorphism with weight strengthened during childhood and adolescence, and weakened during adulthood . This result suggested that the effect of MC4R-rs17782313 on disease risk was dependent on age. In addition, a study showed that rs12970134 increased the risk of T2DM among individuals of European descent , although this effect was not found in our study. In our analysis, rs6567160 reduced the susceptibility to T2DM among individuals ≤ 60 years old but was not associated with the clinical characteristics. However, Carvalho et al. suggested that rs6567160 was associated with a greater postpartum increase in HbA1c in women who had experienced gestational diabetes mellitus than in those who had not . Additionally, rs663129 decreased the risk of T2DM among Han Chinese people. This finding was inconsistent with the discovery of Nikpay et al., which indicated that allele A of rs663129 increased the risk of both coronary artery disease and obesity in individuals of European ancestry . The reason for these inconsistent results may be that the occurrence and development of T2DM are related to a variety of factors, including population, sample size, and environment. Together, these data highlighted the important role of MC4R polymorphisms in the occurrence of T2DM.
Moreover, these selected SNPs in the PDX1 and MC4R genes can affect promoter histone marks, enhancer histone marks, DNAse, proteins bound, motifs changed, NHGRI/EBI GWAS hits, and GRASP QTL hits. Therefore, we presumed that these functions could modify the risk of T2DM by influencing gene expression. The specific mechanisms underlying these effects require further investigation.
There were several limitations in this study. First, this research was performed based on a Chinese Han population. Therefore, further research with subjects of different genetic backgrounds should be conducted to validate our results. Second, selection bias was an unavoidable problem in our research.
In conclusion, our findings demonstrated that the variants in the PDX1 and MC4R genes were related to susceptibility to T2DM in the Chinese Han population. These single polymorphic markers are considered to be new targets in the assessment and prevention of T2DM among Chinese Han people.
Availability of data and material
The datasets generated and/or analyzed during the current study are not publicly available due patient privacy but are available from the corresponding author on reasonable request.
type 2 diabetes mellitus
pancreatic and duodenal homeobox-1
melanocortin receptor 4
95% confidence intervals
Chatterjee S, Davies MJ. Accurate diagnosis of diabetes mellitus and new paradigms of classification. Nat Rev Endocrinol. 2018;14:386–7.
Saeedi P, Petersohn I, Salpea P, Malanda B, Karuranga S, Unwin N, et al. Global and regional diabetes prevalence estimates for 2019 and projections for 2030 and 2045: results from the International Diabetes Federation Diabetes Atlas. Diabetes Res Clin Pract. 2019;157:107843.
Cho NH, Shaw JE, Karuranga S, Huang Y, da Rocha Fernandes JD, Ohlrogge AW, et al. IDF diabetes Atlas: global estimates of diabetes prevalence for 2017 and projections for 2045. Diabetes Res Clin Pract. 2018;138:271–81.
Steinthorsdottir V, Thorleifsson G, Sulem P, Helgason H, Grarup N, Sigurdsson A, et al. Identification of low-frequency and rare sequence variants associated with elevated or reduced risk of type 2 diabetes. Nat Genet. 2014;46:294–8.
Gurevich LE, Proshchina AE, Voronkova IA, Ashevskaya VE, Korosteleva PA. Dolzhansky OV [Differential diagnostic value of the expression of the transcription factor PDX-1 in neuroendocrine and non-neuroendocrine tumors of the pancreas and other organs]. Arkh Patol. 2019;81:11–21.
American Diabetes Association. Standards of medical care in diabetes—2010. Diabetes Care 2010;33(Suppl 1):S11–61.
Du J, Jin T, Cao Y, Chen J, Guo Y, Sun M, et al. Association between genetic polymorphisms of MMP8 and the risk of steroid-induced osteonecrosis of the femoral head in the population of northern China. Medicine (Baltimore). 2016;95:e4794.
Weng J, Macfarlane WM, Lehto M, Gu HF, Shepherd LM, Ivarsson SA, et al. Functional consequences of mutations in the MODY4 gene (IPF1) and coexistence with MODY3 mutations. Diabetologia. 2001;44:249–58.
Manning AK, Hivert MF, Scott RA, Grimsby JL, Bouatia-Naji N, Chen H, et al. A genome-wide approach accounting for body mass index identifies genetic variants influencing fasting glycemic traits and insulin resistance. Nat Genet. 2012;44:659–69.
Mountjoy KG, Mortrud MT, Low MJ, Simerly RB, Cone RD. Localization of the melanocortin-4 receptor (MC4-R) in neuroendocrine and autonomic control circuits in the brain. Mol Endocrinol. 1994;8:1298–308.
Morris AP, Voight BF, Teslovich TM, Ferreira T, Segrè AV, Steinthorsdottir V, et al. Large-scale association analysis provides insights into the genetic architecture and pathophysiology of type 2 diabetes. Nat Genet. 2012;44:981–90.
de Carvalho AM, Shao P, Liu H, Cheng H-L, Zheng Y, Leng J, et al. The MC4R genotype is associated with postpartum weight reduction and glycemic changes among women with prior gestational diabetes: longitudinal analysis. Sci Rep. 2017;7:9654.
Nikpay M, Turner AW, McPherson R. Partitioning the pleiotropy between coronary artery disease and body mass index reveals the importance of low frequency variants and central nervous system-specific functional elements. Circ Genom Precis Med. 2018;11:e002050.
We sincerely thank all participators in this study.
This study was supported by the Bethune-Merck Diabetes Research Foundation (G-X-2019-056), the Clinical Research Project of the First Affiliated Hospital of Xi’an Jiaotong University (XJTU1AF-CRF-2019-007) and the Social Development Science and Technology Project of Shaanxi Province (2016sf-314).
Co-first author: Ning Wang and Rui Tong
Authors and Affiliations
Department of Endocrinology and Second Department of Geriatrics, The First Affiliated Hospital of Xi’an Jiaotong University, 277 West Yanta Road, Xi’an, 710061, Shaanxi, China
Ning Wang, Rui Tong, Jing Xu, Jiaqi Cui, Huan Chen, Yanqi Peng, Sijia Fei, Shujun Yang, Lu Wang, Juanchuan Yao & Wei Cui
Department of Oncology, East Branch of the First Affiliated Hospital of Xi’an Jiaotong University, Xi’an, 710089, Shaanxi, China
Department of Endocrinology, Xianyang Central Hospital, Xianyang, 712000, Shaanxi, China
NW, RT and WC designed this study protocol and drafted the manuscript; JX and JC performed the DNA extraction and genotyping; YT and JP performed the data analysis; HC, YP, SF and SY performed the sample collection and information recording. LW and JY polished the manuscript, WC conceived and supervised the study. All authors read and approved the final manuscript.
This research received approval from the Ethics committee of the Frist Affiliated Hospital of Xi’an Jiaotong University, and conformed to the Declaration of Helsinki. Informed consent was acquired from each participant at recruitment after fully describing our research to them.
Consent for publication
The authors declare that they have no conflict of interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Primer sequences of PDX1 and MC4R for PCR. Table S2. Basic information of candidate SNPs in the study.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
Wang, N., Tong, R., Xu, J. et al.PDX1 and MC4R genetic polymorphisms are associated with type 2 diabetes mellitus risk in the Chinese Han population.
BMC Med Genomics14, 249 (2021). https://doi.org/10.1186/s12920-021-01037-3