A correction has been applied to this article in:
Corrigendum: Exploring the bi-directional relationship and shared genes between depression and stroke via NHANES and bioinformatic analysis
Zhanghuan Yang1†
Qian Zhang
Shifu Li- 1Department of Oncology, Xiangya Cancer Center, Xiangya Hospital, Central South University, Changsha, Hunan, China
- 2Hainan Medical University, Haikou, China
- 3Department of Neurosurgery, Xiangya Hospital, Central South University, Changsha, Hunan, China
- 4National Clinical Research Center for Geriatric Disorders, Central South University, Changsha, Hunan, China
- 5Department of Neurosurgery, The First people’s Hospital of Changde, Changde, China
- 6Department of Neurology, Xiangya Hospital, Central South University, Changsha, Hunan, China
Background: Stroke and depression are the two most common causes of disability worldwide. Growing evidence suggests a bi-directional relationship between stroke and depression, whereas the molecular mechanisms underlying stroke and depression are not well understood. The objectives of this study were to identify hub genes and biological pathways related to the pathogenesis of ischemic stroke (IS) and major depressive disorder (MDD) and to evaluate the infiltration of immune cells in both disorders.
Methods: Participants from the United States National Health and Nutritional Examination Survey (NHANES) 2005–2018 were included to evaluate the association between stroke and MDD. Two differentially expressed genes (DEGs) sets extracted from GSE98793 and GSE16561 datasets were intersected to generate common DEGs, which were further screened out in cytoHubba to identify hub genes. GO, KEGG, Metascape, GeneMANIA, NetworkAnalyst, and DGIdb were used for functional enrichment, pathway analysis, regulatory network analysis, and candidate drugs analysis. ssGSEA algorithm was used to analyze the immune infiltration.
Results: Among the 29706 participants from NHANES 2005–2018, stroke was significantly associated with MDD (OR = 2.79,95% CI:2.26–3.43, p < 0.0001). A total of 41 common upregulated genes and eight common downregulated genes were finally identified between IS and MDD. Enrichment analysis revealed that the shared genes were mainly involved in immune response and immune-related pathways. A protein-protein interaction (PPI) was constructed, from which ten (CD163, AEG1, IRAK3, S100A12, HP, PGLYRP1, CEACAM8, MPO, LCN2, and DEFA4) were screened. In addition, gene-miRNAs, transcription factor-gene interactions, and protein-drug interactions coregulatory networks with hub genes were also identified. Finally, we observed that the innate immunity was activated while acquired immunity was suppressed in both disorders.
Conclusion: We successfully identified the ten hub shared genes linking the IS and MDD and constructed the regulatory networks for them that could serve as novel targeted therapy for the comorbidities.
Introduction
Stroke is the leading cause of death and disability, leading to significant economic losses as a result of its functional impairments (Meschia et al., 2014). Every year more than 795,000 people in the United States (United States) have a stroke, of which 610,000 are first-time strokes, whereas 185,000 patients have previously had a stroke. The vast majority of stroke cases occur as a result of two specific types of stroke (hemorrhagic and ischemic strokes). United States strokes are dominated by ischemic strokes, which account for 87% of all strokes (Barthels and Das, 2020). The prevalence of depression is growing among the general population, typically characterized by anhedonia and the inability to experience pleasure. A variety of somatic symptoms (psychological disturbance, fatigue, and weight fluctuations) as well as cognitive symptoms (poor concentration and negative cognitions) can accompany depression (Bucciarelli et al., 2020). Depression is sub-categorized into major depressive disorder (MDD) and dysthymia. Epidemiologic data indicates that the average lifetime and 12-month prevalence estimates of MDD are 14.6% and 5.5% in high-income and 11.1% and 5.9% in the low-to middle-income countries (Bromet et al., 2011).
Stroke patients are more likely to suffer from depression than the general population. Growing evidence suggests a bi-directional relationship between stroke and depression: 1) depression is associated with an increased risk of stroke (Pan et al., 2011; Wassertheil-Smoller et al., 2018; Harshfield et al., 2020), and 2) depression is particularly prevalent among stroke survivors (Robinson and Jorge, 2016; Das and Rajanikant, 2018). The prevalence of post-stroke depression (PSD) is estimated to be 29% at any time point up to 5 years following a stroke (Hackett and Pickles, 2014). However, the mechanisms underlying the association between depression and stroke are poorly investigated. Biological factors such as alterations in ascending monoamine systems, neuroplasticity, and glutamate neurotransmission and an increasing of pro-inflammatory cytokines were proposed to explain the mechanisms of PSD (Robinson and Jorge, 2016). Moreover, multiple mechanisms may play roles in depression contributing to stroke. First, Smoking, obesity (Ho et al., 2008), and poor health behaviors (i.e., poor diet, physical inactivity, and smoking) (Strine et al., 2008) may increase stroke risk in depression patients. Second, other major comorbidities, such as diabetes (Wesołowska et al., 2018), atherosclerosis (Joynt et al., 2003), and hypertension (Patten et al., 2009), accompanied by depression, are major risk factors for stroke. Finally, the use of antidepressant medication may potentially contribute to the occurrence of stroke events.
Apart from the above-mentioned mechanisms, genetic factors are likely involved in the pathogenesis of depression and stroke. Increased risk for depression in first-degree relatives of depression probands was observed with an estimated odds ratio of 2.84 from a meta-analysis of the highest-quality family studies (Sullivan et al., 2000). The heritability of MDD has been found to be greater in women (42%) than in men (29%) in a Swedish national twin study (Kendler et al., 2006). There are multiple risk factor genes that were thought to participate in the pathogenesis of depression with extremely complex, polygenic, and epistatic inheritance patterns (Zhao et al., 2019a). There is significant evidence that stroke has a hereditary component based on studies of twins, siblings, and families (Humphries and Morgan, 2004). Heritability for all IS is estimated to be 37.9% (Bevan et al., 2012). The heritability of stroke subtypes varies markedly, with 40.3% for large vessels and 32.6% for cardioembolics but lower for cardioembolic small vessels (16.1%). The genetic involvement in the pathogenesis of both stroke and depression as well as the comorbidity frequency is not yet fully established or whether common overlapping genes and biological mechanisms are subserving both disorders.
A common transcription feature may provide new insights into the pathogenesis of depression and stroke. This study aims to identify hub genes and biological pathways related to the pathogenesis of IS and MDD. Furthermore, as increasing evidence points to the involvement of an immune response in both disorders (Beurel et al., 2020; Iadecola et al., 2020), we evaluate the immune cell infiltration and identify the common immune cells.
Materials and methods
Dataset collection and processing
The data used in the present work was downloaded from the National Health and Nutrition Examination Survey (NHANES) (https://www.cdc.gov/nchs/index.htm) and the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/) based on a microarray or RNA-seq dataset of major depressive disorder (MDD) and ischemic stroke (IS). The NHANES is a research project aimed to assess the health and nutritional status of adults and children in the United States, combining interviews and physical examinations to provide vital and health statistics. The GSE98793 microarray profile included 128 MDD whole blood samples and 64 health samples at the platform of GPL570 Affymetrix U133_Plus2.0 Genechips. The effect of two batches in the GSE98793 dataset were removed by applying removed BatchEffect function of the limma package (Ritchie et al., 2015). The GSE76826 dataset is a microarray profile at the platform of GPL17077 Agilent-039494 SurePrint G3 Human GE v2 8 × 60 K Microarray 039381. The GSE16561 microarray profile contained whole blood from 39 IS patients and was compared with 24 healthy control subjects, measured using a GPL6883 Illumina HumanRef-8 v3.0 expression beadchip. The GSE122709 dataset (including 10 peripheral blood mononuclear cells samples of IS patients and five controls) is a RNA-sequencing dataset and measured at GPL20795 HiSeq X Ten. When multiple probes were matched with one gene, the probe with the highest expression values was annotated in the homologous gene symbol based on the annotation information on the platform.
NHANES
Data of 70190 participants were available in NHANES 2005–2018. Age, sex, race or ethnicity, education level, poverty, marital status, smoking status, stroke, body mass index (BMI), waist circumference and diabetes was included as variables in the analysis. Depression was measured using the Patient Health Questionnaire (PHQ-9). Participants with PHQ-9 total scores≥10 were considered as having MDD. After excluding participants with missing data, 29706 participants were included in our analysis. Continuous variables are presented as the mean (standard deviation), and categorical variables are presented as the frequency (percentage). The chi-square test or Student’s t-test were performed to evaluate the differences between the non-exposure and exposure condition on stroke and MDD. Logistic regression models were performed to calculate odds ratios (ORs) for stroke and MDD.
Identification of differentially expressed genes (DEGs)
After the data standardization and normalization of datasets using the normalizeBetweenArrays function in the “limma” R package, a principal component analysis (PCA) was conducted by using the “factoextra” R package. The DEGs between cases and healthy controls were analyzed by using the “limma” R package. The criteria of p-value <0.05 and |log fold change (FC)|> 0.2 were used to screen the DEGs of MDD and controls, and |log FC| > 0.5 were regarded as cut-off criteria for significant DEGs for IS patients and controls. A volcano plot and a heat map plot were performed by using the R software ggplot2 package (Ginestet, 2011) and “ComplexHeatmap” (Gu et al., 2016) to show significant DEGs, respectively.
Screening of communal DEGs of MDD and IS
After having separately identified the DEGs of MDD and IS, we intersected their DEGs to screen out the communal DEGs that may participate in the pathogenesis of the two diseases. Only when the DEGs had the same expression trends in both diseases were these common genes kept. The processes were conducted and visualized using the “ggVennDiagram” package (Gao et al., 2021). The overlapped genes were further shown in two disorders with a heat map from the perspective of logFC and p-value.
Function enrichment analysis
The “clusterProfiler” package (Yu et al., 2012) was used to enrich the biological processes (BP), cellular components (CC), and molecular function (MF) of Gene Ontology (GO) (Gene Ontology Consortium, 2015) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (Kanehisa and Goto, 2000) of common DEGs.
Protein-protein interaction (PPI) network
To detect potential relationships among the DEG-encoded proteins common to both MDD and IS, a protein-protein interaction (PPI) network was constructed using the Search Tool for the Retrieval of Interacting Genes database (STRING, www.string-db.org) (Szklarczyk et al., 2019). Low confidence of 0.15 was set to find more interactions between proteins. The other parameters were set to the default values (i.e., a full STRING network for nerwork type; evidence for meaning of network edges; and all active interaction sources). The contructed network was imported into Cytoscape to be visualized and further analyzed.
Selection and analysis of hub genes
In this work, we used six common algorithms MCC (Maximal Clique Centrality), MNC (Maximum neighborhood component), DMNC (Density of Maximum Neighborhood Component), Degree, Closeness, and Betweenness) in the cytoHubba plugin of Cytoscape to evaluate and identify hub genes. The detailed information about the six algorithms were descripted in previous article (Chin et al., 2014). The relationships among genes were calculated using Pearson’s correlation methods. The GSE76826 and GSE122709 datasets was applied to validate the expression levels of hub genes.
Subsequently, a co-expression network of these hub genes was constructed via GeneMANIA (http://www.genemania.org/) (Warde-Farley et al., 2010), and their potential functional processes were enriched using the Metascape tool (https://metascape.org/) (Zhou et al., 2019a).
DEG-miRNA interaction analysis
NetworkAnalyst (https://www.networkanalyst.ca/) is an online platform that aimed to provide a wide-range for meta-analyzing gene expression data and constructing gene regulatory networks in a user-friendly manner (Zhou et al., 2019b). The miRTarBase database provided comprehensive information on experimentally validated miRNA-target interactions and was used to identify regulatory miRNAs that influenced DEGs at the post-transcriptional level in NetworkAnalyst.
Recognition of transcription factors
Transcription factors influence target genes at a transcriptional level. Using the Binding and Expression Target Analysis Minus algorithm, ENCODE targeted the transcription factor of genes derived from the ChIP-seq data. We adopted the ENCODE to predict regulatory TFs of our hub genes.
Prediction of potential drugs of hub genes
The Drug–Gene Interaction Database (DGIdb) (http://www.dgidb.org/) is an online database for identifying drug-gene interaction by integrating the data from, for examplethe Drug Target Commons, DrugBank, TTD, PharmGKB, and Chembl and so on (Wagner et al., 2016). The common hub genes were imported into the database to search for potential drugs. The drug-gene interactions were visualized by the “ggalluvial” R package.
Immune infiltration analysis
The enrichment for 28 immune infiltrating cells (Bindea et al., 2013) in the MDD and IS was assessed using a single-sample gene set enrichment analysis (ssGSEA) by using the “GSVA” R package (Hänzelmann et al., 2013). The immune cells with the same enrichment trends for both diseases and significant differences between diseases and the healthy controls were identified as the potential immune cells involved in the pathogenesis. The relationships between hub DEGs and immune cells were also constructed.
Statistical analyses
R software (version R-4.1.0) performed all statistical analyses. The Wilcoxon test was used for statistical analysis between two groups. The relationships of genes with genes and genes with immune cells were constructed by using Pearson’s correlation method. A p-value less than 0.05 (p < 0.05) was considered to indicate statistical significance. The significance level is denoted as follows: *p < 0.05, **p < 0.01, and ***p < 0.001.
Results
Association between stroke and MDD
Baseline characteristics and the results of logistic regression analysis for stroke and MDD were shown in Table 1 and Table 2. The results indicated that MDD was significantly associated with an increased risk of stroke. Compared with non-exposure condition, the odds ratios (ORs) with 95% confidence intervals (CIs) for exposure condition between stroke and MDD was 2.79 (2.26,3.43), p < 0.0001.
Table 1
Table 1. Baseline characteristics and odds ratio of participants by stroke levels in NHANES (2005–2018).
Table 2
Table 2. Baseline characteristics and odds ratio of participants by MDD levels in NHANES (2005–2018).
Identification of DEGs and communal DEGs between MDD and IS
The flow diagram for this study is shown in Figure 1. To uncover the interrelationships of IS with MDD, we first analyzed the human gene expression datasets from the GEO database to identify the dysregulated genes that stimulate MDD and IS separate. PCA results showed that there were two distinctive batches in GSE989793 in Supplementary Figure S1A, and the batch effect was removed in Supplementary Figure S1B. A volcano plot showed that a total of 336 DEGs were identified based on the following criteria: |log2FC|>0.2 and a p-value <0.05, including 194 that were upregulated and 142 that were downregulated between MDD patients and healthy controls in Figure 2A. These deregulated genes are shown with a heat map in Figure 2B. The GSE16561 dataset identified 360 upregulated and 295 downregulated genes taken from IS patient peripheral blood with the cut-off of |log2FC|>0.5 and a p-value <0.05. The deregulated genes were presented with a volcano and heat map plot in Figures 2C, D, respectively. The PCA result for IS patients and controls are shown in Supplementary Figure S2.
Figure 1
Figure 1. Workflow of data analysis in our present work. MDD, major depressive disorder; DEGs, differentially expressed genes; ssGSEA, single sample gene set enrichment analysis; TF, transcription factor; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; PPI, protein-protein interaction.
Figure 2
Figure 2. Landscapes of differentially expressed genes (DEGs) in MDD and IS. A volcano plot (A) and heat map (B) show the DEGs in MDD. A volcano plot (C) and heat map (D) show the DEGs in IS. MDD, major depressive disorder; IS, ischemic stroke.
We further overlapped the deregulated genes of MDD and IS with the same expression trends. The Venn diagram showed that 41 common upregulated genes and eight common downregulated genes were finally identified in Figures 3A, B. The differential expression patterns in the two groups were presented with heat map plots from the perspective of logFC and p-value in Figures 3C, D, respectively.
Figure 3
Figure 3. Identification of shared DEGs between MDD and IS. (A) Overlapping the shared upregulated DEGs. (B) Overlapping the shared downregulated DEGs. (C) A heat map indicates the shared DEGs from the perspective of log(Fold change). (D) A heat map indicates the shared DEGs from the perspective of the p-value. MDD, major depressive disorder; IS, ischemic stroke.
Functional enrichment analysis
The biological processes (BP) results showed that the shared genes were mainly enriched in neutrophil activation involved in immune response, defense response to the bacterium, humoral immune response, and reactive oxygen species metabolic process (Figure 4A). In Figure 4B, we observed that these genes were involved in the vesicle lumen, specific granule, and tertiary granule cell component (CC). The molecular functions (MF) of these shared DEGs were enriched in protein heterodimerization activity and serine-related activity in Figure 4C. The KEGG result showed that the primary immunodeficiency, T cell receptor signaling pathway, and antigen processing and presentation pathways were enriched in Figure 4D.
Figure 4
Figure 4. GO and KEGG enrichment analysis for shared DEGs between MDD and IS. (A–C) BP, CC, and MF of GO analysis were enriched for common DEGs. (D) KEGG pathways of common DEGs. BP, biological processes; CC, cellular components; MF, and molecular functions.
Identification and analysis of hub common DEGs
The shared genes were imported into the STRING database to construct a protein-protein network (Supplementary Figure S3). The network was further visualized in Cytoscape in Figure 5. The red represents the upregulated DEGs, while the turquoise represents the downregulated DEGs. The size of the node shows the degree of interaction with other genes.
Figure 5
Figure 5. PPI network showing the protein interaction for shared DEGs between MDD and IS. The red color represents the commonly upregulated genes. The blue color represents the common downregulated genes. The size of the circle indicates the Degree of the node. The PPI network was generated using STRING and visualized in Cytoscape.
Next, we adopted six algorithms in the cytoHubb plugin of Cytoscape to identify the hub genes. The top 20 genes in each method were visualized in Supplementary Figure S4 and listed in Table 3. We then intersected the 20 genes for each method, and ten overlapped genes (CD163, AEG1, IRAK3, S100A12, HP, PGLYRP1, CEACAM8, MPO, LCN2, and DEFA4) denominated as hub communal DEGs were selected in Figure 6A. The detailed descriptions of the hub genes were listed in Table 4, and that of other DEGs were in Supplementary Table S1. The locations of the 10 genes in the corresponding chromosome are presented in Figure 6B. The violin plot showed that the hub shared genes were significantly expressed in MDD and IS with the same trend in Figures 6C, D. The relationships among genes show that most genes were significantly positively related to each other in MDD (Figure 6E) and stroke (Figure 6F).
Table 3
Table 3. The top20 genes identified by six different methods.
Figure 6
Figure 6. Identification of hub genes from PPI network. (A) A Venn diagram shows ten overlapped genes that were screened out by the six methods from the top 20 genes of six methods in Cytohubba plug in Cytoscape. (B) The gene locations of the ten hub genes. (C) and (D) shows the expression levels of the ten hub genes in MDD and IS, respectively. (E) and (F) depicted the correlations of the ten genes with each other in MDD and IS, respectively. Significance level was denoted by *p-value <.05, **p-value <.01, ***p-value <.001.
Table 4
Table 4. The detailed information and descriptions of hub genes.
The diagnostic ablility of the hub genes in MDD (Figure 7A) and IS (Figure 7B) were visualized with receiver operating characteristic curves. The results shows that HP present the greatest diagnostic value with AUC = 0.671 in MDD, while CD163 display the greatest diagnostic value with AUC = 0.965 in IS.
Figure 7
Figure 7. The diagnostic abilities of the ten hub genes in MDD (A) and IS (B) with ROC curve. ROC, receiver operating characteristic.
We further validated the epression of the ten hub genes in other external datasets (Supplementary Figure S5). However, only S100A12 was validated as the common differentially expressed gene of the two disorders, which need to be verified by in vivo or vitro experiments.
In Figure 8A, the co-expression network of hub genes was constructed using the GeneMANIA website. In the complex PPI network, the interaction of the co-expression accounted for 74.83%, physical interactions for 22.14%, and colocalization for 3.04%. In function analysis, these genes were involved in humoral immune response, secretory granule lumen, cell killing, and regulation of inflammatory response, which was almost consistent with the results from the Metascape analysis (Figure 8B; Supplementary Table S2). In Figure 8C and Supplementary Table S3, we also predicted the potential diseases that the hub genes may be involved in through gene-disease association information collected from the DisGeNET database (Piñero et al., 2020) in Metascape. The results showed that these genes participate in intravascular hemolysis, endotoxemia, and bacterial infections.
Figure 8
Figure 8. The co-expression network and biological functions of hub genes. (A) Hub genes and their co-expression genes were analyzed via GeneMANIA. (B) The biological pathways were enriched for hub genes via Metascape. (C) The potential diseases were participated by these hub genes from the DisGeNET database.
Gene regulatory network analysis of Genes-miRNAs and Genes-TFs
TarBase database was utilized to predict the miRNA of hub genes. All hub shared genes were predicted for their interacted miRNA, and a total of 28 miRNA were determined in Figure 9A. In the gene-miRNA interaction network, LCN2 interacted with the most miRNAs with 13 predicted, followed by HP and IRAK3 with 6 miRNAs. hsa-mir-27a-3p were located in a conspicuous place due to interacting with five hub genes.
Figure 9
Figure 9. Construction of the regulatory network for hub genes. (A) gene-miRNAs interactions. (B) gene-TFs interactions. (C) A Sankey diagram indicating the potential drugs for hub genes. The circle nodes represent the hub genes. The square nodes represent miRNAs. The diamond nodes represent TFs. TF, transcription factor.
We also predicted the experimentally validated TFs of hub genes using ENCODE database in Figure 9B. Only four hub genes predicted their regulatory TFs. A total of 61 TFs were identified, and LCN2 also had the most targeting nodes with 45 TFs, followed by HP with 20 TFs.
The Sankey diagram showed the potential drugs that targeted the hub genes from the DGIdb database (Figure 9C). A total of 30 drugs were predicted, and the detailed information were listed in Supplementary Table S4. Of these, 21 drugs targeted MPO; five drugs targeted S100A12, and 2 drugs each targeted HP and ARG1. No potential drugs could be identified for LCN2, DEFA4, PGLYRP1, CEACAM8, CS163, and IRAK3.
Immune cell infiltration analysis
Using the ssGSEA algorithm, we obtained the immune infiltration of 28 immune cells in the MDD group, IS group, and control group. The immune cells with significant differences between cases and the healthy control group and the same trends were regarded as the potential cells. A total of five immune cells among these 28 types of cells, including activated B cell, activated dendritic cell, effector memory CD8 T cell, macrophage, and natural killer cell were identified, among which active B cell and effector memory CD8 T cell were downregulated, whereas other cells were upregulated, implying the innate immunity was activated while acquired immunity was suppressed in the two diseases (Figures 10A, B). Figures 10C, D also showed the strong relationships between the hub gene and immune cells.
Figure 10
Figure 10. Identification of common immune cells between MDD and IS. (A) The abundance of the immune cell in MDD using the ssGSEA method. (B) The abundance of the immune cell in IS using the ssGSEA method. The immune cells with red color indicate the significantly common immune cells. (C) and (D) A heat map visualized the correlations between common immune cells and shared hub genes. Significance level was denoted by *p-value <.05, **p-value <.01, ***p-value <.001.
Discussion
Depression is a global health problem with a high prevalence and the third leading cause of disability globally (Park and Zarate, 2019). The incidence of suicide associated with depression has been increasing and is the 10th leading cause of death in the United States Similarly, stroke, a neurological disorder characterized by blockage of blood vessels, is a major cause of death and disability worldwide (Johnston et al., 2009). Early studies indicated that depression increased the risk of stroke. A prospective longitudinal study showed that a history of depression was associated with an increased risk of stroke by over twofold (Jackson and Mishra, 2013). Compared with participants with stable low/no depressive symptoms, the participants with the stable high and remitted depressive symptoms had a 2.14 and 1.66 elevated hazard risk of stroke, respectively (Gilsanz et al., 2015). In addtion, post-stroke depression (PSD) is one of the common and serious sequelae of stroke. Folstein et al. (1977) first demonstrated that mood disorder is a more specific complication of stroke. Disability, anxiety, stroke severity, depression pre-stroke, and cognitive impairment all play an important role in PSD, according to a meta-analysis (Ayerbe et al., 2013). Although the bi-directional relationship between stroke and depression is recognized, the underlying mechanism remains a provocative and unresolved question. Considering that stroke and depression have genetic roots, as well as their frequent comorbidity, we speculate shared genes and biological pathways for both stroke and depression.
With the rapid development of sequencing technology, for the first time, we explored the shared gene signatures and molecular mechanisms between MDD and IS from transcriptome data. In our study, we observe a total of 41 genes are simultaneously upregulated, and eight genes are downregulated in both MDD and IS. Biological enrichment analysis shows that these common genes are involved in the immune response, cell killing, and defense response to the bacterium. Moreover, the T cell receptor signaling pathway, primary immunodeficiency, malaria, IL-17 signaling pathway, and rap1 signaling pathway are enriched in the KEGG pathway analysis. Ten overlapped genes (CD163, AEG1, IRAK3, S100A12, HP, PGLYRP1, CEACAM8, MPO, LCN2, and DEFA4) denominated as hub communal DEGs are identified. We observe that the ten hub genes participate in the immune response and cell killing processes. The gene-diseases analysis reveals that intravascular hemolysis, endotoxemia, and bacterial infections are correlated with these genes. Furthermore, we construct gene regulatory networks with gene-miRNA, gene-TF, and gene-drugs, which further provide targets for therapeutic interventions. Finally, we depict the immune landscapes for both disorders and found that five immune cells, including activated B cell, activated dendritic cell, effector memory CD8 T cell, macrophage, and natural killer cell, were significantly different in both diseases. Further analysis indicates innate immunity may be activated whereas acquired immunity may be suppressed.
There is consistently and robust evidence supporting the role of inflammation in depression. The inflammatory response in MDD patients was characterized by increased production of complement, chemoattractors, and pro-inflammatory cytokines in peripheral blood cerebrospinal fluid, and post-mortem brain samples (Miller and Raison, 2016). Cytokines, which mediate the innate immune response, including IL-1, tumor necrosis factor (TNF)-alpha, C-reactive protein (CRP), and IL-6, from peripheral blood are considered the most reliable biomarkers of inflammation in patients with depression (Miller et al., 2009). In addition, by inhibiting pro-inflammatory cytokines or their signaling pathways, depressed mood can be improved and conventional antidepressants better tolerated (Kenis and Maes, 2002; Bluthé et al., 2006). Furthermore, by producing anti-inflammatory cytokines (IL-2, IL-4, and IL-10) and/or activating T regulatory (Treg) cells, the effects of immune response were also counter-balanced or compensated in MDD patients (Dowlati et al., 2010). Growing evidence also revealed an intimate relationship between the immune system and all stages of the ischemic cascade, from the acute intravascular events induced by a blockage of the blood supply to the parenchymal process causing brain damage (Iadecola and Anrather, 2011; Endres et al., 2022). A recent review summarized that pro-inflammatory interleukins (IL-1b, IL-6, IL-8, IL-12, IL-15, IL-16, IL-20, IL-18, and IL-23/IL-17) and anti-inflammatory interleukins (IL-2, IL-4, IL-10, IL-13, IL-19, and IL-33) were involved in the pathogenesis of IS (Zhu et al., 2022). As inflammation is common after stroke and depression, immunological processes were proposed as the underlying mechanism triggering PSD (Pascoe et al., 2011). Inflammatory markers such as CRP, ferritin, and neopterin have been linked to PSD development later in life (Becker, 2016). Our GO and KEGG analysis demonstrates that the immune response is enriched in the common DEGs and the ten hub genes of MDD and IS. Moreover, in our immune infiltration analysis, the abundance of active B cells and effector memory CD8 T cell decreases, while that of activated dendritic cell, macrophage, and natural killer cell increases in both disorders, which may provide new insight into the common pathogenesis and immunotherapy for both diseases. Specifically, the IL-17 signaling pathway was also observed in our biological function annotation of shared genes. The IL-17 family is an evolutionarily old cytokine family consisting of six members (IL-17A-F), dominantly produced by immune cells of the adaptive and innate lymphocyte lineages, including CD4+ Th17 cells, CD8+ Tc17 cells, γδT17 cells, MAIT cells, and innate lymphoid cells ILC3 (Majumder and McGeachy, 2021). It has been observed that IL-17 levels are high in the central nervous system (CNS) during inflammatory responses, including IS and MDD. Peripheral blood samples from patients with IS show an increased expression of IL-17 compared with healthy individuals (Kostulas et al., 1999). High plasma levels of IL-17 were also detected in MDD patients (Waisman et al., 2015). Combined with our results, we speculate that the IL-17 signaling pathway plays an important role in the shared mechanisms of MDD and IS.
A PPI network analysis was conducted among the proteins derived from shared DEGs to depict functional and physical interactions between IS and MDD. By integrating with six algorithms (MCC, MNC, DMNC, Degree, Closeness, and Betweenness) in the cytoHubba plugin of Cytoscape, we identify ten hub communal DEGs (CD163, AEG1, IRAK3, S100A12, HP, PGLYRP1, CEACAM8, MPO, LCN2, and DEFA4), which may serve as potential interventional targets.
CD163, the hemoglobin scavenger receptor, is a macrophage-specific protein of the “alternative activation” phenotype and played a major role in dampening the inflammatory response (Moestrup and Møller, 2004). The upregulation of CD163 in monocytes was observed in MDD patients compared with healthy controls (Simon et al., 2021). However, in MDD in vitro experiments, sub-anesthetic doses of ketamine, an antidepressant (Murrough et al., 2013), program human monocytes into M2c-like macrophages (anti-inflammatory phenotype) by inducing high levels of CD163 and MERTK (Nowak et al., 2019). Compared with the CD14+ classical subtype, CD163 expression was more pronounced in CD16+ non-classical and intermediate monocytes after IS and may serve as a potential biomarker of monocyte activation (Greco et al., 2021). Moreover, the percentage of CD163+/CD16+ events 24 h after IS was positively associated with stroke severity and disability. In our analysis, in comparison with the controls, the higher expression of CD163 in both MDD and IS is observed, which may act as a shared risk gene for IS and MDD.
The interleukin receptor-associated kinase (IRAK) family [including IRAK-1, IRAK-2, IRAK-M (IRAK-3), and IRAK-4] are involved in regulating Toll-like receptor (TLR) and interleukin-1 (IL-1) signaling pathways. Interleukin one receptor-associated kinase 3 (IRAK3) is a protein of 596 amino acids with a molecular mass of 68 kDa and is limited to monocytes and macrophages (Wesche et al., 1999). A recent genome-wide association study (GWAS) identified a genome-wide significant locus (rs11465988) in IRAK3 for esketamine efficacy of anti-depression (i.e., percentage change in symptom severity score compared with baseline). The potential roles of IRAK3 in IS have also been discovered recently. The expression levels of IRAK3 that may link natural killer cells to apoptosis were upregulated in IS through bioinformatics analysis (Feng et al., 2022). In experimental stroke mice, IRAK3 has neuroprotective effects, and its deletion can exacerbate neurovascular damages (Lyu et al., 2018). However, our results identify the enhanced expression in both IS and MDD.
Myeloperoxidase (MPO) is a member of the superfamily of heme peroxidases, that is, mainly found in neutrophils and monocytes. High levels of MPO have been detected in the serum of depressive patients in a twin study (Vaccarino et al., 2008). Inhibiting MPO activity and serotonin reuptake may be a potential new approach to MDD with inflammatory syndrome (Soubhye et al., 2014). Moreover, a significant increase in MPO mRNA expression was observed in peripheral blood cells from patients with recurrent depressive disorder (rDD) compared to controls (Gałecki et al., 2012; Talarowska et al., 2015). The expression of MPO was also associated with the risk of IS (Wright et al., 2009). Concentrations of serum MPO are increased after IS and was associated with stroke severity (Palm et al., 2018; Orion et al., 2020). Inhibiting MPO activity increased cell proliferation and improved neurogenesis after IS (Kim et al., 2016; Kim et al., 2019). Although MPO contributes to both IS and MDD, few studies were conducted to explore the potential mechanism of IS complicated with MDD. Our result may provide a bridge linking the two disorders.
Lipocalin-2 (LCN2) is a member of the highly heterogeneous lipocalin family of secretory proteins. The roles of LCN2 in IS and depression have been proved recently (Zhao et al., 2019b; Vichaya et al., 2019). A study demonstrated that the relationship between LCN2 and the process of PSD may be mediated via the P38 MAPK pathway (Wei et al., 2021). Our study provides potential association for LCN2 and comorbidity between depression and IS.
Haptoglobin encoded by HP participates in the process of depression and stroke from the perspective of genetic and proteomic levels (Maes et al., 1993; Kiga et al., 2008; Ijäs et al., 2013). Considering few studies focusing on the intermediate role of haptoglobin between depression and stroke, this study provides new insight and reference for investigating its potential roles in the comorbidity, such as PSD.
Regulatory biomolecules may serve as potential interventional targets in multiple complex illnesses. TFs play a key role in regulating the ratio of transcription, and miRNAs handle gene regulation and RNA silencing at the post-transcription level. Given the crucial roles of the ten hub common genes, we also analyze the TFs–gene, miRNAs-gene, and drugs-gene interaction to find transcriptional, post-transcriptional, and therapeutic regulators. mir-27a-3p, mir-146a-5p, mir-335-5p, and let-7b-5p are identified to be interacting with at least three hub genes. Furthermore, we discovered that TFs (such as ZNF197, KLF9, KLF11, RERE, ARID4B, TFE3, and FOSL1) target LCN2 and HP simultaneously. Finally, 31 candidate drugs were predicted, among which 21 drugs target MPO. Combined with the above-mentioned roles of MPO in both depression and stroke, these drugs may serve as potential theraputics to treat the comorbidities.
Some limitations should be noted in our work. First, although the gene expression profiling from both diseases are derived from the same tissues, there is inadequate information regarding the blood sample collection time for the studies. The disease course of depression and stroke are different. Second, in this study, all the results were acquired by bioinformatic analysis, and we have not conducted any in vivo or in vitro experiments to verify the different expression levels. Hence, the findings should be interpreted with caution.
Conclusion
We performed a bioinformatic analysis to identify overlapping DEGs subserving both MDD and IS. The communal DEGs participate in the immune response and cell killing processes. Furthermore, ten hub DEGs (CD163, AEG1, IRAK3, S100A12, HP, PGLYRP1, CEACAM8, MPO, LCN2, and DEFA4) were screened out based on six algorithms (MCC, MNC, DMNC, Degree, Closeness, and Betweenness). Immune infiltration analysis shows that the innate immunity was activated whereas acquired immunity was suppressed in both diseases. These findings increase our understating of the association of IS with depression at a transcriptional level. The final gene regulatory network may shed light on novel therapeutic targets for both disorders.
Data availability statement
Publicly available datasets were analyzed in this study. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
Conceptualization, ZY, SL, and MH; methodology, ZY, SL, and MH; software, ZY and MH; validation, QZ; formal analysis, QZ; investigation, QZ; writing-original draft preparation, ZY, SL, and MH; writing-review and editing, HC; visualization, DL; supervision, DL.
Funding
This work was supported by the National Natural Science Foundation of China (82101947), the China Postdoctoral Science Foundation (2021TQ0371 and 2021M703636), and the Fundamental Research Funds for the Central Universities of Central South University (2020zzts269).
Acknowledgments
We acknowledge GEO databases for providing their platforms and contributors for uploading their meaningful datasets. Thanks to Jing Zhang (Shanghai Tongren Hospital) for his work in the NHANES database. His outstanding work, nhanesR package and webpage, makes it easier for us to explore the NHANES database.
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.2023.1004457/full#supplementary-material
References
Ayerbe, L., Ayis, S., Wolfe, C. D., and Rudd, A. G. (2013). Natural history, predictors and outcomes of depression after stroke: Systematic review and meta-analysis. Br. J. psychiatry J. Ment. Sci. 202, 14–21. doi:10.1192/bjp.bp.111.107664
Barthels, D., and Das, H. (2020). Current advances in ischemic stroke research and therapies. Biochimica biophysica acta Mol. basis Dis. 1866, 165260. doi:10.1016/j.bbadis.2018.09.012
Becker, K. J. (2016). Inflammation and the silent sequelae of stroke. Neurother. J. Am. Soc. Exp. Neurother. 13, 801–810. doi:10.1007/s13311-016-0451-5
Beurel, E., Toups, M., and Nemeroff, C. B. (2020). The bidirectional relationship of depression and inflammation: Double trouble. Neuron 107, 234–256. doi:10.1016/j.neuron.2020.06.002
Bevan, S., Traylor, M., Adib-Samii, P., Malik, R., Paul, N. L. M., Jackson, C., et al. (2012). Genetic heritability of ischemic stroke and the contribution of previously reported candidate gene and genomewide associations. Stroke 43, 3161–3167. doi:10.1161/strokeaha.112.665760
Bindea, G., Mlecnik, B., Tosolini, M., Kirilovsky, A., Waldner, M., Obenauf, A. C., et al. (2013). Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity 39, 782–795. doi:10.1016/j.immuni.2013.10.003
Bluthé, R. M., Kelley, K. W., and Dantzer, R. (2006). Effects of insulin-like growth factor-I on cytokine-induced sickness behavior in mice. Brain, Behav. Immun. 20, 57–63. doi:10.1016/j.bbi.2005.02.003
Bromet, E., Andrade, L. H., Hwang, I., Sampson, N. A., Alonso, J., de Girolamo, G., et al. (2011). Cross-national epidemiology of DSM-IV major depressive episode. BMC Med. 9, 90. doi:10.1186/1741-7015-9-90
Bucciarelli, V., Caterino, A. L., Bianco, F., Caputi, C. G., Salerni, S., Sciomer, S., et al. (2020). Depression and cardiovascular disease: The deep blue sea of women's heart. Trends Cardiovasc. Med. 30, 170–176. doi:10.1016/j.tcm.2019.05.001
Chin, C. H., Chen, S. H., Wu, H. H., Ho, C. W., Ko, M. T., and Lin, C. Y. (2014). cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst. Biol. 8, S11. doi:10.1186/1752-0509-8-s4-s11
Das, J., and Rajanikant, G. K. (2018). Post stroke depression: The sequelae of cerebral stroke. Neurosci. Biobehav. Rev. 90, 104–114. doi:10.1016/j.neubiorev.2018.04.005
Dowlati, Y., Herrmann, N., Swardfager, W., Liu, H., Sham, L., Reim, E. K., et al. (2010). A meta-analysis of cytokines in major depression. Biol. psychiatry 67, 446–457. doi:10.1016/j.biopsych.2009.09.033
Endres, M., Moro, M. A., Nolte, C. H., Dames, C., Buckwalter, M. S., and Meisel, A. (2022). Immune pathways in etiology, acute phase, and chronic sequelae of ischemic stroke. Circulation Res. 130, 1167–1186. doi:10.1161/circresaha.121.319994
Feng, L., Tian, R., Mu, X., Chen, C., Zhang, Y., Cui, J., et al. (2022). Identification of genes linking natural killer cells to apoptosis in acute myocardial infarction and ischemic stroke. Front. Immunol. 13, 817377. doi:10.3389/fimmu.2022.817377
Folstein, M. F., Maiberger, R., and McHugh, P. R. (1977). Mood disorder as a specific complication of stroke. J. Neurol. Neurosurg. Psychiatry 40, 1018–1020. doi:10.1136/jnnp.40.10.1018
Gałecki, P., Gałecka, E., Maes, M., Chamielec, M., Orzechowska, A., Bobińska, K., et al. (2012). The expression of genes encoding for COX-2, MPO, iNOS, and sPLA2-IIA in patients with recurrent depressive disorder. J. Affect. Disord. 138, 360–366. doi:10.1016/j.jad.2012.01.016
Gao, C. H., Yu, G., and Cai, P. (2021). ggVennDiagram: An intuitive, easy-to-use, and highly customizable R package to generate Venn diagram. Front. Genet. 12, 706907. doi:10.3389/fgene.2021.706907
Gene Ontology Consortium (2015). Gene Ontology Consortium: Going forward. Nucleic acids Res. 43, D1049–D1056. doi:10.1093/nar/gku1179
Gilsanz, P., Walter, S., Tchetgen Tchetgen, E. J., Patton, K. K., Moon, J. R., Capistrant, B. D., et al. (2015). Changes in depressive symptoms and incidence of first stroke among middle-aged and older US adults. J. Am. Heart Assoc. 4, e001923. doi:10.1161/jaha.115.001923
Ginestet, C. (2011). ggplot2: Elegant graphics for data analysis. J. R. Stat. Soc. Ser. A-STATISTICS Soc. 174, 245–246. doi:10.1111/j.1467-985X.2010.00676_9.x
Greco, R., Demartini, C., Zanaboni, A. M., Tumelero, E., Persico, A., Candeloro, E., et al. (2021). CD163 as a potential biomarker of monocyte activation in ischemic stroke patients. Int. J. Mol. Sci. 22, 6712. doi:10.3390/ijms22136712
Gu, Z., Eils, R., and Schlesner, M. (2016). Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinforma. Oxf. Engl. 32, 2847–2849. doi:10.1093/bioinformatics/btw313
Hackett, M. L., and Pickles, K. (2014). Part I: Frequency of depression after stroke: An updated systematic review and meta-analysis of observational studies. Int. J. stroke official J. Int. Stroke Soc. 9, 1017–1025. doi:10.1111/ijs.12357
Hänzelmann, S., Castelo, R., and Guinney, J. (2013). Gsva: Gene set variation analysis for microarray and RNA-seq data. BMC Bioinforma. 14, 7. doi:10.1186/1471-2105-14-7
Harshfield, E. L., Pennells, L., Schwartz, J. E., Willeit, P., Kaptoge, S., Bell, S., et al. (2020). Association between depressive symptoms and incident cardiovascular diseases. Jama 324, 2396–2405. doi:10.1001/jama.2020.23068
Ho, R. C., Niti, M., Kua, E. H., and Ng, T. P. (2008). Body mass index, waist circumference, waist-hip ratio and depressive symptoms in Chinese elderly: A population-based study. Int. J. geriatric psychiatry 23, 401–408. doi:10.1002/gps.1893
Humphries, S. E., and Morgan, L. (2004). Genetic risk factors for stroke and carotid atherosclerosis: Insights into pathophysiology from candidate gene approaches. Lancet Neurology 3, 227–235. doi:10.1016/s1474-4422(04)00708-2
Iadecola, C., and Anrather, J. (2011). The immunology of stroke: From mechanisms to translation. Nat. Med. 17, 796–808. doi:10.1038/nm.2399
Iadecola, C., Buckwalter, M. S., and Anrather, J. (2020). Immune responses to stroke: Mechanisms, modulation, and therapeutic potential. J. Clin. investigation 130, 2777–2788. doi:10.1172/jci135530
Ijäs, P., Saksi, J., Soinne, L., Tuimala, J., Jauhiainen, M., Jula, A., et al. (2013). Haptoglobin 2 allele associates with unstable carotid plaque and major cardiovascular events. Atherosclerosis 230, 228–234. doi:10.1016/j.atherosclerosis.2013.07.008
Jackson, C. A., and Mishra, G. D. (2013). Depression and risk of stroke in midaged women: A prospective longitudinal study. Stroke 44, 1555–1560. doi:10.1161/strokeaha.113.001147
Johnston, S. C., Mendis, S., and Mathers, C. D. (2009). Global variation in stroke burden and mortality: Estimates from monitoring, surveillance, and modelling. Lancet Neurology 8, 345–354. doi:10.1016/s1474-4422(09)70023-7
Joynt, K. E., Whellan, D. J., and O'Connor, C. M. (2003). Depression and cardiovascular disease: Mechanisms of interaction. Biol. psychiatry 54, 248–261. doi:10.1016/s0006-3223(03)00568-7
Kanehisa, M., and Goto, S. (2000). Kegg: Kyoto encyclopedia of genes and genomes. Nucleic acids Res. 28, 27–30. doi:10.1093/nar/28.1.27
Kendler, K. S., Gatz, M., Gardner, C. O., and Pedersen, N. L. (2006). A Swedish national twin study of lifetime major depression. Am. J. psychiatry 163, 109–114. doi:10.1176/appi.ajp.163.1.109
Kenis, G., and Maes, M. (2002). Effects of antidepressants on the production of cytokines. Int. J. Neuropsychopharmacol. 5, 401–412. doi:10.1017/s1461145702003164
Kiga, C., Sakurai, H., Goto, H., Hayashi, K., Shimada, Y., and Saiki, I. (2008). Proteomic identification of haptoglobin as a stroke plasma biomarker in spontaneously hypertensive stroke-prone rats. Life Sci. 83, 625–631. doi:10.1016/j.lfs.2008.08.013
Kim, H., Wei, Y., Lee, J. Y., Wu, Y., Zheng, Y., Moskowitz, M. A., et al. (2016). Myeloperoxidase inhibition increases neurogenesis after ischemic stroke. J. Pharmacol. Exp. Ther. 359, 262–272. doi:10.1124/jpet.116.235127
Kim, H. J., Wei, Y., Wojtkiewicz, G. R., Lee, J. Y., Moskowitz, M. A., and Chen, J. W. (2019). Reducing myeloperoxidase activity decreases inflammation and increases cellular protection in ischemic stroke. J. Cereb. blood flow metabolism official J. Int. Soc. Cereb. Blood Flow Metabolism 39, 1864–1877. doi:10.1177/0271678x18771978
Kostulas, N., Pelidou, S. H., Kivisäkk, P., Kostulas, V., and Link, H. (1999). Increased IL-1beta, IL-8, and IL-17 mRNA expression in blood mononuclear cells observed in a prospective ischemic stroke study. Stroke 30, 2174–2179. doi:10.1161/01.str.30.10.2174
Lyu, C., Zhang, Y., Gu, M., Huang, Y., Liu, G., Wang, C., et al. (2018). IRAK-M deficiency exacerbates ischemic neurovascular injuries in experimental stroke mice. Front. Cell. Neurosci. 12, 504. doi:10.3389/fncel.2018.00504
Maes, M., Scharpé, S., Meltzer, H. Y., and Cosyns, P. (1993). Relationships between increased haptoglobin plasma levels and activation of cell-mediated immunity in depression. Biol. psychiatry 34, 690–701. doi:10.1016/0006-3223(93)90042-c
Majumder, S., and McGeachy, M. J. (2021). IL-17 in the pathogenesis of disease: Good intentions gone awry. Annu. Rev. Immunol. 39, 537–556. doi:10.1146/annurev-immunol-101819-092536
Meschia, J. F., Bushnell, C., Boden-Albala, B., Braun, L. T., Bravata, D. M., Chaturvedi, S., et al. (2014). Guidelines for the primary prevention of stroke: A statement for healthcare professionals from the American heart association/American stroke association. Stroke 45, 3754–3832. doi:10.1161/str.0000000000000046
Miller, A. H., Maletic, V., and Raison, C. L. (2009). Inflammation and its discontents: The role of cytokines in the pathophysiology of major depression. Biol. psychiatry 65, 732–741. doi:10.1016/j.biopsych.2008.11.029
Miller, A. H., and Raison, C. L. (2016). The role of inflammation in depression: From evolutionary imperative to modern treatment target. Nat. Rev. Immunol. 16, 22–34. doi:10.1038/nri.2015.5
Moestrup, S. K., and Møller, H. J. (2004). CD163: A regulated hemoglobin scavenger receptor with a role in the anti-inflammatory response. Ann. Med. 36, 347–354. doi:10.1080/07853890410033171
Murrough, J. W., Iosifescu, D. V., Chang, L. C., Al Jurdi, R. K., Green, C. E., Perez, A. M., et al. (2013). Antidepressant efficacy of ketamine in treatment-resistant major depression: A two-site randomized controlled trial. Am. J. psychiatry 170, 1134–1142. doi:10.1176/appi.ajp.2013.13030392
Nowak, W., Grendas, L. N., Sanmarco, L. M., Estecho, I. G., Arena, A. R., Eberhardt, N., et al. (2019). Pro-inflammatory monocyte profile in patients with major depressive disorder and suicide behaviour and how ketamine induces anti-inflammatory M2 macrophages by NMDAR and mTOR. EBioMedicine 50, 290–305. doi:10.1016/j.ebiom.2019.10.063
Orion, D., von Landenberg, P., Itsekson-Hayosh, Z., Schwammenthal, Y., Tsabari, R., Merzeliak, O., et al. (2020). Plasma myeloperoxidase levels in acute brain ischaemia and high grade carotid stenosis. Eur. J. neurology 27, 1604–1611. doi:10.1111/ene.14279
Palm, F., Pussinen, P. J., Safer, A., Tervahartiala, T., Sorsa, T., Urbanek, C., et al. (2018). Serum matrix metalloproteinase-8, tissue inhibitor of metalloproteinase and myeloperoxidase in ischemic stroke. Atherosclerosis 271, 9–14. doi:10.1016/j.atherosclerosis.2018.02.012
Pan, A., Sun, Q., Okereke, O. I., Rexrode, K. M., and Hu, F. B. (2011). Depression and risk of stroke morbidity and mortality: A meta-analysis and systematic review. Jama 306, 1241–1249. doi:10.1001/jama.2011.1282
Park, L. T., and Zarate, C. A. (2019). Depression in the primary care setting. N. Engl. J. Med. 380, 559–568. doi:10.1056/NEJMcp1712493
Pascoe, M. C., Crewther, S. G., Carey, L. M., and Crewther, D. P. (2011). Inflammation and depression: Why poststroke depression may be the norm and not the exception. Int. J. stroke official J. Int. Stroke Soc. 6, 128–135. doi:10.1111/j.1747-4949.2010.00565.x
Patten, S. B., Williams, J. V., Lavorato, D. H., Campbell, N. R. C., Eliasziw, M., and Campbell, T. S. (2009). Major depression as a risk factor for high blood pressure: Epidemiologic evidence from a national longitudinal study. Psychosom. Med. 71, 273–279. doi:10.1097/PSY.0b013e3181988e5f
Piñero, J., Ramírez-Anguita, J. M., Saüch-Pitarch, J., Ronzano, F., Centeno, E., Sanz, F., et al. (2020). The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic acids Res. 48, D845–D855. doi:10.1093/nar/gkz1021
Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic acids Res. 43, e47. doi:10.1093/nar/gkv007
Robinson, R. G., and Jorge, R. E. (2016). Post-stroke depression: A review. Am. J. psychiatry 173, 221–231. doi:10.1176/appi.ajp.2015.15030363
Simon, M. S., Schiweck, C., Arteaga-Henríquez, G., Poletti, S., Haarman, B. C. M., Dik, W. A., et al. (2021). Monocyte mitochondrial dysfunction, inflammaging, and inflammatory pyroptosis in major depression. Prog. neuro-psychopharmacology Biol. psychiatry 111, 110391. doi:10.1016/j.pnpbp.2021.110391
Soubhye, J., Aldib, I., Prévost, M., Elfving, B., Gelbcke, M., Podrecca, M., et al. (2014). Hybrid molecules inhibiting myeloperoxidase activity and serotonin reuptake: A possible new approach of major depressive disorders with inflammatory syndrome. J. Pharm. Pharmacol. 66, 1122–1132. doi:10.1111/jphp.12236
Strine, T. W., Mokdad, A. H., Dube, S. R., Balluz, L. S., Gonzalez, O., Berry, J. T., et al. (2008). The association of depression and anxiety with obesity and unhealthy behaviors among community-dwelling US adults. General Hosp. psychiatry 30, 127–137. doi:10.1016/j.genhosppsych.2007.12.008
Sullivan, P. F., Neale, M. C., and Kendler, K. S. (2000). Genetic epidemiology of major depression: Review and meta-analysis. Am. J. psychiatry 157, 1552–1562. doi:10.1176/appi.ajp.157.10.1552
Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING v11: Protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic acids Res. 47, D607–D613. doi:10.1093/nar/gky1131
Talarowska, M., Szemraj, J., and Gałecki, P. (2015). Myeloperoxidase gene expression and cognitive functions in depression. Adv. Med. Sci. 60, 1–5. doi:10.1016/j.advms.2014.06.001
Vaccarino, V., Brennan, M. L., Miller, A. H., Bremner, J. D., Ritchie, J. C., Lindau, F., et al. (2008). Association of major depressive disorder with serum myeloperoxidase and other markers of inflammation: A twin study. Biol. psychiatry 64, 476–483. doi:10.1016/j.biopsych.2008.04.023
Vichaya, E. G., Gross, P. S., Estrada, D. J., Cole, S. W., Grossberg, A. J., Evans, S. E., et al. (2019). Lipocalin-2 is dispensable in inflammation-induced sickness and depression-like behavior. Psychopharmacology 236, 2975–2982. doi:10.1007/s00213-019-05190-7
Wagner, A. H., Coffman, A. C., Ainscough, B. J., Spies, N. C., Skidmore, Z. L., Campbell, K. M., et al. (2016). DGIdb 2.0: Mining clinically relevant drug-gene interactions. Nucleic acids Res. 44, D1036–D1044. doi:10.1093/nar/gkv1165
Waisman, A., Hauptmann, J., and Regen, T. (2015). The role of IL-17 in CNS diseases. Acta neuropathol. 129, 625–637. doi:10.1007/s00401-015-1402-7
Warde-Farley, D., Donaldson, S. L., Comes, O., Zuberi, K., Badrawi, R., Chao, P., et al. (2010). The GeneMANIA prediction server: Biological network integration for gene prioritization and predicting gene function. Nucleic acids Res. 38, W214–W220. doi:10.1093/nar/gkq537
Wassertheil-Smoller, S., Qi, Q., Dave, T., Mitchell, B. D., Jackson, R. D., Liu, S., et al. (2018). Polygenic risk for depression increases risk of ischemic stroke: From the stroke genetics network study. Stroke 49, 543–548. doi:10.1161/strokeaha.117.018857
Wei, L., Du, Y., Xie, Y., Yu, X., Chen, H., and Qiu, Y. (2021). Lipocalin-2 regulates hippocampal microglial activation in poststroke depression. Front. aging Neurosci. 13, 798335. doi:10.3389/fnagi.2021.798335
Wesche, H., Gao, X., Li, X., Kirschning, C. J., Stark, G. R., and Cao, Z. (1999). IRAK-M is a novel member of the Pelle/interleukin-1 receptor-associated kinase (IRAK) family. J. Biol. Chem. 274, 19403–19410. doi:10.1074/jbc.274.27.19403
Wesołowska, K., Elovainio, M., Hintsa, T., Jokela, M., Pulkki-Raback, L., Lipsanen, J., et al. (2018). Is the association between depressive symptoms and glucose bidirectional? A population-based study. Health Psychol. official J. Div. Health Psychol. Am. Psychol. Assoc. 37, 603–612. doi:10.1037/hea0000612
Wright, C. B., Moon, Y., Paik, M. C., Brown, T. R., Rabbani, L., Yoshita, M., et al. (2009). Inflammatory biomarkers of vascular risk as correlates of leukoariosis. Stroke 40, 3466–3471. doi:10.1161/strokeaha.109.559567
Yu, G., Wang, L. G., Han, Y., and He, Q. Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters. Omics a J. Integr. Biol. 16, 284–287. doi:10.1089/omi.2011.0118
Zhao, F., Yue, Y., Jiang, H., and Yuan, Y. (2019a). Shared genetic risk factors for depression and stroke. Prog. neuro-psychopharmacology Biol. psychiatry 93, 55–70. doi:10.1016/j.pnpbp.2019.03.003
Zhao, N., Xu, X., Jiang, Y., Gao, J., Wang, F., Xu, X., et al. (2019b). Lipocalin-2 may produce damaging effect after cerebral ischemia by inducing astrocytes classical activation. J. neuroinflammation 16, 168. doi:10.1186/s12974-019-1556-7
Zhou, G., Soufan, O., Ewald, J., Hancock, R. E. W., Basu, N., and Xia, J. (2019). NetworkAnalyst 3.0: A visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic acids Res. 47, W234–W241. doi:10.1093/nar/gkz240
Zhou, Y., Zhou, B., Pache, L., Chang, M., Khodabakhshi, A. H., Tanaseichuk, O., et al. (2019). Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat. Commun. 10, 1523. doi:10.1038/s41467-019-09234-6
Keywords: ischemic stroke, major depressive disorder, bioinformatics, shared genes, immune infiltration, NHANES
Citation: Yang Z, He M, Zhang Q, Li S, Chen H and Liao D (2023) Exploring the bi-directional relationship and shared genes between depression and stroke via NHANES and bioinformatic analysis. Front. Genet. 14:1004457. doi: 10.3389/fgene.2023.1004457
Received: 27 July 2022; Accepted: 13 March 2023;
Published: 31 March 2023.
Edited by:
Thayne Kowalski, Centro Universitário Cesuca, BrazilReviewed by:
Julia do Amaral Gomes, Ministério da Saúde, BrazilGiovanna Giudicelli, Federal University of Rio Grande do Sul, Brazil
Copyright © 2023 Yang, He, Zhang, Li, Chen and Liao. 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: Hua Chen, 936483206@qq.com; Di Liao, liaodi1118@csu.edu.cn
†These authors have contributed equally to this work