Towards stable non-invasive estimation of Interstitial Fluid Pressure incorporating diffusion weighted data
Charles Baker1,2, Nicholas Dowson1, Michael Fay3, Rosalind Jeffree4, and Stephen Rose1,2

1CSIRO, Brisbane, Australia, 2University of Queensland, Brisbane, Australia, 3Genesis Cancer Care, Newcastle, Australia, 4Royal Brisbane and Women's Hospital, Brisbane, Australia

Synopsis

Interstitial Fluid Pressure (IFP) is an important potential reason behind the typically ineluctable progression of late stage glioma in the face of radical treatment. It can be inferred from dynamic contrast enhanced (DCE) MRI images, by voxel-wise fitting of a Tofts model to extract the plasma flux per unit volume and subsequent solving of a pressure equation. Ideally model fitting should be numerically stable and should account for the non-isotropic directionality of fluid flow. This paper describes our approach to achieve this.

Purpose

IFP could explain the typically ineluctable progression of late stage glioma, which despite maximal safe surgical resection and subsequent chemo-radiotherapy has a median survival of 15 months (with all subgroups combined) 1. Direct measurement of IFP using a pressure transducer is invasive, and only the measurement at the point of surgical ingress will remain free of bias arising from tissue movement subsequent to incision. Dynamic contrast enhanced (DCE) MRI enables IFP to be inferred non-invasively by modelling the rate of fluid flux from the region of blood brain barrier breakdown 2, and utilising the empirically established physical properties of tissue 3. The model requires an accurate characterisation of the rate of flux from the blood into tissue and out again, which is numerically stable and ideally should also account for the preferential directionality of fluid flow, which is obtainable from Diffusion Weighted Images (DWI). This paper defines an approach to achieve this.

Methods

A two compartment (plasma and tissue) model has been found to describe the temporal changes of Gadolinium concentration within tissues, , in regions of blood brain breakdown well, if the influx, $$$ K_{in}$$$, and outflux, $$$ k_{out}$$$, constants are allowed to differ as in the extended Tofts model 4:

$$\frac{dC_t}{dt} = K_{in}C_p(t) - k_{out}C_t(t) \qquad\qquad\qquad (1) $$

where is the plasma concentration at time t. The flow of fluid through tissue is governed by both hydrostatic pressure (pushing plasma into the interstitium in the arterioles) and osmotic pressure (extracting plasma from interstitium in venules) 5, which is mediated by the conductivity of the tissue as expressed by Darcy’s Law 6. Where tumours disrupt normal vasculature, the natural equilibrium in pressure is disrupted, resulting in an influx of fluid into the tumour and surrounding regions 3. The fluid flux, $$$ \frac{J_V}{V}$$$, per unit volume from blood into tissue may be computed from the difference between and , after compensating for the osmotic reflection coefficient of Gadolinium, $$$\sigma$$$ 3. Assuming dissipation of excess pressure within the parenchyma is slow, but that steady state has been reached: conservation of mass allows the computation of local pressure, p, from the flux in the region of blood brain barrier breakdown:

$$ \nabla\cdot\frac{J_V}{v}=\nabla\cdot\frac{K_{in}(\mathbf{x}-k_{out}(\mathbf{x}))}{1-\sigma} = \nabla\cdot K\nabla p \qquad\qquad\qquad (2)$$,

where K is the hydraulic conductivity. For a scalar K, $$$\nabla\cdot\nabla$$$ can be encoded as a discrete Laplacian matrix operator and p solved directly. Accounting for local plasma volume, $$$v_p$$$, (1) can be formulated to be optimised over a single non-linear parameter for the time-varying concentration within each voxel, $$$T_i$$$,at location $$$\mathbf{x}_i$$$ over discrete times, $$$t_u$$$:

$$\arg_{k_{out}}\min \arg_{v_p, k_{in}}\min \sum_u \left( T_i(t_u)-[v_p+K_{in}e^{-k_{out}\cdot t}]*C_p(t)\right)^2 \qquad\qquad\qquad (3) $$

while the linear parameters are optimised using a standard non-negative least squares technique 7. The optimisation equation (3) was implemented using an open source dynamic analysis platform 8. Note that $$$K_{in}$$$ should be inversely scaled to account for plasma fraction, $$$v_p$$$. The hydraulic conductivity, K, is not necessarily isotropic as fluid preferentially flows parallel to white matter tracts 9,10. If available, DWI data can be incorporated into the matrix operator in (2), where discrete voxels correspond to distinct rows and columns and element ij encodes the relationship between voxels i and j. Only neighbouring voxels, $$$j\in N_i$$$ have non-zero values. Directional preferences of fluid flow can be encoded as reweightings, so if $$$ K_i^{ij}$$$ is the conductivity of voxel i, in the unit direction $$$(x_i-x_j$$$, the Laplacian operator becomes:

$$\nabla\cdot K\nabla = \left\{ \begin{matrix} -\frac{1}{2}\sum_{j\in N_i} K_i^{ij}+K_j^{ij} & i=j \\ \frac{1}{2}(K_i^{ij}+K_j^{ij}) & i\neq j, j \in N_i \\ 0&\textrm{elsewhere}\\ \end{matrix}\right. \qquad\qquad\qquad (4) $$

Results

Fig. 1 shows the pressure map of a glioma patient overlayed on the post-contrast and FLAIR image of the patient. $$$C_p$$$ was estimated from a voxel-cluster in the sagittal sinus of 15x55second DCE. (3) yielded numerically stable influx parameters that varied smoothly over the spatial volume, without regularisation. This was used as an input to compute pressure using (2) and (4).

Discussion

The measurement of IFP using a numerically stable formulation of the extended Tofts model and incorporating the directional preferences of fluid flow has been demonstrated. IFP corresponds to the enhancing region of oedema in the FLAIR image. Future work will explore inclusion of the fact that flow will not traverse sulci, which is expected to further improve the correspondence with the FLAIR image.

Acknowledgements

This work was partly funded by a grant from the Cancer Council Foundation, Queensland. The funding bodies have not contributed to the study design, the collection, management, analysis and interpretation of data, the writing of final reports or the decision to submit findings for publication. No other authors have potential conflicts of interest to declare.

References

1. Stupp R, Hegi ME, Mason WP, et al. Effects of radiotherapy with concomitant and adjuvant temozolomide versus radiotherapy alone on survival in glioblastoma in a randomised phase III study: 5-year analysis of the EORTC-NCIC trial. Lancet Oncol. 2009;10(5):459-66. doi:10.1016/S1470-2045(09)70025-7.

2. Zhao J, Salmon H, Sarntinoranont M. Effect of heterogeneous vasculature on interstitial transport within a solid tumor. Microvasc. Res. 2007;73(3):224-36. doi:10.1016/j.mvr.2006.12.003.

3. Baxter LT, Jain RK. Transport of fluid and macromolecules in tumors. I. Role of interstitial pressure and convection. Microvasc. Res. 1989;37(1):77-104.

4. Tofts PS. Modeling tracer kinetics in dynamic Gd-DTPA MR imaging. J. Magn. Reson. Imaging 1997;7(1):91-101. doi:10.1002/jmri.1880070113.

5. Starling EH. On the Absorption of Fluids from the Connective Tissue Spaces. J. Physiol. 1896;19(4):312-26.

6. Whitaker S. Flow in porous media I: A theoretical derivation of Darcy’s law. Transp. Porous Media 1986;1(1):3-25. doi:10.1007/BF01036523.

7. Lawson CL, Hanson RJ, Society for Industrial and Applied Mathematics. Solving Least Squares Problems, Society for Industrial and Applied Mathematics 1974.

8. Dowson N, Baker C, Raffelt D, Thomas P, Smith J, Rose S. InsightToolkit Kinetic Analysis (itk::ka) Library. Insight J. (ISSN 2327-770X) 2014. doi:10.4225/08/540E9A7D11EB0.

9. Beaulieu C. The basis of anisotropic water diffusion in the nervous system - a technical review. NMR Biomed. 2002;15(7-8):435-55. doi:10.1002/nbm.782.

10. Rose SE, Chen F, Chalk JB, Zelaya FO, Strugnell WE, Benson M, Semple J, Doddrell DM. Loss of connectivity in Alzheimer’s disease: an evaluation of white matter tract integrity with colour coded MR diffusion tensor imaging. J. Neurol. Neurosurg. Psychiatry 2000;69(4):528-530. doi:10.1136/jnnp.69.4.528.

Figures

IFP overlayed onto the post-constrast and FLAIR image of a glioma patient before surgery. Interstitial fluid pressure is shown on a log scale covering four orders of magnitude. The region of increased pressure corresponds to the region of FLAIR enhancement.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
2524