Abstract

Objective. This study is aimed at analyzing the relationship between leptin (LEP) signaling pathway and type 2 diabetes mellitus (T2DM) and at providing support for molecular genetic research on the pathogenesis of T2DM in Chinese Han population. Methods. A case-control study was designed, including 1092 cases with T2DM and 1092 healthy controls of Chinese Han origin recruited from ten hospitals in Guangdong Province, Southern China. Twenty-three single nucleotide polymorphisms (SNPs) of 15 genes in LEP signaling pathway were genotyped by SNPscan™ kit. The Pearson chi-square test, Cochran-Armitage trend test, MAX3, and logistic regression were applied to analyze the association between single nucleotide polymorphism (SNP) and T2DM; unconditional logistic regression was used to analyze haplotype in LD block; and SNP set analysis based on logistic kernel machine regression was used to analyze pathway. All statistical analysis was performed by SPSS25.0, R2.14, Haploview4.2, SNPStats, and other statistical software packages. Results. In association analysis based on SNP, rs2167270 had statistical significance both in the adjusted and unadjusted covariate dominant model and in the unadjusted covariate overdominant model while it had no significant difference in the adjusted covariate overdominant model. Compared to GG genotype, rs2167270 of AG genotype had statistical significance in both the adjusted and unadjusted covariate codominant models. And rs16147 had statistical significance in robust test, stealth model and overdominant model, and adjusting and unadjusting covariate. This study found linkage disequilibrium existed between rs2167270 and rs4731426 of LEP, rs10889502 and rs17127107 of JAK1, rs2970847 and rs6821591 of PPARGC1A, rs249429 and rs3805486 of PRKAA1, rs1342382 and rs6588640 of PRKAA2, rs3766522 and rs6937 of PRKAB2, rs2970847 and rs6821591 of PRKAG2, and rs6436094 and rs645163 of PRKAG3. There was no positive finding with statistical significance from the unconditional logistic regression of the mentioned genes’ haplotype of LD block. Conclusions. LEP signaling pathway association with T2DM remained to be confirmed in Chinese Han population, although rs2167270 and rs16147 were significantly associated with T2DM.

1. Introduction

Diabetes has become a major public problem around the world [1], with the numbers increasing from 108 million in 1980 to 422 million in 2014 [2], and the number of people with diabetes is expected to increase to 592 million by 2035 [3]. Type 2 diabetes mellitus (T2DM) is the most prevalent form of diabetes mellitus (DM), accounting for about 90-95% of the total number of DM cases worldwide [4].

Insulin resistance is the basic pathophysiological change of type 2 diabetes, almost throughout the entire process of the occurrence and development of type 2 diabetes [5]. Studies have shown that free fatty acid (FFA) and adipocytokines released by adipose tissue can cause insulin resistance through a variety of ways [69].

For a long time, it has been believed that adipose tissue is only involved in the body’s energy storage and nontrembling thermogenesis. With the discovery of leptin (LEP) in 1994 [10], a series of new functions of adipose tissue were revealed. Adipose tissue is not only an energy storage organ of the body but also an endocrine organ [11]. It is known that human adipocytes secrete dozens of adipocytokines and protein factors, such as leptin, adiponectin, resistin, and visfatin. Some of these adipocytokines enter the human blood circulation to regulate the functions of distant target organs, and the other part plays a role of adjacent secretion or autocrine regulation in adjacent tissues or cells or within cells.

Leptin is a secreted protein synthesized by fat cells. It is the coded product of obese gene (ob gene). It exists freely in the blood or binds to leptin-binding protein. It has a wide range of centers such as reducing food intake and increasing energy consumption. And peripheral tissue effects are mediated through leptin receptors [12]. Leptin receptors belong to the class I cytokine receptor family and are widely distributed in central and peripheral tissues. Recent research results have shown that the occurrence and development of obesity, insulin resistance, and type 2 diabetes are closely related to the abnormal function of adipocytokines or their receptors [1315]. Leptin can cause insulin resistance and hyperinsulinemia through a variety of ways, and hyperinsulinemia and insulin resistance can cause excessive secretion of leptin, thus forming a vicious circle, leading to abnormal glucose and lipid metabolism [12, 16].

Single nucleotide polymorphism mainly refers to the polymorphism of DNA sequence caused by the variation of single nucleotide at the genome level. It is the most common type of human heritable variation. Recent studies have shown that there is an association between SNPs in multiple genes and T2DM and its complications. Letonja et al.’s [17] group finds a significant association between the SIRT1 rs7069102 polymorphism and diabetic nephropathy caused by T2DM. Zhang et al.’s [18] group discovers that rs62099905 in the CNDP1 gene was related to the serum glucose level of healthy Chinese Han population. Cui et al.’s [19] study shows that ADIPOQ gene polymorphisms rs10937273, rs1501299, rs182052, rs2241767, and rs266729 are associated with type 2 diabetes. Zhu et al.’s [20] study suggests that the interaction effects of SNP-SNP and SNP-environmental factors were related to T2DM susceptibility, and Khan et al.’s [21] research shows that the genetic variation of the TCF7L2 gene is related to the increased susceptibility to T2DM. Juttada et al.’s [22] group discovers that TCF7L2 gene polymorphism is closely related to a positive family history of diabetes. Bains et al.’s [23] study suggests that LEPrs7799039 and LEPRrs1137101 polymorphism conferred 3.4- and 2.1-fold risk towards the development of T2DM after adjustment of various covariates under a recessive genetic model.

This study intends to analyze the association between a SNP and type 2 diabetes based on the LEP signaling pathway, hoping to discover new mutation sites that may cause insulin resistance, conduct linkage disequilibrium analysis and haplotype-based association analysis, and analyze multiple loci of a single gene. Evaluate the distribution difference of haplotypes formed by points between the case group and the control group; conduct SNP-SNP interaction analysis to find meaningful interaction pairs; use pathway analysis to explore the relationship between LEP signaling pathway and T2DM.

2. Materials and Methods

2.1. Study Population

During November 2011 to October 2013, case-control groups of 1092 patients diagnosed as having T2DM were recruited from 10 hospitals in Zhanjiang, Maoming, Shaoguan, Dongguan, and Shenzhen of Guangdong Province. T2DM was diagnosed according to the criteria of the World Health Organization (WHO) in 1999. The inclusive criteria of the T2DM patient group were as follows: (1) range of ages 20-70; (2) random blood glucose  mmol/l with diabetes symptoms including polydipsia, polyphagia, polyuria, weight loss, itchiness, blurred vision, and other acute metabolic disorders caused by hyperglycemia and fasting plasma glucose  mmol/l without diabetes symptoms or blood glucose  mmol/l with a glucose tolerance test after two hours oral dose; and (3) patients without other serious diseases, such as cardiovascular and cerebrovascular diseases, malignant tumors, and other diseases. The control group comprised 1092 normal persons who were examined in the same hospital at the same time as the cases. The inclusive criteria of the control group were as follows: (1) range of ages 20-70; (2) without a family history of diabetes; and (3) healthy after physical examination including medical history, blood glucose, and other biochemical test results. This study was approved by the Ethics Committee. All the surveys and samples obtained the consents of participants in advance, and the informed consent forms had been legally consented. All the participants were the permanent residents of Chinese Han nationality in Guangdong Province, and there is no blood relationship among them.

2.2. Information Collection and Blood Sample Collection

A standardized survey is aimed at respondents who met the inclusion criteria, including age, gender, place of origin, occupation, disease history, duration of illness, smoking history, family history, complications, and diet, and exercise was conducted by a uniformly trained investigator. At the end of the investigation, the height and weight of the subjects were measured by a height and weight scale, and the blood pressure of the right arm in a sitting position was measured with a benchtop mercury sphygmomanometer, twice, and averaged.

5 ml of peripheral blood was collected by the endocrine nurses in the morning, and the relative clinical biochemical indicators, including fasting plasma glucose (FPG), total cholesterol (TC), triglyceride (TG), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and glycosylated hemoglobin A1c (HbA1c), were detected. The FPG level was measured by the glycosylation method, and the plasma levels of TC, TG, HDL-C, and LDL-C were determined in an enzymatic method. HBA1C was determined by a high-performance liquid chromatography (BIO-RADD-10TM glycated hemoglobin detection system), and other biochemical indicators were detected by an automatic biochemical analyzer. Besides, 4 ml of peripheral blood of the subjects (2 ml per tube) was collected, anticoagulated with EDTA·k2, and stored at -80°C.

This study was approved by the ethics committee, and all the investigations and sampling were subject to the consent of the subjects with signed informed consent.

2.3. Data Collation and Database Establishment

All completed questionnaires were uniformly coded, and all participants’ questionnaire information, physical examination, and clinical biochemical examination results were compiled. We used EpiData3.1 software to build a database with data by double input. The entered data was checked by both manual and computer methods to ensure that the data has no logic errors and no entry errors.

2.4. DNA Extraction

(1)Transfer the blood sample from the anticoagulation tube to a 10 ml centrifuge tube, and mark the centrifuge tube(2)Add 2 ml of cell lysate CL to blood containing anticoagulant, mix upside down, centrifuge at 5000 rpm for 10 min, and discard the supernatant(3)Add 3 ml of cell lysate CL to it, mix upside down, centrifuge at 5000 rpm for 10 min, discard the supernatant, and put the centrifuge tube upside down on clean absorbent paper for 1 min and dry by airing(4)Configure a mixture of buffer FG and proteinase K (, ready for use)(5)Add 1 ml of mixed solution to each sample and vortex immediately until the solution is free of clumps(6)Perform water bath at 65°C for 20 minutes, slowly tilting and inverting and mixing several times during this period(7)Add 1 ml of isopropanol, and mix thoroughly by inversion until filaments or clusters (DNA) appear(8)Centrifuge at 5000 rpm for 15 minutes, discard the supernatant, and place the centrifuge tube upside down on a clean absorbent paper to ensure that the precipitate is in the tube(9)Add 1 ml of 70% ethanol, transfer to a 1.5 ml centrifuge tube, centrifuge at 8000 rpm for 10 min in a small centrifuge, and discard the supernatant(10)Add 1 ml of 70% ethanol, centrifuge at 8000 rpm for 10 min in a small centrifuge, and discard the supernatant(11)Dry the DNA pellet in the centrifuge tube at room temperature until the liquid evaporates(12)Add the corresponding amount of buffer TB according to the amount of DNA; heat at 65°C for 1 hour to dissolve the DNA, during which time the finger flicks several times to aid the solution; and finally centrifuge for 30 s to shake the undissolved precipitate to the bottom of the tube(13)Perform DNA quality detection with a nucleic acid protein content detector, record DNA concentration (ng/μl) and purity (OD260/280 value, OD260/230 value), and finally store the DNA in a refrigerator at -80°C

2.5. Screening and Typing of SNP

The pathway diagram of LEP signaling pathway was obtained from the KEGG database (Kyoto Encyclopedia of Genes and Genomes), and 15 main genes were selected, including LEP, LEPR, JAK (JAK1, JAK2), STAT3, SHP-2 (PTPN11), PPARGC1Alpha (PPARGC1A), AMPK (PRKAA1, PRKAA2, PRKAB2, PRKAG1, PRKAG2, and PRKAG3), alpha-MSH (POMC), and NPY. Haploview (ver.4.2) was used to select tagSNPs from selected genes and their upstream and downstream 5 kb regions, requiring and (SNP MAF and LD data come from the HapMap database, release 27, http://snp.cshl.org/). Then, use FastSNP (http://fastsnp.ibms.sinica.edu.tw/pages/input_CandidateGeneSearch.jsp) to predict the potential function of tagSNP, select tagSNP with higher prediction score, and select 1-2 tagSNP for each gene. Finally, 23 SNPs were selected from 15 genes. These SNPs were distributed on 8 chromosomes, one of which was located downstream of the 3 end, 4 of which were located in the untranslated region of the 3 end, 2 of which were located upstream of the 5 end, one of which was located in the 5untranslated region, 13 of which were located in the intron region, and 2 of which were located in the exon region. In the Chinese population, the minimum MAF was 0.128 (rs3766522) and the maximum MAF was 0.444 (rs6436094). Related information is listed in Table 1.

2.6. Statistical Analysis

Comparisons of all variables between T2DM cases and control subjects were carried out by the test for categorical variables or -test for continuous variables. The Pearson chi-square test, Cochran-Armitage trend test, MAX3, and logistic regression were used to analyze the association between SNP and T2DM; unconditional logistic regression was used to analyze haplotype in LD block; and SNP set analysis based on logistic kernel machine regression was used to analyze pathway. All statistical analysis was performed by SPSS25.0, PLINK 1.07, R 2.14.2, Haploview 4.2, SNPStats, and other statistical software packages.

3. Results

3.1. Equilibrium Test Results of Baseline Data in the Case Group and Control Group

In this study, 1092 cases and 1092 controls were included, and the individuals with the loss rate of SNP were excluded. Eventually, 1067 cases and 1054 controls were included in the following analysis, including 532 males and 535 females in the case group with an average age of 49.86 years and 532 males and 522 females in the control group with an average age of 50.47 years. Differences of the age composition, body mass index (BMI), fasting plasma glucose (FPG), triglyceride level, and low-density cholesterol level between the two groups were statistically significant. The average age, BMI, FPG, and triglyceride level of the case group are higher than those of the control group, as shown in Table 2.

3.2. SNP Typing Results

The typing success rate of 23 SNPs was all above 98%, and the minimum allele frequency was 0.119, and the maximum was 0.474. The Hardy-Weinberg equilibrium test shows that each point satisfies the Hardy-Weinberg equilibrium, as shown in Table 3.

3.3. Allele Association Analysis Results

The results of genotype association analysis are shown in Table 4. There was no statistical difference in the suballelic frequency of each SNP between the case group and the control group (). After adjusting for covariates such as age and BMI, there was no statistical difference in the suballelic frequency between the case group and the control group.

3.4. Genotype Association Analysis Results

The genotype (CC/CT/TT) distribution of rs16147 was statistically significant between the two groups (, ). There was no statistical difference in genotype distribution of other SNPs, as shown in Table 5.

In order to further confirm whether each SNP was associated with T2DM and whether the incidence of disease increases with the increase of risk alleles in the genotype, we conducted the Cochran-Armitage trend test under different genetic models (additive model, codominant model, dominant model, invisible model, and superdominant model), and the results are shown in Table 6.

Since the trend test relies on prespecified scores, different scores correspond to different genetic models. Only when the designated gene model is a true gene model the corresponding trend test is the test with the best power. We further applied the robust test method and compared the results obtained with the methods based on various genetic models (see Table 7 for details). The result of rs16147 was statistically significant (). The other SNPs were not statistically significant in the robust test results.

After adjusting the covariate, the results under the five genetic models are shown in Table 8. rs2167270 still has statistical significance in the dominant model (, , 95% CI: 1.01 1.44), while there is no statistical significance in the superdominant model; in the codominant model, compared to GG genotype, AG genotype (, 95% CI: 1.01 1.47, ) is statistically significant, and compared to GG genotype, AA genotype (, 95% CI: 0.78 1.63, ) is still not statistically significant. rs16147 was still statistically significant in both the stealth model and the superdominant model, in the recessive model (, , 95% CI: 0.56-0.93), and in the superdominant model (, , 95% CI: 1.01-1.44).

3.5. Linkage Disequilibrium Analysis and Association Analysis Based on Haplotype

Linkage is a genetic tendency where genetic markers are inherited together as a result of being near to one another on the same chromosome. Genetic linkage analysis, one of the old study approaches, focuses on genomic regions with large genetic effect that can influence the development of a disease [24, 25].

Haploview 4.2 software was used to analyze the linkage imbalance (LD) between different loci in the same gene, and the results showed that there was linkage imbalance between loci in 8 genes including LEP. Figure 1 shows the composition of LD block of 8 genes including LEP.

The haplotype in LD block was analyzed by using SNPStats software. The analysis results are shown in Table 9. The haplotype of LD block in each gene was not found to be statistically significant.

3.6. SNP-SNP Interaction Results

We analyzed 23 SNPS of 15 genes under the pathways of the LEP pathway by logistic regression model-based SNP-SNP interaction analysis, and the partial results are shown in Table 10. According to the test level of 0.01, 2 pairs of statistically significant results were obtained. The OR value of the interaction term is 1.311, which corresponds to the interaction between rs16147 and rs1044471, and the corresponding genes are NPY and ADIPOR2. The OR value of the interaction term is 1.422, which corresponds to the interaction between rs2167270 and rs5435, and the corresponding genes are LEP and SLC2A4.

We uploaded 8 genes from the LEP signaling pathway to the STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) tool. The interaction between proteins encoded by these genes was analyzed, and the results are shown in Figure 2.

3.7. Pathway Analysis Result

Four kinds of kernel functions linear, linear.weighted, IBS, and IBS.weighted were used to analyze the SNP set based on pathways. The results are shown in Table 11. No statistical significance was found whether covariates were added or not (), nor was the empirical P value obtained by bootstrap method.

4. Discussion

The prevalence of diabetes increases significantly in recent decades, affecting about 6% of adult population globally. Therefore, it is one of the major health care challenges in the world [26]. In China, the incidence of T2D in adults has been increasing over recent decades [27]. T2DM is a major chronic disease that is affected by genetic and environmental factors [28]. The development of T2DM involves multiple factors and signaling pathways, such as hyperglycemia, dyslipidemia, and oxidative damage. The main pathological features are insulin resistance and damage to the function of cells, which is a chronic metabolic disease caused by the interaction of genetic and environmental factors [29, 30].

Heredity has all along been mentioned to play an important role in the development of diabetes [31]. Genome-wide association studies have identified more than 100 independent SNPS associated with T2DM regulation and risk [32]. Several gene polymorphisms of leptin have been reported to be associated with T2DM [2, 33]. For example, the leptin (LEP) G2548A polymorphism has been associated with increased leptin production and plasma secretion from adipocytes [34]. Human LEP is located on chromosome7q31.3, and its translational product is leptin, which plays a decisive role in the regulation of human appetite and results in severe metabolic disorders [35]. Such signalization upregulates and downregulates genes involved in cell differentiation, proliferation, apoptosis, and synthesis of extracellular matrix proteins [36].

In this study, rs2167270 and rs16147 were associated with type 2 diabetes in the unit-point correlation analysis. However, in order to further confirm its association with type 2 diabetes, independent sample verification is needed in Chinese Han population. Linkage disequilibrium analysis and correlation analysis based on single type times get within the eight genes of the pathway between the two respective loci; there is a linkage disequilibrium relationship; each genotype in LD block of haploid analysis has not been having a positive result. The LEP pathway gene SNP-SNP interactions analysis results have statistical significance and illustrate that the SNP-SNP interactions (gene interactions) are part of the genetic structure of type 2 diabetes mellitus. LEP gene SNP in pathway, the statistically significant results in the occurrence of type 2 diabetes, one of more than just a channel or certain genes, plays an important role, may also involve two pathways even based on the joint action of multiple pathways between LEP pathways of SNP set correlation analysis, and has statistically significant results. Considering that a few SNPS were selected in each gene, it is necessary to increase the number of SNPS in follow-up studies.

However, the sample of this study is only the Han population in Guangdong Province. The research group has certain limitations. It is necessary to further verify in other populations to provide new clues for the prevention and treatment strategy of T2DM.

5. Conclusion

In summary, we found that rs2167270 and rs16147 may interact to affect the risk of T2DM but remained to further confirm their association, and independent sample validation is required in the Chinese Han population. According to the logistic kernel machine regression results, we cannot think the LEP signaling pathway is associated with T2DM.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declared no conflicts of interest.

Authors’ Contributions

All the authors made substantial contributions to the conception and design. Haibing Yu, Lin Xu, Hao Liu, Jialu Huang, Yuanlin Ding, and Danli Kong collected the data. Ling Luo, Rong Chen, Jian xu, Chunwen Lin, Weiying Chen, Yuxun Xie, and Huihuang Yang made substantial contributions to the analysis and interpretation of data. Haibing Yu wrote the first draft of the manuscript. All the authors were involved in revising the manuscript critically for important intellectual content. Yuanlin Ding and Danli Kong gave the final approval of the version to be published. All the authors read and approved the final manuscript.

Acknowledgments

We appreciate all authors for their contributions and physicians and participants. The study is supported by the Discipline Construction Project of Guangdong Medical University (Fund No. 4SG21264P), the Basic and Applied Basic Research Foundation of Guangdong Province Regional Joint Fund Project (the Key Project, Fund No. 2020B1515120021), the 2021 Guangdong Provincial Medical Research Fund Project (Fund No. A2021395), the Natural Science Foundation of Guangdong Province Basic and Applied Basic Research Fund (Fund No. 2021A1515010061), the 2020 Zhanjiang City Science and Technology Development Special Fund Competitive Allocation Project (Fund No. 2020A01031), the 2020 Guangdong Medical University Scientific Research Fund Natural Science Key Cultivation Project (Fund No. GDMUZ2020008), the 2019 Guangdong Medical University College Student Innovation Experimental Project (Key Project, Fund No. ZZDG003), the 2019 Characteristic Innovation Projects (Natural Science) of Ordinary Universities in Guangdong Province (Fund No. 2019KTSCX047), the 2020 Dongguan Science and Technology Commissioner Project, and the Young Innovative Talents Projects in Ordinary Universities in Guangdong Province (Fund No. 2018KQNCX088).