- Open Access
Biomarkers from circulating neutrophil transcriptomes have potential to detect unruptured intracranial aneurysms
Journal of Translational Medicinevolume 16, Article number: 373 (2018)
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.
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).
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.
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.
Intracranial aneurysm (IA) rupture is the primary cause of nontraumatic subarachnoid hemorrhage and its devastating sequalae . The risk of rupture can be reduced by elective endovascular or surgical treatment [2, 3]. However, because IAs are almost invariably asymptomatic until rupture [4, 5], unruptured aneurysms are usually detected incidentally in individuals who are prescribed cerebral imaging for other reasons. A blood test capable of identifying individuals who harbor unruptured aneurysms would therefore be a significant advance in the field.
Circulating blood is a dynamic, information-rich tissue that continuously interacts with the aneurysm. In many different disease states, transcriptome profiling has been used to discover panels of differentially expressed genes in the circulating blood that could serve as useful diagnostic markers . This strategy has been successfully applied to complex vascular diseases, including IA [7,8,9]. In a recent case–controlled study, we performed transcriptome profiling on circulating neutrophils from patients with and without IA and discovered an IA-associated RNA expression signature (NCBI Gene Expression Omnibus Accession Number GSE106520) . This signature was characteristic of peripheral neutrophil activation and was able to separate patients with IA from controls in several statistical analyses. Although these results showed that differences in neutrophil expression may be able to distinguish patients with IA, biomarker development from expression data is more complex than the identification of differentially expressed genes. Average expression differences between patients with and without IA cannot themselves be used to predict aneurysm on an individual basis . Rather, classification algorithms that use discrete expression levels of informative transcripts to predict the presence of IA from individual samples are required.
In this study, we sought to test the feasibility of creating classification models of unruptured IAs based on RNA expression of circulating neutrophils. We recruited additional patients with and without unruptured IAs (confirmed on angiography), isolated peripheral blood neutrophils, and performed next-generation RNA sequencing to obtain the neutrophil transcriptomes. Differential expression analysis was used to identify highly significantly differentially expressed transcripts as features for model development. As there was no precedent for the specific type of algorithm best suited for IA detection from neutrophil expression differences, we applied four widely used supervised machine-learning approaches to select a classification model most fitting for our data. Trained classification models were then validated in an independent testing dataset of neutrophil transcriptomes. Furthermore, real-time quantitative polymerase chain reaction (qPCR) was used to corroborate expression differences of model transcripts in neutrophil samples from a new, separate validation cohort of patients with and without IA. Results from this study could lay the groundwork for future, larger efforts towards developing a blood-based IA diagnostic.
This study was approved by the University at Buffalo Health Sciences Institutional Review Board (Study No. 030-474433). Methods were carried out in accordance with the approved protocol and written informed consent was obtained from all subjects. We included individuals who were older than 18 years, spoke English, and had not received previous treatment for IA. We excluded individuals with potentially altered immune systems; including patients who were pregnant, had recently undergone invasive surgery, were undergoing chemotherapy, had a body temperature above 37.78 °C (100 °F), had received solid organ transplants, had autoimmune diseases, and those who were taking prednisone or any other immunomodulating drugs.
Between December 2013 to May 2016, 106 peripheral blood samples were collected from patients undergoing cerebral digital subtraction angiography (DSA) at Gates Vascular Institute in Buffalo, New York: 51 patients had a positive IA diagnosis and 55 had a negative IA diagnosis (controls). DSA imaging was used to confirm IA diagnosis by either positive or negative angiographic presence of IA. Patient medical electronic records were also collected. Reasons for the patients to receive DSA included confirmation of findings from noninvasive imaging of the presence of unruptured IAs, vascular malformations, or carotid stenosis, or follow-up noninvasive imaging of previously detected IAs. Prior to RNA expression analysis, we further excluded patients with other known cerebrovascular malformations or extracranial aneurysms, including abdominal aortic aneurysms. The presence of other cerebrovascular malformations or extracranial aneurysms was recorded from both the patient’s operative report following DSA and their recorded medical history.
During the DSA procedure, 16 mL of blood was drawn from the access catheter in the femoral artery and transferred into two 8 mL, citrated, cell preparation tubes (BD, Franklin Lakes, NJ). Neutrophils were isolated within 1 h of peripheral blood collection, as described elsewhere . Cell preparation tubes were centrifuged at 1700×g for 25 min to separate erythrocytes and neutrophils from mononuclear cells and plasma in the peripheral blood samples via a Ficoll density gradient. Erythrocytes and neutrophils were collected into a 3 mL syringe. Following hypotonic lysis of red blood cells, neutrophils were isolated by centrifugation at 400×g for 10 min and disrupted and stored in TRIzol reagent (Life Technologies, Carlsbad, CA) at − 80 °C until further processing. Neutrophils isolated in this fashion are more than 98% CD66b+ by flow cytometry and contain no contaminating CD14+ monocytes .
Neutrophil RNA was extracted using TRIzol, according to the manufacturer’s instructions. Trace DNA was removed by DNase I (Life Technologies, Carlsbad, CA) treatment. RNA was purified using the RNeasy MinElute Cleanup Kit (Qiagen, Venlo, Limburg, Netherlands) and suspended in RNase-free water. The purity and concentration of RNA in each sample were measured by absorbance at 260 nm on a NanoDrop 2000 spectrophotometer (Thermo Scientific, Waltham, MA), and 200–400 ng of RNA was sent to our university’s Next-Generation Sequencing and Expression Analysis Core facility for further quality control. Precise RNA concentration was measured at the core facility via the Quant-iT RiboGreen Assay (Invitrogen, Carlsbad, CA) with a TBS-380 Fluorometer (Promega, Madison, WI). The quality of the RNA samples was measured with an Agilent 2100 BioAnalyzer RNA 6000 Pico Chip (Agilent, Las Vegas, NV). RNA samples of acceptable purity (260/280 ratio of ≥ 1.9) and integrity (RIN ≥ 5.0) were considered for RNA sequencing.
RNA libraries were constructed using the Illumina TruSeq RNA Library Preparation Kit (Illumina, San Diego, CA). All samples were subjected to 50-cycle, single-read sequencing in a HiSeq 2500 system (Illumina) and demultiplexed using Bcl2Fastq v126.96.36.199 (Illumina). Gene expression analysis was carried out using the Tuxedo Suite [14,15,16,17]. For each sample, short RNA fragment data in FASTQ format was compiled and aligned to the human reference genome (human genome 19—hg19) using TopHat v2.1.13 . To evaluate the quality of RNA sequencing, we performed quality control analysis using FASTQC  and visualized and compared the aggregate quality control data using MultiQC .
Transcript expression levels were calculated from counts using transcripts per million (TPM) normalization for comparison of RNA levels between samples. Since samples were processed in two batches, we performed batch effect correction using ComBat under the default settings in R [20, 21]. This was performed on expression data for all transcripts with an average TPM > 1.0 in at least one of the two groups (see Additional file 1: Table S1 for batch information).
Differential expression analysis
Prior to differential expression analysis, neutrophil transcriptomes were randomly divided into two cohorts. Based on standard convention [22, 23], and to maximize learning while retaining a substantial testing group, 3/4 (75%) of the samples were allocated to a training cohort and 1/4 (25%) was allocated to a testing cohort, each containing half IA and half control samples. Although the investigators could not be blinded to sample class, randomization was performed to unbiasedly allocate samples to the training and testing cohorts. Differential gene expression analysis in the training cohort was carried out using F statistics to assess differential variation in the mean on a transcript-by-transcript basis [24,25,26,27]. Multiple testing correction was performed by using the John Storey method , and q-values were reported for each transcript. Transcripts were considered significantly differentially expressed at an FDR-adjusted p-value (q-value) < 0.05.
We performed gene ontology (GO) term enrichment analysis using the open source Gene Ontology enRIchment anaLysis and visuaLizAtion tool (GORILLA) on all differentially expressed transcripts (q < 0.05) . This was done using a background gene list of previously published neutrophil RNA expression patterns (average fragments per kilo base of transcript per million mapped reads, FPKM > 1.0) of three healthy individuals, described elsewhere . This tool identified GO terms that are enriched in genes with increased or decreased expression in IA compared to the background neutrophil expression using standard hypergeometric statistics. We reported associated GO processes and functions if the enrichment FDR-adjusted p-value (q-value) was < 0.20 (20% FDR).
Feature selection for classification model development
Prior to model training, the set of differentially expressed transcripts was reduced by filtering. We retained only transcripts with an FDR < 0.05 and absolute fold-change ≥ 1.5. To visualize how those transcripts separated samples from patients with and without IAs, we performed principal component analysis (PCA) in R using the prcomp package under the default settings . We also performed post hoc power estimation following a method by Hart et al.  for 15 samples in each group, with an α = 0.05, and a coefficient of variation and counts per million mapped reads of 0.40 and 38, respectively .
Using the selected transcripts, we trained classification models using MATLAB Statistics and Machine Learning Toolbox (MathWorks, Natick, MA) and R bioconductor (https://www.bioconductor.org/). Specifically, we used four algorithms that have been successfully used for disease classification from gene expression data . These machine-learning methods included cosine nearest neighbors (cosine NN) classification , diagonal linear discriminant analysis (DLDA) , nearest shrunken centroids (NSC) classification , and support vector machines (SVM) . Each method was applied to the training cohort separately and evaluated with a leave-one-out (LOO) cross-validation to estimate model performance and prevent overfitting.
The k-nearest neighbor method  with a cosine metric (cosine NN) was employed. The number of neighbors, k, was set as 5 for cosine NN. The resulting model classified test samples by calculating their distance to each training sample. The test sample labels were predicted by choosing the class that was most common among their k-nearest neighbors.
We trained a classifier using diagonal LDA (DLDA), as described elsewhere . This method seeks the linear combination of transcripts that best separates two classes using a diagonal covariance matrix. The linear model coefficients associated with transcripts (discriminant scores) relayed the importance of each transcript to the prediction model . Classification was performed by projecting a test sample onto the maximally separating direction that was determined by discriminant scores and calculating the corresponding posterior probability of IA.
Nearest centroids classification
We used a modification of the nearest centroids technique, called NSC . This method calculates class-specific centroids (standard deviation normalized averages) for each transcript and refines them by eliminating those with variable expression. Classification was performed by comparing the expression of the included model transcripts with the centroids of the two classes and assigning it to the class that it was closest to in squared distance .
The most complex classification algorithm we implemented was SVM . To separate the binary labeled training samples, SVM finds a hyperplane that is maximally distant from samples of either class. A linear kernel was used in model creation. The resulting model classified test samples by mapping them to a higher-dimensional space and making decisions based on their signed distance to the hyperplane.
Model assessment in the training cohort
The performance of each model in the training cohort was estimated using the results of the LOO cross-validation. The model classifications were compared to each patient’s clinical diagnosis from imaging, and the true positives (TP), true negatives (TN), false positives (FP), and false negatives (FN) were counted. Each model’s performance was first assessed by calculating the model’s sensitivity, specificity, and accuracy, as follows:
Based on model predictions, we created receiver operating characteristic (ROC) curves and calculated the area under the ROC curve (AUC) to assess model performance .
Validation of the models in an independent testing cohort
Classification models were independently tested on transcriptomes from the testing cohort. TPM values of these model features were input into the models for classification of IA presence. The classification results were compared to clinical diagnoses to calculate the true sensitivity, specificity, and accuracy for each model. ROC curves were constructed and AUCs were used to assess the performance of each classifier .
Cross-validation over all samples
Because the models were fit using data points from a randomly selected training dataset (n = 30), selection bias may introduce inconsistency in model predictions. To increase the prediction reliability of the models and to create algorithms more generalizable to a broader population, we implemented LOO cross-validation using the expression levels of the 26 selected transcripts from all 40 patients for each model. The LOO cross-validation method essentially retrained the models in 40 different training sets consisting of 39 samples and performed testing on the remaining sample. As before, classification results were used to calculate sensitivity, specificity, and accuracy for each model, as well as to find the AUC of the ROC curve for each modified classifier.
Positive and negative predictive values of the models
On the basis of the cross-validation results over all samples, we assessed the predictive value of the classification models by calculating their positive predictive values (PPV) and negative predictive values (NPV) . PPVs and NPVs were estimated using the following formulas based on Bayes’ theorem [43, 44]:
Validation of expression differences by qPCR in an independent validation cohort
To validate expression differences in the 26 genes in our models, quantitative polymerase chain reaction (qPCR) was performed. Due to limitations in mRNA volume, qPCR was performed on seven model transcripts in 10 additional patients (independent cohort of 5 with IA and 5 controls), as described previously . In brief, oligonucleotide primers were designed with a 60 °C melting temperature and a length of 15–25 nucleotides to produce PCR products with lengths of 50–250 base pairs using Primer3 software  and Primer BLAST (NCBI, Bethesda, MD). The replication efficiency of each primer set was tested by performing qPCR on serial dilutions of cDNA samples (primer sequences, annealing temperatures, efficiencies, and product lengths are shown in Additional file 1: Table S2).
For reverse transcription, first-strand cDNA was generated from total RNA using OmniScript Reverse Transcriptase kit (Qiagen, Venlo, Limburg, Netherlands) according to the manufacturer’s directions. qPCR was run with 10 ng of cDNA in 25 µL reactions in triplicate in Bio-Rad CFX Connect (Bio-Rad, Hercules, California) using ABI SYBR Green I Master Mix (Applied Biosystems, Foster City, California) and gene-specific primers at a concentration of 0.02 μM each. The temperature profile consisted of an initial step of 95 °C for 10 min, followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min, and then a final melting curve analysis from 60 °C to 95 °C over 20 min.
Gene-specific amplification was demonstrated by a single peak using the Bio-Rad dissociation melt curve. Samples were normalized based on GAPDH, 18s rRNA, and GPI expression, which were run in parallel reactions to the genes of interest. These values were used to calculate average fold-change between the two groups using the 2−ΔΔCt method . These values were calculated for each housekeeping gene and averaged. Average fold-change in gene expression measured by qPCR data in the new cohort was then compared to the fold-change calculated from RNA sequencing in the training cohort.
Influence of IA size on the 26 classifier transcripts
To determine if aneurysm size could affect classification results, we dichotomized the entire IA cohort into a “small” (< 5 mm) group and a “large” (≥ 5 mm) group and analyzed expression differences of the 26 classifier transcripts in samples from patients with IA in the two groups separately. The 5 mm aneurysm size cutoff was based on data reported by the PHASES study , which pooled an analysis of six longitudinal investigations [3, 52,53,54,55,56] and found that aneurysms < 5 mm were less likely to rupture than those with larger diameters; IAs between 7 and 10 mm had 2.7 times greater risk of rupture than small IAs (< 5 mm); and IAs > 20 mm had 14.3 times greater risk . We investigated fold-changes in the 26 genes between the “small” IA group and the entire control group (n = 20) and between the “large” IA group and the entire control group (n = 20) to determine if aneurysm size affected their expression pattern.
During the study period, we collected 106 blood samples (51 from patients with IA, 55 from control subjects) as well as angiographic images and medical records data from individuals undergoing cerebral DSA. Of the blood samples collected, 43 (20 from IA patients, 23 from controls) met our criteria and also had neutrophil RNA of sufficient quality and volume for sequencing. A total of 40 patients (20 with IA and 20 controls) were then chosen and randomly divided into a 30-patient training cohort (n = 15 IA and n = 15 control) and a 10-patient testing cohort (n = 5 IA and n = 5 control). Characteristics of the two cohorts are provided in Table 1. These samples were of sufficient quality and had an average 260/280 ratio of 2.02 (range 1.90–2.12) and an average RNA integrity number (RIN) of 6.88 (range 5.2–8.2) (Additional file 1: Table S3). Patients with IAs had only saccular aneurysms that ranged in size (greatest diameter) from 1 to 19 mm. Five patients with IA had multiple aneurysms (Additional file 1: Table S4). A portion of these samples (n = 22) had been previously analyzed in our aforementioned study that investigated neutrophil expression differences between patients with and without IA .
Differential RNA expression in neutrophils from patients with IA vs. controls
RNA sequencing data were used to identify differentially expressed neutrophil transcripts between the 15 patients with IA and 15 controls in the training cohort. Overall, our sequencing experiments had an average of 53.84 million sequences per sample and a 95.4% read mapping rate (% aligned) (Additional file 1: Table S5). The volcano plot in Fig. 1a shows neutrophil expression differences between IA patients and controls in terms of average fold-change in expression and significance level. From 12,775 different transcripts with average TPM > 1.0 across samples in either group, differential expression analysis identified 95 transcripts that were significantly differentially expressed (q < 0.05). Gene set enrichment analysis performed using these 95 differentially expressed transcripts showed that genes with higher levels in the IA group were involved in defense response, leukocyte activation, stem cell maintenance, maintenance of cell number, cell activation, and stem cell development (Table 2). Genes with lower levels in IAs were involved in regulation of glutathione and tetrapyrrole binding.
Selected transcripts for model training
Prior to model training, we performed feature selection by filtering to identify disease-related transcripts and reduce the data dimensionality to facilitate downstream analysis. Our statistical criteria of false discovery rate (FDR) < 0.05 and an absolute fold-change ≥ 1.5 resulted in the retention of the 26 transcripts that are shown in red in Fig. 1a and listed in Table 3. The PCA in Fig. 1b shows that these 26 transcripts could generally discriminate patients with IAs from the controls. The top three principal components represented 47.8% of the variation: PC1 contained 22.4% variation, PC2 contained 15.3% variation, and PC3 contained 10.1% variation. Overall, 60% of the samples from the IA patients and 80% of those from controls could be grouped together by PCA. Furthermore, our post hoc power analysis estimated that in the 15 vs. 15 training dataset a power of 0.78 was achieved in detecting expression differences of at least 1.5 fold at an α = 0.05. Thus, our feature selection criteria resulting in the identification of the 26 transcripts had a power of > 0.78.
Classification models of IA have high performance in training and testing datasets
Using the expression of these 26 transcripts, we trained the aforementioned four classification models (cosine NN, DLDA, NSC, and SVM). Figure 1c shows the sensitivity, specificity, and accuracy of the models, which were estimated from LOO cross-validation. There was moderate performance by each classification method, with accuracies that ranged from 0.50 to 0.73. Evaluation by ROC curve analysis showed a range in AUCs from 0.54 to 0.73 (Fig. 1d) across all methods. In this dataset, DLDA performed the best, with a sensitivity of 0.67, a specificity of 0.80, an accuracy of 0.73, and an AUC of 0.73.
To independently validate the models, we implemented them in neutrophil transcriptomes from 10 patients in the testing cohort. The PCA in Fig. 2a shows that the 26 transcripts could discriminate patients with IAs from controls in the testing cohort as well. Overall, 100% of the samples from IA patients and 80% of those from controls could be grouped together by PCA. In the testing cohort, the models predicted aneurysm status with a range in accuracy from 0.60 to 0.90 (Fig. 2b). The ROC analysis in Fig. 2c shows that model AUCs ranged from 0.62 to 0.80. In this cohort, the DLDA classification model again performed the best, with a sensitivity of 0.80, specificity of 1.0, an accuracy of 0.90, and an AUC of 0.80.
Cross-validation to increase model reliability
To increase model reliability, we applied LOO cross-validation using all patient transcriptomes and refit them in all 40 datasets (all combinations of 39 training samples and 1 testing sample). This analysis revealed that model accuracy ranged from 0.63 to 0.80 (Fig. 3a) and model AUCs ranged from 0.68 to 0.84 (Fig. 3b). Again, the DLDA model performed the best, with a sensitivity of 0.65, specificity of 0.95, accuracy of 0.80, and an AUC of 0.84.
Models have high negative predictive value
Given their range of performance, we wanted to know how useful the models would be at detecting IA. However, their value would be inherently influenced by the prevalence of IA in a given target population. To estimate this, we plotted the positive predictive value (PPV) and negative predictive value (NPV) for each model (Fig. 3c, d) using the sensitivity and specificity reported after the LOO cross-validation in all datasets. The rate of aneurysm prevalence found in the literature is between 3.2%  to 7% , and is highlighted in the blue region on the graphs in Fig. 3c, d. At 5% prevalence, the PPV of the models ranged from 0.10 to 0.41 and NPV ranged from 0.96 to 0.98. The DLDA classifier had the highest PPV (0.41) and NPV (0.98).
Independent validation of expression differences by RT-qPCR
We performed a corroboration study to determine if the differential expression of seven model genes could be detected in a new population of IA and control patients. We used samples from 10 additional patients (5 with IA and 5 controls) from which we collected neutrophil RNA but did not sequence (see Additional file 1: Table S6 for patient information for this cohort and Additional file 1: Table S7 for aneurysm information from the IA patients in this cohort). These samples were used for quantitative qPCR analysis of CD177, CYP1B1, ARMC12, OLAH, CD163, G0S2, and FCRL5, which were selected because they were also differentially expressed in our previous study . Figure 4 shows the qPCR results of this corroborative study in comparison with expression differences obtained from RNA sequencing in the training cohort. Six of the seven genes demonstrated average fold-change in the same direction and of similar magnitude to those in the original cohort. This provides evidence that in ~ 86% (6 of 7) of the tested transcripts the expression differences between patients with IA and controls is consistent. This result may serve as the basis for further development of a qPCR-based diagnostic.
Signals of the classifier genes are stronger for larger IAs
To determine if there was a correlation between aneurysm size and differential expression of the 26 neutrophil classifier transcripts, we separately compared expression levels in patients with “small” and “large” IAs to those in all control patients (n = 20). We used the 5 mm aneurysm size cutoff , which dichotomized the IA cohort (n = 20) into a 10-patient “small” group and a 10-patient “large” group. Figure 5 shows the fold-changes between the “small” group and control samples (green) and the “large” group and control samples (orange) for each of the 26 classifier transcripts. These are compared to the fold-changes between aneurysm and control samples in the entire training cohort as a baseline (solid black line). Expression changes were more pronounced in both the positive and negative direction in patients with larger IAs. On average, the fold-change was 24% greater than the baseline for the “large” group, but 35% lower than baseline for the “small” group.
Neutrophils play a role in the progressive inflammation that typifies IAs . We hypothesized that gene expression patterns in circulating neutrophils could indicate the presence of aneurysm. Previously, we found that patients with unruptured IAs and aneurysm-free controls had significant RNA expression differences in circulating neutrophils . In the present study, we tested the feasibility of developing biomarkers using neutrophil RNA expression levels from blood samples to predict IA presence. Using RNA expression profiling in circulating neutrophils, we identified 26 transcripts that were highly associated with the presence of IA. Machine-learning algorithms were then implemented to develop classification models to test whether these 26 transcripts could predict the presence of an aneurysm.
Classification models of IA using neutrophil RNA expression
To our knowledge, this is the first to demonstrate IA prediction from RNA expression patterns in the blood. The prediction models we trained had a classification accuracy of up to 90% in the test dataset, a level that indicates promise for further investigation of RNA expression biomarkers for IA. Overall, classification by DLDA achieved the best performance in our data. This model consistently had the highest accuracy and AUC over multiple analyses, including cross-validation during model training (accuracy = 0.73, AUC = 0.73), independent model testing (accuracy = 0.90, AUC = 0.80) and cross-validation across the entire dataset (accuracy = 0.80, AUC = 0.84). See Additional file 1: Table S8 for a summary of the performance of all four models.
We suspect that DLDA outperformed the other methods because it best accounted for potential intersample variability of the 26 transcripts. Modeling techniques that broadly survey patterns of gene expression may afford better IA classification [58, 59]. The DLDA method does this by (1) ranking transcripts by importance, giving more weight to the most consistently informative transcripts (unlike non-parametric approaches such as NN); and (2) using information from all transcripts to project test samples to the direction which best separates the classes. In this way, a linear combination of transcripts, which may be individually inconsistent, can generate a stable IA prediction and accommodate for potential intersample variability. Additionally, ignoring correlations between genes as DLDA does, provided a simple model and produced lower misclassification rates than more sophisticated classifiers, such as SVM.
In the current study, classifiers were developed based on 30 transcriptomes that were randomly selected from all available data (n = 40). Randomization was used so we could test the viability of IA prediction in patients who have potentially confounding covariates (comorbidities and demographics). Table 1 shows differences in smoking history between the IA and control groups in the training cohort, which may reflect an established association between the presence of an IA and smoking . It is worth noting that this study was designed differently from our previous investigation . There we identified an 82-transcript expression signature of IA by transcriptome profiling of a cohort-controlled group of patients, whereas in this study the randomly selected training cohort was not cohort-controlled. Yet, even with this difference, 10 of the 26 classifier transcripts (38%) were also a part of the previously discovered 82-transcript signature. These genes include CYP1B1, CD177, ARMC12, OLAH, CD163, ADTRP, VWA8, G0S2, FCRL5, and C1orf226. Notably, in qPCR validation on seven of these genes, six of them (CYP1B1, ARMC12, OLAH, CD163, G0S2, and FCRL5) showed consistent expression differences. These six transcripts may warrant further investigation as they may be most important for IA detection.
Biological role of classifier transcripts
The natural history of IA has been characterized by mounting inflammatory responses accompanied by progressive degradation of the aneurysmal wall [61, 62]. This begins during aneurysm initiation, in which a combination of risk factors and hemodynamic stresses elicit pro-inflammatory changes in smooth muscle cells that lead to overproduction of matrix metalloproteinases (MMPs) that degrade the arterial extracellular matrix [61, 62]. Once the aneurysmal sac forms, recirculating blood in the IA is conducive to inflammatory cell infiltration into the wall, which is also assisted by an increase of chemokines and cytokines in both the aneurysm wall and in the plasma within the aneurysm [63, 64]. Recruited inflammatory infiltrates also produce MMPs that continue to degrade the aneurysm wall and advance its growth and rupture [57, 62], which can ultimately occur when the wall is weakened to the point that it can no longer contain the force of the blood pressure . This is demonstrated by histological analyses and gene expression studies of human aneurysmal tissues, which have found increased matrix degradation proteins, inflammatory processes, and inflammatory cytokines and chemoattractant proteins in the walls of aneurysms [65, 66].
Despite being the most abundant circulating immune cell, the role that neutrophils play in IA pathophysiology is relatively unknown. However, since neutrophils are recruited to sites of injury to coordinate the inflammatory response and attract inflammatory cells (monocytes) during other vascular pathologies , we initially suspected they may play a similar role in IA. Studies suggest that neutrophils reside in intraluminal thrombi that have formed on the wall of the aneurysmal sac [68, 69]. There, besides secreting MMP-9, activated neutrophils can release myeloperoxidase (MPO) and neutrophil gelatinase associated lipocalin (NGAL) that can indirectly promote extracellular matrix degradation and cytotoxicity. Elevated levels of MPO, a peroxidase enzyme secreted during degranulation, provoke pro-inflammatory cell signaling and oxidative stress via increased production of reactive oxygen species . Increased NGAL protects MMP-9 from degradation, thereby increasing its activity and promoting wall degeneration . Interestingly, levels of MPO and NGAL have been shown to be elevated in the peripheral blood of patients with IAs [64, 71], which can act in an autocrine manner to promote activation of circulating neutrophils [72, 73]. In this study, however, we did not observe significantly higher expression of the MPO or NGAL genes in circulating neutrophils, which suggests that the source of these circulating proteins could be the wall itself or other blood cells.
Our data shows that circulating neutrophils from patients with IA are peripherally activated. From gene ontology analysis on all 95 differentially expressed genes (q < 0.05), we found that they have dysregulated inflammatory and defense responses, and increased signaling and response to stimuli. Increased IL-1 signal transduction through receptors IL1R1 and IL1R2 has been shown to play a major role in neutrophil activation [74,75,76]. Increased neutrophil activation was also evidenced through greater expression levels of several CD antigens, namely CD36, CD99L2, CD163, and CD177. Specifically, CD177 is a marker of neutrophil activation that plays a role in migration , and CD99L2 is involved in neutrophil recruitment to inflamed tissues . These findings are consistent with the results of our previous cohort-controlled study , which also showed increased peripheral leukocyte activation in neutrophils from IA patients.
In the subset of the 26 classifier transcripts, neutrophil activation was reflected through the roles of five genes (CD177, IL18R1, PVRL2, PDE9A, and PTGES). Like CD177, IL18R1 has been shown to be involved in neutrophil activation as well as migration via IL-18 signaling . Nectin-2 (PVRL2), a membrane glycoprotein, is involved in cell adhesion , and has been shown to have increased expression in the aneurysm wall tissue . Similarly, PDE9A (a cGMP-specific phosphodiesterase) is also involved in cell adhesion [80, 82] and, as demonstrated by Li et al. , is regulated by two of the most active transcription factors in the IA tissue. Furthermore, lower PTGES expression may be partially responsible for increasing the lifespan of neutrophils, because it is involved in p53-induced apoptosis . We suspect that capturing neutrophil activation responses involved in IA is the reason why the 26-transcript biomarker was able to detect IA.
Besides these five genes, nine other transcripts (CD163, TGS1, CYP26B1, ADTRP, OCLN, OLAH, C1QL1, FCRL5, and IGSF23) in the 26-transcript classifier reflect complex inflammatory processes, albeit not specifically attributed to neutrophil activation. For example, CD163, which is abundant in the walls of IAs (but typically associated with macrophages [85, 86]), has been shown to be increased in neutrophils during sepsis  and thus could be contributing to vascular inflammation in IA. Expression differences of other transcripts, like TGS1 and CYP26B1 (that are differentially expressed in tuberculosis  and juvenile idiopathic arthritis , respectively) could be related to neutrophil responses to intravascular perturbations during IA. Other transcripts—such as ADTRP (expressed by macrophages in coronary artery plaques) , OCLN (increased in activated T-lymphocytes and in whole blood during sepsis) [90, 91], OLAH (increased in peripheral blood mononuclear cells during non-small cell lung cancer) , C1QL1 (a complement component decreed in glioblastoma multiform) , FCRL5 (an immunoglobulin receptor that regulates B cell response to antigen) , and IGSF23 (decreased during the inflammatory response associated with mycoestrogen exposure) —are involved in inflammation but have not been reported to be differentially expressed in neutrophils. The roles of the remaining model transcripts (e.g., C1orf226, LOC100506229, MTRNR2L10) in neutrophils are unknown. Further study into the roles of these transcripts in IA may be warranted, as they could represent novel predictive targets in neutrophil RNA expression.
Taken together, our results suggest that circulating neutrophils are peripherally activated in patients with IA, which leads to a change in their RNA expression profile. We postulate this activation could be caused by contact with inflamed aneurysm tissue that often contains denuded regions and plaque or thrombus [96,97,98]. Alternatively, the activation may be caused by chemokines and cytokines released from the aneurysm. Chalouhi et al.  demonstrated that blood inside IAs contain increased concentrations of the chemokines MCP-1, RANTES, MIG, IP-10, and exotoxin, and chemoattractant cytokines, including IL-8 and IL-17. Either of the above two scenarios may explain why expression differences of the 26 classifier transcripts were more pronounced in patients with larger IAs, since larger IAs provide greater luminal surface area for either contact or release of inflammatory mediators. It would be interesting to conduct a longitudinal study of patients with growing aneurysms to ascertain the relationship between aneurysm development and the effect on gene expression in circulating neutrophils.
Due to our small sample size, the results in this study are rather preliminary. However, we have confidence in the discovered classifier transcripts for the following reasons. (a) The classifier identified patients with IA in an independent testing cohort with 90% accuracy. (b) qPCR confirmed expression differences in an independent validation cohort. (c) Our post hoc power analysis demonstrated > 0.78 power. In the future, we could further increase reliability in the model transcripts by decreasing the number of features in the data and increasing sample size.
Secondly, our classifier was derived from a population of patients who had different rates of demographic factors and comorbidities between aneurysm and control patients. It is unclear whether the presence of these confounding factors contributes to the differential neutrophil expression we detected. We are currently conducting research on larger cohorts by including multiple control groups with different vascular pathologies (including extracranial aneurysms such as abdominal aortic aneurysm) and immunological conditions (both of which were excluded in the current study) to narrow down transcripts specific to IA. Finally, although we collected basic demographics and information about comorbidities, including hypertension and diabetes, there could be others factors, such as the presence of immune-metabolism mediators in the blood that could affect gene expression in circulating neutrophils. Efforts should be made to collect higher fidelity patient health metadata.
We have shown for the first time the feasibility of using blood-based biomarkers to detect unruptured IA. A model consisting of 26 transcripts predicted IA presence with an accuracy of 0.80 and an AUC of 0.84. This biomarker reflects gene expression changes due in part to activation of circulating neutrophils associated with IA. Pending confirmation in larger cohorts, these results suggest an exciting potential to develop blood-based IA tests to screen for IA in high-risk populations.
area under the receiver operating characteristic curve
- cosine NN:
cosine nearest neighbors
diagonal linear discriminant analysis
digital subtraction angiography
false discovery rate
fragments per kilobase million
Gene Ontology enRIchment anaLysis and visuaLizAtion tool
linear discriminant analysis
negative predictive value
nearest shrunken centroids
principal component analysis
positive predictive value
quantitative polymerase chain reaction
receiver operator characteristic
support vector machines
transcripts per million
de Rooij NK, Linn FH, van der Plas JA, Algra A, Rinkel GJ. Incidence of subarachnoid haemorrhage: a systematic review with emphasis on region, age, gender and time trends. J Neurol Neurosurg Psychiatry. 2007;78:1365–72.
Greving JP, Rinkel GJ, Buskens E, Algra A. Cost-effectiveness of preventive treatment of intracranial aneurysms: new data and uncertainties. Neurology. 2009;73:258–65.
Juvela S. Treatment options of unruptured intracranial aneurysms. Stroke. 2004;35:372–4.
Vega C, Kwoon JV, Lavine SD. Intracranial aneurysms: current evidence and clinical practice. Am Fam Physician. 2002;66:601–8.
Keedy A. An overview of intracranial aneurysms. Mcgill J Med. 2006;9:141–6.
Aziz H, Zaas A, Ginsburg GS. Peripheral blood gene expression profiling for cardiovascular disease assessment. Genom Med. 2007;1:105–12.
Jin H, Li C, Ge H, Jiang Y, Li Y. Circulating microRNA: a novel potential biomarker for early diagnosis of intracranial aneurysm rupture a case control study. J Transl Med. 2013;11:296.
Liu L, Gritz D, Parent CA. PKCβII acts downstream of chemoattractant receptors and mTORC2 to regulate cAMP production and myosin II activity in neutrophils. Mol Biol Cell. 2014;25:1446–57.
Sabatino G, Rigante L, Minella D, Novelli G, Della Pepa GM, Esposito G, Albanese A, Maira G, Marchese E. Transcriptional profile characterization for the identification of peripheral blood biomarkers in patients with cerebral aneurysms. J Biol Regul Homeost Agents. 2013;27:729–38.
Tutino VM, Poppenberg KE, Jiang K, Jarvis JN, Sun Y, Sonig A, Siddiqui AH, Snyder KV, Levy EI, Kolega J, Meng H. Circulating neutrophil transcriptome may reveal intracranial aneurysm signature. PLoS ONE. 2018;13:e0191407.
Deng X, Campagne F. Introduction to the development and validation of predictive biomarker models from high-throughput data sets. Methods Mol Biol. 2010;620:435–70.
Jiang K, Sun X, Chen Y, Shen Y, Jarvis JN. RNA sequencing from human neutrophils reveals distinct transcriptional differences associated with chronic inflammatory states. BMC Med Genom. 2015;8:55.
Jarvis JN, Dozmorov I, Jiang K, Frank MB, Szodoray P, Alex P, Centola M. Novel approaches to gene expression analysis of active polyarticular juvenile rheumatoid arthritis. Arthritis Res Ther. 2004;6:R15–32.
Ghosh S, Chan CK. Analysis of RNA-seq data using TopHat and Cufflinks. Methods Mol Biol. 2016;1374:339–61.
Kulahoglu C, Brautigam A. Quantitative transcriptome analysis using RNA-seq. Methods Mol Biol. 2014;1158:71–91.
Pollier J, Rombauts S, Goossens A. Analysis of RNA-seq data with TopHat and Cufflinks for genome-wide expression analysis of jasmonate-treated plants and plant cultures. Methods Mol Biol. 2013;1011:305–15.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7:562–78.
Andrews S (2010) FASTQC. A quality control tool for high throughput sequence data. http://www.bioinformatics.babraham.ac.uk/projects/fastqc/.
Ewels P, Magnusson M, Lundin S, Kaller M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics. 2016;32:3047–8.
Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007;8:118–27.
Leek JT, Johnson WE, Parker HS, Jaffe AE, Storey JD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics. 2012;28:882–3.
Hamori S, Kawai M, Kume T, Murakami Y, Watanabe C. Ensemble learning or deep learning? Application to default risk analysis. J Risk Financ Manag. 2018;11:12.
Paliwal N, Jaiswal P, Tutino VM, Shallwani H, Davies JM, Siddiqui AH, Rai R, Meng H. Outcome prediction of intracranial aneurysm treatment by flow diverters using machine learning. Neurosurg Focus. 2018;45:10.
Auer PL, Doerge RW. Statistical design and analysis of RNA sequencing data. Genetics. 2010;185:405–16.
Cui X, Hwang JT, Qiu J, Blades NJ, Churchill GA. Improved statistical tests for differential gene expression by shrinking variance components estimates. Biostatistics. 2005;6:59–75.
Reeb PD, Steibel JP. Evaluating statistical analysis models for RNA sequencing experiments. Front Genet. 2013;4:178.
Urasoko Y, He XJ, Masao T, Kinoshita Y, Edamoto H, Hatayama K, Asano Y, Tamura K, Mochizuki M. Changes in blood parameters and the expression of coagulation-related genes in lactating Sprague–Dawley rats. J Am Assoc Lab Anim Sci. 2012;51:144–9.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a pratical and powerful approach to multiple testing. J R Stat Soc Ser B (Methodol). 1995;57:289–300.
Eden E, Navon R, Steinfeld I, Lipson D, Yakhini Z. GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinform. 2009;10:48.
Everitt BS, Hothorn T. A handbook of statistical analyses using R. http://www.ecostat.unical.it/tarsitano/Didattica/LabStat2/Everitt.pdf. 2005.
Hart SN, Therneau TM, Zhang Y, Poland GA, Kocher JP. Calculating sample size estimates for RNA sequencing data. J Comput Biol. 2013;20:970–8.
Ye J, Li T, Xiong T, Janardan R. Using uncorrelated discriminant analysis for tissue classification with gene expression data. IEEE/ACM Trans Comput Biol Bioinform. 2004;1:181–90.
Dudoit S, Fridlyand J, Speed TP. Comparison of discrimination methods for the classification of tumors using gene expression data. J Am Stat Assoc. 2002;97:77–87.
Golub TR, Slonim DK, Tamayo P, Huard C, Gaasenbeek M, Mesirov JP, Coller H, Loh ML, Downing JR, Caligiuri MA, et al. Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. Science. 1999;286:531–7.
Tibshirani R, Hastie T, Narasimhan B, Chu G. Diagnosis of multiple cancer types by shrunken centroids of gene expression. Proc Natl Acad Sci USA. 2002;99:6567–72.
Lee Y, Lee CK. Classification of multiple cancer types by multicategory support vector machines using gene expression data. Bioinformatics. 2003;19:1132–9.
Dudoit S, Fridlyand J. A prediction-based resampling method for estimating the number of clusters in a dataset. Genome Biol. 2002;3:Research0036.
Huerta EB, Duval B, Hao J-K. Gene selection for microarray data by a LDA-based genetic algorithm. In: LERIA, Universit´e d’Angers. 2008. pp. 252–63.
Furey TS, Cristianini N, Duffy N, Bednarski DW, Schummer M, Haussler D. Support vector machine classification and validation of cancer tissue samples using microarray expression data. Bioinformatics. 2000;16:906–14.
Florkowski CM. Sensitivity, specificity, receiver-operating characteristic (ROC) curves and likelihood ratios: communicating the performance of diagnostic tests. Clin Biochem Rev. 2008;29(Suppl 1):S83–7.
Hanley JA, McNeil BJ. A method of comparing the areas under receiver operating characteristic curves derived from the same cases. Radiology. 1983;148:839–43.
Ray P, Le Manach Y, Riou B, Houle TT. Statistical evaluation of a biomarker. Anesthesiology. 2010;112:1023–40.
Mercaldo ND, Lau KF, Zhou XH. Confidence intervals for predictive values with an emphasis to case-control studies. Stat Med. 2007;26:2170–83.
Habib A, Alalyani M, Musa IH, Almutheibi MS. Brief review on sensitivity, specificity and predictivities. IOSR J Dent Med Sci. 2015;14:64–8.
Harada K, Fukuyama K, Shirouzu T, Ichinose M, Fujimura H, Kakumoto K, Yamanaga Y. Prevalence of unruptured intracranial aneurysms in healthy asymptomatic Japanese adults: differences in gender and age. Acta Neurochir (Wien). 2013;155:2037–43.
Li MH, Chen SW, Li YD, Chen YC, Cheng YS, Hu DJ, Tan HQ, Wu Q, Wang W, Sun ZK, et al. Prevalence of unruptured cerebral aneurysms in Chinese adults aged 35 to 75 years: a cross-sectional study. Ann Intern Med. 2013;159:514–21.
Rinkel GJ. Intracranial aneurysm screening: indications and advice for practice. Lancet Neurol. 2005;4:122–8.
Vlak MH, Algra A, Brandenburg R, Rinkel GJ. Prevalence of unruptured intracranial aneurysms, with emphasis on sex, age, comorbidity, country, and time period: a systematic review and meta-analysis. Lancet Neurol. 2011;10:626–36.
Rozen S, Skaletsky H. Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000;132:365–86.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−∆∆C(T)) method. Methods. 2001;25:402–8.
Greving JP, Wermer MJ, Brown RD Jr, Morita A, Juvela S, Yonekura M, Ishibashi T, Torner JC, Nakayama T, Rinkel GJ, Algra A. Development of the PHASES score for prediction of risk of rupture of intracranial aneurysms: a pooled analysis of six prospective cohort studies. Lancet Neurol. 2014;13:59–66.
Ishibashi T, Murayama Y, Urashima M, Saguchi T, Ebara M, Arakawa H, Irie K, Takao H, Abe T. Unruptured intracranial aneurysms: incidence of rupture and risk factors. Stroke. 2009;40:313–6.
Morita A, Kirino T, Hashi K, Aoki N, Fukuhara S, Hashimoto N, Nakayama T, Sakai M, Teramoto A, Tominari S, Yoshimoto T. The natural course of unruptured cerebral aneurysms in a Japanese cohort. N Engl J Med. 2012;366:2474–82.
Wermer MJ, van der Schaaf IC, Velthuis BK, Majoie CB, Albrecht KW, Rinkel GJ. Yield of short-term follow-up CT/MR angiography for small aneurysms detected at screening. Stroke. 2006;37:414–8.
Wiebers DO, Whisnant JP, Huston J 3rd, Meissner I, Brown RD Jr, Piepgras DG, Forbes GS, Thielen K, Nichols D, O’Fallon WM, et al. Unruptured intracranial aneurysms: natural history, clinical outcome, and risks of surgical and endovascular treatment. Lancet. 2003;362:103–10.
Yonekura M. Small unruptured aneurysm verification (SUAVe study, Japan)—interim report. Neurol Med Chir (Tokyo). 2004;44:213–4.
Frosen J, Tulamo R, Paetau A, Laaksamo E, Korja M, Laakso A, Niemela M, Hernesniemi J. Saccular intracranial aneurysm: pathology and mechanisms. Acta Neuropathol. 2012;123:773–86.
Deng X, Campagne F. Introduction to the development and validation of predictive biomarker models from high-throughput data sets. In: Bang H, Zhou XK, van Epps HL, Mazumdar M, editors. Statistical methods in molecular biology. Totowa: Humana Press; 2010. p. 435–70.
Jagga ZG, Gupta D. Machine learning for biomarker identification in cancer research—developments toward its clinical application. Pers Med. 2015;12:371–87.
Weir B. Unruptured intracranial aneurysms: a review. J Neurosurg. 2002;96:3–42.
Meng H, Tutino VM, Xiang J, Siddiqui A. High WSS or low WSS? Complex interactions of hemodynamics with intracranial aneurysm initiation, growth, and rupture: toward a unifying hypothesis. AJNR Am J Neuroradiol. 2014;35:1254–62.
Tulamo R, Frosen J, Hernesniemi J, Niemela M. Inflammatory changes in the aneurysm wall: a review. J Neurointerv Surg. 2010;2:120–30.
Chalouhi N, Points L, Pierce GL, Ballas Z, Jabbour P, Hasan D. Localized increase of chemokines in the lumen of human cerebral aneurysms. Stroke. 2013;44:2594–7.
Serra R, Volpentesta G, Gallelli L, Grande R, Buffone G, Lavano A, de Franciscis S. Metalloproteinase-9 and neutrophil gelatinase-associated lipocalin plasma and tissue levels evaluation in middle cerebral artery aneurysms. Br J Neurosurg. 2014. https://doi.org/10.3109/02688697.2014.913777.
Pera J, Korostynski M, Krzyszkowski T, Czopek J, Slowik A, Dziedzic T, Piechota M, Stachura K, Moskala M, Przewlocki R, Szczudlik A. Gene expression profiles in human ruptured and unruptured intracranial aneurysms: what is the role of inflammation? Stroke. 2010;41:224–31.
Tromp G, Weinsheimer S, Ronkainen A, Kuivaniemi H. Molecular basis and genetic predisposition to intracranial aneurysm. Ann Med. 2014;46:597–606.
Strong MJ, Amenta PS, Dumont AS, Medel R. The role of leukocytes in the formation and rupture of intracranial aneurysms. Neuroimmunol Neuroinflammation. 2015;2:8.
Marbacher S, Marjamaa J, Bradacova K, von Gunten M, Honkanen P, Abo-Ramadan U, Hernesniemi J, Niemela M, Frosen J. Loss of mural cells leads to wall degeneration, aneurysm growth, and eventual rupture in a rat aneurysm model. Stroke. 2014;45:248–54.
Strong M, Amenta P, Dumont A, Medel R. The role of leukocytes in the formation and rupture of intracranial aneurysms. Neuro Immunol Inflamm. 2015;2:107–14.
Gounis MJ, Vedantham S, Weaver JP, Puri AS, Brooks CS, Wakhloo AK, Bogdanov AA Jr. Myeloperoxidase in human intracranial aneurysms: preliminary evidence. Stroke. 2014;45:1474–7.
Chu Y, Wilson K, Gu H, Wegman-Points L, Dooley SA, Pierce GL, Cheng G, Pena Silva RA, Heistad DD, Hasan D. Myeloperoxidase is increased in human cerebral aneurysms and increases formation and rupture of cerebral aneurysms in mice. Stroke. 2015;46:1651–6.
Leopold JA. The central role of neutrophil gelatinase-associated lipocalin in cardiovascular fibrosis. Hypertension. 2015;66:20–2.
Stapleton PP, Redmond HP, Bouchier-Hayes DJ. Myeloperoxidase (MPO) may mediate neutrophil adherence to the endothelium through upregulation of CD11B expression—an effect downregulated by taurine. Adv Exp Med Biol. 1998;442:183–92.
Schett G, Dayer JM, Manger B. Interleukin-1 function and role in rheumatic disease. Nat Rev Rheumatol. 2016;12:14–24.
Paulsson JM, Moshfegh A, Dadfar E, Held C, Jacobson SH, Lundahl J. In-vivo extravasation induces the expression of interleukin 1 receptor type 1 in human neutrophils. Clin Exp Immunol. 2012;168:105–12.
Peters VA, Joesting JJ, Freund GG. IL-1 receptor 2 (IL-1R2) and its role in immune regulation. Brain Behav Immun. 2013;32:1–8.
Bai M, Grieshaber-Bouyer R. CD177 modulates human neutrophil migration through activation-mediated integrin and chemoreceptor regulation. Blood. 2017;130:2092–100.
Borregaard N. Neutrophils, from marrow to microbes. Immunity. 2010;33:657–70.
Leung BP, Culshaw S, Gracie JA, Hunter D, Canetti CA, Campbell C, Cunha F, Liew FY, McInnes IB. A role for IL-18 in neutrophil activation. J Immunol. 2001;167:2879–86.
Zen K, Liu Y, McCall IC, Wu T, Lee W, Babbin BA, Nusrat A, Parkos CA. Neutrophil migration across tight junctions is mediated by adhesive interactions between epithelial coxsackie and adenovirus receptor and a junctional adhesion molecule-like protein on neutrophils. Mol Biol Cell. 2005;16:2694–703.
Weinsheimer S, Lenk GM, van der Voet M, Land S, Ronkainen A, Alafuzoff I, Kuivaniemi H, Tromp G. Integration of expression profiles and genetic mapping data to identify candidate genes in intracranial aneurysm. Physiol Genom. 2007;32:45–57.
Miguel LI, Almeida CB, Traina F, Canalli AA, Dominical VM, Saad ST, Costa FF, Conran N. Inhibition of phosphodiesterase 9A reduces cytokine-stimulated in vitro adhesion of neutrophils from sickle cell anemia individuals. Inflamm Res. 2011;60:633–42.
Li Z, Huang G, Chen L, Tan H, Wang Z, Liu L, Shi Y, Yin C, Wang Q. Construction of transcriptional regulatory network in the intracranial aneurysms. Int J Clin Exp Med. 2017;10:1583–91.
Coulombe F, Jaworska J, Verway M, Tzelepis F, Massoud A, Gillard J, Wong G, Kobinger G, Xing Z, Couture C, et al. Targeted prostaglandin E2 inhibition enhances antiviral immunity through induction of type I interferon and apoptosis in macrophages. Immunity. 2014;40:554–68.
Ollikainen E, Tulamo R, Kaitainen S, Honkanen P, Lehti S, Liimatainen T, Hernesniemi J, Niemela M, Kovanen PT, Frosen J. Macrophage infiltration in the saccular intracranial aneurysm wall as a response to locally lysed erythrocytes that promote degeneration. J Neuropathol Exp Neurol. 2018;77:890–903.
Ollikainen E, Tulamo R, Lehti S, Lee-Rueckert M, Hernesniemi J, Niemela M, Yla-Herttuala S, Kovanen PT, Frosen J. Smooth muscle cell foam cell formation, apolipoproteins, and ABCA1 in intracranial aneurysms: implications for lipid accumulation as a promoter of aneurysm wall rupture. J Neuropathol Exp Neurol. 2016;75:689–99.
Groselj-Grenc M, Ihan A, Derganc M. Neutrophil and monocyte CD64 and CD163 expression in critically ill neonates and children with sepsis: comparison of fluorescence intensities and calculated indexes. Mediators Inflamm. 2008;2008:202646.
Liu Y, Tan S, Huang L, Abramovitch RB, Rohde KH, Zimmerman MD, Chen C, Dartois V, VanderVen BC, Russell DG. Immune activation of the host cell induces drug tolerance in Mycobacterium tuberculosis both in vitro and in vivo. J Exp Med. 2016;213:809–25.
Chinetti-Gbaguidi G, Copin C, Derudas B, Vanhoutte J, Zawadzki C, Jude B, Haulon S, Pattou F, Marx N, Staels B. The coronary artery disease-associated gene C6ORF105 is expressed in human macrophages under the transcriptional control of PPARgamma. FEBS Lett. 2015;589:461–6.
Alexander JS, Dayton T, Davis C, Hill S, Jackson TH, Blaschuk O, Symonds M, Okayama N, Kevil CG, Laroux FS, et al. Activated T-lymphocytes express occludin, a component of tight junctions. Inflammation. 1998;22:573–82.
Kangelaris KN, Prakash A, Liu KD, Aouizerat B, Woodruff PG, Erle DJ, Rogers A, Seeley EJ, Chu J, Liu T, et al. Increased expression of neutrophil-related genes in patients with early sepsis-induced ARDS. Am J Physiol Lung Cell Mol Physiol. 2015;308:L1102–13.
Kossenkov AV, Dawany N, Evans TL, Kucharczuk JC, Albelda SM, Showe LC, Showe MK, Vachani A. Peripheral immune cell gene expression predicts survival of patients with non-small cell lung cancer. PLoS ONE. 2012;7:e34392.
Ruiz Esparza-Garrido R, Rodriguez-Corona JM, Lopez-Aguilar JE, Rodriguez-Florido MA, Velazquez-Wong AC, Viedma-Rodriguez R, Salamanca-Gomez F, Velazquez-Flores MA. Differentially expressed long non-coding RNAs were predicted to be involved in the control of signaling pathways in pediatric astrocytoma. Mol Neurobiol. 2017;54:6598–608.
Damdinsuren B, Dement-Brown J, Li H, Tolnay M. B cell receptor induced Fc receptor-like 5 expression is mediated by multiple signaling pathways converging on NF-kappaB and NFAT. Mol Immunol. 2016;73:112–21.
So MY, Tian Z, Phoon YS, Sha S, Antoniou MN, Zhang J, Wu RS, Tan-Un KC. Gene expression profile and toxic effects in human bronchial epithelial cells exposed to zearalenone. PLoS ONE. 2014;9:e96404.
Chyatte D, Bruno G, Desai S, Todor DR. Inflammation and intracranial aneurysms. Neurosurgery. 1999;45:1137–46 (discussion 1146-1137).
Frosen J, Piippo A, Paetau A, Kangasniemi M, Niemela M, Hernesniemi J, Jaaskelainen J. Remodeling of saccular cerebral artery aneurysm wall is associated with rupture: histological analysis of 24 unruptured and 42 ruptured cases. Stroke. 2004;35:2287–93.
Kataoka K, Taneda M, Asai T, Kinoshita A, Ito M, Kuroda R. Structural fragility and inflammatory response of ruptured cerebral aneurysms. A comparative study between ruptured and unruptured cerebral aneurysms. Stroke. 1999;30:1396–401.
Study conception and design: VMT, JNJ, JK, HM. Acquisition of data: VMT, KEP, KJ, HS, AHS, KVS, EIL. Analysis and interpretation of data: VMT, KEP, LL, YS, KJ, JNJ, YS, JK, HM. Drafting of manuscript: VMT, KEP, LL, JNJ, YS, JK, HM. Critical revision: VMT, KEP, LL, HS, KJ, JNJ, YS, KVS, EIL, AHS, KJ, HM. All authors read and approved the final manuscript.
We thank the patients who participated in this study; Ashish Sonig MD MS MCh, Sabareesh Natarajan MD MS, Ning Lin MD, Chandan Krishna MD, Marshall C. Cress MD, Leonardo Rangel-Castilla MD, and Maxim Mokin MD PhD for blood collection; Jonathan Bard MA and Brandon Marzulo MS for sequencing assistance; Jennifer L. Gay CCRP for study protocol management; Paul H. Dressel BFA for preparation of the illustrations; and Debra J Zimmer for editorial assistance.
KEP, LL, HS, KJ, YS, JK—None. VMT—Co-founder: Neurovascular Diagnostics, Inc. JNJ—Principal Investigator: NIH Grant R01-AR-060604. AHS—Financial Interest/Investor/Stock Options/Ownership: Amnis Therapeutics, Apama Medical, BlinkTBI, Inc, Buffalo Technology Partners, Inc., Cardinal Health, Cerebrotech Medical Systems, Inc, Claret Medical, Cognition Medical, Endostream Medical, Ltd, Imperative Care, International Medical Distribution Partners, Rebound Therapeutics Corp., Silk Road Medical, StimMed, Synchron, Three Rivers Medical, Inc., Viseon Spine, Inc. Consultant/Advisory Board: Amnis Therapeutics, Boston Scientific, Canon Medical Systems USA, Inc., Cerebrotech Medical Systems, Inc., Cerenovus, Claret Medical, Corindus, Inc., Endostream Medical, Ltd, Guidepoint Global Consulting, Imperative Care, Integra, Medtronic, MicroVention, Northwest University—DSMB Chair for HEAT Trial, Penumbra, Rapid Medical, Rebound Therapeutics Corp., Silk Road Medical, StimMed, Stryker, Three Rivers Medical, Inc., VasSol, W.L. Gore & Associates. National PI/Steering Committees: Cerenovus LARGE Trial and ARISE II Trial, Medtronic SWIFT PRIME and SWIFT DIRECT Trials, MicroVention FRED Trial & CONFIDENCE Study, MUSC POSITIVE Trial, Penumbra 3D Separator Trial, COMPASS Trial, INVEST Trial.
KVS—Consulting and teaching for Canon Medical Systems Corporation, Penumbra Inc., Medtronic, and Jacobs Institute. Co-Founder: Neurovascular Diagnostics, Inc.
EIL—Intratech Medical Ltd., NeXtGen Biologics. Principal investigator: Medtronic US SWIFT PRIME Trials. Honoraria–Medtronic. Consultant–Pulsar Vascular. Advisory Board-Stryker, NeXtGen Biologics, MEDX. Cognition Medical. Other financial support—Abbott Vascular for carotid training sessions.
HM—Principal investigator: the 3 Brain Aneurysm Foundation grants mentioned below (no applicable grant numbers), NIH Grants R01-NS-091075 and R01-NS-064592. Grant support: Toshiba Medical Systems. Co-founder: Neurovascular Diagnostics, Inc.
Availability of data and materials
The datasets used in the current study are available from the corresponding author upon reasonable request.
Consent for publication
No specific patient information or images are given in the manuscript.
Ethics approval and consent to participate
This study was approved by the University at Buffalo Health Sciences Institutional Review Board (Study No. 030-474433). All methods were carried out in accordance with the approved protocol and informed consent was obtained from all subjects.
This work was supported by the Brain Aneurysm Foundation Carol Harvey Chair of Research Grant (HM), the Brain Aneurysm Foundation Team Cindy-Alcatraz Chair of Research Grant (HM), the Brain Aneurysm Foundation Kristen’s Legacy of Love Chair of Research Grant (HM), and NIH Grant # R01-AR-060604 (JNJ).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.