key: cord-0019032-czsr3iif authors: Zhou, Hui; Huang, Liujia; Liang, Lucong; Chen, Liechun; Zou, Chun; Li, Zhenhua; Li, Rongjie; Jian, Chongdong; Zou, Donghua title: Identification of an miRNA Regulatory Network and Candidate Markers for Ischemic Stroke Related to Diabetes date: 2021-07-07 journal: Int J Gen Med DOI: 10.2147/ijgm.s319503 sha: 480f1a7f1f0bd47bb00654f06f8ea7d73859d12b doc_id: 19032 cord_uid: czsr3iif PURPOSE: Type 2 diabetes mellitus (T2DM) increases the risk of ischemic stroke and poor prognosis. This study aimed to identify molecular mechanisms that are dysregulated in T2DM-associated ischemic stroke and candidate genes that might serve as biomarkers. METHODS: The top 25% variance genes in the GSE21321 and GSE22255 datasets were analyzed for coexpression. The differentially expressed mRNAs (DEmRs) between patients with T2DM or ischemic stroke and controls were analyzed. Then, the union of overlapping coexpressed genes and overlapping DEmRs was analyzed. The miRNAs differentially expressed in T2DM-associated ischemic stroke were also analyzed. CIBERSORT was used to evaluate the levels of infiltration by immune cells in T2DM-associated stroke. RESULTS: Thirteen coexpression modules were identified in T2DM and 10 in ischemic stroke, and 594 module genes were shared between the two conditions. A total of 4452 mRNAs differentially expressed between T2DM patients and controls were identified, as were 2390 mRNAs differentially expressed between ischemic stroke and controls. The 771 union genes were enriched mainly in immune-related biological functions and signaling pathways. UBE2N, TGFB3, EXOSC1, and VIM were identified as candidate markers. In addition, we identified miR-576-3p as having the most regulatory roles in both T2DM and ischemic stroke. Mast cell activation was significantly down-regulated in T2DM but up-regulated in ischemic stroke. CONCLUSION: These findings provide numerous testable hypotheses about the pathways underlying T2DM-associated ischemic stroke, which may help identify therapeutic targets. Diabetes mellitus, the most common form of which is type 2 diabetes (T2DM), is a major public health problem. It is estimated that the number of people in the world with the condition will reach 642 million by 2040. 1, 2 The rising prevalence of diabetes mellitus means rising prevalence of numerous health complications such as retinopathy, nephropathy, neuropathy, cognitive impairment and ischemic stroke. 3 Ischemic stroke is the leading cause of death and disability worldwide, affecting about 30 million people. 4 It is the second most frequent cause of death in the world, after heart disease. 5 T2DM increases the risk of ischemic stroke by nearly 2 to 3 times, thereby increasing the risk of associated mortality. [6] [7] [8] Ischemic stroke can cause more serious micro-and macrovascular damage when it occurs in a diabetes background, profoundly impacting multiple organs and aggravating pathological cascades after stroke. 9 Diabetes can aggravate the acute inflammatory response after ischemic stroke and increase levels of inflammatory factors in the brain. 10 Ischemic stroke can induce or aggravate cerebrovascular damage, especially in diabetic patients under 65 years old. 11 Antihypertensive, antiplatelet, and hypolipidemic agents can reduce the risk of ischemic stroke in diabetic patients, 12 but even with treatment, individuals with diabetes are at higher risk of cerebrovascular accidents than individuals without diabetes. 13 In fact, up to 80% of individuals with diabetes eventually die of macrovascular complications. 14 The clinical symptoms of ischemic stroke appear late in the course of the disease, and given the limited therapeutic options, effective preventive treatments and early diagnostic markers are urgently needed. 15 In addition, more effective risk stratification may lead to better management and even prevention of ischemic stroke. 16 Identifying molecular pathways involved in ischemic stroke means examining gene expression, reflected in levels of mRNAs, but also levels of microRNAs (miRNAs). These regulatory molecules interact with the 3ʹ untranslated region (3ʹUTR) of the target mRNA to inhibit translation. 17 Various stages of cerebral ischemic injury and various types of brain injury involve miRNAs. 18, 19 In fact, ischemic stroke has been linked to dysregulation of miRNAs, which have been proposed as potential therapeutic strategies. 10, 20 The present study aimed to explore molecular mechanisms and potential markers of ischemic stroke in patients with T2DM, based on analysis of differential mRNA and miRNA expression profiles in patients and controls. All data were obtained from the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/ geo). The GSE21321 dataset included mRNA expression profiles derived from blood samples of nine men with T2DM and eight men with no history of T2DM based on the GPL6883 platform. The GSE21321 dataset included miRNA expression profiles derived from blood samples of nine men with T2DM and 10 healthy men with no history of T2DM, based on the GPL10322 platform. The GSE22255 dataset included mRNA expression profiles derived from peripheral blood mononuclear cells from 20 ischemic stroke patients and 20 sex-and age-matched controls, based on the GPL570 platform. The GSE110993 dataset included miRNA expression profiles derived from peripheral blood samples from 20 ischemic stroke patients and 20 matched healthy control subjects according to age and sex, based on the GPL15456 platform. Weighted gene co-expression network analysis (WGCNA), a method that identifies gene coexpression networks based on topological overlap, 20 was performed on the top 25% of genes explaining the observed expression differences in the GPL6883 and GPL570 platforms. Coexpression network modules were constructed using the WGCNA package in the R suite, 21 as follows. Pairwise correlations between genes were used to generate a similarity matrix, then soft threshold power values were calculated to generate a scalefree network topology. The topological overlap matrix (TOM) similarity function 20 was used to convert adjacency values into a TOM matrix, which was used to cluster genes into different modules. The limma package inR 22 was used to identify mRNAs and miRNAs differentially expressed between T2DM patients and controls in the GSE21321 dataset, as well as mRNAs differentially expressed between ischemic stroke patients and controls in the GSE22255 dataset. The DESeq2 package inR 23 was used to identify miRNAs differentially expressed between ischemic stroke patients and controls in the GSE110993 dataset. Differences associated with P < 0.05 were considered significant and included in subsequent analyses. The clusterProfiler package inR 24 was used to examine functional enrichment of differentially expressed mRNAs based on Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). The GO terms included biological processes (BPs), cellular components (CCs), and molecular functions (MFs). Enrichment results were visualized using the ggplot2 https://doi.org/10.2147/IJGM.S319503 International Journal of General Medicine 2021:14 package in R. 25 Gene set enrichment analysis (GSEA) was carried out for mRNAs differentially expressed in both T2DM and ischemic stroke using the fgsea package in R. 26 Enrichment was considered significant if it was associated with P < 0.05. CIBERSORT (https://cibersort.stanford.edu/) was used to assess the levels of infiltration by 22 types of immune cells in T2DM and ischemic stroke. Immune cells expressed as 0 were excluded from the analysis. The limma package in R was used to calculate differences in infiltration levels between T2DM patients and controls, as well as between ischemic stroke patients and controls. We also explored potential correlations between hub gene expression and immune cell infiltration using Pearson correlation analysis. Results associated with P < 0.05 were considered significant. The target genes regulated by differentially expressed miRNAs were predicted using Targetscan (http://www.tar getscan.org/vert_72/). Targets were defined as those with a total context++ score greater than 0.1. These target genes were assessed for their ability to diagnose ischemic stroke in T2DM patients based on the area under the receiver operating characteristic curve (AUC), which was calculated using the pROC package in R. 27 To explore gene expression relationships between T2DM and ischemic stroke, we performed WGCNA ( Figure 1 ). In the case of T2DM, 4658 genes explaining the top 25% in observed expression variance were assigned to 13 coexpression modules based on a soft power threshold β = 20 ( Figure 2A ) and the TOM matrix ( Figure 2B ). In the case of ischemic stroke, 5116 genes were assigned to 10 coexpression modules based on a soft power threshold of β = 20 ( Figures 2C and D) . A total of 594 genes were found to overlap between the coexpression modules in T2DM and ischemic stroke, and these genes were considered to be associated with both conditions. A total of 4452 mRNAs were differentially expressed between T2DM patients and controls ( Figure 3A) , of which 1878 were up-regulated and 2574 down-regulated in the disease. A total of 2390 mRNAs were differentially expressed between ischemic stroke patients and controls ( Figure 3B ), of which 1408 were up-regulated and 982 down-regulated in the condition. We identified 70 genes that were up-regulated in both T2DM and ischemic stroke, and 107 that were down-regulated in both conditions ( Figure 3C) . We defined these genes as potentially associated with T2DM and ischemic stroke. To begin to elucidate biological processes and signaling pathways that might be associated with both T2DM and ischemic stroke, we performed enrichment analysis on the union genes which were common across coexpression modules in the two conditions, as well as on the genes differentially expressed in the two conditions. These genes were enriched for the following GO BPs: neutrophil activation, regulation of leukocyte activation, and T cell activation ( Figure 4A ). The genes were enriched for the following CCs: adherens junction, receptor complex, and membrane region ( Figure 4B ). They were enriched for the following MFs: protein heterodimerization, cell adhesion molecule binding, and cytokine receptor binding ( Figure 4C ). The genes were significantly enriched in the following KEGG signaling pathways: cytokine-cytokine receptor interaction, hematopoietic cell lineage, and Th17 cell differentiation ( Figure 4D ). GSEA showed that genes differentially expressed in T2DM were involved in neutrophil extracellular trap formation and bacterial invasion of epithelial cells ( Figure 4E ). Genes differentially expressed in ischemic stroke were involved in cytokine-cytokine receptor interaction and MAPK signaling ( Figure 4F ). We identified 29 miRNAs that were differentially expressed between T2DM patients and controls and 4446 differentially expressed between ischemic stroke patients and controls. The two sets of miRNAs shared four miRNAs: hsa-miR-299-3p, hsa-miR-320a, hsa-miR-576-3p, and hsa-miR-665. These four miRNAs were predicted to target the mRNAs of 134 union genes. These genes were significantly enriched in 83 KEGG pathways, primarily apoptosis, FoxO signaling, and cytokine-cytokine receptor interaction ( Figure S1 ). The following seven genes in the GSE21321 and GSE22255 datasets were involved in the four signaling pathways mentioned above, and their expression levels were able to predict T2DM and ischemic stroke with an AUC > 0.75 ( Figure 5A ): UBE2N, TGFB3, EXOSC1, VIM, PTGS2, IL10RB, and CXCL3. In addition, the expression of the first four genes was altered in the same direction in both T2DM and ischemic stroke, so they were considered candidate markers ( Figure 5B ). Based on these results, we generated a regulatory network of miRNAs in the two conditions ( Figure 5C ). In this network, miR-576-3p emerged as regulating the most target genes and affecting the most KEGG pathways. The enrichment analysis suggested that genes differentially expressed in T2DM and ischemic stroke were enriched in immune-related functions. Therefore, we compared levels of immune cells in T2DM patients ( Figure 6A ) and ischemic stroke patients ( Figure 6B) . T2DM patients showed a higher proportion of resting mast cells, while ischemic stroke patients showed a higher proportion of monocytes. Compared to controls, T2DM patients showed significantly higher levels of neutrophils, lower levels of CD8+ T cells, and higher levels of activated mast cells ( Figure 6C ). Compared to controls, ischemic stroke patients showed significantly lower levels of resting mast cells ( Figure 6D ). The strongest correlations between differentially expressed genes and immune cell infiltration were the positive correlation of TGS2 and CXCL3 with mast cell activation in T2DM ( Figure 6E) , and the positive correlation of VIM with levels of resting mast cells in ischemic stroke ( Figure 6F ). The rising prevalence of T2DM worldwide means rising risk of ischemic stroke. To enable screening of diabetic patients for stroke risk and to guide appropriate management and treatment strategies, we identified numerous candidate markers and their potential miRNA regulators in T2DM-associated ischemic stroke. The findings from this study will guide future experimental and bioinformatics analyses that may help treat and even prevent ischemic stroke among individuals with diabetes. To maximize the probability of identifying genes associated with T2DM-related ischemic stroke, we searched for genes whose expression was up-or down-regulated in both conditions. Our enrichment analysis suggests that many immune responses may be altered in the two conditions. Weak immune activation is a risk factor for T2DM onset and for T2DM-associated ischemic stroke. 28 In ischemic stroke patients, an excessive number of neutrophils infiltrate ischemic brain tissue, which can lead to systemic inflammation and breakdown of the blood-brain barrier. 29 Acute ischemic stroke patients show an elevated ratio of neutrophils to lymphocytes and worse prognosis. 30 T cells are involved in the late phase of cerebral ischemia, and different T cell subtypes play different roles in ischemic stroke. 31 Several studies have suggested that in T2DM and ischemic stroke, miRNAs regulate the expression of target genes that mediate inflammatory responses, cell proliferation and apoptosis. [35] [36] [37] We identified four miRNAs that may be associated with the occurrence of ischemic stroke among T2DM patients: hsa-miR-299-3p, hsa-miR-320a, hsa-miR-576-3p, and hsa-miR-665. The miR-299-3p target is associated with N-terminal pro brain natural peptide, levels of which can predict in-hospital mortality of patients with acute ischemic stroke. 38, 39 Levels of miR-320 can discriminate between diabetic and non-diabetic patients and are significantly reduced in stroke patients, especially those with good prognosis. 40, 41 The miR-576-3p has been reported to induce interferon production, and gene therapy to restore interferon production may improve prognosis after ischemic stroke. 42, 43 The miR-665 is down- 44, 45 This suggests that the four miRNAs that we identified -especially miR-576-3p, with the most target genes -may be risk factors and markers for the development of ischemic stroke in T2DM patients. Among the target genes of these miRNAs linked to T2DM and ischemic stroke, we identified several that showed particularly high AUCs for predicting disease: UBE2N, TGFB3, PTGS2, IL10RB, EXOSC1, CXCL3, and VIM. In particular, the differential expression of UBE2N, TGFB3, EXOSC1, and VIM was similar in T2DM and ischemic stroke, so we consider them candidate markers. UBE2N, an E2 ubiquitin-conjugating enzyme, is a key enzyme in Parkin-dependent mitophagy. 46 It plays a key role in synaptosomal accumulation of mutant huntingtin and is involved in multiple neurodegenerative disease processes. 47 Transforming growth factor-β (TGF-β) is a multifunctional inflammatory cytokine that is produced by a variety of inflammatory cells, including leukocytes and macrophages. 48 TGF-β and interleukin-10 (IL-10) drive damage to the bloodbrain barrier in ischemic stroke and may have neuroinflammatory effects. 49 Expression of the IL10RB gene has been associated with the risk of ischemic stroke. 50 EXOSC1 has been associated with intracerebral hemorrhage, 51 though it appears never to have been linked to risk of ischemic stroke. CXCL3 is up-regulated in ischemic brain tissue and can predict stroke. 52 Up- DovePress regulation of VIM, which encodes vimentin, can drive injury after ischemic stroke by triggering reactive gliosis and scar formation. 53 The findings from our study require validation in additional, larger samples and in experimental studies. This work should begin to test the ability of miRNAs to This study explored molecular mechanisms that may be dysregulated in ischemic stroke associated with T2DM, and it identified potential diagnostic markers of such stroke. Potential candidate driver genes and regulatory miRNAs were identified, which will help guide further studies into this debilitating complication of T2DM as well as efforts to develop effective therapies. The raw analyses from this study and codes can be obtained from the corresponding author upon reasonable request. The International Journal of General Medicine is an international, peer-reviewed open-access journal that focuses on general and internal medicine, pathogenesis, epidemiology, diagnosis, monitoring and treatment protocols. The journal is characterized by the rapid reporting of reviews, original research and clinical studies across all disease areas. The manuscript management system is completely online and includes a very quick and fair peer-review system, which is all easy to use. Visit http://www.dovepress.com/ testimonials.php to read real quotes from published authors. Post-stroke neovascularization and functional outcomes differ in diabetes depending on severity of injury and sex: potential link to hemorrhagic transformation IDF diabetes atlas: global estimates for the prevalence of diabetes for 2015 and 2040 Diabetes mellitus predicts cognitive impairment in patients with ischemic stroke NLRP3 inflammasome as a potential treatment in ischemic stroke concomitant with diabetes Heart disease and stroke statistics-2018 update: a report from the American heart association The independent effect of type 2 diabetes mellitus on ischemic heart disease, stroke, and death: a population-based study of 13,000 men and women with 20 years of follow-up Diabetic aggravation of stroke and animal models Comprehensive identification of key genes involved in development of diabetes mellitus-related atherogenesis using weighted gene correlation network analysis Cell-based and exosome therapy in diabetic stroke MicroRNA-146b-3p regulates the development and progression of cerebral infarction with diabetes through RAF1/P38MAPK/COX-2 signaling pathway Fuyuan Xingnao decoction promotes angiogenesis through the Rab1/AT1R pathway in diabetes mellitus complicated with cerebral infarction The diagnosis and management of cerebrovascular disease in diabetes Blood pressure and stroke risk among diabetic patients Establishment of clinical prediction model based on the study of risk factors of stroke in patients with type 2 diabetes mellitus Association of diabetic retinopathy with stroke: a systematic review and meta-analysis. Front Neurol Stroke type, etiology, clinical features and prognosis of diabetic patients in Southern China Recent highlights of research on miRNAs as early potential biomarkers for cardiovascular complications of type 2 diabetes mellitus Non-coding RNAs participate in the ischemia-reperfusion injury Association between polymorphisms in microRNAs and ischemic stroke in an Asian population: evidence based on 6,083 cases and 7,248 controls Mesenchymal stem cell-derived exosome miR-542-3p suppresses inflammation and prevents cerebral infarction WGCNA: an R package for weighted correlation network analysis limma powers differential expression analyses for RNA-sequencing and microarray studies Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 clusterProfiler: an R package for comparing biological themes among gene clusters Application of ggplot2 to pharmacometric graphics An algorithm for fast preranked gene set enrichment analysis using cumulative statistic calculation pROC: an open-source package for R and S+ to analyze and compare ROC curves Sex-related disparities in the incidence and outcomes of ischemic stroke among type 2 diabetes patients. A matched-pair analysis using the Spanish national hospital discharge database for years 2016-2018 Harmful neutrophil subsets in patients with ischemic stroke: association with disease severity Association of a high neutrophil-to-lymphocyte ratio with hyperdense artery sign and unfavorable short-term outcomes in patients with acute ischemic stroke Neuroinflammatory mechanisms in ischemic stroke: focus on cardioembolic stroke, background, and therapeutic approaches Age-dependent involvement of gut mast cells and histamine in post-stroke inflammation Local mast cell activation promotes neovascularization Targeting the autonomic nervous system for risk stratification, outcome prediction and neuromodulation in ischemic stroke Type 2 diabetes mellitus and cardiovascular disease: genetic and epigenetic links MicroRNAs as diagnostic and prognostic biomarkers in ischemic stroke-a comprehensive review and bioinformatic analysis. Cells Can miRNAs be considered as diagnostic and therapeutic molecules in ischemic stroke pathogenesis?-Current status The prognostic value of circulating microRNAs in heart failure: preliminary results from a genome-wide expression study N-terminal pro-brain natriuretic peptide levels and short term prognosis in acute ischemic stroke Molecular mechanisms underpinning microparticle-mediated cellular injury in cardiovascular complications associated with diabetes Expression profile of MicroRNAs in young stroke patients The Neat Dance of COVID-19: NEAT1, DANCR, and co-modulated cholinergic RNAs link to inflammation Therapeutic spectrum of interferon-beta in ischemic stroke Serum microRNA profiling and bioinformatics analysis of patients with type 2 diabetes mellitus in a Chinese population miRNAs: roles and clinical applications in vascular disease The ubiquitin-conjugating enzymes UBE2N, UBE2L3 and UBE2D2/3 are essential for Parkin-dependent mitophagy Aged monkey brains reveal the role of ubiquitin-conjugating enzyme UBE2N in the synaptosomal accumulation of mutant huntingtin Platelet expression of transforming growth factor beta 1 is enhanced and associated with cardiovascular prognosis in patients with acute coronary syndrome Neuroinflammatory mechanisms of blood-brain barrier damage in ischemic stroke Association between IL10, IL10RA, and IL10RB SNPs and ischemic stroke with hypertension in Korean population Intracerebral hemorrhage and ischemic stroke of different etiologies have distinct alternatively spliced mRNA profiles in the blood: a pilot RNA-seq study Expression of cytokines and chemokines as predictors of stroke outcomes in acute ischemic stroke Neuroinflammation: friend and foe for ischemic stroke The authors declare that this research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.