The prognostic value of dynamic contrast-enhanced MRI contrast agent transfer constant Ktrans in cervical cancer is explained by plasma flow rather than vessel permeability.

BACKGROUND: The microvascular contrast agent transfer constant Ktrans has shown prognostic value in cervical cancer patients treated with chemoradiotherapy. This study aims to determine whether this is explained by the contribution to Ktrans of plasma flow (Fp), vessel permeability surface-area product (PS), or a combination of both. METHODS: Pre-treatment dynamic contrast-enhanced MRI (DCE-MRI) data from 36 patients were analysed using the two-compartment exchange model. Estimates of Fp, PS, Ktrans, and fractional plasma and interstitial volumes (vp and ve) were made and used in univariate and multivariate survival analyses adjusting for clinicopathologic variables tumour stage, nodal status, histological subtype, patient age, tumour volume, and treatment type (chemoradiotherapy vs radiotherapy alone). RESULTS: In univariate analyses, Fp (HR=0.25, P=0.0095) and Ktrans (HR=0.20, P=0.032) were significantly associated with disease-free survival while PS, vp and ve were not. In multivariate analyses adjusting for clinicopathologic variables, Fp and Ktrans significantly increased the accuracy of survival predictions (P=0.0089). CONCLUSIONS: The prognostic value of Ktrans in cervical cancer patients treated with chemoradiotherapy is explained by microvascular plasma flow (Fp) rather than vessel permeability surface-area product (PS).

Improved accuracy and precision of tracer kinetic parameters by joint fitting to variable flip angle and dynamic contrast enhanced MRI data.

PURPOSE: To improve the accuracy and precision of tracer kinetic model parameter estimates for use in dynamic contrast enhanced (DCE) MRI studies of solid tumors. THEORY: Quantitative DCE-MRI requires an estimate of precontrast T1 , which is obtained prior to fitting a tracer kinetic model. As T1 mapping and tracer kinetic signal models are both a function of precontrast T1 it was hypothesized that its joint estimation would improve the accuracy and precision of both precontrast T1 and tracer kinetic model parameters. METHODS: Accuracy and/or precision of two-compartment exchange model (2CXM) parameters were evaluated for standard and joint fitting methods in well-controlled synthetic data and for 36 bladder cancer patients. Methods were compared under a number of experimental conditions. RESULTS: In synthetic data, joint estimation led to statistically significant improvements in the accuracy of estimated parameters in 30 of 42 conditions (improvements between 1.8% and 49%). Reduced accuracy was observed in 7 of the remaining 12 conditions. Significant improvements in precision were observed in 35 of 42 conditions (between 4.7% and 50%). In clinical data, significant improvements in precision were observed in 18 of 21 conditions (between 4.6% and 38%). CONCLUSION: Accuracy and precision of DCE-MRI parameter estimates are improved when signal models are fit jointly rather than sequentially. Magn Reson Med 76:1270-1281, 2016. © 2015 Wiley Periodicals, Inc.

Automated prognostic pattern detection shows favourable diffuse pattern of FOXP3(+) Tregs in follicular lymphoma.

BACKGROUND: Histopathological prognostication relies on morphological pattern recognition, but as numbers of biomarkers increase, human prognostic pattern-recognition ability decreases. Follicular lymphoma (FL) has a variable outcome, partly determined by FOXP3 Tregs. We have developed an automated method, hypothesised interaction distribution (HID) analysis, to analyse spatial patterns of multiple biomarkers which we have applied to tumour-infiltrating lymphocytes in FL. METHODS: A tissue microarray of 40 patient samples was used in triplex immunohistochemistry for FOXP3, CD3 and CD69, and multispectral imaging used to determine the numbers and locations of CD3(+), FOXP3/CD3(+) and CD69/CD3(+) T cells. HID analysis was used to identify associations between cellular pattern and outcome. RESULTS: Higher numbers of CD3(+) (P=0.0001), FOXP3/CD3(+) (P=0.0031) and CD69/CD3(+) (P=0.0006) cells were favourable. Cross-validated HID analysis of cell pattern identified patient subgroups with statistically significantly different survival (35.5 vs 142 months, P=0.00255), a more diffuse pattern associated with favourable outcome and an aggregated pattern with unfavourable outcome. CONCLUSIONS: A diffuse pattern of FOXP3 and CD69 positivity was favourable, demonstrating ability of HID analysis to automatically identify prognostic cellular patterns. It is applicable to large numbers of biomarkers, representing an unsupervised, automated method for identification of undiscovered prognostic cellular patterns in cancer tissue samples.

Indexed distribution analysis for improved significance testing of spatially heterogeneous parameter maps: application to dynamic contrast-enhanced MRI biomarkers.

PURPOSE: To develop significance testing methodology applicable to spatially heterogeneous parametric maps of biophysical and physiological measurements arising from imaging studies. THEORY: Heterogeneity can confound statistical analyses. Indexed distribution analysis (IDA) transforms a reference distribution, establishing correspondences across parameter maps to which significance tests are applied. METHODS: Well-controlled simulated and clinical K(trans) data from a dynamic contrast-enhanced magnetic resonance imaging study of bevacizumab were analyzed using conventional significance tests of parameter averages, histogram analysis, and IDA. Repeated pretreatment scans provided negative control; a post treatment scan provided positive control. RESULTS: Histogram analysis was insensitive to simulated and known effects. Simulation: conventional analysis identified treatment effect (P ≈ 5 × 10(-4)) and direction, but underestimated magnitude (relative error 67-81%); IDA identified treatment effect (P = 0.001), magnitude, direction, and spatial extent (100% accuracy). Bevacizumab: conventional analysis was sensitive to treatment effect (P = 0.01; 95% confidence interval on K(trans) decrease: 23-37%); IDA was sensitive to treatment effect (P < 0.05; K(trans) decrease approximately 25%), inferred its spatial extent to be 94-96%, and inferred that K(trans) decrease is independent of baseline value, an inference that conventional and histogram analyses cannot make. CONCLUSIONS: In the presence of heterogeneity, IDA can accurately infer the magnitude, direction, and spatial extent of between samples of parametric maps, which can be visualized spatially with respect to the original parameter maps.

Quantifying spatial heterogeneity in dynamic contrast-enhanced MRI parameter maps.

Dynamic contrast-enhanced MRI is becoming a standard tool for imaging-based trials of anti-vascular/angiogenic agents in cancer. So far, however, biomarkers derived from DCE-MRI parameter maps have largely neglected the fact that the maps have spatial structure and instead focussed on distributional summary statistics. Such statistics-e.g., biomarkers based on median values-neglect the spatial arrangement of parameters, which may carry important diagnostic and prognostic information. This article describes two types of heterogeneity biomarker that are sensitive to both parameter values and their spatial arrangement. Methods based on Rényi fractal dimensions and geometrical properties are developed, both of which attempt to describe the complexity of DCE-MRI parameter maps. Experiments using simulated data show that the proposed biomarkers are sensitive to changes that distribution-based summary statistics cannot detect and demonstrate that heterogeneity biomarkers could be applied in the drug trial setting. An experiment using 23 DCE-MRI parameter maps of gliomas-a class of tumour that is graded on the basis of heterogeneity-shows that the proposed heterogeneity biomarkers are able to differentiate between low- and high-grade tumours.

Tumour enhancing fraction (EnF) in glioma: relationship to tumour grade.

The aim of this research was to determine whether the proportion of a tumour that enhances (enhancing fraction, EnF) and changes in EnF with enhancement threshold differ between low and high grade glioma. Forty-four patients (45 gliomas comprising 16 grade II, 5 grade III and 24 grade IV) were studied. Imaging included pre- and post-contrast-enhanced T(1)-weighted sequences and T(1)-weighted DCE-MRI. Thresholded enhancement maps were generated for each tumour by using a range of values of the initial area under the contrast concentration curve (IAUC). A plot of EnF versus threshold value was generated. We examined the relationship between tumour grade and enhancement metrics including: EnF (threshold IAUC > 0 mMol s), EnF (threshold IAUC > 2.5 mMol s), initial slope of the EnF/threshold curve (partial differentialEnF), IAUC, and two previously described signal-intensity-based metrics. EnF, defined as the proportion of tumour showing any enhancement (threshold IAUC > 0 mMol s), showed no difference between low and high grade glioma. All other measures demonstrated significant differences between grade II and IV, and low (grade II) and high grade (grades III/ IV) gliomas (p < 0.01). Two measures, partial differentialEnF and Pronin's measure of enhancement, showed differences between grade III and IV (p < 0.05). No measure separated grade II from III. Metrics which describe the enhancing fraction and its variation with enhancement threshold partial differentialEnF show considerably different behaviour in low and high grade tumours. These observations suggest that these metrics may provide important biological information concerning tumour biology and therapeutic responses and encourage further research to characterise and validate these novel biomarkers.

Glandular function in Sjögren syndrome: assessment with dynamic contrast-enhanced MR imaging and tracer kinetic modeling--initial experience.

PURPOSE: To prospectively use dynamic contrast material-enhanced magnetic resonance (MR) imaging and a tracer kinetic model to compare parotid gland microvascular characteristics in patients who have Sjögren syndrome (SS) with those in healthy volunteers. MATERIALS AND METHODS: The local research ethics committee approved the study, and written informed consent was obtained from all participants. Twenty-one patients (19 women, two men; age range, 31-73 years) with a diagnosis of SS and 11 healthy volunteers (10 women, one man; age range, 41-68 years) underwent three-dimensional T1-weighted dynamic contrast-enhanced MR imaging of the parotid gland at 1.5 T. A voxel-wise tracer kinetic model and a model-free analysis were applied to the dynamic MR data. Parameter medians and standard deviations were computed to summarize gland microvascular characteristics and gland heterogeneity, respectively. Differences were investigated by using multivariate analysis of variance, t, or U tests. Further investigation was performed by using linear discriminant and receiver operating characteristic analyses. RESULTS: Compared with the healthy volunteers, the patients with SS had highly significant elevations (P << .001) in the model-free parameter initial area under the curve and in tracer kinetic model parameters, including transcapillary contrast agent transfer constant (P < .001) and extracellular extravascular volume (P < .001). Gland heterogeneity was significantly greater (P < .001) in the patients with SS. Parameter medians and standard deviations enabled excellent differentiation (areas under receiver operating characteristic curve, 0.96 and 1.00, respectively) between the patients with SS and the healthy volunteers. CONCLUSION: Dynamic contrast-enhanced MR imaging has the potential to be used in clinical settings to quantify microvascular function in SS and to differentiate between patients with and those without SS.

Enhancing fraction predicts clinical outcome following first-line chemotherapy in patients with epithelial ovarian carcinoma.

PURPOSE: To define a simple radiologic biomarker of prognosis in patients with advanced epithelial ovarian carcinoma on first-line chemotherapy. EXPERIMENTAL DESIGN: Twenty-seven patients receiving platinum-based chemotherapy with >2 cm residual disease [International Federation of Gynecology and Obstetrics (FIGO) stages IIIC or IV] after surgery were identified. The proportion of enhancing tumor tissue--the enhancing fraction--was calculated on pre-chemotherapy computed tomography scans at four Hounsfield unit (HU) thresholds and assessed for correlation with CA125 response, Response Evaluation Criteria in Solid Tumors (RECIST) radiologic response, and time to progression. Discriminative power was assessed by leave-one-out discriminant analysis. RESULTS: Pre-chemotherapy residual tumor volume did not correlate with clinical outcome. Pre-chemotherapy enhancing fraction at all thresholds significantly correlated with CA125 response (P < 0.001, rho = 0.553 for 50 HU; P < 0.001, rho = 0.565 for 60 HU; P < 0.001, rho = 0.553 for 70 HU; P = 0.001, rho = 0.516 for 80 HU). Significant correlations were also shown for radiologic response at all thresholds. Enhancing fraction predicted CA125 response with 81.9% to 86.4% specificity and Response Evaluation Criteria in Solid Tumors response with 74.9% to 76.8% specificity at 95% sensitivity (dependent on threshold). Enhancing fraction correlated with time to progression at the 60 HU (P = 0.045, rho = 0.336) and 70 HU (P = 0.042; rho = 0.340) thresholds. CONCLUSION: Pre-chemotherapy enhancing fraction is a simple quantitative radiologic measure. Further evaluation in larger trials is required to confirm the potential of enhancing fraction as a predictive factor, particularly for patients who may benefit from the addition of antiangiogenic therapy.

Imaging intratumor heterogeneity: role in therapy response, resistance, and clinical outcome.

Tumors exhibit genomic and phenotypic heterogeneity, which has prognostic significance and may influence response to therapy. Imaging can quantify the spatial variation in architecture and function of individual tumors through quantifying basic biophysical parameters such as CT density or MRI signal relaxation rate; through measurements of blood flow, hypoxia, metabolism, cell death, and other phenotypic features; and through mapping the spatial distribution of biochemical pathways and cell signaling networks using PET, MRI, and other emerging molecular imaging techniques. These methods can establish whether one tumor is more or less heterogeneous than another and can identify subregions with differing biology. In this article, we review the image analysis methods currently used to quantify spatial heterogeneity within tumors. We discuss how analysis of intratumor heterogeneity can provide benefit over more simple biomarkers such as tumor size and average function. We consider how imaging methods can be integrated with genomic and pathology data, instead of being developed in isolation. Finally, we identify the challenges that must be overcome before measurements of intratumoral heterogeneity can be used routinely to guide patient care.

A statistical framework for analyzing hypothesized interactions between cells imaged using multispectral microscopy and multiple immunohistochemical markers.

BACKGROUND: Multispectral microscopy and multiple staining can be used to identify cells with distinct immunohistochemical (IHC) characteristics. We present here a method called hypothesized interaction distribution (HID) analysis for characterizing the statistical distribution of pair-wise spatial relationships between cells with particular IHC characteristics and apply it to clinical data. MATERIALS AND METHODS: We retrospectively analyzed data from a study of 26 follicular lymphoma patients in which sections were stained for CD20 and YY1. HID analysis, using leave-one-out validation, was used to assign patients to one of two groups. We tested the null hypothesis of no difference in Kaplan-Meier survival curves between the groups. RESULTS: Shannon entropy of HIDs assigned patients to groups that had significantly different survival curves (median survival was 7.70 versus 110 months, P = 0.00750). Hypothesized interactions between pairs of cells positive for both CD20 and YY1 were associated with poor survival. CONCLUSIONS: HID analysis provides quantitative inferences about possible interactions between spatially proximal cells with particular IHC characteristics. In follicular lymphoma, HID analysis was able to distinguish between patients with poor versus good survival, and it may have diagnostic and prognostic utility in this and other diseases.

Quantifying antivascular effects of monoclonal antibodies to vascular endothelial growth factor: insights from imaging.

PURPOSE: Little is known concerning the onset, duration, and magnitude of direct therapeutic effects of anti-vascular endothelial growth factor (VEGF) therapies. Such knowledge would help guide the rational development of targeted therapeutics from bench to bedside and optimize use of imaging technologies that quantify tumor function in early-phase clinical trials. EXPERIMENTAL DESIGN: Preclinical studies were done using ex vivo microcomputed tomography and in vivo ultrasound imaging to characterize tumor vasculature in a human HM-7 colorectal xenograft model treated with the anti-VEGF antibody G6-31. Clinical evaluation was by quantitative magnetic resonance imaging in 10 patients with metastatic colorectal cancer treated with bevacizumab. RESULTS: Microcomputed tomography experiments showed reduction in perfused vessels within 24 to 48 h of G6-31 drug administration (P

Organ-specific effects of oxygen and carbogen gas inhalation on tissue longitudinal relaxation times.

Molecular oxygen has been previously shown to shorten longitudinal relaxation time (T1) in the spleen and renal cortex, but not in the liver or fat. In this study, the magnitude and temporal evolution of this effect were investigated. Medical air, oxygen, and carbogen (95% oxygen/5% CO2) were administered sequentially in 16 healthy volunteers. T1 maps were acquired using spoiled gradient echo sequences (TR=3.5 ms, TE=0.9 ms, alpha=2 degrees/8 degrees/17 degrees) with six acquisitions on air, 12 on oxygen, 12 on carbogen, and six to 12 back on air. Mean T1 values and change in relaxation rate were compared between each phase of gas inhalation in the liver, spleen, skeletal muscle, renal cortex, and fat by one-way analysis of variance. Oxygen-induced T1-shortening occurred in the liver in fasted subjects (P<0.001) but not in non-fasted subjects (P=0.244). T1-shortening in spleen and renal cortex (both P<0.001) were greater than previously reported. Carbogen induced conflicting responses in different organs, suggesting a complex relationship with organ vasculature. Shortening of tissue T1 by oxygen is more pronounced and more complex than previously recognized. The effect may be useful as a biomarker of arterial flow and oxygen delivery to vascular beds.