RNA sequencing analysis between ruptured and un-ruptured brain AVM
Chinese Neurosurgical Journal volume 8, Article number: 13 (2022)
A brain arteriovenous malformation (BAVM) is a tangle of abnormal blood vessels connecting the arteries and veins in the brain and is associated with a higher risk for intracerebral hemorrhage (ICH). RNA sequencing technology has been recently used to investigate the mechanism of diseases owing to its ability to identify the gene changes on a transcriptome-wide level. This study aims to gain insights into the potential mechanism involved in BAVM rupture.
Sixty-five BAVM nidus samples were collected, among which 28 were ruptured and 37 were un-ruptured. Then, next-generation RNA sequencing was performed on all of them to obtain differential expressed genes (DEGs) between the two groups. In addition, bioinformatics analysis was performed to evaluate the involved biological processes and pathways by GO and KEGG analysis. Finally, we performed a univariate Cox regression analysis to obtain the early rupture-prone DEGs.
A total of 951 genes were differentially expressed between the ruptured and un-ruptured BAVM groups, of which 740 genes were upregulated and 211 genes were downregulated in ruptured BAVMs. Then, bioinformatics analysis showed the biological processes and pathways related to the inflammatory processes and extracellular matrix organization were significantly enriched. Meanwhile, some downregulated genes are involved in cell adhesion and genes participating in response to muscle activity and the terms of nervous system development. Finally, one hundred twenty-five genes, many were involved in inflammation, were correlated with the early rupture of BAVMs.
The upregulated genes in the ruptured BAVM group were involved in inflammatory processes and extracellular matrix organization. Some of the downregulated genes participated in cell adhesion and myofibril assembly, indicating the role of enhanced inflammation and reduced inflammation vessel strength in BAVMs rupture.
Brain arteriovenous malformations (BAVMs) are a complex tangle of abnormal arteries and veins located in the brain, and this vascular abnormality lacks an intervening capillary network . The risk of spontaneous intracranial hemorrhage (ICH) caused by the rupture of BAVM was estimated at 2–4% per year and can be as high as 34% under certain conditions, and the ICH could cause catastrophic neurological consequences [2, 3]. Early surgical management can effectively eliminate the cumulative lifetime risk of ICH; however, the higher rate of permanent neurological defects or death makes the treatment of un-ruptured BAVMs increasingly controversial . Currently, there is no specific available medical therapy for the treatment of BAVMs. The investigation of the mechanism of BAVM rupture would shed light on developing new medical therapies, but the underlying mechanism of this process remained largely elusive.
It is likely that many cellular and molecular events played orchestrally, which eventually led to the rupture of BAVM . Instead of examining individual or a few interrelated gene products, RNA sequencing technology enables the investigation of the entire transcriptome simultaneously and profiling the gene signature comprehensively. Recent RNA sequencing studies have been used to explore the pathogenesis and development of BAVM [6, 7]. However, no such studies have looked into the rupture processes of BAVM.
Thus, in this study, we classified the relatively large number of BAVM samples into the un-ruptured and ruptured groups and performed RNA sequencing analysis of these samples to unveil the underlying expression signatures along with the rupture of BAVM.
Patients and samples
All study participants with a defined nidus and evidence of arteriovenous shunting on digital subtraction angiography were recruited and underwent microsurgical resection in the Neurosurgical Department of Beijing Tiantan Hospital. A group of 65 BAVMs nidus was collected after the resection. Then, 31 BAVM nidus were obtained as an independent validation group. All the specimens were snap-frozen and stored in liquid nitrogen. All the participants provided written informed consent, and the institutional review board of Beijing Tiantan Hospital approved this study, Capital Medical University. BAVM patients were dichotomized as ruptured and un-ruptured groups based on the radiological examination (CT or MRI) combined with their medical history. The location of BAVM was defined as the brain lobe to which most part of the BAVM belonged.
RNA extraction, library preparation, and sequencing
Total RNA was extracted from the BAVM nidus using TRIzol Reagent according to the manufacturer’s instructions. Then, RNA purity was determined by NanoPhotometer® spectrophotometer (IMPLEN, Germany). The concentration and integrity were examined by Qubit® RNA Assay Kit with a Qubit® 2.0 Fluorometer (Life Technologies, USA) and RNA Nano 6000 Assay Kit in a Bioanalyzer 2100 system (Agilent Technologies, USA). The detailed method, RNA sequencing data, and RNA sequencing validation using RT-qPCR could be found in our previous studies [7,8,9].
Differential expressed genes (DEGs) and bioinformatics analysis
To obtain the DEGs between the ruptured and un-ruptured BAVM groups, reads containing adaptors or poly-N and low-quality reads were first removed from the raw data. Then, the clean reads were aligned with both the transcript reference and Hg19 RefSeq (RNA sequences, GRCh37) using STAR v2.2.1. The gene expression was calculated as fragments per kilobase of exon per million mapped fragment values with RSEM v1.3.0. DEGs were defined as differential expressed transcripts with P < 0.05 and fold change > 1.50 or < 0.67.
GO and KEGG analysis
For further analyzing the common DEGs, Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted on DAVID (The Database for Annotation, Visualization and Integrated Discovery, https://david.ncifcrf.gov/). An FDR value of 0.05 was set as the cutoff criterion. Gene expression data were analyzed in R v3.5.1.
Cox regression analysis
The univariate Cox regression analysis was applied to select for BAVM early rupture-prone prognostic DEGs. DEGs associated with overall survival (OS) with P < 0.05 were considered statistically significant. The univariate analysis was performed using the R packages “survival” and “surveminer” to identify OS-related DEGs.
The baseline clinical characteristics of the 65 BAVM patients were summarized in Table 1, among which 28 were ruptured BAVMs and 37 were un-ruptured BAVMs. The typical CT\MRI images of un-ruptured and ruptured BAVM were present in Fig. 1.
All of these patients were of Chinese Han ethnicity and had no BAVM family history. There is no significant difference in baseline characteristics between these two groups (Additional file 1: Table S1).
Differential expressed genes analysis
A total of 951 genes were differentially expressed between the ruptured and un-ruptured BAVM groups (P < 0.05; fold change > 1.5 or < 0.67), of which 740 (77.8%) were upregulated and 211 (22.1%) genes were downregulated in ruptured BAVMs (Fig. 2). The top 5 upregulated and downregulated DEGs (based on the fold changes) were summarized in Table 2. Notably, all of the 5 upregulated DEGs (HMOX1, PLA2G7, FUCA1, ACP5, and IFI30) were related to an inflammatory process or oxidative stress injury. Meanwhile, among the top 5 downregulated genes, TNFSF18 and CCL4L1 belong to cytokines; TPH1 encodes a member of the aromatic amino acid hydroxylase family, which plays crucial roles in the biosynthesis of serotonin ; PKHD1L1 could encode a receptor with inducible T lymphocyte expression , and TMEM235 encodes a membrane protein which is over-expressed in the brain.
Functional analysis of the DEGs
Functional analysis of the upregulated genes in ruptured BAVM showed the biological processes and pathways related to the inflammatory processes were significantly enriched (Fig. 3, Additional files 2 and 3: Tables S2 and S3). GO analysis demonstrated that several inflammatory terms, including innate immune response and chemotaxis, were highly enriched. In addition, GO terms about extracellular matrix organization were also highly enriched, comprising integrin-mediated signaling pathway and collagen catabolic process, indicating their essential roles in BAVM rupture. KEGG analysis showed inflammation-related diseases and signaling pathways, such as systemic lupus erythematosus, Toll-like receptor signaling pathway, NF-κB signaling pathway, and complement and coagulation cascades which were enriched (FDR-adjusted P < 0.05). Meanwhile, some downregulated genes were involved in cell adhesion (DSCAM, CNTN2, CERCAM, ITGBL1, and LAMA3) and genes participating in myofibril assembly (CAPN3, MYOZ1, and MYOZ2). GO and KEGG analysis showed the downregulated genes were also enriched in nervous system development, possibly indicating an important parenchymal impact of BAVM rupture (Additional files 4 and 5: Tables S4 and S5).
Early rupture-prone DEGs by univariate Cox regression analysis
One hundred twenty-five genes were correlated with the early rupture of BAVMs, among which 96 were upregulated and 29 were downregulated. GO analysis showed that inflammatory response and innate immune response were significantly enriched (Additional file 6: Table S6). There were no KEGG terms from upregulating genes, and GO and KEGG terms from downregulating genes were significantly enriched.
In the ruptured BAVM group, a large cohort of upregulated genes was involved in inflammatory processes, including recruited leukocytes, enhanced cytokines, and leukocyte migration, contributing to vessel wall damage and weakening . Notably, the innate immune responses were highly upregulated in ruptured BAVMs, indicating the crucial role of the innate immune system in BAVMs rupture processes. The cell components of the innate immune system mainly include monocytes, macrophages, and highly phagocytic neutrophils , among which macrophage load in BAVM nidus was correlated to higher hemorrhage risk by using ferumoxytol-enhanced MRI . Another study found a higher density of M2 macrophages in the perinidal dilated capillary network (PDCN) in the ruptured BAVM, and it postulated M2 macrophages could lead to an inappropriate BBB which was prone to bleed . Toll-like receptors (TLRs) have crucial roles in activating the innate immune system . As a member of the TLR family, TLR4 is mainly expressed on the cell surface of inflammatory cells, endothelial cells, and smooth muscle cells [17, 18]. TLR4 can bind to MyD88, and the binding between these two proteins leads to downstream expression of pro-inflammatory cytokines . It has been shown that the pathway involving TLR4 and MyD88 promotes tissue remodeling and rupture of a brain aneurysm . In this study, the TLR4 signaling pathway and MyD88-dependent Toll-like receptor signaling pathway were upregulated in ruptured BAVMs. Except for innate immune responses, several other inflammation-related pathways were also increased in the ruptured BAVM group. For instance, NF-κB is a transcription factor that could functionally regulate several pro-inflammatory genes. NF-κB p65, a subunit of NF-κB, was found at a higher level in the ruptured BAVM group compared to the un-ruptured BAVM group . Consistently, in our study, the NF-κB signaling pathway was upregulated in ruptured BAVMs. Several DEGs identified in this study were also tested in the peripheral circulation blood . For instance, TNFAIP6 were found both in ruptured BAVM tissue and peripheral blood. This protein plays a crucial role in the protease network associated with inflammation and is also involved in extracellular matrix (ECM) remodeling, cell adhesion, and cell migration . It is indicated that these DEGs had the potential as biomarkers for the BAVM rupture.
Moreover, extracellular matrix organization was also activated in ruptured BAVMs. The genes involved in this process include COL15A1, COL11A1, COL1A1, COL3A1, COL1A2, COL5A1, COL5A2, COL6A1, and COL6A3, all of which encode subtypes of collagen that constitute the ECM basement membrane and interstitial matrix . Matrix metalloproteinases (MMPs) can degrade all kinds of extracellular matrix proteins and can process a number of bioactive molecules. MMP9 has been reported causing damage to the BAVM vessel walls, leading to rupture of the nidus . In our study, MMP2, instead of MMP9, was increased in ruptured BAVMs, possibly indicating its important role in BAVM rupture. In addition, several genes encoding the activities of the central nervous system, such as MOG, ZBTB16, CNTN2, and MAL, were downregulated in ruptured BAVMs. Due to the abnormal vasculature of the blood vessels in BAVM, the perinidal brain parenchyma, which mainly includes the glial cells, astrocytes, and neurons, is highly affected . We hypothesized that destabilizing BAVM which had increased hemorrhage risk would affect more on the brain parenchyma and result in abnormal gene expressions and cell functions.
In conclusion, the upregulated genes in the ruptured BAVM group were involved in inflammatory processes and extracellular matrix organization, while some of the downregulated genes participated in cell adhesion and myofibril assembly, indicating the role of enhanced inflammation and reduced vessel strength in BAVMs rupture.
Availability of data and materials
Please contact the authors for data requests.
Brain arteriovenous malformation
Differential expressed genes
Kyoto Encyclopedia of Genes and Genomes
Horton JC. Arteriovenous malformations of the brain. N Engl J Med. 2007;357(17):1774–5 author reply 5.
Kim H, Sidney S, McCulloch CE, Poon KY, Singh V, Johnston SC, et al. Racial/ethnic differences in longitudinal risk of intracranial hemorrhage in brain arteriovenous malformation patients. Stroke. 2007;38(9):2430–7.
van Beijnum J, van der Worp HB, Buis DR, Al-Shahi Salman R, Kappelle LJ, Rinkel GJ, et al. Treatment of brain arteriovenous malformations: a systematic review and meta-analysis. JAMA. 2011;306(18):2011–9.
Davidson AS, Morgan MK. How safe is arteriovenous malformation surgery? A prospective, observational study of surgery as first-line treatment for brain arteriovenous malformations. Neurosurgery. 2010;66(3):498–504 discussion -5.
Kim H, Pawlikowska L, Chen Y, Su H, Yang GY, Young WL. Brain arteriovenous malformation biology relevant to hemorrhage and implication for therapeutic development. Stroke. 2009;40(3 Suppl):S95–7.
Hauer AJ, Kleinloog R, Giuliani F, Rinkel GJE, de Kort GA, Berkelbach van der Sprenkel JW, et al. RNA-sequencing highlights inflammation and impaired integrity of the vascular wall in brain arteriovenous malformations. Stroke. 2020;51(1):268–74.
Fu W, Huo R, Yan Z, Xu H, Li H, Jiao Y, et al. Mesenchymal behavior of the endothelium promoted by Smad6 downregulation is associated with brain arteriovenous malformation microhemorrhage. Stroke. 2020;51(7):2197–207.
Li H, Nam Y, Huo R, Fu W, Jiang B, Zhou Q, et al. De novo germline and somatic variants convergently promote endothelial-to-mesenchymal transition in simplex brain arteriovenous malformation. Circ Res. 2021;129(9):825–39.
Huo R, Fu W, Li H, Jiao Y, Yan Z, Wang L, et al. RNA sequencing reveals the activation of Wnt signaling in low flow rate brain arteriovenous malformations. J Am Heart Assoc. 2019;8(12):e012746.
Jin B, Ha SE, Wei L, Singh R, Zogg H, Clemmensen B, et al. Colonic motility is improved by the activation of 5-Ht(2b) receptors on interstitial cells of Cajal in diabetic mice. Gastroenterology. 2021;161(2):608–22.e7.
Gekara NO. DNA damage-induced immune response: micronuclei provide key platform. J Cell Biol. 2017;216(10):2999–3001.
Mouchtouris N, Jabbour PM, Starke RM, Hasan DM, Zanaty M, Theofanis T, et al. Biology of cerebral arteriovenous malformations with a focus on inflammation. J Cereb Blood Flow Metab. 2015;35(2):167–75.
Kaur BP, Secord E. Innate immunity. Immunol Allergy Clin N Am. 2021;41(4):535–41.
Zanaty M, Chalouhi N, Starke RM, Jabbour P, Hasan D. Molecular imaging in neurovascular diseases: the use of ferumoxytol to assess cerebral aneurysms and arteriovenous malformations. Top Magn Reson Imaging. 2016;25(2):57–61.
Nakamura Y, Sugita Y, Nakashima S, Okada Y, Yoshitomi M, Kimura Y, et al. Alternatively activated macrophages play an important role in vascular remodeling and hemorrhaging in patients with brain arteriovenous malformation. J Stroke Cerebrovasc Dis. 2016;25(3):600–9.
Beutler BA. TLRs and innate immunity. Blood. 2009;113(7):1399–407.
Akira S, Takeda K. Toll-like receptor signalling. Nat Rev Immunol. 2004;4(7):499–511.
Andonegui G, Zhou H, Bullard D, Kelly MM, Mullaly SC, McDonald B, et al. Mice that exclusively express Tlr4 on endothelial cells can efficiently clear a lethal systemic Gram-negative bacterial infection. J Clin Investig. 2009;119(7):1921–30.
Fang H, Wang PF, Zhou Y, Wang YC, Yang QW. Toll-like receptor 4 signaling in intracerebral hemorrhage-induced inflammation and injury. J Neuroinflammation. 2013;10:27.
Mitsui K, Ikedo T, Kamio Y, Furukawa H, Lawton MT, Hashimoto T. Tlr4 (Toll-like receptor 4) mediates the development of intracranial aneurysm rupture. Hypertension (Dallas, Tex: 1979). 2020;75(2):468–76.
Aziz MM, Takagi Y, Hashimoto N, Miyamoto S. Activation of nuclear factor κB in cerebral arteriovenous malformations. Neurosurgery. 2010;67(6):1669–79 discussion 79-80.
Weinsheimer SM, Xu H, Achrol AS, Stamova B, McCulloch CE, Pawlikowska L, et al. Gene expression profiling of blood in brain arteriovenous malformation patients. Transl Stroke Res. 2011;2(4):575–87.
Chan TC, Li CF, Ke HL, Wei YC, Shiue YL, Li CC, et al. High Tnfaip6 level is associated with poor prognosis of urothelial carcinomas. Urol Oncol. 2019;37(4):293.e11–24.
Manon-Jensen T, Kjeld NG, Karsdal MA. Collagen-mediated hemostasis. J Thromb Haemost. 2016;14(3):438–48.
Li X, Wang R, Wang X, Xue X, Ran D, Wang S. Relevance of Il-6 and Mmp-9 to cerebral arteriovenous malformation and hemorrhage. Mol Med Rep. 2013;7(4):1261–6.
Ng I, Tan WL, Ng PY, Lim J. Hypoxia inducible factor-1alpha and expression of vascular endothelial growth factor and its receptors in cerebral arteriovenous malformations. J Clin Neurosci. 2005;12(7):794–9.
First, neuroinflammation might exhibit a few hours after intracranial hemorrhage, which might partly impact the DEGs involved in inflammation. In this study, we removed the surrounding tissue of BAVM nidus as clean as possible to reduce the above impact, but there were still remnants. Second, further confirmatory studies, including RT-qPCR and immunofluorescent staining, are necessary to validate the DEGs. Third, multivariate Cox regression analysis was further needed to achieve high-risk factors for BAVM rupture when the sample size is large enough.
This study was supported by the project Genomics Platform Construction for Chinese Major Brain Disease-AVM (grant No. PXM2019_026280_000002-AVM, principal investigator: Dr. Cao).
Ethics approval and consent to participate
This study was approved by the institutional review board of Beijing Tiantan Hospital, Capital Medical University (KY2017-035-02).
Consent for publication
All the participants provided written informed consent.
The authors declare that they have no competing interests.
comparison of baseline between ruptured and un-ruptured BAVM group.
GO terms enriched by the up-regulated genes in ruptured BAVMs.
KEGG terms enriched by the up-regulated genes in ruptured BAVMs.
GO terms enriched by the down-regulated genes in ruptured BAVMs.
KEGG terms enriched by the down-regulated genes in ruptured BAVMs.
GO terms enriched by the up-regulated BAVM early rupture-prone genes.
The typical CT\MRI images of un-ruptured and ruptured BAVM.
About this article
Cite this article
Li, H., Yan, Z., Huo, R. et al. RNA sequencing analysis between ruptured and un-ruptured brain AVM. Chin Neurosurg Jl 8, 13 (2022). https://doi.org/10.1186/s41016-022-00282-4
- Brain arteriovenous malformation
- Inflammatory processes
- Cell adhesion
- Myofibril assembly