Characterizing Heterogeneity within Head and Neck Lesions Using Cluster Analysis of Multi-Parametric MRI Data

Oropharyngeal cancer patient stratification using random forest based-learning over high-dimensional radiomic features

To improve risk prediction for oropharyngeal cancer (OPC) patients using cluster analysis on the radiomic features extracted from pre-treatment Computed Tomography (CT) scans. 553 OPC Patients randomly split into training (80%) and validation (20%), were classified into 2 or 3 risk groups by applying hierarchical clustering over the co-occurrence matrix obtained from a random survival forest (RSF) trained over 301 radiomic features. The cluster label was included together with other clinical data to train an ensemble model using five predictive models (Cox, random forest, RSF, logistic regression, and logistic-elastic net). Ensemble performance was evaluated over the independent test set for both recurrence free survival (RFS) and overall survival (OS). The Kaplan-Meier curves for OS stratified by cluster label show significant differences for both training and testing (p val < 0.0001). When compared to the models trained using clinical data only, the inclusion of the cluster label improves AUC test performance from .62 to .79 and from .66 to .80 for OS and RFS, respectively. The extraction of a single feature, namely a cluster label, to represent the high-dimensional radiomic feature space reduces the dimensionality and sparsity of the data. Moreover, inclusion of the cluster label improves model performance compared to clinical data only and offers comparable performance to the models including raw radiomic features.

Pre-treatment functional MRI of breast cancer: T2* evaluation at 3 T and relationship to dynamic contrast-enhanced and diffusion-weighted imaging

PURPOSE

Baseline T2* relaxation time has been proposed as an imaging biomarker in cancer, in addition to Dynamic Contrast-Enhanced (DCE) MRI and diffusion-weighted imaging (DWI) parameters. The purpose of the current work is to investigate sources of error in T2* measurements and the relationship between T2* and DCE and DWI functional parameters in breast cancer.

METHODS

Five female volunteers and thirty-two women with biopsy proven breast cancer were scanned at 3 T, with Research Ethics Committee approval. T2* values of the normal breast were acquired from high-resolution, low-resolution and fat-suppressed gradient-echo sequences in volunteers, and compared. In breast cancer patients, pre-treatment T2*, DCE MRI and DWI were performed at baseline. Pathologically complete responders at surgery and non-responders were identified and compared. Principal component analysis (PCA) and cluster analysis (CA) were performed.

RESULTS

There were no significant differences between T2* values from high-resolution, low-resolution and fat-suppressed datasets (p > 0.05). There were not significant differences between baseline functional parameters in responders and non-responders (p > 0.05). However, there were differences in the relationship between T2* and contrast-agent uptake in responders and non-responders. Voxels of similar characteristics were grouped in 5 clusters, and large intra-tumoural variations of all parameters were demonstrated.

CONCLUSION

Breast T2* measurements at 3 T are robust, but spatial resolution should be carefully considered. T2* of breast tumours at baseline is unrelated to DCE and DWI parameters and contribute towards describing functional heterogeneity of breast tumours.

Linking imaging to omics utilizing image-guided tissue extraction

Imaging provides an insight into biological patho-mechanisms of diseases. However, the link between the imaging phenotype and the underlying molecular processes is often not well understood. Methods such as metabolomics and proteomics reveal detailed information about these processes. Unfortunately, they provide no spatial information and thus cannot be easily correlated with functional imaging. We have developed an image-guided milling machine and unique workflows to precisely isolate tissue samples based on imaging data. The tissue samples remain cooled during the entire procedure, preventing sample degradation. This enables us to correlate, at an unprecedented spatial precision, comprehensive imaging information with metabolomics and proteomics data, leading to a better understanding of diseases.

Phenotypic heterogeneity is commonly observed in diseased tissue, specifically in tumors. Multimodal imaging technologies can reveal tissue heterogeneity noninvasively in vivo, enabling imaging-based profiling of receptors, metabolism, morphology, or function on a macroscopic scale. In contrast, in vitro multiomics, immunohistochemistry, or histology techniques accurately characterize these heterogeneities in the cellular and subcellular scales in a more comprehensive but ex vivo manner. The complementary in vivo and ex vivo information would provide an enormous potential to better characterize a disease. However, this requires spatially accurate coregistration of these data by image-driven sampling as well as fast sample-preparation methods. Here, a unique image-guided milling machine and workflow for precise extraction of tissue samples from small laboratory animals or excised organs has been developed and evaluated. The samples can be delineated on tomographic images as volumes of interest and can be extracted with a spatial accuracy better than 0.25 mm. The samples remain cooled throughout the procedure to ensure metabolic stability, a precondition for accurate in vitro analysis.

Investigating intracranial tumour growth patterns with multiparametric MRI incorporating Gd-DTPA and USPIO-enhanced imaging

High grade and metastatic brain tumours exhibit considerable spatial variations in proliferation, angiogenesis, invasion, necrosis and oedema. Vascular heterogeneity arising from vascular co-option in regions of invasive growth (in which the blood-brain barrier remains intact) and neoangiogenesis is a major challenge faced in the assessment of brain tumours by conventional MRI. A multiparametric MRI approach, incorporating native measurements and both Gd-DTPA (Magnevist) and ultrasmall superparamagnetic iron oxide (P904)-enhanced imaging, was used in combination with histogram and unsupervised cluster analysis using a k-means algorithm to examine the spatial distribution of vascular parameters, water diffusion characteristics and invasion in intracranially propagated rat RG2 gliomas and human MDA-MB-231 LM2-4 breast adenocarcinomas in mice. Both tumour models presented with higher ΔR1 (the change in transverse relaxation rate R1 induced by Gd-DTPA), fractional blood volume (fBV) and apparent diffusion coefficient than uninvolved regions of the brain. MDA-MB-231 LM2-4 tumours were less densely cellular than RG2 tumours and exhibited substantial local invasion, associated with oedema, whereas invasion in RG2 tumours was minimal. These additional features were reflected in the more heterogeneous appearance of MDA-MB-231 LM2-4 tumours on T2 -weighted images and maps of functional MRI parameters. Unsupervised cluster analysis separated subregions with distinct functional properties; areas with a low fBV and relatively impermeable blood vessels (low ΔR1 ) were predominantly located at the tumour margins, regions of MDA-MB-231 LM2-4 tumours with relatively high levels of water diffusion and low vascular permeability and/or fBV corresponded to histologically identified regions of invasion and oedema, and areas of mismatch between vascular permeability and blood volume were identified. We demonstrate that dual contrast MRI and evaluation of tissue diffusion properties, coupled with cluster analysis, allows for the assessment of heterogeneity within invasive brain tumours and the designation of functionally diverse subregions that may provide more informative predictive biomarkers.