Abstract
Background:
Neutrophil extracellular traps (NETs) cause delayed wound closed up in type 2 diabetes mellitus (T2DM), but the specific regulatory mechanism of NETs-related genes (NETs-RGs) in T2DM is unclear.
Methods:
We acquired GSE21321 and GSE15932 datasets from gene expression omnibus (GEO) database. First, differentially expressed genes (DEGs) between T2DM and control samples of GSE21321 dataset were sifted out by differential expression analysis. NETs scores were calculated for all samples in GSE21321 dataset, and key module genes associated with NETs scores were screened by constructing co-expression network. Then, DEGs and key module genes were intersected to yield intersection genes, and candidate genes were identified by constructing a protein protein interaction (PPI) network. Least absolute shrinkage and selection operator (LASSO) regression analysis was implemented on candidate genes to screen out diagnostic genes, and they were subjected to single sample gene set enrichment analysis (ssGSEA). Finally, immune characteristic analysis was carried out, and we constructed the gene-drug and transcription factor (TF)-miRNA-mRNA networks. Besides, we validated the expression of diagnostic genes by quantitative real-time polymerase chain reaction (qRT-PCR).
Results:
In total, 23 candidate genes were gained by PPI analysis. The 5 diagnostic genes, namely, inter-trypsin inhibitor heavy chain 3 (ITIH3), fibroblast growth factor 1 (FGF1), neuron cell adhesion molecule (NRCAM), advanced glycosylation end-product-specific receptor (AGER), and calcium voltage-gated channel subunit alpha1 C (CACNA1C), were identified via LASSO analysis, and they were involved in carboxylic acid transport, axonogenesis, etc. M2 Macrophage, Monocyte, Natural killer (NK) cell, and Myeloid dendritic cells (DC) were remarkably different between T2DM and control samples. Diagnostic genes had the strongest and the most significant positive correlation with B cells. The gene-drug network included CACNA1C-Isradipine, CACNA1C-Benidipine and other relationship pairs. Totally 76 nodes and 44 edges constituted the TF-miRNA-mRNA network, including signal transducer and activator of transcription 1(STAT1) -hsa-miR-3170-AGER, CCCTC-binding factor (CTCF)-hsa-miR-455-5p-CACNA1C, etc. Moreover, qRT-PCR suggested that the expression trends of FGF1 and AGER were in keeping with the results of bioinformatic analysis. FGF1 and AGER were markedly regulated downwards in the T2DM group.
Conclusion:
Through bioinformatic analysis, we identified NETs-related diagnostic genes (ITIH3, FGF1, NRCAM, AGER, CACNA1C) in T2DM, and explored their mechanism of action from different aspects, providing new ideas for the studies related to diagnosis and treatment of T2DM.
1 Introduction
With the progress of society, the incidence and prevalence of type 2 diabetes mellitus (T2DM) has sharply risen due to the change of lifestyle, environmental pollution, mental stress and other factors, and has become a serious global health problem. T2DM is caused by a combination of genetic and environmental factors that affect the metabolism of sugar, fat and protein. Extended lasting hyperglycemia can cause multi-system harm to organs and tissues, such as, eyes (), kidneys (), nerves (Yang et al., 2022), heart and blood vessels (). The pathogenesis of T2DM is complex. There is research indicating a certain correlation between T2DM and Alzheimer’s disease (Surguchov, 2020; ). We know that there are many tissues and organs, such as pancreas, liver, muscle, kidney, heart, brain and adipose tissue, are involved in regulating blood glucose metabolism. Although there are many clinical hypoglycemic drugs that can bring many benefits to patients, they still cannot prevent the emergence of serious complications of T2DM, such as diabetic foot, diabetic nephropathy, coronary heart disease, etc., which not only cause patients great physical and psychological pain, but also a heavy financial burden. Many studies have shown that T2DM is also considered an autoimmune disease (; ; Zhou et al., 2018). Including congenital immune response deficits (neutrophils and macrophages dysfunction) and adaptive immune response dysfunction (T cells) (). Chronic hyperglycemia increases levels of various chemokines like monocyte chemokine protein-1 (MCP-1) and pro-inflammatory cytokines such as Tumour Necrosis Factor alpha (TNF-α), interleukin-1beta (IL-1β), and interleukin-6(IL-6) (). These inflammatory factors can infiltrate macrophages in the pancreas and adipose tissue, leading to insulin resistance and islet cell dysfunction ().
Neutrophil extracellular traps (NETs) are made up of filaments of chromatin DNA coated with granular proteins that can be released by neutrophils to trap microorganisms (; Yang et al., 2020). NETs are part of the innate immune response, which is beneficial to the human body, but in some cases, the disorder of NETs can also cause pathological effects on the human body (Yang et al., 2020; ). The formation of NETs is due to the induction of NETosis (inflammatory cell death mode of neutrophils) under the stimulation of various pathogens, activated platelets, chemokines, phorbol esters, etc., directly activating the protein kinase C (PKC) and Raf-MEK-ERK-MAP kinase pathways (; ). The activation of MAP kinase will initiate the formation of NADPH oxidase complex, rapidly generating reactive oxygen species (ROS). Meanwhile, key proteins in NETosis, neutrophil elastase (NE) and myeloperoxidase (MPO) (; ), contribute to nuclear membrane permeability and further development of chromatin; Peptidylarginine deiminase 4 (PAD4) modifies histones by converting arginine to citrulline, leading to chromatin depolymerization (; ). Chromatin is released outside the cell through membrane pores, ultimately releasing DNA, citrulline histone (citH3), and other intracellular particles to form a cell capture network. Owing to the unspecific action of the emitted enzyme proteins, NETs may bring out an unregulated inflammatory response leading to pathological changes resulting in direct cell damage. NETs also recruit extra pro-inflammatory cytokines, inducing the production of autoantibodies, forming immune complexes, and causing tissue damage (; ). NETs have been demonstrated to join in the pathophysiological processes of various diseases, for instance autoimmune diseases, cancer, chronic inflammation, delayed wound healing (; ), etc. In recent years, relevant literature has reported that NETs are related to the occurrence and development of diabetic nephropathy, diabetic foot ulcers and diabetes-related atherosclerotic diseases in diabetic cases (Xu et al., 2019; ).However, the definite regulatory mechanism of NETs-related genes (NETs-RGs) in T2DM remains indefinite.
In this study the relevant data sets of T2DM, including GSE21321 and GSE15932, was obtained from the Gene Expression Omnibus (GEO) database. NETs-RGs were drawn out from the literature (Wu et al., 2022). Then the diagnostic genes were screened through bioinformatic analyses such as differential expression analysis and weighted gene co-expression network analysis (WGCNA) and least absolute shrinkage and selection operator (LASSO) analysis. Moreover, the gene set enrichment analysis (GSEA), immune profile analysis, construction of the gene-drug network and transcription factor (TF)-miRNA-mRNA network were also performed on diagnostic genes to explore their mechanism of action in T2DM.This research is of great significance for the diagnosis and treatment of T2DM.
2 Materials and methods
2.1 Sources of data
T2DM-related datasets, microarray and RNA-seq of GSE21321 and GSE15932 were acquired from the GEO database (https://www.ncbi.nlm.nih.gov/geo/). Among them, GSE21321 contained blood samples from 9 T2DM and 8 control, and GSE15932 contained peripheral blood samples from 8 T2DM and 8 control. The platform of GSE21321 was GPL6883 Illumina HumanRef-8 v3.0 expression beadchip, and GSE21321 was applied as training set. The platform of GSE15932 was GPL570 [HG-U133 Plus 2] Affymetrix Human Genome U133 Plus 2.0 Array, and GSE15932 was employed as validation set to verify the expression of the diagnostic genes and the diagnostic performance for the disease (Table 1). The 137 NETs-RGs were extracted from the literature (Supplementary Table S1) (Wu et al., 2022).
TABLE 1
| Dataset | T2DM | Control | Platform |
|---|---|---|---|
| GSE21321 | 9 | 8 | GPL6883 |
| GSE15932 | 8 | 8 | GPL570 |
Information on the datasets used in this study.
2.2 Differential expression analysis and WGCNA
In the GSE21321 dataset, the differentially expressed genes (DEGs) between T2DM and control samples were sifted out by limma package (version 3.46.0) (Wang et al., 2021) setting p < 0.05 and |log2FC| > 0.5. Next, NETs scores were calculated for all samples in GSE21321 dataset utilizing the GSVA package (version 1.38.2) (), and the difference of NETs score between T2DM and control samples was compared by rank sum test. Then, to screen for key module genes associated with NETs scores, we constructed a co-expression network using WGCNA (). WGCNA serves as an algorithmic approach tailored for scrutinizing gene expression patterns across a multitude of samples. It possesses the ability to cluster genes and erect modules predicated on akin gene expression patterns, all the while scrutinizing the correlations between said modules and biological traits. Firstly, the samples in the training set were clustered and outliers were eliminated, then this network was constructed based on a soft threshold of R2 over 0.85 and connectivity close to 0. The adjacency matrix was then converted into the topological overlap matrix (TOM). Subsequently, the similar genes were classified into different modules based on dynamic tree cutting, and different colors were employed to express different modules. The NETs score was treated as the trait, Preason was then utilized to explore the correlation between these modules and the traits (p < 0.05), and the most remarkably correlated module was selected as the key module, and then the genes therein were filtered according to |gene significance (GS)| > 0.6 and |module membership (MM)| > 0.8) as key module genes.
2.3 Functional enrichment analysis of intersected genes and framework of a PPI network
We identified the intersection genes by intersecting differentially expressed genes (DEGs) with key module genes. Subsequently, the intersection genes were enrichment analysed using clusterProfiler (version 4.4.4) (Yu et al., 2012) and org.Hs.eg.db (version 3.12.0) () packages to explore the pathways and functions they were related to, encompassing gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) (adj.p < 0.05). To further understand the interaction relationships among intersection genes, PPI network (Confidence = 0.4, the discrete proteins were removed) of intersecting genes was constructed in the STRING database (https://string-db.org) and set degree cutoff = 2, node score cutoff = 0.2, k-core = 0.6, max. depth = 100 to screen candidate genes.
2.4 Acquisition of diagnostic genes
First, in the GSE21321 dataset, LASSO analysis was implemented on candidate genes to screen out diagnostic genes, and we compared the expressions of diagnostic genes between T2DM and control samples in training set and validation set. Then, receiver operating characteristic (ROC) curves were plotted utilizing the pROC package (version 1.17.0.1) (), and area under the curve (AUC) values were calculated to assess the diagnostic ability of diagnostic genes for T2DM. Further, the nomogram containing diagnostic genes was constructed using rms package (version 6.2–0) () to predict risk rates for T2DM, and its predictive value was assessed by ROC curve. Similarly, the diagnostic performance of diagnostic genes for T2DM was verified in the validation set. To probe the relevant pathways involved in diagnostic genes, we conducted gene set enrichment analysis (GSEA) on them by clusterProfiler package (version 3.18.0) (Yu et al., 2012) based on the GO-Biological process (BP) database.
2.5 Immune characteristic analysis
The study calculated the 10 immune cell proportions (scores) for all samples utilizing the quantiseq algorithm of immunedeconv (version 2.0.4). After that, the immune cell proportions were compared between T2DM and control samples by rank sum test. Eventually, spearman correlation analysis between among immune cells, and diagnostic genes and immune cells were carried out.
2.6 Framework of gene-drug and TF-miRNA-mRNA networks
In order to identify small-molecule compounds with potential therapeutic effects on T2DM, we predicted the target drugs associated with diagnostic genes via DrugBank database (https://go.drugbank.com/), and constructed gene-drug network by Cytoscape software. Then, TFs and target miRNAs of diagnostic genes were predicted via NetworkAnalyst (https://www.networkanalyst.ca/) and miRWalk (http://mirwalk.umm.uni-heidelberg.de/) databases, respectively, the TF-mRNA and miRNA-mRNA relationship pairs were merged to yield the TF-miRNA-mRNA regulatory network.
2.7 Quantitative real-time polymerase chain reaction (qRT-PCR)
We collected fresh blood samples of T2DM (n = 10) and control (n = 10) groups from the central hospital of Taiyuan city. All participants have signed informed consent forms and have been agreed by the hospital ethics committee. The characteristic information of these patients is in Supplementary Table S2. We extracted 500 μL of frozen blood from the samples and added TRIzol reagent to extract total RNA. Afterwards, 1 μL of RNA was taken and the concentration was measured using a nanophotometer N50. Reverse transcription of mRNA was performed using Servicebio’s Surscript First Strand cDNA Synthesis Kit. Dilute the reverse transcription product cDNA with ddH2O without RNase/DNase and perform qPCR reaction by diluting it to 5–20 times. GAPDH serves as an internal reference. Eventually, compare the expression of diagnostic genes in T2DM and normal blood. The primer sequence is shown in Supplementary Table S3.
2.8 Statistical analysis
In this study all bioinformatics analyses were conducted using the R package. Differential expression analyses between T2DM and controls were executed via the limma package. Co-expression network construction utilised the WGCNA package. Enrichment analyses were performed in clusterProfiler and org.Hs.eg.db packages. ROC was plotted through pROC package. Nomograms were constructed utilising the rms package. Immunological analysis was implemented in immunedeconv. Differences between groups were compared using wilcoxon test. If not otherwise specified, p < 0.05 indicates statistical significance.
3 Results
3.1 Identification of DEGs and key module
There were 671 DEGs between T2DM and control samples (Figure 1A), and the heat map illustrated the top50 of them (Figure 1B). As seen from the violin plot, the NETs score was remarkably higher in T2DM samples than that in the control samples (Supplementary Figure S1). In Figure 1C, it's evident that the clustering performance of the dataset samples was remarkably strong, indicating no necessity for sample exclusion. A soft threshold of 11 was applied to construct the scale-free network when R2 was approximately 0.85 and the connectivity was close to 0 (Figure 1D). Subsequently, genes were grouped into 11 modules (Figure 1E). Among which, MEturquoise (Cor = −0.63 and p < 0.05) had the strongest correlation with NETs score, so we treated it as the key module (Figure 1F). MEturquoise module contained 3,368 key module genes (|GS| > 0.6 and |MM| > 0.8), and they were positively correlated with the NETs score (Figure 1G).
FIGURE 1
3.2 Functional enrichment and PPI analyses of intersection genes
DEGs and key module genes were intersected to yield 302 intersection genes (Figure 2A). GO items displayed that these genes were enriched to cellular response to amyloid−beta, growth factor activity, transporter complex and so on (Figure 2B). KEGG pathways revealed that intersection genes were primarily related to lipid and atherosclerosis, aldosterone synthesis and secretion, Interleukin 17(IL−17) signaling pathway, etc. (Figure 2C). A total of 165 nodes and 203 edges constituted the PPI network, containing Interleukin 17 Receptor E (IL17RE)-FAM92B, Brain-derived neurotrophic factor (BDNF)- calcium voltage-gated channel subunit alpha1 C (CACNA1C), peroxisome proliferator-activated receptor gamma (PPARG)- Paired box protein 6(PAX6) and other reciprocal relationship pairs (Figure 2D). After analyzing the topological properties of the PPI network, we acquired 23 candidate genes, such as PPARG, acetyl-CoA synthetase 2 (ACSS2), Complement factor I (CFI) (Figure 2E; Supplementary Table S4).
FIGURE 2
3.3 Excellent diagnostic power of diagnostic genes for T2DM
The five diagnostic genes, namely, inter-trypsin inhibitor heavy chain 3 (ITIH3), fibroblast growth factor 1 (FGF1), neuron cell adhesion molecule (NRCAM), advanced glycosylation end-product-specific receptor (AGER), and CACNA1C, were identified when the lambda min = 0.0106 via LASSO analysis (Figure 3A), and they were all downregulated in T2DM samples of the GSE21321 dataset (Supplementary Figure S2). Further,AUC values of ITIH3 (AUC = 0.833), FGF1 (AUC = 0.792), NRCAM (AUC = 0.847), AGER (AUC = 0.764), and CACNA1C (AUC = 0.750) were all greater than 0.7, indicating that diagnostic genes had strong diagnostic power for T2DM (Figure 3B). The nomogram suggested that the diagnostic genes had strong risk prediction ability for T2DM, and the results were further validated by the ROC curve (AUC > 0.7) (Figures 3C,D). Similarly, the diagnostic efficacy of diagnostic genes for T2DM was validated in the validation set (Supplementary Figures S3, S4).
FIGURE 3
3.4 Diagnostic genes were involved in multiple pathways
We applied GSEA to further explore the biological pathways that diagnostic genes were involved in T2DM, and found that all five diagnostic genes were mainly associated with detection of chemical stimulus, detection of chemical stimulus involved in sensory perception, detection of chemical stimulus involved in sensory perception of smell, detection of stimulus involved in sensory perception, pattern specification process, sensory perception of chemical stimulus and sensory perception of smell (Figure 4). This suggests that these genes may influence the onset and development of T2DM through multiple pathways.
FIGURE 4
3.5 Immuno-infiltration analysis
From the box plot, it could be concluded that the proportion of Neutrophil and Tregs in the samples were relatively high (Figure 5A). M2 Macrophage, Monocyte, Natural killer (NK) cell, and Myeloid dendritic cells (DC) were dramatically different between T2DM and control samples (Figure 5B). The all diagnostic genes had positive and the most significant correlation with B cells (Figure 5C). Correlation analysis showed differences in correlations among 9 immune cells, among them, Monocyte was markedly positively correlated with Myeloid DC, while it was remarkably negatively correlated with Regulatory T cells (Tregs) (Figure 5D).
FIGURE 5
3.6 Construction of gene-drug and TF-miRNA-mRNA networks
The gene-drug network included CACNA1C-Isradipine, CACNA1C-Benidipine, AGER-Pyridoxamine, ITIH3-Clozapine, FGF1-Muparfostat and other relationship pairs (Figure 6A). Figure 6B demonstrated the 2D structure of some drugs. A total of 76 nodes and 44 edges constituted the TF-miRNA-mRNA network, including signal transducer and activator of transcription 1(STAT1)-hsa-miR-3170-AGER, CCCTC-binding factor (CTCF)-hsa-miR-455-5p-CACNA1C, GA binding protein transcription factor Alpha (GABPA) -hsa-miR-452-5p-NRCAM, TEA domain transcription factor 3 (TEAD3)-hsa-miR-4749-3p-ITIH3, ZFP2 zinc finger protein (ZFP2)-hsa-miR-4265-AGER et al. (Figure 6C).
FIGURE 6
3.7 Validation of diagnostic genes expressions by qRT-PCR
In comparison with the control group, FGF1 and AGER were markedly regulated downwards in the T2DM group (Figure 7). The expression trends of these two diagnostic genes were consistent with the consequences of bioinformatic analysis.
FIGURE 7
4 Discussion
The pathogenesis of T2DM is multifaceted, primarily characterized by insulin resistance and relative insulin insufficiency. The greatest challenge in our treatment of T2DM remains the significant harm caused by its complications. Immune dysfunction is increasingly recognized as a fundamental pathophysiological mechanism underlying both T2DM and its chronic complications (). Studies have displayed that NETs have a crucial role in the initiation and development of diabetes-associated atherosclerosis (AS) (), diabetic nephropathy (Zheng et al., 2022), diabetic retinopathy () and diabetic foot ulcer (Yang et al., 2023). However, the pathogenesis of NETs and T2DM is still not clear. Through bioinformatics analysis, this study successfully identified and characterized NETs-RGs in T2DM. The diagnostic genes ITIH3, FGF1, NRCAM, AGER, and CACNA1C showed significant association with NETs in T2DM, underscoring their potential as valuable biomarkers for diagnosing T2DM. Furthermore, studies on their mechanisms using immune infiltration analysis, enrichment analysis and regulatory network analysis have provided important insights into the role of NETs-RGs in T2DM. These findings present innovative ideas for the diagnosis and treatment of T2DM patients.
FGF1 is involved in the regulation of glycolipid metabolism. Studies have shown that FGF1 reduces blood sugar by increasing insulin sensitivity, providing a new approach for the treatment of metabolic diseases like fatness, NAFLD and T2DM (; ). It plays a role in coronary heart disease, myocardial ischemia, nerve injury and wound healing (). FGF1 can enhance the healing process of diabetic wounds (). In individuals with diabetes, during the late stages of inflammation in injured tissue, macrophages often remain in an inflammatory state, failing to transition into a reparative phenotype. Consequently, they are unable to secrete factors that facilitate tissue repair, hindering the transition of the wound into the proliferation phase and leading to the development of chronic inflammation (). Our study showed that in contrast with the control group, the proportion of neutrophils in T2DM was remarkably increased, and the proportion of M2 macrophages was significantly reduced, thus promoting the occurrence and development of inflammation. It was also found that excessive neutrophils release NETs in diabetes tissue damage. NETs can activate NLRP 3 inflammasomes in macrophages, thereby releasing IL-1β. These inflammatory cytokines exist in diabetes wounds for much longer than normal, which may prolong the inflammatory phase and inhibit the formation of granulation tissue (). In this study, FGF1 was significantly downregulated in both bioinformatics analysis and qRT-PCR validation in T2DM patients, which may exacerbate abnormal glucose and lipid metabolism and promote the occurrence and development of diabetic complications and delayed wound healing. Therefore, FGF1 may be a latent new aim for the treatment of insulin resistance and T2DM in the future.
AGER is a highly polymorphic gene involved in multiple pathways, such as nuclear factor (NF-κB), protein kinase B (AKT), p38 and mitogen-activated protein (MAP) kinase, which are activated to cause proinflammatory states in the body and may be linked to the growth and development of human autoimmune diseases (), diabetes complications, cancer (), coronary heart disease and lung disease. In this study, we found that intersection genes are trapped in IL-17 signaling pathway, aldosterone synthesis and secretion, cyclic adenosine monophosphate (cAMP), rat sarcoma (RAS), phosphoinositide 3-kinase (PI3K) Akt and mitogen-activated protein kinases (MAPK) signaling pathway, advanced glycation end product (AGE)-receptor for AGE (RAGE) signaling pathway in diabetes complications, which is consistent with literature results. In patients with T2DM, prolonged hyperglycemia leads to the accumulation of advanced glycation end products (AGEs). AGER serves as a receptor for AGEs, initiating AGE-AGER interaction, which in turn activates nuclear factor-kappaB (NF-κB). NF-κB activation triggers the transcription of various pro-inflammatory cytokines such as IL-1α, IL-6, and TNF-α, as well as growth factors and adhesion molecules. This cascade promotes the production of ROS, membrane oxidation, increased membrane permeability, and activation of phospholipase A2 (PLA2). These events collectively contribute to the onset and progression of vascular complications in diabetes (Yuan et al., 2019; ). Research shows that (), there are significantly more macrophages in diabetes plaques, and the expression of AGER is increased in monocytes/macrophages, the key cells in atherosclerosis, so diabetes plaques show more expression of AGER. This study demonstrated that the expression of AGER in T2DM was lower than that in the control group, which may be associated with the decrease of AGER expression due to the excessive accumulation of AGEs caused by long-term hyperglycemia. Therefore, AGER may get trapped in the occurrence and development of T2DM, but its specific mechanism needs further exploration and research.
CACNA1C encodes the L-type voltage-dependent calcium channel. Insulin secretion is regulated by two distinct phases: early-phase secretion and late-phase secretion, each governed by different regulatory mechanisms. Early-phase insulin secretion is initiated by the uptake of glucose via glucose transporter (GLUT) transporters. Subsequently, the activation of potassium channels leads to the influx of Ca2+ through calcium channels, facilitating the release of insulin stored in vesicles through exocytosis (). Glucose-stimulated insulin secretion (GSIS) of pancreatic beta cells plays a key role in maintaining glucose homeostasis and preventing hyperglycemia. It has been shown that Ca2+ flows primarily into beta cells via L-type voltage-dependent Ca2+ (Ca V) channels, which are necessary for GSIS (). This study revealed a significant decrease in the expression of CACNA1C in T2DM patients compared to the control group. This decrease in CACNA1C expression could potentially impair early insulin secretion and GSIS, suggesting it may play a crucial role in the pathogenesis of T2DM. Interestingly, inhibition of excess Ca2+ inflow in pre-diabetic islets by the drug isradipine to prevent Ca2+ induced hypersensitivity may be an option for early intervention, at least to delay the onset of T2DM (; ).
In this study, M2 Macrophage, Monocyte, NK cell, Myeloid DC were found to be significantly different between T2DM and control samples. The positive correlation between diagnostic genes and B cells was the strongest and most significant. T2DM is recognized as a chronic low-grade inflammatory condition, driven by the inflammatory activation of both recruited and resident macrophages, known as the classical pathway. This entails the infiltration and activation of macrophages within adipose tissue, which is closely linked to insulin resistance induced by obesity (). Most metabolic organs are made up of dense cells of the inborn and adaptive immune systems. Macrophages play an important role in atherosclerosis (). M2 Macrophages has anti-inflammatory and immunomodulatory effects, which helps to promote tissue repair, remodeling, healing, etc. (; ; ). This study demonstrated a notable decrease in M2 Macrophages in T2DM patients in comparison with control group, which may be an important reason for delayed wound healing in T2DM. Research has shown that hyperglycemia accelerates the formation of AGEs. Compared with non-diabetes plaques, diabetes plaques have significantly more macrophages. At the same time, diabetes plaques show more expression of AGER (especially in Macrophages). The interaction between AGEs and AGER promotes the occurrence of inflammation and OS, leading to the formation of foam cells in macrophages, which leads to the occurrence and development of diabetes related atherosclerosis (). Therefore, studying the relationship between T2DM and immune cells is particularly important. NK cells are believed to be associated with T2DM by regulating systemic inflammation (; Zitti and Bryceson, 2018). The study reported that the NK cells contribute to inflammation and insulin resistance induced by obesity (). The NK cell activity of T2DM patients is significantly reduced and linearly deteriorates with the degree of hyperglycemia, which in turn disrupts the maintenance of the immune system (). This study demonstrated a notable decrease in NK cell levels among T2DM patients, with consistent results across the board. Consequently, comprehending the functional alterations of NK cells presents novel avenues for addressing T2DM resulting from obesity.
This study was conducted through DrugBank (https://go.drugbank.com/) online databases predict targeted drugs related to diagnostic genes. In our study, we found that Pyridoxamine was one of the drugs searched for by AGER. Pyridoxamine is a derivative of vitamin B6, an inhibitor of advanced glycation and lipid oxidation end products. Studies have shown that pyridoxamine can inhibit the kidney lesions caused by AGEs, and can also significantly reduce plasma cholesterol, triglyceride, creatinine and urinary albumin in diabetic rats, and reduce the content of skin collagen AGEs (). Muparfostat is also known as PI-88, an inhibitor of heparinase (HPSE). HPSE is an endo-β-D-glucuronidase capable of degrading heparan sulfate (HS). Studies have reported that the increase of blood glucose can lead to the upregulation of HPSE expression, and with the prolongation of diabetes, HPSE may cause more serious damage to islet function (; ). At the same time, HS acts in the FGFs signaling pathway. FGFRS mediates FGFs signal transduction through heparan sulfate or Klotho-dependent ways. Long-term hyperglycemia leads to the overexpression of HPSE and the destruction of HS, which in turn affects the FGFs signaling pathway and has a negative impact on blood glucose control (; ). In the future, specific studies can be conducted on the treatment of T2DM with Pyridoxamine and Muparfostat.
In this study, the relationship between T2DM and NETs-RGs was discussed through bioinformatics analysis, and five diagnostic genes related to NETs were identified in T2DM, including ITIH3, FGF1, NRCAM, AGER, and CACNA1C. The mechanism of action of diagnostic genes on T2DM was also explored from different perspectives. It provides a new way to study the diagnosis and treatment of T2DM. However, the method of confidence generation limited our ability to predict the TFs and related regulatory networks of only five diagnostic genes using the online database. The association between certain diagnostic genes, such as ITIH3 and NRCAM, and diabetes remains understudied. Only FGF1 and AGER genes were validated by qRT PCR in this study, while the other three genes were not validated. This may be related to small sample sizes, limited to specific populations, or possible biases during the validation process. As a result, the other three were not validated successfully. The precise regulatory mechanisms of NETs-RGs AGER, FGF1, and CACNA1C in T2DM have not been fully elucidated. Moving forward, we will continue to closely monitor the research progress concerning these diagnostic genes and their roles in T2DM.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The studies involving humans were approved by the Ethics Committee of Taiyuan Central Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author contributions
MH: Writing–review and editing, Writing–original draft. JN: Writing–review and editing, Writing–original draft. HC: Funding acquisition, Writing–review and editing. CG: Writing–review and editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2024.1373807/full#supplementary-material
References
1
AmorimM.MartinsB.FernandesR. (2023). Immune fingerprint in diabetes: ocular surface and retinal inflammation. Int. J. Mol. Sci.24, 9821. 10.3390/ijms24129821
2
ArabpourM.SaghazadehA.RezaeiN. (2021). Anti-inflammatory and M2 macrophage polarization-promoting effect of mesenchymal stem cell-derived exosomes. Int. Immunopharmacol.97, 107823. 10.1016/j.intimp.2021.107823
3
ArakiK.ArakiA.HondaD.IzumotoT.HashizumeA.HijikataY.et al (2019). TDP-43 regulates early-phase insulin secretion via CaV1.2-mediated exocytosis in islets. J. Clin. Investig.129, 3578–3593. 10.1172/jci124481
4
BerbudiA.RahmadikaN.TjahjadiA. I.RuslamiR. (2020). Type 2 diabetes and its impact on the immune system. Curr. Diabetes Rev.16, 442–449. 10.2174/1573399815666191024085838
5
BoniakowskiA. E.KimballA. S.JoshiA.SchallerM.DavisF. M.denDekkerA.et al (2018). Murine macrophage chemokine receptor CCR2 plays a crucial role in macrophage recruitment and regulated inflammation in wound healing. Eur. J. Immunol.48, 1445–1455. 10.1002/eji.201747400
6
BorissoffJ. I.JoosenI. A.VersteylenM. O.BrillA.FuchsT. A.SavchenkoA. S.et al (2013). Elevated levels of circulating DNA and chromatin are independently associated with severe coronary atherosclerosis and a prothrombotic state. Arterioscler. Thromb. Vasc. Biol.33, 2032–2040. 10.1161/atvbaha.113.301627
7
BrinkmannV.ZychlinskyA. (2007). Beneficial suicide: why neutrophils die to make NETs. Nat. Rev. Microbiol.5, 577–582. 10.1038/nrmicro1710
8
BruschiM.BonanniA.PetrettoA.VaglioA.PratesiF.SantucciL.et al (2020). Neutrophil extracellular traps profiles in patients with incident systemic lupus erythematosus and lupus nephritis. J. Rheumatol.47, 377–386. 10.3899/jrheum.181232
9
CaligiuriM. A. (2008). Human natural killer cells. Blood112, 461–469. 10.1182/blood-2007-09-077438
10
CaveolinS. A. (2020). Caveolin: a new link between diabetes and ad. Cell. Mol. Neurobiol.40 (7), 1059–1066. Epub 2020 Jan 23. PMID: 31974905. 10.1007/s10571-020-00796-4
11
CoopeA.TorsoniA. S.VellosoL. A. (2016). Mechanisms in endocrinology: metabolic and inflammatory pathways on the pathogenesis of type 2 diabetes. Eur. J. Endocrinol.174, R175–R187. 10.1530/eje-15-1065
12
DalmasE. (2019). Role of innate immune cells in metabolism: from physiology to type 2 diabetes. Semin. Immunopathol.41, 531–545. 10.1007/s00281-019-00736-5
13
DhounchakS.PoppS. K.BrownD. J.LaybuttD. R.BidenT. J.BornsteinS. R.et al (2021). Heparan sulfate proteoglycans in beta cells provide a critical link between endoplasmic reticulum stress, oxidative stress and type 2 diabetes. PLoS One16, e0252607. 10.1371/journal.pone.0252607
14
DickersonM. T.DadiP. K.ButterworthR. B.NakheA. Y.GraffS. M.ZaborskaK. E.et al (2020). Tetraspanin-7 regulation of L-type voltage-dependent calcium channels controls pancreatic β-cell insulin secretion. J. Physiol.598, 4887–4905. 10.1113/jp279941
15
FerroV.HammondE.FairweatherJ. K. (2004). The development of inhibitors of heparanase, a key enzyme involved in tumour metastasis, angiogenesis and inflammation. Mini Rev. Med. Chem.4, 693–702. 10.2174/1389557043403729
16
HawezA.Al-HaidariA.MadhiR.RahmanM.ThorlaciusH. (2019). MiR-155 regulates PAD4-dependent formation of neutrophil extracellular traps. Front. Immunol.10, 2462. 10.3389/fimmu.2019.02462
17
HirotaT.LevyJ. H.IbaT. (2020). The influence of hyperglycemia on neutrophil extracellular trap formation and endothelial glycocalyx damage in a mouse model of type 2 diabetes. Microcirculation27, e12617. 10.1111/micc.12617
18
ItariuB. K.StulnigT. M. (2014). Autoimmune aspects of type 2 diabetes mellitus - a mini-review. Gerontology60, 189–196. 10.1159/000356747
19
JoS.FangS. (2021). Therapeutic strategies for diabetes: immune modulation in pancreatic β cells. Front. Endocrinol. (Lausanne)12, 716692. 10.3389/fendo.2021.716692
20
JonkerJ. W.SuhJ. M.AtkinsA. R.AhmadianM.LiP.WhyteJ.et al (2012). A PPARγ-FGF1 axis is required for adaptive adipose remodelling and metabolic homeostasis. Nature485, 391–394. 10.1038/nature10998
21
KhalifahR. G.BaynesJ. W.HudsonB. G. (1999). Amadorins: novel post-Amadori inhibitors of advanced glycation reactions. Biochem. Biophys. Res. Commun.257, 251–258. 10.1006/bbrc.1999.0371
22
LangfelderP.HorvathS. (2008). WGCNA: an R package for weighted correlation network analysis. BMC Bioinforma.9, 559. 10.1186/1471-2105-9-559
23
LewisH. D.LiddleJ.CooteJ. E.AtkinsonS. J.BarkerM. D.BaxB. D.et al (2015). Inhibition of PAD4 activity is sufficient to disrupt mouse and human NET formation. Nat. Chem. Biol.11, 189–191. 10.1038/nchembio.1735
24
LiH. D.YouY. K.ShaoB. Y.WuW. F.WangY. F.GuoJ. B.et al (2022). Roles and crosstalks of macrophages in diabetic nephropathy. Front. Immunol.13, 1015142. 10.3389/fimmu.2022.1015142
25
LiX. (2019). The FGF metabolic axis. Front. Med.13, 511–530. 10.1007/s11684-019-0711-y
26
LiuD.YangP.GaoM.YuT.ShiY.ZhangM.et al (2019). NLRP3 activation induced by neutrophil extracellular traps sustains inflammatory response in the diabetic wound. Clin. Sci. (Lond).133, 565–582. 10.1042/cs20180600
27
LiuT. T.LiR.HuoC.LiJ. P.YaoJ.JiX. L.et al (2021a). Identification of CDK2-related immune forecast model and ceRNA in lung adenocarcinoma, a pan-cancer analysis. Front. Cell. Dev. Biol.9, 682002. 10.3389/fcell.2021.682002
28
LiuY.LiuY.DengJ.LiW.NieX. (2021b). Fibroblast growth factor in diabetic foot ulcer: progress and therapeutic prospects. Front. Endocrinol. (Lausanne)12, 744868. 10.3389/fendo.2021.744868
29
LoweW. L.Jr.BoydF. T.ClarkeD. W.RaizadaM. K.HartC.LeRoithD. (1986). Development of brain insulin receptors: structural and functional studies of insulin receptors from whole brain and primary cell cultures. Endocrinology119, 25–35. 10.1210/endo-119-1-25
30
Magaña-GuerreroF. S.Aguayo-FloresJ. E.Buentello-VolanteB.Zarco-ÁvilaK.Sánchez-CisnerosP.Castro-SalasI.et al (2023). Spontaneous neutrophil extracellular traps release are inflammatory markers associated with hyperglycemia and renal failure on diabetic retinopathy. Biomedicines11, 1791. 10.3390/biomedicines11071791
31
MalikP.ChaudhryN.MittalR.MukherjeeT. K. (2015). Role of receptor for advanced glycation end products in the complication and progression of various types of cancers. Biochim. Biophys. Acta1850, 1898–1904. 10.1016/j.bbagen.2015.05.020
32
MarassiM.FadiniG. P. (2023). The cardio-renal-metabolic connection: a review of the evidence. Cardiovasc Diabetol.22 (1), 195. PMID: 37525273; PMCID: PMC10391899. 10.1186/s12933-023-01937-x
33
MasolaV.GambaroG.TibaldiE.OnistoM.AbaterussoC.LupoA. (2011). Regulation of heparanase by albumin and advanced glycation end products in proximal tubular cells. Biochim. Biophys. Acta1813, 1475–1482. 10.1016/j.bbamcr.2011.05.004
34
NienhuisH. L.WestraJ.SmitA. J.LimburgP. C.KallenbergC. G.BijlM. (2009). AGE and their receptor RAGE in systemic autoimmune diseases: an inflammation propagating factor contributing to accelerated atherosclerosis. Autoimmunity42, 302–304. 10.1080/08916930902831746
35
PapayannopoulosV.MetzlerK. D.HakkimA.ZychlinskyA. (2010). Neutrophil elastase and myeloperoxidase regulate the formation of neutrophil extracellular traps. J. Cell. Biol.191, 677–691. 10.1083/jcb.201006052
36
ParkerH.WinterbournC. C. (2012). Reactive oxidants and myeloperoxidase and their involvement in neutrophil extracellular traps. Front. Immunol.3, 424. 10.3389/fimmu.2012.00424
37
PetrettoA.BruschiM.PratesiF.CroiaC.CandianoG.GhiggeriG.et al (2019). Neutrophil extracellular traps (NET) induced by different stimuli: a comparative proteomic analysis. PLoS One14, e0218946. 10.1371/journal.pone.0218946
38
PĭtreM.GaudreaultN.SanturéM.NadeauA.BachelardH. (1999). Isradipine and insulin sensitivity in hypertensive rats. Am. J. Physiol.276, E1038–E1048. 10.1152/ajpendo.1999.276.6.E1038
39
PoznyakA.GrechkoA. V.PoggioP.MyasoedovaV. A.AlfieriV.OrekhovA. N. (2020). The diabetes mellitus-atherosclerosis connection: the role of lipid and glucose metabolism and chronic inflammation. Int. J. Mol. Sci.21, 1835. 10.3390/ijms21051835
40
Rayego-MateosS.Marquez-ExpósitoL.Rodrigues-DiezR.SanzA. B.GuiterasR.DoladéN.et al (2022). Molecular mechanisms of kidney injury and repair. Int. J. Mol. Sci.23, 1542. 10.3390/ijms23031542
41
ReevesE. P.LuH.JacobsH. L.MessinaC. G.BolsoverS.GabellaG.et al (2002). Killing activity of neutrophils is mediated through activation of proteases by K+ flux. Nature416, 291–297. 10.1038/416291a
42
ReinbotheT. M.AlkayyaliS.AhlqvistE.TuomiT.IsomaaB.LyssenkoV.et al (2013). The human L-type calcium channel Cav1.3 regulates insulin release and polymorphisms in CACNA1D associate with type 2 diabetes. Diabetologia56, 340–349. 10.1007/s00125-012-2758-z
43
RobinX.TurckN.HainardA.TibertiN.LisacekF.SanchezJ. C.et al (2011). pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinforma.12, 77. 10.1186/1471-2105-12-77
44
RowishaM.El-BatchM.El ShikhT.El MelegyS.AlyH. (2016). Soluble receptor and gene polymorphism for AGE: relationship with obesity and cardiovascular risks. Pediatr. Res.80, 67–71. 10.1038/pr.2016.55
45
SchmidtA. M. (2017). 22016 ATVB plenary lecture: receptor for advanced glycation endproducts and implications for the pathogenesis and treatment of cardiometabolic disorders: spotlight on the macrophage. Arterioscler. Thromb. Vasc. Biol.37, 613–621. 10.1161/atvbaha.117.307263
46
SiY.LiuJ.ShanW.ZhangY.HanC.WangR.et al (2020). Association of lymphocyte-to-monocyte ratio with total coronary plaque burden in patients with coronary artery disease. Coron. Artery Dis.31, 650–655. 10.1097/mca.0000000000000857
47
SimeonovicC. J.ZiolkowskiA. F.WuZ.ChoongF. J.FreemanC.ParishC. R. (2013). Heparanase and autoimmune diabetes. Front. Immunol.4, 471. 10.3389/fimmu.2013.00471
48
SmithC. K.Vivekanandan-GiriA.TangC.KnightJ. S.MathewA.PadillaR. L.et al (2014). Neutrophil extracellular trap-derived enzymes oxidize high-density lipoprotein: an additional proatherogenic mechanism in systemic lupus erythematosus. Arthritis Rheumatol.66, 2532–2544. 10.1002/art.38703
49
SuR.JinC.BuH.XiangJ.ZhouL.JinC. (2022). Development and validation of an immune-related prognostic signature in cervical cancer. Front. Oncol.12, 861392. 10.3389/fonc.2022.861392
50
SuhJ. M.JonkerJ. W.AhmadianM.GoetzR.LackeyD.OsbornO.et al (2014). Endocrinization of FGF1 produces a neomorphic and potent insulin sensitizer. Nature513, 436–439. 10.1038/nature13540
51
SurguchovA. C. (2020). A new link between diabetes and AD. Cell. Mol. Neurobiol.40 (7), 1059–1066. 10.1007/s10571-020-00796-4
52
TanT. E.WongT. Y. (2023). Diabetic retinopathy: looking forward to 2030. Front. Endocrinol. (Lausanne)13, 1077669. PMID: 36699020; PMCID: PMC9868457. 10.3389/fendo.2022.1077669
53
TheurichS.TsaousidouE.HanssenR.LempradlA. M.MauerJ.TimperK.et al (2017). IL-6/Stat3-Dependent induction of a distinct, obesity-associated NK cell subpopulation deteriorates energy and glucose homeostasis. Cell. Metab.26, 171–184. 10.1016/j.cmet.2017.05.018
54
ThipsawatS. (2021). Early detection of diabetic nephropathy in patient with type 2 diabetes mellitus: a review of the literature. Diab Vasc. Dis. Res.18 (6), 14791641211058856. 10.1177/14791641211058856
55
WangH. W.TangJ.SunL.LiZ.DengM.DaiZ. (2023). Mechanism of immune attack in the progression of obesity-related type 2 diabetes. World J. Diabetes14, 494–511. 10.4239/wjd.v14.i5.494
56
WangY.LiM.StadlerS.CorrellS.LiP.WangD.et al (2009). Histone hypercitrullination mediates chromatin decondensation and neutrophil extracellular trap formation. J. Cell. Biol.184, 205–213. 10.1083/jcb.200806072
57
WangY.WangZ.SunJ.QianY. (2021). Identification of HCC subtypes with different prognosis and metabolic patterns based on mitophagy. Front. Cell. Dev. Biol.9, 799507. 10.3389/fcell.2021.799507
58
WuJ.ZhangF.ZhengX.ZhangJ.CaoP.SunZ.et al (2022). Identification of renal ischemia reperfusion injury subtypes and predictive strategies for delayed graft function and graft survival based on neutrophil extracellular trap-related genes. Front. Immunol.13, 1047367. 10.3389/fimmu.2022.1047367
59
XuH.JiangJ.ChenW.LiW.ChenZ. (2019). Vascular macrophages in atherosclerosis. J. Immunol. Res.2019, 4354786. 10.1155/2019/4354786
60
YangL.LiuQ.ZhangX.LiuX.ZhouB.ChenJ.et al (2020). DNA of neutrophil extracellular traps promotes cancer metastasis via CCDC25. Nature583, 133–138. 10.1038/s41586-020-2394-6
61
YangS.WangS.ChenL.WangZ.ChenJ.NiQ.et al (2023). Neutrophil extracellular traps delay diabetic wound healing by inducing endothelial-to-mesenchymal transition via the hippo pathway. Int. J. Biol. Sci.19 (1), 347–361. 10.7150/ijbs.78046
62
YuG.WangL. G.HanY.HeQ. Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters. Omics16, 284–287. 10.1089/omi.2011.0118
63
YangK.WangY.LiY. W.ChenY. G.XingN.LinH. B. (2022). Progress in the treatment of diabetic peripheral neuropathy. Biomed Pharmacother. 148, 112717. 10.1016/j.biopha.2022.112717
64
YuanT.YangT.ChenH.FuD.HuY.WangJ.et al (2019). New insights into oxidative stress and inflammation during diabetes mellitus-accelerated atherosclerosis. Redox Biol.20, 247–260. 10.1016/j.redox.2018.09.025
65
ZhengF.MaL.LiX.WangZ.GaoR.PengC.et al (2022). Neutrophil extracellular traps induce glomerular endothelial cell dysfunction and pyroptosis in diabetic kidney disease. Diabetes71, 2739–2750. 10.2337/db22-0153
66
ZhouT.HuZ.YangS.SunL.YuZ.WangG. (2018). Role of adaptive and innate immunity in type 2 diabetes mellitus. J. Diabetes Res.2018, 7457269. 10.1155/2018/7457269
67
ZittiB.BrycesonY. T. (2018). Natural killer cells in inflammation and autoimmunity. Cytokine Growth Factor Rev.42, 37–46. 10.1016/j.cytogfr.2018.08.001
Summary
Keywords
neutrophil extracellular traps, type 2 diabetes mellitus, GEO, bioinformatic, diagnosis
Citation
He M, Niu J, Cheng H and Guo C (2024) Identification and validation of diagnostic genes associated with neutrophil extracellular traps of type 2 diabetes mellitus. Front. Genet. 15:1373807. doi: 10.3389/fgene.2024.1373807
Received
12 February 2024
Accepted
20 August 2024
Published
04 September 2024
Volume
15 - 2024
Edited by
Jinhui Liu, Nanjing Medical University, China
Reviewed by
Andrei Surguchov, University of Kansas Medical Center, United States
Geeta Rai, Banaras Hindu University, India
Updates
Copyright
© 2024 He, Niu, Cheng and Guo.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Meifang He, 717120015@qq.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.