Radiology

Open journal

ISSN 2473-0947

Bio-Mathematical Modelling in Tumor Evaluation Via Magnetic Resonance Imaging

Mario Sansone*

Mario Sansone, PhD

Assistant Professor, Department of Electrical Engineering and Information Technologies, University of Naples Federico II, Via Claudio 21, Naples 8013, Italy; E-mail: msansone@unina.it

In the last 25 years, within the field of tumor evaluation using magnetic resonance imaging (MRI), 2 main promising applications have emerged that can provide useful functional information beyond the morphological images: dynamic contrast enhanced MRI (DCE-MRI) and diffusion weighted imaging (DWI or DW-MRI). DCE-MRI involves the administration of a contrast medium that, flowing into the vascular network feeding the tumor, might give information on tissue vascular perfusion which is related to the angiogenesis phenomenon, in turn associated to tumor growth. DWI exploits the water brownian motion that is affected by tissue cellularity and vascularization. The successful application of both DCE-MRI and DWI lays strongly on both technical MRI improvements achieved in the last decades and on adequate bio-mathematical modelling of tissues. In this editorial, I will briefly discuss a few main issues of both applications focussing on bio-mathematical modelling.

DCE-MRI

The main technical issues in DCE-MRI concern the trade-off between time-space resolution and accurate quantification of contrast medium concentration. They both affect the accuracy of bio-mathematical modelling typically used in DCE assessment of the tumor.

In the early years of DCE, radiologists mainly looked at the shape of the time-intensity curve of a single voxel or averaged regions of interest (ROIs); in fact, it had been observed that, after contrast medium administration, a typical aggressive (malignant) tumor showed a wash-in (absorption) and a wash-out (excretion) phase while a benign lesion showed a slow absorption rate. However, this type of analysis was strongly reader-dependent and did not allow to easily and quantitatively  compare subsequent studies of the same subject or studies from different subjects. It was clear that a more quantitative approach was needed.

Several groups (Tofts et al1, Larsson et al2, Brix et al3) proposed to use bio-mathematical models to analyze DCE data. These early models are based mainly on tissue compartmentalization (plasma and extracellular compartments) and the interesting parameters are the kinetics constants of exchange across compartments that are related to the permeability of the vessel walls: in fact, immature, leaky vessels (quickly grown stimulated with the vascular endothelial growing factor (VEGF) by the tumor) should show a higher permeability.4 However, these models are a strong simplification of reality because of limitation of the technical equipments available at that time. In fact, the need for high resolution images prevented fast acquisition of multiple series.

Today, a number of fast pulse sequences of acquisition give the possibility to use more realistic bio-mathematical models: different groups Koh et al5, Schabel et al6, Sourbron7 have proposed more sophisticated models including some aspects not considered by the first generation models: the rough approximation of only 2 compartments within a voxel has been abandoned and a more realistic distribution of compartments is considered instead: moreover, capillary transit-time (the time required for a contrast medium particle to travel across a capillary) and vascular permeability are considered random variables with a certain distribution within the voxel (region) of interest.

The complexity of the model and the measurement uncertainty pose some limits in the evaluation of the kinetics parameters (permeability, vascular fraction, transit time, etc.): the accuracy (repeatability) with which these parameters can be evaluated is intrinsically limited by mathematical constraints.8 Moreover, accurate measurement of the contrast medium concentration within the voxel is still an issue because fast pulse sequences can overcome the time-space resolution trade-off but typically introduce some additional uncertainty in the measurement.9 

Therefore, it is important for the radiologists evaluating the parametric maps to be supported by specialised professionals (such as biomedical engineers) in the interpretation of these maps in order to take into account the quality of the map reconstruction.

DWI 

DWI was born for brain perfusion assessment10 (and has quickly evolved leading to the tractography11 that is not considered further here) and has been quickly becoming an important tool for the evaluation of tumors in other body parts such as prostate12 and breast.13 The diffusion of water (brownian motion) has been studied by many physicists among which Einstein gave an important contribution: in fact, considering the diffusion as a random process he established the relation between the diffusion constant, the time of observation and the path of a water particle. Diffusion within biological tissues depends on tissue organization: in order to evaluate whether there are some preferred directions of motion, the concept of diffusion tensor imaging (roughly, diffusion constants along 3 orthogonal directions) has been introduced.14 

During DWI acquisition random moving spins are affected by magnetic gradients in such a way that fast diffusing water molecules contribute to attenuate the signal from a specific voxel. The magnetic gradients parameters are called b-values. In this way, qualitative maps can be obtained in which voxels with fast diffusing water are represented by darker colors.15 However, this type of images give information on the apparent diffusion coefficient (ADC) because true molecular water diffusion is restricted and hindered by vascular tortuosity.

To overcome this issue, the first bio-mathematical model of water diffusion within biological tissues has been developed by Le Bihan et al16 in the late 1980’s. This model, called intra-voxel-incoherent-motion (IVIM), again considered the tissues within a voxel as organized in compartments (capillaries and extracellular space). In this model the interesting quantities are the vascular fraction, the water diffusion constant and the pseudo-diffusion constant. This latter has been introduced to account for the random orientation of vessels in tissues.

As in the case of DCE-MRI, also in the case of DW-MRI the quantitative analysis is limited by mathematical constraints. Specifically, the b-values used in DW acquisition strongly affect the accuracy of parameters estimation. Very sophisticated solutions have been proposed in the field of tractography14; moreover, some studies have been performed in the field of tumor evaluation17; however, this issue is still under debate and no definitive solution is currently available.

Finally, a very interesting issue is the link between the different models. As the mentioned models describe different measurements of related quantities (liquid flow and diffusion across compartments within tissues) it seems not unreasonable that model parameters should be related across different models. A link between perfusion and IVIM models was first suggested theoretically by Le Bihan et al.18 This link, if definitively shown, might be used to cross-validate DCE and DW studies and possibly to improve the accuracy of parametric maps. However, in the field of tumor studies, it seems that this issue has not been yet addressed specifically: a very small number of studies has been performed in this direction13,19 but they did not show this link in a conclusive manner.

1. Tofts PS, Kermode AG. Measurement of the blood-brain barrier permeability and leakage space using dynamic MR imaging. 1. Fundamental concepts. Magn Reson Med. 1991; 17(2): 357-367. doi: 10.1002/mrm.1910170208

2. Larsson HBW, Tofts PS. Measurement of blood-brain barrier permeability using dynamic Gd-DTPA scanning-a comparison of methods. Magn Reson Med. 1992; 24(1): 174-176. doi: 10.1002/mrm.1910240119

3. Brix G, Semmler W, Port R, Schad LR, Layer G, Lorenz WJ. Pharmacokinetic parameters in CNS Gd-DTPA enhanced MR imaging. J Comput Assist Tomogr. 1991; 15(4): 621-628. Web site. http://journals.lww.com/jcat/abstract/1991/07000/pharmacokinetic_parameters_in_cns_gd_dtpa_enhanced.18.aspx. Accessed August 24, 2016.

4. Alan J, Buckley DL, Parker GJM. Dynamic Contrast-Enhanced Magnetic Resonance Imaging in Oncology. Berlin, Germany: Springer; 2005: 309.

5. Koh TS, Cheong LH, Tan CK, Lim CC. A distributed parameter model of cerebral blood-tissue exchange with account of capillary transit time distribution. Neuroimage. 2006; 30(2): 426-435. Web site. http://search.proquest.com/openview/6ae90c9b0b9898df4dbf5cc177a6c42b/1?pq-origsite=gscholar&cbl=2031077. Accessed August 24, 2016

6. Schabel MC. A unified impulse response model for DCE-MRI. Magn Reson Med. 2012; 68(5): 1632-1646. doi: 10.1002/mrm.24162

7. Sourbron S. A tracer-kinetic field theory for medical imaging. IEEE Trans Med Imaging. 2014; 33(4): 935-946. doi: 10.1109/TMI.2014.2300450

8. Bates DM, Watts DG. Nonlinear Regression Analysis and Its Applications. 2nd ed. USA: John Wiley & Sons, Inc; 1988: 87.

9. Sansone M, Fusco R, Petrillo A. Accuracy of contrast agent quantification in MRI: A comparison between two k-space sampling schemes. Appl Magn Reson. 2015; 46: 1283. doi: 10.1007/s00723-015-0718-8

10. Le Bihan D, Breton E, Lallemand D, Grenier P, Cabanis E, Laval-Jeantet M. MR imaging of intravoxel incoherent motions: application to diffusion and perfusion in neurologic disorders. Radiology. 1986; 161(2): 401-407. Web site. http://www.meteoreservice.com/PDFs/DLBRadiology86.pdf. Accessed August 24, 2016

11. Kinosada Y, Ono M, Okuda Y, et al. MR tractography–visualization of structure of nerve fiber system from diffusion weighted images with maximum intensity projection method. Nihon Igaku Hoshasen Gakkai Zasshi. 1993; 53(2): 171-179. Web site. http://eurekamag.com/research/008/982/008982405.php. Accessed August 24, 2016.

12. Fusco R, Sansone M, Petrillo M, et al. Multiparametric MRI for prostate cancer detection: Preliminary results on quantitative analysis of dynamic contrast enhanced imaging, diffusion-weighted imaging and spectroscopy imaging. Magn Reson Imaging. 34(7): 839-845. doi: 10.1016/j.mri.2016.04.001

13. Fusco R, Sansone M, Filice S, Granata V, Catalano O, Amato DM. Integration of DCE-MRI and DW-MRI quantitative parameters for breast lesion classification. BioMed Research International. 2015; 237863: 12. doi: 10.1155/2015/237863

14. Minati L, Węglarz WP. Physical foundations, models, and methods of diffusion magnetic resonance imaging of the brain: A review. Concepts Magn Reson. 2007; 30A: 278-307. doi: 10.1002/cmr.a.20094

15. Stejskal EO, Tanner JE. Spin diffusion measurements: Spin echoes in the presence of a time-dependent field gradient. J Chem Phys. 1965; 42(1): 288. doi: 10.1063/1.1695690

16. Le Bihan D. Intravoxel incoherent motion imaging using steady-state free precession. Magn Reson Med. 1988; 7(3): 346-351. doi: 10.1002/mrm.1910070312

17. Lemke A, Stieltjes B, Schad LR, Laun FB. Toward an optimal distribution of b values for intravoxel incoherent motion imaging. Magn Reson Imaging. 2011; 29(6): 766-776. doi: 10.1016/j.mri.2011.03.004

18. Le Bihan D, Turner R. The capillary network: A link between IVIM and classical perfusion. Magn Reson Med. 1992; 27(1): 171-178. Web site. http://meteoreservice.com/PDFs/LeBihan92.pdf. Accessed August 24, 2016

19. Sansone M, Fusco R, Petrillo A. An attempt to integrate diffusion weighted and dynamic contrast enhanced MRI. Radiol Open J. 1(1): 31-34. doi: 10.17140/ROJ-1-105

LATEST ARTICLES

Unraveling the Mysteries of Type-A Aortic Dissection Using POCUS/Echocardiography

Syeda Rukh*, Sathyanarayana Machani and Milind Awale

doi.

Blood Sample from the Patient

Hypertriglyceridemia-Induced Pancreatitis: A Case Report and Literature Review

Maarten Bulterys, Melvin Willems* and Agnes Meersman

doi.

From Neck Pain to a Life-Threatening Condition: A Case Report

Floris Vandewoude* and Sören Verstraete

doi.

LATEST ARTICLES