Skip to main content

Identification of upregulated NF-κB inhibitor alpha and IRAK3 targeting lncRNA following intracranial aneurysm rupture-induced subarachnoid hemorrhage



This study was performed to identify genes and lncRNAs involved in the pathogenesis of subarachnoid hemorrhage (SAH) from ruptured intracranial aneurysm (RIA).


Microarray GSE36791 was downloaded from Gene Expression Omnibus (GEO) database followed by the identification of significantly different expressed RNAs (DERs, including lncRNA and mRNA) between patients with SAH and healthy individuals. Then, the functional analyses of DEmRNAs were conducted and weighted gene co-expression network analysis (WGCNA) was also performed to extract the modules associated with SAH. Following, the lncRNA-mRNA co-expression network was constructed and the gene set enrichment analysis (GSEA) was performed to screen key RNA biomarkers involved in the pathogenesis of SAH from RIA. We also verified the results in a bigger dataset GSE7337.


Totally, 561 DERs, including 25 DElncRNAs and 536 DEmRNAs, were identified. Functional analysis revealed that the DEmRNAs were mainly associated with immune response-associated GO-BP terms and KEGG pathways. Moreover, there were 6 modules significantly positive-correlated with SAH. The lncRNA-mRNA co-expression network contained 2 lncRNAs (LINC00265 and LINC00937) and 169 mRNAs. The GSEA analysis showed that these two lncRNAs were associated with three pathways (cytokine-cytokine receptor interaction, neurotrophin signaling pathway, and apoptosis). Additionally, IRAK3 and NFKBIA involved in the neurotrophin signaling pathway and apoptosis while IL1R2, IL18RAP and IL18R1 was associated with cytokine-cytokine receptor interaction pathway. The expression levels of these genes have the same trend in GSE36791 and GSE7337.


LINC00265 and LINC00937 may be implicated with the pathogenesis of SAH from RIA. They were involved in three important regulatory pathways. 5 mRNAs played important roles in the three pathways.

Peer Review reports


Subarachnoid hemorrhage (SAH) is an acute, devastating hemorrhagic stroke accounting for 5% of cerebrovascular strokes [1, 2]. SAH from ruptured intracranial aneurysm (RIA) is a non-traumatic type SAH with destructive central nervous system, and has a high in-hospital mortality (45%) [3,4,5], disability rate (30%), and high morbidity of long-term cognitive impairment (50%) among survivors [6, 7]. Delayed cerebral ischemia and vasospasm following SAH are the primary causes of SAH-induced mortality in intensive care unit [8, 9]. Although drug therapy, such as intrathecal nicardipine and statins, has showed some effect on reducing SAH-induced mortality, the overall clinical outcomes has been unsatisfactory due to drug-related myotoxicity and side effect on liver and kidney [10,11,12]. Therefore, it is imperative to gain comprehensive understanding of the pathogenesis of SAH from RIA, which will greatly promote to develop effective therapeutic strategies against this disease.

The transcriptome profiling analysis provides new opportunities for identifying promising therapeutic targets and exploring the molecular mechanisms of various diseases. Notably, increasing researchers have focused on elaborating the underlying pathological mechanisms of SAH by the transcriptomic analysis. For example, Pera et al suggested that SAH from RIA significantly influences the gene expression profiles of peripheral blood cells and 16 transcriptional biomarkers could differentiate IRA patients and healthy individuals by a microarray analysis [13]. Lai et al performed a miRNA microarray analysis and demonstrated that miR-4320 may be a potentially valuable signatures for the detection of SAH [14]. Additionally, Liang et al reported that the expression of a long non-coding RNA (lncRNA) MEG3 was positively correlated with the severity of SAH [15]. They also showed that lncRNA MEG3 inhibited the neuron activity through the PI3K/Akt signaling pathway. Although several RNA biomarkers have been showed to be associated with SAH, the potential mechanism of SAH from RIA have not been fully understood.

We performed a co-expression analysis to identify potential mRNA and lncRNA signatures involved in the pathogenesis of SAH from RIA. The differential expression analysis was carried out and then subjected to a subsequent functional analyses RNAs makers involved in the pathogenesis of SAH from RIA. This study might provide a deeper insight into the molecular mechanisms of SAH from RIA.


Affrymetrix microarray dataset and standard processing

Microarray dataset GSE36791 (including lncRNA and mRNA) was downloaded from the National Center for Biotechnology Information-Gene Expression Omnibus (NCBI-GEO) database ( [13], it was generated by GPL10558 Illumina HumanHT-12 V4.0 expression beadchip platform (, and includes 61 peripheral blood samples (18 healthy controls and 43 patients with SAH from RIA).

Afterwards, Limma package (version 3.34.0; [16] in R 3.4.1 was employed to preprocess raw data. First, we log2 the expression profile data to transform the gene expression data from a skewed distribution to an approximate normal distribution, and then normalize the data with the median standardization method.

Identification of differentially expressed RNAs (DERs) and functional analyses

Firstly, The lncRNAs and mRNAs in GSE36791 datasetes were re-annotated using the HUGO Gene Nomenclature Committee (HGNC) database, consisting of annotated 4055 lncRNAs and 19,198 protein coding genes [17]. Then, the significantly DERs (, including lncRNA and mRNA) between SAH group and normal controls were identified using Limma package (version 3.34.0, in R3.4.1, with the threshold of false discovery rate (FDR) < 0.05 and |logFC| > 0.5. Afterwards, the bidirectional clustering analysis of DEGs was also conducted using R pheatmap package (version 1.0.8; [18].

Finally, the analysis of Gene Ontology biological processes (GO-BP) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways was carried out for those significant DERs using the web-accessible the Database for Annotation, Visualization and Integrated Discovery (DAVID, version 6.8; software according to the cutoff criterion of FDR < 0.05 [19, 20].

Weighted gene co-expression network analysis (WGCNA)

WGCNA is a tool to construct networks and identify gene clusters or modules of gene with co-expression profiling. Here, we carried out a WGCNA analysis to identify the modules associated with SAH using WGCNA integrated algorithm (version 1.61; in R 3.4.1 [21]. We set that each RNA module contains at least 30 RNA elements (cutHeight = 0.995). Then by calculating the Pearson correlation coefficient between the module eigengene (ME) of each module and each sample group, we selected modules that are significantly related to the disease. The selection criteria are as follows: 1) P value was less than 0.05; 2) the correlation coefficient was higher than that in the control (grey) module.

Construction of lncRNA-mRNA regulatory network

The R cor function ( was employed to compute the Pearson correlation coefficient (PCC) between DElncRNAs and DEmRNAs in significant target modules. Thereby, the DElncRNAs-DEmRNAs co-expression network was established and visualized by Cytoscape software (version 3.6.1; [22].

Gene set enrichment analysis (GSEA) for RNAs in regulatory network

To further screen hub genes involved pathogenesis of SAH from RIA, the KEGG enrichment pathway analysis of genes which exhibited co-expression relationships with DElncRNAs was performed by GSEA (, which is a popular tool for interpreting the effects of collective behavior of genes on observed phenotypes by evaluating the deviation of gene expression between disease groups and healthy control group [23]. The cutoff of P < 0.05 was chosen as a statistically significant threshold.

Expression level analysis for hub gens

In the GEO database, we searched a new data set (GSE7337) containing more samples to verify our results. There were 119 peripheral blood samples from SAH patients and 118 controls. The gene expression profiles were determined by Illumina HumanHT-12v4 BeadChips. We extracted the expression levels of hub genes in these two data sets and compared the expression in different samples.


DERs identification and functional analyses

First, we need to standardize the downloaded expression profile data set (Fig. 1). Totally, 1053 lncRNAs and 18,320 protein coding RNAs were identified by annotating with HGNC database after data standardization. Moreover, 561 DERs were identified, including 25 DElncRNAs (10 up-regulated lncRNAs and 15 down-regulated lncRNAs) and 536 DEmRNAs (229 down-regulated mRNAs and 307 up-regulated mRNAs) according to the criteria of FDR < 0.05 and |logFC| ≥ 0.5 (Supplementary Table S1; Fig. 2a and b), with the criteria of FDR < 0.05 and |logFC| ≥ 0.5. Hierarchical clustering showed the distinct expression profiles of these DERs in the SAH from RIA samples and in controls (Fig. 2c). Additionally, functional enrichment analyses of DEmRNAs showed there were 21 GO-BP terms and 7 KEGG significant pathways. Briefly, these genes were mainly associated with GO-BP terms immune response, defense response, inflammatory response, leukocyte activation and regulation. The KEGG analysis indicated that these genes played essential roles in T cell receptor signaling pathway, NOD-like receptor signaling pathway, Cytokine-cytokine receptor interaction as well as Adipocytokine signaling pathway (Fig. 3).

Fig. 1

The boxplots of the data downloaded from GEO database. Gene expression level was analyzed using log2 transformation. Data before normalization

Fig. 2

Statistics and expression of the significantly differentially expressed genes (DEGs) in SAH samples. a the Volcano plot of the SARs in SAH samples. Green lines indicate false discovery rate (FDR) =0.05 (horizontal line) and |log2[Fold change (FC) of FPKM] | = 0.5 (vertical lines), respectively. b the statistics of the DERs and percentage of lncRNAs and mRNAs, as well as down- and up-regulated DERs. c heatmap of the DERs in 43 SAH samples and 18 control (CTRL) samples

Fig. 3

The ingenuity pathways and GO categories associated with significantly differentially expressed mRNAs (DEmRNAs) in SAH. 7 KEGG pathways and 21 GO biological processes were identified to be associated with 536 significantly DEmRNAs. Black to blue note p value close to 0 and 1.0 respectively. Bubbles indicate –log (p value). The larger the bubble, the higher –log (p value)

Construction and analysis of WGCNA modules associated with SAH traits

All the 561 DERs were selected for the WGCNA module analysis. When we choose the parameters of an adjacency function, we only consider those parameter values that cause the network to only approximately satisfy the scale-free topology, such as R ^ 2 > 0.80. First, we need to consider average connectivity. It can determine whether the network can contain enough information such as module detection. Second, the slope of the regression line should be around − 1. In this research, the first parameter value was 0.8. the threshold of soft was power 20 (Fig. 4a). The mean connectivity reached 1 when soft-thresholding power was 20 (Fig. 4b), revealing WGCNA modules were constructed using the approximate scale-free topology. Then, the dissimilarity coefficients between gene nodes were calculated and the clustering tree was built. Accordingly, a total of nine WGCNA modules were constructed in the gene dendrogram (Fig. 4c), with the criteria of eigengenes number ≥ 100 genes, cut height = 0.995. Then we calculated the correlation between modules and status of each sample: disease or not (RIA or CTRL, Supplementary Table S2). The results revealed that six modules significantly positively correlated with disease, including yellow (P-value = 2e-26), blue (P-value = 9e-40), red (P-value = 5e-21), brown (P-value = 2e-31), black (P-value = 2e-31) and pink (P-value = 3e-43) modules. The list of the 217 DERs (2 lncRNAs and 215 mRNAs) included in the six WGCNA modules is shown in Supplementary Table S3.

Fig. 4

The construction of WGCNA modules in dataset GSE36791. a the correlation coefficient square (r2) of eigengenes (red line) under soft threshold power (red number). b the mean connectivity of eigengenes (red line) under soft threshold power (red number). c WGCNA modules in the gene dendrogram based on the criteria of number of genes > 100, cut Height = 0.995, soft threshold power 20 and r2 = 0.8. D, the heatmap showing the correlation of WGCNA modules with SAH clinical traits. Green and red color note negative and positive correlation, respectively. Red stars note 6 modules showed higher, closer correlation (correlation p value ≤0.05 and coefficient > 0.32) than gray module

Construction of lncRNA and mRNA regulatory network

We calculated the PCC between DElncRNA and DEmRNA in the 6 WGCNA modules and retained the connection pairs with PCC higher than 0.6. A total of 240 pairs meet the requirements (Supplementary Table S4). They were composed of 171 nodes. Then the lncRNA and mRNA regulatory network was constructed (Fig. 5). All the notes including LINC00265 and LINC00937 were upregulated in SAH patients from RIA. We found that there were strong relationships between LINC00265 and two mRNAs (nuclear factor kappa B inhibitor alpha, NFKBIA and interleukin 1 receptor associated kinase 3, IRAK3). In addion, LINC00937 closely interacts with NFKBIA.

Fig. 5

The lncRNA and mRNA regulatory network. Box and circle represent lncRNA and mRNA, respectively. Red color notes the upregulation in SAH relative to control

GSEA for the lncRNAs in regulatory network

At last, we performed the GSEA to identify the KEGG pathways associated with the 169 mRNAs genes connected LINC00265 and LINC00937 in regulatory network as mentioned above. The results showed that a total of four KEGG pathways were associated with LINC00265, including cytokine-cytokine receptor interaction (P-value = 0.0042), neurotrophin signaling pathway (P-value = 0.0226), apoptosis (P-value = 0.0269) and neuroactive ligand receptor interaction (P-value = 0.0309). Five KEGG pathways were associated to LINC00937, including apoptosis (P-value = 0.0036), Toll like receptor signaling pathway (P-value = 0.0178), neurotrophin signaling pathway (P-value = 0.0353), cytokine-cytokine receptor interaction (P-value = 0.0080) and MAPK signaling pathway (P-value = 0.0489) (Table 1). Moreover, three overlapping KEGG pathways (‘Cytokine-Cytokine receptor interaction’, ‘Neurotrophin signaling pathway’, and ‘Apoptosis’) were positively correlated with LINC00265 and LINC00937 (Fig. 6). Two genes (NFKBIA and IRAK3) were overlapped in apoptosis and neurotrophin signaling pathway, respectively. Meanwhile, three overlapped genes, including interleukin 1 receptor type 2 (IL1R2), interleukin 18 receptor accessory protein (IL18RAP) and interleukin 18 receptor 1 (IL18R1) were involved in cytokine-cytokine receptor interaction.

Table 1 KEGG pathways associated with LINC00265 and LINC00937 and the target genes
Fig. 6

Overlapping KEGG pathways associated with both LINC00265 and LINC00937. KEGG pathways were identified using GSEA software and the positive correlations of lncRNA expression with pathway activation were confirmed

Expression level analysis for hub gens

We compared the expression level of the seven genes in different samples. As is shown in Fig. 7, the expression trend is consistent in the two data sets.

Fig. 7

Comparison of expression levels of hub genes in different sample groups. Up: GSE36791. Down: GSE73378. White square: The average expression level of genes in the CTRL sample group. Black square: The average expression level of genes in the SHA patients form RIA sample group


SAH caused by RIA is one of the common critical illnesses in cerebrovascular disease. The disability rate among surviving patients is approximately 30% [24]. Bioinformatics was used to analyze the datasets of the RIA and identify potential biomarkers, which is very helpful for clinical diagnosis. Wang et al. identified six hub genes associated with rupture of intracranial aneurysms by weighted gene co-expression network analysis [25]. Wan et al. identified 4 hub genes. Represent potential biomarkers of SAH or impending the likelihood of IA progression and rupture [26]. hsa-miR-1304, hsa-miR-33b, hsa-miR-125b, and hsa-miR-125a-5p were predicted to take part in the pathogenesis of SAH caused by RIA [27]. At present, lncRNA has gradually received attention. Its functions are complex and can participate in the various stages of regulating gene expression.

In previous research, Pera et al. proved that RIA strongly influences the transcriptional profiles of peripheral blood cells [13]. They found 16 genes and lymphocyte-to-monocyte-and-neutrophil gene expression ratios distinguished RIA patients from normal. In our research, we want to explore some hub lncRNAs related to SHA. Five hundred sixty-one DERs, including 25 DElncRNAs and 536 DEmRNAs were identified between SAH group and normal controls and the functional analyses of DEmRNAs showed that they were dramatically enriched in GO-BP terms of immune response and KEGG pathway of T cell receptor signaling pathway. Moreover, the WGCNA module analysis revealed that six modules significantly positively correlated with SAH, which includes 217 DERs (2 lncRNAs and 215 mRNAs). Additionally, lncRNA and mRNA regulatory network contained 240 lncRNA-mRNA interaction pairs among 169 DEmRNAs and 2 DElncRNAs (LINC00265 and LINC00937). The GSEA analysis revealed that LINC00265 and LINC00937 all played essential roles in cytokine-cytokine receptor interaction, neurotrophin signaling pathway, and apoptosis pathways, and these three pathways were found to be implicated with SAH development and suggested five genes (NFKBIA, IRAK3, IL1R2, IL18RAP, and IL18R1) exhibited strong correlation with them. We compared the expression levels of these genes between different samples in data sets GSE36791 and GSE7337. The two data sets show the same expression trend.

A previous study revealed that RIA was preferentially associated with increased cell apoptosis levels and proposed that apoptosis might decrease the resistant ability of aneurysm wall, leading to its rupture [28]. Overwhelming evidence has demonstrated that early brain injury (EBI) largely contributes to elevated mortality risks within 24–72 h after RIA and neuronal apoptosis is responsible for EBI [29, 30]. Yuksel et al. pointed out that apoptotic cell was observed in cortical, subcortical or hippocampal neurons after SAH and cell apoptosis exerted critical roles in EBI, which suggested that apoptosis may be a research agenda for clarifying the mechanisms of SAH from RIA [31]. Herein, we found that LINC00265/LINC00937-NFKBIA was predominantly related to cell apoptosis. Increasing studies have reported that that activation of NFKB pathway is the main cause of EBI occurring following SAH. Conversely, its inhibition had neuroprotective roles against SAH [32, 33]. Interestingly, Zhang et al. argued that inhibiting TGFβ-activated kinase 1 significantly suppressed NFKB activity, reduced neuronal apoptosis in SAH [34]. In neurons or cerebral tissue, NFKB signaling pathways could activate the anti-apoptotic proteins (Bcl-2 and Bcl-xL) [35]. In addition, NFKB also participates in regulating the inflammatory responses. NFKB activation could enhance gene expression level of many inflammatory factors such as IL-1β, which may accelerate the cell death in the brain [36]. You et al. reported that NFKB was activated in the neurons following SAH and regulated inflammatory genes expression, thereby causing delayed brain injury [37]. Our finding showed that a NFKB inhibitor gene, NFKBIA was increased in blood samples of SAH patients compared with healthy controls. The anti-inflammatory effect of NFKBIA in SAH needs to be investigated in the following analysis. Our results also showed that there were strong associations between LINC00265/LINC00937-NFKBIA and neurotrophin signaling pathway. Notably, the dysgregulation of neurotrophin was correlated with neuron survival and death. For example, the binding of neurotrophin to p75 neurotrophin receptor (p75NTR) is necessary for neuron apoptosis [38]. Recently, the interaction between p75 neurotrophin receptor (p75NTR) and the pro-apoptotic BH3-only protein NIX was identified by Shen et al., who argued that the interaction of p75NTR and protein NIX were crucial for the p75NTR-mediated neuron apoptosis in intracerebral hemorrhage [39]. Many studies have also highlighted that neurotrophic factors were crucial for the neurological recovery following SAH [40, 41]. A research showed that oxidized high-mobility group box1 protein dramatically facilitated brain recovery through increasing the expression of neurotrophin in the late stage of SAH, which implied that neurotrophin pathway activation was a key pathological processe in SAH [42]. Therefore, we inferred that LINC00265/LINC00937-NFKBIA might be implicated with the underlying pathogenesis of SAH from RIA via apoptosis and neurotrophin signaling pathway. IRAK3, another NFKB inhibitor was also upregulated in SAH patients compared to healthy controls. The co-expression analysis indicated that LINC00265 strongly interacts with IRAK3. This gene also participated in neurotrophin signaling pathway and apoptosis. Although few studies explored the potential molecular mechanisms with the involvement of IRAK3 in SAH from IAR and there was no direct evidence involving the underlying roles of LINC00265/LINC00937 in SAH, our results provide a new clue that LINC00265-IRAK3 might be involved in the pathogenesis of SAH from RIA by inhibiting NFKB.


This study identified 2 lncRNAs and 5 genes from SAH patients with RIA, which may enhance our current knowledge on this disease and may provide potential biomarkers of this disease. They may be involved in the pathogenesis of SAH from RIA by activating neurotrophin signaling pathway and apoptosis.

Availability of data and materials

The datasets analysed during the current study are available in the National Center for Biotechnology Information-Gene Expression Omnibus (NCBI-GEO) database ( We downloaded two datasets: GSE36791 and GSE7337.



Subarachnoid hemorrhage


Weighted gene co-expression network analysis


NF-κB inhibitor alpha gene


Endothelial nitric oxide synthase


Long non-coding RNA


Gene Ontology


Kyoto Encyclopedia of Genes and Genomes


Gene Set Enrichment Analysis


  1. 1.

    Diringer MN, Dhar R, Scalfani M, Zazulia AR, Chicoine M, Powers WJ, Derdeyn CP. Effect of high-dose simvastatin on cerebral blood flow and static autoregulation in subarachnoid hemorrhage. Neurocrit Care. 2016;25(1):56–63.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  2. 2.

    Garg K, Sinha S, Kale SS, Chandra PS, Suri A, Singh MM, Kumar R, Sharma MS, Pandey RM, Sharma BS, Mahapatra AK. Role of simvastatin in prevention of vasospasm and improving functional outcome after aneurysmal sub-arachnoid hemorrhage: a prospective, randomized, double-blind, placebo-controlled pilot trial. Br J Neurosurg. 2013;27(2):181–6.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Wesali S, Persson HC, Cederin B, Sunnerhagen KS. Improved survival after non-traumatic subarachnoid haemorrhage with structured care pathways and modern intensive care. Clin Neurol Neurosurg. 2015;138:52–8.

    Article  PubMed  Google Scholar 

  4. 4.

    Banki NM, Kopelnik A, Dae MW, Miss J, Tung P, Lawton MT, Drew BJ, Foster E, Smith W, Parmley WW, Zaroff JG. Acute neurocardiogenic injury after subarachnoid hemorrhage. Circulation. 2005;112(21):3314–9.

    Article  PubMed  Google Scholar 

  5. 5.

    Johnston SC, Selvin S, Gress DR. The burden, trends, and demographics of mortality from subarachnoid hemorrhage. Neurology. 1998;50(5):1413–8.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Suarez JI, Tarr RW, Selman WRJS. Aneurysmal subarachnoid hemorrhage. N Engl J Med. 2006;354(4):387.

    CAS  Article  Google Scholar 

  7. 7.

    Zumofen DW, Roethlisberger M, Achermann R, et al. Factors associated with clinical and radiological status on admission in patients with aneurysmal subarachnoid hemorrhage. Neurosurg Rev. 2018;41(4):1059–69.

    Article  PubMed  Google Scholar 

  8. 8.

    Diringer MN, Zazulia AR. Aneurysmal subarachnoid hemorrhage: strategies for preventing vasospasm in the intensive care unit. Semin Respir Crit Care Med. 2017;38(6):760–7.

    Article  PubMed  Google Scholar 

  9. 9.

    Leclerc JL, Blackburn S, Neal D, Mendez NV, Wharton JA, Waters MF, Doré S. Haptoglobin phenotype predicts the development of focal and global cerebral vasospasm and may influence outcomes after aneurysmal subarachnoid hemorrhage. Proc Natl Acad Sci U S A. 2015;112(4):1155–60.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Hafeez S, Grandhi R. Systematic review of Intrathecal Nicardipine for the treatment of cerebral vasospasm in aneurysmal subarachnoid hemorrhage. Neurocrit Care. 2019;31(2):399–405.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Shen J, Huang KY, Zhu Y, Pan JW, Jiang H, Weng YX, Zhan RY. Effect of statin treatment on vasospasm-related morbidity and functional outcome in patients with aneurysmal subarachnoid hemorrhage: a systematic review and meta-analysis. J Neurosurg. 2017;127(2):291–301.

    Article  PubMed  Google Scholar 

  12. 12.

    Tseng MY, Czosnyka M, Richards H, Pickard JD, Kirkpatrick PJ. Effects of acute treatment with pravastatin on cerebral vasospasm, autoregulation, and delayed ischemic deficits after aneurysmal subarachnoid hemorrhage: a phase II randomized placebo-controlled trial. Stroke. 2005;36(8):1627–32.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Pera J, Korostynski M, Golda S, Piechota M, Dzbek J, Krzyszkowski T, Dziedzic T, Moskala M, Przewlocki R, Szczudlik A, Slowik A. Gene expression profiling of blood in ruptured intracranial aneurysms: in search of biomarkers. J Cereb Blood Flow Metab. 2013;33(7):1025–31.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Lai NS, Zhang JQ, Qin FY, Sheng B, Fang XG, Li ZB. Serum microRNAs are non-invasive biomarkers for the presence and progression of subarachnoid haemorrhage. Biosci Rep. 2017;37(1):BSR20160480.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Liang Z, Chi YJ, Lin GQ, et al. LncRNA MEG3 participates in neuronal cell injury induced by subarachnoid hemorrhage via inhibiting the Pi3k/Akt pathway. Eur Rev Med Pharmacol Sci. 2018;22(9):2824–31.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Carvalho BS, Irizarry RA. A framework for oligonucleotide microarray preprocessing. Bioinformatics. 2010;26(19):2363–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Yates B, Braschi B, Gray KA, Seal RL, Tweedie S, Bruford EA. The HGNC and VGNC resources in 2017. Nucleic Acids Res. 2017;45(D1):D619–d625.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Szekely GJ, Rizzo ML. Hierarchical clustering via joint between-within distances: extending Ward’s minimum variance method. J Classif. 2005;22(2):151–83.

    Article  Google Scholar 

  19. 19.

    da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.

    CAS  Article  Google Scholar 

  20. 20.

    Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. 2003;4(5):P3.

    Article  PubMed  Google Scholar 

  21. 21.

    Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9(1):559.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Carpenter CR, Hussain AM, Ward MJ, Zipfel GJ, Fowler S, Pines JM, Sivilotti MLA. Spontaneous subarachnoid hemorrhage: a systematic review and meta-analysis describing the diagnostic accuracy of history, physical examination, imaging, and lumbar puncture with an exploration of test thresholds. Acad Emerg Med Off J Soc Acad Emerg Med. 2016;23(9):963–1003.

    Article  Google Scholar 

  25. 25.

    Wang Q, Luo Q, Yang Z, Zhao YH, Li J, Wang J, Piao J, Chen X. Weighted gene co-expression network analysis identified six hub genes associated with rupture of intracranial aneurysms. PLoS One. 2020;15(2):e0229308.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Zhao H, Li ST, Zhu J, Hua XM, Wan L. Analysis of peripheral blood cells’ Transcriptome in patients with subarachnoid hemorrhage from ruptured aneurysm reveals potential biomarkers. World Neurosurg. 2019;129:e16–22.

    Article  PubMed  Google Scholar 

  27. 27.

    Bo L, Wei B, Wang Z, et al. Screening of critical genes and MicroRNAs in blood samples of patients with ruptured intracranial aneurysms by bioinformatic analysis of gene expression data. Med Sci Monit. 2017;23:4518–25.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Pentimalli L, Modesti A, Vignati A, Marchese E, Albanese A, di Rocco F, Coletti A, di Nardo P, Fantini C, Tirpakova B, Maira G. Role of apoptosis in intracranial aneurysm rupture. J Neurosurg. 2004;101(6):1018–25.

    Article  PubMed  Google Scholar 

  29. 29.

    Pluta RM, Hansen-Schwartz J, Dreier J, Vajkoczy P, Macdonald RL, Nishizawa S, Kasuya H, Wellman G, Keller E, Zauner A, Dorsch N, Clark J, Ono S, Kiris T, LeRoux P, Zhang JH. Cerebral vasospasm following subarachnoid hemorrhage: time for a new world of thought. Neurol Res. 2009;31(2):151–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Fujii M, Yan J, Rolland WB, Soejima Y, Caner B, Zhang JH. Early brain injury, an evolving frontier in subarachnoid hemorrhage research. Transl Stroke Res. 2013;4(4):432–46.

    Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Yuksel S, Tosun YB, Cahill J, Solaroglu I. Early brain injury following aneurysmal subarachnoid hemorrhage: emphasis on cellular apoptosis. Turk Neurosurg. 2012;22(5):529–33.

    Article  PubMed  Google Scholar 

  32. 32.

    Wang Z, Zuo G, Shi XY, Zhang J, Fang Q, Chen G. Progesterone administration modulates cortical TLR4/NF-κB signaling pathway after subarachnoid hemorrhage in male rats. Mediat Inflamm. 2011;2011:848309.

    CAS  Article  Google Scholar 

  33. 33.

    Hasegawa Y, Suzuki H, Sozen T, et al. Apoptotic mechanisms for neuronal cells in early brain injury after subarachnoid hemorrhage. Acta Neurochir Suppl. 2011;110(Pt 1):43–8.

    Article  PubMed  Google Scholar 

  34. 34.

    Zhang D, Yan H, Li H, Hao S, Zhuang Z, Liu M, Sun Q, Yang Y, Zhou M, Li K, Hang C. TGFβ-activated kinase 1 (TAK1) inhibition by 5Z-7-oxozeaenol attenuates early brain injury after experimental subarachnoid hemorrhage. J Biol Chem. 2015;290(32):19900–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Fann DY, Lim YA, Cheng YL, et al. Evidence that NF-κB and MAPK signaling promotes NLRP Inflammasome activation in neurons following ischemic stroke. Mol Neurobiol. 2018;55(2):1082–96.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Crack PJ, Taylor JM, Ali U, et al. Potential contribution of NF-kappaB in neuronal cell death in the glutathione peroxidase-1 knockout mouse in response to ischemia-reperfusion injury. Stroke. 2006;37(6):1533–8.

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    You WC, Wang CX, Pan YX, Zhang X, Zhou XM, Zhang XS, Shi JX, Zhou ML. Activation of nuclear factor-κB in the brain after experimental subarachnoid hemorrhage and its potential role in delayed brain injury. PLoS One. 2013;8(3):e60290.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  38. 38.

    Pathak A, Stanley EM, Hickman FE, et al. Retrograde degenerative signaling mediated by the p75 neurotrophin receptor requires p150(glued) deacetylation by axonal HDAC1. Dev Cell. 2018;46(3):376–387.e377.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Shen J, Chen X, Li H, Wang Y, Huo K, Ke K. p75 neurotrophin receptor and its novel interaction partner, NIX, are involved in neuronal apoptosis after intracerebral hemorrhage. Cell Tissue Res. 2017;368(1):13–27.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Li T, Liu H, Xue H, Zhang J, Han X, Yan S, Bo S, Liu S, Yuan L, Deng L, Li G, Wang Z. Neuroprotective effects of hydrogen sulfide against early brain injury and secondary cognitive deficits following subarachnoid hemorrhage. Brain Pathol. 2017;27(1):51–63.

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    Stanfill A, Simpson C, Sherwood P, Poloyac S, Crago E, Kim H, Conley Y. A pilot study on the impact of dopamine, serotonin, and brain-derived neurotrophic factor genotype on long-term functional outcomes after subarachnoid hemorrhage. SAGE Open Med. 2017;5:2050312117726725.

    Article  PubMed  PubMed Central  Google Scholar 

  42. 42.

    Tian X, Sun L, Feng D, Sun Q, Dou Y, Liu C, Zhou F, Li H, Shen H, Wang Z, Chen G. HMGB1 promotes neurovascular remodeling via rage in the late phase of subarachnoid hemorrhage. Brain Res. 2017;1670:135–45.

    CAS  Article  PubMed  Google Scholar 

Download references





Author information




WL and DF participated in the design of this study, and they both performed the statistical analysis. ZR carried out the study and collected important background information. QYL drafted the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Qiaoying Li.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no conflict of interest.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1:

Supplementary Table S1. The list of significantly differentially expressed RNAs (DERs) (lncRNA and mRNAs) in SAH samples relative to controls.

Additional file 2:

Supplementary Table S2. The list of sample information.

Additional file 3:

Supplementary Table S3. The list of significantly differentially expressed RNAs (DERs) included in the 6 WGCNA modules.

Additional file 4:

Supplementary Table S4. The list of the lncRNA-mRNA pairt with Pearson correlation coefficients (PCC) higher than 0.6.

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 The Creative Commons Public Domain Dedication waiver ( 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

Verify currency and authenticity via CrossMark

Cite this article

Leng, W., Fan, D., Ren, Z. et al. Identification of upregulated NF-κB inhibitor alpha and IRAK3 targeting lncRNA following intracranial aneurysm rupture-induced subarachnoid hemorrhage. BMC Neurol 21, 197 (2021).

Download citation