Background

Glioblastoma (GBM) is the deadliest and most common brain tumor in adults, with poor survival and response to aggressive therapy. Limited access of drugs to tumor cells is one reason for such grim clinical outcomes. A driving force for therapeutic delivery is interstitial fluid flow (IFF), both within the tumor and in the surrounding brain parenchyma. However, convective and diffusive transport mechanisms are understudied. In this study, we examined the application of a novel image analysis method to measure fluid flow and diffusion in GBM patients.

Methods

Here, we applied an imaging methodology that had been previously tested and validated in vitro, in silico, and in preclinical models of disease to archival patient data from the Ivy Glioblastoma Atlas Project (GAP) dataset. The analysis required the use of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI), which is readily available in the database. The analysis results, which consisted of IFF flow velocity and diffusion coefficients, were then compared to patient outcomes such as survival.

Results

We characterized IFF and diffusion patterns in patients. We found strong correlations between flow rates measured within tumors and in the surrounding parenchymal space, where we hypothesized that velocities would be higher. Analyzing overall magnitudes indicated a significant correlation with both age and survival in this patient cohort. Additionally, we found that neither tumor size nor resection significantly altered the velocity magnitude. Lastly, we mapped the flow pathways in patient tumors and found a variability in the degree of directionality that we hypothesize may lead to information concerning treatment, invasive spread, and progression in future studies.

Conclusions

An analysis of standard DCE-MRI in patients with GBM offers more information regarding IFF and transport within and around the tumor, shows that IFF is still detected post-resection, and indicates that velocity magnitudes correlate with patient prognosis.

Contrast transport models are widely used to quantify blood flow and transport in dynamic contrast-enhanced magnetic resonance imaging. These models analyze the time course of the contrast agent concentration, providing diagnostic and prognostic value for many biological systems. Thus, ensuring accuracy and repeatability of the model parameter estimation is a fundamental concern. In this work, we analyze the structural and practical identifiability of a class of nested compartment models pervasively used in analysis of MRI data. We combine artificial and real data to study the role of noise in model parameter estimation. We observe that although all the models are structurally identifiable, practical identifiability strongly depends on the data characteristics. We analyze the impact of increasing data noise on parameter identifiability and show how the latter can be recovered with increased data quality. To complete the analysis, we show that the results do not depend on specific tissue characteristics or the type of enhancement patterns of contrast agent signal.

Glioblastoma (GBM), a highly aggressive form of brain tumor, is a disease marked by extensive invasion into the surrounding brain. Interstitial fluid flow (IFF), or the movement of fluid within the spaces between cells, has been linked to increased invasion of GBM cells. Better characterization of IFF could elucidate underlying mechanisms driving this invasion in vivo. Here, we develop a technique to noninvasively measure interstitial flow velocities in the glioma microenvironment of mice using dynamic contrast-enhanced magnetic resonance imaging (MRI), a common clinical technique. Using our in vitro model as a phantom "tumor" system and in silico models of velocity vector fields, we show we can measure average velocities and accurately reconstruct velocity directions. With our combined MR and analysis method, we show that velocity magnitudes are similar across four human GBM cell line xenograft models and the direction of fluid flow is heterogeneous within and around the tumors, and not always in the outward direction. These values were not linked to the tumor size. Finally, we compare our flow velocity magnitudes and the direction of flow to a classical marker of vessel leakage and bulk fluid drainage, Evans blue. With these data, we validate its use as a marker of high and low IFF rates and IFF in the outward direction from the tumor border in implanted glioma models. These methods show, for the first time, the nature of interstitial fluid flow in models of glioma using a technique that is translatable to clinical and preclinical models currently using contrast-enhanced MRI.

Contrast transport models are widely used to quantify blood flow and transport in dynamic contrast-enhanced magnetic resonance imaging. These models analyze the time course of the contrast agent concentration, providing diagnostic and prognostic value for many biological systems. Thus, ensuring accuracy and repeatability of the model parameter estimation is a fundamental concern. In this work, we analyze the structural and practical identifiability of a class of nested compartment models pervasively used in analysis of MRI data. We combine artificial and real data to study the role of noise in model parameter estimation. We observe that although all the models are structurally identifiable, practical identifiability strongly depends on the data characteristics. We analyze the impact of increasing data noise on parameter identifiability and show how the latter can be recovered with increased data quality. To complete the analysis, we show that the results do not depend on specific tissue characteristics or the type of enhancement patterns of contrast agent signal.

Background

Dynamic contrast-enhanced MRI (DCE-MRI) parameters have been shown to be biomarkers for treatment response in glioblastoma (GBM). However, variations in analysis and measurement methodology complicate determination of biological changes measured via DCE. The aim of this study is to quantify DCE-MRI variations attributable to analysis methodology and image quality in GBM patients.

Methods

The Extended Tofts model (eTM) and Leaky Tracer Kinetic Model (LTKM), with manually and automatically segmented vascular input functions (VIFs), were used to calculate perfusion kinetic parameters from 29 GBM patients with double-baseline DCE-MRI data. DCE-MRI images were acquired 2-5 days apart with no change in treatment. Repeatability of kinetic parameters was quantified with Bland-Altman and percent repeatability coefficient (%RC) analysis.

Results

The perfusion parameter with the least RC was the plasma volume fraction (v _p ), with a %RC of 53%. The extra-cellular extra-vascular volume fraction (v _e ) %RC was 82% and 81%, for extended Tofts-Kety Model (eTM) and LTKM respectively. The %RC of the volume transfer rate constant (K ^trans ) was 72% for the eTM, and 82% for the LTKM, respectively. Using an automatic VIF resulted in smaller %RCs for all model parameters, as compared to manual VIF.

Conclusions

As much as 72% change in K ^trans (eTM, autoVIF) can be attributable to non-biological changes in the 2-5 days between double-baseline imaging. Poor K ^trans repeatability may result from inferior temporal resolution and short image acquisition time. This variation suggests DCE-MRI repeatability studies should be performed institutionally, using an automatic VIF method and following quantitative imaging biomarkers alliance guidelines.

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) is a routine method to noninvasively quantify perfusion dynamics in tissues. The standard practice for analyzing DCE-MRI data is to fit an ordinary differential equation to each voxel. Recent advances in data science provide an opportunity to move beyond existing methods to obtain more accurate measurements of fluid properties. Here, we developed a localized convolutional function regression that enables simultaneous measurement of interstitial fluid velocity, diffusion, and perfusion in 3D. We validated the method computationally and experimentally, demonstrating accurate measurement of fluid dynamics in situ and in vivo. Applying the method to human MRIs, we observed tissue-specific differences in fluid dynamics, with an increased fluid velocity in breast cancer as compared to brain cancer. Overall, our method represents an improved strategy for studying interstitial flows and interstitial transport in tumors and patients. We expect that our method will contribute to the better understanding of cancer progression and therapeutic response.