- 1Department of Neurosurgery, Sun Yat-sen Memorial Hospital, Sun Yat-sen University, Guangzhou, China
- 2RNA Biomedical Institute, Sun Yat-sen Memorial Hospital, Sun Yat-sen University, Guangzhou, China
- 3Guangdong Province Key Laboratory of Brain Function and Disease, Sun Yat-sen University, Guangzhou, China
Background: Circular RNAs (circRNAs) may involve the formation and rupture of intracranial aneurysms (IA). Inflammation plays a vital role in the development and progression of IA, which can be reflected by aneurysm wall enhancement (AWE) on high-resolution vessel wall magnetic resonance imaging (HR-VWI). This study aims to evaluate the role of circRNAs as the blood inflammatory biomarker for unruptured IA (UIA) patients with AWE on HR-VWI.
Methods: We analyzed the circRNA expression profiles in the peripheral blood samples among subjects from saccular UIA with AWE, UIA without AWE, and healthy controls by the circRNA microarray. The differential expression of hsa_circ_0007990 was assessed. We constructed the hsa_circ_0007990-microRNA-mRNA network and the regulatory axis of hub genes associated with the AWE in UIA.
Results: Eighteen patients harboring saccular UIAs with HR VWI and five healthy controls were included. We found 412 differentially expressed circRNAs between UIA patients and healthy controls by circRNA microarray. Two hundred thirty-one circRNAs were significantly differentially expressed in UIA patients with AWE compared with those without AWE. Twelve upregulated circRNAs were associated with AWE of UIA, including hsa_circ_0007990, hsa_circ_0114507, hsa_circ_0020460, hsa_circ_0053944, hsa_circ_0000758, hsa_circ_0000034, hsa_circ_0009127, hsa_circ_0052793, hsa_circ_0000301 and hsa_circ_0000729. The expression of hsa_circ_0007990 was increased gradually in the healthy control, UIA without AWE, and UIA with AWE confirmed by RT-PCR (P<0.001). We predicted 4 RNA binding proteins (Ago2, DGCR8, EIF4A3, PTB) and period circadian regulator 1 as an encoding protein with hsa_circ_0007990. The hsa_circ_0007990-microRNA-mRNA network containing five microRNAs (miR-4717-5p, miR-1275, miR-150-3p, miR-18a-5p, miR-18b-5p), and 97 mRNAs was constructed. The five hub genes (hypoxia-inducible factor 1 subunit alpha, estrogen receptor 1, forkhead box O1, insulin-like growth factor 1, CREB binding protein) were involved in the inflammatory response.
Conclusion: Differentially expressed blood circRNAs associated with AWE on HR-VWI may be the novel inflammatory biomarkers for assessing UIA patients. The mechanism of hsa_circRNA_0007990 for UIA progression needs to investigate further.
Introduction
Unruptured intracranial aneurysms (UIA) cause subarachnoid hemorrhage once bleeding, resulting in severe poor clinical outcomes. The risk of rupture in UIA significantly decreases after receiving endovascular treatment or surgical clipping in cases of aneurysm growth (1, 2). Inflammation plays a vital role in the development and progression of intracranial aneurysms (IA), which can be reflected by aneurysm wall enhancement (AWE) on high-resolution vessel wall magnetic resonance imaging (HR-VWI) (3–8). Recent studies suggest the association of inflammatory blood biomarkers in the aneurysm sac and peripheral blood with AWE on HR-VWI in patients with UIA (9–13).
Recent evidence shows that circular RNAs (circRNAs), one type of non-coding RNA, involve the formation and rupture of IA (14–19). CircRNAs regulate vascular smooth muscle cells phenotype involving the pathological process of IA (14–16). Differentially expressed circRNAs in the peripheral blood are associated with the aneurysm rupture status in patients with IA and may be used as potential clinical diagnostic biomarkers (20–22). However, the association of circRNAs in the peripheral blood with AWE on HR-VWI for IA remains unclear. In this study, we aim to evaluate the role of circRNAs as the blood inflammatory biomarker for UIA patients with AWE on HR-VWI.
Materials and methods
Patients selection
The institutional review board approved this study at Sun Yat-sen Memorial Hospital, Sun Yat-sen University. All participants signed their written informed consent to participate in this study. Patients with UIA for initial endovascular and surgical treatment from our HR-VWI aneurysm database were selected for this study between March 2018 and March 2021. Patients were older than 18 years old. Saccular UIA was confirmed by digital subtraction angiography and can be identified on MR angiography. Patients had images of HR-VWI without artifacts. We excluded patients with aneurysms in the extracranial or cavernous sinus of the internal carotid artery and with a size less than 3mm. Healthy adults were recruited for this study and had brain MR angiography examinations to exclude the presence of UIA and other intracranial lesions. All participants had no history of pneumonia, heart disease, cancer, autoimmune disease, hematological disease, chronic liver disease, kidney disease, or recent brain surgery. We recorded the data on medical history, laboratory examination, and clinical and radiological characteristics.
HR-VWI study
Patients with UIA underwent brain MR for HR-VWI on a 3.0T MRI scanner. We divided the pattern of AWE into none, focal, and circumferential AWE according to the precontrast and postcontrast HR-VWI T1-weighted sequences described in our previous study. The aneurysm-to-pituitary stalk contrast ratio (CRstalk) was calculated using the maximal signal intensity value of the aneurysm wall (SIwall) and the pituitary stalk (SIstalk) on the 3D postcontrast HR-VWI T1-weighted sequence as CRstalk = SIwall/SIstalk. CRstalk was used as an optimal quantitative analysis of AWE in UIA (8, 12).
Peripheral blood sample and RNA extraction
Before any antiplatelet drug administration or treatment, 2 ml peripheral venous blood was collected in the early morning after overnight fasting from all participants. We extracted total RNA from each peripheral whole blood sample using TRI Reagent BD as described in the manufacturer’s brochure (Molecular Research Center, Cincinnati, OH, United States) (23). RNA quality was assessed by the NanoDrop ND-1000 Spectrophotometer (Thermo Fisher Scientific, Waltham, MA, United States), and RNA with an OD260/OD280 ratio between 1.8 and 2.1 was considered good quality. Agarose gel electrophoresis was used to assess RNA integrity.
RNA labeling and hybridization
The RNase R (Epicentre, Inc.) was used to remove linear RNAs to enrich the circular RNAs. Circular RNAs were amplified, transcribed into cRNA, and purified using an RNeasy Mini Kit (Qiagen) according to the manufacturer’s protocols. The purified cRNA was labeled using Quick Amp Labeling Kit (Agilent) at 40°C for 2 hours. The NanoDrop ND-1000 was used to determine the yield and specific activity of the labeled cRNAs. Then, cDNA samples were hybridized onto the Human circRNA Arrays (V2.0, Arraystar Inc). Hybridization was performed at 65°C for 17 hours in an Agilent Hybridization Oven. Hybridized arrays were washed, fixed, and scanned with a G2505C Agilent Microarray Scanner.
circRNA microarray analysis
We obtained the original data by Agilent Feature Extraction software (version 11.0.1.1). The data were normalized and processed using the R software limma package. The expression level of circRNA was analyzed and compared among groups. False discovery rate (FDR) was calculated to account for multiple testing using R (version 3.6.7; R Foundation for Statistical Computing) (24). The difference in circRNAs expression was considered significant with Fold change (FC) ≥ 1.3 and an FDR-adjusted P value < 0.05.
Quantitative real-time PCR
qRT-PCR was utilized to assess the significantly differentially expressed circRNAs screened from microarray analysis. Total RNA was extracted and used as a template for reverse transcription to first-strand cDNA with SYBR Premix Ex Taq II (Tli RNaseH Plus; TaKaRa, Osaka, Japan). qRT-PCR with a total 10μl reaction volume was performed on the Quanstudio DX (ABI) machine using GAPDH as an endogenous reference for normalization (25). The primer sequences for circRNAs are listed in Supplementary Table 1. Each sample was tested three times. The 2−ΔΔCt method was applied to obtain the relative circRNA expression.
GO and KEGG pathway analysis
We used the circBank database to search microRNAs having a binding site with our target circRNA. Then, we predicted the simultaneous target mRNAs of the circRNA-targeted microRNAs by both TargetScan and miRDB databases. We performed the GO enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis to reveal the biological function of target genes predicted by the circRNA-binding microRNAs using the David (https://david.ncifcrf.gov). The GO enrichment analysis included the biological process, cellular component, and molecular function. Furthermore, the RNA binding proteins with target circRNA were searched by the CircInteractome database. The open reading frame of target circRNA and their potential encoding proteins were also analyzed using the online database.
Construction circRNA-microRNA-mRNA network with hub gene prediction
We identified the crucial circRNA and its targeted microRNAs and genes to construct the circRNA-microRNA-mRNA network using Cytoscape software (version 3.9.1) for visualization. Then, predicted genes from the top ten enriched KEGG pathways from the target circRNA-microRNAs were selected to construct the protein-protein interaction network by the Search Tool for the Retrieval of Interacting Genes (STRING, version 11.5; https://cn.string-db.org/cgi/input). After obtaining gene interaction for those target genes, the top ten hub genes were defined from the calculation by the maximal clique centrality algorithm in the CytoHubba plugin. Finally, we obtained the circRNA-microRNA-mRNA network containing ten key regulatory axes using the Cytoscape.
Statistical analysis
Continuous and categorical variables were analyzed by Student t-test, Mann-Whitney U test, Fisher’s exact, or chi-square test. A P-value less than 0.05 was considered significant for the results. Statistical analyses were performed using SPSS 22.0 software.
Results
CircRNA expression profiles in peripheral blood from UIA
We identified twenty-three subjects in this study, including 18 patients with UIAs and five healthy normal controls. The characteristics of all participants are shown in Table 1. We analyzed the circRNA expression profile of peripheral blood in six samples from UIA with AWE group, three from UIA without AWE, and three from normal controls by circRNA microarray. A total of 13 615 circRNAs were detected in the three groups using the Arraystar human circRNA Microarray V2.0. The differentially expressed peripheral blood circRNAs were generated and assessed by hierarchical clustering, scatter plots, and volcano plots between the UIA group and healthy control group (Figure 1A-C) and between UIA with AWE and UIA without AWE (Figure 1D-F).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1061592/fimmu-13-1061592-HTML-r1/image_m/fimmu-13-1061592-t001.jpg)
Table 1 Clinical characteristics of participants with unruptured intracranial aneurysm and healthy control.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1061592/fimmu-13-1061592-HTML-r1/image_m/fimmu-13-1061592-g001.jpg)
Figure 1 Differentially expressed circRNAs in peripheral blood identified by circRNAs microarray. (A) The hierarchical clustering of differentially expressed circRNAs between the unruptured intracranial aneurysm (UIA) and healthy control groups. Group A represented healthy control; group B represented the UIA group. (B) The scatterplot is used to assess the circRNA expression variation between the UIA and healthy control group of samples. The circRNAs above the top green line and below the green bottom line indicate more than a 1.3-fold change of circRNAs between the two compared samples. (C) Volcano Plots are used to visualize differential expression between the UIA and healthy control groups. The red point in the plot represents the differentially expressed circRNAs with statistical significance. (D) The hierarchical clustering of differentially expressed circRNAs between the UIA with aneurysm wall enhancement (AWE) and UIA without AWE. Group B1 represented UIA with AWE of samples; group B2 represented UIA without AWE of samples. (E) The scatterplot is used for assessing the circRNA expression variation between the UIA with AWE and UIA without AWE. The circRNAs above the top green line and below the green bottom line indicate more than a 1.3-fold change of circRNAs between the two compared samples. (F) Volcano Plots are used for visualizing differential expression between the UIA with AWE and UIA without AWE. The red point in the plot represents the differentially expressed circRNAs with statistical significance. (G) The genomic origin of the differentially expressed circRNAs among three groups of samples (UIA with AWE, UIA without AWE, and healthy controls). (H) A total of 12 significantly differentially expressed peripheral blood circRNAs among the three compared groups (UIA with AWE, UIA without AWE, and healthy controls) were identified by circRNA microarray.
We found 412 circRNAs significantly differentially expressed in the UIA group compared with the healthy control group (fold change ≥1.3, p-value <0.05). Among the differentially expressed circRNAs, 320 were upregulated while 92 were downregulated. We identified 231 circRNAs significantly differentially expressed in the UIA with AWE compared with the UIA without AWE. These included 154 upregulated circRNAs and 77 downregulated circRNAs. Differentially expressed circRNAs were classified into exonic, intronic, antisense, sense overlapping, and intergenic. The genomic origin of the differentially expressed circRNAs among three groups is shown in Figure 1G. Among the 412 differentially expressed circRNAs between the UIA and healthy control groups, 332 (80.58%) circRNAs are from exons. Among the 231 differentially expressed circRNAs between UIA with AWE and UIA without AWE, 183 (79.22%) circRNAs are from exons.
Differentially expressed CircRNAs in UIA with AWE
We found a total of 12 significantly differentially expressed peripheral blood circRNAs among the three compared groups identified by circRNA microarray (Figure 1H). Among these 12 upregulated circRNAs in UIA with AWE, 10 (83.33%) have been identified in circBase, including hsa_circ_0007990, hsa_circ_0114507, hsa_circ_0020460, hsa_circ_0053944, hsa_circ_0000758, hsa_circ_0000034, hsa_circ_0009127, hsa_circ_0052793, hsa_circ_0000301, hsa_circ_0000729. In addition, 9 (75%) of 12 differentially expressed circRNAs consist of exons (Table 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1061592/fimmu-13-1061592-HTML-r1/image_m/fimmu-13-1061592-t002.jpg)
Table 2 A total of 12 significantly differentially expressed peripheral blood circRNAs associated with aneurysm wall enhancement in UIA.
Upregulation of hsa_circ_0007990 in peripheral blood of UIA with AWE
We selected hsa_circ_0007990 to compare the expression levels among the three groups by qRT-PCR because hsa_circ_0007990 has the most remarkable fold change among the differentially expressed exonic circRNA associated with UIA with AWE. The expression of peripheral blood hsa_circ_0007990 was significantly increased in turn among the group of healthy control, UIA without AWE, and UIA with AWE (P<0.001; Figure 2). The level of hsa_circ_0007990 was higher in the UIA group than in the healthy control group (3.38 vs. 1.00, P=0.005). The level of hsa_circ_0007990 was significantly higher in the UIA with AWE than those without AWE (4.19 vs. 1.78, P=0.001).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1061592/fimmu-13-1061592-HTML-r1/image_m/fimmu-13-1061592-g002.jpg)
Figure 2 Upregulated expression of hsa_circ_0007990 associated with aneurysm wall enhancement confirmed by RT-PCR. The data are the mean ± SD, **p < 0.01, ***p < 0.001.
Characteristics and functions of hsa_circ_0007990
Hsa_circRNA_0007990 is derived from exon 8, exon 9, and exon 10 of the PGAP3 gene (NM_033419) located on homo sapiens chromosome 17 with a sequence length of 251 bp, and its circBank ID is hsa_circPGAP3_003, which contains the open reading frame (ORF), microRNA binding sites, and the ability of RNA binding protein (Figure 3). One function of circRNA is as microRNA sponges to control gene transcription by binding with microRNAs. A total of 29 microRNAs having at least two binding sites with hsa_circRNA_0007990 were predicted by the circBank database (Supplementary Table 2). The top ten predicted microRNAs are as follows, including miR-4717-5p, miR-1275, miR-150-3p, miR-18a-5p, miR-18b-5p, miR-3914, miR-3934-5p, miR-4254, miR-433-5p. Furthermore, four RNA binding proteins with hsa_circRNA_0007990 (argonaute-2 [Ago2], DiGeorge syndrome critical region 8 [DGCR8], eukaryotic initiation factor 4A-III [EIF4A3], polypyrimidine tract binding [PTB] protein) were found from the CircInteractome database (Figure 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1061592/fimmu-13-1061592-HTML-r1/image_m/fimmu-13-1061592-g003.jpg)
Figure 3 Schematic diagram of hsa_circ_0007990. The red bar represents the microRNA binding sites, the blue bar represents proteins binding with hsa_circRNA_0007990, and the light green bar represents the open reading frame.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1061592/fimmu-13-1061592-HTML-r1/image_m/fimmu-13-1061592-g004.jpg)
Figure 4 RNA binding proteins with hsa_circ_0007990. According to the CircInteractome database, the 3′start flanking intron upstream of hsa_circRNA_0007990 has a potential binding region with DGCR8 protein, and the 3′start flanking intron downstream of hsa_circRNA_0007990 has a potential binding region with Ago2, EIF4A3, and PTB protein.
We had the nucleotide sequence of hsa_circRNA_0007990 from the CircBase database, and then four ORFs were detected by the online database ORFfinder (https://www.ncbi.nlm.nih.gov/orffinder). Among the four ORFs (Table 3), only the ORF2 has a function of encoding protein named as the period circadian regulator 1 (PER1) superfamily searched from the online database of Conserved Domains (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi).
Gene oncology and KEGG pathway analysis of hsa_circ_0007990
We selected the top five predicted microRNAs binding on hsa_circRNA_0007990 through the CircBank database and searched the target mRNAs of these five microRNAs by both TargetScan and miRDB databases simultaneously (Supplementary Table 3). Then, we performed the GO and KEGG enrichment analysis to analyze the biological function of 961 target genes predicted by hsa_circRNA_0007990. The top ten enriched pathways from KEGG pathways analysis were found as follows: transcriptional misregulation in cancer, prostate cancer, Hippo signaling pathway, signaling pathways regulating pluripotency of stem cells, FoxO signaling pathway, Ras signaling pathway, fluid shear stress and atherosclerosis, focal adhesion, renal cell carcinoma, pathways in cancer (Figure 5A). The top ten enriched terms of biological process, cellular component, and molecular function of GO are shown in Figure 5B-D.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1061592/fimmu-13-1061592-HTML-r1/image_m/fimmu-13-1061592-g005.jpg)
Figure 5 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and Gene Ontology (GO) functional enrichment analysis of the predicted target genes from has_circRNA_0007990. (A) KEGG. (B) GO-molecular functions. (C) GO-biological processes. (D) GO-cell components.
CircRNA–microRNA–mRNA network construction
Ninety-seven targeted genes from the top ten enriched KEGG pathways were selected to construct the protein-protein interaction network by the STRING online tool (Figure 6A and Supplementary Table 4). Then the top ten hub genes were calculated by the MCC method in the CytoHubba plugin, including hypoxia-inducible factor 1 subunit alpha (HIF1A), estrogen receptor 1 (ESR1), forkhead box O1 (FOXO1), insulin-like growth factor 1 (IGF1), CREB binding protein (CREBBP), SMAD family member 2 (SMAD2), Kruppel like factor 4 (KLF4), ATM serine/threonine kinase (ATM), platelet-derived growth factor receptor beta (PDGFRB), protein tyrosine kinase 2 (PTK2) (Figure 6B and Supplementary Table 5). Finally, we obtained the hsa_circRNA_0007990-microRNA-mRNA network using the Cytoscape 3.9.1 software, including five microRNAs and 93 mRNAs containing ten key circRNA-microRNA-mRNA regulatory axes (Figure 6C). The top ten hub genes in the hsa_circRNA_0007990 ceRNA network were associated with inflammation by searching the previous studies (Supplementary Table 6).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1061592/fimmu-13-1061592-HTML-r1/image_m/fimmu-13-1061592-g006.jpg)
Figure 6 The circRNA-miRNA-mRNA regulatory network of hsa_circRNA_0007990. (A) Protein-protein interaction network obtained on STRING website. (B) Relationship network diagram of hub genes from protein-protein interaction network. (C) The reconstructed circRNA-microRNA-hub genes network.
Discussion
Our study showed 412 differentially expressed circRNAs in the peripheral blood between UIA patients and healthy controls by circRNA microarray. Twelve differentially expressed circRNAs were associated with AWE on HR-VMI in UIA patients. Hsa_circ_0007990 had the most remarkable fold change among the differentially expressed exonic circRNAs. The expression of hsa_circ_0007990 was increased gradually in the healthy control, UIA without AWE, and UIA with AWE. We predicted Ago2, DGCR8, EIF4A3, and PTB as the RNA binding proteins and period circadian regulator 1 as an encoding protein with hsa_circ_0007990. The hsa_circ_0007990-microRNA-mRNA network containing five microRNAs (miR-4717-5p, miR-1275, miR-150-3p, miR-18a-5p, miR-18b-5p), and 97 mRNAs was constructed. The ten hub genes, including HIF1A, ESR1, FOXO1, IGF1, CREBBP, SMAD2, KLF4, ATM, PDGFRB, and PTK2, were all associated with the inflammatory response. These results suggest that hsa_circ_0007990 may be a novel inflammatory blood biomarker for assessing UIA patients.
AWE on HR-VWI has been used as an imaging marker for vessel wall inflammation and aneurysm instability in patients with UIA. Recent studies showed the association of inflammatory blood biomarkers with AWE on HR-VWI in UIA patients (9–13). The plasma concentration of anti-inflammatory cytokine interleukin (IL)-10 in the aneurysm sac decreased in UIA with AWE (10). Our previous study showed that the neutrophil-lymphocyte ratio as an inflammatory peripheral blood biomarker was associated with AWE on HR-VWI in UIA patients (12). The up-regulated proinflammatory IL-1β and down-regulated anti-inflammatory IL-1RA in the serum were associated with AWE on HR-VWI in 34 patients with UIA. There was histological inflammation infiltration in aneurysm tissues with AWE. The serum IL-1β, IL-1RA, and IL-1RA/IL-1β ratio were associated with the inflammation-related proteins in aneurysm tissues (13). These results suggested that systematic inflammatory cytokines in the peripheral blood correlated with inflammatory factors in aneurysm tissues with AWE on HR-VWI, which may be a potential biomarker for aneurysm instability in UIA. In this study, we assessed the relationship of blood circRNAs as an inflammatory biomarker with AWE on HR-VWI in UIA patients.
CircRNAs play a key role in forming and rupturing IA. CircRNAs modulated phenotype, proliferation, migration, invasion, and apoptosis of vascular smooth muscle cells in aneurysm tissues (14–16). Inflammation is a vital pathological mechanism in IA rupture. Differentially expressed circRNAs related to IA have been shown to involve the inflammatory pathways and regulate the inflammatory gene expression through the circRNA-microRNA-mRNA axis (17–19). CircRNA_0079586 and circRNA_RanGAP1 in the endothelial cell of IA tissues were involved in the inflammatory response and IA rupture via miR-183-5p/myeloperoxidase and miR-877-3p/myeloperoxidase, respectively (17). CircLIFR in the IA tissues regulated the proliferation, migration, invasion, and apoptosis of vascular smooth muscle cells via the axis of miR-1299 and kinase insert domain receptor (18).
CircRNA hsa_circ_0007990 as a peripheral blood biomarker was associated with vessel wall inflammation on HR-VWI in this study. The ten hub genes in the predicted hsa_circ_0007990-microRNA-mRNA network were all involved in the inflammatory response (Supplementary Table 6). Recent studies suggested that circRNA may be used as a peripheral inflammatory blood biomarker for assessing IA instability. CircRNA expression profile of the peripheral blood in patients with multiple IA has been identified. The predicted dysregulated genes were mainly involved in the inflammatory response (21). Two hundred thirty-five differentially expressed circRNAs were shown in patients between ruptured IA and UIA by a circRNA microarray analysis of the peripheral blood. The mammalian target of rapamycin signaling pathway was one of the most significant pathways predicted by the targeted genes (22). Several differentially expressed circRNAs in IA tissues regulated inflammatory pathways. In the peripheral blood, two of these circRNAs (hsa_circ_0008433 and hsa_circ_0001946) were risk factors for IA rupture (19, 20).
In this study, the hsa_circ_0007990-microRNA-mRNA regulatory network consisted of five microRNAs and 97 mRNAs. CircRNAs play their functional effects via microRNA sponge activity to regulate gene transcription. Functional analysis indicated the top five predicted microRNAs (miR-18a-5p, miR-18b-5p, miR-150-3p, miR-1275, miR-4717-5p) binding on hsa_circRNA_0007990 through the CircBank database. MiR-18b-5p was shown as one of the significantly downregulated plasma microRNA biomarkers in ruptured IA, and its target genes were involved in inflammatory signaling pathways (26). Macrophage-derived miR-18a-5p was a potential urinary marker of early acute kidney injury following renal transplantation, and its several targets were involved in the cell cycle (27). MiR-150-3p was reported to involve in cyclophosphamide-induced lung inflammation and fibrosis via NF-κB and MAPK signaling pathways (28). In addition, five hub genes (HIF1A, ESR1, FOXO1, IGF1, and CREBBP) were identified in the hsa_circ_0007990-microRNA-mRNA regulatory network. HIF1A was upregulated in ruptured IAs and played an essential role in the PI3K-Akt signaling pathway (29, 30). FOXO1 increased in IA tissues and blocked myeloid cell leukemia sequence 1 transcription to involve the proliferation, apoptosis, and inflammation of vascular smooth muscle cells. MiR-124-5p downregulated FoxO1 expression to inhibit the migration and invasion of endothelial cells and inflammatory response (31, 32). IGF-1 alleviated the neuroinflammation and microglial activation induced by aneurysmal subarachnoid hemorrhage (33).
This study has some limitations. The sample size in this study was small with a single center data. The inflammatory response may promote intracranial aneurysm formation, progression, and rupture. In this study, the aneurysm size was significantly larger in patients with AWE than those without AWE (6.15 mm vs. 4.15 mm, P=0.025). Circ_0007990 expression was significantly higher in the UIA with AWE, which also tended to be associated with a larger size of the aneurysm. In one study, hsa_circ_0021001 expression in peripheral blood was not significantly associated with the size of the aneurysm (34). The relationship of differentially expressed circRNAs with a larger size of UIA needs to assess during aneurysm progression induced by inflammation. Three differentially expressed circRNAs (hsa_circ_0008433, hsa_circ_0001946, and hsa_circ_0021001) may be used as the diagnostic blood biomarkers for IA and the risk factors for aneurysms rupture (20, 34). The potential diagnostic value of 12 differentially expressed circRNAs identified from our study for IA diagnosis needs further assessment. The expression of predicted key microRNAs, hub mRNAs, and RNA binding proteins with hsa_circ_0007990 needs further confirmation in larger UIA samples. The exact mechanism of hsa_circ_0007990 involved in the inflammation pathway is still needed for further in vitro and in vivo study.
Conclusions
Differentially expressed blood circRNAs associated with AWE on HR-VWI can be used as novel inflammatory biomarkers for assessing UIA patients. The mechanism of hsa_circRNA_0007990 for UIA progression needs to investigate further.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The studies involving human participants were reviewed and approved by Sun Yat-sen Memorial Hospital, Sun Yat-sen University. The patients/participants provided their written informed consent to participate in this study.
Author contributions
Z-SS participated in the design of the study. Z-SS, X-BW and Y-TW wrote the initial manuscript. Z-SS revised the manuscript. All authors critically reviewed and edited the manuscript, approved the final version, participated in the interpretation and collection of the data. All authors contributed to the article and approved the submitted version.
Funding
Z-SS is supported by the National Natural Science Foundation of China (Grant Numbers: 81873752 and 81720108014) and the Science and Technology Planning Project of Guangzhou City (Grant Number: 201704020166).
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/fimmu.2022.1061592/full#supplementary-material
References
1. van der Kamp LT, Rinkel GJE, Verbaan D, van den Berg R, Vandertop WP, Murayama Y, et al. Risk of rupture after intracranial aneurysm growth. JAMA Neurol (2021) 78:1228–35. doi: 10.1001/jamaneurol.2021.2915
2. Etminan N, de Sousa DA, Tiseo C, Bourcier R, Desal H, Lindgren A, et al. European Stroke organisation (ESO) guidelines on management of unruptured intracranial aneurysms. Eur Stroke J (2022) 7:LXXXI-CVI. doi: 10.1177/23969873221099736
3. Tulamo R, Frösen J, Hernesniemi J, Niemelä M. Inflammatory changes in the aneurysm wall: a review. J Neurointerv Surg (2018) 10:i58–67. doi: 10.1136/jnis.2009.002055.rep
4. Shimonaga K, Matsushige T, Ishii D, Sakamoto S, Hosogai M, Kawasumi T, et al. Clinicopathological insights from vessel wall imaging of unruptured intracranial aneurysms. Stroke (2018) 49:2516–9. doi: 10.1161/STROKEAHA.118.021819
5. Samaniego EA, Roa JA, Hasan D. Vessel wall imaging in intracranial aneurysms. J neurointerv Surg (2019) 11:1105–12. doi: 10.1136/neurintsurg-2019-014938
6. Zhong W, Su W, Li T, Tan X, Chen C, Wang Q, et al. Aneurysm wall enhancement in unruptured intracranial aneurysms: a histopathological evaluation. J Am Heart Assoc (2021) 10:e018633. doi: 10.1161/JAHA.120.018633
7. Bae H, Suh SI, Yoon WK, Roh H, Kim C, Kwon TH. Correlation of aneurysmal wall enhancement of unruptured intracranial aneurysms on high-resolution vessel-wall imaging with clinical indices and surgical findings. Neurosurgery (2021) 89:420–7. doi: 10.1093/neuros/nyab178
8. Wu XB, Zhong JL, Wang SW, Su Y, Chen PS, Li ZJ, et al. Circumferential wall enhancement with contrast ratio measurement in unruptured intracranial aneurysm for aneurysm instability. Brain Behav (2022) 12:e2568. doi: 10.1002/brb3.2568
9. Tutino VM, Zebraski HR, Rajabzadeh-Oghaz H, Waqas M, Jarvis JN, Bach K, et al. Identification of circulating gene expression signatures of intracranial aneurysm in peripheral blood mononuclear cells. Diagn (Basel) (2021) 11:1092. doi: 10.3390/diagnostics11061092
10. Swiatek VM, Neyazi B, Roa JA, Zanaty M, Samaniego EA, Ishii D, et al. Aneurysm wall enhancement is associated with decreased intrasaccular IL-10 and morphological features of instability. Neurosurgery (2021) 89:664–71. doi: 10.1093/neuros/nyab279
11. Tutino VM, Lu Y, Ishii D, Poppenberg KE, Rajabzadeh-Oghaz H, Siddiqui AH, et al. Aberrant whole blood gene expression in the lumen of human intracranial aneurysms. Diagn (Basel) (2021) 11:1442. doi: 10.3390/diagnostics11081442
12. Wu XB, Zhong JL, Wang SW, Su Y, Chen PS, Li ZJ, et al. Neutrophil-to-lymphocyte ratio is associated with circumferential wall enhancement of unruptured intracranial aneurysm. Front Neurol (2022) 13:879882. doi: 10.3389/fneur.2022.879882
13. Liu Q, Zhang Y, Zhu C, Liu W, Ma X, Chen J, et al. Serum IL-1, pyroptosis and intracranial aneurysm wall enhancement: Analysis integrating radiology, serum cytokines and histology. Front Cardiovasc Med (2022) 9:818789. doi: 10.3389/fcvm.2022.818789
14. Chen X, Yang S, Yang J, Liu Q, Li M, Wu J, et al. Circular RNA circDUS2 is a potential biomarker for intracranial aneurysm. Front Aging Neurosci (2021) 13:632448. doi: 10.3389/fnagi.2021.632448
15. Chen X, Yang S, Yang J, Liu Q, Li M, Wu J, et al. The potential role of hsa_circ_0005505 in the rupture of human intracranial aneurysm. Front Mol Biosci (2021) 8:670691. doi: 10.3389/fmolb.2021.670691
16. Wang C, Luo Y, Tang H, Yan Y, Chang X, Zhao R, et al. Hsa_circ_0031608: a potential modulator of VSMC phenotype in the rupture of intracranial aneurysms. Front Mol Neurosci (2022) 15:842865. doi: 10.3389/fnmol.2022.842865
17. Zhang Z, Sui R, Ge L, Xia D. CircRNA_0079586 and circRNA_RanGAP1 are involved in the pathogenesis of intracranial aneurysms rupture by regulating the expression of MPO. Sci Rep (2021) 11:19800. doi: 10.1038/s41598-021-99062-w
18. Zhang H, Zhang B, Chen C, Chen J. Circular RNA circLIFR regulates the proliferation, migration, invasion and apoptosis of human vascular smooth muscle cells via the miR-1299/KDR axis. Metab Brain Dis (2022) 37:253–63. doi: 10.1007/s11011-021-00853-x
19. Huang Q, Huang QY, Sun Y, Wu S. High-throughput data reveals novel circular RNAs via competitive endogenous RNA networks associated with human intracranial aneurysms. Med Sci Monit (2019) 25:4819–30. doi: 10.12659/MSM.917081
20. Huang Q, Sun Y, Huang Q, Zeng Y, Lin S, Huang S, et al. Association between circular RNAs and intracranial aneurysm rupture under the synergistic effect of individual environmental factors. Front Neurol (2021) 12:594835. doi: 10.3389/fneur.2021.594835
21. Cao H, Chen J, Lai X, Liu T, Qiu P, Que S, et al. Circular RNA expression profile in human primary multiple intracranial aneurysm. Exp Ther Med (2021) 21:239. doi: 10.3892/etm.2021.9670
22. Ma Y, Zhang B, Zhang D, Wang S, Li M, Zhao J. Differentially expressed circular RNA profile in an intracranial aneurysm group compared with a healthy control group. Dis Markers (2021) 2021:8889569. doi: 10.1155/2021/8889569
23. Chomczynski P, Wilfinger WW, Eghbalnia HR, Kennedy A, Rymaszewski M, Mackey K. Inter-individual differences in RNA levels in human peripheral blood. PloS One (2016) 11:e0148260. doi: 10.1371/journal.pone.0148260
24. Korthauer K, Kimes PK, Duvallet C, Reyes A, Subramanian A, Teng M, et al. A practical guide to methods controlling false discoveries in computational biology. Genome Biol (2019) 20:118. doi: 10.1186/s13059-019-1716-1
25. Liu Z, Tuo YH, Chen JW, Wang QY, Li S, Li MC, et al. NADPH oxidase inhibitor regulates microRNAs with improved outcome after mechanical reperfusion. J neurointerv Surg (2017) 9:702–6. doi: 10.1136/neurintsurg-2016-012463
26. Supriya M, Christopher R, Indira Devi B, Bhat DI, Shukla D. Circulating microRNAs as potential molecular biomarkers for intracranial aneurysmal rupture. Mol Diagn Ther (2020) 24:351–64. doi: 10.1007/s40291-020-00465-8
27. Connor KL, Teenan O, Cairns C, Banwell V, Thomas RA, Rodor J, et al. Identifying cell-enriched miRNAs in kidney injury and repair. JCI Insight (2020) 5:e140399. doi: 10.1172/jci.insight.140399
28. Abdel-Latif GA, Elwahab AHA, Hasan RA, ElMongy NF, Ramzy MM, Louka ML, et al. A novel protective role of sacubitril/valsartan in cyclophosphamide induced lung injury in rats: impact of miRNA-150-3p on NF-κB/MAPK signaling trajectories. Sci Rep (2020) 10:13045. doi: 10.1038/s41598-020-69810-5
29. Kurki MI, Häkkinen SK, Frösen J, Tulamo R, von und zu Fraunberg M, Wong G, et al. Upregulated signaling pathways in ruptured human saccular intracranial aneurysm wall: an emerging regulative role of toll-like receptor signaling and nuclear factor-κB, hypoxia-inducible factor-1A, and ETS transcription factors. Neurosurgery (2011) 68:1667–75. doi: 10.1227/NEU.0b013e318210f001
30. Pan YB, Lu J, Yang B, Lenahan C, Zhang J, Shao A. Construction of competitive endogenous RNA network reveals regulatory role of long non-coding RNAs in intracranial aneurysm. BMC Neurosci (2021) 22:15. doi: 10.1186/s12868-021-00622-7
31. Huang J, Hong L, Shen B, Zhou Y, Lan J, Peng Y. FOXO1 represses MCL1 transcription to regulate the function of vascular smooth muscle cells in intracranial aneurysm. Exp Brain Res (2022) 240:2861–70. doi: 10.1007/s00221-022-06461-0
32. Wang RK, Sun YY, Li GY, Yang HT, Liu XJ, Li KF, et al. MicroRNA-124-5p delays the progression of cerebral aneurysm by regulating FoxO1. Exp Ther Med (2021) 22:1172. doi: 10.3892/etm.2021.10606
33. Chen TF, Chen KW, Chien Y, Lai YH, Hsieh ST, Ma HY, et al. Dental pulp stem cell-derived factors alleviate subarachnoid hemorrhage-induced neuroinflammation and ischemic neurological deficits. Int J Mol Sci (2019) 20:3747. doi: 10.3390/ijms20153747
Keywords: aneurysm wall enhancement, blood marker, circular RNA, intracranial aneurysm, vessel wall imaging, inflammation
Citation: Wu X-B, Wu Y-T, Guo X-X, Xiang C, Chen P-S, Qin W and Shi Z-S (2022) Circular RNA hsa_circ_0007990 as a blood biomarker for unruptured intracranial aneurysm with aneurysm wall enhancement. Front. Immunol. 13:1061592. doi: 10.3389/fimmu.2022.1061592
Received: 04 October 2022; Accepted: 04 November 2022;
Published: 17 November 2022.
Edited by:
Yingmei Feng, Beijing Youan Hospital, Capital Medical University, ChinaReviewed by:
Haocheng Lu, Southern University of Science and Technology, ChinaLixin Wang, Fudan University, China
Copyright © 2022 Wu, Wu, Guo, Xiang, Chen, Qin and Shi. 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: Zhong-Song Shi, shizhs@mail.sysu.edu.cn
†These authors have contributed equally to this work