Histogram analysis of the microvasculature of intracerebral human and murine glioma xenografts

Histogram analysis of MR quantitative parameters: are they correlated with prognostic factors in prostate cancer?

PURPOSE

To investigate the correlation between quantitative MR parameters and prognostic factors in prostate cancer (PCa).

METHOD

A total of 186 patients with pathologically confirmed PCa who underwent preoperative multiparametric MRI (mpMRI), including synthetic MRI (SyMRI), were enrolled from two medical centers. The histogram metrics of SyMRI [T1, T2, proton density (PD)] and apparent diffusion coefficient (ADC) values were extracted. The Mann‒Whitney U test or Student's t test was employed to determine the association between these histogram features and the prognostically relevant factors. Receiver operating characteristic (ROC) curve analysis was conducted to evaluate the differentiation performance. Spearman's rank correlation coefficients were calculated to determine the correlations between histogram parameters and the International Society of Urological Pathology (ISUP) grade group as well as pathological T stage.

RESULTS

Significant correlations were found between the histogram parameters and the ISUP grade as well as pathological T stage of PCa. Among these histogram parameters, ADC_minimum had the strongest correlation with the ISUP grade (r = - 0.481, p < 0.001), and ADC_Median showed the strongest association with pathological T stage (r = - 0.285, p = 0.008). The ADC_10th percentile exhibited the highest performance in identifying clinically significant prostate cancer (csPCa) (AUC 0.833; 95% CI 0.771-0.883). When discriminating between the status of different prognostically relevant factors, a significant difference was observed between extraprostatic extension-positive and -negative cancers with regard to histogram parameters of the ADC map (10th percentile, 90th percentile, mean, median, minimum) and T1 map (minimum) (p = 0.002-0.032). Moreover, histogram parameters of the ADC map (90th percentile, maximum, mean, median), T2 map (10th percentile, median), and PD map (10th percentile, median) were significantly lower in PCa with perineural invasion (p = 0.009-0.049). The T2 values were significantly lower in patients with seminal vesicle invasion (minimum, p = 0.036) and positive surgical margin (10th percentile, 90th percentile, mean, median, and minimum, p = 0.015-0.025).

CONCLUSION

Quantitative histogram parameters derived from synthetic MRI and ADC maps may have great potential for predicting the prognostic features of PCa.

Radiomics analysis of pre-treatment F-18 FDG PET/CT for predicting response to transarterial radioembolization in liver tumors

Aim: To investigate the relationship between the textural features extracted from pre-treatment fluorine-18 fluorodeoxyglucose positron emission with computed tomography (F-18 FDG PET/CT) and the response to treatment in patients undergoing transarterial radioembolization (TARE) due to primary or metastatic liver tumors. Material and Method: A total of 25 liver lesions from the pre-treatment F-18 PET/CT images of 14 patients were segmented manually. Standard uptake value (SUV) metrics and radiomics features were extracted for each lesion. Metabolic treatment response was determined according to PERCIST criteria in 18F-FDG PET/CT imaging performed 2 months after the treatment. Feature selection was done with recursive feature elimination (RFE). The association between selected features and treatment response was evaluated with logistic regression analysis. Results: Eventually, 13 lesions responded to TARE, while 12 lesions remain stable or progressed. All standard uptake values and 27 out of 30 textural heterogeneity indicators were significantly higher in lesions that responded to treatment. SUVmax, kurtosis and dissimilarity features were selected by the RFE algorithm for the prediction of response to TARE. Logistic regression analysis revealed that all three parameters were significantly associated with treatment outcome. Conclusion: Textural features extracted from pre-treatment F-18 FDG PET/CT in patients undergoing TARE due to liver tumors are promising biomarkers that can be potentially used to predict metabolic treatment response.

Application of Synthetic MRI for Direct Measurement of Magnetic Resonance Relaxation Time and Tumor Volume at Multiple Time Points after Contrast Administration: Preliminary Results in Patients with Brain Metastasis

Objective

The purpose of this study was to investigate the time-dependent effects of contrast medium on multi-dynamic, multi-echo (MDME) sequence in patients with brain metastases.

Materials and Methods

This study included 7 patients with 15 brain metastases who underwent magnetic resonance (MR) examination which included MDME sequences at 1 minute, 10 minutes and 20 minutes after contrast injection. Two volumes of interests, covering an entire tumor (whole tumor) and the enhancing portion of the tumor, were derived from post-contrast synthetic T1-weighted images. Statistical comparisons were performed for three different time delays for histogram parameters of the longitudinal relaxation rate (R₁) and the transverse relaxation rate (R₂), and lesion volumes.

Results

The mean and the median of R₁ and the mean of R₂ in both the whole tumor and the inner enhancing portion were larger on the 10 minutes delayed images than on the 1 minute or 20 minutes delayed images (mean of R₁ in the whole tumor on the 1 minute, 10 minutes, and 20 minutes delayed images: 1.26 ms, 1.39 ms, and 1.37 ms; mean of R₁ in the inner enhancing portion: 1.43 ms, 1.53 ms and 1.44 ms; all p < 0.017). The volumes of the whole tumor and the inner enhancing portion were significantly larger in the 10 minutes and 20 minutes delayed images than on the 1 minute delayed images (all p < 0.017).

Conclusion

Magnetic resonance relaxation times and the volumes of the whole tumor and the inner enhancing portion were measured larger on the 10 minutes or 20 minutes delayed images than on the 1 minute delayed images. The MDME sequence immediately after contrast injection cannot fully reflect the effects of gadolinium-based contrast agent leakage in the tissue.

Histogram analysis of stretched-exponential and monoexponential diffusion-weighted imaging models for distinguishing low and intermediate/high gleason scores in prostate carcinoma

BACKGROUND

Noninvasive measures to evaluate the aggressiveness of prostate carcinoma (PCa) may benefit patients.

PURPOSE

To assess the value of stretched-exponential and monoexponential diffusion-weighted imaging (DWI) for predicting the aggressiveness of PCa.

STUDY TYPE

Retrospective study.

SUBJECTS

Seventy-five patients with PCa.

FIELD STRENGTH

3T DWI examinations were performed using b-values of 0, 500, 1000, and 2000 s/mm² .

ASSESSMENT

The research were based on entire-tumor histogram analysis and the reference standard was radical prostectomy.

STATISTICAL TESTS

The correlation analysis was programmed with Spearman's rank-order analysis between the histogram variables and Gleason grade group (GG). Receiver operating characteristic (ROC) regression was used to analyze the ability of these histogram variables to differentiate low-grade (LG) from intermediate/high-grade (HG) PCa.

RESULTS

The percentiles and mean of apparent diffusion coefficient (ADC) and distributed diffusion coefficient (DDC) were correlated with GG (ρ: 0.414-0.593), while there was no significant relation among α value, skewnesses, and kurtosises with GG (ρ:0.034-0.323). HG tumors (ADC:484 ± 136, 592 ± 139, 670 ± 144, 788 ± 146, 895 ± 141 mm² /s; DDC: 410 ± 142, 532 ± 172, 666 ± 193, 786 ± 196, 914 ± 181 mm² /s) had lower values in the 10^th , 25^th , 50^th , 75^th percentiles and means than LG tumors (ADC: 644 ± 779, 737 ± 84, 836 ± 83, 919 ± 82, 997 ± 107 mm² /s; DDC: 552 ± 82, 680 ± 94, 829 ± 112, 931 ± 106, 1045 ± 100 mm² /s). However, there was no difference between LG and HG tumors in α value (0.671 ± 0.041 vs. 0.633 ± 0.114), kurtosises (ADC 0.09 vs. 0.086; DDC -0.033 vs. -0.317), or skewnesses (ADC -0.036 vs. 0.073; DDC -0.063 vs. 0.136). The above statistics were P < 0.01. ADC10 with AUC = 0.840 and DDC10 with AUC = 0.799 were similar in discriminating between LG and HG PCa at P < 0.05.

DATA CONCLUSION

Histogram variables of DDC and ADC may predict the aggressiveness of PCa, while α value does not. The abilities of ADC10 and DDC10 to discriminate LG from HG tumors were similar, and both better than their respective means.

LEVEL OF EVIDENCE

3 Technical Efficacy: Stage 1 J. MAGN. RESON. IMAGING 2018;48:491-498.

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.

MRI Evaluation of Non-Necrotic T2-Hyperintense Foci in Pediatric Diffuse Intrinsic Pontine Glioma

BACKGROUND AND PURPOSE

The conventional MR imaging appearance of diffuse intrinsic pontine glioma suggests intralesional histopathologic heterogeneity, and various distinct lesion components, including T2-hypointense foci, have been described. Here we report the prevalence, conventional MR imaging semiology, and advanced MR imaging features of non-necrotic T2-hyperintense foci in diffuse intrinsic pontine glioma.

MATERIALS AND METHODS

Twenty-five patients with diffuse intrinsic pontine gliomas were included in this study. MR imaging was performed at 3T by using conventional and advanced MR imaging sequences. Perfusion (CBV), vascular permeability (v_e, K^trans), and diffusion (ADC) metrics were calculated and used to characterize non-necrotic T2-hyperintense foci in comparison with other lesion components, namely necrotic T2-hyperintense foci, T2-hypointense foci, peritumoral edema, and normal brain stem. Statistical analysis was performed by using Kruskal-Wallis and Wilcoxon rank sum tests.

RESULTS

Sixteen non-necrotic T2-hyperintense foci were found in 12 tumors. In these foci, ADC values were significantly higher than those in either T2-hypointense foci (P = .002) or normal parenchyma (P = .0002), and relative CBV values were significantly lower than those in either T2-hypointense (P = .0002) or necrotic T2-hyperintense (P = .006) foci. Volume transfer coefficient values in T2-hyperintense foci were lower than those in T2-hypointense (P = .0005) or necrotic T2-hyperintense (P = .0348) foci.

CONCLUSIONS

Non-necrotic T2-hyperintense foci are common, distinct lesion components within diffuse intrinsic pontine gliomas. Advanced MR imaging data suggest low cellularity and an early stage of angioneogenesis with leaky vessels resulting in expansion of the extracellular space. Because of the lack of biopsy validation, the underlying histoarchitectural and pathophysiologic changes remain unclear; therefore, these foci may correspond to a poorly understood biologic event in tumor evolution.

DEMARCATE: Density-based magnetic resonance image clustering for assessing tumor heterogeneity in cancer

Tumor heterogeneity is a crucial area of cancer research wherein inter- and intra-tumor differences are investigated to assess and monitor disease development and progression, especially in cancer. The proliferation of imaging and linked genomic data has enabled us to evaluate tumor heterogeneity on multiple levels. In this work, we examine magnetic resonance imaging (MRI) in patients with brain cancer to assess image-based tumor heterogeneity. Standard approaches to this problem use scalar summary measures (e.g., intensity-based histogram statistics) that do not adequately capture the complete and finer scale information in the voxel-level data. In this paper, we introduce a novel technique, DEMARCATE (DEnsity-based MAgnetic Resonance image Clustering for Assessing Tumor hEterogeneity) to explore the entire tumor heterogeneity density profiles (THDPs) obtained from the full tumor voxel space. THDPs are smoothed representations of the probability density function of the tumor images. We develop tools for analyzing such objects under the Fisher-Rao Riemannian framework that allows us to construct metrics for THDP comparisons across patients, which can be used in conjunction with standard clustering approaches. Our analyses of The Cancer Genome Atlas (TCGA) based Glioblastoma dataset reveal two significant clusters of patients with marked differences in tumor morphology, genomic characteristics and prognostic clinical outcomes. In addition, we see enrichment of image-based clusters with known molecular subtypes of glioblastoma multiforme, which further validates our representation of tumor heterogeneity and subsequent clustering techniques.

Improving tumour heterogeneity MRI assessment with histograms

By definition, tumours are heterogeneous. They are defined by marked differences in cells, microenvironmental factors (oxygenation levels, pH, VEGF, VPF and TGF-α) metabolism, vasculature, structure and function that in turn translate into heterogeneous drug delivery and therapeutic outcome. Ways to estimate quantitatively tumour heterogeneity can improve drug discovery, treatment planning and therapeutic responses. It is therefore of paramount importance to have reliable and reproducible biomarkers of cancerous lesions' heterogeneity. During the past decade, the number of studies using histogram approaches increased drastically with various magnetic resonance imaging (MRI) techniques (DCE-MRI, DWI, SWI etc.) although information on tumour heterogeneity remains poorly exploited. This fact can be attributed to a poor knowledge of the available metrics and of their specific meaning as well as to the lack of literature references to standardised histogram methods with which surrogate markers of heterogeneity can be compared. This review highlights the current knowledge and critical advances needed to investigate and quantify tumour heterogeneity. The key role of imaging techniques and in particular the key role of MRI for an accurate investigation of tumour heterogeneity is reviewed with a particular emphasis on histogram approaches and derived methods.

MRI measurements of vessel calibre in tumour xenografts: comparison with vascular corrosion casting

Vessel size index (R_v, μm) has been proposed as a quantitative magnetic resonance imaging (MRI) derived imaging biomarker in oncology, for the non-invasive assessment of tumour blood vessel architecture and vascular targeted therapies. Appropriate pre-clinical evaluation of R_v in animal tumour models will improve the interpretation and guide the introduction of the biomarker into clinical studies. The objective of this study was to compare R_v measured in vivo with vessel size measurements from high-resolution X-ray computed tomography (μCT) of vascular corrosion casts measured post mortem from the same tumours, with and without vascular targeted therapy. MRI measurements were first acquired from subcutaneous SW1222 colorectal xenografts in mice following treatment with 0 (n = 6), 30 (n = 6) or 200 mg/kg (n = 3) of the vascular disrupting agent ZD6126. The mice were then immediately infused with a low viscosity resin and, following polymerisation and maceration of surrounding tissues, the resulting tumour vascular casts were dissected and subsequently imaged using an optimised μCT imaging approach. Vessel diameters were not measurable by μCT in the 200 mg/kg group as the high dose of ZD6126 precluded delivery of the resin to the tumour vascular bed. The mean R_v for the three treatment groups was 24, 23 and 23.5 μm respectively; the corresponding μCT measurements from corrosion casts from the 0 and 30 mg/kg cohorts were 25 and 28 μm. The strong association between the in vivo MRI and post mortem μCT values supports the use of R_v as an imaging biomarker in clinical trials of investigational vascular targeted therapies.