Interstitial fluid pressure correlates with intravoxel incoherent motion imaging metrics in a mouse mammary carcinoma model

The impact of tumor microenvironment: unraveling the role of physical cues in breast cancer progression

Metastasis accounts for the vast majority of breast cancer-related fatalities. Although the contribution of genetic and epigenetic modifications to breast cancer progression has been widely acknowledged, emerging evidence underscores the pivotal role of physical stimuli in driving breast cancer metastasis. In this review, we summarize the changes in the mechanics of the breast cancer microenvironment and describe the various forces that impact migrating and circulating tumor cells throughout the metastatic process. We also discuss the mechanosensing and mechanotransducing molecules responsible for promoting the malignant phenotype in breast cancer cells. Gaining a comprehensive understanding of the mechanobiology of breast cancer carries substantial potential to propel progress in prognosis, diagnosis, and patient treatment.

Hydraulic conductivity of human cancer tissue: A hybrid study

Background

Elevated tumor tissue interstitial fluid pressure (IFP) is an adverse biomechanical biomarker that predicts poor therapy response and an aggressive phenotype. Advances in functional imaging have opened the prospect of measuring IFP non-invasively. Image-based estimation of the IFP requires knowledge of the tissue hydraulic conductivity (K), a measure for the ease of bulk flow through the interstitium. However, data on the magnitude of K in human cancer tissue are not available.

Methods

We measured the hydraulic conductivity of tumor tissue using modified Ussing chambers in surgical resection specimens. The effect of the tumor microenvironment (TME) on K was investigated by quantifying the collagen content, cell density, and fibroblast density of the tested samples using quantitative immune histochemistry. Also, we developed a computational fluid dynamics (CFD) model to evaluate the role of K on interstitial fluid flow and drug transport in solid tumors.

Results

The results show that the hydraulic conductivity of human tumor tissues is very limited, ranging from approximately 10-15 to 10-14 m2/Pa∙s. Moreover, K values varied significantly between tumor types and between different samples from the same tumor. A significant inverse correlation was found between collagen fiber density and hydraulic conductivity values. However, no correlation was detected between K and cancer cell or fibroblast densities. The computational model demonstrated the impact of K on the interstitial fluid flow and the drug concentration profile: higher K values led to a lower IFP and deeper drug penetration.

Conclusions

Human tumor tissue is characterized by a very limited hydraulic conductivity, representing a barrier to effective drug transport. The results of this study can inform the development of realistic computational models, facilitate non-invasive IFP estimation, and contribute to stromal targeting anticancer therapies.

Interstitial fluid pressure as an emerging biomarker in solid tumors

The physical microenvironment of cancer is characterized by elevated stiffness and tissue pressure, the main component of which is the interstitial fluid pressure (IFP). Elevated IFP is an established negative predictive and prognostic parameter, directly affecting malignant behavior and therapy response. As such, measurement of the IFP would allow to develop strategies aimed at engineering the physical microenvironment of cancer. Traditionally, IFP measurement required the use of invasive methods. Recent progress in dynamic and functional imaging methods such as dynamic contrast enhanced (DCE) magnetic resonance imaging and elastography, combined with numerical models and simulation, allows to comprehensively assess the biomechanical landscape of cancer, and may help to overcome physical barriers to drug delivery and immune cell infiltration. Here, we provide a comprehensive overview of the origin of elevated IFP, and its role in the malignant phenotype. Also, we review the methods used to measure IFP using invasive and imaging based methods, and highlight remaining obstacles and potential areas of progress in order to implement IFP measurement in clinical practice.

Non-Invasive Monitoring of Increased Fibrotic Tissue and Hyaluronan Deposition in the Tumor Microenvironment in the Advanced Stages of Pancreatic Ductal Adenocarcinoma

Simple Summary

Pancreatic ductal adenocarcinoma (PDAC) is a deadly disease with a poor prognosis. A better understanding of the tumor microenvironment may help better treat the disease. Magnetic resonance imaging may be a great tool for monitoring the tumor microenvironment at different stages of tumor evolution. Here, we used multi-parametric magnetic resonance imaging techniques to monitor underlying pathophysiologic processes during the advanced stages of tumor development and correlated with histologic measurements.

Abstract

Pancreatic ductal adenocarcinomas are characterized by a complex and robust tumor microenvironment (TME) consisting of fibrotic tissue, excessive levels of hyaluronan (HA), and immune cells. We utilized quantitative multi-parametric magnetic resonance imaging (mp-MRI) methods at 14 Tesla in a genetically engineered KPC (Kras^LSL-G12D/+, Trp53^LSL-R172H/+, Cre) mouse model to assess the complex TME in advanced stages of tumor development. The whole tumor, excluding cystic areas, was selected as the region of interest for data analysis and subsequent statistical analysis. Pearson correlation was used for statistical inference. There was a significant correlation between tumor volume and T2 (r = −0.66), magnetization transfer ratio (MTR) (r = 0.60), apparent diffusion coefficient (ADC) (r = 0.48), and Glycosaminoglycan-chemical exchange saturation transfer (GagCEST) (r = 0.51). A subset of mice was randomly selected for histological analysis. There were positive correlations between tumor volume and fibrosis (0.92), and HA (r = 0.76); GagCEST and HA (r = 0.81); and MTR and CD31 (r = 0.48). We found a negative correlation between ADC low-b (perfusion) and Ki67 (r = −0.82). Strong correlations between mp-MRI and histology results suggest that mp-MRI can be used as a non-invasive tool to monitor the tumor microenvironment.

Comparison of different ROI analysis methods for liver lesion characterization with simplified intravoxel incoherent motion (IVIM)

This study investigated the impact of different ROI placement and analysis methods on the diagnostic performance of simplified IVIM-DWI for differentiating liver lesions. 1.5/3.0-T DWI data from a respiratory-gated MRI sequence (b = 0, 50, 250, 800 s/mm²) were analyzed in patients with malignant (n = 74/54) and benign (n = 35/19) lesions. Apparent diffusion coefficient ADC = ADC(0,800) and IVIM parameters D₁' = ADC(50,800), D₂' = ADC(250,800), f₁' = f(0,50,800), f₂' = f(0,250,800), and D*' = D*(0,50,250,800) were calculated voxel-wise. For each lesion, a representative 2D-ROI, a 3D-ROI whole lesion, and a 3D-ROI from "good" slices were placed, including and excluding centrally deviating areas (CDA) if present, and analyzed with various histogram metrics. The diagnostic performance of 2D- and 3D-ROIs was not significantly different; e.g. AUC (ADC/D₁'/f₁') were 0.958/0.902/0.622 for 2D- and 0.942/0.892/0.712 for whole lesion 3D-ROIs excluding CDA at 1.5 T (p > 0.05). For 2D- and 3D-ROIs, AUC (ADC/D₁'/D₂') were significantly higher, when CDA were excluded. With CDA included, AUC (ADC/D₁'/D₂'/f₁'/D*') improved when low percentiles were used instead of averages, and was then comparable to the results of average ROI analysis excluding CDA. For lesion differentiation the use of a representative 2D-ROI is sufficient. CDA should be excluded from ROIs by hand or automatically using low percentiles of diffusion coefficients.

Understanding diffusion-weighted MRI analysis: Repeatability and performance of diffusion models in a benign breast lesion cohort

Diffusion-weighted MRI (DWI) is an important tool for oncology research, with great clinical potential for the classification and monitoring of breast lesions. The utility of parameters derived from DWI, however, is influenced by specific analysis choices. The purpose of this study was to critically evaluate repeatability and curve-fitting performance of common DWI signal representations, for a prospective cohort of patients with benign breast lesions. Twenty informed, consented patients with confirmed benign breast lesions underwent repeated DWI (3 T) using: sagittal single-shot spin-echo echo planar imaging, bipolar encoding, TR/TE: 11,600/86 ms, FOV: 180 x 180 mm, matrix: 90 x 90, slices: 60 x 2.5 mm, iPAT: GRAPPA 2, fat suppression, and 13 b-values: 0-700 s/mm² . A phase-reversed scan (b = 0 s/mm² ) was acquired for distortion correction. Voxel-wise repeat-measures coefficients of variation (CoVs) were derived for monoexponential (apparent diffusion coefficient [ADC]), biexponential (intravoxel incoherent motion: f, D, D*) and stretched exponential (α, DDC) across the parameter histograms for lesion regions of interest (ROIs). Goodness-of-fit for each representation was assessed by Bayesian information criterion. The volume of interest (VOI) definition was repeatable (CoV 13.9%). Within lesions, and across both visits and the cohort, there was no dominant best-fit model, with all representations giving the best fit for a fraction of the voxels. Diffusivity measures from the signal representations (ADC, D, DDC) all showed good repeatability (CoV < 10%), whereas parameters associated with pseudodiffusion (f, D*) performed poorly (CoV > 50%). The stretching exponent α was repeatable (CoV < 12%). This pattern of repeatability was consistent over the central part of the parameter percentiles. Assumptions often made in diffusion studies about analysis choices will influence the detectability of changes, potentially obscuring useful information. No single signal representation prevails within or across lesions, or across repeated visits; parameter robustness is therefore a critical consideration. Our results suggest that stretched exponential representation is more repeatable than biexponential, with pseudodiffusion parameters unlikely to provide clinically useful biomarkers.

Statistical Evaluation of Different Mathematical Models for Diffusion Weighted Imaging of Prostate Cancer Xenografts in Mice

Purpose

To evaluate fitting quality and repeatability of four mathematical models for diffusion weighted imaging (DWI) during tumor progression in mouse xenograft model of prostate cancer.

Methods

Human prostate cancer cells (PC-3) were implanted subcutaneously in right hind limbs of 11 immunodeficient mice. Tumor growth was followed by weekly DWI examinations using a 7T MR scanner. Additional DWI examination was performed after repositioning following the fourth DWI examination to evaluate short term repeatability. DWI was performed using 15 and 12 b-values in the ranges of 0-500 and 0-2000 s/mm², respectively. Corrected Akaike information criteria and F-ratio were used to evaluate fitting quality of each model (mono-exponential, stretched exponential, kurtosis, and bi-exponential).

Results

Significant changes were observed in DWI data during the tumor growth, indicated by ADC_m, ADC_s, and ADC_k. Similar results were obtained using low as well as high b-values. No marked changes in model preference were present between the weeks 1−4. The parameters of the mono-exponential, stretched exponential, and kurtosis models had smaller confidence interval and coefficient of repeatability values than the parameters of the bi-exponential model.

Conclusion

Stretched exponential and kurtosis models showed better fit to DWI data than the mono-exponential model and presented with good repeatability.

Evaluation of metronomic chemotherapy response using diffusion and dynamic contrast-enhanced MRI

PURPOSE

To investigate the feasibility of using diffusion MRI (dMRI) and dynamic contrast-enhanced (DCE) MRI to evaluate the treatment response of metronomic chemotherapy (MCT) in the 4T1 mammary tumor model of locally advanced breast cancer.

METHODS

Twelve Balb/c mice with metastatic breast cancer were divided into treated and untreated (control) groups. The treated group (n = 6) received five treatments of anti-metabolite agent 5-Fluorouracil (5FU) in the span of two weeks. dMRI and DCE-MRI were acquired for both treated and control groups before and after MCT. Immunohistochemically staining and measurements were performed after the post-MRI measurements for comparison.

RESULTS

The control mice had significantly (p<0.005) larger tumors than the MCT treated mice. The DCE-MRI analysis showed a decrease in contrast enhancement for the control group, whereas the MCT mice had a more stable enhancement between the pre-chemo and post-chemo time points. This confirms the antiangiogenic effects of 5FU treatment. Comparing amplitude of enhancement revealed a significantly (p<0.05) higher enhancement in the MCT tumors than in the controls. Moreover, the MCT uptake rate was significantly (p<0.001) slower than the controls. dMRI analysis showed the MCT ADC values were significantly larger than the control group at the post-scan time point.

CONCLUSION

dMRI and DCE-MRI can be used as potential biomarkers for assessing the treatment response of MCT. The MRI and pathology observations suggested that in addition to the cytotoxic effect of cell kills, the MCT with a cytotoxic drug, 5FU, induced changes in the tumor vasculature similar to the anti-angiogenic effect.

Magnetic resonance imaging biomarkers for pulsed focused ultrasound treatment of pancreatic ductal adenocarcinoma

BACKGROUND

The robust fibroinflammatory stroma characteristic of pancreatic ductal adenocarcinoma (PDA) impedes effective drug delivery. Pulsed focused ultrasound (pFUS) can disrupt this stroma and has improved survival in an early clinical trial. Non-invasive methods to characterize pFUS treatment effects are desirable for advancement of this promising treatment modality in larger clinical trials.

AIM

To identify promising, non-invasive pre-clinical imaging methods to characterize acute pFUS treatment effects for in vivo models of PDA.

METHODS

We utilized quantitative magnetic resonance imaging methods at 14 tesla in three mouse models of PDA (subcutaneous, orthotopic and transgenic - Kras^LSL-G12D/+, Trp53^LSL-R172H/+, Cre or “KPC”) to assess immediate tumor response to pFUS treatment (VIFU 2000 Alpinion Medical Systems; 475 W peak electric power, 1 ms pulse duration, 1 Hz, duty cycle 0.1%) vs sham therapy, and correlated our results with histochemical data. These pFUS treatment parameters were previously shown to enhance tumor permeability to chemotherapeutics. T1 and T2 relaxation maps, high (126, 180, 234, 340, 549) vs low (7, 47, 81) b-value apparent diffusion coefficient (ADC) maps, magnetization transfer ratio (MTR) maps, and chemical exchange saturation transfer (CEST) maps for the amide proton spectrum (3.5 parts per million or “ppm”) and the glycosaminoglycan spectrum (0.5-1.5 ppm) were generated and analyzed pre-treatment, and immediately post-treatment, using ImageJ. Animals were sacrificed immediately following post-treatment imaging. The whole-tumor was selected as the region of interest for data analysis and subsequent statistical analysis. T-tests and Pearson correlation were used for statistical inference.

RESULTS

Mean high-b value ADC measurements increased significantly with pFUS treatment for all models. Mean glycosaminoglycan CEST and T2 measurements decreased significantly post-treatment for the KPC group. Mean MTR and amide CEST values increased significantly for the KPC group. Hyaluronic acid focal intensities in the treated regions were significantly lower following pFUS treatment for all animal models. The magnetic resonance imaging changes observed acutely following pFUS therapy likely reflect: (1) Sequelae of variable degrees of microcapillary hemorrhage (T1, MTR and amide CEST); (2) Lower PDA glycosaminoglycan content and associated water content (glycosaminoglycan CEST, T2 and hyaluronic acid focal intensity); and (3) Improved tumor diffusivity (ADC) post pFUS treatment.

CONCLUSION

T2, glycosaminoglycan CEST, and ADC maps may provide reliable quantitation of acute pFUS treatment effects for patients with PDA.

Intravoxel Incoherent Motion of Colon Cancer Liver Metastases for the Assessment of Response to Antiangiogenic Treatment: Results from a Pilot Study

OBJECTIVE

This study was aimed at evaluating the intravoxel incoherent motion (IVIM) parameter alterations of liver metastases of colorectal carcinoma (CRC) during antiangiogenic bevacizumab combination therapy.

METHODS

Twenty-five patients with CRC liver metastases treated with bevacizumab in combination with FOLFOX-or-FOLFIRI protocols were enrolled in the study. MRI was performed using a 1.5-tesla scanner pre-treatment (PT) and at 3, 6, and 9 months of therapy. Routine abdominal MRI sequences and an IVIM-DWI (diffusion-weighted imaging) sequence were obtained. The IVIM-DWI sequence was executed with 16 b-values varying from 0 to 1,400 s/mm2. The mean values of apparent diffusion coefficient (ADC), true diffusion (D), pseudodiffusion (D*), and perfusion fraction (f) of each metastasis were obtained for all b-values, and the time-related changes were recorded to analyze the chronologic responses to antiangiogenic therapy. The RECIST 1.1 criteria were used for the evaluation of treatment response.

RESULTS

The diameters of the metastases diminished significantly at 9 months when compared with PT (p = 0.03). The D (p = 0.10) and ADC (p = 0.21) values of the metastases increased at 9 months of therapy. D* was the highest at 3 months (p =0.24); it decreased at 6 (p =0.97) and 9 months (p =0.87) of therapy. The f value had peaked at 3 months (p =0.51) and started to decrease thereafter. At 6 months, f decreased to the lowest values (p =0.12).

CONCLUSION

IVIM parameters, particularly the perfusion fraction, may quantitatively reflect the response to antiangiogenic treatment. The antiangiogenic response manifests after 3 months of therapy before the RECIST-related response.

Advances in Diffusion and Perfusion MRI for Quantitative Cancer Imaging

PURPOSE OF REVIEW

This article is to review recent technical developments and their clinical applications in cancer imaging quantitative measurement of cellular and vascular properties of the tumors.

RECENT FINDINGS

Rapid development of fast Magnetic Resonance Imaging (MRI) technologies over last decade brought new opportunities in quantitative MRI methods to measure both cellular and vascular properties of tumors simultaneously.

SUMMARY

Diffusion MRI (dMRI) and dynamic contrast enhanced (DCE)-MRI have become widely used to assess the tissue structural and vascular properties, respectively. However, the ultimate potential of these advanced imaging modalities has not been fully exploited. The dependency of dMRI on the diffusion weighting gradient strength and diffusion time can be utilized to measure tumor perfusion, cellular structure, and cellular membrane permeability. Similarly, DCE-MRI can be used to measure vascular and cellular membrane permeability along with cellular compartment volume fractions. To facilitate the understanding of these potentially important methods for quantitative cancer imaging, we discuss the basic concepts and recent developments, as well as future directions for further development.

Repeatability of radiomics and machine learning for DWI: Short-term repeatability study of 112 patients with prostate cancer

PURPOSE

To evaluate repeatability of prostate DWI-derived radiomics and machine learning methods for prostate cancer (PCa) characterization.

METHODS

A total of 112 patients with diagnosed PCa underwent 2 prostate MRI examinations (Scan1 and Scan2) performed on the same day. DWI was performed using 12 b-values (0-2000 s/mm² ), post-processed using kurtosis function, and PCa areas were annotated using whole mount prostatectomy sections. A total of 1694 radiomic features including Sobel, Kirch, Gradient, Zernike Moments, Gabor, Haralick, CoLIAGe, Haar wavelet coefficients, 3D analogue to Laws features, 2D contours, and corner detectors were calculated. Radiomics and 4 feature pruning methods (area under the receiver operator characteristic curve, maximum relevance minimum redundancy, Spearman's ρ, Wilcoxon rank-sum) were evaluated in terms of Scan1-Scan2 repeatability using intraclass correlation coefficient (ICC)(3,1). Classification performance for clinically significant and insignificant PCa with Gleason grade groups 1 versus >1 was evaluated by area under the receiver operator characteristic curve in unseen random 30% data split.

RESULTS

The ICC(3,1) values for conventional radiomics and feature pruning methods were in the range of 0.28-0.90. The machine learning classifications varied between Scan1 and Scan2 with % of same class labels between Scan1 and Scan2 in the range of 61-81%. Surface-to-volume ratio and corner detector-based features were among the most represented features with high repeatability, ICC(3,1) >0.75, consistently high ranking using all 4 feature pruning methods, and classification performance with area under the receiver operator characteristic curve >0.70.

CONCLUSION

Surface-to-volume ratio and corner detectors for prostate DWI led to good classification of unseen data and performed similarly in Scan1 and Scan2 in contrast to multiple conventional radiomic features.

Non-Invasive Monitoring of Stromal Biophysics with Targeted Depletion of Hyaluronan in Pancreatic Ductal Adenocarcinoma

Pancreatic ductal adenocarcinoma (PDA) is characterized by a pronounced fibroinflammatory stromal reaction consisting of inordinate levels of hyaluronan (HA), collagen, immune cells, and activated fibroblasts that work in concert to generate a robust physical barrier to the perfusion and diffusion of small molecule therapeutics. The targeted depletion of hyaluronan with a PEGylated recombinant human hyaluronidase (PEGPH20) lowers interstitial gel–fluid pressures and re-expands collapsed intratumoral vasculature, improving the delivery of concurrently administered agents. Here we report a non-invasive means of assessing biophysical responses to stromal intervention with quantitative multiparametric magnetic resonance imaging (MRI) at 14 Tesla (T). We found that spin-spin relaxation time T2 values and glycosaminoglycan chemical exchange saturation transfer (GagCEST) values decreased at 24 h, reflecting depletion of intratumoral HA content, and that these parameters recovered at 7 days concurrent with replenishment of intratumoral HA. This was also reflected in an increase in low-b apparent diffusion coefficient (ADC) at 24 h, consistent with improved tumor perfusion that again normalized at 7 days after treatment. Phantom imaging suggests that the GagCEST signal is driven by changes in HA versus other glycosaminoglycans. Thus, multiparametric magnetic resonance imaging (MRI) can be used as a non-invasive tool to assess therapeutic responses to targeted stromal therapy in PDA and likely other stroma-rich solid tumors that have high levels of hyaluronan and collagen.

Evaluation of pancreatic tumor development in KPC mice using multi-parametric MRI

Background

Pancreatic ductal adenocarcinoma (PDA) is a fatal disease with very poor prognosis. Development of sensitive and noninvasive methods to monitor tumor progression in PDA is a critical and unmet need. Magnetic resonance imaging (MRI) can noninvasively provide information regarding underlying pathophysiological processes such as necrosis, inflammatory changes and fibrotic tissue deposition.

Methods

A genetically engineered KPC mouse model that recapitulates human PDA was used to characterize disease progression. MR measures of T₁ and T₂ relaxation times, magnetization transfer ratio (MTR), diffusion and chemical exchange saturation transfer were compared in two separate phases i.e. slow and rapid growth phase of tumor. Fibrotic tissue accumulation was assessed histologically using Masson’s trichrome staining. Pearson correlation coefficient (r) was computed to assess the relationship between the fibrotic tissue accumulation and different MR parameters.

Results

There was a negative correlation between amide proton transfer signal intensity and tumor volume (r = − 0.63, p = 0.003) in the slow growth phase of the tumor development. In the terminal stage of rapid growth phase of the tumor development MTR was strongly correlated with tumor volume (r = 0.62, p = 0.008). Finally, MTR was significantly correlated with % fibrosis (r = 0.87; p < 0.01), followed by moderate correlation between tumor volume (r = 0.42); T₁ (r = − 0.61), T₂ (r = − 0.61) and accumulation of fibrotic tissue.

Conclusions

Here we demonstrated, using multi-parametric MRI (mp-MRI), that MRI parameters changed with tumor progression in a mouse model of PDA. Use of mp-MRI may have the potential to monitor the dynamic changes of tumor microenvironment with increase in tumor size in the transgenic KPC mouse model of pancreatic tumor.

Electronic supplementary material

The online version of this article (10.1186/s40644-018-0172-6) contains supplementary material, which is available to authorized users.

Accurate IVIM model-based liver lesion characterisation can be achieved with only three b-value DWI

OBJECTIVE

The objective of this study was to evaluate a simplified intravoxel incoherent motion (IVIM) approach of diffusion-weighted imaging (DWI) with four b-values for liver lesion characterisation at 1.5 T.

METHODS

DWI data from a respiratory-gated MRI sequence with b = 0, 50, 250, 800 s/mm² were retrospectively analysed in 173 lesions and 40 healthy livers. The apparent diffusion coefficient ADC = ADC(0,800) and IVIM-based parameters D₁' = ADC(50,800), D₂' =ADC(250,800), f₁', f₂', D*', ADC_low = ADC(0,50), and ADC_diff=ADC_low-D₂' were calculated voxel-wise without fitting procedures. Differences between lesion groups were investigated.

RESULTS

Focal nodular hyperplasias were best discriminated from all other lesions by f₁' with an area under the curve (AUC) of 0.989. Haemangiomas were best discriminated by D₁' (AUC of 0.994). For discrimination between malignant and benign lesions, ADC(0,800) and D₁' were best suited (AUC of 0.915 and 0.858, respectively). Discriminatory power was further increased by using a combination of D₁' and f₁'.

CONCLUSION

IVIM parameters D and f approximated from three b-values provided more discriminatory power between liver lesions than ADC determined from two b-values. The use of b = 0, 50, 800 s/mm² was superior to that of b = 0, 250, 800 s/mm². The acquisition of four instead of three b-values has no further benefit for lesion characterisation.

KEY POINTS

• Diffusion and perfusion characteristics are assessable with only three b-values. • Association of b = 0, 50, 800 s/mm²is superior to b = 0, 250, 800 s/mm². • A fourth acquired b-value has no benefit for differential diagnosis. • For liver lesion characterisation, simplified IVIM analysis is superior to ADC determination. • Simplified IVIM approach guarantees numerically stable, voxel-wise results and short acquisition times.

Identification of Potential MR-Derived Biomarkers for Tumor Tissue Response to 177Lu-Octreotate Therapy in an Animal Model of Small Intestine Neuroendocrine Tumor

Magnetic resonance (MR) methods enable noninvasive, regional tumor therapy response assessment, but associations between MR parameters, underlying biology, and therapeutic effects must be investigated. The aim of this study was to investigate response assessment efficacy and biological associations of MR parameters in a neuroendocrine tumor (NET) model subjected to radionuclide treatment. Twenty-one mice with NETs received 177Lu-octreotate at day 0. MR experiments (day -1, 1, 3, 8, and 13) included T2-weighted, dynamic contrast-enhanced (DCE) and diffusion-weighted imaging (DWI) and relaxation measurements (T1/T2*). Tumor tissue was analyzed using proteomics. MR-derived parameters were evaluated for each examination day and for different radial distances from the tumor center. Response assessment efficacy and biological associations were evaluated using feature selection and protein expression correlations, respectively. Reduced tumor growth rate or shrinkage was observed until day 8, followed by reestablished growth in most tumors. The most important MR parameter for response prediction was DCE-MRI-derived pretreatment signal enhancement ratio (SER) at 40% to 60% radial distance, where it correlated significantly also with centrally sampled protein CCD89 (association: DNA damage and repair, proliferation, cell cycle arrest). The second most important was changed diffusion (D) between day -1 and day 3, at 60% to 80% radial distance, where it correlated significantly also with peripherally sampled protein CATA (association: oxidative stress, proliferation, cell cycle arrest, apoptotic cell death). Important information regarding tumor biology in response to radionuclide therapy is reflected in several MR parameters, SER and D in particular. The spatial and temporal information provided by MR methods increases the sensitivity for tumor therapy response.

Investigating Low-Velocity Fluid Flow in Tumors with Convection-MRI

Several distinct fluid flow phenomena occur in solid tumors, including intravascular blood flow and interstitial convection. Interstitial fluid pressure is often raised in solid tumors, which can limit drug delivery. To probe low-velocity flow in tumors resulting from raised interstitial fluid pressure, we developed a novel MRI technique named convection-MRI, which uses a phase-contrast acquisition with a dual-inversion vascular nulling preparation to separate intra- and extravascular flow. Here, we report the results of experiments in flow phantoms, numerical simulations, and tumor xenograft models to investigate the technical feasibility of convection-MRI. We observed a significant correlation between estimates of effective fluid pressure from convection-MRI with gold-standard, invasive measurements of interstitial fluid pressure in mouse models of human colorectal carcinoma. Our results show how convection-MRI can provide insights into the growth and responsiveness to vascular-targeting therapy in colorectal cancers.Significance: A noninvasive method for measuring low-velocity fluid flow caused by raised fluid pressure can be used to assess changes caused by therapy. Cancer Res; 78(7); 1859-72. ©2018 AACR.

Optimization of prostate MRI acquisition and post-processing protocol: a pictorial review with access to acquisition protocols

The aim of this review article is to provide insight into the optimization of 1.5-Testla (T) and 3-T prostate magnetic resonance imaging (MRI). An approach for optimization of data quantification, especially diffusion-weighted imaging (DWI), is provided. Benefits and limitations of various pulse sequences are discussed. Importable MRI protocols and access to imaging datasets is provided. Careful optimization of prostate MR acquisition protocol allows the acquisition of high-quality prostate MRI using clinical 1.5-T/3-T MR scanners with an overall acquisition time < 15 min.

Intravoxel incoherent motion (IVIM) histogram biomarkers for prediction of neoadjuvant treatment response in breast cancer patients

OBJECTIVE

To examine the prognostic capabilities of intravoxel incoherent motion (IVIM) metrics and their ability to predict response to neoadjuvant treatment (NAT). Additionally, to observe changes in IVIM metrics between pre- and post-treatment MRI.

METHODS

This IRB-approved, HIPAA-compliant retrospective study observed 31 breast cancer patients (32 lesions). Patients underwent standard bilateral breast MRI along with diffusion-weighted imaging before and after NAT. Six patients underwent an additional IVIM-MRI scan 12-14 weeks after initial scan and 2 cycles of treatment. In addition to apparent diffusion coefficients (ADC) from monoexponential decay, IVIM mean values (tissue diffusivity D_t, perfusion fraction f_p, and pseudodiffusivity D_p) and histogram metrics were derived using a biexponential model. An additional filter identified voxels of highly vascular tumor tissue (VTT), excluding necrotic or normal tissue. Clinical data include histology of biopsy and clinical response to treatment through RECIST assessment. Comparisons of treatment response were made using Wilcoxon rank-sum tests.

RESULTS

Average, kurtosis, and skewness of pseudodiffusion D_p significantly differentiated RECIST responders from nonresponders. ADC and D_t values generally increased (∼70%) and VTT% values generally decreased (∼20%) post-treatment.

CONCLUSION

D_p metrics showed prognostic capabilities; slow and heterogeneous pseudodiffusion offer poor prognosis. Baseline ADC/D_t parameters were not significant predictors of response. This work suggests that IVIM mean values and heterogeneity metrics may have prognostic value in the setting of breast cancer NAT.

Quantitative assessment of hyperacute cerebral infarction with intravoxel incoherent motion MR imaging: Initial experience in a canine stroke model

PURPOSE

To evaluate the feasibility of intravoxel incoherent motion (IVIM) for the measurement of diffusion and perfusion parameters in hyperacute strokes.

MATERIALS AND METHODS

An embolic ischemic model was established with an autologous thrombus in 20 beagles. IVIM imaging was performed on a 3.0 Tesla platform at 4.5 h and 6 h after embolization. Ten b values from 0 to 900 s/mm² were fitted with a bi-exponential model to extract perfusion fraction f, diffusion coefficient D, and pseudo-diffusion coefficient D*. Additionally, the apparent diffusion coefficient (ADC) was calculated using the mono-exponential model with all the b values. Statistical analysis was performed using the pairwise Student's t test and Pearson's correlation test.

RESULTS

A significant decrease in f and D was observed in the ischemic area when compared with those in the contralateral side at 4.5 h and 6 h after embolization (P < 0.01 for all). No significant difference was observed in D* between the two sides at either time point (P = 0.086 and 0.336, respectively). In the stroke area, f at 6 h was significantly lower than that at 4.5 h (P = 0.016). A significantly positive correlation was detected between ADC and D in both stroke and contralateral sides at 4.5 h and 6 h (P < 0.001 for both). Significant correlation between ADC and f was only observed in the contralateral side at 4.5 h and 6 h (P = 0.019 and 0.021, respectively).

CONCLUSION

IVIM imaging could simultaneously evaluate the diffusion and microvascular perfusion characteristics in hyperacute strokes.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. MAGN. RESON. IMAGING 2017;46:550-556.

Intravoxel incoherent motion diffusion-weighted MRI of the abdomen: The effect of fitting algorithms on the accuracy and reliability of the parameters

PURPOSE

To evaluate the influence of fitting methods on the accuracy and reliability of intravoxel incoherent motion (IVIM) parameters, with a particular emphasis on the constraint function.

MATERIALS AND METHODS

Diffusion-weighted (DW) imaging data were analyzed using IVIM-based full-fitting (simultaneous fit of all parameters) and segmented-fitting (step-by-step fit of each parameter), each with and without the constraint function, to estimate the molecular diffusion coefficient (D_slow ), perfusion fraction (f), and flow-related diffusion coefficient (D_fast ). Computational simulations were performed at variable signal-to-noise ratios to evaluate the relative error (RE) and coefficient of variation (CV) of the estimated IVIM parameters. DW imaging of the abdomen was performed twice at 1.5 Tesla using nine b-values (0-900 s/mm² ) in 12 health volunteers (6 men and 6 women; mean age: 30 years). The measurement repeatability of IVIM parameters in the liver and the pancreas was evaluated using the within-subject coefficient of variation (_w CV).

RESULTS

In simulations, full-fitting without the constraint function yielded the largest RE (P < 0.001 for D_slow and f; P ≤ 0.044 for D_fast ) and CV (P ≤ 0.033 for D_slow and f; P ≤ 0.473 for D_fast ) for IVIM parameters among all four algorithms. In volunteer imaging, full-fitting without the constraint function also resulted in the poorest repeatability for D_slow (_w CV, 17.12%-65.45%) and f (_w CV, 19.35%-42.84%) in the liver and pancreas, while the other algorithms had similar repeatability values (_w CV, 4.05%-11.99% for D_slow and 9.65%-18.66% for f). Measurement repeatability of D_fast (_w CV, 29.52%-85.01%) was the poorest among the IVIM parameters.

CONCLUSION

For accurate and reliable measurement of IVIM parameters, segmented fitting or full-fitting with the constraint function should be used for IVIM-based analysis of DW imaging.

LEVEL OF EVIDENCE

3 Technical Efficacy: Stage 2 J. MAGN. RESON. IMAGING 2017;45:1637-1647.

Intravoxel incoherent motion model-based analysis of diffusion-weighted magnetic resonance imaging with 3 b-values for response assessment in locoregional therapy of hepatocellular carcinoma

Purpose

The aim of this study was to evaluate an intravoxel incoherent motion (IVIM) model–based analysis of diffusion-weighted imaging (DWI) for assessing the response of hepatocellular carcinoma (HCC) to locoregional therapy.

Patients and methods

Respiratory-gated DWI (b=0, 50, and 800 s/mm²) was retrospectively analyzed in 25 patients who underwent magnetic resonance imaging at 1.5 T before and 6 weeks following the first cycle of transarterial chemoembolization therapy, transarterial ethanol-lipiodol embolization therapy, and transarterial radioembolization therapy. In addition to the determination of apparent diffusion coefficient, ADC(0,800), an estimation of the diffusion coefficient, D′, and the perfusion fraction, f′, was performed by using a simplified IVIM approach. Parameters were analyzed voxel-wise. Tumor response was assessed in a central slice by using a region of interest (ROI) covering the whole tumor. HCCs were categorized into two groups, responders and nonresponders, according to tumor size changes on first and second follow ups (if available) and changes of contrast-enhanced region on the first follow up.

Results

In total, 31 HCCs were analyzed: 17 lesions were assigned to responders and 14 were to nonresponders. In responders, ADC(0,800) and D′ were increased after therapy by ~30% (P=0.00004) and ~42% (P=0.00001), respectively, whereas f′ was decreased by ~37% (P=0.00094). No significant changes were found in nonresponders. Responders and nonresponders were better differentiated by changes in D′ than by changes in ADC(0,800) (area under the curve =0.878 vs 0.819 or 0.714, respectively).

Conclusion

In patients with HCCs undergoing embolization therapy, diffusion changes were better reflected by D′ than by conventional ADC(0,800), which is influenced by counteracting perfusion changes as assessed by f′.

Intravoxel Incoherent Motion Diffusion-Weighted MR Imaging for Prediction of Early Arterial Blood Flow Stasis in Radioembolization of Breast Cancer Liver Metastases

PURPOSE

To retrospectively evaluate predictive value of intravoxel incoherent motion (IVIM) diffusion-weighted imaging (DWI) for early arterial blood flow stasis during transarterial radioembolization (TARE) of liver dominant breast metastases (LdBM).

MATERIALS AND METHODS

Preinterventional 1.5T DWI (b0, b1, b2 = 0, 50, 800 s/mm(2)) data for 28 liver lobes of 18 female patients treated by resin-based radioembolization (10 bilobar and 8 unilobar treatments) were analyzed. Apparent diffusion coefficient (ADC) (0, 800) and an estimation of the true diffusion coefficient D' and of the perfusion fraction f' were calculated for the 2 largest metastases. Response rate at 3 months and survival were analyzed. Procedures without full dose application because of early stasis were assigned to group A (n = 15), and procedures with full dose application were assigned to group B (n = 13).

RESULTS

Metastases in group A showed significantly lower f' (0.035 ± 0.018 vs 0.076 ± 0.015, P < .0001) and a trend toward lower ADC(0, 800) with values given in 10(-6) mm(2)/s (1,066 ± 141 vs 1,189 ± 176, P = .051); no group difference was shown for D'. Groups were best discriminated by weighted mean f' values of the 2 largest metastases with accuracy of 100%. Mean tumor diameter before and after TARE was 51 mm ± 18 and 50 mm ± 24 in group A and 47 mm ± 27 and 48 mm ± 32 for group B. Imaging response did not differ between groups (P = .545). Overall survival did not differ significantly between group A (230 d) and B (155 d) (P = .124).

CONCLUSIONS

Perfusion-sensitive IVIM parameter f' may predict early blood flow stasis in patients undergoing TARE for LdBM. Determination of this parameter before intervention may increase awareness of the interventionalist and increase safety of microsphere administration.

The value of intravoxel incoherent motion model-based diffusion-weighted imaging for outcome prediction in resin-based radioembolization of breast cancer liver metastases

Purpose

To evaluate prognostic values of clinical and diffusion-weighted magnetic resonance imaging-derived intravoxel incoherent motion (IVIM) parameters in patients undergoing primary radioembolization for metastatic breast cancer liver metastases.

Subjects and methods

A total of 21 females (mean age 54 years, range 43–72 years) with liver-dominant metastatic breast cancer underwent standard liver magnetic resonance imaging (1.5 T, diffusion-weighted imaging with b-values of 0, 50, and 800 s/mm²) before and 4–6 weeks after radioembolization. The IVIM model-derived estimated diffusion coefficient D’ and the perfusion fraction f’ were evaluated by averaging the values of the two largest treated metastases in each patient. Kaplan–Meier and Cox regression analyses for overall survival (OS) were performed. Investigated parameters were changes in f’- and D’-values after therapy, age, sex, Eastern Cooperative Oncology Group (ECOG) status, grading of primary tumor, hepatic tumor burden, presence of extrahepatic disease, baseline bilirubin, previous bevacizumab therapy, early stasis during radioembolization, chemotherapy after radioembolization, repeated radioembolization and Response Evaluation Criteria in Solid Tumors (RECIST) response at 6-week follow-up.

Results

Median OS after radioembolization was 6 (range 1.5–54.9) months. In patients with therapy-induced decreasing or stable f’-values, median OS was significantly longer than in those with increased f’-values (7.6 [range 2.6–54.9] vs 2.6 [range 1.5–17.4] months, P<0.0001). Longer median OS was also seen in patients with increased D’-values (6 [range 1.6–54.9] vs 2.8 [range 1.5–17.4] months, P=0.008). Patients with remission or stable disease (responders) according to RECIST survived longer than nonresponders (7.2 [range 2.6–54.9] vs 2.6 [range 1.5–17.4] months, P<0.0001). An ECOG status ≤1 resulted in longer median OS than >1 (7.6 [range 2.6–54.9] vs 1.7 [range 1.5–4.5] months, P<0.0001). Pretreatment IVIM parameters and the other clinical characteristics were not associated with OS. Classification by f’-value changes and ECOG status remained as independent predictors of OS on multivariate analysis, while RECIST response and D’-value changes did not predict survival.

Conclusion

Following radioembolization of breast cancer liver metastases, early changes in the IVIM model-derived perfusion fraction f’ and baseline ECOG score were predictive of patient outcome, and may thus help to guide treatment strategy.

Evaluation of a Simplified Intravoxel Incoherent Motion (IVIM) Analysis of Diffusion-Weighted Imaging for Prediction of Tumor Size Changes and Imaging Response in Breast Cancer Liver Metastases Undergoing Radioembolization: A Retrospective Single Center Analysis

To investigate the value of a simplified intravoxel incoherent motion (IVIM) analysis for evaluation of therapy-induced tumor changes and response of breast cancer liver metastases (mBRC) undergoing radioembolization.In 21 females (mean age 54 years, range 43-72) with mBRC tumor size changes and response evaluation criteria in solid tumors (RECIST) response to 26 primary radioembolization procedures were analyzed. Standard 1.5-T liver magnetic resonance imaging including respiratory-gated diffusion-weighted imaging (DWI) with b0 = 0 s/mm, b1 = 50 s/mm, b2 = 800 s/mm before and 6 weeks after each treatment was performed. In addition to the apparent diffusion coefficient (ADC)(0,800), the estimated diffusion coefficient D' and the perfusion fraction f' were determined using a simplified IVIM approach. For each radioembolization, the 2 largest treated metastases (if available) were analyzed. Lesions were categorized according to size changes into group A (reduction of longest diameter [LD]) and group B (LD increase) after 3 months. Radioembolization procedures were further categorized into "response" (partial response and stable disease) and "nonresponse" (progressive disease) according to RECIST after 3 months. ADC and D' are given in 10 mm/s.Forty-five metastases were analyzed. Thirty-two lesions were categorized as A; 13 as B. Before therapy, group A lesions showed significantly larger f'-values than B (P = 0.001), but ADC(0,800) and D' did not differ. After therapy, in group A lesions the ADC(0,800)- and D'-values increased and f' decreased (P < 0.0001); in contrast in group B lesions f' increased (P = 0.001). Groups could be differentiated by preinterventional f' and by changes of D' and f' between pre and postinterventional imaging (area under the curve [AUC] of 0.903, 0.747 and 1.0, respectively).Preinterventional parameters did not differ between responders and nonresponders according to RECIST. ADC(0,800)- and D'-values showed a larger increase in responders compared with nonresponders (P = 0.013 and P = 0.001, respectively). After therapy f'-values decreased significantly in responders (P = 0.001). Good to excellent prediction of long-term RECIST response was possible by therapy-induced changes in LD, D', and f' (AUC 0.903, 0.879, and 0.867, respectively).A simplified IVIM model-based analysis of early post-treatment DWI can deliver additional information on tumor size changes and long-term RECIST response after radioembolization of mBRC. The estimated perfusion fraction f' is better suited for response assessment than the conventional ADC(0,800) or D'. This can be useful to guide further treatment strategy.

Stroke assessment with intravoxel incoherent motion diffusion-weighted MRI

Intravoxel incoherent motion (IVIM) diffusion-weighted MRI can simultaneously measure diffusion and perfusion characteristics in a non-invasive way. This study aimed to determine the potential utility of IVIM in characterizing brain diffusion and perfusion properties for clinical stroke. The multi-b-value diffusion-weighted images of 101 patients diagnosed with acute/subacute ischemic stroke were retrospectively evaluated. The diffusion coefficient D, representing the water apparent diffusivity, was obtained by fitting the diffusion data with increasing high b-values to a simple mono-exponential model. The IVIM-derived perfusion parameters, pseudodiffusion coefficient D*, vascular volume fraction f and blood flow-related parameter fD*, were calculated with the bi-exponential model. Additionally, the apparent diffusion coefficient (ADC) was fitted according to the mono-exponential model using all b-values. The diffusion parameters for the ischemic lesion and normal contralateral region were measured in each patient. Statistical analysis was performed using the paired Student t-test and Pearson correlation test. Diffusion data in both the ischemic lesion and normal contralateral region followed the IVIM bi-exponential behavior, and the IVIM model showed better goodness of fit than the mono-exponential model with lower Akaike information criterion values. The paired Student t-test revealed significant differences for all diffusion parameters (all P < 0.001) except D* (P = 0.218) between ischemic and normal areas. For all patients in both ischemic and normal regions, ADC was significantly positively correlated with D (both r = 1, both P < 0.001) and f (r = 0.541, P < 0.001; r = 0.262, P = 0.008); significant correlation was also found between ADC and fD* in the ischemic region (r = 0.254, P = 0.010). For all pixels within the region of interest from a representative subject in both ischemic and normal regions, ADC was significantly positively correlated with D (both r = 1, both P < 0.001), f (r = 0.823, P < 0.001; r = 0.652, P < 0.001) and fD* (r = 0.294, P < 0.001; r = 0.340, P < 0.001). These findings may have clinical implications for the use of IVIM imaging in the assessment and management of acute/subacute stroke patients. Copyright © 2016 John Wiley & Sons, Ltd.

Fitting methods for intravoxel incoherent motion imaging of prostate cancer on region of interest level: Repeatability and gleason score prediction

PURPOSE

To evaluate different fitting methods for intravoxel incoherent motion (IVIM) imaging of prostate cancer in the terms of repeatability and Gleason score prediction.

METHODS

Eighty-one patients with histologically confirmed prostate cancer underwent two repeated 3 Tesla diffusion-weighted imaging (DWI) examinations performed using 14 b-values in the range of 0-500 s/mm² and diffusion time of 19.004 ms. Mean signal intensities of regions-of-interest were fitted using five different fitting methods for IVIM as well as monoexponential, kurtosis, and stretched exponential models. The fitting methods and models were evaluated in the terms of fitting quality [Akaike information criteria (AIC)], repeatability, and Gleason score prediction. Tumors were classified into three groups (3 + 3, 3 + 4, > 3 + 4). Machine learning algorithms were used to evaluate the performance of the combined use of the parameters. Simulation studies were performed to evaluate robustness of the fitting methods against noise.

RESULTS

Monoexponential model was preferred over IVIM based on AIC. The "pseudodiffusion" parameters demonstrated low repeatability and clinical value. Median "pseudodiffusion" fraction values were below 8.00%. Combined use of the parameters did not outperform the monoexponential model.

CONCLUSION

Monoexponential model demonstrated the highest repeatability and clinical values in the regions-of-interest based analysis of prostate cancer DWI, b-values in the range of 0-500 s/mm² . Magn Reson Med 77:1249-1264, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Evaluation of breast cancer using intravoxel incoherent motion (IVIM) histogram analysis: comparison with malignant status, histological subtype, and molecular prognostic factors

PURPOSE

To examine heterogeneous breast cancer through intravoxel incoherent motion (IVIM) histogram analysis.

MATERIALS AND METHODS

This HIPAA-compliant, IRB-approved retrospective study included 62 patients (age 48.44 ± 11.14 years, 50 malignant lesions and 12 benign) who underwent contrast-enhanced 3 T breast MRI and diffusion-weighted imaging. Apparent diffusion coefficient (ADC) and IVIM biomarkers of tissue diffusivity (Dt), perfusion fraction (fp), and pseudo-diffusivity (Dp) were calculated using voxel-based analysis for the whole lesion volume. Histogram analysis was performed to quantify tumour heterogeneity. Comparisons were made using Mann-Whitney tests between benign/malignant status, histological subtype, and molecular prognostic factor status while Spearman's rank correlation was used to characterize the association between imaging biomarkers and prognostic factor expression.

RESULTS

The average values of the ADC and IVIM biomarkers, Dt and fp, showed significant differences between benign and malignant lesions. Additional significant differences were found in the histogram parameters among tumour subtypes and molecular prognostic factor status. IVIM histogram metrics, particularly fp and Dp, showed significant correlation with hormonal factor expression.

CONCLUSION

Advanced diffusion imaging biomarkers show relationships with molecular prognostic factors and breast cancer malignancy. This analysis reveals novel diagnostic metrics that may explain some of the observed variability in treatment response among breast cancer patients.

KEY POINTS

• Novel IVIM biomarkers characterize heterogeneous breast cancer. • Histogram analysis enables quantification of tumour heterogeneity. • IVIM biomarkers show relationships with breast cancer malignancy and molecular prognostic factors.

Comparison of fitting methods and b-value sampling strategies for intravoxel incoherent motion in breast cancer

PURPOSE

To compare fitting methods and sampling strategies, including the implementation of an optimized b-value selection for improved estimation of intravoxel incoherent motion (IVIM) parameters in breast cancer.

METHODS

Fourteen patients (age, 48.4 ± 14.27 years) with cancerous lesions underwent 3 Tesla breast MRI examination for a HIPAA-compliant, institutional review board approved diffusion MR study. IVIM biomarkers were calculated using "free" versus "segmented" fitting for conventional or optimized (repetitions of key b-values) b-value selection. Monte Carlo simulations were performed over a range of IVIM parameters to evaluate methods of analysis. Relative bias values, relative error, and coefficients of variation (CV) were obtained for assessment of methods. Statistical paired t-tests were used for comparison of experimental mean values and errors from each fitting and sampling method.

RESULTS

Comparison of the different analysis/sampling methods in simulations and experiments showed that the "segmented" analysis and the optimized method have higher precision and accuracy, in general, compared with "free" fitting of conventional sampling when considering all parameters. Regarding relative bias, IVIM parameters fp and Dt differed significantly between "segmented" and "free" fitting methods.

CONCLUSION

IVIM analysis may improve using optimized selection and "segmented" analysis, potentially enabling better differentiation of breast cancer subtypes and monitoring of treatment.

Ductal carcinoma in situ of the breasts: review of MR imaging features

The incidence of ductal carcinoma in situ (DCIS) has increased over the past few decades and now accounts for over 20% of newly diagnosed cases of breast cancer. Although the detection of DCIS has increased with the advent of widespread mammography screening, it is essential to have a more accurate assessment of the extent of DCIS for successful breast conservation therapy. Recent studies evaluating the detection of DCIS with magnetic resonance (MR) imaging have used high spatial resolution techniques and have increasingly been performed to screen a high-risk population as well as to evaluate the extent of disease. This work has shown that MR imaging is the most sensitive modality currently available for identifying DCIS and is more accurate than mammography in evaluating the extent of DCIS. MR imaging is particularly sensitive for identifying high-grade and intermediate-grade DCIS. DCIS may have variable morphologic features on MR images, with non-mass enhancement morphology being the most common manifestation. Less commonly, DCIS may also manifest as a mass on MR images, in which case it is most likely to be irregular. The kinetics of DCIS are also variable, with fast uptake and a plateau curve reported as the most common kinetic pattern. Additional MR imaging tools such as diffusion-weighted imaging and quantitative kinetic analysis combined with the benefit of high field strength, such as 3 T, may increase the sensitivity and specificity of breast MR imaging in the detection of DCIS.

Early radiation-induced changes evaluated by intravoxel incoherent motion in the major salivary glands

PURPOSE

To investigate the potential of intravoxel incoherent motion (IVIM) MRI for early evaluation of irradiated major salivary glands.

MATERIALS AND METHODS

Thirty-four patients with head-neck cancer were included in a prospective study. All patients underwent three serial IVIM-MRI: before, half-way through, and at the end of radiotherapy (RT). Apparent diffusion coefficient (ADC), ADClow derived in the low b-value range, perfusion fraction f, and pure diffusion coefficient D were estimated. Pretreatment values and early changes of diffusion parameters were correlated with parotid mean dose (Dmean ) and volume reduction after RT.

RESULTS

Changes in diffusion parameters over time were all significant (P < 0.001 for ADC, ADClow , and D, P = 0.003 for f). Variations of ADC, ADClow , and f were not correlated with Dmean (P = 0.089, P = 0.252 and P = 0.884, respectively), whereas a significant relationship was found between changes in D and Dmean (r = 0.197 with CI95% = 0.004-0.375, P = 0.046). Pretreatment f and Dmean were the best independent predictors for the percentage shrinkage (P = 0.0003 and 0.0597 respectively; R(2) = 0.391).

CONCLUSION

Early changes of irradiated major salivary glands can be noninvasively evaluated by IVIM-MRI. Perfusion-related coefficients in conjunction with dosimetric information increase our capability to predict the change in parotid volume and hence, if further validated, guide treatment strategy in RT.

High-field small animal magnetic resonance oncology studies

This review focuses on the applications of high magnetic field magnetic resonance imaging (MRI) and spectroscopy (MRS) to cancer studies in small animals. High-field MRI can provide information about tumor physiology, the microenvironment, metabolism, vascularity and cellularity. Such studies are invaluable for understanding tumor growth and proliferation, response to treatment and drug development. The MR techniques reviewed here include (1)H, (31)P, chemical exchange saturation transfer imaging and hyperpolarized (13)C MRS as well as diffusion-weighted, blood oxygen level dependent contrast imaging and dynamic contrast-enhanced MRI. These methods have been proven effective in animal studies and are highly relevant to human clinical studies.

Assessment of diffusion parameters by intravoxel incoherent motion MRI in head and neck squamous cell carcinoma

The objectives of this study were to assess the diffusion parameters derived from intravoxel incoherent motion (IVIM) MRI in head and neck squamous cell carcinoma (HNSCC) and to investigate the agreement between different methods of tumor delineation and two numerical methods to extract the perfusion fraction f. Thirty-seven untreated patients with histopathologically confirmed primary HNSCC were included retrospectively in the study. The entire volume of the primary tumor was outlined on diffusion-weighted images using co-registered morphological images as a guide to the tumor location. Apparent diffusion coefficient (ADC) and IVIM diffusion parameters were estimated considering the largest tumor section as well as the entire tumor volume. A bi-exponential fit was implemented to extract f, D (pure diffusion coefficient) and D* (pseudo-diffusion coefficient). A second simplified method, based on an asymptotic extrapolation, was used to determine f. The agreement between ADC and IVIM diffusion parameters derived from the delineation of single and multiple slices, and between the two f estimations, was assessed by Bland-Altman plots. The inter-slice variability of ADC and IVIM diffusion parameters was evaluated. The Kruskal-Wallis test was used to investigate whether the tumor location had a statistically significant influence on the values of the parameters. Comparing the tumor delineation methods, a better accordance was found for ADC and D, with a mean percentage difference of less than 2%. Larger discrepancies were found for f and D*, with mean differences of 4.5% and 5.5%, respectively. When comparing the two f estimation methods, small mean differences were found (<3.5%), suggesting that the two methods may be considered as equivalent for the assessment of f in our patient population. The observed ADC and IVIM diffusion parameters were dependent on the anatomic site of the lesion, carcinoma of the nasopharynx showing more homogeneous and dissimilar estimations than other HNSCCs.

Intravoxel incoherent motion diffusion-weighted MRI at 3.0 T differentiates malignant breast lesions from benign lesions and breast parenchyma

PURPOSE

To study the differentiation of malignant breast lesions from benign lesions and fibroglandular tissue (FGT) using apparent diffusion coefficient (ADC) and intravoxel incoherent motion (IVIM) parameters.

MATERIALS AND METHODS

This retrospective study included 26 malignant and 14 benign breast lesions in 35 patients who underwent diffusion-weighted MRI at 3.0T and nine b-values (0-1000 s/mm(2) ). ADC and IVIM parameters (perfusion fraction fp , pseudodiffusion coefficient Dp , and true diffusion coefficient Dd ) were determined in lesions and FGT. For comparison, IVIM was also measured in 16 high-risk normal patients. A predictive model was constructed using linear discriminant analysis. Lesion discrimination based on ADC and IVIM parameters was assessed using receiver operating characteristic (ROC) and area under the ROC curve (AUC).

RESULTS

In FGT of normal subjects, fp was 1.1 ± 1.1%. In malignant lesions, fp (6.4 ± 3.1%) was significantly higher than in benign lesions (3.1 ± 3.3%, P = 0.0025) or FGT (1.5 ± 1.2%, P < 0.001), and Dd ((1.29 ± 0.28) × 10(-3) mm(2) /s) was lower than in benign lesions ((1.56 ± 0.28) × 10(-3) mm(2) /s, P = 0.011) or FGT ((1.86 ± 0.34) × 10(-3) mm(2) /s, P < 0.001). A combination of Dd and fp provided higher AUC for discrimination between malignant and benign lesions (0.84) or FGT (0.97) than ADC (0.72 and 0.86, respectively).

CONCLUSION

The IVIM parameters provide accurate identification of malignant lesions.

MR diffusion tensor imaging of normal kidneys

PURPOSE

To assess the feasibility of diffusion tensor imaging (DTI) of normal kidneys and the influence of hydration state.

MATERIALS AND METHODS

Ten healthy volunteers underwent renal DTI after fasting for 12 hours and 4 hours, without fasting, and following water diuresis. Medullary and cortical apparent diffusion coefficient (ADC) and fractional anisotropy (FA) values were measured and compared in the four different states of hydration. DTI was performed with a 3T magnetic resonance imaging (MRI) system using fat-saturated single-shot spin-echo echo planar imaging sequence.

RESULTS

ADC of normal cortex (2.387 ± 0.081 × 10(-3) mm(2) /s) was significantly higher (t = 20.126, P = 0) than that of medulla (1.990 ± 0.063 × 10(-3) mm(2) /s). The FA value of normal cortex (0.282 ± 0.017) was significantly lower (t = -42.713, P = 0) than that of medulla (0.447 ± 0.022). The ADC and FA values of the left renal cortex (2.404 ± 0.082 × 10(-3) mm(2) /s, 0.282 ± 0.017) and medulla (2.002 ± 0.081 × 10(-3) mm(2) /s, 0.452 ± 0.024) were not significantly different (P > 0.05) from those of right renal cortex (2.369 ± 0.080 × 10(-3) mm(2) /s, 0.283 ± 0.018) and medulla (1.978 ± 0.039 × 10(-3) mm(2) /s, 0.443 ± 0.019). Values for ADC (×10(-3) mm(2) /s) and FA in the 12-hour fasting, 4-hour fasting, nonfasting, and water diuresis states were 2.372 ± 0.095 and 0.278 ± 0.018, 2.387 ± 0.081 and 0.282 ± 0.017, 2.416 ± 0.051 and 0.279 ± 0.023, 2.421 ± 0.068, and 0.270 ± 0.021, respectively, in cortex, 1.972 ± 0.084 and 0.438 ± 0.014, 1.990 ± 0.063 and 0.447 ± 0.022, 2.021 ± 0.081 and 0.450 ± 0.031, 2.016 ± 0.076 and 0.449 ± 0.028, respectively, in medulla. The ADC and FA values in different hydration states were not significantly different (P > 0.05).

CONCLUSION

DTI of normal kidneys is feasible with reproducible ADC and FA values independent of hydration states.

Diffusion-weighted imaging with acquisition of three b-values for response evaluation of neuroendocrine liver metastases undergoing selective internal radiotherapy

OBJECTIVES

To evaluate diffusion-weighted MRI with acquisition of three b-values and calculation of fractioned ADCs for response evaluation of neuroendocrine liver metastases undergoing selective internal radiotherapy (SIRT).

METHODS

Ten consecutive patients with neuroendocrine liver metastases underwent MRI before and following SIRT. Diffusion-weighted imaging included acquisition of the b-values 0, 50 and 800 s/mm(2) and calculation of ADC(50,800), ADC(0,50) and ADC(0,800) maps. According to therapy response, lesions were categorised into group A [≥20% reduction of the longest diameter (LD) in comparison to baseline MRI] and group B (<20% reduction of the LD).

RESULTS

Twelve out of 31 metastases were categorised as group A and 19 out of 31 metastases were categorised as group B. Pretherapeutic values of ADC(0,800) and ADC(50,800) did not differ significantly between the two groups; however, ADC(0,50) was 32% lower in group A (P = 0.049). ADC(0,800) and ADC(50,800) increased significantly after therapy in both groups, however, group differences were not statistically significant. Conversely, the increase in ADC(0,50) was about a factor of 7 larger in group A than in group B (P = 0.023).

CONCLUSIONS

Our study showed that the ADC(0,50) is a promising biomarker for response assessment of neuroendocrine liver metastases following SIRT.

KEY POINTS

• Diffusion-weighted MRI offers new information about neuroendocrine hepatic metastases. • Evaluation of perfusion and diffusion components requires fractioned apparent diffusion coefficients (ADCs). • Perfusion effects represented by ADC (0.50) can be observed in neuroendocrine metastases. • Pretherapeutic ADC (0.50) was significantly lower in metastases with a response ≥20%. • Such biomarkers may help evaluate liver metastases in patients undergoing therapy.