Semi-automatic tumor segmentation of rectal cancer based on functional magnetic resonance imaging.

Background and purpose: Tumor delineation is required both for radiotherapy planning and quantitative imaging biomarker purposes. It is a manual, time- and labor-intensive process prone to inter- and intraobserver variations. Semi or fully automatic segmentation could provide better efficiency and consistency. This study aimed to investigate the influence of including and combining functional with anatomical magnetic resonance imaging (MRI) sequences on the quality of automatic segmentations. Materials and methods: T2-weighted (T2w), diffusion weighted, multi-echo T2*-weighted, and contrast enhanced dynamic multi-echo (DME) MR images of eighty-one patients with rectal cancer were used in the analysis. Four classical machine learning algorithms; adaptive boosting (ADA), linear and quadratic discriminant analysis and support vector machines, were trained for automatic segmentation of tumor and normal tissue using different combinations of the MR images as input, followed by semi-automatic morphological post-processing. Manual delineations from two experts served as ground truth. The Sørensen-Dice similarity coefficient (DICE) and mean symmetric surface distance (MSD) were used as performance metric in leave-one-out cross validation. Results: Using T2w images alone, ADA outperformed the other algorithms, yielding a median per patient DICE of 0.67 and MSD of 3.6 mm. The performance improved when functional images were added and was highest for models based on either T2w and DME images (DICE: 0.72, MSD: 2.7 mm) or all four MRI sequences (DICE: 0.72, MSD: 2.5 mm). Conclusion: Machine learning models using functional MRI, in particular DME, have the potential to improve automatic segmentation of rectal cancer relative to models using T2w MRI alone.

A comparison of methods for fully automatic segmentation of tumors and involved nodes in PET/CT of head and neck cancers.

Target volume delineation is a vital but time-consuming and challenging part of radiotherapy, where the goal is to deliver sufficient dose to the target while reducing risks of side effects. For head and neck cancer (HNC) this is complicated by the complex anatomy of the head and neck region and the proximity of target volumes to organs at risk. The purpose of this study was to compare and evaluate conventional PET thresholding methods, six classical machine learning algorithms and a 2D U-Net convolutional neural network (CNN) for automatic gross tumor volume (GTV) segmentation of HNC in PET/CT images. For the latter two approaches the impact of single versus multimodality input on segmentation quality was also assessed. 197 patients were included in the study. The cohort was split into training and test sets (157 and 40 patients, respectively). Five-fold cross-validation was used on the training set for model comparison and selection. Manual GTV delineations represented the ground truth. Tresholding, classical machine learning and CNN segmentation models were ranked separately according to the cross-validation Sørensen-Dice similarity coefficient (Dice). PET thresholding gave a maximum mean Dice of 0.62, whereas classical machine learning resulted in maximum mean Dice scores of 0.24 (CT) and 0.66 (PET; PET/CT). CNN models obtained maximum mean Dice scores of 0.66 (CT), 0.68 (PET) and 0.74 (PET/CT). The difference in cross-validation Dice between multimodality PET/CT and single modality CNN models was significant (p ≤ 0.0001). The top-ranked PET/CT-based CNN model outperformed the best-performing thresholding and classical machine learning models, giving significantly better segmentations in terms of cross-validation and test set Dice, true positive rate, positive predictive value and surface distance-based metrics (p ≤ 0.0001). Thus, deep learning based on multimodality PET/CT input resulted in superior target coverage and less inclusion of surrounding normal tissue.

Hypoxia and perfusion in breast cancer: simultaneous assessment using PET/MR imaging.

OBJECTIVES: Hypoxia is associated with poor prognosis and treatment resistance in breast cancer. However, the temporally variant nature of hypoxia can complicate interpretation of imaging findings. We explored the relationship between hypoxia and vascular function in breast tumours through combined 18F-fluoromisonidazole (18 F-FMISO) PET/MRI, with simultaneous assessment circumventing the effect of temporal variation in hypoxia and perfusion. METHODS: Women with histologically confirmed, primary breast cancer underwent a simultaneous 18F-FMISO-PET/MR examination. Tumour hypoxia was assessed using influx rate constant Ki and hypoxic fractions (%HF), while parameters of vascular function (Ktrans, kep, ve, vp) and cellularity (ADC) were derived from dynamic contrast-enhanced (DCE) and diffusion-weighted (DW)-MRI, respectively. Additional correlates included histological subtype, grade and size. Relationships between imaging variables were assessed using Pearson correlation (r). RESULTS: Twenty-nine women with 32 lesions were assessed. Hypoxic fractions > 1% were observed in 6/32 (19%) cancers, while 18/32 (56%) tumours showed a %HF of zero. The presence of hypoxia in lesions was independent of histological subtype or grade. Mean tumour Ktrans correlated negatively with Ki (r = - 0.38, p = 0.04) and %HF (r = - 0.33, p = 0.04), though parametric maps exhibited intratumoural heterogeneity with hypoxic regions colocalising with both hypo- and hyperperfused areas. No correlation was observed between ADC and DCE-MRI or PET parameters. %HF correlated positively with lesion size (r = 0.63, p = 0.001). CONCLUSION: Hypoxia measured by 18F-FMISO-PET correlated negatively with Ktrans from DCE-MRI, supporting the hypothesis of perfusion-driven hypoxia in breast cancer. Intratumoural hypoxia-perfusion relationships were heterogeneous, suggesting that combined assessment may be needed for disease characterisation, which could be achieved using simultaneous multimodality imaging. KEY POINTS: • At the tumour level, hypoxia measured by 18F-FMISO-PET was negatively correlated with perfusion measured by DCE-MRI, which supports the hypothesis of perfusion-driven hypoxia in breast cancer. • No associations were observed between 18F-FMISO-PET parameters and tumour histology or grade, but tumour hypoxic fractions increased with lesion size. • Intratumoural hypoxia-perfusion relationships were heterogeneous, suggesting that the combined hypoxia-perfusion status of tumours may need to be considered for disease characterisation, which can be achieved via simultaneous multimodality imaging as reported here.

Hyperpolarized 13C MRI of Tumor Metabolism Demonstrates Early Metabolic Response to Neoadjuvant Chemotherapy in Breast Cancer.

Purpose: To compare hyperpolarized carbon 13 (13C) MRI with dynamic contrast material-enhanced (DCE) MRI in the detection of early treatment response in breast cancer. Materials and Methods: In this institutional review board-approved prospective study, a woman with triple-negative breast cancer (age, 49 years) underwent 13C MRI after injection of hyperpolarized [1-carbon 13 {13C}]-pyruvate and DCE MRI at 3 T at baseline and after one cycle of neoadjuvant therapy. The 13C-labeled lactate-to-pyruvate ratio derived from hyperpolarized 13C MRI and the pharmacokinetic parameters transfer constant (K trans) and washout parameter (k ep) derived from DCE MRI were compared before and after treatment. Results: Exchange of the 13C label between injected hyperpolarized [1-13C]-pyruvate and the endogenous lactate pool was observed, catalyzed by the enzyme lactate dehydrogenase. After one cycle of neoadjuvant chemotherapy, a 34% reduction in the 13C-labeled lactate-to-pyruvate ratio resulted in correct identification of the patient as a responder to therapy, which was subsequently confirmed via a complete pathologic response. However, DCE MRI showed an increase in mean K trans (132%) and mean k ep (31%), which could be incorrectly interpreted as a poor response to treatment. Conclusion: Hyperpolarized 13C MRI enabled successful identification of breast cancer response after one cycle of neoadjuvant chemotherapy and may improve response prediction when used in conjunction with multiparametric proton MRI.Published under a CC BY 4.0 license.

Imaging breast cancer using hyperpolarized carbon-13 MRI.

Our purpose is to investigate the feasibility of imaging tumor metabolism in breast cancer patients using 13C magnetic resonance spectroscopic imaging (MRSI) of hyperpolarized 13C label exchange between injected [1-13C]pyruvate and the endogenous tumor lactate pool. Treatment-naïve breast cancer patients were recruited: four triple-negative grade 3 cancers; two invasive ductal carcinomas that were estrogen and progesterone receptor-positive (ER/PR+) and HER2/neu-negative (HER2-), one grade 2 and one grade 3; and one grade 2 ER/PR+ HER2- invasive lobular carcinoma (ILC). Dynamic 13C MRSI was performed following injection of hyperpolarized [1-13C]pyruvate. Expression of lactate dehydrogenase A (LDHA), which catalyzes 13C label exchange between pyruvate and lactate, hypoxia-inducible factor-1 (HIF1α), and the monocarboxylate transporters MCT1 and MCT4 were quantified using immunohistochemistry and RNA sequencing. We have demonstrated the feasibility and safety of hyperpolarized 13C MRI in early breast cancer. Both intertumoral and intratumoral heterogeneity of the hyperpolarized pyruvate and lactate signals were observed. The lactate-to-pyruvate signal ratio (LAC/PYR) ranged from 0.021 to 0.473 across the tumor subtypes (mean ± SD: 0.145 ± 0.164), and a lactate signal was observed in all of the grade 3 tumors. The LAC/PYR was significantly correlated with tumor volume (R = 0.903, P = 0.005) and MCT 1 (R = 0.85, P = 0.032) and HIF1α expression (R = 0.83, P = 0.043). Imaging of hyperpolarized [1-13C]pyruvate metabolism in breast cancer is feasible and demonstrated significant intertumoral and intratumoral metabolic heterogeneity, where lactate labeling correlated with MCT1 expression and hypoxia.

Decoding the Interdependence of Multiparametric Magnetic Resonance Imaging to Reveal Patient Subgroups Correlated with Survivals.

Glioblastoma is highly heterogeneous in microstructure and vasculature, creating various tumor microenvironments among patients, which may lead to different phenotypes. The purpose was to interrogate the interdependence of microstructure and vasculature using perfusion and diffusion imaging and to investigate the utility of this approach in tumor invasiveness assessment. A total of 115 primary glioblastoma patients were prospectively recruited for preoperative magnetic resonance imaging (MRI) and surgery. Apparent diffusion coefficient (ADC) was calculated from diffusion imaging, and relative cerebral blood volume (rCBV) was calculated from perfusion imaging. The empirical copula transform was applied to ADC and rCBV voxels in the contrast-enhancing tumor region to obtain their joint distribution, which was discretized to extract second-order features for an unsupervised hierarchical clustering. The lactate levels of patient subgroups, measured by MR spectroscopy, were compared. Survivals were analyzed using Kaplan-Meier and multivariate Cox regression analyses. The results showed that three patient subgroups were identified by the unsupervised clustering. These subtypes showed no significant differences in clinical characteristics but were significantly different in lactate level and patient survivals. Specifically, the subtype demonstrating high interdependence of ADC and rCBV displayed a higher lactate level than the other two subtypes (P = .016 and P = .044, respectively). Both subtypes of low and high interdependence showed worse progression-free survival than the intermediate (P = .046 and P = .009 respectively). Our results suggest that the interdependence between perfusion and diffusion imaging may be useful in stratifying patients and evaluating tumor invasiveness, providing overall measure of tumor microenvironment using multiparametric MRI.

Characterizing tumor invasiveness of glioblastoma using multiparametric magnetic resonance imaging.

OBJECTIVE: The objective of this study was to characterize the abnormalities revealed by diffusion tensor imaging (DTI) using MR spectroscopy (MRS) and perfusion imaging, and to evaluate the prognostic value of a proposed quantitative measure of tumor invasiveness by combining contrast-enhancing (CE) and DTI abnormalities in patients with glioblastoma. METHODS: Eighty-four patients with glioblastoma were recruited preoperatively. DTI was decomposed into isotropic (p) and anisotropic (q) components. The relative cerebral blood volume (rCBV) was calculated from the dynamic susceptibility contrast imaging. Values of N-acetylaspartate, myoinositol, choline (Cho), lactate (Lac), and glutamate + glutamine (Glx) were measured from multivoxel MRS and normalized as ratios to creatine (Cr). Tumor regions of interest (ROIs) were manually segmented from the CE T1-weighted (CE-ROI) and DTI-q (q-ROI) maps. Perfusion and metabolic characteristics of these ROIs were measured and compared. The relative invasiveness coefficient (RIC) was calculated as a ratio of the characteristic radii of CE-ROI and q-ROI. The prognostic significance of RIC was tested using Kaplan-Meier and multivariate Cox regression analyses. RESULTS: The Cho/Cr, Lac/Cr, and Glx/Cr in q-ROI were significantly higher than CE-ROI (p = 0.004, p = 0.005, and p = 0.007, respectively). CE-ROI had significantly higher rCBV values than q-ROI (p < 0.001). A higher RIC was associated with worse survival in a multivariate overall survival (OS) model (hazard ratio [HR] 1.40, 95% confidence interval [CI] 1.06-1.85, p = 0.016) and progression-free survival (PFS) model (HR 1.55, 95% CI 1.16-2.07, p = 0.003). An RIC cutoff value of 0.89 significantly predicted shorter OS (median 384 vs 605 days, p = 0.002) and PFS (median 244 vs 406 days, p = 0.001). CONCLUSIONS: DTI-q abnormalities displayed higher tumor load and hypoxic signatures compared with CE abnormalities, whereas CE regions potentially represented the tumor proliferation edge. Integrating the extents of invasion visualized by DTI-q and CE images into clinical practice may lead to improved treatment efficacy.

Low perfusion compartments in glioblastoma quantified by advanced magnetic resonance imaging and correlated with patient survival.

BACKGROUND AND PURPOSE: Glioblastoma exhibits profound intratumoral heterogeneity in perfusion. Particularly, low perfusion may induce treatment resistance. Thus, imaging approaches that define low perfusion compartments are crucial for clinical management. MATERIALS AND METHODS: A total of 112 newly diagnosed glioblastoma patients were prospectively recruited for maximal safe resection. The apparent diffusion coefficient (ADC) and relative cerebral blood volume (rCBV) were calculated from diffusion and perfusion imaging, respectively. Based on the overlapping regions of lowest rCBV quartile (rCBVL) with the lowest ADC quartile (ADCL) and highest ADC quartile (ADCH) in each tumor, two low perfusion compartments (ADCH-rCBVL and ADCL-rCBVL) were identified for volumetric analysis. Lactate and macromolecule/lipid levels were determined from multivoxel MR spectroscopic imaging. Progression-free survival (PFS) and overall survival (OS) were analyzed using Kaplan-Meier's and multivariate Cox regression analyses, to evaluate the effects of compartment volume and lactate level on survival. RESULTS: Two compartments displayed higher lactate and macromolecule/lipid levels compared to contralateral normal-appearing white matter (each P < 0.001). The proportion of the ADCL-rCBVL compartment in the contrast-enhancing tumor was associated with a larger infiltration on FLAIR (P < 0.001, rho = 0.42). The minimally invasive phenotype displayed a lower proportion of the ADCL-rCBVL compartment than the localized (P = 0.031) and diffuse phenotypes (not significant). Multivariate Cox regression showed higher lactate level in the ADCL-rCBVL compartment was associated with worsened survival (PFS: HR 2.995, P = 0.047; OS: HR 4.974, P = 0.005). CONCLUSIONS: Our results suggest that the ADCL-rCBVL compartment may potentially indicate a clinically measurable resistant compartment.

Multi-parametric and multi-regional histogram analysis of MRI: modality integration reveals imaging phenotypes of glioblastoma.

OBJECTIVES: Integrating multiple imaging modalities is crucial for MRI data interpretation. The purpose of this study is to determine whether a previously proposed multi-view approach can effectively integrate the histogram features from multi-parametric MRI and whether the selected features can offer incremental prognostic values over clinical variables. METHODS: Eighty newly-diagnosed glioblastoma patients underwent surgery and chemoradiotherapy. Histogram features of diffusion and perfusion imaging were extracted from contrast-enhancing (CE) and non-enhancing (NE) regions independently. An unsupervised patient clustering was performed by the multi-view approach. Kaplan-Meier and Cox proportional hazards regression analyses were performed to evaluate the relevance of patient clustering to survival. The metabolic signatures of patient clusters were compared using multi-voxel spectroscopy analysis. The prognostic values of histogram features were evaluated by survival and ROC curve analyses. RESULTS: Two patient clusters were generated, consisting of 53 and 27 patients respectively. Cluster 2 demonstrated better overall survival (OS) (p = 0.007) and progression-free survival (PFS) (p < 0.001) than Cluster 1. Cluster 2 displayed lower N-acetylaspartate/creatine ratio in NE region (p = 0.040). A higher mean value of anisotropic diffusion in NE region was associated with worse OS (hazard ratio [HR] = 1.40, p = 0.020) and PFS (HR = 1.36, p = 0.031). The seven features selected by this approach showed significantly incremental value in predicting 12-month OS (p = 0.020) and PFS (p = 0.022). CONCLUSIONS: The multi-view clustering method can provide an effective integration of multi-parametric MRI. The histogram features selected may be used as potential prognostic markers. KEY POINTS: • Multi-parametric magnetic resonance imaging captures multi-faceted tumor physiology. • Contrast-enhancing and non-enhancing tumor regions represent different tumor components with distinct clinical relevance. • Multi-view data analysis offers a method which can effectively select and integrate multi-parametric and multi-regional imaging features.

Optoacoustic Imaging Detects Hormone-Related Physiological Changes of Breast Parenchyma.

PURPOSE: Optoacoustic imaging with ultrasound (OPUS) can assess in-vivo perfusion/oxygenation through surrogate measures of oxy, deoxy and total hemoglobin content in tissues. The primary aim of our study was to evaluate the ability of OPUS to detect physiological changes in the breast during the menstrual cycle and to determine qualitative/quantitative metrics of normal parenchymal tissue in pre-/post-menopausal women. The secondary aim was to assess the technique's repeatability. MATERIALS AND METHODS: We performed a prospective ethically approved study in volunteers using OPUS (700, 800 and 850 nm wavelengths) in the proliferative/follicular and secretory phase of the menstrual cycle. Regions of interest (ROIs) were drawn on the most superficial region of fibroglandular tissue and same-day intra-observer repeatability was assessed. We used t-tests to interrogate differences in the OPUS measurements due to hormonal changes and interclass correlation coefficients/Bland-Altman plots to evaluate the repeatability of mean ROI signal intensities. RESULTS: 22 pre-menopausal and 8 post-menopausal volunteers were recruited. 21 participants underwent repeatability examinations. OPUS intensity values were significantly higher (p < 0.0001) at all excitation wavelengths in the secretory compared to the proliferative/follicular phase. Post-menopausal volunteers showed similar optoacoustic values to the proliferative/follicular phase of pre-menopausal volunteers. The repeatability of the technique was comparable to other handheld ultrasound modalities. CONCLUSION: OPUS detects changes in perfusion/vascularity related to the menstrual cycle and menopausal status of breast parenchyma.

Intratumoral Heterogeneity of Glioblastoma Infiltration Revealed by Joint Histogram Analysis of Diffusion Tensor Imaging.

BACKGROUND: Glioblastoma is a heterogeneous disease characterized by its infiltrative growth, rendering complete resection impossible. Diffusion tensor imaging (DTI) shows potential in detecting tumor infiltration by reflecting microstructure disruption. OBJECTIVE: To explore the heterogeneity of glioblastoma infiltration using joint histogram analysis of DTI, to investigate the incremental prognostic value of infiltrative patterns over clinical factors, and to identify specific subregions for targeted therapy. METHODS: A total of 115 primary glioblastoma patients were prospectively recruited for surgery and preoperative magnetic resonance imaging. The joint histograms of decomposed anisotropic and isotropic components of DTI were constructed in both contrast-enhancing and nonenhancing tumor regions. Patient survival was analyzed with joint histogram features and relevant clinical factors. The incremental prognostic values of histogram features were assessed using receiver operating characteristic curve analysis. The correlation between the proportion of diffusion patterns and tumor progression rate was tested using Pearson correlation. RESULTS: We found that joint histogram features were associated with patient survival and improved survival model performance. Specifically, the proportion of nonenhancing tumor subregion with decreased isotropic diffusion and increased anisotropic diffusion was correlated with tumor progression rate (P = .010, r = 0.35), affected progression-free survival (hazard ratio = 1.08, P < .001), and overall survival (hazard ratio = 1.36, P < .001) in multivariate models. CONCLUSION: Joint histogram features of DTI showed incremental prognostic values over clinical factors for glioblastoma patients. The nonenhancing tumor subregion with decreased isotropic diffusion and increased anisotropic diffusion may indicate a more infiltrative habitat and potential treatment target.

In Quest of the Alanine R3 Radical: Multivariate EPR Spectral Analyses of X-Irradiated Alanine in the Solid State.

The amino acid l-α-alanine is the most commonly used material for solid-state electron paramagnetic resonance (EPR) dosimetry, due to the formation of highly stable radicals upon irradiation, with yields proportional to the radiation dose. Two major alanine radical components designated R1 and R2 have previously been uniquely characterized from EPR and electron-nuclear double resonance (ENDOR) studies as well as from quantum chemical calculations. There is also convincing experimental evidence of a third minor radical component R3, and a tentative radical structure has been suggested, even though no well-defined spectral signature has been observed experimentally. In the present study, temperature dependent EPR spectra of X-ray irradiated polycrystalline alanine were analyzed using five multivariate methods in further attempts to understand the composite nature of the alanine dosimeter EPR spectrum. Principal component analysis (PCA), maximum likelihood common factor analysis (MLCFA), independent component analysis (ICA), self-modeling mixture analysis (SMA), and multivariate curve resolution (MCR) were used to extract pure radical spectra and their fractional contributions from the experimental EPR spectra. All methods yielded spectral estimates resembling the established R1 spectrum. Furthermore, SMA and MCR consistently predicted both the established R2 spectrum and the shape of the R3 spectrum. The predicted shape of the R3 spectrum corresponded well with the proposed tentative spectrum derived from spectrum simulations. Thus, results from two independent multivariate data analysis techniques strongly support the previous evidence that three radicals are indeed present in irradiated alanine samples.

Haralick texture features from apparent diffusion coefficient (ADC) MRI images depend on imaging and pre-processing parameters.

In recent years, texture analysis of medical images has become increasingly popular in studies investigating diagnosis, classification and treatment response assessment of cancerous disease. Despite numerous applications in oncology and medical imaging in general, there is no consensus regarding texture analysis workflow, or reporting of parameter settings crucial for replication of results. The aim of this study was to assess how sensitive Haralick texture features of apparent diffusion coefficient (ADC) MR images are to changes in five parameters related to image acquisition and pre-processing: noise, resolution, how the ADC map is constructed, the choice of quantization method, and the number of gray levels in the quantized image. We found that noise, resolution, choice of quantization method and the number of gray levels in the quantized images had a significant influence on most texture features, and that the effect size varied between different features. Different methods for constructing the ADC maps did not have an impact on any texture feature. Based on our results, we recommend using images with similar resolutions and noise levels, using one quantization method, and the same number of gray levels in all quantized images, to make meaningful comparisons of texture feature results between different subjects.

Autodelineation of cervical cancers using multiparametric magnetic resonance imaging and machine learning.

BACKGROUND: Tumour delineation is a challenging, time-consuming and complex part of radiotherapy planning. In this study, an automatic method for delineating locally advanced cervical cancers was developed using a machine learning approach. MATERIALS AND METHODS: A method for tumour segmentation based on image voxel classification using Fisher?s Linear Discriminant Analysis (LDA) was developed. This was applied to magnetic resonance (MR) images of 78 patients with locally advanced cervical cancer. The segmentation was based on multiparametric MRI consisting of T2- weighted (T2w), T1-weighted (T1w) and dynamic contrast-enhanced (DCE) sequences, and included intensity and spatial information from the images. The model was trained and assessed using delineations made by two radiologists. RESULTS: Segmentation based on T2w or T1w images resulted in mean sensitivity and specificity of 94% and 52%, respectively. Including DCE-MR images improved the segmentation model?s performance significantly, giving mean sensitivity and specificity of 85?93%. Comparisons with radiologists? tumour delineations gave Dice similarity coefficients of up to 0.44. CONCLUSION: Voxel classification using a machine learning approach is a flexible and fully automatic method for tumour delineation. Combining all relevant MR image series resulted in high sensitivity and specificity. Moreover, the presented method can be extended to include additional imaging modalities.

Dynamic 2-Deoxy-2-[18F]Fluoro-D-Glucose Positron Emission Tomography for Chemotherapy Response Monitoring of Breast Cancer Xenografts.

PURPOSE: Non-invasive response monitoring can potentially be used to guide therapy selection for breast cancer patients. We employed dynamic 2-deoxy-2-[18F]fluoro-D-glucose positron emission tomography ([18F]FDG PET) to evaluate changes in three breast cancer xenograft lines in mice following three chemotherapy regimens. PROCEDURES: Sixty-six athymic nude mice bearing bilateral breast cancer xenografts (two basal-like and one luminal-like subtype) underwent three 60 min [18F]FDG PET scans. Scans were performed prior to and 3 and 10 days after treatment with doxorubicin, paclitaxel, or carboplatin. Tumor growth was monitored in parallel. A pharmacokinetic compartmental model was fitted to the tumor uptake curves, providing estimates of transfer rates between the vascular, non-metabolized, and metabolized compartments. Early and late standardized uptake values (SUVE and SUVL, respectively); the rate constants k 1, k 2, and k 3, and the intravascular fraction v B were estimated. Changes in tumor volume were used as a response measure. Multivariate partial least-squares regression (PLSR) was used to assess if PET parameters could model tumor response and to identify PET parameters with the largest impact on response. RESULTS: Treatment responders had significantly larger perfusion-related parameters (k 1 and k 2) and lower metabolism-related parameter (k 3) than non-responders 10 days after the start of treatment. These findings were further supported by the PLSR analysis, which showed that k 1 and k 2 at day 10 and changes in k 3 explained most of the variability in response to therapy, whereas SUVL and particularly SUVE were of lesser importance. CONCLUSIONS: Overall, rate parameters related to both tumor perfusion and metabolism were associated with tumor response. Conventional metrics of [18F]FDG uptake such as SUVE and SUVL apparently had little relation to tumor response, thus necessitating full dynamic scanning and pharmacokinetic analysis for optimal evaluation of chemotherapy-induced changes in breast cancers.

Cluster analysis of dynamic contrast enhanced MRI reveals tumor subregions related to locoregional relapse for cervical cancer patients.

BACKGROUND: Solid tumors are known to be spatially heterogeneous. Detection of treatment-resistant tumor regions can improve clinical outcome, by enabling implementation of strategies targeting such regions. In this study, K-means clustering was used to group voxels in dynamic contrast enhanced magnetic resonance images (DCE-MRI) of cervical cancers. The aim was to identify clusters reflecting treatment resistance that could be used for targeted radiotherapy with a dose-painting approach. MATERIAL AND METHODS: Eighty-one patients with locally advanced cervical cancer underwent DCE-MRI prior to chemoradiotherapy. The resulting image time series were fitted to two pharmacokinetic models, the Tofts model (yielding parameters Ktrans and νe) and the Brix model (ABrix, kep and kel). K-means clustering was used to group similar voxels based on either the pharmacokinetic parameter maps or the relative signal increase (RSI) time series. The associations between voxel clusters and treatment outcome (measured as locoregional control) were evaluated using the volume fraction or the spatial distribution of each cluster. RESULTS: One voxel cluster based on the RSI time series was significantly related to locoregional control (adjusted p-value 0.048). This cluster consisted of low-enhancing voxels. We found that tumors with poor prognosis had this RSI-based cluster gathered into few patches, making this cluster a potential candidate for targeted radiotherapy. None of the voxels clusters based on Tofts or Brix parameter maps were significantly related to treatment outcome. CONCLUSION: We identified one group of tumor voxels significantly associated with locoregional relapse that could potentially be used for dose painting. This tumor voxel cluster was identified using the raw MRI time series rather than the pharmacokinetic maps.

Classification of dynamic contrast enhanced MR images of cervical cancers using texture analysis and support vector machines.

Dynamic contrast enhanced MRI (DCE-MRI) provides insight into the vascular properties of tissue. Pharmacokinetic models may be fitted to DCE-MRI uptake patterns, enabling biologically relevant interpretations. The aim of our study was to determine whether treatment outcome for 81 patients with locally advanced cervical cancer could be predicted from parameters of the Brix pharmacokinetic model derived from pre-chemoradiotherapy DCE-MRI. First-order statistical features of the Brix parameters were used. In addition, texture analysis of Brix parameter maps was done by constructing gray level co-occurrence matrices (GLCM) from the maps. Clinical factors and first- and second-order features were used as explanatory variables for support vector machine (SVM) classification, with treatment outcome as response. Classification models were validated using leave-one-out cross-model validation. A random value permutation test was used to evaluate model significance. Features derived from first-order statistics could not discriminate between cured and relapsed patients (specificity 0%-20%, p-values close to unity). However, second-order GLCM features could significantly predict treatment outcome with accuracies (~70%) similar to the clinical factors tumor volume and stage (69%). The results indicate that the spatial relations within the tumor, quantified by texture features, were more suitable for outcome prediction than first-order features.