- Research article
- Open Access
Construction of competitive endogenous RNA network reveals regulatory role of long non-coding RNAs in intracranial aneurysm
BMC Neuroscience volume 22, Article number: 15 (2021)
Rupture of intracranial aneurysm (IA) is the main cause of devastating subarachnoid hemorrhage, which urges our understanding of the pathogenesis and regulatory mechanisms of IA. However, the regulatory roles of long non-coding RNAs (lncRNAs) in IA is less known.
We processed the raw SRR files of 12 superficial temporal artery (STA) samples and 6 IA samples to count files. Then the differentially expressed (DE) mRNAs, miRNAs, and lncRNAs between STAs and IAs were identified. The enrichment analyses were performed using DEmRNAs. Next, a lncRNA-miRNA-mRNA regulatory network was constructed using integrated bioinformatics analysis. In summary, 341 DElncRNAs, 234 DEmiRNAs, and 2914 DEmRNAs between the STA and IA. The lncRNA-miRNA-mRNA regulatory network of IA contains 91 nodes and 146 edges. The subnetwork of hub lncRNA PVT1 was extracted. The expression level of PVT1 was positively correlated with a majority of the mRNAs in its subnetwork. Moreover, we found that several mRNAs (CCND1, HIF1A, E2F1, CDKN1A, VEGFA, COL1A1 and COL5A2) in the PVT1 subnetwork served as essential components in the PI3K-Akt signaling pathway, and that some of the non-coding RNAs (ncRNAs) (PVT1, HOTAIR, hsa-miR-17, hsa-miR-142, hsa-miR-383 and hsa-miR-193b) interacted with these mRNAs.
Our annotations noting ncRNA’s role in the pathway may uncover novel regulatory mechanisms of ncRNAs and mRNAs in IA. These findings provide significant insights into the lncRNA regulatory network in IA.
Rupture of intracranial aneurysm (IA) is the main cause of subarachnoid hemorrhage (SAH), leading to exceedingly high mortality and morbidity [1, 2]. The prevalence of IA in the general population is reported to be 3.2% . Approximately 85% of non-traumatic SAH is caused by IA rupture . Moreover, although unruptured IAs (UIs) are typically asymptomatic, UIs are being detected more and more, and have become an important health issue . Since surgical clipping and intravascular coiling are invasive procedures with potentially serious complications, the effective management of UIs remains a challenge . A deeper understanding of the pathogenesis of IA helps to find more effective treatment for IA.
Recently, an increasing number of studies have focused on the competitive endogenous RNA (ceRNA) network’s impact in cardiac hypertrophy, type 2 diabetes mellitus, and various types of cancer [5,6,7,8]. The hypothesis of ceRNA described a molecular regulatory mechanism for post-transcriptional regulation. Several studies revealed that lncRNAs can function as microRNA (miRNA) sponges in a ceRNA regulatory network, and further regulate mRNA expression levels . LncRNAs are non-coding RNA of greater than 200 nucleotides in length, which were reported to play a crucial role in the biological processes . However, the potential role of the ceRNA regulatory network in IA formation remains unclear.
In this study, we identified differentially expressed lncRNAs, miRNAs, and mRNAs between STAs and IAs. Then, we developed a lncRNA-miRNA-mRNA regulatory network. Previous study found that knockdown of lncRNA PVT1 could reverse the murine angiotensin II-induced abdominal aortic aneurysm (AAA) associated alterations, including inhibition of vascular smooth muscle cell (VSMC) apoptosis and ECM disruption . However, the regulatory mechanism of lncRNA PVT1 in development of IA was not clear. Thus, lncRNA PVT1 was selected for further analysis. The correlations between PVT1 and mRNAs in the PVT1 subnetwork were analyzed. After a functional enrichment analysis, we found that mRNAs in the PVT1 subnetwork were annotated into the PI3K-Akt signaling pathway. In addition, we found that several ncRNAs were connected with these mRNAs, and we annotated these ncRNAs to the PI3K-Akt signaling pathway, which may uncover novel regulatory mechanisms of ncRNAs and mRNAs in IA.
Materials and methods
Patients and RNA-seq data processing
The SRR files (SRR1819905–SRR1819922), including 12 superficial temporal artery (STA) samples and 6 intracranial aneurysm (IA) samples, were retrieved from GEO NCBI  for single ended data, which were converted to ‘fastq’ format data using sratoolkit (version 2.8.2). The reference genome of human (GRCh38.p13) was downloaded from Ensembl (https://asia.ensembl.org/). STAR (version 2.4.0i) , a highly efficient RNAseq alignment tool. This was used to align the ‘fastq’ format data with the human reference genome, using the custom default parameters. The SAM files from STAR alignment were converted to the BAM format by using SAMtools (version 1.3.1) . The BAM files were converted to counts files by using HTSeq .We generated “.csv” files, each consisting of all gene counts for that particular sample. Furthermore, we combined all of the “.csv” files and obtained a single file with sample names depicted as columns, and gene names depicted as rows. Then, the counts files were normalized using the edgER package . LncRNAs and mRNAs were annotated using the Ensembl database . Regarding the microRNA data in GSE66240, the raw probe-level data was preprocessed through the robust multi-array average (RMA) algorithm in the Affy package . For genes corresponding to multiple probes, we used the average probe value as the expression level . The missing data in these gene expression matrices were imputed with the k-Nearest Neighbor (KNN) method (k = 10) . As all the data were retrieved from the GEO database, the approval from the local Ethics Committee was not needed.
Identification of DElncRNAs, DEmiRNAs and DEmRNAs
Counts files, containing expression of lncRNAs and mRNAs, were normalized by the edgER package. Using the limma package, DElncRNAs, DEmiRNAs, and DEmRNAs were identified. lncRNAs, miRNAs, and mRNAs with |Log2(fold change)|> 1 and adjusted P-values < 0.05 were considered differentially expressed lncRNAs, miRNAs, and mRNAs.
Functional enrichment analysis
Both the Gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed using the Database for Annotation, Visualization, and Integrated Discovery platform (DAVID 6.8, https://david.ncifcrf.gov/) . DEmRNAs were utilized for GO and KEGG pathway enrichment analyses. The pathways and GO terms, with corrected P-values < 0.05 using the Benjamini method, were considered significant categories.
Construction of a ceRNA regulatory network
Figure 1 depicts the flow chart of the ceRNA network construction. Firstly, the lncRNA-miRNA potential interaction file (Highly conserved microRNA families) was downloaded from miRcode . MiRcode provides a whole transcriptome human microRNA target prediction based on comprehensive GENCODE gene annotation, including > 10,000 long non-coding RNA genes. Secondly, DElncRNAs were put into the miRcode database and lncRNA-miRNA pairs were identified. Then, we eliminated miRNAs that expressed no difference between the STA and IA tissues. Thirdly, miRNA-targeted mRNAs were identified using three miRNA reference databases: miRDB, miRTarBase, and TargetScan, which contain predictive or experimental validated miRNA-targeted mRNAs [22,23,24]. Only mRNAs appearing in the three databases were defined as miRNA-targeted mRNAs for increasing the prediction reliability. Furthermore, targeted mRNAs that expressed no difference between STA and IA tissues were filtered out. The ceRNA network was constructed and viewed using Cytoscape (http://www.cytoscape.org/). The relationship between PVT1 and mRNA expression level was analyzed with Pearson correlation. The heatmaps based on DEmRNAs, DEmiRNAs and DElncRNAs included in the ceRNA regulatory network were generated by “pheatmap” package.
GraphPad Prism (version 6.0, GraphPad Software, San Diego, CA, USA) and R language (3.4.0) were used for statistical analysis. A P-value < 0.05 was considered statistically significant.
Differentially expressed RNAs in IA
Compared with the superficial temporal artery (STA) tissues, a total 2914 differentially expressed mRNA, 234 miRNA and 341 lncRNA were identified in intracranial aneurysm tissues. Of these, 201 (58.8%) lncRNAs, 1.807 (62%) mRNAs and 10 (4.3%) miRNAs were upregulated, whereas 141 (41.2%) lncRNAs, 1.107 (38%) mRNAs, and 224 (95.7%) miRNAs were downregulated in the tissue of intracranial aneurysm, when compared with STA tissues.
Functional enrichment analysis of DEmRNAs
GO enrichment and KEGG pathway analyses were performed to explore the potential functions of the 2.914 differentially expressed mRNAs. In the GO enrichment analysis, a total of 378 enriched GO terms have been identified in the Biological Process (BP). And the top 10 significantly enriched terms are shown in Fig. 2a. Moreover, the DEmRNAs were primarily enriched in collagen and extracellular matrix-related BPs, such as “extracellular matrix organization”, “collagen catabolic process”, “collagen fibril organization”, and “cell adhesion”. Moreover, we also found that these DEmRNAs were also enriched in “angiogenesis”, “skeletal system development”, “muscle contraction”, and “inflammatory response”, which indicates that these BPs may play roles in IA formation. In addition, a total 53 enriched pathways were identified after the KEGG pathway analysis. The top 10 enriched pathways are shown in Fig. 2b. And “vascular smooth muscle contraction” is closely associated with IA formation in these pathways. Furthermore, the enriched “focal adhesion”, “ECM-receptor interaction”, and “PI3K-Akt signaling pathway” indicated that the PI3K-Akt signaling pathway may play an important role in IA, which was further analyzed in this study.
Construction of a ceRNA regulatory network in IA
A lncRNA‐miRNA‐mRNA regulatory network of IA was constructed (Fig. 3) for exploring the regulatory mechanism of IA. The mRNA-miRNA and lncRNA‐miRNA relationship pairs (Additional file 1: Table S1, Additional file 2: Table S2) were combined into the ceRNA network, following a negative regulation pattern. The positively co-expressed mRNA-miRNA and lncRNA‐miRNA pairs were excluded. Finally, the network was constructed with 91 nodes (60 mRNAs, 9 miRNAs and 22 lncRNAs) and 146 edges. The degree of a node is defined as the total number of edges connecting it in the network, with an increased degree indicating that the node is highly connected, and likely to be a hub gene. The heatmaps based on DEmRNAs, DEmiRNAs and DElncRNAs included in the ceRNA regulatory network were shown in Fig. 4. Among the DEmiRNAs in the network, hsa-miR-17 has the highest degree (degree = 52), indicating that 52 nodes directly connect to it (Fig. 5). Moreover, other DEmiRNAs that connect lncRNAs and mRNAs in the network are also shown (Fig. 5). The highest degree DElncRNAs were PVT1, NEAT1 and KCNQ1OT1 (degree = 8).
The components of PVT1 subnetwork
PVT1 was the core lncRNA, and had the highest degree in the network. A previous study found that knockdown of lncRNA PVT1 reversed the murine angiotensin II-induced AAA associated alterations, such as attenuation of aortic diameter dilation and marked adventitial thickening . To explore the regulatory mechanism of PVT1 in IA, we focused on the PVT1 subnetwork. The components of the PVT1 subnetwork are shown in Fig. 6a. The predicted or validated binding sites of hsa-miR-17, hsa-miR-143, hsa-miR-152, hsa-miR-761 and hsa-miR-383 with PVT1 are shown (Fig. 6b–e).
Interaction between lncRNA and mRNA in the ceRNA network
LncRNA can interact with mRNAs indirectly following post-transcriptional regulatory mechanisms. We analyzed the correlation between PVT1 and mRNA’s expression levels in its subnetwork, which showed that there were strong positive correlations. There were 52 mRNAs interacting indirectly with PVT1. The correlation coefficients between 38.5% (20/52) mRNAs and PVT1 were more than 0.6 (Additional file 3: Table S3), and the correlation coefficients between 90.4% (47/52) mRNAs and PVT1 were more than 0.3 (Additional file 3: Table S3). Scatter plots of the top 10 mRNAs with the highest correlation coefficients are shown in Fig. 7.
The non-coding RNAs-involved PI3K-AKT signaling pathway
The KEGG pathway analysis revealed that the PI3K-AKT signaling pathway was significantly enriched, indicating that the PI3K-AKT signaling pathway could play an important role in IA formation. In addition, we found that several mRNAs in PVT1 subnetwork were curial genes in the PI3K-AKT signaling pathway, such as CCND1, HIF1A, E2F1, CDKN1A, VEGFA, COL1A1, and COL5A2. Moreover, these mRNAs were directly connected with hsa-miR-17, hsa-miR-152, hsa-miR-383, and hsa-miR-193b, and these miRNAs were connected with PVT1 and HOTAIR (Fig. 8a). hsa-miR-17, hsa-miR-143, hsa-miR-193b and hsa-miR-383 and their predicted or validated binding sites with these mRNAs were shown (Fig. 8b–i). hsa-miR-17, hsa-miR-143, hsa-miR-193b, hsa-miR-761, hsa-miR-152 and their predicted or validated binding sites with LncRNA HOTAIR was shown in Fig. 8j–m. To explore how these non-coding RNAs (ncRNAs) regulate the PI3K-AKT signaling pathway and further affect IA formation, we mapped these ncRNAs to the PI3K-AKT signaling pathway (Fig. 9). Regulation of E2F transcription factor 1 (E2F1) by topoisomerase 2-binding protein 1 (TOPBP1) involves the PI3K-AKT signaling pathway, which can mediate both cell proliferation and apoptosis . Hypoxia inducible Factor 1 Subunit Alpha (HIF1A) is regulated by PI3K-AKT-MTOR signaling pathway, and further promotes angiogenesis . CDKN1A and CCND1 can regulate cell cycle, which are also part of the PI3K-AKT signaling pathway . VEGFA can bind to VEGFR1/2 and further activate PI3K-AKT signaling pathway .
LncRNA are a type of non-coding RNA with more than 200 bp transcript length. Various studies have found that lncRNAs play important roles in different types of diseases. Large amounts of lncRNAs have been annotated, however, the functions of them, especially in diseases, still need to be further explored. Many previous studies focused on ncRNA-mediated regulation of mRNAs, but recent studies have indicated that lncRNAs could regulate mRNAs via miRNAs and further form a regulatory ceRNA network. The ceRNA hypothesis refers that the lncRNAs and mRNAs compete to bind miRNAs and further regulate each other’s expression levels. Recently, various studies have focused on the functions of ceRNA networks in cardiac hypertrophy, type 2 diabetes mellitus, and various types of cancer [5,6,7,8]. However, the potential role of the ceRNA network in IA formation remains unclear.
In this study, all of the DEmiRNAs in the ceRNA network were downregulated in IAs relative to STAs. According to ceRNA theory, the down-regulated lncRNAs and mRNAs were excluded. Previous studies found that downregulation of miR-143/145 cluster could upregulate their targeted gene KLF5, and further promote proliferation and migration of VSMCs . In our study, we found that miR-143 interacted with lncRNA PVT1 and HOTAIR, and that miR-143 could target COL1A1 and COL5A2 and further affect the PI3K-Akt signaling pathway. Thus, we may provide a potential axis PVT1/HOTAIR-miR-143-COL1A1/COL5A2 as an upstream regulator of the PI3K-Akt signaling pathway.
The PI3K-Akt signaling pathway is activated by various types of cellular stimuli or toxic insults, and regulates fundamental cellular functions such as transcription, translation, proliferation, growth, and survival. Dysregulation in the PI3K-Akt signaling pathway is implicated in various diseases such as cancer, type 2 diabetes mellitus, and IA [7, 30, 31]. Previous research reported that activation of PI3K-Akt signaling pathway promoted proliferation of vascular smooth muscle cells which were found to be involved in aortic aneurysms . Inhibition of VEGF/PI3K/Akt signaling pathway mediated by miR-195 could suppress formation of abdominal aortic aneurysm . In addition, the PI3K-Akt signaling pathway was less studied in IA formation. Sun et al. found that bone marrow mesenchymal stem cells-derived exosomes supressed the activation of the PI3k/Akt/NF-κB signaling pathway and maintained Th17/Treg balance, which could inhibite development of IA. Recently, various studies revealed that various lncRNAs affect the PI3K-Akt signaling pathway by targeting miRNAs, and further regulating mRNAs in the pathway [34, 35].
In our study, we found that PVT1 and HOTAIR were involved in the regulation of PI3K-Akt signaling pathway by targeting different miRNAs and regulating upstream/downstream effectors (COL1A1, COL5A2, VEGFA, CCND1, HIF1A, E2F1 and CDKN1A). Several members of collagen family, such as COL1A1, COL5A2, were overexpressed in IAs, which were consistent with previous studies [36, 37]. And these two genes were involved in collagen formation that could result in the ECM remodeling of intracranial blood vessels, which promotes formation of IA . VEGFA, related to aSAH, was reported to be associated with thrombospondin and vascular endothelial growth factor . A previous study revealed that HIF1A may be associated with rupture of human saccular IA wall . Another research found that HIF1A was pivotal for the development of AAA . HIF1A could upregulate expression levels of MMP-2 and MMP-9, and further promote aneurysmal progression . Down-regulation of PVT1 could reduce proinflammatory cytokines, MMP-2 and MMP-9, increase TIMP-1, and further reverse angiotensin II-induced AAA-associated alterations in mice . In this study, we annotated these ncRNAs to PI3K-Akt signaling pathway. These findings establish novel connections among lncRNAs, miRNAs, and mRNAs in PI3K-Akt signaling pathway, which play an important role in IA.
In this study, we also found that there were strong correlations between PVT1 expression and several genes’ expression levels (Fig. 7 and Additional file 3: Table S3), including DNMT1, E2F1 et al. Jin et al. reported that PVT1 could recruit DNMT1 via EZH2 to miR-18b-5p DNA promoter and inhibited the miR-18b-5p transcription via DNA methylation . Previous study found that PVT1 could regulate E2F1 expression levels in tumor development . These genes might also play roles in development of IA, which needs more verification in future research. The potential axis PVT1/HOTAIR-miR-143-COL1A1/COL5A2 needs further verification in future research, and the function of the axis in formation of IA needs to be explored.
First, we identified differentially expressed lncRNAs, miRNAs and mRNAs between STAs and IAs. Functional enrichment analysis of differentially expressed mRNAs was performed. We further developed a lncRNA–miRNA–mRNA regulatory network. Then, the subnetworks of each core miRNAs were extracted and analyzed. Moreover, PVT1, the hub lncRNA in the ceRNA network, was reported to be related to aneurysm. The subnetwork of PVT1 was extracted and analyzed. Furthermore, the correlations of expression levels between PVT1 and mRNAs in its subnetwork were analyzed. In addition, we found that several mRNAs in the PVT1 subnetwork, including CCND1, HIF1A, E2F1, CDKN1A, VEGFA, COL1A1, and COL5A2, were an important part of the PI3K-Akt signaling pathway. We found that several ncRNAs were connected with these mRNAs. Then, we annotated these ncRNAs to the PI3K-Akt signaling pathway. This study may provide a comprehensive view of the underlying mechanisms of gene regulation and interaction in IA. Additionally, this study also revealed that ncRNAs involved in the PI3K-Akt signaling pathway could play an important role in IA formation.
Availability of data and materials
The expression data analyzed in this study are available in GEO public repository (GSE66240).
Long noncoding RNA
Superficial temporal artery
Competitive endogenous RNA
Robust multi-array average (RMA)
Kyoto Encyclopedia of Genes and Genomes
Abdominal aortic aneurysm
E2F transcription factor 1
Topoisomerase 2-binding protein 1
Vlak MH, Algra A, Brandenburg R, Rinkel GJ. Prevalence of unruptured intracranial aneurysms, with emphasis on sex, age, comorbidity, country, and time period: a systematic review and meta-analysis. Lancet Neurol. 2011;10(7):626–36.
Westerlaan HE, van Dijk JM. Intracranial aneurysms in patients with subarachnoid hemorrhage: CT angiography as a primary examination tool for diagnosis–systematic review and meta-analysis. Radiology. 2011;258(1):134–45.
van Gijn J, Kerr RS, Rinkel GJ. Subarachnoid haemorrhage. Lancet (London, England). 2007;369(9558):306–18.
Lylyk P, Ferrario A, Pasbon B, Miranda C, Doroszuk G. Buenos Aires experience with the Neuroform self-expanding stent for the treatment of intracranial aneurysms. J Neurosurg. 2005;102(2):235–41.
Wang JD, Zhou HS, Tu XX, He Y, Liu QF, Liu Q, Long ZJ. Prediction of competing endogenous RNA coexpression network as prognostic markers in AML. Aging. 2019;11(10):3333–47.
Gao S, Chen Y, Xu B, Yu C, Yue M, Tan X, Zhang J, Feng C, Song C, Ai B, et al. Identification and analysis of a key long non-coding RNAs (lncRNAs)-associated module reveal functional lncRNAs in cardiac hypertrophy. J Cell Mol Med. 2018;22(2):892–903.
Lin Z, Li X, Zhan X. Construction of competitive endogenous RNA network reveals regulatory role of long non-coding RNAs in type 2 diabetes mellitus. J Cell Mol Med. 2017;21(12):3204–13.
Pan H, Pan J, Song S, Ji L, Lv H, Yang Z. Identification and development of long non-coding RNA-associated regulatory network in colorectal cancer. J Cell Mol Med. 2019;23(8):5200–10.
Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language? Cell. 2011;146(3):353–8.
Wahlestedt C. Targeting long non-coding RNA to therapeutically upregulate gene expression. Nat Rev Drug Discovery. 2013;12(6):433–46.
Zhang Z, Zou G, Chen X, Lu W, Liu J, Zhai S, Qiao G. Knockdown of lncRNA PVT1 inhibits vascular smooth muscle cell apoptosis and extracellular matrix disruption in a murine abdominal aortic aneurysm model. Mol Cells. 2019;42(3):218–27.
Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30(1):207–10.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Batut P, Chaisson M, Gingeras TR. STAR: ultrafast universal RNA-seq aligner. Bioinformatics (Oxford, England). 2013;29(1):15–21.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. The Sequence Alignment/Map format and SAMtools. Bioinformatics (Oxford, England). 2009;25(16):2078–9.
Anders S, Pyl PT, Huber W. HTSeq–a Python framework to work with high-throughput sequencing data. Bioinformatics (Oxford, England). 2015;31(2):166–9.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics (Oxford, England). 2010;26(1):139–40.
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics (Oxford, England). 2003;4(2):249–64.
Li W, Li K, Zhao L, Zou H. Bioinformatics analysis reveals disturbance mechanism of MAPK signaling pathway and cell cycle in Glioblastoma multiforme. Gene. 2014;547(2):346–50.
Garcia-Laencina PJ, Abreu PH, Abreu MH, Afonoso N. Missing data imputation on the 5-year survival prediction of breast cancer patients with unknown discrete values. Comput Biol Med. 2015;59:125–33.
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.
Jeggari A, Marks DS, Larsson E. miRcode: a map of putative microRNA target sites in the long non-coding transcriptome. Bioinformatics (Oxford, England). 2012;28(15):2062–3.
Agarwal V, Bell GW. Nam JW. Bartel DP: Predicting effective microRNA target sites in mammalian mRNAs. eLife; 2015. p. 4.
Chou CH, Shrestha S, Yang CD, Chang NW, Lin YL, Liao KW, Huang WC, Sun TH, Tu SJ, Lee WH, et al. miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. 2018;46(D1):D296-d302.
Wong N, Wang X. miRDB: an online resource for microRNA target prediction and functional annotations. Nucleic acids Res. 2015;43:D146-152.
Liu K, Paik JC, Wang B, Lin FT, Lin WC. Regulation of TopBP1 oligomerization by Akt/PKB for cell survival. EMBO J. 2006;25(20):4795–807.
Zhang L, Chen C, Duanmu J, Wu Y, Tao J, Yang A, Yin X, Xiong B, Gu J, Li C, et al. Cryptotanshinone inhibits the growth and invasion of colon cancer by suppressing inflammation and tumor angiogenesis through modulating MMP/TIMP system, PI3K/Akt/mTOR signaling and HIF-1alpha nuclear translocation. Int Immunopharmacol. 2018;65:429–37.
Bollaert E, de Rocca SA, Demoulin JB. The HMG box transcription factor HBP1: a cell cycle inhibitor at the crossroads of cancer signaling pathways. Cell Mol Life Sci. 2019;76(8):1529–39.
Sadremomtaz A, Mansouri K, Alemzadeh G, Safa M, Rastaghi AE, Asghari SM. Dual blockade of VEGFR1 and VEGFR2 by a novel peptide abrogates VEGF-driven angiogenesis, tumor growth, and metastasis through PI3K/AKT and MAPK/ERK1/2 pathway. Biochim Biophys Acta. 2018;1862(12):2688–700.
Xu J, Yan S, Tan H, Ma L, Feng H, Han H, Pan M, Yu L, Fang C. The miR-143/145 cluster reverses the regulation effect of KLF5 in smooth muscle cells with proliferation and contractility in intracranial aneurysm. Gene. 2018;679:266–73.
Li XG, Wang YB. SRPK1 gene silencing promotes vascular smooth muscle cell proliferation and vascular remodeling via inhibition of the PI3K/Akt signaling pathway in a rat model of intracranial aneurysms. CNS Neurosci Ther. 2019;25(2):233–44.
Martini M, De Santis MC, Braccini L, Gulluni F, Hirsch E. PI3K/AKT signaling pathway and cancer: an updated review. Ann Med. 2014;46(6):372–83.
Liu C, Su T, Li F, Li L, Qin X, Pan W, Feng F, Chen F, Liao D, Chen L. PI3K/Akt signaling transduction pathway is involved in rat vascular smooth muscle cell proliferation induced by apelin-13. Acta Biochim Biophys Sin (Shanghai). 2010;42(6):396–402.
Ma X, Yao H, Yang Y, Jin L, Wang Y, Wu L, Yang S, Cheng K. miR-195 suppresses abdominal aortic aneurysm through the TNF-α/NF-κB and VEGF/PI3K/Akt pathway. Int J Mol Med. 2018;41(4):2350–8.
Zhang W, Zhang Y, Xi S. Upregulation of lncRNA HAGLROS enhances the development of nasopharyngeal carcinoma via modulating miR-100/ATG14 axis-mediated PI3K/AKT/mTOR signals. Artif Cells Nanomed Biotechnol. 2019;47(1):3043–52.
Chen S, Chen H, Yu C, Lu R, Song T, Wang X, Tang W, Gao Y. Long noncoding RNA myocardial infarction associated transcript promotes the development of thoracic aortic by targeting microRNA-145 via the PI3K/Akt signaling pathway. J Cell Biochem. 2019;120(9):14405–13.
Bekelis K, Kerley-Hamilton JS, Teegarden A, Tomlinson CR, Kuintzle R, Simmons N, Singer RJ, Roberts DW, Kellis M, Hendrix DA. MicroRNA and gene expression changes in unruptured human cerebral aneurysms. J Neurosurg. 2016;125(6):1390–9.
Leeper N, Raiesdana A, Kojima Y, Chun H, Azuma J, Maegdefessel L. MicroRNA-26a is a novel regulator of vascular smooth muscle cell function. J Cell Physiol. 2011;226:1035–43.
Lopes KP, Vinasco-Sandoval T, Vialle RA, Paschoal FM Jr, Bastos V, Bor-Seng-Shu E, Teixeira MJ, Yamada ES, Pinto P, Vidal AF, et al. Global miRNA expression profile reveals novel molecular players in aneurysmal subarachnoid haemorrhage. Sci Rep. 2018;8(1):8786.
Kurki MI, Hakkinen SK, Frosen J, Tulamo R, Fraunberg M, Wong G, Tromp G, Niemela M, Hernesniemi J, Jaaskelainen JE, et al. Upregulated signaling pathways in ruptured human saccular intracranial aneurysm wall: an emerging regulative role of Toll-like receptor signaling and nuclear factor-kappaB, hypoxia-inducible factor-1A, and ETS transcription factors. Neurosurgery. 2011;68(6):1667–75.
Tsai SH, Huang PH, Hsu YJ, Peng YJ, Lee CH, Wang JC, Chen JW, Lin SJ. Inhibition of hypoxia inducible factor-1alpha attenuates abdominal aortic aneurysm progression through the down-regulation of matrix metalloproteinases. Sci Rep. 2016;6:28612.
Jin L, Cai Q, Wang S, Wang S, Wang J, Quan Z. Long noncoding RNA PVT1 promoted gallbladder cancer proliferation by epigenetically suppressing miR-18b-5p via DNA methylation. Cell Death Dis. 2020;11(10):871.
Saleembhasha A, Mishra S. Long non-coding RNAs as pan-cancer master gene regulators of associated protein-coding genes: a systems biology approach. PeerJ. 2019;7:e6388.
We would like to express our thanks to the National Library of Medicine for giving us the privilege to access the raw data of various GEO series.
This work was funded by the National Natural Science Foundation of China (Grant No: 81870916), National Natural Science Foundation of China (Grant No: 81701144).
Ethical approval and consent to participate
All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and, or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. As personal identifying information was not included in GEO database, thus the informed consent was not required in this study.
Consent for publication
The authors declare no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Pan, YB., Lu, J., Yang, B. et al. Construction of competitive endogenous RNA network reveals regulatory role of long non-coding RNAs in intracranial aneurysm. BMC Neurosci 22, 15 (2021). https://doi.org/10.1186/s12868-021-00622-7
- Bioinformatics analysis
- Intracranial aneurysm
- Competitive endogenous RNA
- Long non-coding RNA
- PI3K-Akt signaling pathway