How much microvascular anatomy is in T1-DCE MRI? – A computerized analysis of meningioma microvasculature correlated to kinetic parameters applying the extended Tofts model
Vera Catharina Keil1, Kanishka Hiththetiya2, Gerrit H. Gielen3, Matthias Simon4, Bogdan Pintea4, Anna Vogelgesang1, Juergen Gieseke1,5, Burkhard Maedler5, Hans Heinz Schild1, and Dariusch Reza Hadizadeh1

1Department of Radiology, Universitätsklinikum Bonn, Bonn, Germany, 2Center for Pathology, Universitätsklinikum Bonn, Bonn, Germany, 3Department of Neuropathology, Universitätsklinikum Bonn, Bonn, Germany, 4Clinic for Neurosurgery and Stereotaxy, Universitätsklinikum Bonn, Bonn, Germany, 5Philips Healthcare, Best, Netherlands

Synopsis

Tissue perfusion is co-defined by anatomical factors of the microvasculature. If and how kinetic parameters of T1w dynamic-contrast enhanced (T1-DCE) MRI fit into this anatomical perfusion model, is a topic of on-going discussion. Based on the extended Tofts model (ETK) we therefore performed a focal analysis of kinetic parameters in meningioma and surgically retrieved precisely corresponding tissue specimens. Their microvasculature underwent multimodal computerized analysis. Kinetic parameters were found to correlate poorly and inconsistently with the microvascular anatomy on both an inter- and intra-individual level.

Target Audience

Scientists interested in the physiological and histological basis of T1-DCE MRI as well as clinicians applying T1-DCE in neuroradiology.

Purpose

To elucidate if kinetic parameters of T1-DCE MRI based on the extended Tofts model (ETK) have a foundation in the microvascular anatomy of meningioma [1].

Introduction

In T1-DCE MRI, kinetic parameters, such as transfer constants Ktrans and kep but also the extracellular-extravascular (ve) and the plasma volume fractions (vp), relate to vascular density and surface area. They are presumed to characterize tissue microvasculature and cellularity [2]. T1-DCE MRI is thus increasingly clinically applied to deduce characteristics of tissue composition. However, few studies investigate if this correlation with microvascular anatomy actually exists [3,4].

Bearing in mind that models to derive kinetic parameters are known to be vulnerable to factors such as vascular permeability, global vascularization and perfusion, correlations depend on both tissue and applied model [5,6]. Meningioma is a usually highly perfused and vascularized extra-axial tumor entity that seems ideal to elucidate whether and respectively to what extent ETK-derived kinetic parameters are related to microvascular anatomy based on precisely matching tissue specimens [1,7].

Methods

19 patients with confirmed meningiomas (15 WHO I, 3 WHO III, 1 WHO II) underwent a whole-brain T1-DCE MRI [3.0T Achieva TX (Philips Healthcare, Best/NL); 8-channel head coil; axial orientation; 36 slices; 1.57x1.6x3mm acq. voxel size; dual FA technique (5°/15°; TE 2.3ms); 50 dynamics at 5.3s each (FA 8°; 1.7ms); contrast: 0.1mmol/kg BW Gadobutrol (Bayer Healthcare, Leverkusen/D); 24ml saline flush; flow rate 3ml/s].

The ETK, defined by C(t) = vp* Ca(t) + Ktrans*e−tkep* Ca(t), was applied to derive Ktrans, kep, ve and vp (fig. 1, Intellispace Portal 5.0, Philips Healthcare) [1]. Arterial vessels from internal and external carotid artery feed a meningioma, but cannot be distinguished on imaging. However, as meningiomas certainly drain into large venous sinuses the ROI for vascular input function (VIF) calculation was placed in the superior sagittal sinus. VIF considered hematocrit and contrast relaxivity 5.0 L/(mmol*s)].

In each patient 1 to 5 biopsy regions were marked (5x5mm) and tissue based kinetic parameters were measured. T1-DCE MRI data sets were exported to neuronavigation software (Brainlab, Feldkirchen/D). The marked biopsy specimens were collected during total resection of the meningioma (64 specimens in total). These were fixated, sliced and immunohistochemically stained with CD34 (marker for endothelial cells). The slices were scanned with tissue quantification software (Tissue Studio, Definiens, fig. 2). Vessels were divided into three size groups: small<2mm², intermediate<5mm², large>5mm². Software was trained to reliably detect vessels with and without lumen. For subsequent automated analysis the tissue was marked, then divided into 0.6*0.6mm ROI grids. Representative ROIs were checked for correct registration of vessel wall, lumen and size segmentation (fig. 3). Data retrieved was: 1. vessel density (vessels/mm²), 2. relative area of vasculature in total slice area (%), 3. number of vessels by size group, 4. distribution of vessel sizes per slice and 5. mean vessel wall thickness.

All data was split into subgroups of vessels with and without lumen. Kinetic parameter data (Ktrans, kep, ve, vp) was inter- and intra-individually analyzed for correlation with the microvascular data (mixed linear model; SAS9.4, SAS Institute Inc., Cary/USA.)

Results

The vascular fraction ranged from 0.5 to 14.9% of the total tissue area. It did not significantly correlate with any of the kinetic parameters on an inter-individual basis (e.g.: vp: p=0.954, correlation r=0.01; Ktrans: p=0.642; correlation r=0.09). Vascular size, relative fractions of vascular sizes and thickness of vascular lumina were also not significantly correlated to any of the kinetic parameters.

The inter-individual variability of kinetic parameters was expectedly high (Ktrans: 0.28±0.29/min; kep: 0.99±0.79/min; vp: 36.3±57.0%; ve: 27.9±21.9%). Therefore an intra-individual correlation was performed despite the low number of 3 to 5 biopsy points per patient. This rendered inconsistent results (e.g.: vp vs. vascular fraction of total tissue: r=-0.99 to +0.89; mean r: 0.11; Ktrans vs. vascular fraction of total tissue: r= -0.99 to+0.99; mean r: 0.04).

Discussion

Quite different from initial expectations kinetic parameters of T1-DCE MRI did not correlate with any of the microvascular tissue factors anylyzed in this multimodal setting. Pronounced inter-individual variability may not explain this lack of correlation, as intra-individual results were highly inconsistent. Future studies applying different kinetic and VIF models are needed to rule out a bias from ill fitting of the models applied in this piece of research.

Conclusion

Kinetic parameters of T1-DCE MRI cannot directly be deduced from the microvascular architecture in meningioma. Clinical users of T1-DCE MRI should therefore be careful to derive anatomical tissue aspects from their kinetic data.

Acknowledgements

Thanks to the technical staff of both MRI and neuropathology unit and to Mr. H. Peusens (Philips Healthcare) for technical advice.

References

[1] Tofts PS et al. JMRI. 1999 Sep; 10(3):223-32.

[2] Tofts PS and Kermode AG. MRM. 1991 Feb; 17 (2):357-67.

[3] Jia ZZ et al. EurJRadiol. 2015 Sep; 84(9):1805-9.

[4] Van Niekerk CG et al. Eur Radiol. 2014 Oct; 24(10):2597-605.

[5] Sourbron SP and Buckley DL. MRM. 2011 Sep; 66(3):735-45.

[6] Cuenod CA and Balvay D. DiagnIntervImaging. 2013 Dec; 94(12):1187-204.

[7] Mawrin C and Perry A. J Neurooncol. 2010 Sep; 99(3):379-91.

Figures

Fig. 1 ROI-based analysis of kinetic parameters (red dot) on color maps of kinetic parameters (here l.t.r. Ktrans, vp, ve) was applied to define tissue for neurosurgical retrieval and histological matching.

Fig. 2 The kinetic parameter ROI matched tissue specimens were fixated, sliced and CD34-stained (a). The tissue slice areas for software analysis were manually defined (b).

Fig. 3 Microvasculature was quantified based on high resolution grids (a, vessel endothelium in brown). After training and empirical definition of identification thresholds, software could reliably distinguish between vessel walls (b, lilac) and lumen (green) as well as separate vessels into small (c, yellow), intermediate (orange) and large (red) sizes. Vascular data was then compared to kinetic parameters.



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