Whole blood transcriptome biomarkers of unruptured intracranial aneurysm

Modulation of the Immunological Milieu in Acute Aneurysmal Subarachnoid Hemorrhage: The Potential Role of Monocytes Through CXCL10 Secretion

Emerging evidence indicates that aneurysmal subarachnoid hemorrhage (aSAH) elicits a response from both innate and adaptive immune systems. An upregulation of CD8 + CD161 + cells has been observed in the cerebrospinal fluid (CSF) after aSAH, yet the precise role of these cells in the context of aSAH is unkown. CSF samples from patients with aSAH and non-aneurysmal SAH (naSAH) were analyzed. Single-cell RNA sequencing (scRNAseq) was performed on CD8 + CD161 + sorted samples from aSAH patients. Cell populations were identified using "clustering." Gene expression levels of ten previously described genes involved in inflammation were quantified from aSAH and naSAH samples using RT-qPCR. The study focused on the following genes: CCL5, CCL7, APOE, SPP1, CXCL8, CXCL10, HMOX1, LTB, MAL, and HLA-DRB1. Gene clustering analysis revealed that monocytes, NK cells, and T cells expressed CD8 + CD161 + in the CSF of patients with aSAH. In comparison to naSAH samples, aSAH samples exhibited higher mRNA levels of CXCL10 (median, IQR = 90, 16-149 vs. 0.5, 0-6.75, p = 0.02). A trend towards higher HMOX1 levels was also observed in aSAH (median, IQR = 12.6, 9-17.6 vs. 2.55, 1.68-5.7, p = 0.076). Specifically, CXCL10 and HMOX1 were expressed by the monocyte subpopulation. Monocytes, NK cells, and T cells can potentially express CD8 + CD161 + in patients with aSAH. Notably, monocytes show high levels of CXCL10. The elevated expression of CXCL10 in aSAH compared to naSAH indicates its potential significance as a target for future studies.

Dysregulated Genes and Signaling Pathways in the Formation and Rupture of Intracranial Aneurysm

Intracranial aneurysm (IA) has the potential to rupture. Despite scientific advances, we are still not in a position to screen patients for IA and identify those at risk of rupture. It is critical to comprehend the molecular basis of disease to facilitate the development of novel diagnostic strategies. We used transcriptomics to identify the dysregulated genes and understand their role in the disease biology. In particular, RNA-Seq was performed in tissue samples of controls, unruptured IA, and ruptured IA. Dysregulated genes (DGs) were identified and analyzed to understand the functional aspects of molecules. Subsequently, candidate genes were validated at both transcript and protein level. There were 314 DGs in patients with unruptured IA when compared to control samples. Out of these, SPARC and OSM were validated as candidate molecules in unruptured IA. PI3K-AKT signaling pathway was found to be an important pathway for the formation of IA. Similarly, 301 DGs were identified in the samples of ruptured IA when compared with unruptured IAs. CTSL was found to be a key candidate molecule which along with Hippo signaling pathway may be involved in the rupture of IA. We conclude that activation of PI3K-AKT signaling pathway by OSM along with up-regulation of SPARC is important for the formation of IA. Further, regulation of Hippo pathway through PI3K-AKT signaling results in the down-regulation of YAP1 gene. This along with up-regulation of CTSL leads to further weakening of aneurysm wall and its subsequent rupture.

Diagnosis of Unruptured Intracranial Aneurysm by High-Performance Serum Metabolic Fingerprints

Unruptured intracranial aneurysm (UIA) is a high-risk cerebrovascular saccular dilatation, the effective medical management of which depends on high-performance diagnosis. However, most UIAs are diagnosed incidentally during neurovascular imaging modalities, which are time-consuming and harmful (e.g., radiation). Serum metabolic fingerprints is a promising alternative for early diagnosis of UIA. Here, nanoparticle enhanced laser desorption/ionization mass spectrometry is applied to obtain high-performance UIA-specific serum metabolic fingerprints. Diagnostic performance with an area-under-the-curve (AUC) of 0.842 (95% confidence interval (CI): 0.783-0.891) is achieved by the constructed machine learning (ML) model, including ML algorithm selection and feature selection. Lactate, glutamine, homoarginine, and 3-methylglutaconic acid are identified as the metabolic biomarker panel, which showed satisfactory diagnosis (AUC of 0.812, 95% CI: 0.727-0.897) and effective growth risk assessment (p<0.05, two-tailed t-test) of UIAs. This work aims to promote the diagnostics of UIAs and metabolic biomarker screening for medical management.

Transcriptomic Studies on Intracranial Aneurysms

Intracranial aneurysm (IA) is a relatively common vascular malformation of an intracranial artery. In most cases, its presence is asymptomatic, but IA rupture causing subarachnoid hemorrhage is a life-threating condition with very high mortality and disability rates. Despite intensive studies, molecular mechanisms underlying the pathophysiology of IA formation, growth, and rupture remain poorly understood. There are no specific biomarkers of IA presence or rupture. Analysis of expression of mRNA and other RNA types offers a deeper insight into IA pathobiology. Here, we present results of published human studies on IA-focused transcriptomics.

A machine learning-derived gene signature for assessing rupture risk and circulatory immunopathologic landscape in patients with intracranial aneurysms

Background

Intracranial aneurysm (IA) is an uncommon but severe subtype of cerebrovascular disease, with high mortality after aneurysm rupture. Current risk assessments are mainly based on clinical and imaging data. This study aimed to develop a molecular assay tool for optimizing the IA risk monitoring system.

Methods

Peripheral blood gene expression datasets obtained from the Gene Expression Omnibus were integrated into a discovery cohort. Weighted gene co-expression network analysis (WGCNA) and machine learning integrative approaches were utilized to construct a risk signature. QRT-PCR assay was performed to validate the model in an in-house cohort. Immunopathological features were estimated using bioinformatics methods.

Results

A four-gene machine learning-derived gene signature (MLDGS) was constructed for identifying patients with IA rupture. The AUC of MLDGS was 1.00 and 0.88 in discovery and validation cohorts, respectively. Calibration curve and decision curve analysis also confirmed the good performance of the MLDGS model. MLDGS was remarkably correlated with the circulating immunopathologic landscape. Higher MLDGS scores may represent higher abundance of innate immune cells, lower abundance of adaptive immune cells, and worse vascular stability.

Conclusions

The MLDGS provides a promising molecular assay panel for identifying patients with adverse immunopathological features and high risk of aneurysm rupture, contributing to advances in IA precision medicine.

RNA Expression Signatures of Intracranial Aneurysm Growth Trajectory Identified in Circulating Whole Blood

After detection, identifying which intracranial aneurysms (IAs) will rupture is imperative. We hypothesized that RNA expression in circulating blood reflects IA growth rate as a surrogate of instability and rupture risk. To this end, we performed RNA sequencing on 66 blood samples from IA patients, for which we also calculated the predicted aneurysm trajectory (PAT), a metric quantifying an IA's future growth rate. We dichotomized dataset using the median PAT score into IAs that were either more stable and more likely to grow quickly. The dataset was then randomly divided into training (n = 46) and testing cohorts (n = 20). In training, differentially expressed protein-coding genes were identified as those with expression (TPM > 0.5) in at least 50% of the samples, a q-value < 0.05 (based on modified F-statistics with Benjamini-Hochberg correction), and an absolute fold-change ≥ 1.5. Ingenuity Pathway Analysis was used to construct networks of gene associations and to perform ontology term enrichment analysis. The MATLAB Classification Learner was then employed to assess modeling capability of the differentially expressed genes, using a 5-fold cross validation in training. Finally, the model was applied to the withheld, independent testing cohort (n = 20) to assess its predictive ability. In all, we examined transcriptomes of 66 IA patients, of which 33 IAs were "growing" (PAT ≥ 4.6) and 33 were more "stable". After dividing dataset into training and testing, we identified 39 genes in training as differentially expressed (11 with decreased expression in "growing" and 28 with increased expression). Model genes largely reflected organismal injury and abnormalities and cell to cell signaling and interaction. Preliminary modeling using a subspace discriminant ensemble model achieved a training AUC of 0.85 and a testing AUC of 0.86. In conclusion, transcriptomic expression in circulating blood indeed can distinguish "growing" and "stable" IA cases. The predictive model constructed from these differentially expressed genes could be used to assess IA stability and rupture potential.

Profiling of Circulating Gene Expression Reveals Molecular Signatures Associated with Intracranial Aneurysm Rupture Risk

BACKGROUND

Following detection, rupture risk assessment for intracranial aneurysms (IAs) is critical. Towards molecular prognostics, we hypothesized that circulating blood RNA expression profiles are associated with IA risk.

METHODS

We performed RNA sequencing on 68 blood samples from IA patients. Here, patients were categorized as either high or low risk by assessment of aneurysm size (≥ 5 mm = high risk) and Population, Hypertension, Age, Size, Earlier subarachnoid hemorrhage, Site (PHASES) score (≥ 1 = high risk). Modified F-statistics and Benjamini-Hochberg false discovery rate correction was performed on transcripts per million-normalized gene counts. Protein-coding genes expressed in ≥ 50% of samples with a q value < 0.05 and an absolute fold-change ≥ 2 were considered significantly differentially expressed. Bioinformatics in Ingenuity Pathway Analysis was performed to understand the biology of risk-associated expression profiles. Association was assessed between gene expression and risk via Pearson correlation analysis. Linear discriminant analysis models using significant genes were created and validated for classification of high-risk cases.

RESULTS

We analyzed transcriptomes of 68 IA patients. In these cases, 31 IAs were large (≥ 5 mm), while 26 IAs had a high PHASES score. Based on size, 36 genes associated with high-risk IAs, and two were correlated with the size measurement. Alternatively, based on PHASES score, 76 genes associated with high-risk cases, and nine of them showed significant correlation to the score. Similar ontological terms were associated with both gene profiles, which reflected inflammatory signaling and vascular remodeling. Prediction models based on size and PHASES stratification were able to correctly predict IA risk status, with > 80% testing accuracy for both.

CONCLUSIONS

Here, we identified genes associated with IA risk, as quantified by common clinical metrics. Preliminary classification models demonstrated feasibility of assessing IA risk using whole blood expression.

Canonical NF-κB signaling pathway and GRO-α/CXCR2 axis are activated in unruptured intracranial aneurysm patients

Activation of the nuclear factor kappa-B (NF-κB) stimulates the production of pro-inflammatory molecules involved in the formation of intracranial aneurysms (IA). The study aimed to assess the NF-κB p65 subunit and the GRO-α chemokine and its receptor CXCR2 concentrations in unruptured intracranial aneurysm patients (UIA, n = 25) compared to individuals without vascular changes in the brain (n = 10). It was also analyzed whether tested proteins are related to the size and number of aneurysms. Cerebrospinal fluid (CSF) and serum protein levels were measured using the ELISA method. Median CSF and serum NF-κB p65 concentrations were significantly lower, while median CSF GRO-α and CXCR2 concentrations were significantly higher in UIA patients compared to the control group. CSF and serum NF-κB p65 concentrations negatively correlated with the number of aneurysms. In UIA patients the median GRO-α concentration was two-fold and CXCR2 almost four-fold higher in CSF compared to the serum value. CSF GRO-α concentration positively correlated with the size of aneurysms.Significantly decreased CSF NF-κB p65 and significantly increased CSF GRO-α and its CXCR2 receptor concentrations in UIA patients compared to the control group may altogether suggest that the canonical NF-κB signaling pathway is activated and its target pro-inflammatory genes are highly expressed in UIA patients. However, to unequivocally assess the involvement of the classical NF-κB pathway with the participation of the NF-κB p65 subunit and the GRO-α/CXCR2 axis in the formation of IA, further in vivo model studies are needed.

A two-stage hybrid gene selection algorithm combined with machine learning models to predict the rupture status in intracranial aneurysms

An IA is an abnormal swelling of cerebral vessels, and a subset of these IAs can rupture causing aneurysmal subarachnoid hemorrhage (aSAH), often resulting in death or severe disability. Few studies have used an appropriate method of feature selection combined with machine learning by analyzing transcriptomic sequencing data to identify new molecular biomarkers. Following gene ontology (GO) and enrichment analysis, we found that the distinct status of IAs could lead to differential innate immune responses using all 913 differentially expressed genes, and considering that there are numerous irrelevant and redundant genes, we propose a mixed filter- and wrapper-based feature selection. First, we used the Fast Correlation-Based Filter (FCBF) algorithm to filter a large number of irrelevant and redundant genes in the raw dataset, and then used the wrapper feature selection method based on the he Multi-layer Perceptron (MLP) neural network and the Particle Swarm Optimization (PSO), accuracy (ACC) and mean square error (MSE) were then used as the evaluation criteria. Finally, we constructed a novel 10-gene signature (YIPF1, RAB32, WDR62, ANPEP, LRRCC1, AADAC, GZMK, WBP2NL, PBX1, and TOR1B) by the proposed two-stage hybrid algorithm FCBF-MLP-PSO and used different machine learning models to predict the rupture status in IAs. The highest ACC value increased from 0.817 to 0.919 (12.5% increase), the highest area under ROC curve (AUC) value increased from 0.87 to 0.94 (8.0% increase), and all evaluation metrics improved by approximately 10% after being processed by our proposed gene selection algorithm. Therefore, these 10 informative genes used to predict rupture status of IAs can be used as complements to imaging examinations in the clinic, meanwhile, this selected gene signature also provides new targets and approaches for the treatment of ruptured IAs.

Identification of the key immune-related genes in aneurysmal subarachnoid hemorrhage

Subarachnoid hemorrhage (SAH) is a major cause of death and morbidity worldwide, often due to rupture of intracranial aneurysms (IAs). Immune infiltration and inflammatory activation play key roles in the process of aneurysmal SAH (aSAH). This study aimed to elaborate the immune infiltration and identify related biomarkers both in blood and tissue samples from patients with aSAH. Expression data of aSAH and healthy control samples were obtained from gene expression omnibus (GEO) database. Overall, a blood sample dataset GSE36791 and a tissue sample dataset GSE122897 were included. Differentially expressed genes (DEGs) between aSAH and healthy samples were explored. We applied GO biological and Gene Set Enrichment Analyses (GSEA) processes to access the functional enrichment. Then feature elimination algorithms based on random forest were used to screen and verify the biomarkers of aSAH. We performed three computational algorithms including Cell type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT), Microenvironment Cell Populations-counter (MCPcounter), and xcell to evaluate the immune cell infiltration landscape to identify the unique infiltration characteristics associated with rupturing. We found 2,220 DEGs (856 upregulated and 1,364 downregulated) in the original dataset. Functional analysis revealed most of these genes are enriched in immunological process, especially related with neutrophil response. Similar signaling pathway enrichment patterns were observed in tissue sample dataset and ClueGo. Analysis of immune microenvironment infiltration suggested neutrophils were abnormally upregulated in aSAH compared with those in the control group. Key gene SRPK1 was then filtered based on feature elimination algorithms, and transcription factor (TF) ZNF281 is assumed to participate in immunomodulation by regulating expression of SRPK1. Several immunomodulators such as CXCR1 and CXCR2 also appear to be involved in the progression of aSAH. In the present study, we performed a comprehensive stratification and quantification of the immune infiltration status of aSAH. By exploring the potential mechanism for aSAH based on several computational algorithms, key genes including SRPK1 and ZNF281 were filtered. This study may be of benefit to patients who are at high risk of suffering aSAH which allows for early diagnosis and potential therapy.

Intracranial aneurysm wall (in)stability-current state of knowledge and clinical perspectives

Intracranial aneurysm (IA), a local outpouching of cerebral arteries, is present in 3 to 5% of the population. Once formed, an IA can remain stable, grow, or rupture. Determining the evolution of IAs is almost impossible. Rupture of an IA leads to subarachnoid hemorrhage and affects mostly young people with heavy consequences in terms of death, disabilities, and socioeconomic burden. Even if the large majority of IAs will never rupture, it is critical to determine which IA might be at risk of rupture. IA (in)stability is dependent on the composition of its wall and on its ability to repair. The biology of the IA wall is complex and not completely understood. Nowadays, the risk of rupture of an IA is estimated in clinics by using scores based on the characteristics of the IA itself and on the anamnesis of the patient. Classification and prediction using these scores are not satisfying and decisions whether a patient should be observed or treated need to be better informed by more reliable biomarkers. In the present review, the effects of known risk factors for rupture, as well as the effects of biomechanical forces on the IA wall composition, will be summarized. Moreover, recent advances in high-resolution vessel wall magnetic resonance imaging, which are promising tools to discriminate between stable and unstable IAs, will be described. Common data elements recently defined to improve IA disease knowledge and disease management will be presented. Finally, recent findings in genetics will be introduced and future directions in the field of IA will be exposed.

Isolation of RNA from Acute Ischemic Stroke Clots Retrieved by Mechanical Thrombectomy

Mechanical thrombectomy (MT) for large vessel acute ischemic stroke (AIS) has enabled biologic analyses of resected clots. While clot histology has been well-studied, little is known about gene expression within the tissue, which could shed light on stroke pathophysiology. In this methodological study, we develop a pipeline for obtaining useful RNA from AIS clots. A total of 73 clot samples retrieved by MT were collected and stored in RNALater and in 10% phosphate-buffered formalin. RNA was extracted from all samples using a modified Chemagen magnetic bead extraction protocol on the PerkinElmer Chemagic 360. RNA was interrogated by UV–Vis absorption and electrophoretic quality control analysis. All samples with sufficient volume underwent traditional qPCR analysis and samples with sufficient RNA quality were subjected to next-generation RNA sequencing on the Illumina NovaSeq platform. Whole blood RNA samples from three patients were used as controls, and H&E-stained histological sections of the clots were used to assess clot cellular makeup. Isolated mRNA was eluted into a volume of 140 µL and had a concentration ranging from 0.01 ng/µL to 46 ng/µL. Most mRNA samples were partially degraded, with RNA integrity numbers ranging from 0 to 9.5. The majority of samples (71/73) underwent qPCR analysis, which showed linear relationships between the expression of three housekeeping genes (GAPDH, GPI, and HPRT1) across all samples. Of these, 48 samples were used for RNA sequencing, which had moderate quality based on MultiQC evaluation (on average, ~35 M reads were sequenced). Analysis of clot histology showed that more acellular samples yielded RNA of lower quantity and quality. We obtained useful mRNA from AIS clot samples stored in RNALater. qPCR analysis could be performed in almost all cases, while sequencing data could only be performed in approximately two-thirds of the samples. Acellular clots tended to have lower RNA quantity and quality.

A Two-Gene-Based Diagnostic Signature for Ruptured Intracranial Aneurysms

Background: Ruptured intracranial aneurysm (IA) is a disease with high mortality. Despite the great progress in treating ruptured IA, methods for risk assessment of ruptured IA remain limited. Methods: In this study, we aim to develop a robust diagnostic model for ruptured IA. Gene expression profiles in blood samples of 18 healthy persons and 43 ruptured IA patients were obtained from the Gene Expression Omnibus (GEO). Differential expression analysis was performed using limma Bioconductor package followed by functional enrichment analysis via clusterProfiler Bioconductor package. Immune cell compositions in ruptured IA and healthy samples were assessed through the CIBERSORT tool. Protein-protein interaction (PPI) was predicted based on the STRING database. Logistic regression model was used for the construction of predictive model for distinguishing ruptured IA and healthy samples. Real-time quantitative polymerase chain reaction (RT-qPCR) was performed to validate the gene expression between the ruptured IA and healthy samples. Results: A total of 58 differentially expressed genes (DEGs) were obtained for ruptured IA patients compared with healthy controls. Functional enrichment analysis showed that the DEGs were enriched in biological processes related to neutrophil activation, neutrophil degranulation, and cytokine-cytokine receptor interaction. Notably, immune analysis results proved that the rupture of IA might be related to immune cell distribution. We further identified 24 key genes as hub genes using the PPI networks. The logistic regression model trained based on the 24 key genes ultimately retained two genes, i.e., IL2RB and CCR7, which had great potential for risk assessment for rupture of IA. The RT-qPCR further validated that compared with the healthy samples, the expression levels of IL2RB and CCR7 were decreased in ruptured IA samples. Conclusions: This study might be helpful for cohorts who have a high risk of ruptured IA for early diagnosis and prevention of the disease.

Aberrant Whole Blood Gene Expression in the Lumen of Human Intracranial Aneurysms

The rupture of an intracranial aneurysm (IA) causes devastating hemorrhagic strokes. Yet, most IAs remain asymptomatic and undetected until they rupture. In the search for circulating biomarkers of unruptured IAs, we previously performed transcriptome profiling on whole blood and identified an IA-associated panel of 18 genes. In this study, we seek to determine if these genes are also differentially expressed within the IA lumen, which could provide a mechanistic link between the disease and the observed circulating gene expression patterns. To this end, we collected blood from the lumen of 37 IAs and their proximal parent vessels in 31 patients. The expression levels of 18 genes in the lumen and proximal vessel were then measured by quantitative polymerase chain reaction. This analysis revealed that the expression of 6/18 genes (CBWD6, MT2A, MZT2B, PIM3, SLC37A3, and TNFRSF4) was significantly higher in intraluminal blood, while the expression of 3/18 genes (ST6GALNAC1, TCN2, and UFSP1) was significantly lower. There was a significant, positive correlation between intraluminal and proximal expression of CXCL10, MT2A, and MZT2B, suggesting local increases of these genes is reflected in the periphery. Expression of ST6GALNAC1 and TIFAB was significantly positively correlated with IA size, while expression of CCDC85B was significantly positively correlated with IA enhancement on post-contrast MRI, a metric of IA instability and risk. In conclusion, intraluminal expression differences in half of the IA-associated genes observed in this study provide evidence for IA tissue-mediated transcriptional changes in whole blood. Additionally, some genes may be informative in assessing IA risk, as their intraluminal expression was correlated to IA size and aneurysmal wall enhancement.

Identification of Circulating Gene Expression Signatures of Intracranial Aneurysm in Peripheral Blood Mononuclear Cells

Peripheral blood mononuclear cells (PBMCs) play an important role in the inflammation that accompanies intracranial aneurysm (IA) pathophysiology. We hypothesized that PBMCs have different transcriptional profiles in patients harboring IAs as compared to IA-free controls, which could be the basis for potential blood-based biomarkers for the disease. To test this, we isolated PBMC RNA from whole blood of 52 subjects (24 with IA, 28 without) and performed next-generation RNA sequencing to obtain their transcriptomes. In a randomly assigned discovery cohort of n = 39 patients, we performed differential expression analysis to define an IA-associated signature of 54 genes (q < 0.05 and an absolute fold-change ≥ 1.3). In the withheld validation dataset, these genes could delineate patients with IAs from controls, as the majority of them still had the same direction of expression difference. Bioinformatics analyses by gene ontology enrichment analysis and Ingenuity Pathway Analysis (IPA) demonstrated enrichment of structural regulation processes, intracellular signaling function, regulation of ion transport, and cell adhesion. IPA analysis showed that these processes were likely coordinated through NF-kB, cytokine signaling, growth factors, and TNF activity. Correlation analysis with aneurysm size and risk assessment metrics showed that 4/54 genes were associated with rupture risk. These findings highlight the potential to develop predictive biomarkers from PBMCs to identify patients harboring IAs.