Circulating neutrophil transcriptome may reveal intracranial aneurysm signature

Screening key genes for intracranial aneurysm rupture using LASSO regression and the SVM-RFE algorithm

Background

Although an intracranial aneurysm (IA) is widespread and fatal, few drugs can be used to prevent its rupture. This study explored the molecular mechanism and potential targets of IA rupture through bioinformatics methods.

Methods

The gene expression matrices of GSE13353, GSE122897, and GSE15629 were downloaded. Differentially expressed genes (DEGs) were screened using the limma package. Functional enrichment analysis was performed, and a PPI network was constructed. Furthermore, candidate key genes were identified using the least absolute shrinkage and selection operator (LASSO) regression model, support vector machine-recursive feature elimination (SVM-RFE) analysis, and PPI network analysis. ROC analysis was conducted to further verify the diagnostic value of the key genes.

Results

A total of 334 DEGs were screened, including 175 upregulated genes and 159 downregulated genes. Further functional analysis suggested that the DEGs were enriched in inflammation and immune response pathways. Fourteen hub genes were identified using the two algorithms. The PPI networks of the hub genes were analyzed using the Cytoscape plugin CytoNCA to obtain two key genes (IL10 and Integrin α5 (ITGA5)). The ROC curve analysis showed that the AUC values of IL10 and ITGA5 were 0.801, and 0.786, respectively. In addition, the two key genes were significantly positively correlated with macrophages and Treg (T) cells. The immune score and ESTIMATE score of the ruptured IA group were significantly higher than those of the unruptured IA group.

Conclusion

The increase in IL-10 and ITGA5 may weaken the vascular wall by promoting inflammation in blood vessels and immune cells, which could have a harmful effect on the rupture of IAs.

Single-cell atlas of peripheral blood by CyTOF revealed peripheral blood immune cells metabolic alterations and neutrophil changes in intracranial aneurysm rupture

Previous studies have found that the peripheral immune environment is closely related to the occurrence and development of intracranial aneurysms. However, it remains unclear how the metabolism of peripheral blood mononuclear cells (PBMCs) and the composition of polymorphonuclear leukocytes (PMNs) changes in the process of intracranial aneurysm rupture. This study utilized cytometry by time of flight technology to conduct single-cell profiling analysis of PBMCs and PMNs from 72 patients with IAs. By comparing the expression differences of key metabolic enzymes in PBMCs between patients with ruptured intracranial aneurysms (RIAs) and unruptured intracranial aneurysms, we found that most PBMCs subsets from RIA group showed upregulation of rate-limiting enzymes related to the glycolytic pathway. By comparing the composition of PMNs, it was found that the proinflammatory CD101+HLA DR+ subsets were increased in the RIA group, accompanied by a decrease in the anti-inflammatory polymorphonuclear myeloid-derived suppressor cells. In conclusion, this study showed the changes in the peripheral immune profile of RIAs, which is helpful for our understanding of the mechanisms underlying peripheral changes and provides a direction for future related research.

MicroCT and Histological Analysis of Clot Composition in Acute Ischemic Stroke : A Comparative Study of MT-Retrieved Clots and Clot Analogs

PURPOSE

Assessing clot composition on prethrombectomy computed tomography (CT) imaging may help in stroke treatment planning. In this study we seek to use microCT imaging of fabricated blood clots to understand the relationship between CT radiographic signals and the biological makeup.

METHODS

Clots (n = 10) retrieved by mechanical thrombectomy (MT) were collected, and 6 clot analogs of varying RBC composition were made. We performed paired microCT and histological image analysis of all 16 clots using a ScanCo microCT 100 (4.9 µm resolution) and standard H&E staining (imaged at 40×). From these data types, first order statistic (FOS) radiomics were computed from microCT, and percent composition of RBCs (%RBC) was computed from histology. Polynomial and linear regression (LR) were used to build statistical models based on retrieved thrombus microCT and %RBC that were evaluated for their ability to predict the %RBC of clot analogs from mean HU. Correlation analyses of microCT FOS with composition were completed for both retrieved clots and analogs.

RESULTS

The LR model fits relating MT-retrieved clot %RBC with mean (R2 = 0.625, p = 0.006) and standard deviation (R2 = 0.564, p < 0.05) in HUs on microCT were significant. Similarly, LR models relating analog histological %RBC to analog protocol %RBC (R2 = 0.915, p = 0.003) and mean HUs on microCT (R2 = 0.872, p = 0.007) were also significant. When the LR model built using MT-retrieved clots was used to predict analog %RBC from mean HUs, significant correlation was observed between predictions and actual histological %RBC (R2 = 0.852, p = 0.009). For retrieved clots, significant correlations were observed for energy and total energy with %RBC and %FP (|R| > 0.7, q < 0.01). Analogs further demonstrated significant correlation between FOS energy, total energy, variance and %WBC (|R| > 0.9, q < 0.01).

CONCLUSION

MicroCT can be used to build models that predict AIS clot composition from routine CT parameters and help us to better understand radiomic signatures associated with clot composition and first pass outcomes. In future work, such observations can be used to better infer clot composition and inform thrombectomy prognostics from pretreatment CTs.

Metabolic alterations of peripheral blood immune cells and heterogeneity of neutrophil in intracranial aneurysms patients

BACKGROUND

Intracranial aneurysms (IAs) represent a severe cerebrovascular disease that can potentially lead to subarachnoid haemorrhage. Previous studies have demonstrated the involvement of peripheral immune cells in the formation and progression of IAs. Nevertheless, the impact of metabolic alterations in peripheral immune cells and changes in neutrophil heterogeneity on the occurrence and progression of IAs remains uncertain.

METHODS

Single-cell Cytometry by Time-of-Flight (CyTOF) technology was employed to profile the single-cell atlas of peripheral blood mononuclear cells (PBMCs) and polymorphonuclear cells (PMNs) in 72 patients with IAs. In a matched cohort, metabolic shifts in PBMC subsets of IA patients were investigated by contrasting the expression levels of key metabolic enzymes with their respective counterparts in the healthy control group. Simultaneously, compositional differences in peripheral blood PMNs subsets between the two groups were analysed to explore the impact of altered heterogeneity in neutrophils on the initiation and progression of IAs. Furthermore, integrating immune features based on CyTOF analysis and clinical characteristics, we constructed an aneurysm occurrence model and an aneurysm growth model using the random forest method in conjunction with LASSO regression.

RESULTS

Different subsets exhibited distinct metabolic characteristics. Overall, PBMCs from patients elevated CD98 expression and increased proliferation. Conversely, CD36 was up-regulated in T cells, B cells and monocytes from the controls but down-regulated in NK and NKT cells. The comparison also revealed differences in the metabolism and function of specific subsets between the two groups. In terms of PMNs, the neutrophil landscape within patients group revealed a pronounced shift towards heightened complexity. Various neutrophil subsets from the IA group generally exhibited lower expression levels of anti-inflammatory functional molecules (IL-4 and IL-10). By integrating clinical and immune features, the constructed aneurysm occurrence model could precisely identify patients with IAs with high prediction accuracy (AUC = 0.987). Furthermore, the aneurysm growth model also exhibited superiority over ELAPSS scores in predicting aneurysm growth (lower prediction errors and out-of-bag errors).

CONCLUSION

These findings enhanced our understanding of peripheral immune cell participation in aneurysm formation and growth from the perspectives of immune metabolism and neutrophil heterogeneity. Moreover, the predictive model based on CyTOF features holds the potential to aid in diagnosing and monitoring the progression of human IAs.

Gene expression profiles of ischemic stroke clots retrieved by mechanical thrombectomy are associated with disease etiology

BACKGROUND

Determining stroke etiology is crucial for secondary prevention, but intensive workups fail to classify ~30% of strokes that are cryptogenic.

OBJECTIVE

To examine the hypothesis that the transcriptomic profiles of clots retrieved during mechanical thrombectomy are unique to strokes of different subtypes.

METHODS

We isolated RNA from the clots of 73 patients undergoing mechanical thrombectomy. Samples of sufficient quality were subjected to 100-cycle, paired-end RNAseq, and transcriptomes with less than 10 million unique reads were excluded from analysis. Significant differentially expressed genes (DEGs) between subtypes (defined by the Trial of Org 10 172 in Acute Stroke Treatment) were identified by expression analysis in edgeR. Gene ontology enrichment analysis was used to study the biologic differences between stroke etiologies.

RESULTS

In all, 38 clot transcriptomes were analyzed; 6 from large artery atherosclerosis (LAA), 21 from cardioembolism (CE), 5 from strokes of other determined origin, and 6 from cryptogenic strokes. Among all comparisons, there were 816 unique DEGs, 174 of which were shared by at least two comparisons, and 20 of which were shared by all three. Gene ontology analysis showed that CE clots reflected high levels of inflammation, LAA clots had greater oxidoreduction and T-cell processes, and clots of other determined origin were enriched for aberrant platelet and hemoglobin-related processes. Principal component analysis indicated separation between these subtypes and showed cryptogenic samples clustered among several different groups.

CONCLUSIONS

Expression profiles of stroke clots were identified between stroke etiologies and reflected different biologic responses. Cryptogenic thrombi may be related to multiple etiologies.

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.

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.

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.

Machine Learning and Intracranial Aneurysms: From Detection to Outcome Prediction

Machine learning (ML) is a rapidly rising research tool in biomedical sciences whose applications include segmentation, classification, disease detection, and outcome prediction. With respect to traditional statistical methods, ML algorithms have the potential to learn and improve their predictive performance when fed with large data sets without the need of being specifically programmed. In recent years, this technology has been increasingly applied for tackling clinical issues in intracranial aneurysm (IA) research. Several studies attempted to provide reliable models for enhanced aneurysm detection. Convolutional neural networks trained with variable degrees of human interaction on data from diverse imaging modalities showed high sensitivity in aneurysm detection tasks, also outperforming expert image analysis. Algorithms were also shown to differentiate ruptured from unruptured IAs, with however limited clinical relevance. For prediction of rupture and stability assessment, ML was preliminarily shown to achieve better performance compared to conventional statistical methods and existing risk scores. ML-based complication and functional outcome prediction in the event of SAH have been more extensively reported, in contrast with periprocedural outcome investigation in unruptured IA patients. ML has the potential to be a game changer in IA patient management. Currently clinical translation of experimental results is limited.

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.

Characterization of Long Non-coding RNA Signatures of Intracranial Aneurysm in Circulating Whole Blood

BACKGROUND AND OBJECTIVE

Long non-coding RNAs (lncRNAs) may serve as biomarkers for complex disease states, such as intracranial aneurysms. In this study, we investigated lncRNA expression differences in the whole blood of patients with unruptured aneurysms.

METHODS

Whole blood RNA from 67 subjects (34 with aneurysm, 33 without) was used for next-generation RNA sequencing. Differential expression analysis was used to define a signature of intracranial aneurysm-associated lncRNAs. To estimate the signature's ability to classify aneurysms and to identify the most predictive lncRNAs, we implemented a nested cross-validation pipeline to train classifiers using linear discriminant analysis. Ingenuity pathway analysis was used to study potential biological roles of differentially expressed lncRNAs, and lncRNA ontology was used to investigate ontologies enriched in signature lncRNAs. Co-expression correlation analysis was performed to investigate associated differential protein-coding messenger RNA expression.

RESULTS

Of 4639 detected lncRNAs, 263 were significantly different (p < 0.05) between the two groups, and 84 of those had an absolute fold-change ≥ 1.5. An eight-lncRNA signature (q < 0.35, fold-change ≥ 1.5) was able to separate patients with and without aneurysms on principal component analysis, and had an estimated accuracy of 70.9% in nested cross-validation. Bioinformatics analyses showed that networks of differentially expressed lncRNAs (p < 0.05) were enriched for cell death and survival, connective tissue disorders, carbohydrate metabolism, and cardiovascular disease. Signature lncRNAs shared ontologies that reflected regulation of gene expression, signaling, ubiquitin processing, and p53 signaling. Co-expression analysis showed correlations with messenger RNAs related to inflammatory responses.

CONCLUSIONS

Differential expression in whole blood lncRNAs is detectable in patients harboring aneurysms, and reflects expression/signaling regulation, and ubiquitin and p53 pathways. Following validation in larger cohorts, these lncRNAs may be potential diagnostic targets for aneurysm detection by blood testing.

Inflammatory Biomarkers and Intracranial Hemorrhage after Endovascular Thrombectomy

BACKGROUND

Intracranial hemorrhage after endovascular thrombectomy is associated with poorer prognosis compared with those who do not develop the complication. Our study aims to determine predictors of post-EVT hemorrhage - more specifically, inflammatory biomarkers present in baseline serology.

METHODS

We performed a retrospective review of consecutive patients treated with EVT for acute large vessel ischemic stroke. The primary outcome of the study is the presence of ICH on the post-EVT scan. We used four definitions: the SITS-MOST criteria, the NINDS criteria, asymptomatic hemorrhage, and overall hemorrhage. We identified nonredundant predictors of outcome using backward elimination based on Akaike Information Criteria. We then assessed prediction accuracy using area under the receiver operating curve. Then we implemented variable importance ranking from logistic regression models using the drop in Naegelkerke R2 with the exclusion of each predictor.

RESULTS

Our study demonstrates a 6.3% SITS (16/252) and 10.0% NINDS (25/252) sICH rate, as well as a 19.4% asymptomatic (49/252) and 29.4% (74/252) overall hemorrhage rate. Serologic markers that demonstrated association with post-EVT hemorrhage were: low lymphocyte count (SITS), high neutrophil count (NINDS, overall hemorrhage), low platelet to lymphocyte ratio (NINDS), and low total WBC (NINDS, asymptomatic hemorrhage).

CONCLUSION

Higher neutrophil counts, low WBC counts, low lymphocyte counts, and low platelet to lymphoycyte ratio were baseline serology biomarkers that were associated with post-EVT hemorrhage. Our findings, particularly the association of diabetes mellitus and high neutrophil, support experimental data on the role of thromboinflammation in hemorrhagic transformation of large vessel occlusions.

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.

Transcriptional firing represses bactericidal activity in cystic fibrosis airway neutrophils

Neutrophils are often considered terminally differentiated and poised for bacterial killing. In chronic diseases such as cystic fibrosis (CF), an unexplained paradox pits massive neutrophil presence against prolonged bacterial infections. Here, we show that neutrophils recruited to CF airways in vivo and in an in vitro transmigration model display rapid and broad transcriptional firing, leading to an upregulation of anabolic genes and a downregulation of antimicrobial genes. Newly transcribed RNAs are mirrored by the appearance of corresponding proteins, confirming active translation in these cells. Treatment by the RNA polymerase II and III inhibitor α-amanitin restores the expression of key antimicrobial genes and increases the bactericidal capacity of CF airway neutrophils in vitro and in short-term sputum cultures ex vivo. Broadly, our findings show that neutrophil plasticity is regulated at the site of inflammation via RNA and protein synthesis, leading to adaptations that affect their canonical functions (i.e., bacterial clearance).

Endogenous animal models of intracranial aneurysm development: a review

The pathogenesis and natural history of intracranial aneurysm (IA) remains poorly understood. To this end, animal models with induced cerebral vessel lesions mimicking human aneurysms have provided the ability to greatly expand our understanding. In this review, we comprehensively searched the published literature to identify studies that endogenously induced IA formation in animals. Studies that constructed aneurysms (i.e., by surgically creating a sac) were excluded. From the eligible studies, we reported information including the animal species, method for aneurysm induction, aneurysm definitions, evaluation methods, aneurysm characteristics, formation rate, rupture rate, and time course. Between 1960 and 2019, 174 articles reported endogenous animal models of IA. The majority used flow modification, hypertension, and vessel wall weakening (i.e., elastase treatment) to induce IAs, primarily in rats and mice. Most studies utilized subjective or qualitative descriptions to define experimental aneurysms and histology to study them. In general, experimental IAs resembled the pathobiology of the human disease in terms of internal elastic lamina loss, medial layer degradation, and inflammatory cell infiltration. After the early 2000s, many endogenous animal models of IA began to incorporate state-of-the-art technology, such as gene expression profiling and 9.4-T magnetic resonance imaging (MRI) in vivo imaging, to quantitatively analyze the biological mechanisms of IA. Future studies aimed at longitudinally assessing IA pathobiology in models that incorporate aneurysm growth will likely have the largest impact on our understanding of the disease. We believe this will be aided by high-resolution, small animal, survival imaging, in situ live-cell imaging, and next-generation omics technology.

Whole blood transcriptome biomarkers of unruptured intracranial aneurysm

BACKGROUND

The rupture of an intracranial aneurysm (IA) causes devastating subarachnoid hemorrhages, yet most IAs remain undiscovered until they rupture. Recently, we found an IA RNA expression signature of circulating neutrophils, and used transcriptome data to build predictive models for unruptured IAs. In this study, we evaluate the feasibility of using whole blood transcriptomes to predict the presence of unruptured IAs.

METHODS

We subjected RNA from peripheral whole blood of 67 patients (34 with unruptured IA, 33 without IA) to next-generation RNA sequencing. Model genes were identified using the least absolute shrinkage and selection operator (LASSO) in a random training cohort (n = 47). These genes were used to train a Gaussian Support Vector Machine (gSVM) model to distinguish patients with IA. The model was applied to an independent testing cohort (n = 20) to evaluate performance by receiver operating characteristic (ROC) curve. Gene ontology and pathway analyses investigated the underlying biology of the model genes.

RESULTS

We identified 18 genes that could distinguish IA patients in a training cohort with 85% accuracy. This SVM model also had 85% accuracy in the testing cohort, with an area under the ROC curve of 0.91. Bioinformatics reflected activation and recruitment of leukocytes, activation of macrophages, and inflammatory response, suggesting that the biomarker captures important processes in IA pathogenesis.

CONCLUSIONS

Circulating whole blood transcriptomes can detect the presence of unruptured IAs. Pending additional testing in larger cohorts, this could serve as a foundation to develop a simple blood-based test to facilitate screening and early detection of IAs.

Identification and validation of key genes mediating intracranial aneurysm rupture by weighted correlation network analysis

Background

Rupture of intracranial aneurysm (IA) is the leading cause of subarachnoid hemorrhage. However, there are few pharmacological therapies available for the prevention of IA rupture. Therefore, exploring the molecular mechanisms which underlie IA rupture and identifying the potential molecular targets for preventing the rupture of IA is of vital importance.

Methods

We used the Gene Expression Omnibus (GEO) datasets GSE13353, GSE15629, and GSE54083 in our study. The 3 datasets were merged and normalized. Differentially expressed gene (DEG) screening and weighted correlation network analysis (WGCNA) were conducted. The co-expression patterns between ruptured IA samples and unruptured IA samples were compared. Then, the DEGs were mapped into the whole co-expression network of ruptured IA samples, and a DEG co-expression network was generated. Molecular Complex Detection (MCODE) (http://baderlab.org/Software/MCODE) was used to identify key genes based on the DEG co-expression network. Finally, key genes were validated using another GEO dataset (GSE122897), and their potential diagnostic values were shown using receiver operating characteristic (ROC) analysis.

Results

In our study, 49 DEGs were screened while 8 and 6 gene modules were detected based on ruptured IA samples and unruptured IA samples, respectively. Pathways associated with inflammation and immune response were clustered in the salmon module of ruptured IA samples. The DEG co-expression network with 35 nodes and 168 edges was generated, and 14 key genes were identified based on this DEG co-expression network. The gene with the highest degree in the key gene cluster was CXCR4. All key genes were validated using GSE122897, and they all showed the potential diagnostic value in predicting IA rupture.

Conclusions

Using a weighted gene co-expression network approach, we identified 8 and 6 modules for ruptured IA and unruptured IA, respectively. After that, we identified the hub genes for each module and key genes based on the DEG co-expression network. All these key genes were validated by another GEO dataset and might serve as potential targets for pharmacological therapies and diagnostic markers in predicting IA rupture. Further studies are needed to elucidate the detailed molecular mechanisms and biological functions of these key genes which underlie the rupture of IA.

Classification models using circulating neutrophil transcripts can detect unruptured intracranial aneurysm

Background

Intracranial aneurysms (IAs) are dangerous because of their potential to rupture. We previously found significant RNA expression differences in circulating neutrophils between patients with and without unruptured IAs and trained machine learning models to predict presence of IA using 40 neutrophil transcriptomes. Here, we aim to develop a predictive model for unruptured IA using neutrophil transcriptomes from a larger population and more robust machine learning methods.

Methods

Neutrophil RNA extracted from the blood of 134 patients (55 with IA, 79 IA-free controls) was subjected to next-generation RNA sequencing. In a randomly-selected training cohort (n = 94), the Least Absolute Shrinkage and Selection Operator (LASSO) selected transcripts, from which we constructed prediction models via 4 well-established supervised machine-learning algorithms (K-Nearest Neighbors, Random Forest, and Support Vector Machines with Gaussian and cubic kernels). We tested the models in the remaining samples (n = 40) and assessed model performance by receiver-operating-characteristic (ROC) curves. Real-time quantitative polymerase chain reaction (RT-qPCR) of 9 IA-associated genes was used to verify gene expression in a subset of 49 neutrophil RNA samples. We also examined the potential influence of demographics and comorbidities on model prediction.

Results

Feature selection using LASSO in the training cohort identified 37 IA-associated transcripts. Models trained using these transcripts had a maximum accuracy of 90% in the testing cohort. The testing performance across all methods had an average area under ROC curve (AUC) = 0.97, an improvement over our previous models. The Random Forest model performed best across both training and testing cohorts. RT-qPCR confirmed expression differences in 7 of 9 genes tested. Gene ontology and IPA network analyses performed on the 37 model genes reflected dysregulated inflammation, cell signaling, and apoptosis processes. In our data, demographics and comorbidities did not affect model performance.

Conclusions

We improved upon our previous IA prediction models based on circulating neutrophil transcriptomes by increasing sample size and by implementing LASSO and more robust machine learning methods. Future studies are needed to validate these models in larger cohorts and further investigate effect of covariates.

Genetic polymorphisms and transcription profiles associated with intracranial aneurysm: a key role for NOTCH3

Intracranial aneurysm (IA) incidence is about 1~2%. However, the specific mechanisms of IA onset and development need further study. Our objective was to discover novel IA-related genes to determine possible etiologies further. We performed next-generation sequencing on nineteen Chinese patients with familial IA and one patient with sporadic IA. We obtained mRNA expression data of 129 samples from Gene Expression Omnibus (GEO) and made statistical computing to discover differentially expressed genes (DEGs). The screened IA-related gene NOTCH3 was determined by bioinformatic data mining. We verified the IA-related indicators of NOTCH3. Association was found between IA and the NOTCH3 SNPs rs779314594, rs200504060 and rs2285981. Levels of NOTCH3 mRNA were lower in IA tissue than in control tissue, but higher in peripheral blood neutrophils from IA patients than in neutrophils from controls. Levels of NOTCH3 protein were lower in IA tissue than in cerebral artery tissue. NOTCH3 also decreased the expression of angiogenesis factors in human umbilical vein endothelial cells. Variation in NOTCH3 and alteration of its expression in cerebral artery or neutrophils may contribute to IA. Our findings also describe a bioinformatic-experimental approach that may prove useful for probing the pathophysiology of other complex diseases.

Identification of potential key pathways, genes and circulating markers in the development of intracranial aneurysm based on weighted gene co-expression network analysis

Background: Intracranial aneurysm (IA) is a disease resulted from weak brain control, characterized by local expansion or dilation of brain artery. This study aimed to construct a gene co-expression network by Weighted Gene Correlation Network Analysis (WGCNA) to explore the potential key pathways and genes for the development of IA.Method: Six IA-related gene expression data sets were downloaded from the Gene Expression Omnibus (GEO) database for identifying differentially expressed genes (DEGs). WGCNA was used to identify modules associated with IA. Functional enrichment analysis was used to explore the potential biological functions. ROC analysis was used to find markers for predicting IA.Results: Purple, greenyellow and yellow modules were significantly associated with unruptured intracranial aneurysms, while blue and turquoise modules were significantly associated with ruptured intracranial aneurysms. Functional modules significantly related to IA were enriched in Ribosome, Glutathione metabolism, cAMP signalling pathway, Lysosome, Glycosaminoglycan degradation and other pathways. CD163, FCEREG, FPR1, ITGAM, NLRC4, PDG, and TYROBP were up-regulated ruptured intracranial aneurysms and serum, these genes were potential circulating markers for predicting IA rupture.Conclusions: Potential IA-related key pathways, genes and circulating markers were identified for predicting IA rupture by WGCNA analysis.

Human primary neutrophil mRNA does not contaminate human resolving macrophage mRNA after efferocytosis

The ingestion of apoptotic corpses by macrophages, a process called efferocytosis, is a crucial step in inflammation resolution, since it alters macrophage phenotype toward a pro-resolving profile to foil inflammation and to favor tissue repair. Up to now, the resolving macrophages remain poorly characterized, especially in humans. Global investigations, like RNA sequencing, would be very helpful to unravel some features of these elusive cells. Nonetheless, these inquiries may be challenging in a single-species model, since the fate of ingested mRNA remains unknown and may hinder any subsequent mRNA investigations in the phagocyte. A full human model consisting of primary human neutrophil and primary human monocyte-derived macrophage co-culture was set up several decades ago to mimic in vitro the efferocytosis process. However, to our knowledge, this model has not been characterized as a suitable model to perform global mRNA investigations. Indeed, the extent of ingested neutrophil mRNA contamination has not been assessed in resolving macrophages. This work answers to this crucial question. Indeed, based on the protocols presented in this article, we demonstrate that neutrophil mRNA is severely degraded and is not able to cross-contaminate resolving macrophage mRNA, contrary to apoptotic human peripheral blood derived mononuclear cell (PBMC) or apoptotic leukemic Jurkat cell mRNA. Moreover, this allogenic co-culture system does not favor neither neutrophil activation nor macrophage pro-inflammatory cytokine release. Collectively, we highlight that this model of primary human neutrophil and primary human monocyte-derived macrophage co-culture is the best model for mRNA investigations in human resolving macrophages to help improving our knowledge on these crucial cells.

The feasibility of developing biomarkers from peripheral blood mononuclear cell RNAseq data in children with juvenile idiopathic arthritis using machine learning approaches

Background

The response to treatment for juvenile idiopathic arthritis (JIA) can be staged using clinical features. However, objective laboratory biomarkers of remission are still lacking. In this study, we used machine learning to predict JIA activity from transcriptomes from peripheral blood mononuclear cells (PBMCs). We included samples from children with Native American ancestry to determine whether the model maintained validity in an ethnically heterogeneous population.

Methods

Our dataset consisted of 50 samples, 23 from children in remission and 27 from children with an active disease on therapy. Nine of these samples were from children with mixed European/Native American ancestry. We used 4 different machine learning methods to create predictive models in 2 populations: the whole dataset and then the samples from children with exclusively European ancestry.

Results

In both populations, models were able to predict JIA status well, with training accuracies > 74% and testing accuracies > 78%. Performance was better in the whole dataset model. We note a high degree of overlap between genes identified in both populations. Using ingenuity pathway analysis, genes from the whole dataset associated with cell-to-cell signaling and interactions, cell morphology, organismal injury and abnormalities, and protein synthesis.

Conclusions

This study demonstrates it is feasible to use machine learning in conjunction with RNA sequencing of PBMCs to predict JIA stage. Thus, developing objective biomarkers from easy to obtain clinical samples remains an achievable goal.

Epigenetic landscapes suggest that genetic risk for intracranial aneurysm operates on the endothelium

BACKGROUND

Genetics play an important role in intracranial aneurysm (IA) pathophysiology. Genome-wide association studies have identified several single nucleotide polymorphisms (SNPs) that are linked to IA but how they affect disease pathobiology remains poorly understood. We used Encyclopedia of DNA Elements (ENCODE) data to investigate the epigenetic landscapes surrounding genetic risk loci to determine if IA-associated SNPs affect functional elements that regulate gene expression and if those SNPs are most likely to impact a specific type of cells.

METHODS

We mapped 16 highly significant IA-associated SNPs to linkage disequilibrium (LD) blocks within the human genome. Within these regions, we examined the presence of H3K4me1 and H3K27ac histone marks and CCCTC-binding factor (CTCF) and transcription-factor binding sites using chromatin immunoprecipitation-sequencing (ChIP-Seq) data. This analysis was conducted in several cell types relevant to endothelial (human umbilical vein endothelial cells [HUVECs]) and inflammatory (monocytes, neutrophils, and peripheral blood mononuclear cells [PBMCs]) biology. Gene ontology analysis was performed on genes within extended IA-risk regions to understand which biological processes could be affected by IA-risk SNPs. We also evaluated recently published data that showed differential methylation and differential ribonucleic acid (RNA) expression in IA to investigate the correlation between differentially regulated elements and the IA-risk LD blocks.

RESULTS

The IA-associated LD blocks were statistically significantly enriched for H3K4me1 and/or H3K27ac marks (markers of enhancer function) in endothelial cells but not in immune cells. The IA-associated LD blocks also contained more binding sites for CTCF in endothelial cells than monocytes, although not statistically significant. Differentially methylated regions of DNA identified in IA tissue were also present in several IA-risk LD blocks, suggesting SNPs could affect this epigenetic machinery. Gene ontology analysis supports that genes affected by IA-risk SNPs are associated with extracellular matrix reorganization and endopeptidase activity.

CONCLUSION

These findings suggest that known genetic alterations linked to IA risk act on endothelial cell function. These alterations do not correlate with IA-associated gene expression signatures of circulating blood cells, which suggests that such signatures are a secondary response reflecting the presence of IA rather than indicating risk for IA.

Common functional alterations identified in blood transcriptome of autoimmune cholestatic liver and inflammatory bowel diseases

Primary biliary cholangitis (PBC), primary sclerosing cholangitis (PSC), and inflammatory bowel diseases (IBDs), including Crohn’s disease (CD) and ulcerative colitis (UC), are heterogeneous chronic autoimmune diseases that may share underlying pathogenic mechanisms. Herein, we compared simultaneously analyzed blood transcriptomes from patients with PBC, PSC, and IBD. Microarray-based measurements were conducted using RNA isolated from whole blood samples from 90, 45, 95 and 93 patients with PBC, PSC, CD, and UC, respectively, and 47 healthy controls. Expression levels of selected transcripts were analyzed by quantitative reverse-transcribed PCR using an independent cohort of 292, 71 and 727 patients with PBC, PSC, and IBD, respectively. Of 4026, 2650 and 4967 probe sets differentially expressed (adjusted p-value < 0.05) in samples from patients with PBC, PSC, and IBD, respectively, compared with healthy controls, 1946 were common to all three comparisons. Functional analyses indicated that most terms enriched for genes differentially expressed in PBC, PSC, and IBD patients compared with healthy controls were related to mitochondrial function, the vesicle endomembrane system, and GTPase-mediated processes. This study indicates that microarray-based profiling of blood gene expression supports research into the molecular mechanisms underlying disease, rather than being useful for selection of diagnostic biomarkers for use in clinical practice.

Biomarkers from circulating neutrophil transcriptomes have potential to detect unruptured intracranial aneurysms

BACKGROUND

Intracranial aneurysms (IAs) are dangerous because of their potential to rupture and cause deadly subarachnoid hemorrhages. Previously, we found significant RNA expression differences in circulating neutrophils between patients with unruptured IAs and aneurysm-free controls. Searching for circulating biomarkers for unruptured IAs, we tested the feasibility of developing classification algorithms that use neutrophil RNA expression levels from blood samples to predict the presence of an IA.

METHODS

Neutrophil RNA extracted from blood samples from 40 patients (20 with angiography-confirmed unruptured IA, 20 angiography-confirmed IA-free controls) was subjected to next-generation RNA sequencing to obtain neutrophil transcriptomes. In a randomly-selected training cohort of 30 of the 40 samples (15 with IA, 15 controls), we performed differential expression analysis. Significantly differentially expressed transcripts (false discovery rate < 0.05, fold change ≥ 1.5) were used to construct prediction models for IA using four well-known supervised machine-learning approaches (diagonal linear discriminant analysis, cosine nearest neighbors, nearest shrunken centroids, and support vector machines). These models were tested in a testing cohort of the remaining 10 neutrophil samples from the 40 patients (5 with IA, 5 controls), and model performance was assessed by receiver-operating-characteristic (ROC) curves. Real-time quantitative polymerase chain reaction (PCR) was used to corroborate expression differences of a subset of model transcripts in neutrophil samples from a new, separate validation cohort of 10 patients (5 with IA, 5 controls).

RESULTS

The training cohort yielded 26 highly significantly differentially expressed neutrophil transcripts. Models using these transcripts identified IA patients in the testing cohort with accuracy ranging from 0.60 to 0.90. The best performing model was the diagonal linear discriminant analysis classifier (area under the ROC curve = 0.80 and accuracy = 0.90). Six of seven differentially expressed genes we tested were confirmed by quantitative PCR using isolated neutrophils from the separate validation cohort.

CONCLUSIONS

Our findings demonstrate the potential of machine-learning methods to classify IA cases and create predictive models for unruptured IAs using circulating neutrophil transcriptome data. Future studies are needed to replicate these findings in larger cohorts.