Next Article in Journal
Research in the Field of Drug Design and Development
Next Article in Special Issue
The Development of a Smart Magnetic Resonance Imaging and Chemical Exchange Saturation Transfer Contrast Agent for the Imaging of Sulfatase Activity
Previous Article in Journal
Investigation of the Effect on the Albumin Binding Moiety for the Pharmacokinetic Properties of 68Ga-, 205/206Bi-, and 177Lu-Labeled NAPamide-Based Radiopharmaceuticals
Previous Article in Special Issue
Endogenous and Exogenous Antioxidants as Agents Preventing the Negative Effects of Contrast Media (Contrast-Induced Nephropathy)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Multi-Model High Spatial Resolution Method for Analysis of DCE MRI Data: Insights from Vestibular Schwannoma Responses to Antiangiogenic Therapy in Type II Neurofibromatosis

1
Division of Informatics, Imaging and Data Sciences, School of Health Sciences, Faculty of Biology, Medicine and Health, University of Manchester, Manchester M13 9PL, UK
2
Geoffrey Jefferson Brain Research Centre, University of Manchester, Manchester M13 9PL, UK
3
Department of Neurosurgery, Manchester Centre for Clinical Neurosciences, Salford Royal NHS Foundation Trust, Manchester Academic Health Science Centre, Manchester M13 9NT, UK
4
Division of Neuroscience and Experimental Psychology, School of Biological Sciences, Faculty of Biology, Medicine and Health, University of Manchester, Manchester M13 9PL, UK
5
Wolfson Molecular Imaging Centre, University of Manchester, 27 Palatine Road, Manchester M20 3LJ, UK
6
Department of Neuroradiology, Manchester Centre for Clinical Neurosciences, Salford Royal NHS Foundation Trust, Manchester Academic Health Science Centre, Manchester M13 9NT, UK
7
Division of Cardiovascular Sciences, School of Medical Sciences, Faculty of Biology Medicine and Health, University of Manchester, Manchester M13 9PL, UK
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Pharmaceuticals 2023, 16(9), 1282; https://doi.org/10.3390/ph16091282
Submission received: 3 July 2023 / Revised: 1 September 2023 / Accepted: 5 September 2023 / Published: 11 September 2023
(This article belongs to the Special Issue Next-Generation Contrast Agents for Medical Imaging)

Abstract

:
This study aimed to develop and evaluate a new DCE-MRI processing technique that combines LEGATOS, a dual-temporal resolution DCE-MRI technique, with multi-kinetic models. This technique enables high spatial resolution interrogation of flow and permeability effects, which is currently challenging to achieve. Twelve patients with neurofibromatosis type II-related vestibular schwannoma (20 tumours) undergoing bevacizumab therapy were imaged at 1.5 T both before and at 90 days following treatment. Using the new technique, whole-brain, high spatial resolution images of the contrast transfer coefficient (Ktrans), vascular fraction (vp), extravascular extracellular fraction (ve), capillary plasma flow (Fp), and the capillary permeability-surface area product (PS) could be obtained, and their predictive value was examined. Of the five microvascular parameters derived using the new method, baseline PS exhibited the strongest correlation with the baseline tumour volume (p = 0.03). Baseline ve showed the strongest correlation with the change in tumour volume, particularly the percentage tumour volume change at 90 days after treatment (p < 0.001), and PS demonstrated a larger reduction at 90 days after treatment (p = 0.0001) when compared to Ktrans or Fp alone. Both the capillary permeability-surface area product (PS) and the extravascular extracellular fraction (ve) significantly differentiated the ‘responder’ and ‘non-responder’ tumour groups at 90 days (p < 0.05 and p < 0.001, respectively). These results highlight that this novel DCE-MRI analysis approach can be used to evaluate tumour microvascular changes during treatment and the need for future larger clinical studies investigating its role in predicting antiangiogenic therapy response.

1. Introduction

Dynamic contrast-enhanced MRI (DCE-MRI) is a medical imaging technique that can be used to assess tissue perfusion and permeability. It requires an injection of a gadolinium-based-contrast agent (GBCA) into the blood, whilst acquiring a series of dynamic MRI images as the contrast agent passes through the tissue of interest. By analysing the dynamic changes in GBCA concentration in the tissue, pharmacokinetic modelling can then be used to derive various parameters such as the volume of distribution, the rate of blood flow, and the permeability of the tissue, in turn, providing valuable information on drug delivery and efficacy. This technique is particularly useful for evaluating the effectiveness of anti-cancer therapies, as it can provide insight into how efficiently a drug is penetrating the tumour microenvironment and whether the tumour is responding to treatment [1,2,3,4,5,6,7,8].
Traditionally, a key limitation in DCE-MRI was achieving both simultaneous high spatial (HS) and high temporal (HT) resolution during data acquisition. High temporal resolution, especially during the arterial phase of the GBCA time-course, is critical if accurate quantification of kinetic parameters such as the fractional plasma volume (vp), transfer constant (Ktrans), and the fractional volume of extravascular extracellular space (ve) is to be achieved [9,10,11,12,13,14]. Usage of DCE-MRI as a clinical tool, however, also requires high spatial resolution and whole-brain coverage, especially when lesions are small, heterogenous, or widespread throughout the brain. A DCE-MRI technique that can provide both high spatiotemporal resolution and whole-brain coverage for quantitative microvascular analysis is therefore highly desirable [15,16,17], and in earlier works, it was demonstrated that a novel dual-temporal resolution (DTR) DCE-MRI-based data construction technique termed the LEGATOS (level and rescale the gadolinium contrast concentration curves of high temporal to high spatial) method could be used to provide this [18,19].
Previous studies employing DTR DCE-MRI and the LEGATOS technique have used the extended Tofts model (ETM) [18,19] to derive tissue-validated high spatial resolution estimates of Ktrans, vp, and ve [20,21]. A key limitation of the ETM, however, is that the derived parameter Ktrans is a hybrid parameter reflecting both capillary plasma flow (Fp) and the capillary permeability-surface area product (PS) [22]. The ETM is also unable to accurately measure low-level permeability within normal-appearing brain tissue, which is of increasing importance in neuro-oncology and may help predict both prognosis and treatment response [23,24,25].
Recent studies have demonstrated that in addition to kinetic parameter estimation, absolute cerebral blood flow (CBF) estimates can also be derived from T1-weighted DCE-MRI through the use of the microsphere model and the maximum gradient approach [26,27]. In this technical proof of concept study, we demonstrate that incorporating ETM and the established ‘early time points’ method for absolute cerebral blood flow (ET-CBF) quantification [26] with the LEGATOS technique allows simultaneous evaluation of both tissue blood flow and permeability-based microvascular metrics [28]. In addition, we demonstrate that the integration of LEGATOS with the ETM and a previously described established hybrid first-pass Patlak plot (FP-PP) model [24] allows not only for high spatial resolution assessment of permeability effects within tumours but also high spatial resolution interrogation of low-level permeability effects within normal-appearing brain regions.
As a disease model to evaluate the presented multi-kinetic model technique, we retrospectively analysed previously acquired dual-injection DTR DCE-MRI data from a cohort of patients with neurofibromatosis type 2 (NF2) related vestibular schwannoma (VS) undergoing antiangiogenic (bevacizumab) therapy [4]. NF2 is a dominantly inherited tumour predisposition syndrome, the hallmark of which is the development of VS arising from the vestibulocochlear nerve bilaterally [29,30,31]. The management of these sometimes rapidly progressive VS can present a significant treatment challenge [32,33], and affected patients often develop, alongside bilateral VS, other tumours such as multiple meningiomas [34]. Trials of the anti-vascular endothelial growth factor (VEGF) antibody bevacizumab (Avastin©) in NF2 patients have demonstrated a reduction of VS growth rate [35] and regression in some cases [36], but there is variation in treatment results, with reported tumour response rates of 40–60% [35,36,37]. Previous studies have demonstrated that microvascular parameters derived from DCE-MRI and the ETM have possible treatment predictive value [4,35], and within this study, we sought to evaluate the applicability of the presented multi-model DCE-MRI analysis approach to this disease model. In particular, we sought to evaluate if this novel multi-model high spatial DCE-MRI analysis technique could be used to differentiate treatment-related changes in microvascular flow and permeability within the tumour microenvironment and identify parameters that may allow for better treatment response prediction in future clinical trials.

1.1. Multi-Model LEGATOS Analysis Theory

1.1.1. Construction of a 4D High Spatiotemporal Resolution GBCA Concentration-Volume Using the LEGATOS Method

The use of DTR DCE-MRI and the LEGATOS method for deriving whole-brain, high spatial resolution microvascular parameters have been previously described [20,21]. Low-dose high temporal resolution (LDHT) data is initially acquired during the early (arterial) phase of the contrast agent bolus arrival when rapid changes in signal intensity are expected. A second full-dose high spatial resolution (FDHS) dataset is then acquired during the later phases of the bolus arrival when the signal changes are slower and steadier. By combining these two sets of data, the LEGATOS method can improve the accuracy of quantitative analysis of tissue perfusion dynamics, permitting high spatial resolution assessment of tissue microvascular parameters.
The LEGATOS method consists of two key steps. In the first step, a concatenated DTR 4D GBCA concentration volume is constructed containing a high temporal resolution ‘arterial’ phase followed by a later low temporal, high spatial resolution ‘parenchymal’ phase. The 4D LDHT DCE images are first co-registered to a baseline image frame of the FDHS DCE series to obtain a 4D LDHTaligned volume. The signal intensity-time curves from this 4D LDHTaligned and the 4D HS dynamic image volumes are then converted to GBCA concentration-time curves, before combining the HT and HS resolution series to construct a 4D GBCA concentration volume. In the second step, the LDHTaligned arterial phase of each voxel concentration curve is re-scaled using a derived voxel-wise calibration ratio (ratiocalib) to increase the spatial resolution of the derived kinetic parameter maps. Further details on the LEGATOS reconstruction method and derivation of the calibration ratio are provided in prior publications [20,21].

1.1.2. Use of LEGATOS with the ETM and ET-CBF Model to Derive High Spatial Resolution Ktrans, vp, ve, and Fp Maps

The most widely used pharmacokinetic model in brain DCE-MRI is the ETM [18,19], which can be defined as follows:
C t ( t ) = v p C p ( t ) + K trans   0 t C p ( τ ) exp ( K trans v e ( t τ ) ) d τ ,
where Ct is the tissue contrast agent concentration and Cp is the GBCA concentration in the vascular plasma space. In the convolution integral, t is considered a constant and τ is the variable. Through use of the LEGATOS technique in combination with the ETM high spatial resolution, estimates of Ktrans, vp, and ve could be provided.
ET-CBF is a methodology for absolute CBF estimation, which uses data derived from a LDHT T1W DCE-MRI acquisition, and is based on the microsphere model and ET (early time points before the contrast agent has left the tissue) strategy [26]. The microsphere model can be written mathematically as follows:
C t ( t ) = f 0 t C b ( t ) d t with t ETW ,
where the tissue GBCA concentration, Ct(t), is equal to the amount of GBCA delivered to 1 mL of tissue by time t; f is the absolute CBF (mL min−1 mL−1); and Cb(t’) is the arterial blood concentration at time t′. ETW is the early time window, i.e., the time window that meets the microsphere prerequisite. Using the microsphere model and Equation (2), low spatial resolution early time points absolute cerebral blood flow (CBFET) maps can be derived from the arterial concentration-time curve of a low dose high temporal DCE-MRI acquisition.
High spatial resolution CBFET (CBFET-HS) maps can also be derived from DTR DCE-MRI data through the use of the microsphere model with LEGATOS. CBFET-HS maps can be obtained either directly from the concatenated 4D high spatiotemporal resolution GBCA concentration volume generated using the LEGATOS approach described above or by using the calibration ratio maps generated as part of the LEGATOS procedure to rescale (see Equation (3)) the low spatial resolution CBFET (CBFET-HT) derived from a 4D LDHTaligned volume (see Appendix A for details):
CBF ET - HS = CBF ET - HT ratio calib .

1.1.3. Theoretical Derivation of the Capillary Permeability-Surface Area Product (PS) from Derived Ktrans and CBFET Values

Estimation of PS is based on the following relationship [38]:
K trans = F p PS / ( F p + PS ) ,
where Fp is plasma flow. Fp relates to blood flow by the following:
F p = CBF ET ( 1 Hct ) ,
where Hct is the blood haematocrit and CBFET is the absolute cerebral blood flow estimated using the early time points method.
From Equation (4), we have the following:
PS = K trans / ( 1 K trans / F p ) .
Therefore, PS can be calculated from the known Ktrans and Fp using Equation (6). Using Equation (4), Ktrans is dominated by the smallest of PS and Fp [38].
K trans / PS = F p / ( F p + PS ) = R Fp .
Ktrans/PS is the ratio of plasma flow to the sum of Fp and PS or RFp, ranging from 0 to 1:
when Fp >> PS, RFp ≈ 1;
when Fp = PS, RFp = 0.5;
when Fp << PS, RFp becomes a small positive value (≈ Fp/PS).
The RFp maps can be useful in identifying areas of high perfusion but inadequate permeability (High RFp) and high permeability but inadequate perfusion (Low RFp), respectively, thereby providing insight into tumour vascular heterogeneity [39]. High spatial resolution PS maps can also be obtained from estimates of Ktrans and Fp derived using the LEGATOSETM and CBFET-HS techniques described above, respectively.

1.1.4. Use of LEGATOS-Patlak for Measuring Ktrans within Normal-Appearing Brain Regions

The Patlak model describes a highly perfused two-compartment tissue assuming unidirectional transport from the plasma into the extravascular extracellular space [6,40,41,42,43,44]. The GBCA concentration in tissue (Ct) is given by the following:
C t ( t ) = v p C p ( t ) + K trans 0 t C p ( τ ) d τ .
Dividing both sides of Equation (8) with Cp(t), one obtains the Patlak plot (Equation (9)) where the slope represents Ktrans and the intercept represents vp. The abscissa has the units of time, but this is not laboratory time; it is concentration-stretched time and will be referred to hereafter as tstretch [41].
C t ( t ) C p ( t ) = v p + K trans   0   t C p ( τ ) d τ C p ( t ) .
One issue in the conventional Patlak plot method is that experimental errors are distorted when a non-linear model is transformed to a linear one [45]. This can introduce non-uniform (distorted) noise into the fitting procedure and reduce the reliability of the estimates of Ktrans and vp. A previous study proposed a “hybrid FP-PP” method that combines a first-pass analytical approach [46,47] with the Patlak plot to address these limitations [24]. Both computer simulation and in vivo studies demonstrated improved reliability in vp and Ktrans estimates with the hybrid method [24], and further details on this hybrid FP-PP method can be found in the included reference [24].
The “hybrid FP-PP” model, however, requires high temporal resolution DCE-MRI data at the expense of limited spatial resolution and/or volume coverage. High spatial resolution assessment of low-level permeability in normal-appearing brain regions can be achieved by fitting the LEGATOS-generated 4D high spatiotemporal resolution GBCA concentration volumes with the hybrid FP-PP model. vp estimates obtained from LEGATOS-ETM (instead of from the FP analysis) can also be used as the known vp in the modified Patlak plot linear regression analysis [24] to produce the only free-fitting parameter, Ktrans.
A flowchart in Figure 1 illustrates the key steps in the multi-model LEGATOS integrated pharmacokinetic analysis for the derivation of high spatial resolution microvascular biomarkers within both tumour and normal-appearing brain regions (see Appendix B for a list of abbreviations used in the dual-injection dual-temporal resolution DCE-MRI pharmacokinetic analysis).

2. Results

2.1. High Spatial Multi-Model Assessment of Perfusion and Permeability Parameters within Both Tumour and Normal-Appearing Brain

Representative pre-treatment early time points absolute CBFET maps for a patient with NF2-related VS and meningioma are shown in Figure 2 and demonstrate how high spatial resolution absolute CBFET-HS maps can be achieved using the early time points model with LEGATOS. Figure 2A shows representative pre-treatment images from an NF2 patient with a small convexity meningioma (white arrow). This meningioma was not clearly defined within the native ET-CBF map (left) due to the low spatial resolution and partial volume effects, but it is much better demonstrated within the reconstructed HS-CBF map and the LEGATOS-ETM derived high spatial resolution vp map (right column).
Figure 2B shows representative images from a patient with bilateral NF2-related VS. The CBFET-HS map offered superior visualization of intratumoural heterogeneity within the left-sided VS compared to the CBFET-HT map and closely resembled the LEGATOS-ETM derived high spatial resolution tumour vp map as expected. Compared to the unreconstructed, low spatial resolution CBFET images, the CBFET-HS maps offered superior visualization of both the smaller right-sided VS and intratumoural heterogeneity in blood flow. The CBFET-HS maps spatially corresponded with the pattern of vascularisation seen in the LEGATOS-ETM-derived high spatial resolution vp maps.
Representative high spatial resolution Ktrans maps provided by the LEGATOS method, using either the FP-PP approach or ETM are shown in Figure 3A. Both maps showed comparably high Ktrans values within the pre-treatment meningioma, with median tumour Ktrans values of 0.025 and 0.017 min−1 for ETM and FP-PP models, respectively. The Patlak analysis underestimated Ktrans of the meningioma because the backflux of GBCA from the extravascular extracellular space to the plasma was too large to ignore in the tumour tissue. Compared to the LEGATOSETM method, however, integration of the hybrid FP-PP model with LEGATOS (LEGATOSFP-PP) permitted high spatial resolution assessment of low-level Ktrans in the normal-appearing brain.
Descriptive Ktrans statistics obtained from NAGM/NAWM in the 12 patients are shown in Table 1. Segmented NAGM displayed non-significantly higher mean (p = 0.30) and significantly higher median (p = 0.03) and max (p = 0.002) Ktrans values compared to segmented NAWM (paired t-test). Pearson correlation analysis (Figure 3B) shows that in this NF2 patient cohort, significant positive correlations between tumour volume and NAGM/NAWM Ktrans values were observed. Mean Ktrans measured from both the NAGM or NAWM segments showed a significant correlation with both average VS volume size (p ≤ 0.02) and total VS volume (p ≤ 0.01). Kendall rank correlation test also shows significant positive monotonic correlations of NAGM Ktrans values with both average and total VS volume size (p = 0.04 and p = 0.05, respectively).

2.2. High Spatial Evaluation of Changes in Tumour Microvascular Parameters during Antiangiogenic Therapy

Representative pre- and post-bevacizumab high spatial resolution microvascular parameter maps (Ktrans, PS, Fp, and RFp) are shown in Figure 4A for an NF2 patient with a large VS. Pre-treatment, there was intratumoural heterogeneity in the derived perfusion and permeability metrics with distinct tumour regions showing either high Fp or high Ktrans and PS, respectively. Such heterogeneity was also evident in derived RFp maps (ratio of plasma flow to the sum of Fp and PS, Fp/(Fp + PS)), with regions showing high RFp (permeability limited) and low RFp (perfusion limited), respectively. At 90 days post-bevacizumab therapy, reductions in VS volume are also seen along with substantial reductions in both Ktrans and PS. There is also an observed corresponding increase in RFp, with an increase in the tumour subregion displaying high perfusion but inadequate permeability.
In Figure 4B, the histogram distributions before and after 90 days of treatment, calculated for the same tumour, are shown. After 90 days of bevacizumab treatment, Ktrans and PS histograms showed a shift towards lower values, with a median decrease of 21% in Ktrans and a 33% decrease in PS. RFp histograms conversely showed a shift towards higher values, with a median increase of 10%. The overall tumour volume was reduced by 28%, and tumour volume loss was principally seen in voxels with Fp in the range of 0.2–0.8 min−1 (see arrowed red lines on Fp histogram), with 2638 out of 3005 and 1588 out of 2163 tumour voxels being in this range at day 0 and day 90 post-treatment, respectively.
Post-treatment changes in each high spatial resolution microvascular parameter across the twelve imaged patients are shown in Table 2. Responding tumours displayed significant post-treatment reductions in Ktrans (p ≤ 0.001) and PS (p = 0.0002) and significant increases in RFp (p = 0.004). Post-treatment decreases in Fp (p = 0.08), vp (p = 0.07), and ve (p = 0.56) were also seen but these did not reach statistical significance. Non-responding VS also displayed a significant increase in RFp at day 90 (p = 0.04). Decreases in mean Ktrans and PS were not significant (p > 0.05); however, in contrast to responding tumours, non-responding VS displayed significant increases in ve (p = 0.01).

2.3. Differences between Responding and Non-Responding Tumours in Baseline (Pre-Treatment) High Spatial Resolution Microvascular Parameters

Differences in pre-treatment high spatial microvascular parameters between responding and non-responding VS are shown in Figure 5. Responding tumours showed significantly higher pre-treatment PS (p = 0.045) and ve (p < 0.001) values. Responding tumours also displayed higher Ktrans than non-responding VS but these differences (p = 0.07) did not reach statistical significance. Of the five parameters, baseline ve correlated most strongly with percentage volume change at day 90 (p < 0.001).
Table 3 shows associations between the estimated pre-treatment tumour median Ktrans, Fp, PS, vp, and ve and tumour volumetric parameters (the baseline tumour volume, tumour volume change, or percentage tumour volume change at day 90) assessed by linear regression analysis. Both PS and Ktrans (a hybrid flow-permeability parameter) showed a significant correlation with baseline tumour volume (p = 0.03 and p = 0.05, respectively) and a non-significant trend of correlation with tumour volume change at day 90 (p = 0.06 and p = 0.12, respectively). On the other hand, ve, an oedema-associated parameter [48], showed a significant correlation with tumour volume change and percentage tumour volume change at day 90 (p = 0.04 and p = 0.0002, respectively) and a non-significant trend of correlation with baseline tumour volume (p = 0.08). Fp and vp did not show correlations with any of the three tumour volumetric parameters.
Univariate logistic regression analysis revealed that, of the five microvascular parameters derived using the new method, ve estimates (AUC = 0.896, p = 0.02, Table 4) showed the greatest ability in the prediction of 90-day response in NF2-related vestibular schwannoma. vp was the only variable with a p value greater than 0.50. In multivariate analysis, the backward selection procedure starts with the five potential predictors and eliminates the least significant variable at each step. The model with three variables (ve, PS, and Fp) demonstrated a higher AUC, sensitivity, and specificity than the models with fewer variables.

3. Discussion

We have described a novel processing technique for DTR DCE-MRI, which, using the LEGATOS method with multi-model kinetic analysis, allows for the derivation of accurate, high spatial resolution, whole-brain coverage kinetic parameter maps. In contrast to the use of the extended Tofts model alone, this new approach allows for the separate derivation of tumoural blood flow and permeability, as well as a high spatial resolution assessment of low-level permeability in a normal-appearing brain. We evaluated our new analysis method through an in vivo study of patients with neurofibromatosis type 2 (NF2)-related vestibular schwannoma undergoing antiangiogenic (bevacizumab) therapy.
In a previous longitudinal DCE-MRI study of NF2-related VS, it was demonstrated that pre-treatment, low spatial resolution tumoural Ktrans estimates had the potential ability to predict later VS volume response to anti-VEGF therapy [4]. As a parameter, however, Ktrans reflects regional capillary blood flow, capillary endothelial permeability, and the surface area of the capillary endothelial membrane [22], and the extent to which changes in each of these physiological variables changed during treatment was unknown. From earlier in vivo studies, the primary action of VEGF inhibitors, such as bevacizumab, was thought to be the rapid reduction in capillary endothelial membrane permeability [49,50], and it was hypothesized that the reduction in Ktrans observed following treatment might be indicative of this mechanism. However, an alternative interpretation was that changes in tumour blood flow could also affect measured Ktrans, and to address this question, our study utilized a high spatial resolution multi-kinetic model analysis technique that permits simultaneous estimation of both Fp and PS within the tumour microvasculature.
Our results demonstrated that the capillary permeability-surface area product (PS) showed more pronounced treatment-related changes than Ktrans, with responding tumours showing both significantly higher baseline PS and significant reductions in PS at 90 days post-treatment compared to non-responding VS. On the other hand, our findings revealed that Fp alone was not predictive of tumour volume response to bevacizumab treatment. There was no difference in pre-treatment Fp between responding and non-responding tumour groups, and in neither group was there a significant change in Fp during treatment observed. Instead, across all VS studied, there was an increase in the relative scale of perfusion to permeability (higher RFp) during treatment. Overall, these imaging observations suggest that the changes observed in Ktrans following bevacizumab treatment are primarily driven by a reduction in capillary endothelial membrane permeability, rather than hypothesized changes in tumour blood flow. Such a finding is in line with previous in vivo murine models evaluating anti-VEGF therapy response [4,51,52], which have shown early (< 24 h) vascular normalization with reductions in both vascular permeability and the surface area of vessels; this is to our knowledge the first time that microvascular flow and permeability changes following anti-VEGF therapy in NF2-related VS have been differentiated through an in vivo human imaging analysis.
This study demonstrated that, of the five microvascular parameters investigated (Ktrans, PS, Fp, vp, and ve), baseline ve showed the most significant association with tumour volume response to bevacizumab treatment. These findings are supported by previous studies on the predictive value of an apparent diffusional coefficient (ADC) and a native longitudinal relaxation rate (R1N) in relation to tumour volume response to the VEGF inhibition [4,35,53]. High values of ADC and low values of R1N are in keeping with a larger extravascular extracellular space, increased levels of interstitial free fluid, and likely higher capillary permeability. Within our study, we found that baseline ve had a strong correlation with percent volume reduction at day 90, which is consistent with Plotkin et al., who found that high baseline ADC values were associated with tumour volume reduction at 3 months [35]. Although high baseline ve was a strong correlate of tumour shrinkage, tumour median values of ve showed only a non-significant percentage decrease (−3.1% ± 21.4%, p = 0.56) in the ‘responder’ group and a significant percentage increase (19.3% ± 17.0%, p = 0.01) in the ‘non-responder’ group at day 90. A multivariate model combining pre-treatment ve, PS, and Fp demonstrated high sensitivity and specificity for the prediction of volumetric response at 90 days, and these results, whilst preliminary, suggest that initial volumetric response to anti-VEGF therapy appears more likely in VSs that show both high levels of permeability and intratumoural oedema [4,54]. Future larger studies are, however, required to better understand the role of ve and PS in anti-VEGF therapy prediction.
Alongside demonstrating that the presented multi-model analysis approach permits differentiation of flow and permeability effects within the tumour microenvironment, we demonstrated that integration of the hybrid FP-PP model with LEGATOS also permits high spatial resolution assessment of low-level permeability within the normal-appearing brain. In keeping with previous literature demonstrating higher blood volume, higher vascular surface area, and higher transfer constant estimates in grey compared to white matter [15,55,56], our study data demonstrated higher Ktrans estimates within NAGM compared to NAWM. In this NF2 patient cohort, a significant positive correlation between total VS disease burden and Ktrans values in NAGM/NAWM was observed. In patients with VS, there is emerging evidence that there can be effects on the brain remote to the tumour, with previous diffusion and functional MRI studies demonstrating widespread changes in activity networks, grey matter volume, and white matter fibre integrity in auditory and non-auditory regions in patients with these tumours [57,58,59,60,61]. To date, however, changes in DCE-MRI parameters within the normal-appearing brain of patients with NF2 have not been evaluated and further large, detailed studies are required to better understand the pathophysiology of the observed Ktrans changes and their relationship to the tumour burden in these patients.
A key limitation of the presented study is that the number of patient participants was small. The primary aim of this study was to develop and evaluate the in vivo applicability of a high spatial resolution multi-kinetic model DCE-MRI approach rather than robustly evaluate predictors of NF2-related VS antiangiogenic response through a large clinical trial. Patients with NF2-related VS undergoing bevacizumab therapy were used as a disease model to test this new analysis approach due to both the previously reported predictive potential of DCE-MRI in this tumour group and the purported vascular normalization mechanism of anti-VEGF therapy that has been demonstrated in pre-clinical studies, which could be interrogated through our DCE-MRI approach. Larger future prospective studies incorporating this multi-model DCE-MRI analysis should be performed to better understand the predictive potential of each derived parameter. A second limitation is that the DTR DCE-MRI protocol and analysis technique used in the study has relatively long scan (> 10 min) and data processing times. Using the multi-model LEGATOS integrated pharmacokinetic analysis technique described here, the total computer processing time (excluding image registration time) to obtain the comprehensive panel of microvascular biomarkers with high spatial resolution and whole-brain coverage was approximately six hours. The potential impact of MRI data noise on the reliability of derived microvascular parameters also potentially remains a challenge. Further research and optimization are required to overcome these limitations and establish the technique’s practical value in routine clinical practice.

4. Materials and Methods

4.1. Patients

As a disease model to evaluate this new multi-model DCE-MRI analysis technique, previously acquired dual-injection, DTR DCE-MRI data in twelve consecutive patients with NF2-related VS were analysed. All patients were undergoing antiangiogenic therapy with the anti-VEGF antibody, bevacizumab (Avastin©), and across all patients there were twenty VSs, with four patients having undergone previous surgical resection of a VS. All patients had been recruited and scanned as part of an earlier described study investigating MR imaging predictors of antiangiogenic response [4]. Patients with standard contraindications to MR imaging or contrast agent administration were excluded. Only patients with proven NF2 and at least one VS demonstrating a high growth rate (defined as an annual increase of 4 mm or greater in the maximal transverse diameter) were included. Detailed study inclusion/exclusion criteria along with bevacizumab dose information are detailed in prior publications [4]. Ethical approval was in place for this study (REC reference 13/NW/0247) and all recruited participants had previously consented to later analysis of their MRI data.
Patients underwent DCE-MRI on 2 occasions: pre-treatment (day 0) and 3 months (day 90) following treatment. Tumours (VSs) were classified according to volumetric response/reduction over the 90-day bevacizumab treatment period with response defined as a volume reduction exceeding 0.125 cm3 and/or a relative volume decrease exceeding 5%. Further details on how volume response was defined are provided in the included reference [4].

4.2. MRI Data Acquisition

All patients were imaged at 1.5 T (Philips Achieva, Best, The Netherlands) and dual-injection DTR DCE-MRI data were collected using a 3D spoiled gradient-recalled echo sequence as described previously [62]. Using a power injector, a macrocyclic GBCA (gadoterate meglumine; Dotarem, Guerbet S.A., Villepinte, France) was administered at a rate of 3 mL/s as an intravenous bolus, followed by a chaser of 20 mL/s of 0.9% saline at the same rate. A high temporal (Δt = 1.0 s) but low spatial (voxel size = 2.5 × 2.5 × 6.35 mm3) resolution DCE-MRI series (NFrame = 300) with a low-dose fixed-volume (3 mL) of GBCA was performed for a total scan duration of 5.1 min. Subsequently, a high spatial (voxel size = 1 × 1 × 2 mm3) but low temporal (Δt = 10.7 s) resolution DCE-MRI series with a full dose (0.2 mL/kg·weight—3 mL dose of pre-bolus) of GBCA was performed for a total scan duration of 10.6 min (NFrame = 60). Variable flip angle (a = 2°, 8°, 15°, and 20°) acquisitions were performed prior to the LDHT and FDHS DCE series for baseline longitudinal relaxation rate estimation.

4.3. Image Processing

A 4D high spatiotemporal resolution GBCA concentration volume was generated from DTR DCE-MRI using the LEGATOS technique. The LEGATOS reconstructed concentration volume was subsequently fitted with the ETM, ET-CBF, and hybrid FP-PP models to derive high spatial resolution estimates of Ktrans, vp, ve, CBFET, Fp, RFp, and PS, and to detect low-level permeability (Ktrans) in normal-appearing brain regions.
As outlined above, low spatial resolution CBFET maps were derived from the arterial phase of the 4D LDHTaligned volume (which is the arterial phase of the concatenated 4D DTR GBCA concentration volume generated in key step I of the LEGATOS approach). High spatial resolution CBFET maps were obtained by rescaling the low spatial resolution CBFET with the calibration ratio maps generated as part of the LEGATOS procedure. CBFET-HS and CBFET-HT maps could therefore be compared as part of this study.
Quantifying DCE-MRI-derived kinetic parameters through modelling requires a suitable vascular input function (VIF), and previous studies have demonstrated that the usage of the superior sagittal sinus (SSS) as a surrogate input function provides a good approximation [63,64]. This VIF measurement method uses a semi-automatic extraction method to identify voxels within the SSS that display maximum enhancement during the first pass of the GBCA bolus, as has been previously described [65]. For the LEGATOS reconstruction, a combined vascular input function is constructed through concatenation of the GBCA concentration-time curve, Cp (t), from the LDHT series with the later parenchymal phase of the dose-calibrated GBCA concentration-time curve from the FDHS-DCE series [4,6]. The amplitude of the Cp (t) from the FDHS series is scaled down to match the LDHT-derived Cp (t) using the dose calibration ratio prior to concatenation. Further details on the derivation of the VIF for this analysis can be found in the included reference [20].
SPM (statistical parametric mapping) [66] was used for all image co-registration and for segmentation of the 3D T1W MRI data into normal-appearing grey matter (NAGM), white matter (NAWM), and cerebrospinal fluid (CSF) probability maps. After re-alignment/re-slicing of each map to a baseline FDHS DCE frame, masks of NAGM and NAWM were generated using a probability cut-off of 0.95. Only supratentorial NAGM/NAWM regions were included in each mask for quantitative analysis of Ktrans within normal-appearing brain regions.

4.4. Statistical Analysis

High spatial resolution tumour microvascular parameters derived using the above multi-model (ETM, ET-CBF, FP-PP) approach were compared both before and 90 days after antiangiogenic therapy. Median tumour values for Ktrans, PS, Fp, RFp, ve, and vp were calculated for each visit for the 20 VSs across the 12 NF2 patients. The group mean and standard deviation of these tumour median values were compared across the two visits using a paired t-test, with responding (N = 12) and non-responding (N = 8) VSs analysed separately. Histograms of the fitted voxelwise microvascular parameters (Ktrans, PS, Fp, RFp, ve, and vp) before and 90 days after treatment were also compared.
Estimates of Ktrans within NAGM and NAWM across the twelve patients were compared using descriptive statistics and paired t-tests. Estimates of mean Ktrans measured from the NAGM/NAWM segmentations were also compared with VS volume using scatterplot analysis and reported as Pearson’s and Kendall’s correlation coefficients (r, rs, and τ, respectively). The tumour volume used in the correlation analysis was defined as either (1) the average size across both VSs (or the size of the single VS in patients previously having undergone a VS resection) or (2) the cumulative tumour volume across both/single VSs.
The hypothesis that pre-treatment tumour median values of Ktrans, PS, Fp, vp, and ve do not differ between responding and non-responding VS was tested using the unpaired Student’s t-test. Binary logistic regression (S-Plus, version 6.1; Insightful, Seattle, WA, USA) was also performed to assess the ability of these imaging biomarkers to predict tumour response at day 90. The predictive performance of the model was assessed by receiver operator characteristics (ROC) analysis and calculation of the area under the curve (AUC). Univariate and multivariate analyses were performed. A backward elimination variable selection method was used in the multivariate analysis to obtain the best model.
Linear regression analysis was also performed to assess the relationship between imaging parameters at baseline and tumour volumetric parameters (the baseline tumour volume, tumour volume change, or percentage tumour volume change at day 90). A p value less than 0.05 was considered to indicate a statistically significant difference.

5. Conclusions

This study proposed a new DTR DCE-MRI processing technique for high spatial resolution interrogation of flow and permeability effects within both the tumour microenvironment and the associated normal-appearing brain regions. Evaluation in a cohort of patients with NF2-related VS undergoing anti-VEGF therapy demonstrated that there was an association between baseline tumour volume and low-level permeability (Ktrans) changes within the normal-appearing brain and that, within the tumour microenvironment, this new approach allowed for concomitant evaluation of blood flow and permeability changes, with the tumoural capillary permeability-surface area product demonstrating the most pronounced reduction at 90 days. In a preliminary study of volumetric response predictors, baseline ve showed the strongest correlation with the change in tumour volume during treatment, and a multivariate model combining ve, PS, and Fp demonstrated high sensitivity and specificity for the prediction of volume reduction at 90 days. These results highlight the utility of this novel DCE-MRI analysis approach in evaluating tumour microvascular changes during treatment and the need for future larger studies investigating its role in predicting antiangiogenic therapy response.

Author Contributions

Conceptualization, K.-L.L. and X.Z.; methodology, K.-L.L., X.Z., D.L. and A.J.; software, K.-L.L. and X.Z.; validation, D.L., X.Z., D.J.C., I.D., A.T.K. and A.J.; formal analysis, K.-L.L. and X.Z.; investigation, D.L., K.-L.L. and X.Z.; resources, X.Z., I.D. and A.J.; data curation, K.-L.L., X.Z. and I.D.; writing—original draft preparation, K.-L.L., D.L. and X.Z.; writing—review and editing, K.-L.L., D.L., X.Z., T.C. and A.J.; visualization, D.L. and X.Z.; supervision, D.J.C., A.T.K., T.C. and A.J.; project administration, X.Z.; funding acquisition, A.J. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by funding from Cancer Research UK (CRUK) and the Engineering and Physical Sciences Research Council (ESPRC) through the Cancer Imaging Centres grant (C8742/A18097).

Institutional Review Board Statement

This study was conducted in accordance with the Declaration of Helsinki and approved by the National Research Ethics Committee North-West (REC reference 13/NW/0247).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study and all recruited participants had previously consented for later analysis of their MRI data.

Data Availability Statement

The datasets generated during and/or analysed during the current study are available from the corresponding author upon reasonable request.

Acknowledgments

We thank David Russell for his participation in patient recruitment; superintendent (research) radiographer Amy Watkins for all her help with scanning the patient participants; and Sha Zhao for all his assistance with performing quality assurance and Eurospin Phantom T1 Studies.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Measurement of CBFET has been previously described and is based on averaging the concentration of early time points during the first pass of a low-dose compact bolus [26]. For dual-injection, dual-temporal resolution DCE-MRI, high spatial resolution CBFET (CBFET-HS) maps can also be derived through rescaling the low spatial resolution CBFET (CBFET-HT) derived from the 4D LDHTaligned volume with the calibration ratio maps generated as part of the LEGATOS procedure. Assuming that the GBCA concentration curve in each HS pixel has a similar first pass shape as its concatenated LDHTaligned arterial phase, the ratio of the HS arterial phase, Ct-HS(t), over the LDHTaligned arterial phase, Ct-HT(t), can be used to convert the CBF estimate derived from the LDHTaligned arterial phase (CBFET-HT) to the CBF of the concatenated HS pixel, CBFET-HS, that is,
CBF ET - HS / CBF ET - HT = C t - HS ( t ) / C t - HT ( t ) , t arterial   phase .
Based on the above, the calibration ratio, ratiocalib, can be defined as follows:
ratio calib = C t - HS ( t adj ) / C t - HT ( t adj ) = CBF ET - HS / CBF ET - HT ,
where tadj is the ending time point of the arterial phase. In the proposed LEGATOS method, Ct-HS(tadj) was calculated as the mean concentration of the five HS frames following the HT arterial phase and Ct-HT(tadj) as the mean concentration of several ending frames of the HT arterial phase series.
The derived 3D calibration ratio map from Equation (A2) can then be used to convert a low spatial resolution CBF map derived from the LDHTaligned arterial phase to high spatial resolution:
CBF ET - HS = CBF ET - HT ratio calib .

Appendix B

List of abbreviations used in dual-injection dual-temporal resolution DCE-MRI pharmacokinetic analysis.
SymbolDefinitionUnits
CBFCerebral blood flowmL min−1 mL−1
CBFETAbsolute cerebral blood flow generated using early time points methodmL min−1 mL−1
CBFET-HSHigh spatial resolution estimates of CBFETmL min−1 mL−1
CBFET-HTLow spatial resolution estimates of CBFETmL min−1 mL−1
Cb(t)Concentration of contrast medium in arterial blood at time tmmol
CpConcentration of contrast medium in arterial blood plasma at time tmmol
Ct(t)Concentration of contrast medium in the voxel at time tmmol
DCE-MRIDynamic contrast-enhanced MRI
DTRDual-temporal resolution
ETMThe extended Tofts model
ET-CBFThe ‘early time points’ method for absolute cerebral blood flow quantification
FP-PPThe hybrid first-pass Patlak plot method
ETWThe early time window, i.e., the time window that meets the microsphere prerequisite.
FDHSFull-dose high spatial resolution
FpPlasma flowmL min−1 mL−1
HSHigh spatial (HS) resolution
HTHigh temporal (HT) resolution
KtransVolume transfer constant between blood plasma and extravascular extracellular spacemin−1
LDHTLow-dose high temporal resolution
LDHTaligned4D LDHT DCE images co-registered to subsequent HS DCE series in dual-injection DTR DCE-MRI
LEGATOSA DTR DCE-MRI data construction technique: the level and rescale the gadolinium contrast concentration curves of high temporal to high spatial
RFpThe ratio of Fp to the sum of Fp and PSNone
PSCapillary permeability–surface area product per unit mass of tissuemL min−1 mL−1
veVolume of the extravascular extracellular space per unit volume of tissuenone
vpFractional blood plasma volumenone

References

  1. O’Connor, J.P.; Jackson, A.; Parker, G.J.; Roberts, C.; Jayson, G.C. Dynamic contrast-enhanced MRI in clinical trials of antivascular therapies. Nat. Rev. Clin. Oncol. 2012, 9, 167–177. [Google Scholar] [CrossRef] [PubMed]
  2. Essig, M.; Anzalone, N.; Combs, S.E.; Dorfler, A.; Lee, S.K.; Picozzi, P.; Rovira, A.; Weller, M.; Law, M. MR imaging of neoplastic central nervous system lesions: Review and recommendations for current practice. AJNR Am. J. Neuroradiol. 2012, 33, 803–817. [Google Scholar] [CrossRef] [PubMed]
  3. Jain, R. Measurements of tumor vascular leakiness using DCE in brain tumors: Clinical applications. NMR Biomed. 2013, 26, 1042–1049. [Google Scholar] [CrossRef] [PubMed]
  4. Li, K.-L.; Djoukhadar, I.; Zhu, X.; Zhao, S.; Lloyd, S.; McCabe, M.; McBain, C.; Evans, D.G.; Jackson, A. Vascular biomarkers derived from dynamic contrast-enhanced MRI predict response of vestibular schwannoma to antiangiogenic therapy in type 2 neurofibromatosis. Neuro-oncology 2016, 18, 275–282. [Google Scholar] [CrossRef] [PubMed]
  5. Lewis, D.; Roncaroli, F.; Agushi, E.; Mosses, D.; Williams, R.; Li, K.-L.; Zhu, X.; Hinz, R.; Atkinson, R.; Wadeson, A.; et al. Inflammation and vascular permeability correlate with growth in sporadic vestibular schwannoma. Neuro-oncology 2019, 21, 314–325. [Google Scholar] [CrossRef] [PubMed]
  6. Larsson, H.B.; Courivaud, F.; Rostrup, E.; Hansen, A.E. Measurement of brain perfusion, blood volume, and blood-brain barrier permeability, using dynamic contrast-enhanced T(1)-weighted MRI at 3 tesla. Magn. Reson. Med. 2009, 62, 1270–1281. [Google Scholar] [CrossRef] [PubMed]
  7. Van Dijken, B.R.J.; van Laar, P.J.; Smits, M.; Dankbaar, J.W.; Enting, R.H.; van der Hoorn, A. Perfusion MRI in treatment evaluation of glioblastomas: Clinical relevance of current and future techniques. J. Magn. Reson. Imaging 2019, 49, 11–22. [Google Scholar] [CrossRef] [PubMed]
  8. Lewis, D.; Donofrio, C.A.; O’leary, C.; Li, K.-L.; Zhu, X.; Williams, R.; Djoukhadar, I.; Agushi, E.; Hannan, C.J.; Stapleton, E.; et al. The microenvironment in sporadic and neurofibromatosis type II–related vestibular schwannoma: The same tumor or different? A comparative imaging and neuropathology study. J. Neurosurg. 2021, 134, 1419–1429. [Google Scholar] [CrossRef]
  9. Jackson, A.; O’Connor, J.P.; Parker, G.J.; Jayson, G.C. Imaging Tumor Vascular Heterogeneity and Angiogenesis using Dynamic Contrast-Enhanced Magnetic Resonance Imaging. Clin. Cancer Res. 2007, 13, 3449–3459. [Google Scholar] [CrossRef]
  10. O’Connor, J.P.B.; Jackson, A.; Parker, G.J.M.; Jayson, G.C. DCE-MRI biomarkers in the clinical evaluation of antiangiogenic and vascular disrupting agents. Br. J. Cancer 2007, 96, 189–195. [Google Scholar] [CrossRef]
  11. Ingrisch, M.; Sourbron, S.; Morhard, D.; Ertl-Wagner, B.; Kumpfel, T.; Hohlfeld, R.; Reiser, M.; Glaser, C. Quantification of perfusion and permeability in multiple sclerosis: Dynamic contrast-enhanced MRI in 3D at 3T. Investig. Radiol. 2012, 47, 252–258. [Google Scholar] [CrossRef] [PubMed]
  12. Rose, C.J.; O’Connor, J.P.; Cootes, T.F.; Taylor, C.J.; Jayson, G.C.; Parker, G.J.; Waterton, J.C. Indexed distribution analysis for improved significance testing of spatially heterogeneous parameter maps: Application to dynamic contrast-enhanced MRI biomarkers. Magn. Reson. Med. 2014, 71, 1299–1311. [Google Scholar] [CrossRef] [PubMed]
  13. Mallio, C.A.; Rovira, À.; Parizel, P.M.; Quattrocchi, C.C. Exposure to gadolinium and neurotoxicity: Current status of preclinical and clinical studies. Neuroradiology 2020, 62, 925–934. [Google Scholar] [CrossRef] [PubMed]
  14. Gulani, V.; Calamante, F.; Shellock, F.G.; Kanal, E.; Reeder, S.B. Gadolinium deposition in the brain: Summary of evidence and recommendations. Lancet Neurol. 2017, 16, 564–570. [Google Scholar] [CrossRef] [PubMed]
  15. Cramer, S.; Simonsen, H.; Frederiksen, J.; Rostrup, E.; Larsson, H. Abnormal blood–brain barrier permeability in normal appearing white matter in multiple sclerosis investigated by MRI. Neuroimage Clin. 2014, 4, 182–189. [Google Scholar] [CrossRef] [PubMed]
  16. Kershaw, L.E.; Buckley, D.L. Precision in measurements of perfusion and microvascular permeability withT1-weighted dynamic contrast-enhanced MRI. Magn. Reson. Med. 2006, 56, 986–992. [Google Scholar] [CrossRef] [PubMed]
  17. Barnes, S.R.; Ng, T.S.C.; Montagne, A.; Law, M.; Zlokovic, B.V.; Jacobs, R. Optimal acquisition and modeling parameters for accurate assessment of low Ktrans blood-brain barrier permeability using dynamic contrast-enhanced MRI. Magn. Reson. Med. 2015, 75, 1967–1977. [Google Scholar] [CrossRef] [PubMed]
  18. Tofts, P.S. Modeling tracer kinetics in dynamic Gd-DTPA MR imaging. J. Magn. Reson. Imaging 1997, 7, 91–101. [Google Scholar] [CrossRef]
  19. Fritz-Hansen, T.; Rostrup, E.; Sørndergaard, L.; Ring, P.B.; Amtorp, O.; Larsson, H.B.W. Capillary transfer constant of Gd-DTPA in the myocardium at rest and during vasodilation assessed by MRI. Magn. Reson. Med. 1998, 40, 922–929. [Google Scholar] [CrossRef]
  20. Li, K.L.; Lewis, D.; Coope, D.J.; Roncaroli, F.; Agushi, E.; Pathmanaban, O.N.; King, A.T.; Zhao, S.; Jackson, A.; Cootes, T.; et al. The LEGATOS technique: A new tissue-validated dynamic contrast-enhanced MRI method for whole-brain, high-spatial resolution parametric mapping. Magn. Reson. Med. 2021, 86, 2122–2136. [Google Scholar] [CrossRef]
  21. Lewis, D.; McHugh, D.J.; Li, K.-L.; Zhu, X.; Mcbain, C.; Lloyd, S.K.; Jackson, A.; Pathmanaban, O.N.; King, A.T.; Coope, D.J. Detection of early changes in the post-radiosurgery vestibular schwannoma microenvironment using multinuclear MRI. Sci. Rep. 2021, 11, 15712. [Google Scholar] [CrossRef] [PubMed]
  22. Tofts, P.S.; Brix, G.; Buckley, D.L.; Evelhoch, J.L.; Henderson, E.; Knopp, M.V.; Larsson, H.B.; Lee, T.-Y.; Mayr, N.A.; Parker, G.J.; et al. Estimating kinetic parameters from dynamic contrast-enhanced t1-weighted MRI of a diffusable tracer: Standardized quantities and symbols. J. Magn. Reson. Imaging 1999, 10, 223–232. [Google Scholar] [CrossRef]
  23. Sasi, S.D.; Gupta, R.K.; Patir, R.; Ahlawat, S.; Vaishya, S.; Singh, A. A comprehensive evaluation and impact of normalization of generalized tracer kinetic model parameters to characterize blood-brain-barrier permeability in normal-appearing and tumor tissue regions of patients with glioma. Magn. Reson. Imaging 2021, 83, 77–88. [Google Scholar] [CrossRef] [PubMed]
  24. Li, K.L.; Zhu, X.; Zhao, S.; Jackson, A. Blood-brain barrier permeability of normal-appearing white matter in patients with vestibular schwannoma: A new hybrid approach for analysis of T1 -W DCE-MRI. J. Magn. Reson. Imaging 2017, 46, 79–93. [Google Scholar] [PubMed]
  25. Sengupta, A.; Agarwal, S.; Gupta, P.K.; Ahlawat, S.; Patir, R.; Gupta, R.K.; Singh, A. On differentiation between vasogenic edema and non-enhancing tumor in high-grade glioma patients using a support vector machine classifier based upon pre and post-surgery MRI images. Eur. J. Radiol. 2018, 106, 199–208. [Google Scholar] [CrossRef]
  26. Li, K.-L.; Lewis, D.; Jackson, A.; Zhao, S.; Zhu, X. Low-dose T1W DCE-MRI for early time points perfusion measurement in patients with intracranial tumors: A pilot study applying the microsphere model to measure absolute cerebral blood flow. J. Magn. Reson. Imaging 2018, 48, 543–557. [Google Scholar] [CrossRef] [PubMed]
  27. Larsson, H.B.; Hansen, A.E.; Berg, H.K.; Rostrup, E.; Haraldseth, O. Dynamic contrast-enhanced quantitative perfusion measurement of the brain usingT1-weighted MRI at 3T. J. Magn. Reson. Imaging 2008, 27, 754–762. [Google Scholar] [CrossRef]
  28. Jackson, A. Imaging microvascular structure with contrast enhanced MRI. Br. J. Radiol. 2003, 76, S159–S173. [Google Scholar] [CrossRef]
  29. Evans, D.R. Neurofibromatosis type 2 (NF2): A clinical and molecular review. Orphanet J. Rare Dis. 2009, 4, 16. [Google Scholar] [CrossRef]
  30. Evans, D.G.; Howard, E.; Giblin, C.; Clancy, T.; Spencer, H.; Huson, S.M.; Lalloo, F. Birth incidence and prevalence of tumor-prone syndromes: Estimates from a UK family genetic register service. Am. J. Med. Genet. A 2010, 152, 327–332. [Google Scholar] [CrossRef]
  31. Stivaros, S.M.; Stemmer-Rachamimov, A.O.; Alston, R.; Plotkin, S.R.; Nadol, J.B.; Quesnel, A.; O’Malley, J.; Whitfield, G.A.; McCabe, M.G.; Freeman, S.R.; et al. Multiple synchronous sites of origin of vestibular schwannomas in neurofibromatosis Type 2. J. Med. Genet. 2015, 52, 557–562. [Google Scholar] [CrossRef] [PubMed]
  32. Lloyd, S.K.; Evans, D.G. Neurofibromatosis type 2 (NF2): Diagnosis and management. Handb. Clin. Neurol. 2013, 115, 957–967. [Google Scholar] [PubMed]
  33. Evans, D.G.; Stivaros, S.M. Multifocality in neurofibromatosis type 2. Neuro-Oncology 2014, 17, 481–482. [Google Scholar] [CrossRef] [PubMed]
  34. Hannan, C.J.; Hammerbeck-Ward, C.; Pathmanaban, O.N.; Smith, M.J.; Rutherford, S.A.; Lloyd, S.K.; Freeman, S.R.M.; Wallace, A.J.; King, A.T.; Evans, D.G.R. Multiple Meningiomas as a Criterion for the Diagnosis of Neurofibromatosis Type 2 and Other Tumor Predisposition Syndromes. Neurosurgery 2022, 90, 793–799. [Google Scholar] [CrossRef] [PubMed]
  35. Plotkin, S.R.; Merker, V.L.; Halpin, C.; Jennings, D.; McKenna, M.J.; Harris, G.J.; Barker, F.G., 2nd. Bevacizumab for progressive vestibular schwannoma in neurofibromatosis type 2: A retrospective review of 31 patients. Otol. Neurotol. 2012, 33, 1046–1052. [Google Scholar] [CrossRef] [PubMed]
  36. Mautner, V.-F.; Nguyen, R.; Kutta, H.; Fuensterer, C.; Bokemeyer, C.; Hagel, C.; Friedrich, R.E.; Panse, J. Bevacizumab induces regression of vestibular schwannomas in patients with neurofibromatosis type 2. Neuro-oncology 2010, 12, 14–18. [Google Scholar] [CrossRef] [PubMed]
  37. Plotkin, S.R.; Stemmer-Rachamimov, A.O.; Barker, F.G., 2nd; Halpin, C.; Padera, T.P.; Tyrrell, A.; Sorensen, A.G.; Jain, R.K.; di Tomaso, E. Hearing improvement after bevacizumab in patients with neurofibromatosis type 2. N. Engl. J. Med. 2009, 361, 358–367. [Google Scholar] [CrossRef] [PubMed]
  38. Sourbron, S.P.; Buckley, D.L. On the scope and interpretation of the Tofts models for DCE-MRI. Magn. Reson. Med. 2011, 66, 735–745. [Google Scholar] [CrossRef]
  39. Chang, Y.-C.C.; Ackerstaff, E.; Tschudi, Y.; Jimenez, B.; Foltz, W.; Fisher, C.; Lilge, L.; Cho, H.; Carlin, S.; Gillies, R.J.; et al. Delineation of Tumor Habitats based on Dynamic Contrast Enhanced MRI. Sci. Rep. 2017, 7, 9746. [Google Scholar] [CrossRef]
  40. Patlak, C.S.; Blasberg, R.G. Graphical evaluation of blood-to-brain transfer constants from multiple-time uptake data. Generalizations. J. Cereb. Blood Flow Metab. 1985, 5, 584–590. [Google Scholar]
  41. Ewing, J.R.; Knight, R.A.; Nagaraja, T.N.; Yee, J.S.; Nagesh, V.; Whitton, P.A.; Li, L.; Fenstermacher, J.D. Patlak plots of Gd-DTPA MRI data yield blood-brain transfer constants concordant with those of14C-sucrose in areas of blood-brain opening. Magn. Reson. Med. 2003, 50, 283–292. [Google Scholar] [CrossRef] [PubMed]
  42. Durukan, A.; Marinkovic, I.; Strbian, D.; Pitkonen, M.; Pedrono, E.; Soinne, L.; Abo-Ramadan, U.; Tatlisumak, T. Post-ischemic blood–brain barrier leakage in rats: One-week follow-up by MRI. Brain Res. 2009, 1280, 158–165. [Google Scholar] [CrossRef] [PubMed]
  43. Abo-Ramadan, U.; Durukan, A.; Pitkonen, M.; Marinkovic, I.; Tatlisumak, E.; Pedrono, E.; Soinne, L.; Strbian, D.; Tatlisumak, T. Post-ischemic leakiness of the blood–brain barrier: A quantitative and systematic assessment by Patlak plots. Exp. Neurol. 2009, 219, 328–333. [Google Scholar] [CrossRef] [PubMed]
  44. Taheri, S.; Gasparovic, C.; Shah, N.J.; Rosenberg, G.A. Quantitative measurement of blood-brain barrier permeability in human using dynamic contrast-enhanced MRI with fast T1 mapping. Magn. Reson. Med. 2011, 65, 1036–1042. [Google Scholar] [CrossRef] [PubMed]
  45. Chang, L.C.; Koay, C.G.; Basser, P.J.; Pierpaoli, C. Linear least-squares method for unbiased estimation of T1 from SPGR signals. Magn. Reson. Med. 2008, 60, 496–501. [Google Scholar] [CrossRef] [PubMed]
  46. Li, K.L.; Zhu, X.P.; Waterton, J.; Jackson, A. Improved 3D quantitative mapping of blood volume and endothelial permeability in brain tumors. J. Magn. Reson. Imaging 2000, 12, 347–357. [Google Scholar] [CrossRef] [PubMed]
  47. Li, K.L.; Zhu, X.P.; Checkley, D.R.; Tessier, J.J.L.; Hillier, V.F.; Waterton, J.C.; Jackson, A. Simultaneous mapping of blood volume and endothelial permeability surface area product in gliomas using iterative analysis of first-pass dynamic contrast enhanced MRI data. Br. J. Radiol. 2003, 76, 39–51. [Google Scholar] [CrossRef] [PubMed]
  48. Batchelor, T.T.; Sorensen, A.G.; di Tomaso, E.; Zhang, W.-T.; Duda, D.G.; Cohen, K.S.; Kozak, K.R.; Cahill, D.P.; Chen, P.-J.; Zhu, M.; et al. AZD2171, a Pan-VEGF Receptor Tyrosine Kinase Inhibitor, Normalizes Tumor Vasculature and Alleviates Edema in Glioblastoma Patients. Cancer Cell 2007, 11, 83–95. [Google Scholar] [CrossRef]
  49. Pishko, G.L.; Muldoon, L.L.; Pagel, M.A.; Schwartz, D.L.; Neuwelt, E.A. Vascular endothelial growth factor blockade alters magnetic resonance imaging biomarkers of vascular function and decreases barrier permeability in a rat model of lung cancer brain metastasis. Fluids Barriers CNS 2015, 12, 5. [Google Scholar] [CrossRef]
  50. Gerstner, E.R.; Emblem, K.E.; Chang, K.; Vakulenko-Lagun, B.; Yen, Y.F.; Beers, A.L.; Dietrich, J.; Plotkin, S.R.; Catana, C.; Hooker, J.M.; et al. Bevacizumab Reduces Permeability and Concurrent Temozolomide Delivery in a Subset of Patients with Recurrent Glioblastoma. Clin. Cancer Res. 2020, 26, 206–212. [Google Scholar] [CrossRef]
  51. Wong, H.K.; Lahdenranta, J.; Kamoun, W.S.; Chan, A.W.; McClatchey, A.I.; Plotkin, S.R.; Jain, R.K.; di Tomaso, E. Anti–Vascular Endothelial Growth Factor Therapies as a Novel Therapeutic Approach to Treating Neurofibromatosis-Related Tumors. Cancer Res. 2010, 70, 3483–3493. [Google Scholar] [CrossRef] [PubMed]
  52. Navis, A.C.; Hamans, B.C.; Claes, A.; Heerschap, A.; Jeuken, J.W.; Wesseling, P.; Leenders, W.P. Effects of targeting the VEGF and PDGF pathways in diffuse orthotopic glioma models. J. Pathol. 2011, 223, 626–634. [Google Scholar] [CrossRef] [PubMed]
  53. Buemi, F.; Guzzardi, G.; Del Sette, B.; Sponghini, A.P.; Matheoud, R.; Soligo, E.; Trisoglio, A.; Carriero, A.; Stecco, A. Apparent diffusion coefficient and tumor volume measurements help stratify progression-free survival of bevacizumab-treated patients with recurrent glioblastoma multiforme. Neuroradiol. J. 2019, 32, 241–249. [Google Scholar] [CrossRef] [PubMed]
  54. Cazzador, D.; Astolfi, L.; Daloiso, A.; Tealdo, G.; Simoni, E.; Mazzoni, A.; Zanoletti, E.; Marioni, G. Tumor Microenvironment in Sporadic Vestibular Schwannoma: A Systematic, Narrative Review. Int. J. Mol. Sci. 2023, 24, 6522. [Google Scholar] [CrossRef] [PubMed]
  55. Varatharaj, A.; Liljeroth, M.; Darekar, A.; Larsson, H.B.; Galea, I.; Cramer, S.P. Blood-brain barrier permeability measured using dynamic contrast-enhanced magnetic resonance imaging: A validation study. J. Physiol. 2018, 597, 699–709. [Google Scholar] [CrossRef] [PubMed]
  56. Heye, A.K.; Thrippleton, M.J.; Armitage, P.A.; Hernández, M.d.C.V.; Makin, S.D.; Glatz, A.; Sakka, E.; Wardlaw, J.M. Tracer kinetic modelling for DCE-MRI quantification of subtle blood–brain barrier permeability. Neuroimage 2016, 125, 446–455. [Google Scholar] [CrossRef] [PubMed]
  57. Deng, X.; Liu, L.; Luo, J.; Liu, L.; Hui, X.; Feng, H. Research on the Mechanism of Cognitive Decline in Patients with Acoustic Neuroma. Front. Neurosci. 2022, 16, 933825. [Google Scholar] [CrossRef] [PubMed]
  58. Deng, X.; Liu, L.; Zhen, Z.; Chen, Q.; Liu, L.; Hui, X. Cognitive decline in acoustic neuroma patients: An investigation based on resting-state functional magnetic resonance imaging and voxel-based morphometry. Front. Psychiatry 2022, 13, 968859. [Google Scholar] [CrossRef]
  59. Deng, X.; Liu, L.; Li, J.; Yao, H.; He, S.; Guo, Z.; Sun, J.; Liu, W.; Hui, X. Brain structural network to investigate the mechanism of cognitive impairment in patients with acoustic neuroma. Front. Aging Neurosci. 2022, 14, 970159. [Google Scholar] [CrossRef]
  60. Wang, X.; Xu, P.; Li, P.; Wang, Z.; Zhao, F.; Gao, Z.; Xu, L.; Luo, Y.-J.; Fan, J.; Liu, P. Alterations in gray matter volume due to unilateral hearing loss. Sci. Rep. 2016, 6, 25811. [Google Scholar] [CrossRef]
  61. Kurtcan, S.; Alkan, A.; Kilicarslan, R.; Bakan, A.A.; Toprak, H.; Aralasmak, A.; Aksoy, F.; Kocer, A. Auditory Pathway Features Determined by DTI in Subjects with Unilateral Acoustic Neuroma. Clin. Neuroradiol. 2016, 26, 439–444. [Google Scholar] [CrossRef] [PubMed]
  62. Li, K.-L.; Buonaccorsi, G.; Thompson, G.; Cain, J.R.; Watkins, A.; Russell, D.; Qureshi, S.; Evans, D.G.; Lloyd, S.K.; Zhu, X.; et al. An improved coverage and spatial resolution-using dual injection dynamic contrast-enhanced (ICE-DICE) MRI: A novel dynamic contrast-enhanced technique for cerebral tumors. Magn. Reson. Med. 2012, 68, 452–462. [Google Scholar] [CrossRef] [PubMed]
  63. Koenig, M.; Klotz, E.; Luka, B.; Venderink, D.J.; Spittler, J.F.; Heuser, L. Perfusion CT of the brain: Diagnostic approach for early detection of ischemic stroke. Radiology 1998, 209, 85–93. [Google Scholar] [CrossRef] [PubMed]
  64. Keil, V.C.; Mädler, B.; Gieseke, J.; Fimmers, R.; Hattingen, E.; Schild, H.H.; Hadizadeh, D.R. Effects of arterial input function selection on kinetic parameters in brain dynamic contrast-enhanced MRI. Magn. Reson. Imaging 2017, 40, 83–90. [Google Scholar] [CrossRef] [PubMed]
  65. Lewis, D.; Zhu, X.; Coope, D.J.; Zhao, S.; King, A.T.; Cootes, T.; Jackson, A.; Li, K.-L. Surrogate vascular input function measurements from the superior sagittal sinus are repeatable and provide tissue-validated kinetic parameters in brain DCE-MRI. Sci. Rep. 2022, 12, 8737. [Google Scholar] [CrossRef] [PubMed]
  66. Ashburner, J.; Friston, K. Multimodal Image Coregistration and Partitioning—A Unified Framework. NeuroImage 1997, 6, 209–217. [Google Scholar] [CrossRef]
Figure 1. Flowchart showing the key steps in the multi-model LEGATOS integrated kinetic analysis. Blue highlights the key steps of the dual-injection, dual-temporal resolution DCE-MRI acquisition; red highlights the key steps of the LEGATOS data reconstruction to generate a high spatiotemporal resolution GBCA concentration volume; lilac highlights the key kinetic models used; and black highlights the derived microvascular parameters. ET-CBF = ‘early time points’ method for absolute cerebral blood flow quantification; HS = high spatial resolution; FDHS = full-dose high spatial resolution DCE-MRI data; Ktrans = volume transfer constant between blood plasma and extravascular extracellular space; LDHTaligned = low-dose high temporal resolution DCE-MRI data co-registered to the FDHS DCE series; Fp = plasma flow, where Fp = CBFET·(1 − Hct); FP-PP = the hybrid first-pass Patlak plot model; GBCA= gadolinium-based contrast agent; LEGATOS = level and rescale the gadolinium contrast concentration curves of high temporal to high spatial DCE-MRI; PS = permeability-surface area product; ve = volume of the extravascular extracellular space per unit volume of tissue; vp = fractional blood plasma volume.
Figure 1. Flowchart showing the key steps in the multi-model LEGATOS integrated kinetic analysis. Blue highlights the key steps of the dual-injection, dual-temporal resolution DCE-MRI acquisition; red highlights the key steps of the LEGATOS data reconstruction to generate a high spatiotemporal resolution GBCA concentration volume; lilac highlights the key kinetic models used; and black highlights the derived microvascular parameters. ET-CBF = ‘early time points’ method for absolute cerebral blood flow quantification; HS = high spatial resolution; FDHS = full-dose high spatial resolution DCE-MRI data; Ktrans = volume transfer constant between blood plasma and extravascular extracellular space; LDHTaligned = low-dose high temporal resolution DCE-MRI data co-registered to the FDHS DCE series; Fp = plasma flow, where Fp = CBFET·(1 − Hct); FP-PP = the hybrid first-pass Patlak plot model; GBCA= gadolinium-based contrast agent; LEGATOS = level and rescale the gadolinium contrast concentration curves of high temporal to high spatial DCE-MRI; PS = permeability-surface area product; ve = volume of the extravascular extracellular space per unit volume of tissue; vp = fractional blood plasma volume.
Pharmaceuticals 16 01282 g001
Figure 2. Use of the LEGATOS technique with the ‘early time points’ absolute cerebral blood flow (ET-CBF) method to generate high spatial resolution CBFET estimates. (A) Representative pre-treatment images from an NF2 patient with multiple meningiomas including a right parasagittal convexity meningioma (white arrow). (B) Representative images from a patient with bilateral NF2-related VS. Note the imaging artifact from the left-sided bone-anchored hearing aid within the post-contrast T1-weighted image. The low spatial resolution CBFET map (CBFET-HT) derived from the low dose high temporal resolution DCE data is rescaled by the LEGATOS calibration ratio map to generate the high spatial resolution CBFEHS map. T1W + C = T1-weighted image post-contrast.
Figure 2. Use of the LEGATOS technique with the ‘early time points’ absolute cerebral blood flow (ET-CBF) method to generate high spatial resolution CBFET estimates. (A) Representative pre-treatment images from an NF2 patient with multiple meningiomas including a right parasagittal convexity meningioma (white arrow). (B) Representative images from a patient with bilateral NF2-related VS. Note the imaging artifact from the left-sided bone-anchored hearing aid within the post-contrast T1-weighted image. The low spatial resolution CBFET map (CBFET-HT) derived from the low dose high temporal resolution DCE data is rescaled by the LEGATOS calibration ratio map to generate the high spatial resolution CBFEHS map. T1W + C = T1-weighted image post-contrast.
Pharmaceuticals 16 01282 g002
Figure 3. (A) Representative post-contrast T1W image and high spatial resolution Ktrans maps in an NF2 patient with multiple meningioma. Data were fitted using either the hybrid FP-PP approach (LEGATOSFP-PP, middle) or ETM (LEGATOSETM, right). Note the comparatively high Ktrans values seen within the right-frontal and left-parasagittal meningioma on both maps. Compared to LEGATOSETM the LEGTOSFP-PP permits high spatial resolution assessment of low-level vascular permeability within the normal-appearing brain. (B) Correlation between mean NAGM/NAWM Ktrans and VS volume. Top row: correlation between average VS volume and mean NAGM (left)/NAWM (right) Ktrans values; bottom row: correlation between cumulative VS volume and mean NAGM (left)/NAWM (right) Ktrans values. NAGM = normal-appearing grey matter; NAWM = normal-appearing white matter. T1W + C = T1-weighted image post-contrast.
Figure 3. (A) Representative post-contrast T1W image and high spatial resolution Ktrans maps in an NF2 patient with multiple meningioma. Data were fitted using either the hybrid FP-PP approach (LEGATOSFP-PP, middle) or ETM (LEGATOSETM, right). Note the comparatively high Ktrans values seen within the right-frontal and left-parasagittal meningioma on both maps. Compared to LEGATOSETM the LEGTOSFP-PP permits high spatial resolution assessment of low-level vascular permeability within the normal-appearing brain. (B) Correlation between mean NAGM/NAWM Ktrans and VS volume. Top row: correlation between average VS volume and mean NAGM (left)/NAWM (right) Ktrans values; bottom row: correlation between cumulative VS volume and mean NAGM (left)/NAWM (right) Ktrans values. NAGM = normal-appearing grey matter; NAWM = normal-appearing white matter. T1W + C = T1-weighted image post-contrast.
Pharmaceuticals 16 01282 g003
Figure 4. (A) High spatial resolution tumour Ktrans, PS, Fp, and RFp (= Fp/(Fp + PS)) maps imaged at day 0 (baseline) and day 90 in an NF2 patient with left-sided VS. (B) Histograms of tumour voxel values for Ktrans, Fp, vp, PS, RFp, and ve at baseline (the solid line) and 90 days post-bevacizumab treatment (dashed line), calculated in the same tumour as in panel A. Tumour volume size was 6.01 cm3 (3005 voxels) at day 0 and 4.33 cm3 (2163 voxels) at day 90. The overall tumour volume was reduced by 28%, and tumour volume loss was principally seen in voxels with Fp in the range of 0.2 to 0.8 min−1 (red arrow lines) and vp in the range of 0.02 to 0.10 in the histogram. Tumour mean ± SD (median) for each of the 6 parameters measured at day 0 and day 90 are shown in the panel corresponding to that parameter, respectively. T1W + C = T1-weighted image post-contrast. The x-axis of each histogram represents the range of parameter values divided into bins, and the y-axis represents frequency (how many voxels fall within each bin).
Figure 4. (A) High spatial resolution tumour Ktrans, PS, Fp, and RFp (= Fp/(Fp + PS)) maps imaged at day 0 (baseline) and day 90 in an NF2 patient with left-sided VS. (B) Histograms of tumour voxel values for Ktrans, Fp, vp, PS, RFp, and ve at baseline (the solid line) and 90 days post-bevacizumab treatment (dashed line), calculated in the same tumour as in panel A. Tumour volume size was 6.01 cm3 (3005 voxels) at day 0 and 4.33 cm3 (2163 voxels) at day 90. The overall tumour volume was reduced by 28%, and tumour volume loss was principally seen in voxels with Fp in the range of 0.2 to 0.8 min−1 (red arrow lines) and vp in the range of 0.02 to 0.10 in the histogram. Tumour mean ± SD (median) for each of the 6 parameters measured at day 0 and day 90 are shown in the panel corresponding to that parameter, respectively. T1W + C = T1-weighted image post-contrast. The x-axis of each histogram represents the range of parameter values divided into bins, and the y-axis represents frequency (how many voxels fall within each bin).
Pharmaceuticals 16 01282 g004
Figure 5. Predictive value of pre-treatment tumour median values of microvascular parameters. (AE) Comparison of the pre-treatment (day 0) microvascular parameters (Ktrans, PS, Fp, vp, and ve) between the responder (res) and non-responder (non) VS groups. Boxplots show mean (dot on the box plot), median (bar within the box), upper and lower quartiles (box limits), and extreme values (whiskers). The p values are calculated using a two-sided Student t-test. (F) Correlation between pre-treatment tumour median ve and percentage change in tumour volume at day 90; volumed0 = baseline tumour volume; volumed90 = tumour volume at day 90.
Figure 5. Predictive value of pre-treatment tumour median values of microvascular parameters. (AE) Comparison of the pre-treatment (day 0) microvascular parameters (Ktrans, PS, Fp, vp, and ve) between the responder (res) and non-responder (non) VS groups. Boxplots show mean (dot on the box plot), median (bar within the box), upper and lower quartiles (box limits), and extreme values (whiskers). The p values are calculated using a two-sided Student t-test. (F) Correlation between pre-treatment tumour median ve and percentage change in tumour volume at day 90; volumed0 = baseline tumour volume; volumed90 = tumour volume at day 90.
Pharmaceuticals 16 01282 g005
Table 1. Ktrans statistics estimated from segmented normal-appearing grey matter (NAGM) and white matter (NAWM) in twelve patients with NF2.
Table 1. Ktrans statistics estimated from segmented normal-appearing grey matter (NAGM) and white matter (NAWM) in twelve patients with NF2.
Ktrans Statistics of NAGM/NAWM SegmentsNAGM Segment Statistics
Mean ± SD
NAWM Segment Statistics
Mean ± SD
p Value
Paired t-Test
Mean (min−1)0.00089 ± 0.00110.00062 ± 0.000300.30
Median (min−1)0.00047 ± 0.000370.00028 ± 0.000340.03
SD (min−1)0.0030 ± 0.00260.0012 ± 0.000210.04
Max (min−1)0.050 ± 0.0340.014 ± 0.00580.002
Min (min−1)−0.023 ± 0.0074−0.015 ± 0.00400.002
Significant p values are shown in bold type.
Table 2. Mean and standard deviation of the tumour median Ktrans, Fp, PS, RFp, vp, and ve estimated at day 0 and day 90 from 20 VSs in 12 NF2 patients.
Table 2. Mean and standard deviation of the tumour median Ktrans, Fp, PS, RFp, vp, and ve estimated at day 0 and day 90 from 20 VSs in 12 NF2 patients.
Microvascular ParameterDay 0
Mean ± SD (Median)
Day 90
Mean ± SD (Median)
p Value
Paired t-Test
Ktrans (min−1)
Res (N = 12)0.121 ± 0.023 (0.121)0.083 ± 0.031 (0.082)0.001
Non (N = 8)0.095 ± 0.037 (0.100)0.086 ± 0.025 (0.096)0.18
All (N = 20)0.111 ± 0.059 (0.123)0.085 ± 0.028 (0.093)0.0008
PS (mL min−1 mL−1)
Res (N = 12)0.169 ± 0.039 (0.174)0.109 ± 0.043 (0.106)0.0002
Non (N = 8)0.125 ± 0.053 (0.128)0.107 ± 0.032 (0.118)0.10
All (N = 20)0.151 ± 0.087 (0.162)0.108 ± 0.038 (0.117)0.0001
Fp (mL min−1 mL−1)
Res (N = 12)0.514 ± 0.114 (0.496)0.414 ± 0.103 (0.423)0.08
Non (N = 8)0.485 ± 0.175 (0.470)0.541 ± 0.139 (0.579)0.36
All (N = 20)0.502 ± 0.270 (0.477)0.465 ± 0.132 (0.452)0.37
RFp (%)
Res (N = 12)0.752 ± 0.065 (0.744)0.792 ± 0.048 (0.788)0.004
Non (N = 8)0.784 ± 0.063 (0.790)0.826 ± 0.045 (0.820)0.04
All (N = 20)0.765 ± 0.114 (0.761)0.806 ±0.049 (0.799)0.0003
vp (%)
Res (N = 12)0.047 ± 0.012 (0.043)0.038 ± 0.009 (0.040)0.07
Non (N = 8)0.045 ± 0.019 (0.040)0.045 ± 0.010 (0.047)0.94
All (N = 20)0.046 ± 0.025 (0.042)0.040 ± 0.010 (0.041)0.13
ve (%)
Res (N = 12)0.519 ± 0.047 (0.523)0.498 ± 0.103 (0.457)0.56
Non (N = 8)0.431 ± 0.042 (0.423)0.511 ± 0.063 (0.500)0.01
All (N = 20)0.484 ± 0.259 (0.500)0.504 ± 0.196 (0.488)0.43
Abbreviations: Res = responders, Non = non-responders. Significant p values are shown in bold type.
Table 3. Correlations between tumour median Ktrans, Fp, PS, vp, and ve estimated at day 0 and tumour volumetric parameters.
Table 3. Correlations between tumour median Ktrans, Fp, PS, vp, and ve estimated at day 0 and tumour volumetric parameters.
Linear Regression Analysis
(N = 20)
Tumour Volume
(cm3; Day 0)
Tumour Volume Change
(cm3; Day 90)
Percentage Tumour Volume Change
(%; Day 90)
Ktrans (min−1)R2 = 0.19
p = 0.05
R2 = 0.13
p = 0.12
R2 = 0.08
p = 0.23
PS (mL min−1 mL−1)R2 = 0.24
p = 0.03
R2 = 0.18
p = 0.06
R2 = 0.10
p = 0.18
Fp (mL min−1 mL−1)R2 = 0.01
p = 0.71
R2 = 0.01
p = 0.67
R2 = 0.01
p = 0.66
vp (%)R2 = 0.02
p = 0.55
R2 = 0.03
p = 0.48
R2 = 0.00
p = 0.90
ve (%)R2 = 0.16
p = 0.08
R2 = 0.21
p = 0.04
R2 = 0.56
p < 0.001
Bold type indicates statistically significant p ≤ 0.05.
Table 4. Area under the receiver operator curve, sensitivity, and specificity of the binomial regression model for the prediction of tumour response using univariate (A) and multivariate analysis (B).
Table 4. Area under the receiver operator curve, sensitivity, and specificity of the binomial regression model for the prediction of tumour response using univariate (A) and multivariate analysis (B).
Prediction of ResponseAUC-ROC
(p Value)
SensitivitySpecificityOverall
Classification
A. Univariate analysis
ve (%)0.896
(0.024)
0.8300.8750.850
PS (mL min−1 mL−1)0.708
(0.10)
0.9200.5000.750
Ktrans (min−1)0.688
(0.11)
0.920.3750.700
Fp (mL min−1 mL−1)0.667
(0.50)
1.000.1250.650
vp (%)0.615
(0.61)
---
B. Multivariate analysis with backward selection
Step 1
ve + PS + Ktrans + Fp + vp
0.948
(0.031; 0.92; 0.93; 0.69; 0.83)
0.8300.8750.850
Step 2
ve + PS + Fp + vp
0.948
(0.030; 0.31; 0.70; 0.78)
0.8300.8750.850
Step 3
ve + PS + Fp
0.948
(0.032; 0.35; 0.70)
0.8300.8750.850
Step 4
ve + PS
0.938
(0.036; 0.37)
0.8300.7500.800
Step 5
ve
0.896
(0.024)
0.8300.8750.850
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, K.-L.; Lewis, D.; Zhu, X.; Coope, D.J.; Djoukhadar, I.; King, A.T.; Cootes, T.; Jackson, A. A Novel Multi-Model High Spatial Resolution Method for Analysis of DCE MRI Data: Insights from Vestibular Schwannoma Responses to Antiangiogenic Therapy in Type II Neurofibromatosis. Pharmaceuticals 2023, 16, 1282. https://doi.org/10.3390/ph16091282

AMA Style

Li K-L, Lewis D, Zhu X, Coope DJ, Djoukhadar I, King AT, Cootes T, Jackson A. A Novel Multi-Model High Spatial Resolution Method for Analysis of DCE MRI Data: Insights from Vestibular Schwannoma Responses to Antiangiogenic Therapy in Type II Neurofibromatosis. Pharmaceuticals. 2023; 16(9):1282. https://doi.org/10.3390/ph16091282

Chicago/Turabian Style

Li, Ka-Loh, Daniel Lewis, Xiaoping Zhu, David J. Coope, Ibrahim Djoukhadar, Andrew T. King, Timothy Cootes, and Alan Jackson. 2023. "A Novel Multi-Model High Spatial Resolution Method for Analysis of DCE MRI Data: Insights from Vestibular Schwannoma Responses to Antiangiogenic Therapy in Type II Neurofibromatosis" Pharmaceuticals 16, no. 9: 1282. https://doi.org/10.3390/ph16091282

APA Style

Li, K. -L., Lewis, D., Zhu, X., Coope, D. J., Djoukhadar, I., King, A. T., Cootes, T., & Jackson, A. (2023). A Novel Multi-Model High Spatial Resolution Method for Analysis of DCE MRI Data: Insights from Vestibular Schwannoma Responses to Antiangiogenic Therapy in Type II Neurofibromatosis. Pharmaceuticals, 16(9), 1282. https://doi.org/10.3390/ph16091282

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop