|Home | About | Journals | Submit | Contact Us | Français|
Brain arteriovenous malformations (BAVMs) are an important cause of intracranial hemorrhage (ICH) in young adults. Gene expression profiling of blood has led to the identification of stroke biomarkers, and may help identify BAVM biomarkers and illuminate BAVM pathogenesis. It is unknown whether blood gene expression profiles differ between 1) BAVM patients and healthy controls, or 2) unruptured and ruptured BAVM patients at presentation. We characterized blood transcriptional profiles in 60 subjects (20 unruptured BAVM, 20 ruptured BAVM, and 20 healthy controls) using Affymetrix whole genome expression arrays. Expression differences between groups were tested by ANOVA, adjusting for potential confounders. Genes with absolute fold change ≥ 1.2 (false discovery rate corrected p ≤ 0.1) were selected as differentially expressed and evaluated for over-representation in KEGG biological pathways (p ≤ 0.05). Twenty-nine genes were differentially expressed between unruptured BAVM patients and controls, including 13 which may be predictive of BAVM. Patients with ruptured BAVM compared to unruptured BAVM differed in expression of 1490 genes, with over-representation of genes in 8 pathways including MAPK, VEGF, Wnt signaling and several inflammatory pathways. These results suggest clues to the pathogenesis of BAVM and/or BAVM rupture and point to potential biomarkers or new treatment targets.
Brain arteriovenous malformations (BAVMs) are a tangle of abnormal vessels directly shunting blood from the arterial to venous circulation, resulting in high-flow lesions prone to rupture . BAVMs account for ~2% of all hemorrhagic strokes each year, and are an important cause of intracerebral hemorrhage (ICH) in young adults. At present, only a few robust risk predictors exist, including deep venous drainage, deep location, and ICH presentation . Approximately 50% of patients present initially with ICH, and ICH presentation is the strongest predictor of future hemorrhage [2–4]. However, these risk factors do not help those patients who present unruptured with other clinical symptoms (e.g., seizure, headache, focal neurologic deficit) or who are otherwise asymptomatic. Thus, identification of biomarkers that predict BAVM and/or ICH risk would help to stratify risk in BAVM patients for optimal patient management and represent a major advance in clinical care.
Gene (mRNA) expression profiles of tissue and peripheral blood can reflect pathological changes in the body and brain. Numerous studies have demonstrated that even closely related vascular diseases (e.g., stroke subtypes , intracranial aneurysm  and abdominal aortic aneurysm [7, 8]), have very different gene signatures, and even individuals with the same disease but different severity (e.g., coronary artery disease) may be distinguished [9, 10]. For many diseases, especially those which affect the vasculature, tissue specimens are usually only obtained after a patient has become symptomatic and treated. Peripheral blood continuously interacts with the tissues of the body and different disease states may trigger specific changes in blood cell gene expression which could be measured at earlier stages of disease . For example, a signature set of 41 genes identified by microarray analysis of patient blood was shown to identify asymptomatic thoracic aortic aneurysm patients with high accuracy (~80%) . Therefore, it seems likely that patients harboring BAVM (or those at increased risk for future bleeding) will display a distinct gene expression profile in their blood, which could be used for better risk stratification and patient selection for invasive treatment, especially for unruptured BAVM patients. Prior data from our group suggests that inflammation plays an important role in the pathogenesis of BAVM [12–14]. Thus, the rationale for using peripheral whole blood to measure gene expression in BAVM patients is that inflammatory cells (present in whole blood) may act as reporters, contributing to the distinct gene expression profiles that reflect the presence or severity of disease states associated with BAVM.
The purpose of this study was to perform a preliminary gene expression analysis of whole blood samples from unruptured BAVM patients and healthy controls to characterize the expression profile for BAVM. Additionally, we examined gene expression profiles in the blood of unruptured BAVM patients compared to ruptured BAVM patients to identify genes that are regulated in ICH. In this discovery-driven study, we demonstrate initial evidence that blood gene expression profiling may be a useful tool to identify genes that are regulated in BAVM disease or ICH presentation.
Forty BAVM patients (20 ruptured and 20 unruptured at presentation) were recruited at the University of California, San Francisco (UCSF), and diagnosed using standardized guidelines . Ruptured BAVM at initial presentation was defined as evidence of new intracranial blood on computed tomography or magnetic resonance imaging regardless of clinical presentation. Unruptured BAVM diagnosed by angiography was defined as all other presentations without evidence of new bleeding, including seizure, focal ischemic deficit, headache, apparently unrelated symptoms or asymptomatic, incidental discovery. Twenty controls were recruited from healthy ambulatory volunteers at Stanford University, with no history of cardiovascular or cerebrovascular disease. All subjects provided written, informed consent and blood specimens for genetic studies, in full accordance with the declaration of Helsinki principles. The study was approved by the Institutional Review Boards at UCSF, University of California at Davis, and Stanford University.
After informed consent, whole blood (15 mL) was drawn via venipuncture from each subject into PAXgene tubes (PreAnalytiX, Hilden, Germany) and processed as previously described . Total RNA preparation included RNA from all cells in whole blood: polymorphonuclear cells (neutrophils, basophils, and eosinophils), mononuclear cells (lymphocytes and macrophages/monocytes), platelet precursors and red blood cell precursors.
Preparation of biotin-labeled cDNA from 50 ng total RNA and hybridization to Affymetrix GeneChip® Human Genome U133 Plus 2.0 arrays was performed at the University of California at Davis Genomics facility as previously described . Microarrays were hybridized and scanned in batches over the course of 10 weeks, with approximately equal representation of control, unruptured BAVM, and ruptured BAVM samples processed each week. For data analysis, microarrays were grouped by batch week (n=10).
Probe-level data were preprocessed using the robust multi-array analysis (RMA) normalization method and analyzed using Partek Genomics Suite (Partek Inc., St. Louis, MO) and GeneSpring GX (version 10.0) (Agilent Technologies, Inc., Santa Clara, CA). Expression values in each sample were baseline-corrected relative to the median expression value from all samples.
Our primary analysis compared expression profiles from unruptured BAVM with that of healthy controls to identify genes that are differentially expressed in BAVM. We included only unruptured BAVM in this analysis because confounders such as the hemorrhagic event itself and treatment effects may hinder the identification of genes specifically related to BAVM. Secondary analysis compared expression profiles between unruptured and ruptured BAVM to identify genes related to disease severity. Three types of statistical analyses were used to facilitate the interpretation of gene expression data: 1) analysis of variance (ANOVA) adjusting for covariates to identify differentially expressed genes between groups, 2) hierarchical clustering to observe gene expression patterns across subjects, and 3) prediction analysis to determine if a specific set of genes can accurately distinguish BAVM cases from controls. In addition, to provide a functional interpretation of rupture-related genes, we performed pathway analysis.
We performed ANOVA with equal variances for each probeset with adjustment for group and microarray batch (reduced model). We also evaluated a full multivariable model including additional adjustments for potential confounders including age at blood collection, gender, and race/ethnicity (white/non-white). In a secondary analysis, we compared unruptured to ruptured BAVM. Gene expression may also be influenced by timing of sample collection or treatment effects, thus we performed sensitivity analyses to explore the effects of these potential confounders on gene expression in BAVM patients, defining these continuous variables as log10 days from BAVM presentation or treatment to blood draw.
To adjust for multiple testing, the Benjamini-Hochberg false discovery rate (FDR) was used . Given our relatively small sample size (n=20 per group), we chose to use a FDR corrected significance threshold p ≤ 0.1, and differentially expressed genes were selected as those with an absolute fold change (|FC|) ≥ 1.2. This liberal cutoff, which is commonly used in hypothesis-generating, discovery-driven gene expression studies [17, 18], was chosen to avoid declaring too many false negatives, while accepting that 10% of the genes identified as significant may be false positives. Descriptive statistics including Student’s t-test or χ2 test were performed using Intercooled Stata 10 statistical software (StataCorp, TX).
To visualize the expression profiles that distinguish 1) unruptured BAVM patients from healthy controls or 2) unruptured BAVM patients from ruptured BAVM patients, the differentially expressed genes identified from the two comparisons were organized using the hierarchical clustering method (Pearson’s centered distance metric and average linkage rule). This method groups co-expressed genes into clusters based on similarity of expression values among genes and sample groups. Clustering was performed for genes identified as differentially expressed in the unruptured BAVM patients vs. healthy control analysis (FDR ≤ 0.1, |FC| ≥ 1.2). For the unruptured vs. ruptured BAVM analysis, the number of differentially expressed genes at the chosen significance threshold (FDR ≤ 0.1, |FC| ≥ 1.2) was very large; hence we chose to perform the cluster analysis using the subset of genes with the largest fold change values from the full multivariate ANOVA model (FDR ≤ 0.1, |FC| ≥ 2.0).
To identify a minimum set of genes (represented among the 54,675 probesets tested) that can best distinguish unruptured BAVM patients from healthy controls, we performed prediction analysis using Prediction Analysis of Microarrays software which uses the ‘nearest shrunken centroids’ method . The accuracy of the classifier (e.g., how well the classifier correctly predicts BAVM disease [sensitivity] and control status [specificity]) is then assessed using 10-fold cross-validation .
To provide a functional interpretation of the genes related to BAVM rupture, we performed pathway analysis using the web-based gene set analysis toolkit (WebGestalt) . Differentially expressed genes identified in the full multivariate ANOVA model (FDR ≤ 0.1 and |FC| ≥ 1.2) were evaluated for over representation in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (Fisher’s exact p ≤ 0.05). The reference probesets (n=17,730) included those with median expression values ≥1 in at least 1 sample. Significant KEGG pathways with >5 genes observed in the differentially expressed gene list were reported. Since genes that play a functional role in disease pathogenesis may have small expression differences that are difficult to detect in a small sample with confounding factors, we also present pathway results for the larger set of differentially expressed genes identified in the reduced ANOVA model (FDR ≤ 0.1 and |FC| ≥ 1.2).
Demographic and clinical characteristics for BAVM patients and controls are summarized in Table 1. All BAVM patients had blood specimens collected after initial presentation; most (n=18, 90%) ruptured BAVM patients were treated either by embolization, AVM resection, radiosurgery, or some combination prior to blood collection. There were no significant differences between controls and unruptured BAVM patients regarding gender (p = 1.000) or race/ethnicity (p = 0.749). However, controls were significantly younger than BAVM patients at time of blood collection (p = 0.001). Unruptured BAVM patients initially presented with headache (n=5, 25%), seizure (n=9, 45%), focal neurologic deficit (n=1, 5%), or incidentally (n=5, 25%). There were no significant differences between unruptured and ruptured BAVM patient groups regarding age at diagnosis, gender or race/ethnicity (Table 1). Ruptured patients were significantly younger at time of blood collection than unruptured patients (p = 0.037), presented with smaller BAVM lesion size (p = 0.032), and had a greater percentage with deep venous drainage (p = 0.001). Spetzler-Martin grades (surgical risk score) were similar in ruptured (range 1 – 4, mean 2.6, SD 0.8) and unruptured BAVM patients (range 1 – 5, mean 2.7, SD 1.3).
The primary analysis compared blood gene expression profiles between 20 unruptured BAVM patients and 20 healthy controls to identify genes showing consistent expression differences. A total of 33 mRNA probesets (representing 29 known genes) were identified as differentially expressed in unruptured BAVM vs. control blood samples after adjustment for microarray batch (FDR ≤ 0.1, |FC| ≥ 1.2) (Supplementary Table 1). The observed absolute fold change in blood expression for these genes ranged from 1.4 to 2.1, with a mean of 1.4. There were 31 mRNA probesets (28 known genes) with decreased expression in unruptured BAVM including the most differentially expressed gene IL1RAP (|FC| = 2.1). Only 2 mRNA probesets had increased expression in the unruptured BAVM group: one of which represented a known gene, GPR174 (|FC| = 1.6)
Hierarchical cluster analysis of the 33 differentially expressed mRNA probesets showed segregation of unruptured BAVM patients and controls into two main branches, with some exceptions (Figure 1). The dendrogram shows good segregation of unruptured BAVM cases and controls related to genes on the far left and far right of the cluster, while it is difficult to distinguish 9 individuals in the middle of the cluster. The two BAVM patients that cluster with the controls (left side of dendrogram) had blood specimens collected >5 years after initial presentation, while 85% of the unruptured BAVM patients had blood collection <5 years after initial presentation. Hence, this longer time between BAVM presentation and blood collection may partially explain why gene expression profiles for these 2 patients are more similar to controls and suggests their lesions may be less biologically active (ie. more stable). Unruptured patients tended to have a greater proportion of genes down-regulated (Figure 1, blue color), whereas controls had a greater proportion of the same genes up-regulated (Figure 1, red color).
We next examined the effect of possible confounding factors (e.g., age, gender and race/ethnicity) on gene expression. Among the unruptured BAVM cases and controls, we identified 3 age-related, 60 gender-related, and 0 race/ethnicity-related genes (FDR ≤ 0.1, |FC| ≥ 1.2). None of the age or gender-related genes overlapped with the 29 unruptured BAVM-related genes. However, when age, gender and race/ethnicity were included in the full multivariate ANOVA model, no genes were significantly differentially expressed in unruptured BAVM patients compared to controls (FDR ≤ 0.1).
We performed prediction analysis to determine the minimum number of genes needed to accurately distinguish unruptured BAVM patients from controls. A minimum of 47 out of 54,675 mRNA probesets, representing 33 distinct known genes, best differentiated unruptured BAVM patients from controls (Table 2) with a sensitivity of 70% (correctly predict unruptured BAVM) and specificity of 80% (correctly predict controls), and a cross-validated overall misclassification rate of 25%. Prediction model results were consistent with the reduced ANOVA model, as 13 out of 47 BAVM predictor probesets were differentially expressed among unruptured BAVM patients vs. controls (FDR ≤ 0.1, |FC| ≥ 1.2).
In a secondary analysis, we examined the blood RNA profiles for 20 unruptured and 20 ruptured BAVM patients to identify genes reflecting disease severity. In the reduced ANOVA model, 2809 differentially expressed mRNA probesets (FDR ≤ 0.1, |FC| ≥ 1.2), representing 1490 known genes (Supplementary Table 2), were identified. The observed fold change in blood expression for these genes ranged from 1.2 to 4.7, with a mean of 1.4. After further adjustments for age, gender, and race/ethnicity, 896 differentially expressed mRNA probesets remained associated (FDR ≤ 0.1, |FC| ≥ 1.2), representing 493 known genes. Of these, 49 mRNA probesets had fold change values ≥ 2.0 (Supplementary Table 3). Most of these mRNA probesets (42/49) had higher expression in ruptured BAVM patients including probesets for ARG1, TNFAIP6, IL1R1, IL1RAP, IL18R1, and GPR97, while 7 mRNA probesets had decreased expression including probesets for NOG and LEF1. Only three genes (IL1RAP, RGL4, and IDS) were differentially expressed in both the primary (unruptured BAVM vs. control) and secondary (unruptured BAVM vs. ruptured BAVM) analyses.
Hierarchical cluster analysis of the 49 mRNA probesets that had the largest differences between ruptured and unruptured BAVM patients (full multivariate model, FDR ≤ 0.1, |FC| ≥ 2.0) revealed clusters of co-expressed genes as shown in the gene expression heat map (Figure 2). The dendrogram shows good segregation of 11 ruptured (far left cluster) and 13 unruptured (far right cluster) patients. Ruptured patients tended to have a greater proportion of genes up-regulated (Figure 2, red color), whereas unruptured patients had a greater proportion of the same genes down-regulated (Figure 2, blue color).
Timing of blood collection and treatment itself can be potential confounders of blood gene expression levels, especially with respect to our secondary analysis comparing unruptured to ruptured BAVM patients. To assess whether the gene expression profiles attributed to BAVM rupture were confounded by blood collection or treatment time, we looked at the overlap between genes differentially expressed between ruptured and unruptured patients, in samples with different blood collection times, and patients with different treatment times since rupture. There were 32 mRNA probesets with significant changes in expression associated with blood collection time (FDR ≤ 0.1), including only 1 (representing MAP2K6) that was also differentially expressed between ruptured and unruptured BAVM patients (FDR ≤ 0.1, |FC| ≥ 1.2). These results suggest that the expression levels of a subset of genes may be influenced by the blood collection time after clinical presentation. However, these genes are different from those identified as BAVM rupture-related genes. There were no significant differences in mRNA expression associated with treatment time (one-way ANOVA FDR ≤ 0.1) for the 18 ruptured BAVM patients treated prior to blood collection.
To identify potentially novel functional pathways that influence BAVM rupture, we performed exploratory pathway analysis for the 493 genes that were differentially expressed between unruptured and ruptured BAVM patients in the full multivariate ANOVA model. BAVM rupture-related genes were over represented in the following KEGG biological pathways: natural killer cell mediated cytotoxicity (p = 0.01), T cell receptor signaling (p = 0.01), antigen processing and presentation (p = 0.02) and Wnt signaling (p = 0.05) (Table 3). In addition, the following KEGG metabolic and genetic information processing pathways were significant: Aminoacyl-tRNA biosynthesis (p < 0.001), N-Glycan biosynthesis (p < 0.001), and Arginine and proline metabolism (p = 0.03) (Table 3).
When we expanded the rupture-related gene set to include the 1490 genes that were differentially expressed between ruptured and unruptured BAVM patients in the reduced ANOVA model, we observed 8 KEGG biological pathways with gene over representation (p ≤ 0.05), summarized in Table 3. Three of these pathways were consistent with the multivariate ANOVA results including natural killer cell mediated cytotoxicity (p < 0.001), T cell receptor signaling (p = 0.004), and Wnt signaling (p = 0.05). Five additional pathways were identified only in the expanded gene set including MAPK signaling (p = 0.004), VEGF signaling (p = 0.02), GnRH signaling (p = 0.04), and the inflammatory pathways Fc epsilon RI signaling (p = 0.01) and B cell receptor signaling (p = 0.04). Six of the BAVM rupture-related genes are also in the Notch signaling pathway; however, this finding was not statistically significant (p = 0.06). BAVM rupture-related genes in the expanded gene set were also over represented in 10 metabolic or genetic information processing pathways, 2 which were consistent with the multivariate ANOVA results, and two additional human disease pathways described in Table 3.
Our study is the first to examine blood mRNA expression patterns in BAVM patients to identify genes regulated in BAVM disease. We identified 33 probesets differentially expressed between unruptured BAVM and controls, including 13 which may be predictive of unruptured BAVM. However, the misclassification rate (25%) using these 33 probesets was too high to correctly predict presence of unruptured BAVM for clinical purposes. The expression profiles of ruptured BAVM patients at presentation were markedly different from unruptured BAVM, with many more genes (n=493) showing differential expression. Functional interpretation of these data via pathway analysis suggests novel biological pathways may contribute to the pathogenesis of BAVM rupture.
Three differentially expressed genes (IL1RAP, RGL4, and IDS) overlapped between the primary (unruptured BAVM vs. control) and secondary (unruptured BAVM vs. ruptured BAVM) analysis and are all novel candidate genes for BAVM. We observed down-regulation of these genes in unruptured BAVM patients relative to controls while ruptured BAVM patients had an up-regulation of these three genes relative to unruptured BAVM. This may reflect temporal changes in expression related to disease progression that we are unable to address with the present data. Interleukin-1 receptor accessory protein (IL1RAP) encodes a protein that forms a complex with interleukin-1 to induce synthesis of acute phase and pro-inflammatory proteins during infection, tissue damage, or stress. Interestingly, we have previously reported association of polymorphisms in a related cytokine, interleukin-1-beta (IL1B), with susceptibility of both BAVM and ICH presentation .
In this study we observed a down-regulation of protein C (PROC) in the blood of unruptured BAVM patients compared to controls. PROC is involved in thrombin generation and has been identified as a potential blood biomarker for the presence of aortic aneurysm . In addition, mutations in the PROC gene are known to be associated with venous thrombosis in individuals with inherited protein C deficiency . While the role of the venous system in BAVM pathogenesis has not been defined, studies have suggested that venous thrombosis may be involved [25, 26].
In our study, several inflammatory and angiogenesis-related genes were differentially expressed in unruptured BAVM patient blood compared to controls including ANGPT2, ITGAV, TIE1, and TEK (TIE2), though these findings were not statistically significant at FDR ≤ 0.1. These results support previous genome-wide expression studies conducted in BAVM tissue that identified up-regulation of VEGFA, ENG, ANGPT2, ITGAV, FLT1 (VEGFR1), and MMP9, and down-regulation of TIE1, TEK (TIE2), and ANGPT1 [27, 28].
Many of the BAVM rupture-related genes identified in this study also function in inflammation. These findings are consistent with recent evidence that suggests a role for inflammatory cell involvement in BAVM pathogenesis, as both neutrophils and macrophages/microglia are present in the vascular wall and intervening stroma of unruptured, non-embolized AVMs . ICH leads to white blood cell infiltration, inflammatory cytokine production, and secondary inflammatory gene activation during apoptosis [29, 30]. Studies of human white blood cell gene expression patterns have also indicated an inflammatory reaction in blood after ischemic stroke . In this study, we identified an over representation of rupture-related genes in several inflammatory-related pathways that have not been previously implicated in BAVM pathogenesis including: natural killer cell mediated cytotoxicity, antigen processing and presentation, T and B cell receptor signaling, and the Fc epsilon RI signaling pathways. Examples of genes that function in inflammation that were over expressed in ruptured BAVM patients included: Fas (TNF receptor superfamily, member 6) (FAS), toll-like receptor 10 (TLR10), and tumor necrosis factor, alpha-induced protein 6 (TNFAIP6). It is also noteworthy that ruptured BAVM patients over expressed interleukin-1 receptor, type 1 (IL1R1), an important mediator involved in many cytokine-induced immune and inflammatory responses and reportedly up-regulated in human perihematomal tissue . We previously reported increased expression of interleukin-6 (IL6) in ruptured BAVM tissue. Consistent with these studies, we observed a higher level of IL6 mRNA expression in the blood of ruptured BAVM patients compared to unruptured BAVM patients; however, this finding was not statistically significant (|FC|=1.10, p=0.52). Blood levels of IL6 may differ from BAVM tissue because of temporal changes in IL6 expression related to BAVM progression, which may not have been detected in these cases at the time of sampling. In addition, lower fold change values in the blood may be difficult to measure by microarray.
This study also supports a role for VEGF, MAPK, Wnt, and Notch signaling pathways in BAVM or progression to hemorrhage. Prior evidence that these signaling pathways contribute to the pathogenesis of vascular malformations comes mainly from studies of animal and human tissues. For example, endothelial cells derived from BAVMs have aberrant angiogenic characteristics including up-regulation of VEGF pathway genes . The MAPK pathway has been implicated in cerebral cavernous malformation, a cerebrovascular disease with lesions characterized by dilated, thin-walled, leaky vessels . Wnt signaling influences vascular development, often in coordination with Notch signaling . Recent studies of Notch signaling in human BAVM and a BAVM mouse model suggest Notch activation may promote the development and even maintenance of BAVM [35, 36]. Therefore, our mRNA expression profiling data suggest that blood from BAVM patients may act as a reporter of dysregulated VEGF, MAPK, Wnt, and Notch pathway signaling in brain blood vessels. However, future studies will be needed to elucidate whether coordinated expression of genes in these pathways is relevant to blood vessel stability, progression to rupture, or other factors.
The relationship of tissue gene expression profiles to those observed in circulating blood cells is not well understood because of limitations in the availability and quality of matched tissues . Further, interpretation of expression data is complicated by tissue heterogeneity, in which expression levels reflect an average signal from numerous cell types. Despite these limitations, genome-wide gene expression studies of brain tissue and blood in vervet monkeys demonstrated brain transcriptional variation is reproduced in peripheral blood . Additionally, blood gene expression patterns mirror gene expression in atherosclerotic artery walls and correlate with the extent of coronary artery disease in humans . These studies support the use of blood as a surrogate tissue to identify peripheral (ie. systemic) gene expression changes that reflect, at least to some extent, gene expression at the local level.
Our study has several limitations which could affect the interpretation of results. First, we are unable to distinguish gene expression differences due to the hemorrhagic response from those due to factors that predispose to hemorrhage because gene expression measurements are based on blood collection at a single time point, after the hemorrhagic event or treatment has occurred. Timing of blood collection, or treatment type or other variables could lead to acute effects (local or systemic) that influence gene expression. It is possible that intracranial hemorrhage may exaggerate the expression response (up or down regulation) of genes that underlie BAVM rupture. However, our study suggests that the gene expression response is consistent with known gene function. For example, ruptured patients had increased expression of TNFAIP6, which can be induced by TNF-alpha and interleukin-1, and is involved in extracellular matrix stability and cell migration that is important in the protease network associated with inflammation. Neither TNF-alpha nor interleukin-1 were differentially expressed in ruptured BAVM patients, which may reflect temporal changes in expression or hemorrhage effects that we are unable to detect in the current study. However, several inflammatory genes, including IL1R1 and IL-6, also displayed similar direction of expression response (increased expression) in ruptured BAVM patients, which is consistent with previous studies of hemorrhage[12, 30]. We attempted to control for confounding effects by adjusting the analysis for age, gender, and race/ethnicity, and performed additional sensitivity analyses for timing of blood collection and treatment. However, a larger study would be needed to further adjust for the effect of additional factors, such as BAVM size. In addition, a subset of the 20 unruptured BAVM patients may actually be at high risk for bleeding and may experience a future hemorrhagic event. None of the unruptured BAVM patients in this study has suffered a hemorrhagic event to date.
We identified relatively few differentially expressed genes and observed low fold change values in blood gene expression, on the order of 1.5 to 2.0 fold, which may be more difficult to measure using microarrays . While the significance threshold varies substantially across blood expression profiling studies (ranging from no multiple testing correction to stringent FDR < 0.05), several studies with comparable sample size have reported similar numbers of genes and magnitudes of fold change values [5, 8, 38]. Isolation of specific cell types, such as neutrophils, may yield higher expression levels. Despite these limitations, we identified different expression profiles between unruptured BAVM cases and controls, and found that prediction modeling generally supported these findings. The differential expression profiles observed in this cohort, albeit at a liberal significance threshold (FDR ≤ 0.1), suggest new hypotheses for future BAVM studies. Technical validation of these results has not been performed, although Affymetrix microarray data has been shown to agree well with alternative methods . For rare diseases such as BAVM, the availability of biological replication cohorts is limited. Thus, future studies of larger BAVM cohorts with fewer confounding variables will be necessary to determine whether blood gene expression profiling will ultimately aid in the identification of BAVM predictors.
In conclusion, blood RNA profiling identified gene expression differences in patients who harbor BAVM. In addition, BAVM patients who present with ruptured BAVM also display gene expression profiles that differ from unruptured BAVM patients. These results may provide clues to the pathogenesis of BAVM and/or rupture and point to new pathways for investigation. Future studies coupling BAVM tissue (e.g., endothelial cell) gene expression data together with peripheral blood gene expression data from the same patients will provide a more complete picture of which genes are involved in BAVM.
We thank patients who participated in this study and members of the UCSF Brain AVM project for assistance with patient recruitment, technical support, and data management. This study was supported by the National Institutes of Health (NIH) Grants P01 NS044155 and R01 NS034949 (WLY), supplement 3R01NS034949-14S1 (WLY, FRS), K23 NS058357 (HK), NIH/NIGMS T32 GM08440 (SW), NINDS and ARRA NS056302 (FRS); the American Heart Association Bugher Foundation Center for Stroke Prevention Research Grant (FRS, WLY); and The Aneurysm and AVM Foundation Award (WLY).