Skip to main content
Download PDF

Association analysis of SOCS3, JAK2 and STAT3 gene polymorphisms and genetic susceptibility to type 2 diabetes mellitus in Chinese population

Diabetology & Metabolic Syndrome volume 14, Article number: 4 (2022) Cite this article

Abstract

Aim

The association of polymorphisms in the three genes of SOCS3, JAK2 and STAT3 with genetic susceptibility to type 2 diabetes mellitus (T2DM) was explored, and its interaction with environmental factors such as hypertension and triglycerides was analyzed.

Methods

The Hardy–Weinberg balance test was used to analyze the random balance of genes in the population. The analysis of the association of SNPs with T2DM was performed using Pearson’s chi-square test. Haplotype frequency distribution, SNPs-SNPs interaction and environmental factors were analyzed by chi-square test and logistic regression.

Results

The genotype distribution of SNPs rs2280148 of the SOCS3 gene was statistically significant. The allele frequency distribution of SNPs (rs4969168/rs2280148) was statistically different. After covariate correction, the SOCS3 gene locus (rs4969168) showed an association with T2DM in additive model, while the rs2280148 locus showed an association with T2DM in all three models. The locus (rs10974914/rs10815157) allele and genotype frequency distribution of JAK2 were statistically significant. After covariate correction, two SNPs in the gene showed association with T2DM in both additive and recessive models. The distribution of genotype frequencies of SNPs rs1053005 locus in gene STAT3 was statistically significant between the two groups. In recessive genetic models, rs1053005 locus polymorphisms was associated with T2DM. Haplotype S3 (G G)/S 4 (G T) of the SOCS3 gene as well as haplotype J2 (A G)/J 3 (G C) of the JAK2 gene were closely associated with T2DM. There was an interaction between SNPs rs4969168 and SNPs rs2280148 in the SOCS3 gene. There was an interaction between the SOCS3, JAK2 and STAT3 genes and hypertension/triglycerides.

Conclusion

The SOCS3 and JAK2 genes may be associated with T2DM in the Chinese population, in which SNPs carrying the A allele (rs4969168)/G allele (rs2280148)/C allele (rs10815157) have a reduced risk of T2DM. Haplotype S3 (G G)/S 4 (G T) of the SOCS3 gene and haplotype J2 (A G)/J 3 (G C) of the JAK2 gene may be influencing factor for T2DM. The interaction between SNPs rs4969168 and SNPs rs2280148 increases the risk of T2DM. Hypertension and triglycerides may interact with SNPs of T2DM susceptibility genes.

Introduction

Diabetes has quadrupled in the past 30 years and is the ninth leading cause of death [1]. Currently, approximately 1 in 11 adults (20–79 years) worldwide have diabetes (463 million people), and 90% of them have type 2 diabetes mellitus (T2DM) [2]. Diabetes is expected to become the seventh leading cause of death in developing countries by 2030 due to a rapid rise in the incidence of diabetes [3]. The epidemic of T2DM is mainly found in Asia, and is most pronounced in countries with middle economic level, such as China [2, 3]. It is one of the major non-communicable health problems in the world [4]. T2DM is mainly manifested as the deficiency of the body’s own islet function or insulin dysfunction, which is a multi-factor disease with the joint action of genetic factors and environmental factors [5,6,7]. Genetic susceptibility is one of the main influencing factors of T2DM, which is characterized by a higher concordance rate in identical twins and racial differences in the familial aggregation [8,9,10]. Genetic susceptibility to T2DM has led to an intensive study for nearly past few decades. However, the specific genetic architecture of T2DM risk and the underlying mechanism by which the risk gene of T2DM remains poorly understood. Symptoms of T2DM are often accompanied by inflammatory responses [11]. In the inflammatory response, the JAK/STAT signaling pathway plays an important role and affects the insulin signaling pathway under the SOCS3 (suppressor of cytokine signaling3) gene. Therefore, we hypothesized that the SOCS3, JAK2(Janus kinase 2), and STAT3 (Signal transducer and activator of transcription 3) genes involved in the JAK/STAT signaling pathway may have an association for genetic susceptibility to T2DM.

SOCS3 is a cellular signaling inhibitory protein with a relative molecular mass of 24.7 kd, consisting of 225 amino acids. It plays an important role in cell signaling, with functions such as regulating cell physiological processes, regulating inflammation, inhibiting proliferation of cancer cells and metastasis [12]. Recent studies have found that the overexpression of SOCS3 is closely associated with insulin resistance and blood glucose levels in diabetic patients. SOCS3 as a inhibitory factor of insulin signaling interacts with the insulin receptor or induces the degradation of signaling proteins that in turn suppresses the conduction of insulin signaling [13, 14]. At the same time, it can lead to insulin resistance by participating in the conduction processes of other cytokines (TNF-α, IL6, and PPAR) signaling pathways [14]. Insulin can affect β cell function, enhance insulin gene transcription, and promote or inhibit β cell apoptosis [15]. In β cells, the SOCS3 gene is able to regulate insulin levels by regulating leptin and inhibiting the conduction of JAK/STAT signal channel [16]. Thus, we speculate that the SOCS3 gene may be a risk factor for T2DM.

The JAK2 and STAT3 genes are key genes in the JAK/STAT signaling pathway. The molecule of JAK2 contains two kinase domains: the C-terminal domain contains the structures associated with the tyrosine protein kinase and has the catalytic activity, while the N-end domain lacks the amino acids associated with the catalytic phosphorylation reactions and thus has no tyrosine kinase activity [17]. The JAK2 gene polymorphisms was found to be closely associated with the metabolic syndrome in T2DM and related complications [18]. Chen et al. [19] found that the overexpressed JAK2 gene can increase the activity of the insulin promoter, with opposite effects when the JAK2 gene is silenced or poorly expressed. Amyloid-β protein can lead to insulin resistance as well as the development of T2DM through the activation of the JAK2/STAT3/SOCS1 pathway [20]. We also found that Lep-STAT3 signaling pathway is closely associated with metabolic diseases such as hyperglycemia, hyperinsulinemia and glucose tolerance [21]. It has been found that the JAK2 and STAT3 genes play an important role in the pathogenesis of T2DM by participating in the JAK/STAT signaling pathway.

T2DM is a complex disease controlled by polygenes, and environmental factors also have an important role in the pathogenesis [22]. Simply analyzing the association of only a single gene or individual environmental factors with T2DM often leads to the neglect of some factors with weak effects. Therefore, we analyzed the association of the three genes (SOCS3/JAK2/STAT3) and their SNPs-SNPs interaction with T2DM in the Chinese Guangdong population and analyzed the role of environmental factors involved in it. The aim was to explore the relationship between genetic polymorphisms in SOCS3, JAK2 and STAT3 and the genetic susceptibility to T2DM.

Materials and methods

Study population

This study was a case control study. In 2013, 599 T2DM patients (30–70 years) were recruited from endocrinology departments of eight hospitals in Guangdong Province, China. A total of 677 healthy subjects who were enrolled in the same physical examination centers of the 8 hospitals were selected as the normal control (NC) group. Both groups were Han people in Guangdong province and had no genetic relationship. The consistency of genetic background is high in the population. The project has been approved by the Ethics committee. All subjects were been informed of the specific study background and signed informed consent. We used the T2DM diagnostic criteria of Diabetes proposed by American Diabetes Association (ADA): (1) HbA1c (glycated hemoglobin A1c) ≥ 6.5%; (2) FPG (fasting blood-glucose) ≥ 7.0 mmol/L; and (3) Blood glucose was 2 h after glucose load ≥ 11.1 mmol. Control group inclusion criteria: (1) Age 30–70; (2) No family history of diabetes; (3) No other serious diseases, such as cardiovascular and cerebrovascular diseases; and (4) No abnormalities in biochemical results such as physical examination and blood glucose test.

Information inquiry

Subjects who met the inclusion criteria were given a questionnaire and physical examination. A survey scale suitable for the purpose of the study was designed. The survey mainly includes four aspects: demographic information, lifestyle information, clinical examination data and history of drug taking. Standardized questionnaire survey was conducted by uniformly trained investigators, including general information (age, height, weight, sex, native place and occupation, etc.) and medical history (disease history, disease course, smoking history, family history, complications, diet and exercise status, etc.). Calculate Body Mass Index [BMI = weight (kg)/height (m2)]. The blood pressure of the right arm was measured with a desktop mercury sphygmomanometer twice at an interval of 10 min, and the average value was taken. The diagnostic criteria for hypertension were referred to the 1999WHO/ISH guidelines for hypertension: (1) Systolic blood pressure (SBP ≥ 140 mmHg); and (2) Diastolic blood pressure (DBP ≥ 90 mmHg).

Blood sample collection

Professional endocrinology nurses collected 9 ml of peripheral blood in the morning from all subjects with at least 8 h of fasting. 5 ml of blood samples were used to detect clinical biochemical indicators. Including triglycerides (TG), fasting blood glucose (FPG), high density lipoprotein cholesterol (HDL-C), total cholesterol (TC), low density lipoprotein cholesterol (LDL-C) and glycated hemoglobin A1c (HbA1c). Glycosylations method was used for determination of FPG level and enzyme method was used for determination of TC, TG, HDL-C and LDL-C. HbA1c was determined by high performance liquid chromatography (BIO-Radd-10 TM Hb A1c detection system). The other biochemical parameters were determined by automatic biochemical analyzer (Bio-RAD Corporation, USA).

Extracted DNA

The remaining 4 ml blood sample was used to extract DNA. Ethylene Diamine Tetraacetic acid-K2 (EDTA-K2) was used for adequate anticoagulation treatment of blood samples. After digestion by Protease K, DNA was extracted by salting out method. Store the extract in a refrigerator at − 80℃. For specific extraction steps, we referred to the corresponding technical instructions:

  1. (1)

    The blood samples were transferred from the anticoagulant tube to the 10 ml centrifuge tube and marked on the centrifuge tube.

  2. (2)

    2 ml cell lysis solution was added to the blood samples containing anticoagulant, and centrifuged for 10 min (5000 r/min, r = 7.1 cm) after mixing.

  3. (3)

    Remove the liquid from the top of the centrifuge tube, then invert it on a clean blotting paper for 1 min to dry.

  4. (4)

    A mixture of buffer and protease K is configured in time (1 mL buffer + 10 μL protease K solution).

  5. (5)

    Add 1 mL mixture to each sample and mix thoroughly.

  6. (6)

    Take a water bath at 65℃ for 20 min, add 1 mL isopropyl alcohol, and mix thoroughly until filamentous or clumps appear.

  7. (7)

    Centrifugation was performed for 15 min (5000 r/min, r = 7.1 cm), and the upper liquid was discarded.

  8. (8)

    Add 1 mL 70% ethyl alcohol, transfer it into a 1.5 mL centrifuge tube, centrifuge it for 10 min (8000 r/min, r = 4.2 cm), and the upper liquid was discarded. Perform two operations

  9. (9)

    The tube is dried at room temperature until the liquid in the DNA precipitate evaporates.

  10. (10)

    The corresponding elution buffer was added according to the mass of DNA, and then the mixture was bathed at 65℃ for 1 h. During this time, flick the finger several times to help blend and centrifuge for 30 s.

  11. (11)

    DNA quality was tested by nucleic acid content detector, and DNA concentration (A260/280 and A260/230) was read to judge the purity of DNA. Finally, DNA was stored in a refrigerator at − 80 ℃.

SNPs selection and genotyping

We selected tagSNPs [MAF (minor allele frequency) > 0.05, R2 ≥ 0.8] from Upstream and downstream 5 kb region of target gene by Haplotype Mapping (HapMap) public database and Hapoloview (ver.4.2). TagSNPs with high prediction scores were screened by FastSNP (function analysis and selection tool for single nucleotide polymorphisms, http://fastsnp.ibms.sinica.edu.tw). Then six Tag SNPs (rs4969168/rs2280148, rs10974914/rs10815157 and rs1053005, rs957970) were selected from the three genes (SOCS3/JAK2/STAT3). The TaqMan® SNPs probe method was used for typing. The basic principle of this technique is to use the high specificity of ligase chain reaction to realize the recognition of SNPs alleles. Then non-specific sequences of different lengths were introduced into the end of the ligator probe and ligase chain reaction was used to obtain the ligator products of different lengths corresponding to the loci. Universal primers labeled with fluorescence were used to PCR amplify the ligands, and electrophoresis was used to separate the amplified products. Finally, the genotypes of the loci were obtained by analysis of electrophoresis patterns. The corresponding sequences are shown in Table 1.

Table 1 Six sequences corresponding to the SNPs loci
Full size table

The specific steps are as follows:

  1. (1)

    Items prepared: crushed ice, DNA samples and typing reagents. (Genotype Assay* and Universal PCR Master Mix*).

  2. (2)

    The DNA samples were defrosted at low temperatures. Combine the mixture in proportion. (20 × Genotyping assay = 0.25*N*1.1 μl, 2 × universal PCR Master Mix = 2.5*N*1.1 μl, ddH2O = 4.5*N *1.1 μl).

  3. (3)

    Add mixture (4.5 μl) and DNA (1-20 ng) to Orifice.

  4. (4)

    Design allele typing template (AD file) and PCR amplification template (AQ file). The PCR reaction condition was 95℃ for 10 min. 40 PCR cycles were performed: denaturation at 92 °C for 15 s, annealing at 60 °C for 1 min.

  5. (5)

    Raw data were collected and analyzed on a 7900HT fluorescent quantitative PCR instrument. The fluorescence labeling curves of PCR products and the corresponding SNPs locus/allele information were obtained.

*Reagents purchased from Applied Biosystems.

In the typing test, 10% of the samples were randomly selected from the case group and the control group respectively, and the experiment was repeated. The results of the two typing experiments were compared, and the false positive rate was less than 5%, indicating that the results were reliable.

Statistical analysis

Quantitative measurement data were expressed as mean ± standard deviation, and qualitative data were expressed as case number or percentage. Student’s t-test was used to analyze the difference between the two groups of continuous quantitative data. Genotypes and allele frequencies were compared by Chi-square test. The HWE (Hardy–Weinberg equilibrium) test of genotype distribution was performed by χ2 test of goodness of fit. LD block haplotype analysis was performed by unconditional Logistic regression analysis. Allele association analysis, genotype association analysis and Cochran-Armitage trend test were used to analyze the relationship between SNPs and T2DM. Haplotype frequency distribution, SNPs-SNPs interaction, and environmental factors were analyzed by chi-square test and logistic regression. P-value < 0.05 was considered as statistical significance.

Results

Baseline data

In order to prevent a small number of subjects with poor genotyping results from bias to the overall analysis, individuals with a loss rate of over 20% SNPs genotyping were excluded. Finally, there were 599 cases and 677 normal controls (NC) selected for further study. Among them, there were 304 male (50.75%) in T2DM group, with the average age of 58.23 years, while there were 424 male (62.23%) in NC group, with the average age of 58.12 years. The BMI, FPG and TG of NC group were lower than those of T2DM group (P < 0. 001). The baseline data is shown in Table 2.

Table 2 Comparison of general data and biochemical indexes between the two groups
Full size table

Analyzing of the allele and genotype frequency

Six candidate SNPs loci in the three genes were tested by Hardy-Winberg equilibrium test (α = 0.01). The rs957970 locus was removed (P < 0.001). The remaining loci indicate that the frequency distribution of genotype in the object population conforms to the law of genetic balance and is representative of the population. See Table 3.

Table 3 Results of Hardy–Winberg equilibrium test
Full size table

The allele/genotype distribution of the T2DM patients and normal persons is showed in Tables 4, 5. Table 4 showed that the allelic distribution of rs4969168, rs2280148, rs10974914 and rs10815157 were statistically significant in the two groups (P < 0.05). After adjusting for covariates (age, sex, BMI), the difference between the two groups was still statistically significant (P < 0.05). Chi-square test was performed for two loci in SOCS3, JAK2 and STAT3 genes. Table 5 shows that three candidate genes may be associated with T2DM.

Table 4 The results of association analysis between alleles and T2DM
Full size table
Table 5 The results of association analysis between genotype and T2DM
Full size table

Genetic model analysis

To further confirm the relationship between polymorphisms and T2DM at each gene locus, we performed a univariate logistic regression analysis under different genetic models (additive model, dominant model, and recessive model). Suppose A and a are two alleles at the locus, and A is the risk allele. The corresponding genotypes were AA, Aa and aa, and the sample size was N. Allele frequency calculation of each gene: the frequency of A gene is P(A) = (NAa + 2NAA)/2 N. The association between SNP alleles or genotypes and T2DM was analyzed based on additive, dominant and recessive gene models. Genotype corresponding encoding (additive model: AA = 2, Aa = 1, aa = 0; Dominant model: AA = 1, AA = 1, AA = 0; The implicit model is AA = 1, AA = 0, AA = 0).

To make the analysis results more reliable, four covariates (age, sex, BMI) were adjusted under three models. Where locus rs2280148 was statistically different in three models, while locus rs4969168 associated with T2DM in additive and dominant models. Two SNPs loci (rs10974914 and rs10815157) in the JAK2 gene were significantly associated with T2DM in the additive and recessive models (P < 0.05). In the recessive model, the rs105300 locus was associated with T2DM [P = 0.002, OR (95% CI) = 1.91(1.28–2.86)], as detailed in Table 6.

Table 6 Comparison of genetic models between T2DM and NC
Full size table

Trend tests are based on predefined scores, which correspond to different genetic models. Only when the genetic model developed is the data optimal model, the corresponding trend testing effect is the greatest. The test results show the high reliability of the model analysis, and the results are detailed in Table 7.

Table 7 Analysis results of three robustness test methods
Full size table

Haplotype analysis of JAK2 and STAT3 genes

The SNPs rs4969168 and SNPs rs2280148 in SOCS3 gene were constructed haplotypes, and the four haplotypes were named S1-S4. The SNPs rs10974914 and SNPs rs10815157 in JAK2 gene were constructed haplotypes, and the four haplotypes were named J1-J4. Chi-square test showed that haplotypes S3 (GG) and S4 (GT) had statistically significant differences in haplotype frequency distribution between the case group and the control group (P < 0.001). Same thing with J2 (AG) and J3 (GC). According to their frequency distribution, haplotypes of S3 (GG) and J3 (GC) may be the protective factors for T2DM, while haplotypes of J2 (AG) and S4 (GT) may be the risk factors for T2DM. See Table 8 for details.

Table 8 Results of SOCS3 and JAK2 gene haplotype analysis
Full size table

Analysis of SNPs-SNPs interaction among SOCS3, JAK2 and STAT3 genes

The interaction analysis between SNPs-SNPs was performed using multivariate nonconditional Logistic regression models as well as the PLINK software. The results of the analysis revealed an interaction between rs4969168 and rs2280148 in the SOCS3 gene (P < 0.001). This is a statistically significant SNPs-SNPs interaction pair (α = 0.05). After adjusting for the four covariates (age, sex, BMI), an interaction remained between rs4969168 and rs2280148 (P < 0.001). See Tables 9, 10 for details.

Table 9 The interaction of the loci between the genes
Full size table
Table 10 Interaction of each gene locus after covariate correction
Full size table

Interaction between the SNPs site and hypertension/triglycerides

The results showed that the rs10815157 locus interacts with hypertension in the recessive model [P = 0.03, OR = 0.37, 95% CI = (0.16–0.89)] as a protective factor for T2DM. The rs1053005 locus interacts with hypertension as a risk factor for T2DM [P = 0.01, OR = 3.23, 95% CI = (1.43–7.29)]. The SNPs rs4969168 locus interacts with triglycerides in the additive model, while the SNPs rs10974914 locus interacts with triglycerides in both the additive and recessive models. In the additive model, both SNPs rs4969168 and SNP rs10974914 interact with triglycerides. See Tables 11, 12 for details.

Table 11 Interaction analysis of SNPs with hypertension
Full size table
Table 12 Interaction analysis of SNPs with triglycerides
Full size table

Discussion

T2DM is a complex disease caused by the interaction of environmental and genetic factors [23]. In complex diseases, the effect of a single gene is weak, and generally a variety of factors are involved in the occurrence and development of the disease [24]. Haplotype analysis of multiple genes and interactions between genes provide the possibility of finding such a weak effect [25]. Due to genetic heterogeneity among different populations, we selected 1276 people from Guangdong Han population to study the association between single nucleotide polymorphisms (SNPs) of gene loci and genetic susceptibility to T2DM.

This study found that the allele frequency and genotype frequency of rs2280148, the SNP locus in SOCS3 gene, were statistically significant different between the T2DM group and the control group. The "GG" genotype at rs2280148 reduced the risk of T2DM, indicating that the "GG" genotype was resistant to T2DM. It was closely related to T2DM in all three models (P < 0.05). Obesity is a risk factor for T2DM. Matthew E Talbert also found that SNPs rs2280148 in SOCS3 gene is closely related to various types of obesity in the American population [26]. Studies have shown that single nucleotide polymorphisms in SOCS3 gene are closely related to obesity and BMI [27, 28]. The above evidences prove that rs2280148 is correlated with T2DM. At the same time, SNP rs2280148 locus is located in the UTR3 of SOCS3, which may be involved in protein transcription and affect the incidence of T2DM. The specific biological mechanism needs further downstream functional verification. The alleles and genotype frequencies of rs10974914 and rs10815157 in JAK2 gene showed statistically significant differences between the diabetic group and the control group. The "AA" genotype at rs10974914 significantly increased the risk of T2DM, while the "C" allele at rs10815157 was the resistance factor for T2DM. Logistic regression analysis showed that two SNPs in JAK2 gene were associated with T2DM in the additive model and recessive model. Zhang [29] found that the overexpression of JAK2 gene could increase the activity of insulin promoter and affect the incidence of T2DM indirectly. There are studies that have reported that Aβ induces hepatic insulin resistance by activating JAK2/STAT3/SOCS-1 signaling pathway and have implications on resolving insulin resistance and T2DM [30]. Therefore, we believe that JAK2 gene may be associated with T2DM in Guangdong Han population and participates in a variety of physiological responses as a member of the JAK/STAT inflammatory signaling pathway.

The analysis results show that haplotypes of S3 (GG) and J3 (GC) may be the protective factors for T2DM, while haplotypes of J2 (AG) and S4 (GT) may be the risk factors for T2DM. However, in STAT3 gene, only one locus in the analysis process met Hardy-Winberg equilibrium test, which did not meet the conditions of haplotype analysis. If forced analysis was carried out, the results would be unstable. After covariate correction, the SNPs-SNPs(rs4969168-rs2280148) interaction of SOCS3 was statistically significant (OR = 2.450, P < 0.001). These results indicate that there may be a synergistic effect between the two loci of SOCS3 gene, leading to both SNPs loci associated with T2DM and risk factors for T2DM. However, further functional and epigenetic verification is needed to determine the type of interaction between the two SNPs loci. In recessive model, the SNPs rs10815157 interacted with hypertension and was a protective factor of T2DM, indicating that the prevalence of diabetes in hypertensive population with homozygous "C" risk allele was lower than that in hypertensive population without this risk gene. In the recessive model, the SNPs rs1053005 also interacts with hypertension and is a risk factor for T2DM, indicating that hypertensive patients carrying "CC" homozygous gene are more likely to develop diabetes than those without this allele. Martinez et al. [31] showed that SOCS3 gene and hypertension factors jointly act on the population, leading to the onset of T2DM. The SNPs rs4969168 was shown to interact with triglycerides in the additive model, and SNPs rs1053005 was shown to interact with triglycerides in the dominant model. The interaction between the two loci and triglycerides both a risk factor for T2DM. At the same time, the SNPs rs10974914 had an interaction with triglycerides in both additive and recessive models, and this interaction was a factor of resistance to T2DM. Torisu et al. [32] showed that SOCS3 gene and triglyceride could act together to cause the onset of T2DM. The influencing factors of T2DM do not act on individuals alone, but are a combination of multiple factors, and under the action of time effect, genetic factors and environmental factors jointly lead to the occurrence of the disease. The exact mechanism of action is not clear, it may have weak effects on different protein expression and epigenetics aspect, and act on individuals through its combined effect, which is worthy of further study.

Conclusions

In conclusion, we conducted a preliminary study on the association between SOCS3, JAK2, and STAT3 gene polymorphisms and T2DM in China. The related loci of genes (SOCS3, JAK2) may play a key role in the occurrence and development of T2DM. Gene polymorphisms may provide clues for further study of the pathogenesis of T2DM and search for new drug targets. The different specificity of genetic heterogeneity in T2DM among different populations will be an interesting research direction in the future. Our study used gene haplotype analysis, SNP-SNP interaction analysis, and SNP-environment interaction analysis. We were able to detect moderate and weak susceptibility to T2DM. However, the experiment also has some defects and inadequacies. The sample size and the number of SNPs loci of candidate genes are relatively small, and downstream functional verification has not been verified. As a result, the evidence of biological function is insufficient. At the same time, the association between the six SNPs loci selected in this study and T2DM has rarely been reported. Therefore, a larger sample size, multi-regional, multi-ethnic population study needs to be further carried out.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

References

  1. Zheng Y, Ley SH, Hu FB. Global aetiology and epidemiology of type 2 diabetes mellitus and its complications. Nat Rev Endocrinol. 2018;14(2):88–98.

    Article  Google Scholar 

  2. International Diabetes Federation. IDF diabetes atlas 9th edition 2019. Diabetes Atlas[EB/OL]. http://www.diabetesatlas.org/. Accessed 3 Oct 2021.

  3. Lu Y, Xu J, Zhao W, et al. Measuring self-care in persons with type 2 diabetes: a systematic review. Eval Health Prof. 2016;39(2):131–84.

    Article  Google Scholar 

  4. Raaijmakers LG, Martens MK, Bagchus C, et al. Correlates of perceived self-care activities and diabetes control among Dutch type 1 and type 2 diabetics. J Behav Med. 2015;38(3):450–9.

    Article  Google Scholar 

  5. Kim K, Kim ES, Rhee SY, et al. Clinical characteristics and risk factors for retinal diabetic neurodegeneration in type 2 diabetes. Acta Diabetol. 2017;54:993–9.

    Article  CAS  Google Scholar 

  6. Zhou Z, Mahdi A, Tratsiakovich Y, et al. Erythrocytes from patients with type 2 diabetes induce endothelial dysfunction via arginase I. J Am Coll Cardiol. 2018;72:769–80.

    Article  CAS  Google Scholar 

  7. Zhang QJ, Yang CC, Zhang SY, et al. Alteration of NPY in hypothalamus and its correlation with leptin and ghrelin during the development of T2DM in a rat model. SpringerPlus. 2016;5(1):1913.

    Article  Google Scholar 

  8. Zhang X, Gao Y, Zhou Z, et al. Familial clustering of diabetic retinopathy in Chongqing, China, type 2 diabetic patients. Eur J Ophthalmol. 2010;20(5):911–8.

    Article  Google Scholar 

  9. Weires MB, Tausch B, Haug PJ, et al. Familiality of diabetes mellitus. Exp Clin Endocrinol Diabetes. 2007;115(10):634–40.

    Article  CAS  Google Scholar 

  10. Lee JW, Brancati FL, Yeh HC. Trends in the prevalence of type 2 diabetes in Asians versus whites: results from the United States National Health Interview Survey, 1997–2008. Diabetes Care. 2011;34(2):353–7.

    Article  Google Scholar 

  11. Cucak H, Grunnet LG, Rosendahl A. Accumulation of M1-like macrophages in type 2 diabetic islets is followed by a systemic shift in macrophage polarization. J Leukoc Biol. 2014;95(1):149–60.

    Article  Google Scholar 

  12. Haan S, Ferguson P, Sommer U, et al. Tyrosine phosphorylation disrupts elongin interaction and accelerates SOCS3 degradation. J Biol Chem. 2003;278(34):31972–9.

    Article  CAS  Google Scholar 

  13. Emanuelli B, Peraldi P, Filloux C, et al. SOCS-3 is an insulin-induced negative regulator of insulin signaling. J Biol Chem. 2000;275(21):15985–91.

    Article  CAS  Google Scholar 

  14. Rui L, Yuan M, Frantz D, et al. SOCS-1 and SOCS-3 block insulin signaling by ubiquitin-mediated degradation of IRS1 and IRS2. J Biol Chem. 2002;277(44):42394–8.

    Article  CAS  Google Scholar 

  15. Rakatzi I, Ramrath S, Ledwig D, et al. A novel insulin analog with unique properties: LysB3, GluB29 insulin induces prominent activation of insulin receptor substrate 2, but marginal phosphorylation of insulin receptor substrate 1. Diabetes. 2003;52(9):2227–38.

    Article  CAS  Google Scholar 

  16. Seufert J. Leptin effects on pancreatic beta-cell gene expression and function. Diabetes. 2004;53(Suppl 1):S152-8.

    Article  CAS  Google Scholar 

  17. Bole-Feysot C, Goffin V, Edery M, et al. Prolactin (PRL) and its receptor: actions, signal transduction pathways and phenotypes observed in PRL receptor knockout mice. Endocr Rev. 1998;19(3):225–68.

    Article  CAS  Google Scholar 

  18. Kou ZZ, Li CY, Tang J, et al. Down-regulation of insulin signaling is involved in painful diabetic neuropathy in type 2 diabetes. Pain Phys. 2013;16(2):E71-83.

    Google Scholar 

  19. Zhang Y, Zhou B, Deng B, et al. Amyloid-beta induces hepatic insulin resistance in vivo via JAK2. Diabetes. 2013;62(4):1159–66.

    Article  CAS  Google Scholar 

  20. Ezzat S, Zheng L, Florez JC, et al. The cancer-associated FGFR4-G388R polymorphism enhances pancreatic insulin secretion and modifies the risk of diabetes. Cell Metab. 2013;17(6):929–40.

    Article  CAS  Google Scholar 

  21. W Jialing, Q Bo. Regulation of extinction on retinal JAK2-STAT, TM-1 signaling in DM rats. Trad Chi Med. 2011;10.

  22. Alberti KG, Zimmet PZ. Diabetes: a look to the future. Lancet Diabetes Endocrinol. 2014;2(1):e1-2.

    Article  Google Scholar 

  23. Liu T, Yang W, Pang S, et al. Functional genetic variants within the SIRT2 gene promoter in type 2 diabetes mellitus. Diabetes Res Clin Pract. 2018;137:200–7.

    Article  CAS  Google Scholar 

  24. Lander ES, Schork NJ. Genetic dissection of complex traits. Science. 1994;265:2037–48.

    Article  CAS  Google Scholar 

  25. Belsky J, Hartman S. Gene-environment interaction in evolutionary perspective: differential susceptibility to environmental influences. World Psychiatry. 2014;13:87–9.

    Article  Google Scholar 

  26. Talbert ME, Langefeld CD, Ziegler J, et al. Polymorphisms near SOCS3 are associated with obesity and glucose homeostasis traits in Hispanic Americans from the insulin resistance atherosclerosis family study. Hum Genet. 2009;125:153–62.

    Article  CAS  Google Scholar 

  27. Tang W, Zou J, Chen X, et al. Association of two polymorphisms within and near SOCS3 gene with obesity in three nationalities in Xinjiang province of China. Acta Pharmacol Sin. 2011;32:1381–6.

    Article  CAS  Google Scholar 

  28. Santillán-Benítez J, Mendieta-Zerón H, Gómez-Oliván L, et al. JAK2, STAT3 and SOCS3 gene expression in women with and without breast cancer. Gene. 2014;547:70–6.

    Article  Google Scholar 

  29. Zhang Y, Zhou B, Deng B, et al. Amyloid-β induces hepatic insulin resistance in vivo via JAK2. Diabetes. 2013;62:1159–66.

    Article  CAS  Google Scholar 

  30. Zhang Y, Zhou B, Zhang F, et al. Amyloid-β induces hepatic insulin resistance by activating JAK2/STAT3/SOCS-1 signaling pathway. Diabetes. 2012;61:1434–43.

    Article  CAS  Google Scholar 

  31. Martinez JA, Milagro FI, Claycombe KJ, et al. Epigenetics in adipose tissue, obesity, weight loss, and diabetes. Adv Nutr. 2014;5(1):71–81.

    Article  CAS  Google Scholar 

  32. Torisu T, Sato N, Yoshiga D, et al. The dual function of hepatic SOCS3 in insulin resistance in vivo. Genes Cells. 2007;12(2):143–54.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

We appreciate all authors for their contributions and physicians and participants. The study was funded by the Discipline Construction Project of Guangdong Medical University (4SG21276P), the Medical Scientific Research Foundation of Guangdong Province (Grant No. A2021395), the Natural Science Foundation of Basic and Applied Basic Research Foundation of Guangdong Province (Grant No. 2021A1515010061), the Undergraduate Innovation Experiment Project of Guangdong Medical University (Grant No. ZZDG003), the Natural Science Key Cultivation Project of Scientific Research Fund of Guangdong Medical University (Grant No. GDMUZ2020008), the Zhanjiang Science and Technology Development Special Fund Competitive Allocation Project (Grant No. 2020A01031), the Guangdong Provincial Colleges and Universities Characteristic Innovation Project (Grant No. 2019KTSCX047), the Basic and Applied Basic Research Foundation of Guangdong Province Regional Joint Fund Project (Grant No. 2020B1515120021) and the Young Innovative Talents Project of Guangdong Universities (Grant No. 2018KQNCX088).

Author information

Author notes

  1. Yang Zhang and Chunwen Lin have the same contribution

Authors and Affiliations

  1. Department of Epidemiology and Medical Statistics, School of Public Health, Guangdong Medical University, Dongguan, China

    Yang Zhang, Chunwen Lin, Rong Chen, Ling Luo, Jialu Huang, Hao Liu, Weiying Chen, Jian Xu, Haibing Yu & Yuanlin Ding

  2. Affiliated Hospital of Guangdong Medical University, Zhanjiang, China

    Yang Zhang

Authors
  1. Yang Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Chunwen Lin
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Rong Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Ling Luo
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Jialu Huang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Hao Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Weiying Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Jian Xu
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Haibing Yu
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Yuanlin Ding
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

All the authors made substantial contributions to the conception and design. HY and CL collected the data. RC and JH made substantial contributions to the analysis and interpretation of data. CL wrote the first draft of the manuscript. All the authors were involved in revising the manuscript critically for important intellectual content. YD gave final approval of the version to be published. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Haibing Yu or Yuanlin Ding.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Medical Ethics Committee of the Affiliated Hospital of Guangdong Medical (No. PJ2012079, China). All the surveys and samples were obtained with the consents of participants in advance, and the informed consent forms were legally consented.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, Y., Lin, C., Chen, R. et al. Association analysis of SOCS3, JAK2 and STAT3 gene polymorphisms and genetic susceptibility to type 2 diabetes mellitus in Chinese population. Diabetol Metab Syndr 14, 4 (2022). https://doi.org/10.1186/s13098-021-00774-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13098-021-00774-w

Keywords