Magnetic resonance imaging biomarkers for clinical routine assessment of microvascular architecture in glioma

Abstract

Knowledge about the topological and structural heterogeneity of the microvasculature is important for diagnosis and monitoring of glioma. A vessel caliber and type-dependent temporal shift in the magnetic resonance imaging signal forms the basis for vascular architecture mapping. This study introduced a clinically feasible approach for assessment of vascular pathologies in gliomas using vascular architecture mapping. Sixty consecutive patients with known or suspected gliomas were examined using vascular architecture mapping as part of the routine magnetic resonance imaging protocol. Maps of microvessel radius and density, which adapted to the vasculature-dependent temporal shift phenomenon, were calculated using a costume-made software tool. Microvessel radius and density were moderately to severely elevated in a heterogeneous, inversely correlated pattern within high-grade gliomas. Additionally, three new imaging biomarkers were introduced: Microvessel type indicator allowing differentiation between supplying arterial and draining venous microvasculature in high-grade gliomas. Vascular-induced bolus peak time shift may presumably be sensitive for early neovascularization in the infiltration zone. Surprisingly, curvature showed significant changes in peritumoral vasogenic edema which correlated with neovascularization in the tumor core of high-grade gliomas. These new magnetic resonance imaging biomarkers give insights into complexity and heterogeneity of vascular changes in glioma; however, histological validations in more well-defined patient populations are required.

Keywords

Magnetic resonance imaging human brain perfusion microvasculature imaging biomarker

Introduction

Gliomas comprise a group of tumors that originate from glial cells in the central nervous system (CNS). High-grade gliomas (HGGs) include anaplastic glioma (WHO grade III) and glioblastoma (GBM, WHO grade IV), which are the most common and aggressive form of glioma with a median survival of only 14 months.¹ GBMs are among the most vascularized of all solid tumors, and vascular proliferation is a pathological hallmark of HGGs.² Tumor vessel formation, however, is highly defective, resulting in an abundant and aberrant vasculature with vessels of variable diameter, heterogeneous distribution, and increased blood–brain barrier (BBB) permeability.³ Neovascularization in gliomas is of high clinical relevance as it correlates with biological aggressiveness, degree of malignancy, clinical recurrence, and the postoperative survival.^4,5

Magnetic resonance imaging (MRI) including measures of perfusion and blood volume is the modality of choice for the non-invasive examination of brain tumors. Computer simulations and experiments in animals and humans have demonstrated that more advanced MRI techniques can estimate microvascular vessel caliber, thereby providing further insight into tissue microvascularity.^6–9 The physical basis for MR-based vessel caliber estimation is the different sensitivity of gradient-echo (GE) and spin echo (SE) MRI to magnetic susceptibility, i.e. the degree of magnetization induced in a material in response to an applied magnetic field.^6–8 The highly susceptibility-sensitive GE signals are dominated by both microscopic and macroscopic vessels, whereas SE signals are predominantly sensitive to microscopic vessels (radius < 10 µm).^6,10,11

Kiselev et al.⁸ were the first who reported about a “loop” in the parametric plot of the relaxation rates for GE ( $R 2, GE$ ) versus SE ( $R_{2, SE}^{3 / 2}$ ). They found that the shape of the loop and the direction of its passage may differ between normal brain tissue and tumor tissue. Xu et al.¹² demonstrated by considering a tree model of microvasculature, that the direction of the loop is influenced mainly by the relative arterial and venous blood volume, as well as by the tracer bolus dispersion. These findings were confirmed from Hsu et al.¹³ by using a dual contrast agent injections approach. Emblem et al.¹⁴ recently detected an additional vasculature-depended effect on the MR signals: a temporal shift between GE and SE signals, which influences the formation of the vascular hysteresis loop (VHL). Depending on the hemodynamic properties of a tissue, the different sensitivities of the GE and SE signals to microscopic and macroscopic vessels result in a temporal shift between the respective MRI signal readouts. They demonstrated that patients with recurrent GBM, who responded to anti-angiogenic therapy with reduced tumor vessel calibers and improved hemodynamic efficiency and relative oxygen saturation, survived longer than patients without these responses. This promising technique, however, has yet not found its way into clinical routine. Reasons for this could be the lack of (i) clinically compatible MRI sequence protocols; (ii) imaging biomarkers considering the temporal shift phenomenon; and (iii) software tools for post-processing.

The aims of this study were threefold: (i) the development of an MRI sequence protocol for vascular architecture mapping (VAM) which is compatible to and practicable for clinical diagnosis; (ii) the introduction of new MR imaging biomarkers for the assessment of vascular pathologies in gliomas by development of novel evaluation strategies of VAM data; (iii) testing of the approach regarding feasibility and usefulness for the assessment of the microvascular architecture in a consecutive series of glioma patients in clinical routine.

Materials and methods

Patients

The study was approved by the Institutional Review Boards of the University of Erlangen and the University Clinic of St. Pölten. Analysis of patient data was performed in accordance with the Declaration of Helsinki. Written consent was obtained from all subjects. A total of 60 consecutive patients (24 women, 36 men; mean age ± standard deviation, 57.6 ± 15.1 years) with untreated or recurrent glioma were examined with the VAM technique as part of the routine MRI protocol between July and December 2015. Details about patient demographics and diagnosis are summarized in Table 1. From the total of 60 consecutive patients examined with the VAM technique in this study, four patients (6.7%) were excluded due to insufficient data quality: Severe patient motion between or during the MR scans were detected in three patients (patient no. 31, 42, and 54 in Table 1) with severe neurological deficits. In one patient (patient no. 39 in Table 1) with very low body weight, the bolus signal in the SE-EPI dynamic susceptibility contrast (DSC) perfusion MRI was too weak for evaluation. Calculation of VAM biomarker maps was successfully performed for the remaining 56 patients: 31 patients (14 women, 17 men; 62.0 ± 10.7 years) suffered from GBM WHO grade IV, 14 patients (7 women, 7 men; 48.8 ± 17.2 years) suffered from an anaplastic glioma WHO grade III (10 astrocytoma, 3 oligoastrocytoma, 1 oligodendroglioma), and 11 patients (2 women, 9 men; 52.8 ± 9.6 years) had a low-grade glioma WHO grade II (nine astrocytoma and two oligodendroglioma), respectively. All lesions were proven by histopathologic analysis.

Table 1.

Patient demographics and diagnosis.

No.	Age, sex	Diagnosis, localization	No.	Age, sex	Diagnosis, localization
1	62, m	GBM, r temporal	31	72, m	GBM, r temporal
2	70, f	GBM, l frontal	32	44, f	Astrocytoma II, l temporal
3	61,m	Astrocytoma III, r frontal	33	62, m	GBM, r frontal
4	54, m	GBM, r parietal	34	66, f	GBM, l frontal
5	51, m	GBM, r frontal	35	54, m	Astrocytoma II, l frontal
6	60, f	GBM, r parietal	36	45, m	GBM, l temporal
7	69, f	GBM, l parietal	37	63, m	GBM, l temporo-parietal
8	72, f	GBM, l temporal	38	44, m	Astrocytoma III, r frontal
9	57, m	GBM, l frontal	39	77, f	GBM, l parieto-occipital
10	34, m	Oligoastrocytoma III, r frontal	40	44, m	Astrocytoma II, r temporal
11	73, f	GBM, l frontal	41	80, f	GBM, l temporal
12	40, m	GBM, r temporal	42	54, m	GBM, r frontal
13	43, m	Astrocytoma II, l temporal	43	46, f	Astrocytoma III, r parietal
14	51, f	Astrocytoma II, r frontal	44	31, f	Astrocytoma III, r fronto-parietal
15	31, m	Astrocytoma III, l frontal	45	46, m	GBM, bifrontal
16	36, f	Oligoastrocytoma III, l temporo-parietal	46	51, m	GBM, r temporal
17	68, m	Astrocytoma II, l parieto-occipital	47	69, m	Astrocytoma II, r temporal
18	66, f	GBM, l fronto-temporal	48	55, m	Oligodendroglioma II, r parietal
19	54, f	GBM, l fronto-temporal	49	64, f	GBM, l frontal
20	63, f	GBM, r frontal	50	60, m	GBM, l frontal
21	48, m	GBM, l fronto-parietal	51	60, m	Oligodendroglioma III, r frontal
22	60, m	Oligodendroglioma II, l frontal	52	53, m	Oligoastrocytoma III, l frontal
23	67, f	Astrocytoma III, l frontal	53	76, f	Astrocytoma III, l temporal
24	79, f	Astrocytoma III, r temporal	54	71, m	GBM, l temporal
25	75, m	GBM, l temporal	55	43, m	Astrocytoma II, l frontal
26	69, m	GBM, r temporal	56	68, m	GBM, l frontal
27	49, f	GBM, l frontal	57	77, f	GBM, r fronto-parietal
28	34, f	Astrocytoma III, l fronto-temporal	58	61, m	GBM, r temporal
29	32, m	Astrocytoma III, r frontal	59	62, f	GBM, r temporal
30	51, m	Astrocytoma II, r temporal	60	83, m	GBM, l temporal

m: male; f: female; GBM: glioblastoma; r: right; l: left.

MR imaging

All MRI studies were performed on a 3 Tesla clinical scanner (Tim Trio, Siemens AG, Erlangen, Germany) equipped with a standard eight-channel head coil. The MRI protocol included an axial fluid-attenuated inversion-recovery (FLAIR) sequence (TR/TE/TI: 5000/460/1800 ms; in-plane resolution: 0.45 × 0.45 mm, slice thickness: 3 mm; 48 slices; 2 averages), a single-shot diffusion-weighted echo-planar imaging (DW-EPI) sequence (TR/TE: 5300/98 ms; in-plane resolution: 1.2 × 1.2 mm, slice thickness: 4 mm; 29 slices; 4 averages; parallel imaging using generalized autocalibrating partially parallel acquisition (GRAPPA) factor of 2, b-values of 0 and 1000 s/mm²), and a pre-contrast enhanced T1-weighted GE sequences (TR/TE: 250/2.8 ms; in-plane resolution: 0.5 × 0.5 mm, slice thickness: 4 mm; 29 slices; 2 averages).

For vessel architecture mapping, we used a dual contrast agent injections approach in order to obtain DSC perfusion MRI data using SE and GE EPI sequences with high signal-to-noise ratio (SNR), high spatial resolution, and coverage of the whole brain. However, patient motions between the two scans as well as differences in the time to first-pass peak may significantly affect the data evaluation.¹³ To minimize the probability of these sources of error, we used the following strategies:

To prevent patient motion, special caution was made to properly fix the head and provide clear and repeated patient instructions before and during the MRI examination. Due to the fact that patient motions during or between the SE- and GE-EPI DSC perfusion scans may significantly affect a pixel-based analysis, we used the approach described by Hsu et al.¹³ During data preprocessing, dynamic masks for the contours from the SE-EPI and GE-EPI perfusion raw data were created. By matching the contours, we checked for patient motion between and during the scans.¹³ A mismatch of more than one pixel was defined as “severe patient motion” and the patient was excluded from the analysis. This was necessary in three patients as mentioned above. For the remaining patients we did not observe significant motions, i.e. contour mismatch was not more than one pixel, and thus no image coregistration was applied during data analysis.

To prevent time differences of the first-pass peak signal between the contrast agent boluses, a peripheral pulse unit (PPU) was fitted to a finger of the patients in order to monitor heart rate and cardiac cycle. Special attention was paid to perform the two injections at the same heart rate and exactly at the same phase of the cardiac cycle (at PPU’s peak systole signal). Due to the fact that timing differences between the SE- and GE-EPI DSC perfusion scans may significantly affect a pixel-based analysis, we additionally checked for any error in injection timing. During data preprocessing the arterial input functions (AIFs) for the SE- and GE-EPI DSC perfusion data were verified for correct timing, i.e. were checked for differences in the time to peak (TTP) between the two AIFs. For all 56 patients, we did not observe differences between the SE- and GE-AIFs.

Geometric parameters (section orientation, angulation, field-of-view, etc.) and measurement parameters were chosen identical for the SE-EPI and GE-EPI DSC perfusion MRI sequences (TR: 1740 ms; in-plane resolution: 1.8 × 1.8 mm, note: SE-EPI were spatially interpolated to a 128 × 128 matrix to obtain this resolution; slice thickness: 4 mm; 29 slices; GRAPPA: 2; and 60 dynamic measurements) except for TE (33 ms for SE-EPI and 22 ms for GE-EPI). Both DSC perfusion MRI examinations were performed with administration of 0.1 mmol/kg-bodyweight gadoterate meglumine (Dotarem, Guerbet) at a rate of 4 ml/s using an MR-compatible injector (Spectris, Medrad) and with utmost caution for the injection time (see above). A 20 ml bolus of saline was injected subsequently at the same rate. The first DSC MRI was obtained by the SE-EPI technique since SE-EPI DSC perfusion MRI is less sensitive to CA leakage.^15,16 Post-contrast-enhanced T1-weighted images in axial, sagittal, and coronal orientation were performed between two DSC perfusion MRI acquisitions. The total acquisition time of these three sequences was 8 min. Thereafter, the second DSC MRI was obtained using the GE-EPI technique during a bolus injection under the same conditions as for the SE-EPI DSC perfusion MRI acquisition.

MRI data processing

VAM analysis was performed using custom-made in-house software developed using MatLab (The MathWorks, Inc., Natick, MA). SE- and GE-EPI DSC perfusion MRI were loaded into the software tool and firstly checked for motion artefacts during or between the acquisitions by using the contour approach described above and by Hsu et al.¹³ Next, AIFs for the SE- and GE-EPI DSC perfusion, MRI data were automatically identified, respectively, using a method described previously.¹⁷

ΔR(t)_2,SE and ΔR(t)_2,GE were calculated from the signals of the SE and GE perfusion data using the following equation

Δ R (t) 2, XE = - \frac{1}{TE XE} \cdot \ln \frac{S (t) XE}{S 0, XE}

(1)

where XE stands for SE or GE. TE is the echo time, S₀ is the baseline (prebolus) signal, and S(t) is the signal during the first bolus passage of the corresponding sequence. S₀ was determined as the mean of the signals from the 4th to the 15th dynamic volume.¹³ For S(t), the respective DSC perfusion MR data were corrected for remaining contrast agent extravasation^18,19 and a minimum search between the first and the second bolus passage (caused by recirculation) was performed as described previously¹² in order to separate the boluses. The truncated ΔR(t)_2, _SE and ΔR(t)_2, _GE curves of the first bolus were fitted to a previously described gamma-variate function²⁰ and used for calculation of the

Δ R 2, GE

versus

Δ R_{2, SE}^{3 / 2}

diagram, which we termed VHL, for each voxel in the brain covered by the volume of interest of the perfusion MRIs (Figure 1).

Figure 1.

Calculation of the vascular hysteresis loops. (a) Conventional MRIs of a 54-year-old female patient suffering from a recurrent glioma. (b) The principle of the VAM technique by using a dual contrast agent injections approach. (c) Calculation of $R 2, GE$ versus $R_{2, SE}^{3 / 2}$ diagrams, which we termed vascular hysteresis loops (VHLs).

MR imaging biomarkers of microvascular architecture

The VHL of each voxel was evaluated by the following five parameters: (i) microvessel radius (R_U) and (ii) density (N_U), which were adapted to the temporal shift phenomenon; as well as (iii) microvessel type indicator (MTI) as the signed area of the VHL; (iv) vascular-induced bolus peak time shift (VIPS) between the SE- and GE-EPI first passage bolus signal courses; and (v) the curvature (Curv) of the long-axis of the VHL.

Specifically, the new versions of the known parameters upper limit of microvessel radius (R_U) and upper limit of microvessel density (N_U)²¹ were adapted to the temporal shift phenomenon between the SE- and GE-EPI DSC perfusion MRI signal time courses. With this in mind, we calculated in a first step a maximum Q-index

Q \max = \max [Δ R 2, GE] / \max [Δ R_{2, SE}^{3 / 2}]

(2)

for each VHL, which is equivalent to the calculation of the slope of the diagonal of the enveloping rectangle of the VHL (Figure 2(a)). Q_max takes into account that the maxima of the

Δ R 2, GE

and

Δ R_{2, SE}^{3 / 2}

values, i.e. max[

Δ R 2, GE

] and max[

Δ R_{2, SE}^{3 / 2}

], are separated in the VHL when a vascular architecture-induced temporal shift of the first passage bolus signals between the dynamic SE- and GE-perfusion MRI exists. In the next step, R_U and N_U were calculated using the following equations²¹

R U = (\frac{CBV \cdot ADC \cdot b^{3}}{2 π \cdot Q_{\max}^{3}})^{1 / 2} and N u = \frac{Q \max}{b} \cdot (\frac{{CBV}^{2}}{4 π^{2} \cdot ADC \cdot {\bar{R}}^{4}})^{1 / 3}

(3)

where ADC is the apparent diffusion coefficient which was calculated from the DW-EPI MRI data. CBV is the absolute cerebral blood volume which was calculated from the GE-EPI DSC perfusion MRI data. As described by Jensen et al.,²¹

00 A 0 \bar{R}

is the mean vessel lumen radius (

\bar{R} \approx 3.0 μ m

) and b is a numerical constant (b = 1.6781). The MTI parameter was defined as the area of the VHL signed with the rotational direction of the VHL, i.e. a clockwise VHL-direction was identified with a plus-sign, and a counterclockwise VHL-direction was identified with a minus-sign. MTI values were calculated by subtraction of the areas under the ascending branch of the VHL minus the descending branch of the VHL (Figure 2(b)). The maximum

Δ R_{2, SE}^{3 / 2}

value was defined as reversal point. The VIPS parameter was introduced to determine the vascular architecture-depended time shift between the bolus peak signals of SE- and GE-EPI perfusion. The VIPS parameter quantifies the temporal shift of the peak signal of the first passage bolus signal between the SE- and GE-EPI perfusion which is induced by the vascular architecture.¹⁴ VIPS values were calculated by subtraction of the time to SE-EPI perfusion peak signal minus the time to GE-EPI perfusion peak signal. In other words, VIPS is the difference in TTP between the SE- and GE-EPI perfusion: VIPS = time(max[SE]) – time(max[GE]) = TTP(SE) – TTP(GE), i.e. a negative VIPS value meant that the SE-EPI perfusion signal peaks earlier than the GE-EPI perfusion signal and vice versa (Figure 2(c)).

Figure 2.

Calculation of the VAM biomarker. (a) Calculation of the upper limits for microvessel radius (R_U) and density (N_U) adapted to the temporal shift phenomenon. Calculation of (b) the microvessel type indicator (MTI), (c) vascular-induced peak shift (VIPS), and (d) the curvature (Curv) of the VHL. The color codes are depicted at the right image margins.

The Curv parameter was introduced to assess the shape of the VHL deviating from a symmetric ellipse which is associated with differences in the width of the GE- and SE-EPI perfusion boluses. The Curv values were calculated by fitting of the long-axis of the VHL with a quadratic polynomial followed by determination of the second derivative (Figure 2(d)).

Quantitative analysis

For quantitative analysis of the VAM approach in glioma patients, the mean MR imaging biomarker value within five regions of interest (ROI) was calculated. These ROIs were manually defined by a radiologist and an MR physicist in consensus for each patient based on features seen in the FLAIR and the contrast-enhanced T1-weighted images, respectively. ROIs were located in: (i) tumor core (for GBM: enhancing tumor region on CE T1-weighted MRI and for glioma WHO° III and II: hyperintense tumor region on FLAIR images); (ii) immediate peritumoral non-edematous tissue (predominantly gray matter); (iii) peritumoral edema (predominantly white matter); (iv) ipsilateral and (v) contralateral normal appearing brain tissue (iNAB and cNAB) of predominantly white matter. For GBM, we additionally defined ROIs in necrotic regions.

Statistics

Software (SPSS 14, IBM, Chicago, IL, USA) was used for statistical evaluation. Differences in VAM biomarkers between WHO grades were determined using analysis of variance (ANOVA). Additionally, a Dunnett-T3 test was used as post hoc procedure to be consistent with the assumption that homogeneity of variance was not met as well as for correction for multiple comparisons. Homogeneity of variance was tested using the Levene’s test. VAM biomarker values in lesions were compared with related values in peritumoral regions (non-edematous and edema) and NAB (ipsi- and contralateral) using a Wilcoxon signed-rank test. The Pearson’s correlation coefficient (R) was calculated as a measure of the strength of the linear relationship between two variables. P-values less than 0.05 were considered to indicate significance.

Results

Qualitative interpretation of VAM biomarker maps

Maps of microvessel radius R_U and microvessel density N_U showed heterogeneous patterns of mildly to severely increased levels in anaplastic glioma and GBM. Areas with severely increased N_U showed rather mildly increased to normal R_U and vice versa, i.e. maps of R_U and N_U provided complementary information about the heterogeneity of the microvasculature in high-grade glioma (Figures 2(a), 3(c) and (d), 4(c) and (d)). In LGG, however, R_U was found to be only moderately increased and N_U values were even slightly decreased (Figure 5(c) and (d)).

Maps of the MTI were interpreted based on the findings of Emblem et al.¹⁴ and Xu et al.¹² The VHL in a voxel of the VAM data set transverses in the counterclockwise direction if the vascular system in this voxel contains venule- and capillary-like vessel components,¹⁴ i.e. in case of a relatively higher venous blood volume,¹² whereas the VHL transverses in clockwise direction if vascular system consist of arterioles and capillaries,¹⁴ i.e. in case of relatively higher arterial blood volume.¹² Accordingly to our definition for the calculation of MTI values, a counterclockwise VHL-direction was associated with a negative signed VHL-area, i.e. a negative MTI, and vice-versa for the clockwise VHL-direction. In the MTI maps, negative MTI values were assigned to cool colors and positive to warm colors, respectively. Consequently, maps of MTI enabled differentiation between supplying arterial (areas with warm colors in Figures 2(b), 3(e), and 4(c)) and draining venous microvasculature (areas with cool colors in Figures 2(b), 3(e), and 4(c)) within high-grade glioma. More specifically, a voxel with high-arterial blood volume, for example shows a high positive MTI value (MTI > + 60) and an orange to red color in the MTI maps; a capillary voxel shows a very low MTI value around 0 ( − 5 < MTI < + 5) and a black color in the MTI maps; and a voxel with high venous blood volume shows a high, negative MTI value (MTI < − 60) and a purple color in the MTI maps. In line with the findings for LGG in the R_U and N_U maps, MTI maps showed no substantial changes compared to high-grade glioma and cNAB (Figure 5(e)).

Figure 3.

New VAM biomarkers for a 51-year-old male patient (no. 5 in Table 1) suffering from a recurrent glioblastoma. Conventional MR images of (a) contrast-enhanced T1-weighted and (b) FLAIR sequences. Maps of (c) the upper limit for microvessel radius (R_U) and (d) density (N_U), (e) the microvessel type indicator (MTI), the (f) vascular-induced peak shift (VIPS), (g) the curvature (Curv), and (h) the cerebral blood volume (CBV). The color codes are depicted at the right image margins.

Figure 4.

New VAM biomarkers for a 34-year-old male patient (no. 10 in Table 1) suffering from an oligoastrocytoma (WHO°III). Conventional MR images of (a) contrast-enhanced T1-weighted and (b) FLAIR sequences. Maps of (c) the upper limit for microvessel radius (R_U) and (d) density (N_U), (e) the microvessel type indicator (MTI), the (f) vascular-induced peak shift (VIPS), (g) the curvature (Curv), and (h) the cerebral blood volume (CBV). The color codes are depicted at the right image margins.

Figure 5.

New VAM biomarkers for a 68-year-old male patient (no. 17 in Table 1) suffering from an astrocytoma (WHO°II). Conventional MR images of (a) contrast-enhanced T1-weighted and (b) FLAIR sequences. Maps of (c) the upper limit for microvessel radius (R_U) and (d) density (N_U), (e) the microvessel type indicator (MTI), the (f) vascular-induced peak shift (VIPS), (g) the curvature (Curv), and (h) the cerebral blood volume (CBV). The color codes are depicted at the right image margins.

The interpretation of the maps of the vascular-induced peak shift (VIPS) was also based on the findings of Emblem et al.¹⁴ In areas with slower blood flow velocities, e.g. venules and capillaries, the SE-EPI perfusion signal peaks earlier than the GE-EPI perfusion signal. Based on our definition, this corresponded to a negative VIPS value and vice-versa for arteriole-dominated microvasculature with faster blood flow velocities. In the VIPS maps, similarly to MTI, negative VIPS values were assigned to cool colors and positive to warm colors, respectively. VIPS provided additional information to MTI about microvessel type especially at the tumor periphery. Presumably, VIPS is more sensitive to early neovascularization activity in areas in the infiltration zone characterized by a blood supply still within the normal volume range. VIPS, as a time-based parameter, is less sensitive to blood volume unlike MTI. The changes in MTI and VIPS in combination extended the contrast-enhancing tumor areas.

The Curv parameter showed almost exclusively negative values. Curv values in the tumor core of HGG (WHO grade III and IV) were similar compared to ipsi- and contralateral NAB. Curv was moderately decreased, i.e. showed moderately larger negative values, in peritumoral non-edematous tissue and was severely decreased, i.e. showed severely larger negative values, in peritumoral vasogenic edema. LGG, however, showed slightly decreased Curv values in the tumor core and in the peritumoral structures (non-edematous tissue and edema) when compared to ipsi- and contralateral NAB.

Quantitative evaluation of VAM biomarkers

ROI-based quantitative analysis revealed that R_U, N_U, and Curv were significantly increased (P-value range: < 0.001 to 0.013) and MTI and VIPS were significantly decreased (all P-values < 0.001) in the tumor core of GBM compared to peritumoral regions (non-edematous and edema) and ipsi- and contralateral NAB (red boxes in Figure 6). For anaplastic tumor core (HGGs), we found similar results for N_U and MTI only (P-value range: 0.001 to 0.019). R_U and Curv were significantly increased to peritumoral edema (P = 0.004 and 0.001) and ipsi- and contralateral NAB (P-value range: 0.001 to 0.038) but not to peritumoral non-edematous tissue (yellow boxes in Figure 6). VIPS was significantly decreased to ipsi- and contralateral NAB only (P = 0.008 and 0.019). Curv values in edema were significantly decreased compared to all other evaluated brain regions for GBM and anaplastic glioma WHO grade III (P = 0.001 to 0.01). For LGG, however, R_U was increased as well as N_U and Curv were decreased in the tumor core compared to ipsi- and contralateral NAB (green boxes in Figure 6), but these changes did not reach significance (P-value range: 0.05 to 0.131). All other parameter showed also no significant differences between brain regions in LGG. Interestingly, we found positive linear correlations between MTI and VIPS values in tumor core and the Curv values in peritumoral vasogenic edema of GBM (MTI: R = 0.756, P < 0.001; and VIPS: R = 0.608, P < 0.001) and anaplastic glioma WHO grade III (MTI: R = 0.737, P = 0.003; and VIPS: R = 0.558, P = 0.038). This might be interpreted as a relationship between the extent of neovascularization in the tumor and microvascular changes in edema, i.e. possibly severity of edema. LGGs showed no significant correlations.

Figure 6.

Box-whisker plots (whiskers: minimum, maximum; box: 25th to 75th percentile; line: mean value) for (top-down) the upper limit for microvessel radius (R_U) and density (N_U), the microvessel type indicator (MTI), the vascular-induced peak shift (VIPS), and the curvature (Curv) in (left-right) necrosis (only for glioblastoma), tumor core; immediate peritumoral non-edematous tissue, peritumoral edema; ipsilateral and contralateral normal appearing brain tissue (iNAB and cNAB) of predominantly white matter. Note: green boxes for low-grade glioma WHO grade II, yellow boxes for anaplastic glioma WHO grade III, and red boxes for glioblastoma WHO grade IV.

When compared the VAM biomarkers in tumor core between WHO grades, significant differences were found between LGG and GBM for R_U (P = 0.015), N_U, MTI, and VIPS (for all P < 0.001) as well as between glioma WHO grade II and III for N_U (P < 0.001). Curv in peritumoral edema was significantly higher in LGG compared to glioma WHO grade III (P = 0.032) and IV (P = 0.018). Table 2 gives an overview of the quantitative evaluation of the VAM biomarkers.

Table 2.

Overview of the quantitative evaluation of the VAM biomarkers.

	WHO°	Necrosis	Tumor	Peritumor	Edema	iNAB	cNAB
R_U	II	n.a.	50 ± 29	38 ± 10	32 ± 8	30 ± 9	30 ± 8
[µm]	III	n.a.	75 ± 42	53 ± 28	37 ± 16	30 ± 9	28 ± 5
	IV	31 ± 12	92 ± 59	48 ± 16	46 ± 18	31 ± 8	29 ± 7
N_U	II	n.a.	129 ± 55	178 ± 47	98 ± 26	171 ± 34	179 ± 26
[mm ⁻²]	III	n.a.	408 ± 190	212 ± 70	74 ± 18	173 ± 26	171 ± 22
	IV	57 ± 18	400 ± 167	230 ± 62	73 ± 21	194 ± 36	174 ± 28
MTI	II	n.a.	−1.2 ± 2.2	− 0.4 ± 1.3	− 0.2 ± 0.7	0.0 ± 1.5	0.1 ± 1.3
[s^−5/2]	III	n.a.	−33.5 ± 47.2	−2.9 ± 5.6	−0.3 ± 0.6	−0.1 ± 1.2	−0.2 ± 0.6
	IV	−2.8 ± 2.6	−74.8 ± 57.6	1.2 ± 5.7	-0.8 ± 2.1	0.1 ± 2.9	0.2 ± 1.2
VIPS	II	n.a.	−0.3 ± 0.9	−0.1 ± 0.6	−0.2 ± 0.8	0.0 ± 0.6	0.0 ± 0.6
[s]	III	n.a.	−1.9 ± 2.4	−0.6 ± 0.9	−0.9 ± 1.0	0.0 ± 0.5	−0.1 ± 0.3
	IV	−0.9 ± 0.5	−3.7 ± 2.2	0.1 ± 0.9	−0.8 ± 1.1	−0.1 ± 0.6	0.0 ± 0.5
Curv	II	n.a.	−28 ± 14	−31 ± 21	−32 ± 19	−24 ± 19	−20 ± 16
	III	n.a.	−14 ± 12	−27 ± 31	−72 ± 49	−27 ± 16	−28 ± 15
	IV	−17 ± 11	−16 ± 12	−32 ± 27	−65 ± 55	−21 ± 12	−22 ± 13

R_U: upper limit for microvessel radius; N_U: sdensity; MTI: microvessel type indicator; VIPS: vascular-induced peak shift; Curv: curvature; iNAB and cNAB: ipsilateral and contralateral normal appearing brain tissue.

Discussion

In this study, we developed an approach for assessment of microvascular architecture in glioma which is compatible to and applicable in clinical routine. By using two separate DSC perfusion MRI sequences (GE-EPI and SE-EPI) in combination with a dual contrast agent injections approach,¹³ we were able to ensure that the DSC perfusion MRI sequence (GE-EPI), which is essential for clinical brain MRI, was kept unchanged regarding spatial and temporal resolution. We additionally introduced new VAM biomarkers (MTI, VIPS, and Curv) and adapted known parameters (R_U and N_U)²¹ to the temporal shift phenomenon by development of a dedicated software tool. Finally, our VAM approach was tested in 60 consecutive glioma patients and an initial quantitative analysis was performed.

The compatibility and applicability of our approach were demonstrated since VAM data of good quality were obtained from over 94% of the consecutively recruited patient cohort. Furthermore, it was well accepted by both the staff and the patients. The use of the PPU, for example, did not increase the complexity of the procedure in clinical routine. The positioning of the PPU took only few seconds and had the side benefit that the patients kept still than without a PPU. The VAM biomarker maps provided a complementary insight into the heterogeneous structure of the microvasculature pathologies in glioma as consequence to angiogenetic activities initiated by the lesions. Microvessel radius R_U and density N_U were found to be inversely related, i.e. microvasculature in tumor areas showed vessels with small lumen and high density or with larger lumen and lower density. Statistical evaluation, however, did not reveal a significant negative correlation between R_U and N_U in the tumor core because the biomarker values for the tumor core were averaged over one large ROI during this initial quantitative analysis. Consequently, the heterogeneity of the vasculature within the tumor core as demonstrated in the R_U and N_U maps for high-grade glioma was not considered. MTI and VIPS maps color-coded with warm to cool colors enable differentiation of vessel type (arterioles/capillaries/venules) and might thus allow for detection of early neovascularization activity in the tumor periphery, respectively. MTI values are sensitive to both the blood volume and the vascular architecture-induced temporal shift between the SE- and GE-perfusion MRI signals. In lesions with marked hyperperfusion, the MTI maps are dominated by the strongly increased blood volume. VIPS, on the other hand, is a time-based parameter and is, therefore, less sensitive to blood volume. Based on our findings, we suppose that VIPS might be more sensitive to minor changes occurring during the early stages of neovascularization such as vascular co-option by infiltrative tumor cells.

The features of the Curv parameter, however, were surprising and unexpected. Curv values in peritumoral vasogenic edema of high-grade glioma (WHO grade III and IV) were significantly decreased compared to normal brain regions (iNAB and cNAB) as well as to tumor. Furthermore, the Curv values in peritumoral edema of HGGs correlated with both the MTI and the VIPS values in the tumor core of these lesions. These findings may be an indication for an interrelationship between the degree of neovascularization in the tumor, the severity of breakdown of the tight endothelial junctions that make up the BBB, and changes in tissue structure or microvasculature in peritumoral vasogenic edema.³ The Curv parameter might be a marker for the severity of vasogenic edema. However, there is no verified scientific basis for an interrelationship. Correlations with DTI data, histologic parameters, and Monte Carlo simulations for VHL formation could be informative.

Comparisons regarding WHO grade revealed significant differences in the mechanism of glioma-associated neovascularization between low- and high-grade glioma. While LGGs grow along preexisting vessels via vascular co-option, HGGs start to generate own tumor vessels in the process of growing. Vascular co-option is the first mechanism by which gliomas achieve their vasculature, e.g. in LGGs or in the infiltration zone of HGGs. Vascular co-option is followed by the development of new vessels in HGGs via four additional mechanisms of neovascularization which are known as angiogenesis, vasculogenesis, vascular mimicry, and transdifferentiation.¹ Our results are in line with the known fact that neovascularization is a key characteristic that distinguishes low- and HGGs, i.e. it is absent in LGGs, but prolific in HGGs, providing the tumor with ample blood supply for exponential growth.²²

This study has several limitations. First of all, the dual contrast agent injections approach in combination with two separate SE and GE DSC acquisitions is significantly more sensitive to both patient motions and variations in injection timing compared to the single GE-SE double echo EPI sequence, which was used in the most previous studies. Our developed strategies help to minimize the probability of these sources of error but cannot eliminate them entirely, however. A comparison and validation of our motion correction scheme using brain contours¹³ with available freeware like statistical parametric mapping (SPM)²³ are necessary. Furthermore, the injection of a double dose of contrast agent with two separate injections is required.

However, the application of separate SE and GE DSC sequences allowing the acquisition of VAM data with high SNR, high spatial resolution, and coverage of the whole brain because the entire EPI train is exclusively available for GE- or SE-EPI DSC data collection, respectively. For the combined, parallel GE-SE perfusion the EPI acquisition is shared between GE and SE perfusion data. This limits coverage and resolution or SNR. Both, a high in-plane resolution and the coverage of the whole brain are mandatory for routine GE-EPI perfusion. These features are especially important for detection of small or multi-centric lesions. The combined parallel GE-SE perfusion sequence used in previous studies,^14,24–26 however, does not meet these requirements due to insufficient spatial coverage (too few or too thick slices) or resolution (too small matrix size). This may necessitate performing the routine GE perfusion in addition or in a separate study, which consequently requires a further contrast agent injection. Additionally, most studies published so far used a double dose for the combined GE-SE perfusion^8,14,24,25 and/or applied a second injection for DCE perfusion^14,24,26 or correction of leakage effects.²⁵ This resulted in the application of up to a triple-dose. The efforts to upgrade the combined GE-SE perfusion sequence with the promising Simultaneous MultiSlice (SMS) technique for full coverage²⁷ may help to overcome these limitations. The interpolation to the target resolution, which potentially influences the magnitude of the measured susceptibility effect, and the short TE, which reduces contrast agent sensitivity and hence SNR, of the SE-EPI perfusion sequence are a further limitations. An increasing of the acquisition matrix would dramatically increase TR and consequently the temporal resolution. This is also true for a longer TE. The accuracy of the hysteresis quantification, however, is critically dependent on a high temporal resolution. Finally, validations of the biomarkers are required to enable rational interpretation.

In Conclusion, this study introduced new MRI biomarkers for clinical routine assessment of microvascular architecture in glioma. The VAM biomarker maps provided insights into the complexity and heterogeneity of vascular changes in brain tumors. However, investigations in more well-defined patient populations, histological validations, and more costly evaluations are required considering the heterogeneity of the vasculature especially in high-grade glioma.

Footnotes

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the German Research Foundation (Deutsche Forschungsgemeinschaft – DFG; Grant Numbers STA 1331/3-1 and DO 721/8-1) and by the ELAN program (Erlanger Leistungsbezogene Anschubfinanzierung und Nachwuchsförderung; Grant Number 14-05-21-1-Stadlbauer).

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Authors’ contributions

AS, GH, AD, MB, and KR conceived the study. AS and MZ involved in analysis of data. AS and MZ contributed to the development of the MRI biomarkers of microvascular architecture. AS, GH, AD, and KR involved interpretation of data. AS, GH, and SO contributed to acquisition of data. All authors reviewed the manuscript critically for intellectual content, and read and approved the final manuscript.

References

Hardee

Zagzag

. Mechanisms of glioma-associated neovascularization. Am J Pathol 2012; 181: 1126–1141.

Takano

Yamashita

Ohneda

. Molecular therapeutic targets for glioma angiogenesis. J Oncol 2010; 2010: 1–11.

Dvorak

. How tumors make bad blood vessels and stroma. Am J Pathol 2016; 162: 1747–1757.

Daumas-Duport

Tucker

Kolles

. Oligodendrogliomas. Part II: a new grading system based on morphological and imaging criteria. J Neurooncol 1997; 34: 61–78.

Burger

Vogel

Green

. Glioblastoma multiforme and anaplastic astrocytoma. Pathologic criteria and prognostic implications. Cancer 1985; 56: 1106–1111.

Boxerman

Hamberg

Rosen

. MR contrast due to intravascular magnetic susceptibility perturbations. Magn Reson Med 1995; 34: 555–566.

Dennie

Mandeville

Boxerman

. NMR imaging of changes in vascular morphology due to tumor angiogenesis. Magn Reson Med 1998; 40: 793–799.

Kiselev

Strecker

Ziyeh

. Vessel size imaging in humans. Magn Reson Med 2005; 53: 553–563.

Weisskoff

Zuo

Boxerman

. Microscopic susceptibility variation and transverse relaxation: theory and experiment. Magn Reson Med 1994; 31: 601–610.

10.

Fisel

Ackerman

Buxton

. MR contrast due to microscopically heterogeneous magnetic susceptibility: numerical simulations and applications to cerebral physiology. Magn Reson Med 1991; 17: 336–347.

11.

Kennan

Zhong

Gore

. Intravascular susceptibility contrast mechanisms in tissues. Magn Reson Med 1994; 31: 9–21.

12.

Kiselev

Möller

. Dynamic hysteresis between gradient echo and spin echo attenuations in dynamic susceptibility contrast imaging. Magn Reson Med 2013; 69: 981–991.

13.

Hsu

Y-Y

Yang

W-S

Lim

K-E

. Vessel size imaging using dual contrast agent injections. J Magn Reson Imaging 2009; 30: 1078–84.

14.

Emblem

Mouridsen

Bjornerud

. Vessel architectural imaging identifies cancer patient responders to anti-angiogenic therapy. Nat Med 2013; 19: 1178–1183.

15.

Essig

Wenz

Scholdei

. Dynamic susceptibility contrast-enhanced echo-planar imaging of cerebral gliomas. Effect of contrast medium extravasation. Acta Radiol 2002; 43: 354–359.

16.

Carroll

Horowitz

Shin

. Quantification of cerebral perfusion using the ‘bookend technique’: an evaluation in CNS tumors. Magn Reson Imaging 2008; 26: 1352–1359.

17.

Bjornerud

Emblem

. A fully automated method for quantitative cerebral hemodynamic analysis using DSC-MRI. J Cereb Blood Flow Metab 2010; 30: 1066–1078.

18.

Boxerman

Schmainda

Weisskoff

. Relative cerebral blood volume maps corrected for contrast agent extravasation significantly correlate with glioma tumor grade, whereas uncorrected maps do not. Am J Neuroradiol 2006; 27: 859–867.

19.

Boxerman

Prah

Paulson

. The role of preload and leakage correction in gadolinium-based cerebral blood volume estimation determined by comparison with MION as a criterion standard. Am J Neuroradiol 2012; 33: 1081–1087.

20.

Ducreux

Buvat

Meder

. Perfusion-weighted MR imaging studies in brain hypervascular diseases: comparison of arterial input function extractions for perfusion measurement. AJNR Am J Neuroradiol 2006; 27: 1059–1069.

21.

Jensen

Inglese

. Microvessel density estimation in the human brain by means of dynamic contrast-enhanced echo-planar imaging. Magn Reson Med 2006; 56: 1145–1150.

22.

Pisapia

Greenfield

. Malignant transformation in glioma steered by an angiogenic switch: defining a role for bone marrow-derived cells. Cureus 2016; 8: e471.

23.

Friston

Ashburner

. Statistical parametric mapping. Funct Neuroimaging Tech Found 1995; 15: 361–370.

24.

Sorensen

Batchelor

Zhang

. A ‘vascular normalization index’ as potential mechanistic biomarker to predict survival after a single dose of cediranib in recurrent glioblastoma patients. Cancer Res 2009; 69: 5296–5300.

25.

Schmainda

Rand

Joseph

. Characterization of a first-pass gradient-echo spin-echo method to predict brain tumor grade and angiogenesis. Am J Neuroradiol 2004; 25: 1524–1532.

26.

Batchelor

Gerstner

Emblem

. Improved tumor oxygenation and survival in glioblastoma patients who show increased blood perfusion after cediranib and chemoradiation. Proc Natl Acad Sci U S A 2013; 110: 19059–19064.

27.

Eichner

Jafari-Khouzani

Cauley

. Slice accelerated gradient-echo spin-echo dynamic susceptibility contrast imaging with blipped CAIPI for increased slice coverage. Magn Reson Med 2013; 72: 770–778.