Longitudinal Risk of Intracranial Hemorrhage in Patients With Arteriovenous Malformation of the Brain Within a Defined Population
Background and Purpose— Accurate estimates for risk and rates of intracranial hemorrhage (ICH) in the natural course of patients harboring brain arteriovenous malformation (BAVM) are needed to provide a quantitative basis for planning clinical trials to evaluate interventional strategies and to help guide practice management.
Methods— We identified patients with BAVM at the Kaiser Permanente Northern California health maintenance organization and documented their clinical course. The influences of age at diagnosis, gender, race-ethnicity, ICH at presentation, venous draining pattern, and BAVM size on ICH subsequent to presentation were studied using the multivariate Cox proportional hazards model and Kaplan-Meier curves.
Results— We identified 790 patients with BAVM (51% female; 63% white; mean age±SD at diagnosis: 38±19 years) between 1961 and 2001. Patients who presented with ICH experienced a higher rate of subsequent ICH than those who presented without ICH under multivariate analysis (hazard ratio, 3.6; 95% CI, 1.1 to 11.9; P<0.032). The effect was similar across race-ethnicity and gender. This difference in ICH rates was greatest in the first year (7% versus 3% per year) and converged over time. The effect of subsequent ICH on functional status was similar to that of the initial ICH.
Conclusions— Presentation with ICH was the most important predictor of future ICH, confirming previous studies. Future ICH had similar impact on functional outcome as incident ICH. Intervention to prevent ICH would be of potentially greater benefit to patients presenting with ICH, although the advantage decreases over time.
Brain arteriovenous malformations (BAVMs) are an important cause of hemorrhagic stroke in young adults.1–4 Treatment decisions can be challenging because of expensive complex treatment (resection, embolization, radiosurgery) and a relative lack of unbiased natural history data.3,4 Accurate estimates for risk and rates of ICH in the natural course of patients harboring BAVM are needed to provide a quantitative basis for planning clinical trials to evaluate interventional strategies and to help guide practice management.
A controversy exists concerning what determines the risk of new hemorrhage if the lesion is left untreated.5,6 A primary unresolved issue is whether clinical presentation with intracranial hemorrhage (ICH) increases the risk of subsequent hemorrhage6,7 or does not affect future bleeding risk.8,9
Our hypothesis was that initial clinical presentation with ICH is associated with a higher rate of subsequent ICH in the natural course before any treatment is undertaken. We present data obtained from a large, community-based, metropolitan health maintenance organization (HMO) that, compared with earlier reports, is less encumbered by selection bias inherent in pure tertiary care patient samples. Additionally, we examine an index of functional outcome from subsequent ICH. These more widely generalizable data confirm earlier reports and make it possible to carefully balance the risks of intervention with conservative medical management.
Patients and Methods
We identified patients with BAVM at Kaiser Permanente Northern California (KPNC) and documented their clinical course using a multimodal approach using several databases from inpatient and outpatient services, radiology reports, and chart review. KPNC is an HMO currently with ≈3.3 million members in a 14-county region. It covers a population representative of the sociodemographic characteristics of the Bay Area, but underrepresents residents at both extremes of the socioeconomic continuum.10
Ascertainment included an initial screen using ICD-9 code 747.81 to identify all diagnoses of intracranial vascular malformation, which includes BAVM. We searched KPNC’s inpatient admissions and outpatient care database. The ICD-9 codes were assigned by coders based on diagnostic notes recorded by physicians during the encounter in the emergency room, outpatient clinic, or hospital admission (including an extramural claims database). Unique identifiers prevented recounting from multiple patient visits.
Beginning with calendar year 1998, complete text reports for all radiological procedures were available. The radiology text files were searched using CPT-4 codes for various neuroradiological procedures and screened for various permutations of specified text strings germane to BAVM, arterial-venous malformation, angioma, etc.
The initial visit for each patient was taken as the time of index date, ie, detection of a case as a study patient from one of the various sources. A trained medical records analyst (MRA) generated a list of all the encounters after the index date and the respective sites for medical care. In cases in which date of diagnosis preceded the index date, the MRA reviewed the records with the first mention of BAVM and all subsequent encounters. Using a close-ended data collection sheet, the MRA reviewed photocopies of initial evaluation, follow-up, and treatment for BAVM in the patient records. Lesion size and morphological data were taken from either progress notes or radiology reports.
The abstracts were reviewed by 1 or more of the study physicians (V.S., W.L.Y., S.C.J.) to determine or confirm the underlying diagnosis. A standardized definition was used as previously described:11 an abnormal tangle of vessels that results in arteriovenous shunting, excluding vein of Galen AVM, cavernous malformations, anomalous venous drainage, ruptured arterial aneurysms, dural arteriovenous fistulae, venous malformations, venous varices, or any of the other rarer types of cerebrovascular anomalies. The diagnosis of BAVM was based on information from multiple sources including clinical history, neuroimaging, and, in a few cases, pathology.
The rate of ICH in the subsequent period after BAVM diagnosis, before any intervention, was assessed using survival analysis. The period at risk was defined from the date of first BAVM diagnosis to the date of an event, ie, onset of (first subsequent) ICH, or censoring because of either initiation of first BAVM treatment (surgery, embolization, or radiosurgery) or loss to follow-up (using date of last available follow-up). Cases were censored at first treatment to account for any possible change in natural history risk incurred by intervention for BAVM.
Patient and BAVM characteristics examined included age at diagnosis, gender, race-ethnicity (white versus nonwhite), ICH at presentation (ICH versus no ICH), venous drainage pattern (exclusively deep versus other), and BAVM size (largest dimension <3 cm versus ≥3 cm). Because morphological characteristics had to be abstracted from the medical record text, only deep venous drainage and size were chosen, which are components of the Spetzler-Martin surgical risk scoring system.12 All variables were coded following published standardized guidelines.11
The respective associations between these variables and subsequent ICH were compared using univariate and multivariate Cox proportional hazards models. Stratified analyses were performed to assess statistical interaction by age, gender, and race-ethnicity (white versus nonwhite).
In constructing the multivariate model, we began by including all selected variables into a full model. In a backward stepwise approach, individual variables were subsequently removed from the model singly or in combination to determine the most parsimonious model.
Cumulative rate of ICH was displayed using Kaplan-Meier curves. Differences between hemorrhage-free survival for patients who presented with and without ICH were analyzed using the log rank test. Cumulative and annualized ICH estimates at specific time points were derived from the Kaplan-Meier survival tables. We also performed a restricted survival analysis involving only those cases ascertained at the time of diagnosis at KPNC and compared these results with those obtained from the entire cohort.
We analyzed change in functional status using Rankin score before and after hospital admission for BAVM during initial presentation, as well as before and after admission for subsequent ICH. Rankin score was not part of the original medical record, but was estimated by the chart reviewer from the medical record for the hospital admission. Use of Rankin score from chart abstraction has been previously described.13 We dichotomized Rankin score to >2 (poor outcome) and ≤2 (good outcome), and we used paired (McNemar) χ2 tests.
We identified 793 patients at KPNC with BAVM diagnosed between 1961 and 2001; 3 patients without information regarding initial presentation were excluded from further analysis. Among all BAVM patients identified, 367 (47%) presented with an ICH. Using a hierarchical system,11 the remaining presentations were 190 (24%) with seizure, 114 (14%) with headache, and 119 (15%) with other symptoms, including incidental ones not definitively attributable to the BAVM.
The distribution of patient and BAVM characteristics risks and their respective univariate associations with subsequent ICH under the Cox proportional hazards model are presented in Table 1. Patients who initially presented with ICH had a higher rate of ICH before first treatment (hazard ratio, 2.5; 95% CI, 1.5 to 4.0; P<0.001). This result persisted even after adjusting for age at diagnosis, race-ethnicity, gender, BAVM size, and venous drainage (hazard ratio, 3.6; 95% CI, 1.1 to 11.9; P=0.032).
Including year of diagnosis and use of catheter angiography into the model separately did not alter the estimate. In the final multivariate model, only presentation with ICH demonstrated an independent relationship with subsequent ICH (Table 2). Increased age at diagnosis was marginally significant (P=0.056). An alternative multivariate model excluding BAVM size resulted in attenuation of the hazard ratio estimate of ICH presentation to 2.6 (95% CI, 1.4 to 5.1), but an increase in significance (P) from 0.032 to 0.004.
A higher proportion of patients from the ICH group were nonwhite (38% versus 25%, P<0.001), but other demographic characteristics were similar. Stratified analyses found no difference in hazard rate ratios of presentation with ICH across race-ethnicity or gender. Further examination of race-ethnic strata was limited by small sample sizes for other subgroups: black (9.0%), Latino (11.5%), Native American (0.1%), and unknown/unrecorded (7.6%).
Follow-up data were available for 350 (95%) of the 367 patients who presented with ICH and 410 (97%) of the 423 patients who presented without ICH (Table 3). In the group that presented with ICH, 32 (9%) had a second ICH before treatment, with 11 of these occurring by the end of the first year; the median time to hemorrhage was 2.4 years. Among patients who presented without ICH, 33 (8%) had an ICH, 8 of the ICH events occurred by the first year, and the median time to hemorrhage was 3.7 years.
The cumulative ICH rates calculated using the Kaplan-Meier method are presented in Table 3 for specific times after diagnosis and are depicted graphically in Figure 1. Compared with the non-ICH group, patients in the ICH group experienced a higher cumulative rate of subsequent ICH (log rank test, P<0.001) and overall had a shorter time at risk and a shorter duration between diagnosis and subsequent ICH. The annualized rates of ICH differed most during the first 2 years but gradually converged over subsequent years (Figure 2).
To address potential ascertainment bias, we restricted the survival analysis to the 572 patients ascertained at time of diagnosis and followed-up prospectively. Overall, the hazard rate ratio for presentation with ICH in this subset of cases was not different from the total cohort in either univariate or multivariate analyses.
Patients who presented with ICH could be assigned Rankin scores more often (345/367) than those with non-ICH presentation (54/423), because they were more likely to have been admitted to the hospital and have associated neurological status detailed in the medical record. All (100%) of the patients who initially presented with ICH had a Rankin score of ≤2 before admission for the initial ICH (Figure 3). This proportion decreased to 67% at discharge (P<0.001), indicating that 33% of the patients had a new deficit after the ICH. From this group, among those who proceeded to have a second ICH, the proportion with a Rankin score of ≤2 decreased from 97% before admission for the second ICH to 66% at discharge (P<0.05), similar to the proportion of patients with new deficits observed for initial ICH. In the non-ICH presentation group, the proportion of patients with a Rankin score of ≤2 decreased minimally from 94% before admission to 89% at discharge. However, at the time of subsequent ICH, the proportion of patients with a Rankin score of ≤2 decreased from 100% before admission to 56% at discharge, similar in magnitude to the change seen in the ICH presentation group.
This study demonstrates that ICH at initial presentation is the most important predictor of ICH in the natural course of patients harboring a BAVM. It also confirms and elaborates upon previously reported findings (Figures 1 and 2⇑). Our findings offer the first corraborative results from a large community-based, rather than a pure referral, population. Our findings also reconcile the apparent disagreement between studies suggesting differing effects of incident ICH on subsequent ICH: we found an initial difference in ICH rates that diminishes over time and then shows a longer-term trend for the ICH rates to converge (Figure 2). Our examination of the change in Rankin score offers indirect evidence that subsequent ICH is functionally significant to the patient, and the magnitude of the decrease in functional status is similar to that observed for the initial ICH (Figure 3).
Our results confirm and extend previous retrospective observational studies that estimate general risk of spontaneous bleeding to be in the range of 2% to 6% per year,3,14 probably higher in the first year.15,16 Taken together, the evidence seems to indicate that surgical intervention, if decided on, would be more appropriate sooner rather than later given the steep shape of the hemorrhage-free survival curve during the initial period, especially for patients presenting with ICH (Figures 1 and 2⇑). However, this decision must be weighed against the risk of surgical treatment itself.14
The increase in ICH rate was not modified by race-ethnicity or gender and was not biased by the potentially selective ascertainment among cases diagnosed outside of KPNC and ascertained later along the course of their disease. We could not examine the impact of race-ethnicity in detail because of relatively small sample sizes of race-ethnic groups other than white.
To demonstrate an independent effect of ICH presentation on subsequent ICH, we excluded 1 or more BAVM characteristics associated with hemorrhage from the full multivariate model but could not arrive at a more parsimonious model. Because only 30% of all cases had information on lesion size, we also considered a model excluding BAVM size to see the effect with a more complete dataset. This alternative model also demonstrated the independent effect of presentation with ICH.
Patient selection for surgical treatment could introduce bias into the observed ICH rates in several ways, including introduction of unmeasured covariates. Addition of 2 critical components of the Spetzler-Martin surgical risk score to our multivariate model did not change the nature of the relationship between presentation and subsequent ICH.
Actual patient images were not available to obtain expert neuroradiological adjudication to confirm lesion size or detailed morphological characteristics. The descriptors we report were abstracted from the medical record and were based on published guidelines.11 Therefore, they serve primarily to verify the population that we present is comparable to populations previously reported. There are a number of important considerations of specificity, sensitivity, and interobserver agreement in the evaluation of BAVMs1,17 that the current study design cannot accommodate.
Our hazard ratio estimate was lower than that reported by Mast et al (adjusted hazard rate ratio: 13.9).6 This may simply reflect the different BAVM populations seen at an HMO as compared with a tertiary care center, from which Mast et al derived their cases. Two large studies failed to show an effect of ICH presentation on subsequent ICH. One was geographically and temporally dissimilar from the present setting,9 and the other included a majority who were undergoing a course of embolization and, therefore, did not address pure natural history.
We used Rankin score as a surrogate for functional status to address the question of whether the consequences of a subsequent ICH were worse than an initial ICH. Our results offer a preliminary evaluation of functional changes in subsequent ICH, suggesting that subsequent ICH in the natural course was associated with a decline in functional status, similar to the decline associated with an initial ICH.
The data for Rankin scoring may be influenced by the time elapsed between the ICH and the outcome evaluation; there is often considerable improvement over time after BAVM hemorrhage. Our data represent the functional outcome at hospital discharge for the ICH event. Therefore, our data reflect a short-term outcome, which may overestimate the eventual long-term loss of function.
Some cases may have been missed in our ascertainment, eg, as a result of fatal first ICH or ICD-9 miscoding. However, for this analysis, it was more important that we did not miss subsequent ICH events, and our ability to track inpatient, outpatient, and extramural visits allowed us to detect all subsequent events in our cohort. There are at least subtle effects of the time period of diagnosis on our results, because evaluation and management have changed considerably over the past several decades. Nonetheless, adjusting our analysis for year of ascertainment did not affect our conclusions.
We are aware of 2 large BAVM cohorts being assembled.2,18,19 Data from these sources will be helpful in validating prediction rules derived from our results, which may be of use in formulating treatment decision algorithms and can assist in planning clinical trials.
The authors thank Carolyn Salazar, Barbara Rowe, and Mike Sorel of KPNC; Carroll Schreibman, Broderick Belenson, Mary Nelson Tran, PhD, and the UCSF Brain Arteriovenous Malformation Study Project (http://avm.ucsf.edu). Supported in part by PHS grants RO1 NS34949 and K24 NS02091 (W.L.Y.) and a UCSF School of Medicine Centers of Excellence grant, and K08 NS42359 (A.X.H.) and K23 NS42072 (V.S.) from the National Institutes of Health/National Institute of Neurological Disorders and Stroke.
- Received December 5, 2003.
- Accepted February 2, 2004.
Al-Shahi R, Warlow C. A systematic review of the frequency and prognosis of arteriovenous malformations of the brain in adults. Brain. 2001; 124: 1900–1926.
Stapf C, Mast H, Sciacca RR, Berenstein A, Nelson PK, Gobin YP, Pile-Spellman J, Mohr JP. The New York Islands AVM Study. Design, study progress, and initial results. Stroke. 2003.
Ogilvy CS, Stieg PE, Awad I, Brown RD Jr, Kondziolka D, Rosenwasser R, Young WL, Hademenos G. Recommendations for the management of intracranial arteriovenous malformations: a statement for healthcare professionals from a special writing group of the Stroke Council, American Stroke Association. Circulation. 2001; 103: 2644–2657.
Crawford PM, West CR, Chadwick DW, Shaw MD. Arteriovenous malformations of the brain: natural history in unoperated patients. J Neurol Neurosurg Psychiatry. 1986; 49: 1–10.
Stefani MA, Porter PJ, terBrugge KG, Montanera W, Willinsky RA, Wallace MC. Large and deep brain arteriovenous malformations are associated with risk of future hemorrhage. Stroke. 2002; 33: 1220–1224.
Joint Writing Group of the Technology Assessment Committee American Society of Interventional and Therapeutic Neuroradiology; Joint Section on Cerebrovascular Neurosurgery, a section of American Association of Neurological Surgeons and Congress of Neurological Surgeons; and Section of Stroke and the Section of Interventional Neurology of the American Academy of Neurology. Reporting terminology for brain arteriovenous malformation clinical and radiographic features for use in clinical trials. Stroke. 2001; 32: 1430–1442.
Al-Shahi R, Pal N, Lewis SC, Bhattacharya JJ, Sellar RJ, Warlow CP. Observer agreement in the angiographic assessment of arteriovenous malformations of the brain. Stroke. 2002; 33: 1501–1508.
Al-Shahi R, Bhattacharya JJ, Currie DG, Papanastassiou V, Ritchie V, Roberts RC, Sellar RJ, Warlow CP. Scottish Intracranial Vascular Malformation Study (SIVMS). Evaluation of methods, ICD-10 coding, and potential sources of bias in a prospective, population-based cohort. Stroke. 2003; 34: 1156–1162.
Al-Shahi R, Bhattacharya JJ, Currie DG, Papanastassiou V, Ritchie V, Roberts RC, Sellar RJ, Warlow CP. Prospective, population-based detection of intracranial vascular malformations in adults: the Scottish Intracranial Vascular Malformation Study (SIVMS). Stroke. 2003; 34: 1163–1169.