Skip to main content

Identification of co-expressed central genes and transcription factors in acute myocardial infarction and diabetic nephropathy

Abstract

Background

Acute myocardial infarction (AMI) and diabetic nephropathy (DN) are common clinical co-morbidities, but they are challenging to manage and have poor prognoses. There is no research on the bioinformatics mechanisms of comorbidity, and this study aims to investigate such mechanisms.

Methods

We downloaded the AMI data (GSE66360) and DN datasets (GSE30528 and GSE30529) from the Gene Expression Omnibus (GEO) platform. The GSE66360 dataset was divided into two parts: the training set and the validation set, and GSE30529 was used as the training set and GSE30528 as the validation set. After identifying the common differentially expressed genes (DEGs) in AMI and DN in the training set, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses and protein–protein interaction (PPI) network construction were performed. A sub-network graph was constructed by MCODE, and 15 hub genes were screened by the Cytohubba plugin. The screened hub genes were validated, and the 15 screened hub genes were subjected to GO, KEGG, Gene MANIA analysis, and transcription factor (TF) prediction. Finally, we performed TF differential analysis, enrichment analysis, and TF and gene regulatory network construction.

Results

A total of 46 genes (43 up-regulated and 3 down-regulated) were identified for subsequent analysis. GO functional analysis emphasized the presence of genes mainly in the vesicle membrane and secretory granule membrane involved in antigen processing and presentation, lipopeptide binding, NAD + nucleosidase activity, and Toll-like receptor binding. The KEGG pathways analyzed were mainly in the phagosome, neutrophil extracellular trap formation, natural killer cell-mediated cytotoxicity, apoptosis, Fc gamma R-mediated phagocytosis, and Toll-like receptor signaling pathways. Eight co-expressed hub genes were identified and validated, namely TLR2, FCER1G, CD163, CTSS, CLEC4A, IGSF6, NCF2, and MS4A6A. Three transcription factors were identified and validated in AMI, namely NFKB1, HIF1A, and SPI1.

Conclusions

Our study reveals the common pathogenesis of AMI and DN. These common pathways and hub genes may provide new ideas for further mechanistic studies.

Peer Review reports

Introduction

The heart and kidneys are vital organs in the body for maintaining circulatory and internal environmental homeostasis. There is cross-talk between the functions of the two organs, and acute or chronic injury to either organ may cause adaptive changes or adverse reactions in the two organs, resulting in cardiorenal syndrome [1]. Acute myocardial infarction (AMI) is also one of the most common types of acute kidney injury. While diabetes is one of the major risk factors for atherosclerotic vascular disease, it promotes circulatory disorders in the renal and coronary arteries, often concurrently or serially, causing CHD or diabetic nephropathy (DN). Cardiovascular disease mediates nearly 50–80% of deaths in patients with diabetes [2, 3], while diabetic nephropathy is a widespread complication of diabetes and a major cause of end-stage renal disease (ESRD) [4]. Although AMI and DN are both circulatory disorders with common co-morbidities, the molecular mechanisms explaining the coexistence of these two diseases are unclear. Studying common transcriptional signatures may provide new insights into the common pathogenesis of AMI and DN.

AMI is inextricably associated with DN, and studies have shown that the 10-year risk of death is also increased by 10% for mild renal impairment compared to normal renal function in AMI patients. Severe renal impairment has more than doubled the risk of death at 1 year after AMI, and this increased risk persisted, with creatinine clearance becoming a risk factor beyond other risk factors and second only to age by 10 years [5]. Emergency percutaneous coronary intervention (PCI) is widely acknowledged as an effective treatment for AMI to save the endangered myocardium and reduce mortality and disability. However, contrast-induced nephropathy (CIN) is a common complication in patients with PCI, and hyperglycemia increases the risk of CIN after PCI [6], is associated with a lengthened hospital stay and later adverse outcomes [7, 8], and is associated with a significantly increased risk of CIN in AMI patients undergoing emergency PCI [9, 10], which is more pronounced in patients with DN. Diabetes increases the risk of macrovascular and microvascular complications, including AMI and DN. Studies have shown that the annual cases of AMI and ESRD in Singaporean diabetic patients are projected to increase by 50% by 2050 [11, 12], and it is assumed that this risk may exist in many countries around the world. It is urgent to analyze the co-morbidity mechanisms of AMI and DN from a bioinformatics perspective, which is why we conducted the study.

Materials and methods

Data source

We downloaded the expression profile dataset of endothelial cells isolated from the circulation of AMI patients (GSE66360), the transcriptome dataset of glomerular specimens from DN patients (GSE30528), and the transcriptome dataset of tubular specimens from DN patients (GSE30529) from the GEO (http://www.ncbi.nlm.nih.gov/geo/) database. GSE66360 includes endothelial cells isolated from the circulation of 49 AMI and 50 healthy controls; we selected the first 13 controls and 21 AMI samples as the training set and the last 37 controls and 28 AMI samples as the validation set. GSE30528 includes transcriptomic data from glomerular specimens of 13 controls and 9 DNs as the validation set. GSE30529 includes 12 controls and 10 DNs of renal tubular transcriptome data as the training set. The GSE30528 and GSE30529 datasets were matched for gender, and while the difference in age was significant in GSE30528 (p = 3.2E-02), the difference in gender was not in GSE30529 (p = 0.149). The age range for both sexes in GSE66360 was consistent. The flow chart of this paper is shown in Fig. 1.

Fig. 1
figure 1

Study flow chart

Differentially expressed genes

We cleaned the data, removed the probes without gene annotation, averaged the duplicate probes, normalized and log2-transformed the data, and performed difference analysis for the GSE66360 training set and GSE30529 using the "limma" R package [13]. DEGs were screened based on an adjusted p-value of 0.05 and a |log FC| of > 0.585. The "pheatmap" R package and the "ggplot2" R package were used for volcano plotting and heat map plotting, respectively. Then, we used the "VennDiagram" R package to generate Venn diagrams of up-and-down-regulated genes for the intersecting DEGs in both datasets. Due to the batch effect, different datasets cannot be merged directly, so in our study, we used the fetch intersection method to enhance the robustness of DEG analysis. After removing genes with opposite expression trends, we obtained common DEGs for the two diseases.

Enrichment analyses of common DEGs

GO functional enrichment analysis mainly included three aspects: gene biological process (BP), cellular component (CC), and molecular function (MF). KEGG pathway enrichment is the integration of gene function and genomics as a holistic network to study gene function annotation and gene pathway enrichment. Common DEGs were submitted to enrichment analysis with an adjusted p-value of 0.05 as a screening criterion using the clusterProfiler package [14]. Results are presented using the Ggplot2 package (https://ggplot2.tidyverse.org).

PPI network graph construction, model analysis, and candidate gene screening

The PPI network was constructed by applying the String [15] database (http://string-db.org) with a minimum required interaction score of medium confidence 0.4 as the screening criterion and hiding disconnected nodes in the network. The network information was imported into Cytoscape [16] for visualization, the sub-network maps were constructed by the MCODE module, and the top 15 candidate core genes were screened by the Cytohubba plugin's Degree module. CytoHubba is a tool designed to assist in network topology analysis. Its degree algorithm works by calculating the number of connections that exist between nodes. The more connections a node has in the network, the higher the degree of the node, and proteins with a high degree of degree are more likely to be key proteins.

GO, KEGG, and GeneMANIA analyses of candidate core genes

We performed GO and KEGG analyses of the top 15 candidate core genes screened by the Cytohubba plugin using a p-value of 0.05 as the filtering criterion, and GeneMANIA (http://www.genemania.org/) was used for candidate core gene correlation gene analysis [17].

Connectivity map (CMap) analysis of candidate core genes

The CMap [18] database online resource (https://clue.io) was used to identify drug candidates. This gene expression database uses the gene-expression differences in human cells after treatment with different disruptors to create a database of biological applications where disruptor, gene expression, and disease are interrelated. We submitted the previously screened top 15 candidate core genes to the CMAP website for small molecule drugs that may improve the prognosis of AMI and DN comorbidity. The correlation score between drugs and DEGs is expressed on a scale of -100 to 100, with negative values representing interfering gene expression patterns, indicating that this interference has potential therapeutic effects. We selected < -80 as the cutoff value to screen drug candidates and selected the top 10 to be displayed.

Prediction and verification of transcription factors

We performed TF enrichment analysis of candidate core genes using TRRUST version 2 online software (https://www.grnpedia.org/trrust/) and imported the obtained TF enrichment results into Cytoscape to obtain the regulatory network maps of TFs and genes. We used the limma package for differential analysis of TF and visualization of TF expression data in the GSE66360 training set and GSE30528 dataset.

Validation and analysis of hub genes

Both external datasets, GSE66360 (AMI) validation set and GSE30528 (DN), were used to validate the expression of candidate hub genes. The data were analyzed using the limma package, the data were log2 processed and corrected, the data were compared between groups with a p-value of < 0.05 as the criterion for significant differences, and the visualization of the validation results was performed with the Ggpubr package. Candidate genes that passed validation were considered pivotal genes.

Evaluation of immune infiltrating cells

The fraction of co-expressed genes in AMI circulating endothelial cells and DN renal tubular cells in 28 immune cell types was evaluated separately using the single sample gene set enrichment analysis (ssGSEA) method. The ssGSEA analysis was completed using the "GSVA" and "GSEABase" packages in R, and the results were visualized using the ggplot2 package. The ssGSEA is a widely used technique for analyzing immune cell infiltration, where it uses gene expression data to compare each sample to a specific set of immune cell genes. This comparison results in an immune cell content score, which is calculated using the GSVA function. The higher the score, the higher the immune cell content in the sample. The scoring and grouping data are then combined and corrected, and a heat map is generated as the final output.

Results

Identification of DEGs

Difference analysis and visualization of the GSE66360 training set and the GSE30529 dataset were performed using an adjusted p-value of 0.05 and a |log FC| of > 0.585 as filtering criteria (Fig. 2A–D), and 922 DEGs were obtained in the GSE66360 validation set and 719 differential genes in the GSE30529. A total of 43 up-regulated genes and 3 down-regulated genes were obtained after taking the intersection of the two disease DEGs (Fig. 2E, F).

Fig. 2
figure 2

Volcano and heat plots for difference analysis and co-morbidity intersection Venn diagrams. Heat plot and volcano plot for difference analysis of the GSE66360 training set (A, B); heat plot and volcano plot for difference analysis of the GSE30529 data set (C, D); Venn plots of up-regulated DEG intersections and Venn plots of down-regulated DEG intersections for the GSE66360 training set and the GSE30529 data set (E, F)

GO and KEGG enrichment analysis of common DEGs

The GO functional enrichment analysis of the DEGs with co-morbidities enriched a total of 167 results, and we selected the top 10 results with more significant enrichment of each part to show (Fig. 3 and Supplementary Material S1). KEGG pathway enrichment analysis enriched a total of 17 pathway results that matched the adjusted P of 0.05 (Fig. 3 and Supplementary Material S1).

Fig. 3
figure 3

Bar and bubble plots of GO/KEGG enrichment analysis of common DEGs and circle plots of GO analysis; GO functional analysis of common DEGs (A, B); circle plots of GO analysis (C); KEGG pathway analysis of common DEGs (D, E)

PPI network construction and candidate hub gene screening

A PPI network of common DEGs with 27 nodes and 109 edge interaction matches (medium confidence = 0.4; PPI enrichment p-value of 1.0e-16) was constructed (Fig. 4A-B), and interestingly, all genes in this network graph were upregulated DEGs. Two interaction subnetworks (Fig. 4C–D) were constructed using the MCODE plugin (degree cutoff = 2, node score cutoff = 0.2, set K-core = 2, and max depth = 100) (subnetwork 1: 8 nodes, 25 edges; subnetwork 2: 7 nodes, 15 edges). Fifteen candidate hub genes were filtered using the Degree module of the cytoHubba plugin (Fig. 4E).

Fig. 4
figure 4

PPI network map and sub-network construction of common DEGs and GeneMANIA co-expression network map of candidate core genes. PPI network map (A, B); MCODE module sub-network map (C–D); screening of candidate core genes for cytoHubba plug-in (E); GeneMANIA co-expression gene network diagram (F)

GO and KEGG analysis of candidate core genes

The GO functional enrichment analysis of the DEGs with candidate core genes enriched a total of 86 results, and we selected the top 10 results with more significant enrichment of each part to show. GO functional analysis emphasized the presence of genes mainly in the vesicle membrane and secretory granule membrane involved in antigen processing and presentation, lipopeptide binding, NAD + nucleosidase activity, and Toll-like receptor binding (Fig. 5A, B). KEGG pathway enrichment analysis enriched a total of 14 pathway results that matched the adjusted P of 0.05, and the KEGG pathways analyzed were mainly in the phagosome, neutrophil extracellular trap formation, natural killer cell-mediated cytotoxicity, apoptosis, Fc gamma R-mediated phagocytosis, and Toll-like receptor signaling pathways (Fig. 5C, D). We also used the "circlize" package to plot the chord plot, which can more visually show the interrelationships between genes and gene functions and pathways and share some commonalities, so this chart is ideal for comparing similarities between datasets or different data sets, with nodes distributed around a circle and points connected by arcs to show the relationships, and then assigning weights to each connection (Fig. 5E, F).

Fig. 5
figure 5

Functional and pathway analysis of candidate core gene enrichment. GO function-enriched bar and bubble plots (A, B); bar plot and bubble plot of KEGG pathway enrichment (C–D); GO chord plot (E); KEGG chord plot (F)

Interaction network of candidate core genes and their co-expression genes

We analyzed the interaction network and associated functions of these candidate core genes on the GeneMANIA website. These genes exhibited a complex interaction network with 11.43% physical interactions, 80.09% co-expression, 0.80% prediction, 6.30% co-localization, 2.87% genetic interactions, and 1.37% pathways. The biological functions of this interaction network are mainly involved in phagocytosis, regulation of interleukin-6 production, pattern recognition receptor signaling pathway, positive regulation of defense response, cell surface, regulation of innate immune response, and toll-like receptor signaling pathway (Fig. 4F and Supplementary Material S2).

CMap analysis of candidate core genes

We submitted candidate core genes to the CMap database to identify potential small-molecule chemical drugs for AMI and DN co-morbidity core genes. In order to have a more complete understanding of the drug for the disease, we tentatively set the identification score at -80. We selected the top ten small molecule drugs with an identification score close to -100 and well-defined known chemical structures that are more suitable for experimental and clinical applications. The top 10 potential therapeutic compounds (raltegravir, desoxypeganine, phenytoin, CAY-10577, actarit, OMDM-2, GW-311616, clofibrate, chloroquine, and phylloquinone) (Table 1) and the detailed structures of these compounds are shown in Fig. 6.

Table 1 Results of CMap
Fig. 6
figure 6

To reveal promising chemical drugs, 15 candidate core genes were submitted to the linkage map (CMap) (https://clue.io/), and the top 10 small molecule compounds with the highest enrichment fraction were identified as potentially promising drugs that could have therapeutic effects on AMI and DN comorbidity. 2D structure diagram of small molecule chemical drugs: A raltegravir; B Desoxypeganine; C Phenytoin; D CAY-10577; E Actarit; F OMDM-2; G GW-311616; H Clofibrate; I Chloroquine; J Phylloquinone

TF Exploration and Certification

We explored the TFs using the TRRUST database (Fig. 7). Thereafter, we validated the expression of these TFs in the training set dataset. Three of the six TFs were highly expressed in GSE66360 (Fig. 7A-F), namely HIF1A, NFKB1, and SPI1. SPI1 was highly expressed in the GSE30529 control (Fig. 7G–H). Subsequently, a total of eight hub genes (TLR2, FCGR3B, CD163, CTSS, MNDA, CXCR4, NCF2, CD1D) and six TFs (SPI1, SP1, HIF1A, YY1, RELA, NFKB1) constructed a network diagram of TFs and hub genes (Fig. 7I).

Fig. 7
figure 7

Exploration and validation of TFs and their regulatory networks: A–F Validation of the six relevant TFs predicted using the TRRUST database is expressed in the GSE66360 training set dataset, with HIF1A, NFKB1, and SPI1 highly expressed in the AMI group. G, H TFs were expressed in the GSE30529 dataset, and SPI1 was highly expressed in the DN control group. Statistical significance was determined when the p-value < 0.05. AMI, acute myocardial infarction; DN, diabetic nephropathy. *p < 0.05; **p < 0.01; *** < 0.001. I Regulatory Network of TF: Hub genes are highlighted in red, while TF is highlighted in yellow

Verification of candidate core genes

We validated the expression of candidate core genes in the GSE66360 validation set (AMI) and the GSE30538 dataset (DN). The expression of twelve candidate pivotal genes, TLR2, FCGR2A, FCER1G, CD163, CXCR4, CTSS, CLEC4A, IGSF6, NCF2, FCGR3B, MS4A6A, and GZMB, showed an increasing trend in the GSE6630 validation set (Fig. 8), and the candidate genes TLR2, FCER1G, CD163, TLR1, CTSS, CLEC4A, IGSF6, NCF2, MS4A6A, and CD1D, which were ten genes, showed an increasing trend in the GSE30528 dataset (Fig. 9). Eight AMI and DN co-expressed central genes were finally identified, namely TLR2, FCER1G, CD163, CTSS, CLEC4A, IGSF6, NCF2, and MS4A6A.

Fig. 8
figure 8

Candidate core gene expression levels in GSE66360 A mean t-test was used to compare the differences in core gene expression between the experimental and control groups. Statistical significance was determined when the p-value < 0.05. AMI, acute myocardial infarction; control, healthy control group. *p < 0.05; **p < 0.01; *** < 0.001. ns, p > 0.05

Fig. 9
figure 9

Candidate core gene expression levels in GSE30529. A mean t-test was used to compare the differences in core gene expression between experimental and control groups. Statistical significance was determined when the p-value < 0.05. DN, diabetic nephropathy; Control, healthy control group. *p < 0.05; **p < 0.01; *** < 0.001. ns, p > 0.05

Relationship between co-expressed genes and immune cells

Immune cell content in samples from the AMI training set (GSE66360) and DN training set (GSE30529) was analyzed by gene expression profiling using the GSEABase package and GSVA package. Immunocell heat maps of immune cell expression in the samples were obtained by scoring the immune cell content (Fig. 10A, D). Next, a differential analysis of immune cells was performed to compare the expression of immune cells in the control and test groups (Fig. 10B, E). Finally, we performed immune cell correlation analysis using the Spearman method to explore potential associations between co-expressed genes and immune cells (Fig. 10C, F). Among the eight co-expressed genes and immune correlation analysis in DN, CD163 was associated with Type 2 T helper cells, Type 1 T helper cells, plasmacytoid DC, MDSC, mast cell, macrophage, immature B cell, and activated B. CLEC4A, which is positively correlated with activated CD4 T cells, immature B cells, and regulatory T cells. CTSS was positively correlated with Type 17 T helper cells. FCER1G was positively correlated with activated B cells, activated CD8 T cells, activated dendritic cells, CD56bright NKC, effector memory CD4 T cells, and gamma cells. IGFS6 was positively correlated with activated B cells, activated CD4 T cells, activated CD8 T cells, activated dendritic cells, CD56bright NKC, effector delta T cells, macrophages, monocytes, neutrophils, plasmacytoid DC, T follicular helper cells, Activated CD4 T cells, activated CD8 T cells, central memory CD4 T cells, effector memory CD4 T cells, gamma delta T cells, immature B cells, MDSC, mast cell macrophages, plasmacytoid DC, regulatory T cells, T follicular helper cells, and Type 1 T helper cells. MS4A6A is positively associated with activated B cells, activated CD4 T cells, activated CD8 T cells, effector memory CD4 T cells, MDSC, mast cells, NKC, plasmacytoid DC, and Type 1 T helper cells. NCF2 is positively correlated with CD56-bright natural killer cells, gamma delta T cells, MDSC, NKC, plasmacytoid DC, and Type 1 T helper cells. TLR2 was positively correlated with activated B cells, activated CD4 T cells, activated CD8 T cells, activated dendritic cells, effector memory CD8 T cells, gamma delta T cells, plasmacytoid dendritic cells, and T follicular helper cells. In the immune correlation analysis of AMI, eight co-expressed genes were positively correlated with CD56 dim NKC, CD56 bright NKC, activated DC, eosinophils, gamma delta T cells, and immature dendritic cells, macrophages, mast cells, MDSC, monocytes, NKC, natural killer T cells, neutrophils, plasmacytoid dendritic cells, regulatory T cells, T follicular helper cells, and Type 1 T helper cells; negatively correlated with central memory CD4 T cell, effector memory CD4 T cell, and central and effector memory CD8 T cells.

Fig. 10
figure 10

Correlation between DEGs and co-expressed genes with immune cells (A-C: GSE30529 dataset; D-F: GSE66360 training set). A The correlation between DEGs and immune cells in GSE30529; B Differential analysis of immune cells, Activated B cell, Activated CD4 T cell, Activated CD8 T cell, Activated dendritic cell, CD56bright natural killer cell, Gamma delta T cell, Immature B cell, MDSC, Monocyte, NKC, Mast cell, Regulatory T cell, T follicular helper cell, Type 1 T helper cell, Type 2 T helper cell, Central memory CD4 T cells were highly expressed in the experimental group, and central memory CD8 T cells were highly expressed in the control group; C The correlation between co-expressed genes and immune cells (GSE30529); D The correlation between DEGs and immune cells in GSE66360; E Differential analysis of immune cells: activated CD4 T cells, activated dendritic cells, CD56 dim natural killer cell, eosinophils,gamma delta T cell, immature dendritic cell, MDSC, macrophage, mast cell, monocyte, natural killer T cell, natural Neutrophil, Plasmacytoid dendritic cell, Regulatory T cell, T follicular helper cell, and Type 1 T helper cell were highly expressed in the experimental group, and central memory CD4 T cells were highly expressed in the control group; F, the correlation between co-expressed genes and immune cells (GSE66360 training set)

Discussion

CHD is the most common cause of death worldwide, causing 8.9 million deaths in 2019 alone [19]. AMI is a potentially fatal presentation of CHD and is a main component of CHD deaths [20], and its complications can cause a serious economic and social burden. DN is essentially a chronic microvascular complication of diabetes and a major cause of ESRD, the overall prevalence of DN in China is 8.95% [21]. AMI and DN co-morbidities are very commonly seen clinically and have a poor prognosis, and bioinformatics studies of their co-morbidities have not been reported. Our study provides some contributions to the identification of novel biomarkers of susceptibility to diabetic nephropathy and AMI co-morbidity and useful potential therapeutic targets. Transcriptome difference analysis relies on three main R packages: DESeq2, edgeR, and limma. These packets extract large differences in DGEs, thus enhancing the identification of DEGs using intersection analysis methods. However, DESeq2 and edgeR are suitable for analyzing small data samples, and the false positive rate increases significantly with large samples. In this study, we used the most commonly used limma package for analysis to improve the accuracy of DEGs and enrichment analysis.

In this study, we report the following bioinformatics studies on co-morbidities: We identified and validated eight co-morbid genes, namely TLR2, FCER1G, CD163, CTSS, CLEC4A, IGSF6, NCF2, and MS4A6A. TLR2 is a gene encoding a Toll-like receptor family of proteins, and intervention in the TLR2 signaling pathway may be beneficial in the treatment of DN [22]. TLR2 activates inflammatory responses by activating the TLR2-Myd88-NF-κB pathway, also activates immune cells to promote apoptosis [23], and is associated with the severity of coronary artery disease, plaque instability, and myocardial infarction in myocardial injury [24]. High expression of FCER1G in the plasma of patients with plaque rupture [25] may be associated with atheromatous plaque instability by unknown mechanisms. CD163 is expressed only in monocytes and macrophages, and levels of the soluble form of CD163 (sCD163) are elevated in low-grade inflammation [26], which may be involved in chronic inflammatory responses and possibly in AMI and DN development through the induction of foam cell formation and local inflammatory responses. CTSS degrades damaged or unwanted proteins in lysosomes, regulates antigen presentation in atherosclerosis [27], and mediates phagosome involvement in atherosclerosis progression through macrophages [28]. CTSS is also a biomarker of chronic kidney disease, and as the glomerular filtration rate decreases, CTSS increases, exacerbating inflammation-associated endothelial dysfunction [29]. CTSS is closely associated with diabetes and complications (cardiovascular disease, retinopathy, diabetic nephropathy, and other microangiopathies) [30], impairs endothelial function through PAR2 activation, and promotes podocyte reduction [31]. CLEC4A encodes C-type lectin, which is involved in cell adhesion, intercellular signaling, glycoprotein renewal, and inflammatory and immune responses (https://www.ncbi.nlm.nih.gov/gene/). IGSF6, immunoglobulin superfamily member 6, may be a plasma membrane component involved in cross-modal signaling. NCF2 may be involved in iron death in cells from DN patients [32]. MS4A6A is recognized to be associated with aging and neurodegenerative diseases, and bioinformatic analysis has revealed its high expression in DN, but the exact mechanism is unknown.

GO function is enriched in the vesicle membrane and secretory membrane and is involved in antigen processing and presentation, lipopeptide binding, NAD+ nucleosidase activity, and Toll-like receptor binding. Atherosclerosis is a chronic inflammatory vascular disease [33]; both AMI and DN are essentially atherosclerotic diseases; and most immune cells are involved in the atherosclerotic process [34]. Macrophages are considered the main immune cell type in atherosclerotic plaques [35], acting as antigen-presenting cells and presenting antigens to naive CD8 + T cells and naive CD4 + T cells [36]. CD8 + killer T cells induce apoptosis and necrosis of target cells, exacerbating progression in atherosclerotic plaques [37], while CD4T cell subsets can influence atherosclerotic progression through immune activation or immunosuppression [38]. Macrophages can also be involved in atherosclerotic plaque formation through toll-like receptor activation by oxLDL [39]. Toll-like receptor 4 (TLR4), a cellular inflammatory factor, can induce the expression and release of acquired immune-induced inflammatory factors by initiating natural immunity and stimulating the production of cytokines and chemokines through the corresponding signaling pathways. release, which in turn induces DN and atherogenesis [40]. Macrophage infiltration around renal tubular epithelial cells is considered a hallmark of DN, and DN or lipotoxicity promotes the increased release of extracellular vesicles, which were found to induce apoptosis of injured renal tubular epithelial cells through a death receptor 5 (DR5)-dependent process [41]. Nicotinamide adenine dinucleotide (NAD+) is a carrier of high-energy electrons generated by oxidative metabolism in vivo and serves as a cofactor for most oxidative enzymes with redox and signaling functions. NAD+ dynamic balance has a key role in cellular catabolism and anabolism and is also an important cofactor for redox reactions in the cytoplasm and mitochondria, regulating a variety of cellular processes, including the cellular stress response. It is also an important cofactor in the cytoplasm and mitochondrial redox reactions, regulating various cellular processes including cellular stress response, mitochondrial homeostasis, and calcium signaling. The kidney is second only to the heart in the number of mitochondria and has a high oxidative metabolism and a high demand for NAD+ . NAD+ is a determinant of aging [42], and alterations in NAD+ levels have been detected in diabetic nephropathy in studies suggesting a possible involvement in the development of diabetic nephropathy, but no assessment of NAD+ distribution in the kidney is available [43]. Studies have shown that NAD+ supplementation reduces myocardial infarction size during cardiac ischemia through the longevity protein Sirt1. This process, moreover, relies on the longevity protein Sirt1 to act [44]. A follow-up study found that NAD + supplementation with schizomycin inhibited longevity protein activity but still reduced myocardial infarction in mice, suggesting that the cardioprotective effect of NAD+ supplementation is not only dependent on longevity protein. Further studies confirmed [45] that NAD+ supplementation activates intracellular glycolysis and thus generates large amounts of energy, which acts to protect the myocardium and reduce myocardial infarction through glycolysis in the presence of sufficient glucose in the body.

KEGG pathway enrichment focuses on the phagosome, neutrophil extracellular trap formation, natural killer cell-mediated cytotoxicity, apoptosis, Fc gamma R-mediated phagocytosis, and Toll-like receptor signaling pathways. Pathway enrichment analysis suggests that FcγR-mediated phagocytosis and Toll-like receptor signaling pathways are involved in the pathogenesis of AMI and DN. FcγR, a receptor for the Fc portion of IgG, activates macrophages by mediating the activation of the mitogen-activated protein kinase signaling pathway by low-density lipoprotein immune complexes (LDL-ICs) [46]. Macrophage polarization toward the M1 phenotype increases the risk of AMI pathogenesis and rupture and also the formation of atherosclerotic plaques by taking up LDL-ICs and converting them into foam cells [47]. FcγR activates and promotes phagocytosis of macrophages, induces foam cell aggregation, and promotes the development of atherosclerotic plaques [23], inducing and exacerbating AMI and DN. Neutrophil The formation of extracellular traps (NETs) in granulocytes can be triggered by various stimuli induced by contact with surface receptors (cytokine receptors, Fcγ receptors, toll-like receptors, complement receptors, etc.) and adenosine receptors, causing an increase in intracellular Ca2 + concentration, which activates protein kinase c (PKC), leading to reactive oxygen species (ROS) through NADPH oxidase and/or mitochondrial activation production, resulting in organismal damage [48].

Meanwhile, we screened small molecule drugs for co-morbidities through the CMap (https://clue.io/) website and selected the top 10 small molecule drugs, namely raltegravir, desoxypeganine, phenytoin, CAY-10577, actarit, OMDM-2, GW-311616, clofibrate, chloroquine, and phylloquinone, which may become new approaches for the treatment of AMI and DN co-morbidities. As the first integrase inhibitor approved for the treatment of HIV infection, Raltegravir is gradually being discovered to have effects other than anti-HIV. Studies have shown that Raltegravir may be effective in treating multiple myeloma by promoting DNA damage-induced apoptosis [49]. Another study showed that Raltegravir was effective in treating human T-cell lymphoblastic virus-1 (HTLV-1)-associated myelopathy/tropical spastic paraplegia (HAM/TSP). Although it did not significantly reduce the high proviral load, it produced therapeutic effects on this chronic progressive disease by reducing the expression of CD4+ CD25+ T cells, antibody-secreting B cells, activated CD8+ T cells, and IFN-γ [50]. We speculate that Raltegravir may play an anti-arterial inflammatory role by regulating inflammatory cells and inflammatory factors, but no relevant reports have been reported to date. Desoxypeganine (DOP) is a natural alkaloid. It is not only a cholinesterase inhibitor but also a selective inhibitor of monoamine oxidase A. Oral DOP is well tolerated and safe for smoking and alcohol cessation [51]. Cholinesterase plays a key role in lipid metabolism and is also a predisposing factor for metabolic diseases [52]. Cholinesterase activity is usually higher in diabetic patients than in hyperlipidemia patients [53]. Galantamine, a cholinesterase inhibitor, can reduce inflammation and insulin resistance in metabolic syndrome and improve diabetic kidney damage and renal function [54]. The cholinesterase inhibitor pyridostigmine can enhance cardiac parasympathetic tone and induce an anti-inflammatory response in the heart, thereby promoting cardiac recovery after AMI [55]. High-density lipoprotein (HDL) cholesterol is negatively associated with coronary heart disease, and phenytoin, a drug that can raise HDL levels [56], may be related to HDL's involvement in atherosclerosis mechanisms. Actarit is a drug for the treatment of rheumatoid arthritis, and its main target of action is carbonic anhydrase II [57]. This target is very important in the kidney and osteoclasts, and its main function is to maintain acid-base balance in the body as well as normal bone remodeling [58]. At the same time, Actarit's research in the treatment of heart disease and other diseases is also being repositioned [57]. OMDM-2 is a cannabinoid transporter inhibitor. Angiotensin II releases endocannabinoids by activating the AT1 receptor [59]. The mechanism of the cardiovascular system's response to endocannabinoids is complex, and there are still many unsolved mysteries. It has been suggested that cannabinoid-1 receptor (CB1R) activation is associated with oxidative low-density lipoprotein formation and inflammatory response [60], while CB2R agonists attenuate the inflammatory response of human coronary endothelial cells caused by CB1R activation [61]. Overactivation of CB1R in the renal endocannabinoid system contributes to DN development. Moreover, CB1R blockers have been shown to improve renal insufficiency in mouse models of diabetes [62]. In podocytes, CB1R deletion alleviates glomerular and tubular dysfunction in mouse models of diabetic nephropathy [63]. Animal experiments have shown that neutrophil elastase (NE) is involved in the formation of foam cells and the progression of atherosclerosis, and the specific neutrophil elastase inhibitor GW-311616 significantly reduces atherosclerosis [64]. Clofibrate is a peroxisome proliferator-activated receptor-alpha (PPARα) agonist that reduces the inflammatory response and modulates the renin-angiotensin system. One week of Clofibrate treatment reduced angiotensin II and the angiotensin-converting enzyme, AT1, while bradykinin, ACE-2, and AT2 receptors were elevated, nitric oxide levels increased, and fibrosis and structural damage caused by myocardial infarction improved [65]. Clofibrate also reduced glomerular transforming growth factor (TGF)-β1-induced increased albumin permeability in DN rats [66] or may be used for early treatment of DN. Chloroquine (CQ) is an autophagy inhibitor that inhibits lipid deposition-related kidney injury mediated by autophagy deficiency in DN [67]. Studies have shown that hydroxychloroquine has a role in reducing the risk of cardiovascular disease events in patients with traditional risk factors [68], while chloroquine diphosphate promotes caspase-3 and caspase-9 activity and increases AMI cardiomyocyte apoptosis [69]. Vitamin K is an essential bioactive substance required by the body. According to its main structure, it is divided into phylloquinones (K1) and menaquinones (K2). Phylloquinone is mainly derived from leafy green vegetables and can be converted into menaquinone-4 (MK4) in vivo. It can effectively inhibit lipid peroxidation and its mediated iron death. It can effectively prevent tissue damage in the model of hepatic and renal ischemia–reperfusion injury. It may be one of the oldest natural iron-resistant quinones [70]. K2 is potent in preventing arterial calcification, reducing the risk of diabetes, and improving kidney function [71]. In general, desoxypeganine, actarit, OMDM-2, GW-311616, clofibrate, and phylloquinones may play a greater role in the treatment of DN and coronary heart disease.

We also constructed a network map consisting of eight hub genes (TLR2, FCGR3B, CD163, CTSS, MNDA, CXCR4, NCF2, and CD1D) and six transcription factors (SPI1, SP1, HIF1A, YY1, RELA, and NFKB1) using the TRRUST dataset and examined the TF expression in the training set. A total of three AMI-related transcription factors (HIF1A, NFKB1, and SPI1) and one DN-related TF (SPI1) were identified, suggesting that SPI1 may play opposite roles in the two diseases, a conclusion that needs to be further confirmed in external data. To validate transcription factor expression, we initially constructed myocardial infarction models and diabetic nephropathy models to confirm whether transcription factor expression in the models was consistent with bioinformatics predictions. The procedure begins with bioinformatics analysis that screens for conserved transcription factor sequences, transcription factor core regions, signal peptides, protein glycosylation checkpoints, phosphorylation checkpoints, protein kinase binding checkpoints, and protein affinity/hydrophobicity. Transcription factors are thoroughly examined for their overall transcriptional activation using a dual luciferase detection system. Gene expression is further verified by overexpressing or silencing transcription factors, and the regulatory relationship between transcription factors and target genes is elucidated in relation to phenotype. Protein-DNA interactions are used to identify the DNA sequence to which the transcription factor binds, and chromatin immunoprecipitation (TF-CHIP) is used to determine whether the transcription factor binds to gene promoter sequences. This process will be further refined and completed in future work.

Our findings contribute to the clarification of the mechanisms underlying the association between AMI and DN. However, the study still has some shortcomings. Our findings need further confirmation with external data; core genes need further experimental validation in vitro models; and some of the genes enriched are not yet reported in the literature on AMI or DN. Supplementing these shortcomings will be the focus of our future studies.

Conclusion

We identified and validated eight co-expressed central genes and three TFs of AMI and DN. The preliminary screening of small-molecule chemotherapeutics for co-morbidities provided new directions for co-morbidity understanding and research and new ideas for treatment. TLR2, FCER1G, CD163, CTSS, CLEC4A, IGSF6, NCF2, and MS4A6A are eight co-expressed central genes that were enriched in an analysis focusing on NAD+ nucleosidase activity and Toll-like receptor binding in the phagosome. neutrophil extracellular trap formation, natural killer cell-mediated cytotoxicity, apoptosis, Fc gamma R-mediated phagocytosis, and Toll-like receptor signaling pathways.

Availability of data and materials

The current dataset analyzed in the study can be found in the GEO database.

GSE66360: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE66360.

GSE30528: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE30528.

GSE30529: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE30529.

References

  1. Miura T, Kuno A, Tanaka M. Diabetes modulation of the myocardial infarction-acute kidney injury axis. Am J Physiol Heart Circ Physiol. 2022;322(3):H394–405. https://doi.org/10.1152/ajpheart.00639.2021. (Epub 2022 Jan 28 PMID: 35089809).

    Article  CAS  PubMed  Google Scholar 

  2. Kenny HC, Abel ED. Heart failure in type 2 diabetes mellitus. Circ Res. 2019;124(1):121–41. https://doi.org/10.1161/CIRCRESAHA.118.311371.PMID:30605420;PMCID:PMC6447311.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Wang ZV, Hill JA. Diabetic cardiomyopathy: catabolism driving metabolism. Circulation. 2015;131(9):771–3. https://doi.org/10.1161/CIRCULATIONAHA.115.015357. Epub 2015 Jan 30. PMID: 25637626; PMCID: PMC4351172.

    Article  PubMed  PubMed Central  Google Scholar 

  4. GBD Chronic Kidney Disease Collaboration. Global, regional, and national burden of chronic kidney disease, 1990–2017: A systematic analysis for the global burden of disease stud. Lancet (London, Engl). 2020;2017(395):709. https://doi.org/10.1016/S0140-6736(20)30045-3.

    Article  Google Scholar 

  5. Smith GL, Masoudi FA, Shlipak MG, Krumholz HM, Parikh CR. Renal impairment predicts long-term mortality risk after acute myocardial infarction. J Am Soc Nephrol. 2008;19(1):141–50. https://doi.org/10.1681/ASN.2007050554. Epub 2007 Nov 14. PMID: 18003773; PMCID: PMC2391037.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Stolker JM, McCullough PA, Rao S, Inzucchi SE, Spertus JA, Maddox TM, Masoudi FA, Xiao L, Kosiborod M. Pre-procedural glucose levels and the risk for contrast-induced acute kidney injury in patients undergoing coronary angiography. J Am Coll Cardiol. 2010;55(14):1433–40. https://doi.org/10.1016/j.jacc.2009.09.072. (PMID: 20359592).

    Article  CAS  PubMed  Google Scholar 

  7. Abe M, Morimoto T, Akao M, et al. Relation of contrast-induced nephropathy to long-term mortality after percutaneous coronary intervention. Am J Cardiol. 2014;114:362–8. https://doi.org/10.1016/j.amjcard.2014.05.009. (PMID: 24927973).

    Article  PubMed  Google Scholar 

  8. Tsai TT, Patel UD, Chang TI, et al. Contemporary incidence, predictors, and outcomes of acute kidney injury in patients undergoing percutaneous coronary interventions: insights from the NCDR Cath-PCI registry. JACC Cardiovasc Interv. 2014;7:1–9. https://doi.org/10.1016/j.jcin.2013.06.016. (PMID: 24456715).

    Article  PubMed  PubMed Central  Google Scholar 

  9. Ueda H, Yamada T, Masuda M, Okuyama Y, Morita T, Furukawa Y, Koji T, Iwasaki Y, Okada T, Kawasaki M, Kuramoto Y, Naito T, Fujimoto T, Komuro I, Fukunami M. Prevention of contrast-induced nephropathy by bolus injection of sodium bicarbonate in patients with chronic kidney disease undergoing emergent coronary procedures. Am J Cardiol. 2011;107(8):1163–7. https://doi.org/10.1016/j.amjcard.2010.12.012. (Epub 2011 Feb 23 PMID: 21349483).

    Article  CAS  PubMed  Google Scholar 

  10. Marenzi G, Assanelli E, Campodonico J, Lauri G, Marana I, De Metrio M, Moltrasio M, Grazi M, Rubino M, Veglia F, Fabbiocchi F, Bartorelli AL. Contrast volume during primary percutaneous coronary intervention and subsequent contrast-induced nephropathy and mortality. Ann Intern Med. 2009;150(3):170–7. https://doi.org/10.7326/0003-4819-150-3-200902030-00006. (PMID: 19189906).

    Article  PubMed  Google Scholar 

  11. Tan JK, Salim NNM, Lim GH, Chia SY, Thumboo J, Bee YM. Trends in diabetes-related complications in Singapore, 2013–2020: A registry-based study. PLoS ONE. 2022;17(10): e0275920. https://doi.org/10.1371/journal.pone.0275920.PMID:36219616;PMCID:PMC9553054.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Tan KW, Dickens BSL, Cook AR. Projected burden of type 2 diabetes mellitus-related complications in Singapore until 2050: a Bayesian evidence synthesis. BMJ Open Diabetes Res Care. 2020;8(1): e000928. https://doi.org/10.1136/bmjdrc-2019-000928. (PMID: 32184203; PMCID: PMC7076230).

    Article  PubMed  PubMed Central  Google Scholar 

  13. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7) Epub 2015 Jan 20. PMID: 25605792; PMCID: PMC4402510.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7. https://doi.org/10.1089/omi.2011.0118. Epub 2012 Mar 28. PMID: 22455463; PMCID: PMC3339379.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P, Jensen LJ, Mering CV. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13. https://doi.org/10.1093/nar/gky1131.PMID:30476243;PMCID:PMC6323986.

    Article  CAS  PubMed  Google Scholar 

  16. Smoot ME, Ono K, Ruscheinski J, Wang PL, Ideker T. Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2011;27(3):431–2. https://doi.org/10.1093/bioinformatics/btq675. Epub 2010 Dec 12. PMID: 21149340; PMCID: PMC3031041.

    Article  CAS  PubMed  Google Scholar 

  17. Mostafavi S, Ray D, Warde-Farley D, Grouios C, Morris Q. GeneMANIA: a real-time multiple association network integration algorithm for predicting gene function. Genome Biol. 2008;9 Suppl 1(Suppl 1):S4. https://doi.org/10.1186/gb-2008-9-s1-s4. Epub 2008 Jun 27. PMID: 18613948; PMCID: PMC2447538.

  18. Musa A, Ghoraie LS, Zhang SD, Glazko G, Yli-Harja O, Dehmer M, Haibe-Kains B, Emmert-Streib F. A review of connectivity map and computational approaches in pharmacogenomics. Brief Bioinform. 2018;19(3):506–23. https://doi.org/10.1093/bib/bbw112.Erratum.In:BriefBioinform.2017Sep1;18(5):903.PMID:28069634;PMCID:PMC5952941.

    Article  CAS  PubMed  Google Scholar 

  19. Desai R, Mishra V, Chhina AK, Jain A, Vyas A, Allamneni R, Lavie CJ, Sachdeva R, Kumar G. Cardiovascular Disease Risk Factors and Outcomes of Acute Myocardial Infarction in Young Adults: Evidence from Two Nationwide Cohorts in the United States a Decade Apart. Curr Probl Cardiol. 2023:101747. https://doi.org/10.1016/j.cpcardiol.2023.101747. Epub ahead of print. PMID: 37087077.

  20. Virani SS, Alonso A, Benjamin EJ, et al. Heart Disease and Stroke Statistics-2020 Update: A Report From the American Heart Association. Circulation. 2020;141(9):e139–596. https://doi.org/10.1161/CIR.0000000000000757. (Epub 2020 Jan 29 PMID: 31992061).

    Article  PubMed  Google Scholar 

  21. Zhang XX, Kong J, Yun K. Prevalence of diabetic nephropathy among patients with type 2 diabetes mellitus in China: a meta-analysis of observational studies. J Diabetes Res. 2020;3(2020):2315607. https://doi.org/10.1155/2020/2315607.PMID:32090116;PMCID:PMC7023800.

    Article  Google Scholar 

  22. Wada J, Makino H. Innate immunity in diabetes and diabetic nephropathy. Nat Rev Nephrol. 2016;12(1):13–26. https://doi.org/10.1038/nrneph.2015.175. (Epub 2015 Nov 16 PMID: 26568190).

    Article  CAS  PubMed  Google Scholar 

  23. Zhang Q, Liu H, Zhang M, Liu F, Liu T. Identification of co-expressed central genes and transcription factors in atherosclerosis-related intracranial aneurysm. Front Neurol. 2023;2(14):1055456. https://doi.org/10.3389/fneur.2023.1055456.PMID:36937519;PMCID:PMC10017537.

    Article  Google Scholar 

  24. Qi B, Chen JH, Tao L, Zhu CM, Wang Y, Deng GX, Miao L. Integrated Weighted Gene Co-expression Network Analysis Identified That TLR2 and CD40 Are Related to Coronary Artery Disease. Front Genet. 2021;26(11): 613744. https://doi.org/10.3389/fgene.2020.613744.PMID:33574831;PMCID:PMC7870792.

    Article  Google Scholar 

  25. Xiao S, Zhou Y, Wu Q, Liu Q, Chen M, Zhang T, Zhu H, Liu J, Yin T, Pan D. FCER1G and PTGS2 serve as potential diagnostic biomarkers of acute myocardial infarction based on integrated bioinformatics analyses. DNA Cell Biol. 2021;40(8):1064–75. https://doi.org/10.1089/dna.2020.6447. (Epub 2021 Jun 10 PMID: 34115526).

    Article  CAS  PubMed  Google Scholar 

  26. Møller HJ, Frikke-Schmidt R, Moestrup SK, Nordestgaard BG, Tybjærg-Hansen A. Serum soluble CD163 predicts risk of type 2 diabetes in the general population. Clin Chem. 2011;57(2):291–7. https://doi.org/10.1373/clinchem.2010.154724. (Epub 2010 Nov 24 PMID: 21106861).

    Article  CAS  PubMed  Google Scholar 

  27. Borghese F, Clanchy FI. CD74: an emerging opportunity as a therapeutic target in cancer and autoimmune disease. Expert Opin Ther Targets. 2011;15(3):237–51. https://doi.org/10.1517/14728222.2011.550879. (Epub 2011 Jan 6 PMID: 21208136).

    Article  CAS  PubMed  Google Scholar 

  28. Wang H, Jiang H, Cheng XW. Cathepsin S are involved in human carotid atherosclerotic disease progression, mainly by mediating phagosomes: bioinformatics and in vivo and vitro experiments. PeerJ. 2022;8(10): e12846. https://doi.org/10.7717/peerj.12846.PMID:35186462;PMCID:PMC8833225.

    Article  Google Scholar 

  29. Sonkar SK, Singh PK, Chandra S, Sonkar GK, Bhosale V, Sharma S. Cathepsin S as an early biomarker for cardiovascular disease in chronic kidney disease patients. J Bras Nefrol. 2022;44(3):329–35. https://doi.org/10.1590/2175-8239-JBN-2021-0135. PMID: 35023538; PMCID: PMC9518627.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Smyth P, Sasiwachirangkul J, Williams R, Scott CJ. Cathepsin S (CTSS) activity in health and disease - A treasure trove of untapped clinical potential. Mol Aspects Med. 2022;88: 101106. https://doi.org/10.1016/j.mam.2022.101106. (Epub 2022 Jul 19 PMID: 35868042).

    Article  CAS  PubMed  Google Scholar 

  31. Kumar Vr S, Darisipudi MN, Steiger S, Devarapu SK, Tato M, Kukarni OP, Mulay SR, Thomasova D, Popper B, Demleitner J, Zuchtriegel G, Reichel C, Cohen CD, Lindenmeyer MT, Liapis H, Moll S, Reid E, Stitt AW, Schott B, Gruner S, Haap W, Ebeling M, Hartmann G, Anders HJ. Cathepsin S Cleavage of Protease-Activated Receptor-2 on Endothelial Cells Promotes Microvascular Diabetes Complications. J Am Soc Nephrol. 2016;27(6):1635–49. https://doi.org/10.1681/ASN.2015020208. Epub 2015 Nov 13. PMID: 26567242; PMCID: PMC4884104.

    Article  PubMed  Google Scholar 

  32. Ma LL, Bai Y, Liu WH, Diao ZL. Bioinformatics analysis of potential key ferroptosis-related genes involved in tubulointerstitial injury in patients with diabetic nephropathy. Ren Fail. 2023;45(1):2199095. https://doi.org/10.1080/0886022X.2023.2199095.PMID:37038746;PMCID:PMC10101677.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Chistiakov DA, Melnichenko AA, Myasoedova VA, Grechko AV, Orekhov AN. Mechanisms of foam cell formation in atherosclerosis. J Mol Med (Berl). 2017;95(11):1153–65. https://doi.org/10.1007/s00109-017-1575-8. (Epub 2017 Aug 7 PMID: 28785870).

    Article  CAS  PubMed  Google Scholar 

  34. Roy P, Orecchioni M, Ley K. How the immune system shapes atherosclerosis: roles of innate and adaptive immunity. Nat Rev Immunol. 2022;22(4):251–65. https://doi.org/10.1038/s41577-021-00584-1. Epub 2021 Aug 13. PMID: 34389841; PMCID: PMC10111155.

    Article  CAS  PubMed  Google Scholar 

  35. Willemsen L, de Winther MP. Macrophage subsets in atherosclerosis as defined by single-cell technologies. J Pathol. 2020;250(5):705–14. https://doi.org/10.1002/path.5392. Epub 2020 Mar 11. PMID: 32003464; PMCID: PMC7217201.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Fang F, Xiao C, Li C, Liu X, Li S. Tuning macrophages for atherosclerosis treatment. Regen Biomater. 2022;10:rbac103. https://doi.org/10.1093/rb/rbac103. PMID: 36683743; PMCID: PMC9845526.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Kyaw T, Tipping P, Toh BH, Bobik A. Killer cells in atherosclerosis. Eur J Pharmacol. 2017;5(816):67–75. https://doi.org/10.1016/j.ejphar.2017.05.009. (Epub 2017 May 5 PMID: 28483458).

    Article  CAS  Google Scholar 

  38. Saigusa R, Winkels H, Ley K. T cell subsets and functions in atherosclerosis. Nat Rev Cardiol. 2020;17(7):387–401. https://doi.org/10.1038/s41569-020-0352-5. Epub 2020 Mar 16. PMID: 32203286; PMCID: PMC7872210.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Edfeldt K, Swedenborg J, Hansson GK, Yan ZQ. Expression of toll-like receptors in human atherosclerotic lesions: a possible pathway for plaque activation. Circulation. 2002;105(10):1158–61 (PMID: 11889007).

    Article  CAS  PubMed  Google Scholar 

  40. Lin M, Yiu WH, Wu HJ, Chan LY, Leung JC, Au WS, Chan KW, Lai KN, Tang SC. Toll-like receptor 4 promotes tubular inflammation in diabetic nephropathy. J Am Soc Nephrol. 2012;23(1):86–102. https://doi.org/10.1681/ASN.2010111210. Epub 2011 Oct 21. PMID: 22021706; PMCID: PMC3269929.

    Article  CAS  PubMed  Google Scholar 

  41. Jiang WJ, Xu CT, Du CL, Dong JH, Xu SB, Hu BF, Feng R, Zang DD, Meng XM, Huang C, Li J, Ma TT. Tubular epithelial cell-to-macrophage communication forms a negative feedback loop via extracellular vesicle transfer to promote renal inflammation and apoptosis in diabetic nephropathy. Theranostics. 2022;12(1):324–39. https://doi.org/10.7150/thno.63735.PMID:34987648;PMCID:PMC8690920.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Nacarelli T, Lau L, Fukumoto T, Zundell J, Fatkhutdinov N, Wu S, Aird KM, Iwasaki O, Kossenkov AV, Schultz D, Noma KI, Baur JA, Schug Z, Tang HY, Speicher DW, David G, Zhang R. NAD+ metabolism governs the proinflammatory senescence-associated secretome. Nat Cell Biol. 2019;21(3):397–407. https://doi.org/10.1038/s41556-019-0287-4. Epub 2019 Feb 18. PMID: 30778219; PMCID: PMC6448588.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Zapata-Pérez R, Wanders RJA, van Karnebeek CDM, Houtkooper RH. NAD+ homeostasis in human health and disease. EMBO Mol Med. 2021;13(7):e13943. https://doi.org/10.15252/emmm.202113943. Epub 2021 May 27. PMID: 34041853; PMCID: PMC8261484.

  44. Yamamoto T, Byun J, Zhai P, Ikeda Y, Oka S, Sadoshima J. Nicotinamide mononucleotide, an intermediate of NAD+ synthesis, protects the heart from ischemia and reperfusion. PLoS ONE. 2014;9(6): e98972. https://doi.org/10.1371/journal.pone.0098972.PMID:24905194;PMCID:PMC4048236.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Nadtochiy SM, Wang YT, Nehrke K, Munger J, Brookes PS. Cardioprotection by nicotinamide mononucleotide (NMN): Involvement of glycolysis and acidic pH. J Mol Cell Cardiol. 2018;121:155–62. https://doi.org/10.1016/j.yjmcc.2018.06.007. Epub 2018 Jun 26. PMID: 29958828; PMCID: PMC6103815.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Virella G, Muñoz JF, Galbraith GM, Gissinger C, Chassereau C, Lopes-Virella MF. Activation of human monocyte-derived macrophages by immune complexes containing low-density lipoprotein. Clin Immunol Immunopathol. 1995;75(2):179–89. https://doi.org/10.1006/clin.1995.1069. (PMID: 7704977).

    Article  CAS  PubMed  Google Scholar 

  47. Muhammad S, Chaudhry SR, Dobreva G, Lawton MT, Niemelä M, Hänggi D. Vascular macrophages as therapeutic targets to treat intracranial aneurysms. Front Immunol. 2021;8(12): 630381. https://doi.org/10.3389/fimmu.2021.630381.PMID:33763073;PMCID:PMC7982735.

    Article  Google Scholar 

  48. Wigerblad G, Kaplan MJ. Neutrophil extracellular traps in systemic autoimmune and autoinflammatory diseases. Nat Rev Immunol. 2022:1–15. https://doi.org/10.1038/s41577-022-00787-0. Epub ahead of print. PMID: 36257987; PMCID: PMC9579530.

  49. Akcora-Yildiz D, Gonulkirmaz N, Ozkan T, Beksac M, Sunguroglu A. HIV-1 integrase inhibitor raltegravir promotes DNA damage-induced apoptosis in multiple myeloma. Chem Biol Drug Des. 2023;102(2):262–70. https://doi.org/10.1111/cbdd.14237. (Epub 2023 Apr 24 PMID: 37094820).

    Article  CAS  PubMed  Google Scholar 

  50. Enose-Akahata Y, Billioux BJ, Azodi S, Dwyer J, Vellucci A, Ngouth N, Nozuma S, Massoud R, Cortese I, Ohayon J, Jacobson S. Clinical trial of raltegravir, an integrase inhibitor. HAM/TSP Ann Clin Transl Neurol. 2021;8(10):1970–85. https://doi.org/10.1002/acn3.51437. Epub 2021 Sep 25. PMID: 34562313; PMCID: PMC8528465.

    Article  CAS  PubMed  Google Scholar 

  51. Algorta J, Pena MA, Alvarez A, Maraschiello C, Maruhn D, Windisch M, Mucke HA. Randomized, crossover, single-blind, placebo-controlled, human pharmacology clinical trial with desoxypeganine, a new cholinesterase and selective MAO-A inhibitor: multiple-dose pharmacokinetics. Int J Clin Pharmacol Ther. 2009;47(7):483–90 (PMID: 19640356).

    CAS  PubMed  Google Scholar 

  52. Furtado-Alle L, Tureck LV, de Oliveira CS, Hortega JVM, Souza RLR. Butyrylcholinesterase and lipid metabolism: Possible dual role in metabolic disorders. Chem Biol Interact. 2023;383 Epub 2023 Aug 25. PMID: 37634560.

    Article  CAS  PubMed  Google Scholar 

  53. Pavković E, Simeon V, Reiner E, Sućić M, Lipovac V. Serum paraoxonase and cholinesterase activities in individuals with lipid and glucose metabolism disorders. Chem Biol Interact. 1993;87(1–3):179–82. https://doi.org/10.1016/0009-2797(93)90040-6. (PMID: 8393741).

    Article  PubMed  Google Scholar 

  54. Meng Q, Ma J, Suo L, Pruekprasert N, Chakrapani P, Cooney RN. Galantamine improves glycemic control and diabetic nephropathy in Leprdb/db mice. Sci Rep. 2023;13(1):15544. https://doi.org/10.1038/s41598-023-42665-2.PMID:37731032;PMCID:PMC10511534.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Bandoni RL, BricherChoque PN, Dellê H, de Moraes TL, Porter MHM, da Silva BD, Neves GA, Irigoyen MC, De Angelis K, Pavlov VA, Ulloa L, Consolim-Colombo FM. Cholinergic stimulation with pyridostigmine modulates a heart-spleen axis after acute myocardial infarction in spontaneous hypertensive rats. Sci Rep. 2021;11(1):9563. https://doi.org/10.1038/s41598-021-89104-8.PMID:33953291;PMCID:PMC8099899.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Glueck CJ. Nonpharmacologic and pharmacologic alteration of high-density lipoprotein cholesterol: therapeutic approaches to prevention of atherosclerosis. Am Heart J. 1985;110(5):1107–15. https://doi.org/10.1016/0002-8703(85)90225-x. (PMID: 2865887).

    Article  CAS  PubMed  Google Scholar 

  57. Ghislat G, Rahman T, Ballester PJ. Identification and validation of carbonic anhydrase II as the first target of the anti-inflammatory drug actarit. Biomolecules. 2020;10(11):1570. https://doi.org/10.3390/biom10111570.PMID:33227945;PMCID:PMC7699199.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. McMahon C, Will A, Hu P, Shah GN, Sly WS, Smith OP. Bone marrow transplantation corrects osteopetrosis in the carbonic anhydrase II deficiency syndrome. Blood. 2001;97(7):1947–50. https://doi.org/10.1182/blood.v97.7.1947. (PMID: 11264157).

    Article  CAS  PubMed  Google Scholar 

  59. Mińczuk K, Baranowska-Kuczko M, Krzyżewska A, Schlicker E, Malinowska B. Cross-Talk between the (Endo)Cannabinoid and Renin-Angiotensin Systems: Basic Evidence and Potential Therapeutic Significance. Int J Mol Sci. 2022;23(11):6350. https://doi.org/10.3390/ijms23116350.PMID:35683028;PMCID:PMC9181166.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Latif Z, Garg N. The impact of marijuana on the cardiovascular system: a review of the most common cardiovascular events associated with marijuana use. J Clin Med. 2020;9(6):1925. https://doi.org/10.3390/jcm9061925.PMID:32575540;PMCID:PMC7355963.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Ghosh M, Naderi S. Cannabis and cardiovascular disease. Curr Atheroscler Rep. 2019;21(6):21. https://doi.org/10.1007/s11883-019-0783-9. (PMID: 30980200).

    Article  PubMed  Google Scholar 

  62. Hinden L, Tam J. Do endocannabinoids regulate glucose reabsorption in the kidney? Nephron. 2019;143(1):24–7. https://doi.org/10.1159/000494512. (Epub 2019 Jan 11 PMID: 30636250).

    Article  CAS  PubMed  Google Scholar 

  63. Jourdan T, Park JK, Varga ZV, Pálóczi J, Coffey NJ, Rosenberg AZ, Godlewski G, Cinar R, Mackie K, Pacher P, Kunos G. Cannabinoid-1 receptor deletion in podocytes mitigates both glomerular and tubular dysfunction in a mouse model of diabetic nephropathy. Diabetes Obes Metab. 2018;20(3):698–708. https://doi.org/10.1111/dom.13150. (Epub 2017 Dec 3 PMID: 29106063).

    Article  CAS  PubMed  Google Scholar 

  64. Wen G, An W, Chen J, Maguire EM, Chen Q, Yang F, Pearce SWA, Kyriakides M, Zhang L, Ye S, Nourshargh S, Xiao Q. Genetic and pharmacologic inhibition of the neutrophil elastase inhibits experimental atherosclerosis. J Am Heart Assoc. 2018;7(4): e008187. https://doi.org/10.1161/JAHA.117.008187.PMID:29437605;PMCID:PMC5850208.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Ibarra-Lara L, Sánchez-Aguilar M, Del Valle-Mondragón L, Soria-Castro E, Cervantes-Pérez LG, Pastelín-Hernández G, Sánchez-Mendoza A. Clofibrate improves myocardial ischemia-induced damage through regulation of renin-angiotensin system and favours a pro-vasodilator profile in left ventricle. J Pharmacol Sci. 2020;144(4):218–28. https://doi.org/10.1016/j.jphs.2020.09.005. (Epub 2020 Sep 18 PMID: 33070841).

    Article  CAS  PubMed  Google Scholar 

  66. Luo P, Zhou Y, Chang HH, Zhang J, Seki T, Wang CY, Inscho EW, Wang MH. Glomerular 20-HETE, EETs, and TGF-beta1 in diabetic nephropathy. Am J Physiol Renal Physiol. 2009;296(3):F556-63. https://doi.org/10.1152/ajprenal.90613.2008. Epub 2009 Jan 7. PMID: 19129258; PMCID: PMC2660192.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Han Y, Xiong S, Zhao H, Yang S, Yang M, Zhu X, Jiang N, Xiong X, Gao P, Wei L, Xiao Y, Sun L. Lipophagy deficiency exacerbates ectopic lipid accumulation and tubular cells injury in diabetic nephropathy. Cell Death Dis. 2021;12(11):1031. https://doi.org/10.1038/s41419-021-04326-y.PMID:34718329;PMCID:PMC8557213.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Chen CH, Chen HA, Liao HT, Chen CH. Hydroxychloroquine exposure reduces the risk of cardiovasular disease events in patients with hypertension or diabetes mellitus. Clin Exp Rheumatol. 2023;41(11):2230–8. https://doi.org/10.55563/clinexprheumatol/y7psev. (Epub 2023 May 24 PMID: 37246776).

    Article  PubMed  Google Scholar 

  69. Yan L, Guo N, Cao Y, Zeng S, Wang J, Lv F, Wang Y, Cao X. miRNA-145 inhibits myocardial infarction-induced apoptosis through autophagy via Akt3/mTOR signaling pathway in vitro and in vivo. Int J Mol Med. 2018;42(3):1537–47. https://doi.org/10.3892/ijmm.2018.3748. Epub 2018 Jun 28. PMID: 29956747; PMCID: PMC6089768.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Mishima E, Ito J, Wu Z, Nakamura T, Wahida A, Doll S, Tonnus W, Nepachalovich P, Eggenhofer E, Aldrovandi M, Henkelmann B, Yamada KI, Wanninger J, Zilka O, Sato E, Feederle R, Hass D, Maida A, Mourão ASD, Linkermann A, Geissler EK, Nakagawa K, Abe T, Fedorova M, Proneth B, Pratt DA, Conrad M. A non-canonical vitamin K cycle is a potent ferroptosis suppressor. Nature. 2022;608(7924):778–83. https://doi.org/10.1038/s41586-022-05022-3. Epub 2022 Aug 3. PMID: 35922516; PMCID: PMC9402432.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Halder M, Petsophonsakul P, Akbulut AC, Pavlic A, Bohan F, Anderson E, Maresz K, Kramann R, Schurgers L. Vitamin K: Double bonds beyond coagulation insights into differences between vitamin K1 and K2 in health and disease. Int J Mol Sci. 2019;20(4):896. https://doi.org/10.3390/ijms20040896.PMID:30791399;PMCID:PMC6413124.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We greatly acknowledge all the researchers and subjects for providing public data. Animals are not involved in our present study.

Funding

None.

Author information

Authors and Affiliations

Authors

Contributions

Conception and design: X.Z. and B.L.; data collection and analysis: Y.C., W.X., and Z.H.; manuscript writing: B.L.; manuscript revising: Y.Z. and Y.D. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yi Zhang or Yan Ding.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

This study was conducted based on published studies and the GEO public database and did not require additional ethical approval or informed consent for participation.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, B., Zhao, X., Xie, W. et al. Identification of co-expressed central genes and transcription factors in acute myocardial infarction and diabetic nephropathy. BMC Med Genomics 17, 134 (2024). https://doi.org/10.1186/s12920-024-01906-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12920-024-01906-7

Keywords