Correlations between Perfusion MR Imaging Cerebral Blood Volume, Microvessel Quantification, and Clinical Outcome Using Stereotactic Analysis in Recurrent High-Grade Glioma

BACKGROUND AND PURPOSE: Quantifying MVA rather than MVD provides better correlation with survival in HGG. This is attributed to a specific “glomeruloid” vascular pattern, which is better characterized by vessel area than number. Despite its prognostic value, MVA quantification is laborious and clinically impractical. The DSC-MR imaging measure of rCBV offers the advantages of speed and convenience to overcome these limitations; however, clinical use of this technique depends on establishing accurate correlations between rCBV, MVA, and MVD, particularly in the setting of heterogeneous vascular size inherent to human HGG. MATERIALS AND METHODS: We obtained preoperative 3T DSC-MR imaging in patients with HGG before stereotactic surgery. We histologically quantified MVA, MVD, and vascular size heterogeneity from CD34-stained 10-μm sections of stereotactic biopsies, and we coregistered biopsy locations with localized rCBV measurements. We statistically correlated rCBV, MVA, and MVD under conditions of high and low vascular-size heterogeneity and among tumor grades. We correlated all parameters with OS by using Cox regression. RESULTS: We analyzed 38 biopsies from 24 subjects. rCBV correlated strongly with MVA (r = 0.83, P < .0001) but weakly with MVD (r = 0.32, P = .05), due to microvessel size heterogeneity. Among samples with more homogeneous vessel size, rCBV correlation with MVD improved (r = 0.56, P = .01). OS correlated with both rCBV (P = .02) and MVA (P = .01) but not with MVD (P = .17). CONCLUSIONS: rCBV provides a reliable estimation of tumor MVA as a biomarker of glioma outcome. rCBV poorly estimates MVD in the presence of vessel size heterogeneity inherent to human HGG.

A lthough the degree of angiogenesis in glioma generally correlates with malignant potential, histologic outcome studies have shown that the specific method of quantifying angiogenesis can greatly impact the ability to predict survival. [1][2][3][4][5][6][7][8][9][10] Most glioma studies have measured angiogenesis by using vessel number counts or MVD. [11][12][13][14] Both precedent and the relative ease of MVD calculation have made this technique a logical first choice to study outcome; however, while MVD shows good predictive power for non-central nervous system tumors, the correlation with glioma survival has been marginal at best. [1][2][3][4][5][6] Specific characteristics of glioma angiogenesis are likely to account for the poor correlations with survival. HGGs are often characterized as having large-lumen microvessels that are few in number, large in area, and linked with aggressive behavior and poor prognosis. 4,6-8 These large-lumen microvessels can be composed of either a single lumen or multiple lumina, the latter having a "glomeruloid" appearance that is a histologic hallmark of grade IV tumors. MVD poorly characterizes the morphometric complexity of these large-lumen microvessels, because area and number are often inversely related. 1,4,[6][7][8]14,15 In contrast, the MVA of the tumor vessels may provide a more salient clinical biomarker of glioma outcome. This assertion is supported by studies showing that MVA correlates more strongly with survival compared with MVD, even in tumors with similar histologic grades. 1,4,[8][9][10] Despite this evidence, histologically quantifying MVA remains a research tool because the technique is time-consuming, labor-intensive, and thus impractical for routine clinical use. 6,7 For Ͼ2 decades, DSC-MR imaging has been used to study glioma angiogenesis. [16][17][18][19][20][21] On the basis of the indicator dilution theory, DSC-MR imaging calculates rCBV to estimate tissue blood volume. 18,21,22 Because DSC-MR imaging is semiautomated, it offers advantages of speed and convenience that overcome limitations of histologic vessel quantification. Because measurements can also be performed serially and noninvasively without surgical biopsy, DSC-MR imaging has the potential to lower cost, morbidity, and sampling error to better provide measures of outcome as well as treatment efficacy in patients with glioma. The clinical translation of DSC-MR imaging as a feasible alternative to histologic quantification requires correlations between rCBV, MVD, and MVA. To date, nearly all DSC-MR imaging correlative data are from animal models by using only MVD. 21,[23][24][25][26][27] Although MVA likely represents a more clinically relevant parameter, 1,4,6-10 no human studies have been done and only 1 group has reported MVA correlation with rCBV in the rat 9L gliosarcoma model. 22,28 Basing the clinical translation of DSC-MR imaging primarily on animal data has several distinct limitations. Animal models often lack the broad variations in vessel growth pattern and size that are observed in human HGG, [6][7][8]21,22,29 and it is unclear how the heterogeneity of human tumors impacts rCBV correlations with MVD and MVA. Additionally, largearea vessels in human HGG cause MVD and MVA to be inversely related [6][7][8][9][10] ; this feature would directly challenge the consensus from animal studies that rCBV correlates well with both measures. [21][22][23][24][25][26][27][28] Finally, the rationale for using rCBV to predict glioma outcome has not been reconciled with extensive clinical data comparing histologically derived MVD and MVA with survival. Specifically, whether rCBV more strongly reports MVD versus MVA directly impacts the utility of rCBV as a clinical predictor of glioma outcome. 4,6-10 That these relationships are currently unproved in humans might explain why DSC-MR imaging has not yet gained significant traction as a part of clinical standard of care, with the exception of select academic institutions. 30 To address these issues, we present data from a pilot study correlating localized coregistered rCBV with histologically derived MVD and MVA from stereotactic biopsies in human HGG. We hypothesized that rCBV would show stronger correlation with tissue MVA compared with MVD due to vascular patterns that cause vessel size heterogeneity. Similarly, we hypothesized that rCBV and MVA would show similar correlations with survival in a cohort of patients with HGG, whereas MVD would not. Thus, our aims in this study were the following: 1) to compare the strength of rCBV correlation with MVD and MVA, 2) to describe how vascular heterogeneity (ie, size variation) inherent in human HGG impacts these correlations; and 3) to investigate the associations between rCBV, MVD, and MVA as biomarkers of survival.

Subject Recruitment
We recruited patients with previous HGG who were to undergo new operations. We studied them during preoperative imaging performed for the stereotactic surgery procedure. We documented previous treatment (including chemoradiation therapy and/or previous surgical resection) and excluded subjects with an estimated glomerular filtration rate Ͻ60 mg/min/1.72 m 2 . We obtained approval from the institutional review board and written informed consent from each subject.

Preoperative DSC-MR Imaging Protocol
We used the same 3T MR imaging magnet (Sigma HDx; GE Healthcare, Milwaukee, Wisconsin) for all patients. All imaging was performed within 1 day of surgery. An intravenous injection of 0.1 mmol/kg of Gd-based contrast agent (gadodiamide or gadobenate dimeglumines, preload dose) was administered before DSC-MR imaging on the basis of a previously published optimized protocol. 31,32 The DSC-MR imaging was performed by using a 0.05-mmol/kg Gdagent bolus delivered at a rate of 3-5 mL/s. This dosage, at high field strength, provides strong signal intensity to noise while minimizing patient contrast load. 33,34 The DSC-MR imaging sequence parameters were the following: gradient-echo echo-planar imaging with TR/ TE/flip angle, 1500 -2000 ms/20 ms/60°; FOV, 24 ϫ 24 cm; matrix, 128 ϫ 128; 5-mm section; no intersection gap. The total Gd agent dose was 0.15 mmol/kg of body weight.

Intraoperative Tissue Specimen Collection
Each neurosurgeon collected 2-3 tissue specimens from each tumor by using stereotactic surgical localization, following the smallest possible diameter craniotomies to minimize brain shift. [34][35][36] We chose not to use maximum or minimum rCBV values to guide biopsy targeting, to avoid biasing the correlation between rCBV, MVA, and MVD. Instead, biopsy targets were collected by using standard clinical criteria, such that specimens were collected in pseudorandom fashion from different poles of the enhancing lesion periphery while avoiding any necrotic regions. Biopsies were performed without knowledge of the DSC-MR imaging analyses. The locations of the biopsy sites and the neuronavigational coordinates were recorded at the time of surgical tissue collection as screen captures and were later coregistered with the MR imaging datasets to enable localized measurement of rCBV at the sites for each corresponding tissue sample. Multiple biopsy targets in the same patient were separated by a minimum of 2 cm. The neurosurgeon visually validated stereotactic imaging locations with corresponding intracranial anatomic landmarks, such as vascular structures and ventricle margins, before recording specimen locations.

Coregistration of Stereotactic and DSC-MR Imaging Datasets
Stereotactic and DSC-MR imaging were transferred to an off-line Linux-based workstation for rigid-body coregistration by using statistical parametric mapping (SPM5, Wellcome Department of Imaging Neuroscience, University College London, UK). DSC-MR imaging coregistration enabled region-of-interest placement for rCBV calculation within previously recorded stereotactic locations of corresponding surgical tissue specimens, as previously described. [33][34][35] All regions of interest were also inspected to ensure exclusion of visible large vessels.

DSC-MR Imaging Data Postprocessing and rCBV Calculation
DSC-MR imaging was analyzed by using an in-house Matlab-based program (MathWorks, Natick, Massachusetts) based on previously published methods. 18,31,32,34 In short, we generated whole-brain CBV maps by integrating the first-pass ⌬R2*(t), using baseline subtraction to correct residual T2/T2*-weighted effects. 31,32 Following CBV normalization to contralateral normal gray and white matter as previously described, 34 we calculated mean rCBV for each tissue specimen from a 3 ϫ 3 voxel (ϳ0.8 cm 2 ) region of interest centrally placed within corresponding stereotactic locations. All analyses were performed without knowledge of tissue analysis.

Tissue Specimen Preparation and Histopathologic Diagnosis
All stereotactic biopsy specimens were fixed in 10% formalin and paraffin-embedded. CD34 immunohistochemistry was performed on 10-m sections by using the Ventana Nexus immunostainer (Ventana Medical Systems, Tucson, Arizona). A prediluted monoclonal primary antibody (Mouse Anti-Human CD34, clone QBEnd/10, Ventana Medical Systems) was used with standard cell conditioning. The reaction product was visualized with diaminobenzidine. Negative control sections were processed identically with omission of the primary antibodies. 5,9 Corresponding sections were stained with hematoxylin-eosin per standard protocol. Specimens were diagnosed histologically according to the 2007 World Health Organization classification. 37

Tissue Microvascular Parameter Quantification
Whole-section slides from each specimen were imaged and digitized by using an Axiovert 200M light microscope, Achroplan 20x Objective (Zeiss, Gottingen, Germany) and AxioVision Mosaix 4.8 (Zeiss). Whole-section quantitative microvascular analysis was performed by using the computer-assisted image analysis software AxioVision Automeasure 4.8 (Zeiss). Microvessel parameters were defined as previously described on the basis of CD34-stained profiles, as per the clinical protocol at our institution as well as techniques from multiple previous studies. 1,6,7,10,22 The mean smallest microvessel diameter was determined by using the ferret minimum function. The total area (square micrometers) of each section was also calculated and used to normalize respective microvascular parameters. Tissue analyses were supervised by 2 neuropathologists.

Statistical Analysis
We performed Pearson correlations between rCBV and microvascular parameters for all corresponding stereotactic biopsy samples, as demonstrated in Fig 1. We chose to perform linear correlations based on prior reports, which would facilitate the comparison of correlations from our study with those previously performed. 22,38 For clinical outcome analyses, we determined OS from the date of surgery to death for those patients who had died at the time of analysis. For each patient, rCBV, MVA, and MVD were determined across biopsy locations, and the maximum values (1 value per patient for each variable) were correlated with OS by using Cox regression survival analyses for all patients. The Cox regression analysis takes into account follow-up for all patients, not just for outcomes from which patients died. Patients who have not died contribute data until their last date known to be alive. To assess the influence of tumor grade (III or IV), we compared rCBV, MVA, and MVD between grade III and grade IV tumor samples by using the Fisher exact test. We also performed a Cox regression analysis, comparing tumor grade with OS. A biostatistician performed all analyses, and statistical significance was indicated by P Ͻ .05.

Study Subjects and Tissue Specimens
We enrolled 27 subjects and collected a total of 75 tissue specimens, of which 54 contained sufficient tissue for both histopathologic diagnosis and microvascular quantification. Three samples were excluded due to poor CD34 staining quality, primarily seen as diffuse background staining. Also, 13 samples were diagnosed as having radiation effect without evidence of tumor and were therefore excluded. Final analyses were based on 38 HGG tissue specimens from 24 patients containing either grade III (n ϭ 17) or grade IV (n ϭ 21) tumor. The average tissue specimen volume was approximately 0.2 cm 3 with an average tissue cross-sectional area of 1.61 ϫ 10 7 m 2 . All patients had undergone previous surgical diagnosis for their primary tumor. The Table summarizes the original histopathologic diagnosis and the type and timing of prior chemotherapy and radiation treatments. In our study, the average postprocessing time for rCBV calculation was approximately 5 minutes, compared with approximately 45 minutes for histologic multiparametric vascular quantification. One grade III tumor sample demonstrated glomeruloid vessels. Because this patient received prior multimodality therapy (including radiation), the presence of glomeruloid vessels could not be entirely attributed to malignant progression; thus, on the basis of clinical diagnostic criteria at our institution by 2 board-certified neuropathologists, the tumor was not histopathologically upstaged to grade IV but rather maintained as grade III. There were no significant differences in rCBV (P ϭ .22), MVA (P ϭ .16), or MVD (P ϭ .42) between grade III and grade IV tumor samples.

rCBV Correlation with MVA and MVD
Within the entire specimen group (N ϭ 38), the mean MVA value was 0.0735 Ϯ 0.062 with a range of 0.008 -0.26 and the mean MVD value was 8.40E-4 Ϯ 1.11E-3 m Ϫ2 with a range of 5.07E-5-4.62E-3. The mean rCBV value was 1.49 Ϯ 0.97 with a range from 0.38 to 4.64. rCBV correlated strongly with MVA (r ϭ 0.83; P Ͻ .001; 95% confidence interval, r ϭ 0.69 -0.90). Measures of rCBV weakly correlated with MVD over the entire group (N ϭ 38, r ϭ 0.32; P ϭ .05; 95% confidence interval, r ϭ 0.00 -0.58). Confidence intervals for rCBV correlations between MVA and MVD did not overlap. To test for possible intrasubject grouping effects of multiple tissue samples collected from a single subject, we performed a bootstrap analysis in which rCBV/MVA correlations were computed with random selection of 1 specimen per subject. The correlation remained consistent with the original analysis (average r ϭ 0.78, standard error ϭ 0.002).

Effect of Heterogeneity of Vascular Size on rCBV Correlation with Microvessel Area and Number
The rCBV correlation with MVD in our study differed from that of previous animal reports. To explain this difference, we examined the variation in vascular size within tissue samples. We postulated that because heterogeneity of vascular size is more commonly observed in human HGG than in animal models, greater heterogeneity might negatively impact the rCBV correlation with MVD in our dataset. To test this effect, we separated specimen sections into subgroups based on a vascular size heterogeneity index, representing the ratio between the SD of the vessel diameter within each sample and the mean vessel diameter. A larger SD relative to vessel size suggests greater heterogeneity of vessel size. 6,7 We selected the index value of 1.0 to yield separation between high (Ͼ1) and low (Յ1) vessel size heterogeneity. Vascular heterogeneity impacted MVD correlations with rCBV; sections with low vascular size heterogeneity showed better correlation between rCBV and MVD (n ϭ 20, r ϭ 0.56, P ϭ .01) compared with the high vascular size heterogeneity group (n ϭ 18, r ϭ 0.20, P ϭ .42). MVA correlation with rCBV remained strong for both high (r ϭ 0.81, P Ͻ .001) and low (r ϭ 0.82, P Ͻ .001) vessel size heterogeneity groups Fig 2.

Microvessel Parameter and rCBV Correlation with OS
At the time of analysis, 12 of 24 subjects had died. One patient who had died had received bevacizumab before imaging and surgery. Bevacizumab has been shown to reduce tissue MVA, MVD, and rCBV values. 22,24,39 Because rCBV represents an imaging measure of tissue microvasculature both in the ab-sence and presence of antiangiogenic therapy, 22 this patient was not excluded from correlations between rCBV, MVA, and MVD as measures of vascularity. However, because bevacizumab does influence the use of MVA, MVD, and rCBV as biomarkers of malignant potential, we excluded this patient from correlations with outcome. Of the patients surviving at the time of analysis, the median time of clinical follow-up was 9.3 months. Cox regression analysis of 23 subjects showed that survival was associated with both MVA (P ϭ .01) and rCBV (P ϭ .02) but did not reach statistical significance with MVD (P ϭ .17). Among this cohort of patients with HGG, histologic grade (III or IV) was associated with survival (P Ͻ .05).

Discussion
Although histologic grade undoubtedly impacts prognosis, multiple histologic studies have shown that MVA represents an independent predictor of glioma outcome beyond grade alone. Birner et al 8 characterized the microvessel morphometry in 114 GBM (grade IV) tumors, excluding lower grade tumors from analysis. They classified GBM tumors as having the predominance of either glomeruloid vessels (representative of high MVA) or delicate thin-lumen capillary-type vessels (representative of low MVA). Their data showed that GBM tumors lacking glomeruloid vessels were associated with longer patient survival, suggesting that tumoral microvasculature predicts survival differences among tumors with identical histologic grade. Although glomeruloid vessels are considered a histologic hallmark of grade IV tumors, elevated MVA can be observed in grade III tumors as a result of large-area singlelumen microvessels that are not technically considered "glomeruloid." 1,4,8,10 In fact, Korkolopoulou et al 4,10 reported that grade III tumors, even in the absence of glomeruloid vessels, could exhibit MVA that exceeded that of many GBM tumors and supported the findings of Birner et al 8 by showing that elevated MVA was an independent prognostic biomarker of glioma outcome. That MVD may not accurately report the presence of large-area microvessels offers an explanation for why MVD more weakly correlates with outcome compared with MVA. 1,4,[8][9][10] On the basis of prior histologic studies, it is evident that MVA represents a useful prognostic marker of glioma. Despite its potential value, however, MVA calculation is clinically impractical because of the time-consuming and labor-intensive nature of deriving quantitative measures. Because the calculation of rCBV is both fast and noninvasive compared with that of MVA, a primary goal of the current study was to establish whether rCBV represents a reliable surrogate of MVA in the clinical setting. We used careful stereotactic coregistration to show that rCBV correlates highly with MVA and generally poorly with MVD. We also showed that the average postprocessing time for rCBV calculation was approximately 5 minutes, compared with approximately 45 minutes for histologic multiparametric microvascular quantification. Knowledge of these histologic correlations provides justification for the use of rCBV as a clinically feasible alternative to tissue MVA quantification as a predictor of glioma outcome.
The results from this study can be compared with those of previous animal work by Pathak et al,22 who perfused rat 9L gliosarcoma brains with a latex compound (Microfil; Flow Tech, South Windsor, Connecticut) to mark microvessels for histologic area measurement. Although our human tissue specimens demonstrated a comparatively greater range and heterogeneity of vessel sizes (Յ125 m), the human (r ϭ 0.83) and animal (r ϭ 0.74) rCBV correlations with MVA remained similarly high, suggesting resilience in the clinical setting. The rCBV correlation with MVA can be explained by the Delesse principle, which states that the "areal attenuation of profiles on sections is an unbiased estimator of the volume attenuation of structures." 22 were substantially thinner. This difference may have improved our correlation. Another difference between the 2 studies is that Pathak et al used MION contrast agent to prevent contrast leakage effects. Because MION is contraindicated in humans, we used preload contrast administration, which robustly minimizes T1-weighted leakage effects, regardless of concomitant steroid therapy. 18,31,32 In our study, rCBV weakly correlated with MVD within a general population of heterogeneously sized microvessels; however, when the degree of heterogeneity was minimized, rCBV correlation increased. This "conditional correlation" with MVD likely explains the discordance with results from numerous animal tumor models, which generally demonstrate more homogeneous ranges of vessel size compared with human glioma. 21,[23][24][25][26][27][28][29] Thus, rCBV seems to represent a sufficient measure of MVD (ie, vessel number) in experimental animal models but a poor measure of MVD in the clinical setting. Our results challenge previous clinical reports of high rCBV correlation with MVD in human HGG. 38,41 These discrepancies likely reflect methodologic differences both in imaging-tissue correlation and statistics. For instance, Haris et al 38 did not use spatial coregistration, so regions of rCBV did not necessarily spatially correlate with those of MVD measurement. They chose to match values by rank-order grouping, assuming that high-rCBV locations corresponded to high-MVD locations; this correspondence artificially elevated Pearson correlations. Stereotactic coregistration between imagingtissue locations in our study precluded the need for rank-order grouping to match values and provided the ability to directly compare rCBV correlations with multiple histologically derived vascular parameters.
Based on conventional theory, the prognostic value of rCBV relates primarily to its prediction of glioma grade, which, in turn, drives clinical outcome. [16][17][18][19] However, recent studies have proposed an alternative hypothesis, that rCBV predicts outcome independent of histologic grade. 20,42,43 Law et al 20 and Hirai et al 43 separately reported that rCBV could distinguish survival differences among tumors that were diagnosed histologically as the same grade; however, unlike histology-based studies that correlate MVA and grade within the same tissue sample, none of the imaging-based studies to date have measured rCBV from locations that spatially corresponded with the areas of histologic diagnosis. Thus, these imaging-based studies could not exclude the possibility of sampling errors and histologic misgrading that might influence survival differences within their cohorts. 44 Our study provides a direct histologic microvascular correlate to rCBV, which helps to support the alternative hypothesis that rCBV can predict survival beyond grade alone. 20,42,43 Specifically, we validate the hypothesis that rCBV reflects tissue MVA, which permits the unification of otherwise separate sources of literature: those studies correlating survival with histologic MVA quantification 1,4,[8][9][10] and those with rCBV measurement. 16,20,42,43 Unlike previous studies that restricted analyses to tumors of identical grade, the size of our cohort was too small to separately assess survival in grade III and grade IV tumors; however, our purpose in this study was not to replicate what has been shown in prior studies. Rather, the goals of our correlation analyses were to uniquely illustrate the following: 1) rCBV and MVA measures, taken from the same Fig 2. A, Measures of rCBV show poor correlation with MVD under conditions of high vessel size heterogeneity (n ϭ 18, black data series) but good correlation (red data series) when limiting analysis to homogeneous vessel size (n ϭ 20, red series). This "conditional correlation" suggests that rCBV is an unreliable measure of MVD in human HGG, which commonly demonstrates high vessel size heterogeneity. B, MVA strongly correlates with rCBV for both high (n ϭ 18, black series) and low (n ϭ 20, red series) vessel sizeϪheterogeneity groups, suggesting the clinical robustness of rCBV to quantify vessel area in human HGG. Of note, the 6 rightmost data points were measured from 6 different patients.
tumor locations, correlated with each other; 2) rCBV and MVA measures, taken from the same tumor locations, correlated similarly with survival; and 3) MVD correlated relatively poorly with both rCBV and survival.
We derived rCBV from a gradient echoϪbased DSC-MR imaging sequence, which maintains sensitivity to microvessels of all sizes. Other groups have reported spin-echoϪbased rCBV values, which have a maximal sensitivity to small-diameter microvessels (Ͻ25 m). 16,45,46 In our tissue specimen analysis, many of the large-area microvessels, especially in the GBM, often reached diameters of 75-125 m. Because spinecho rCBV values might under-represent these vessel size ranges, 45,46 our choice of gradient-echo DSC-MR imaging seemed appropriate for this study. Future studies by using dual-echo gradient-echo/spin-echo techniques or microvessel tortuosity metrics might help confirm the present findings and also provide insight about microvessel-size heterogeneities. 46,47 We used the region-of-interest technique primarily for correlative purposes, to enable localized rCBV measurements at locations that corresponded to MVA and MVD quantification from multiple stereotactic biopsies. In our study, the biopsy locations dictated the region-of-interest locations. We do not advocate the application of this method to prospective analysis of DSC-MR imaging data in the absence of biopsy because this would require user-determined random selection of region-of-interest locations. We, instead, would advise consideration of more automated methods of rCBV analysis, such as maximum rCBV or histogram techniques, though comparative studies are likely required, determining which of these methods is most predictive of survival. 20,48 There are some small measurement differences between rCBV and microvessel area. First, rCBV represented a summation of physiologic and biophysical processes, such as flow dynamics and contrast-tissue interactions, which may have conveyed in vivo microcirculation information differently from histologic structures in isolation. 31,32,47 For example, the presence of hyalinized vessels from prior radiation might confound the correlation between rCBV and microvessel parameters. In general, however, we focused analysis on tissue samples with predominant histologic evidence of tumor, and we excluded from analysis those tissue samples showing predominant radiation effect that might lead to hyalinized vessels. This greatly reduced the presence of hyalinized vessels in our cohort. Second, we performed microvascular analysis from a single representative section for each specimen. Future studies using 3D techniques to quantify microvascular parameters in the entire tissue specimen might improve correlations. 49 Third, we calculated rCBV by normalizing lesion DSC-MR imaging measures with contralateral normal parenchyma, though the tumor specimen microvascular parameters were not similarly normalized. Fourth, we measured rCBV in vivo before surgery and compared it with ex vivo tissue analysis. We, therefore, cannot exclude the effects of anesthesia or surgery on microvessel contractility or possible distortions in microvasculature during tissue treatment (ie, formalin fixation, sectioning, and staining). Future in vivo microscopy studies in humans might refine correlations. 39 In this study, patients were scanned with either gadodiamide or gadobenate contrast agents, depending on the clini-cal protocol used by our institution at the time of enrollment. These agents have differing T1 relaxivities, which, theoretically, might become evident through T1-weighted leakage effects and errors in rCBV measurement. 18,31,32,34,46,50 To a large extent, we have greatly minimized these T1-weighted leakage errors through use of preload dosing. 31,34,46 Compared with gadodiamide, gadobenate demonstrates approximately twice the in vivo T1 relaxivity at equivalent doses, which suggests that there would be a greater T1 leakage effect with gadobenate. 50 However, at the same time, the preloading dose of gadobenate may be more effective as well. 31,45 This relationship may be confounded by gadobenate protein binding, which can offset the degree of extravascular contrast extravasation.
Despite the use of different contrast agents in our study, rCBV ultimately demonstrated high correlation with MVA (r ϭ 0.83, P Ͻ .0001). This suggests that preload dosing and DSC-MR imaging are robust techniques and that any residual differences in rCBV estimation that may exist between the contrast agents are unlikely to be clinically significant.
We recognize potential limitations regarding coregistration accuracy between DSC measures and stereotactic locations during open resection. 31,34 Possible discrepancies between imaging region-of-interest size and tissue specimen volumes might have confounded correlation, though we conservatively chose 3 ϫ 3 voxel regions of interest to roughly correlate with target tissue specimen volumes and to help optimize DSC localization. Image distortions and brain shift following craniotomy could also lead to misregistration errors. To compensate, neurosurgeons used small craniotomy sizes to minimize brain shift and also used visually validated stereotactic image location with intracranial neuroanatomic landmarks to help correct for random brain shifts. Rigid-body coregistration of stereotactic and DSC-MR imaging also helped reduce possible geometric distortions. 51 Overall, our experience suggests combined misregistration at approximately 1-2 mm from both brain shift and registration technique, which is similar to that from previous studies by using stereotactic needle biopsy. 35,36 Conclusions We used stereotactic coregistration to demonstrate rCBV as an accurate measure of MVA and propose this correlation as a likely factor in the strong prediction by rCBV of clinical outcome. In contrast, MVD may be a less useful predictor of patient outcome in the clinical setting, compared with MVA or rCBV, and weakly correlates with rCBV due to vascular patterns that cause vessel size heterogeneity.