The prognostic values and immune characteristics of polo-like kinases (PLKs) family: A pan-cancer multi-omics analysis.

Background: In the realm of tumor-targeted therapeutics, Polo-like kinases (PLKs) are a significant group of protein kinases that were found recently as being related to tumors. However, the significance of PLKs in pan-cancer remains systematically studied. Methods and materials: We integrated multi-omics data to comprehensively investigate the expression patterns of the PLK family across various cancer types. Subsequently, study examined the associations between tumor mutation burden (TMB), microsatellite instability (MSI), immune subtype classification, immune infiltration, tumor microenvironment scores, immune checkpoint gene expression, and the PLKs expression profiles within various tumor types. Furthermore, using our mRNA sequencing data (TRUCE01) and four bladder cancer (BLCA) cohorts (GSE111636, GSE176307, and IMvigor210), We examined the correlation between the expression level of PLK and immunotherapy effectiveness. Next, Gene set enrichment analysis (GSEA) was evaluated to find that potentially enriched PLK signaling pathways. Utilizing TIMER 2.0, we conducted an immune infiltration analysis underlying transcriptome expression, copy number variations (CNV), or somatic mutations of PLKs in BLCA. Finally, mRNA expression validation of PLK1/3/4 by real-time PCR within 10 paired BLCA tissues, protein expression verification through the Human Protein Atlas (HPA), and PLK4 in vitro cytological studies have been employed in BLCA. Results: The expression of most of the PLK family members exhibits variation between cancerous tissues and adjacent normal tissues across various cancer species. Furthermore, the expression of PLKs demonstrates a significant association with immunotyping, infiltration of immune cell, tumor mutational burden (TMB), microsatellite instability (MSI), immunological checkpoint gene activity and therapeutic effectiveness in pan-tumor tissues. Additional investigation into the correlation between the PLK family and BLCA has revealed that the expression of the PLK genes holds substantial significance in the biological processes of BLCA. Furthermore, a notable association has been observed between the copy number variation, variant status, and the degree of certain immunological cell infiltration. Of note, the expression validation and in vitro phenotypic experiments have demonstrated that PLK4 has a significant function in promoting the BLCA cell proliferation, migration, and invasion. Conclusion: Collectively, based on various databases, our results highlight the involvement of PLK gene family in the formation of different types of tumors and identify PLK-related genes that may be used for therapy.

Identification and validation of a dysregulated TME-related gene signature for predicting prognosis, and immunological properties in bladder cancer.

Background: During tumor growth, tumor cells interact with their tumor microenvironment (TME) resulting in the development of heterogeneous tumors that promote tumor occurrence and progression. Recently, there has been extensive attention on TME as a possible therapeutic target for cancers. However, an accurate TME-related prediction model is urgently needed to aid in the assessment of patients' prognoses and therapeutic value, and to assist in clinical decision-making. As such, this study aimed to develop and validate a new prognostic model based on TME-associated genes for BC patients. Methods: Transcriptome data and clinical information for BC patients were extracted from The Cancer Genome Atlas (TCGA) database. Gene Expression Omnibus (GEO) and IMvigor210 databases, along with the MSigDB, were utilized to identify genes associated with TMEs (TMRGs). A consensus clustering approach was used to identify molecular clusters associated with TMEs. LASSO Cox regression analysis was conducted to establish a prognostic TMRG-related signature, with verifications being successfully conducted internally and externally. Gene ontology (GO), KEGG, and single-sample gene set enrichment analyses (ssGSEA) were performed to investigate the underlying mechanisms. The potential response to ICB therapy was estimated using the Tumor Immune Dysfunction and Exclusion (TIDE) algorithm and Immunophenoscore (IPS). Additionally, it was found that the expression level of certain genes in the model was significantly correlated with objective responses to anti-PD-1 or anti-PD-L1 treatment in the IMvigor210, GSE111636, GSE176307, or Truce01 (registration number NCT04730219) cohorts. Finally, real-time PCR validation was performed on 10 paired tissue samples, and in vitro cytological experiments were also conducted on BC cell lines. Results: In BC patients, 133 genes differentially expressed that were associated with prognosis in TME. Consensus clustering analysis revealed three distinct clinicopathological characteristics and survival outcomes. A novel prognostic model based on nine TMRGs (including C3orf62, DPYSL2, GZMA, SERPINB3, RHCG, PTPRR, STMN3, TMPRSS4, COMP) was identified, and a TMEscore for OS prediction was constructed, with its reliable predictive performance in BC patients being validated. MultiCox analysis showed that the risk score was an independent prognostic factor. A nomogram was developed to facilitate the clinical viability of TMEscore. Based on GO and KEGG enrichment analyses, biological processes related to ECM and collagen binding were significantly enriched among high-risk individuals. In addition, the low-risk group, characterized by a higher number of infiltrating CD8+ T cells and a lower burden of tumor mutations, demonstrated a longer survival time. Our study also found that TMEscore correlated with drug susceptibility, immune cell infiltration, and the prediction of immunotherapy efficacy. Lastly, we identified SERPINB3 as significantly promoting BC cells migration and invasion through differential expression validation and in vitro phenotypic experiments. Conclusion: Our study developed a prognostic model based on nine TMRGs that accurately and stably predicted survival, guiding individual treatment for patients with BC, and providing new therapeutic strategies for the disease.

Construction and Validation of a Prognostic Model Based on Pyroptosisrelated Genes in Bladder Cancer.

BACKGROUND: Bladder cancer (BCa) is a highly prevalent disease with a poor prognosis. There is no better forecasting method for it yet. Current studies demonstrate that pyroptosis is involved in the development and progression of various cancers. METHODS: This study employed bioinformatics techniques to analyze the data of BCa patients obtained from the TCGA and GEO databases in order to construct a prognostic risk model. The TCGA dataset was used for the training set, and the multiple external datasets (including GSE13507, GSE31684, GSE48075, IMvigor210, and GSE32894) were applied as the validation sets. Prognostic-associated pyroptosis genes screened by univariate Cox regression analysis were utilized to construct the lasso Cox regression model. GO and KEGG analysis results identified the selected genes that are primarily involved in the inflammation and cell death processes. The related patients were grouped into low- and high-risk groups. Kaplan-Meier survival analysis was performed to compare survival differences between the risk groups. The accuracy of this risk prediction model was assessed by ROC. We also applied the Human Protein Atlas (HPA) to detect the protein expression of these genes. Subsequently, qRT-PCR was performed to verify the expression of these model genes. RESULTS: There are 29 pyroptosis-related genes with significant expression differences between BCa and corresponding adjacent tissues, and 11 genes (SH2D2A, CHMP4C, MRFAP1L1, GBP2, EHBP1, RAD9A, ANXA1, TMEM109, HEYL, APOL2, ORMDL1) were picked by univariate and LASSO Cox regression analysis. Immunological cell infiltration and ssGSEA results further indicated that the low and high-risk groups were substantially correlated with the immune status of BCa patients. According to TCGA and multiple external datasets, Kaplan-Meier survival curves showed the overall survival rate of the high-risk group to be decreased. ROC curves showed the model established to be accurate and reliable. Moreover, the HPA database also demonstrated the verification of the modeled genes' expression in BCa and normal bladder tissue using the HPA database. qRT-PCR results also suggested the up-regulated EHBP1 and down-regulated RAD9A mRNA expression levels to be confirmed in 15 pairs of BCa and corresponding adjacent tissues. CONCLUSION: This study presents the development and validation of a novel gene signature associated with pyroptosis, which holds the potential for predicting patient outcomes in BCa and providing insights into the immune microenvironment of BCa.

Construction and validation of a metabolism-associated gene signature for predicting the prognosis, immune landscape, and drug sensitivity in bladder cancer.

Tumor Metabolism is strongly correlated with prognosis. Nevertheless, the prognostic and therapeutic value of metabolic-associated genes in BCa patients has not been fully elucidated. First, in this study, metabolism-related differential expressed genes DEGs with prognostic value in BCa were determined. Through the consensus clustering algorithm, we identified two molecular clusters with significantly different clinicopathological features and survival prognosis. Next, a novel metabolism-related prognostic model was established. Its reliable predictive performance in BCa was verified by multiple external datasets. Multivariate Cox analysis exhibited that risk score were independent prognostic factors. Interestingly, GSEA enrichment analysis of GO, KEGG, and Hallmark gene sets showed that the biological processes and pathways associated with ECM and collagen binding in the high-risk group were significantly enriched. Notely, the model was also significantly correlated with drug sensitivity, immune cell infiltration, and immunotherapy efficacy prediction by the wilcox rank test and chi-square test. Based on the 7 immune infiltration algorithm, we found that Neutrophils, Myeloid dendritic cells, M2 macrophages, Cancer-associated fibroblasts, etc., were more concentrated in the high-risk group. Additionally, in the IMvigor210, GSE111636, GSE176307, or our Truce01 (registration number NCT04730219) cohorts, the expression levels of multiple model genes were significantly correlated with objective responses to anti-PD-1/anti-PD-L1 immunotherapy. Finally, the expression of interested model genes were verified in 10 pairs of BCa tissues and para-carcinoma tissues by the HPA and real-time fluorescent quantitative PCR. Altogether, the signature established and validated by us has high predictive power for the prognosis, immunotherapy responsiveness, and chemotherapy sensitivity of BCa.

Withdrawal notice to: "Pan-cancer analyses of clinical prognosis, immune infiltration, and immunotherapy efficacy for TRPV family using multi-omics data" [Heliyon 9, (2023) 16897].

[This corrects the article DOI: 10.1016/j.heliyon.2023.e16897.].

Pan-cancer analyses of clinical prognosis, immune infiltration, and immunotherapy efficacy for TRPV family using multi-omics data.

Background: Transient receptor potential cation channel subfamily V (TRPV) play an essential in cancer initiation, progression, and treatment. TRPV expression alteration are shown relate to multiple cancers prognosis and treatment of cancers but are less-studied in pan-cancer. In this study, we characterize the clinical prediction value of TRPV at pan-cancer level. Methods: Several databases were used to examine the transcript expression difference in tumor vs. normal tissue, copy-number variant (CNV) and single nucleotide polymorphisms (SNP) mutation of each TRPV members in pan-cancer, including The Cancer Genome Atlas (TCGA) and cBioPortal. We performed K-M survival curve and univariate Cox regression analyses to identify survival and prognosis value of TRPV. CellMiner were selected to explore drug sensitivity. We also analyzed association between tumor mutation burden (TMB), microsatellite instability (MSI), tumor immune microenvironment and TRPV family genes expression. Moreover, we investigated the relationship between TRPVs expression and effectiveness of immunotherapy in multiple cohorts, including one melanoma (GSE78220), one renal cell carcinoma (GSE67501), and three bladder cancer cohorts (GSE111636, IMvigor210, GSE176307 and our own sequencing dataset (TRUCE-01)), and further analyzed the changes of TRPVs expression before and after treatment (tislelizumab combined with nab-paclitaxel) of bladder cancer. Next, we made a special effort to investigate and study biological functions of TRPV in bladder cancer using gene set enrichment analysis (GSEA), and conducted immune infiltration analysis with TRPVs family genes expression, copy number or somatic mutations of bladder cancer by TIMER 2.0. Finally, real-time PCR and protein expression validation of TRPVs within 10 paired cancer and para-carcinoma tissue samples, were also performed in bladder cancer. Results: Only TRPV2 expression was lower in most cancer types among TRPV family genes. All TRPVs were correlated with survival changes. Amplification was the significant gene alternation in all TRPVs. Next, analysis between TRPVs and clinical traits showed that TRPVs were related to pathologic stage, TNM stage and first course treatment outcome. Moreover, TRPV expression was highly correlated with MSI and TMB. Immunotherapy is a research hotspot at present, our result showed the significant association between TRPVs expression and immune infiltration indicated that TRPV expression alternation could be used to guide prognosis. In addition, we also discovered that the expression level of TRPV1/2/3/4/6 was positively or negatively correlated with objective responses to anti-PD-1/PD-L1 across multiple immunotherapy cohort. Further analysis of drug sensitivity showed the value to treatment. Based on the above analysis, we next focused on TRPV family in bladder cancer. The result demonstrated TRPV also played an important role in bladder cancer. Finally, qPCR assay verified our analysis in bladder cancer. Conclusion: Our study firstly revealed expression and genome alternation of TRPV in pan-cancer. TRPV could be used to predict prognosis or instructing treatment of human cancers, especially bladder cancer.

IER3IP1 is critical for maintaining glucose homeostasis through regulating the endoplasmic reticulum function and survival of ß cells.

Recessive mutations in IER3IP1 (immediate early response 3 interacting protein 1) cause a syndrome of microcephaly, epilepsy, and permanent neonatal diabetes (MEDS). IER3IP1 encodes an endoplasmic reticulum (ER) membrane protein, which is crucial for brain development; however, the role of IER3IP1 in ß cells remains unknown. We have generated two mouse models with either constitutive or inducible IER3IP1 deletion in ß cells, named IER3IP1-ßKO and IER3IP1-ißKO, respectively. We found that IER3IP1-ßKO causes severe early-onset, insulin-deficient diabetes. Functional studies revealed a markedly dilated ß-cell ER along with increased proinsulin misfolding and elevated expression of the ER chaperones, including PDI, ERO1, BiP, and P58IPK. Islet transcriptome analysis confirmed by qRT-PCR revealed decreased expression of genes associated with ß-cell maturation, cell cycle, and antiapoptotic genes, accompanied by increased expression of antiproliferation genes. Indeed, multiple independent approaches further demonstrated that IER3IP1-ßKO impaired ß-cell maturation and proliferation, along with increased condensation of ß-cell nuclear chromatin. Inducible ß-cell IER3IP1 deletion in adult (8-wk-old) mice induced a similar diabetic phenotype, suggesting that IER3IP1 is also critical for function and survival even after ß-cell early development. Importantly, IER3IP1 was decreased in ß cells of patients with type 2 diabetes (T2D), suggesting an association of IER3IP1 deficiency with ß-cell dysfunction in the more-common form of diabetes. These data not only uncover a critical role of IER3IP1 in ß cells but also provide insight into molecular basis of diabetes caused by IER3IP1 mutations.

Effect of artemisinin on improving islet function in rats fed with maternal high-fat diet.

BACKGROUND: Maternal high-fat diet (HFD) is a detrimental factor in developing glucose intolerance, obesity, and islet dysfunction. However, the effect of artemisinin on maternal HFD and whether it is related to the alterations of islet function is seldom studied since artemisinin treatments not only attenuate insulin resistance (IR) and restore islet ß cell function in Diabetes mellitus type 2. METHODS: Female rats were randomly fed a HFD (45% kcal from fat), HFD + artemisinin, or a regular chow diet (RCD) before pregnancy and during gestation. Glucose metabolism and the ß cell phenotypes were assessed. RESULTS: Maternal HFD increased islet load in female rats, proliferation of pancreatic ß cells, increased insulinogen, and decreased insulin secretion response to high glucose stimulation with delayed insulin release, increased fasting glucose, and glucose area under the curve compared with the general diet group. HFD inhibited expression of Foxo1 and PAX6 in female rats. Under the effect of both HFD and pregnancy, islet load was further increased, insulinogen was further increased, and fasting insulin level and fasting glucose were higher than RCD fed general-pregnancy group. ALDH1a3 transdifferentiation and PAX6, Foxo1, and PDX1 expression were increased in islets of high-fat pregnant rats. When adding artemisinin in HFD treated pregnant rats, islet function was significantly improved. CONCLUSIONS: Intervention with artemisinin in maternal HFD resulted in reduced islet size, decreased number of ß-cells and improved islet microcirculation, insulin processing shear process, decreased insulinogen/insulin ratio, and restored islet function through increased expression of PC1/3.

A Novel Nonsense INS Mutation Causes Inefficient Preproinsulin Translocation Into the Endoplasmic Reticulum.

Preproinsulin (PPI) translocation across the membrane of the endoplasmic reticulum (ER) is the first and critical step of insulin biosynthesis. Inefficient PPI translocation caused by signal peptide (SP) mutations can lead to ß-cell failure and diabetes. However, the effect of proinsulin domain on the efficiency of PPI translocation remains unknown. With whole exome sequencing, we identified a novel INS nonsense mutation resulting in an early termination at the 46th residue of PPI (PPI-R46X) in two unrelated patients with early-onset diabetes. We examined biological behaviors of the mutant and compared them to that of an established neonatal diabetes causing mutant PPI-C96Y. Although both mutants were retained in the cells, unlike C96Y, R46X did not induce ER stress or form abnormal disulfide-linked proinsulin complexes. More importantly, R46X did not interact with co-expressed wild-type (WT) proinsulin in the ER, and did not impair proinsulin-WT folding, trafficking, and insulin production. Metabolic labeling experiments established that, despite with an intact SP, R46X failed to be efficiently translocated into the ER, suggesting that proinsulin domain downstream of SP plays an important unrecognized role in PPI translocation across the ER membrane. The study not only expends the list of INS mutations associated with diabetes, but also provides genetic and biological evidence underlying the regulation mechanism of PPI translocation.

Biological behaviors of mutant proinsulin contribute to the phenotypic spectrum of diabetes associated with insulin gene mutations.

Insulin gene mutation is the second most common cause of neonatal diabetes (NDM). It is also one of the genes involved in maturity-onset diabetes of the young (MODY). We aim to investigate molecular behaviors of different INS gene variants that may correlate with the clinical spectrum of diabetes phenotypes. In this study, we concentrated on two previously uncharacterized MODY-causing mutants, proinsulin-p.Gly44Arg [G(B20)R] and p.Pro52Leu [P(B28)L] (a novel mutant identified in one French family), and an NDM causing proinsulin-p.(Cys96Tyr) [C(A7)Y]. We find that these proinsulin mutants exhibit impaired oxidative folding in the endoplasmic reticulum (ER) with blocked ER export, ER stress, and apoptosis. Importantly, the proinsulin mutants formed abnormal intermolecular disulfide bonds that not only involved the mutant proinsulin, but also the co-expressed WT-proinsulin, forming misfolded disulfide-linked proinsulin complexes. This impaired the intracellular trafficking of WT-proinsulin and limited the production of bioactive mature insulin. Notably, although all three mutants presented with similar defects in folding, trafficking, and dominant negative behavior, the degrees of these defects appeared to be different. Specifically, compared to MODY mutants G(B20)R and P(B28)L that partially affected folding and trafficking of co-expressed WT-proinsulin, the NDM mutant C(A7)Y resulted in an almost complete blockade of the ER export of WT-proinsulin, decreasing insulin production, inducing more severe ER stress and apoptosis. We thus demonstrate that differences in cell biological behaviors among different proinsulin mutants correlate with the spectrum of diabetes phenotypes caused by the different INS gene mutations.

Misfolded proinsulin impairs processing of precursor of insulin receptor and insulin signaling in ß cells.

Insulin resistance in classic insulin-responsive tissues is a hallmark of type 2 diabetes (T2D). However, the pathologic significance of ß-cell insulin resistance and the underlying mechanisms contributing to defective insulin signaling in ß cells remain largely unknown. Emerging evidence indicates that proinsulin misfolding is not only the molecular basis of mutant INS-gene-induced diabetes of youth (MIDY) but also an important contributor in the development and progression of T2D. However, the molecular basis of ß-cell failure caused by misfolded proinsulin is still incompletely understood. Herein, using Akita mice expressing diabetes-causing mutant proinsulin, we found that misfolded proinsulin abnormally interacted with the precursor of insulin receptor (ProIR) in the endoplasmic reticulum (ER), impaired ProIR maturation to insulin receptor (IR), and decreased insulin signaling in ß cells. Importantly, using db/db insulin-resistant mice, we found that oversynthesis of proinsulin led to an increased proinsulin misfolding, which resulted in impairments of ProIR processing and insulin signaling in ß cells. These results reveal for the first time that misfolded proinsulin can interact with ProIR in the ER, impairing intracellular processing of ProIR and leading to defective insulin signaling that may contribute to ß-cell failure in both MIDY and T2D.-Liu, S., Li, X., Yang, J., Zhu, R., Fan, Z., Xu, X., Feng, W., Cui, J., Sun, J., Liu, M. Misfolded proinsulin impairs processing of precursor of insulin receptor and insulin signaling in ß cells.

Age- and Sex-Associated Impacts of Body Mass Index on Stroke Type Risk: A 27-Year Prospective Cohort Study in a Low-Income Population in China.

The relationship between body mass index (BMI) and stroke type has remained controversial despite studies demonstrating that BMI is related to stroke risk, especially in specific groups. We assessed the age- and sex-associated impacts of BMI on stroke type in a low-income, poorly educated population in China. The association of BMI with stroke type was estimated using Cox regression analyses in this prospective cohort study, after adjusting for sex, age, education level, hypertension, diabetes, smoking, and alcohol drinking status. During the follow-up period, 638 stroke cases occurred among the 3,906 participants included in this prospective study. For men aged <65 years, being overweight was an independent predictor of all stroke subtypes, compared with normal-weight individuals; the hazard ratios (HRs) and 95% confidence intervals (CIs) were 1.98 (1.52-2.58) for total stroke, 1.69 (1.22-2.33) for ischemic stroke, and 3.62 (2.09-6.25) for hemorrhagic stroke, all P < 0.001. Being underweight was also an independent predictor of hemorrhagic stroke (HR, 5.10; 95%CI, 1.80-14.50, P = 0.002). For women <65-years-old, being overweight was a risk factor for total (HR, 1.38; 95% CI, 1.01-1.89; P = 0.044) and hemorrhagic strokes (HR, 2.06; 95% CI, 1.00-4.28; P = 0.050); obesity was a risk factor for total (HR, 2.47; 95% CI, 1.60-3.82) and ischemic strokes (HR, 2.53; 95% CI, 1.54-4.15), all P < 0.001. These findings suggest that weight management should be a high priority for substantially reducing the heavy burden of strokes in rural China among both men and women <65-years-old; men<65-years-old should maintain their weight within a reasonable range.

Determinants of Developing Stroke Among Low-Income, Rural Residents: A 27-Year Population-Based, Prospective Cohort Study in Northern China.

Although strokes are the leading cause of death and disability in many countries, China still lacks long-term monitoring data on stroke incidence and risk factors. This study explored stroke risk factors in a low-income, rural population in China. The study population was derived from the Tianjin Brain Study, a population-based stroke monitoring study that began in 1985. This study documented the demographic characteristics, past medical histories, and personal lifestyles of the study participants. In addition, physical examinations, including measurements of blood pressure (BP), height, and weight, were performed. Hazard ratios (HRs) were estimated for the risk factors for all subtypes of stroke using multivariate Cox regression analyses. During the study with mean following-up time of 23.16 years, 3906 individuals were recruited at baseline, and during 27 years of follow-up, 638 strokes were documented. The multivariate Cox regression analyses revealed a positive correlation between age and stroke incidence. Limited education was associated with a 1.9-fold increase in stroke risk (lowest vs. highest education level). Stroke risk was higher among former smokers than among current smokers (HR, 1.8 vs. 1.6; both, P < 0.05). Moreover, stroke risk was significantly associated with sex (HR, 1.8), former alcohol drinking (HR, 2.7), baseline hypertension (HR, 3.1), and overweight (HR, 1.3). In conclusion, this study identified uncontrollable (sex and age) and controllable (education, smoking, alcohol drinking, hypertension, and overweight) risk factors for stroke in a low-income, rural population in China. Therefore, it is critical to control BP and weight effectively, advocate cessation of smoking/alcohol drinking, and enhance the education level in this population to prevent increase in the burden of stroke in China.

Association of the Transcription Factor 7 Like 2 (TCF7L2) Polymorphism With Diabetic Nephropathy Risk: A Meta-Analysis.

It is assumed that genetic factors may participate in the development of diabetic nephropathy (DN). The association between TCF7L2 gene polymorphism and DN risk is still unclear. To evaluate the relationship, we performed this meta-analysis. Eligible relevant studies were searched and selected from PubMed, Embase, and ISI Web of Science. Summary effect estimates were derived using a random effects model, with attention to study quality and publication bias. Ethnical approval was not necessary, because this meta-analysis was based on published articles, and did not involve patient consent. A total of 7 studies were identified. Analysis of all studies indicated significant association between TCF7L2 gene polymorphism and DN risk (odds ratio [OR] = 1.31, 95% confidence interval (CI) = 1.10-1.56, Pheterogeneity < 0.00001, P = 0.002). Subgroup analysis showed similar results in Asian (OR = 1.33, 95% CI = 1.10-1.62, Pheterogeneity = 0.03, P = 0.004), in Caucasian (OR = 2.27, 95% CI = 1.78-2.90, Pheterogeneity = 0.17, P < 0.00001), in rs7903146 mutation (OR = 1.61, 95% CI = 1.25-2.07, Pheterogeneity < 0.00001, P = 0.0002), However, no association was observed in Negroid (OR = 1.10, 95% CI = 0.90-1.35, Pheterogeneity < 00001, P = 0.36). Our results suggest that TCF7L2 gene polymorphism may contribute to the risk of DN. However, more studies should be launched in the future.