Characterization of the diffusion signal of breast tissues using multi-exponential models.

PURPOSE: Restriction spectrum imaging (RSI) decomposes the diffusion-weighted MRI signal into separate components of known apparent diffusion coefficients (ADCs). The number of diffusion components and optimal ADCs for RSI are organ-specific and determined empirically. The purpose of this work was to determine the RSI model for breast tissues. METHODS: The diffusion-weighted MRI signal was described using a linear combination of multiple exponential components. A set of ADC values was estimated to fit voxels in cancer and control ROIs. Later, the signal contributions of each diffusion component were estimated using these fixed ADC values. Relative-fitting residuals and Bayesian information criterion were assessed. Contrast-to-noise ratio between cancer and fibroglandular tissue in RSI-derived signal contribution maps was compared to DCE imaging. RESULTS: A total of 74 women with breast cancer were scanned at 3.0 Tesla MRI. The fitting residuals of conventional ADC and Bayesian information criterion suggest that a 3-component model improves the characterization of the diffusion signal over a biexponential model. Estimated ADCs of triexponential model were D1,3 = 0, D2,3 = 1.5 × 10-3 , and D3,3 = 10.8 × 10-3 mm2 /s. The RSI-derived signal contributions of the slower diffusion components were larger in tumors than in fibroglandular tissues. Further, the contrast-to-noise and specificity at 80% sensitivity of DCE and a subset of RSI-derived maps were equivalent. CONCLUSION: Breast diffusion-weighted MRI signal was best described using a triexponential model. Tumor conspicuity in breast RSI model is comparable to that of DCE without the use of exogenous contrast. These data may be used as differential features between healthy and malignant breast tissues.

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/mm2 . A phase-reversed scan (b = 0 s/mm2 ) 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.

Modeling the diffusion-weighted imaging signal for breast lesions in the b = 200 to 3000 s/mm2 range: quality of fit and classification accuracy for different representations.

PURPOSE: To evaluate different non-Gaussian representations for the diffusion-weighted imaging (DWI) signal in the b-value range 200 to 3000 s/mm2 in benign and malignant breast lesions. METHODS: Forty-three patients diagnosed with benign (n = 18) or malignant (n = 25) tumors of the breast underwent DWI (b-values 200, 600, 1200, 1800, 2400, and 3000 s/mm2 ). Six different representations were fit to the average signal from regions of interest (ROIs) at different b-value ranges. Quality of fit was assessed by the corrected Akaike information criterion (AICc), and the Friedman test was used for assessing representation ranks. The area under the curve (AUC) of receiver operating characteristic curves were used to evaluate the power of derived parameters to differentiate between malignant and benign lesions. The lesion ROI was divided in central and peripheral parts to assess potential effect of heterogeneity. Sensitivity to noise-floor correction was also evaluated. RESULTS: The Padé exponent was ranked as the best based on AICc, whereas 3 models (kurtosis, fractional, and biexponential) achieved the highest AUC = 0.99 for lesion differentiation. The monoexponential model at bmax = 600 s/mm2 already provides AUC = 0.96, with considerably shorter acquisition time and simpler analysis. Significant differences between central and peripheral parts of lesions were found in malignant lesions. The mono- and biexponential models were most stable against varying degrees of noise-floor correction. CONCLUSION: Non-Gaussian representations are required for fitting of the DWI curve at high b-values in breast lesions. However, the added clinical value from the high b-value data for differentiation of benign and malignant lesions is not clear.

Stromal Collagen Content in Breast Tumors Correlates With In Vivo Diffusion-Weighted Imaging: A Comparison of Multi b-Value DWI With Histologic Specimen From Benign and Malignant Breast Lesions.

BACKGROUND: Increased deposition and reorientation of stromal collagen fibers are associated with breast cancer progression and invasiveness. Diffusion-weighted imaging (DWI) may be sensitive to the collagen fiber organization in the stroma and could provide important biomarkers for breast cancer characterization. PURPOSE: To understand how collagen fibers influence water diffusion in vivo and evaluate the relationship between collagen content and the apparent diffusion coefficient (ADC) and the signal fractions of the biexponential model using a high b-value scheme. STUDY TYPE: Prospective. SUBJECTS/SPECIMENS: Forty-five patients with benign (n = 8), malignant (n = 36), and ductal carcinoma in situ (n = 1) breast tumors. Lesions and normal fibroglandular tissue (n = 9) were analyzed using sections of formalin-fixed, paraffin-embedded tissue stained with hematoxylin, erythrosine, and saffron. FIELD STRENGTH/SEQUENCE: MRI (3T) protocols: Protocol I: Twice-refocused spin-echo echo-planar imaging with: echo time (TE) 85 msec; repetition time (TR) 9300/11600 msec; matrix 90 × 90 × 60; voxel size 2 × 2 × 2.5 mm3 ; b-values: 0 and 700 s/mm2 . Protocol II: Stejskal-Tanner spin-echo echo-planar imaging with: TE: 88 msec; TR: 10600/11800 msec, matrix 90 × 90 × 60; voxel size 2 × 2 × 2.5 mm3 ; b-values [0, 200, 600, 1200, 1800, 2400, 3000] s/mm2 . ASSESSMENT: Area fractions of cellular and collagen content in histologic sections were quantified using whole-slide image analysis and compared with the corresponding DWI parameters. STATISTICAL TESTS: Correlations were assessed using Pearson's r. Univariate analysis of group median values was done using the Mann-Whitney U-test. RESULTS: Collagen content correlated with the fast signal fraction (r = 0.67, P < 0.001) and ADC (r = 0.58, P < 0.001) and was lower (P < 0.05) in malignant lesions than benign and normal tissues. Cellular content correlated inversely with the fast signal fraction (r = -0.67, P < 0.001) and ADC (r = -0.61, P < 0.001) and was different (P < 0.05) between malignant, benign, and normal tissues. DATA CONCLUSION: Our findings suggest stromal collagen content increases diffusivity observed by MRI and is associated with higher ADC and fast signal fraction of the biexponential model. LEVEL OF EVIDENCE: 3 Technical Efficacy Stage: 3 J. Magn. Reson. Imaging 2020;51:1868-1878.

Accuracy of breast cancer lesion classification using intravoxel incoherent motion diffusion-weighted imaging is improved by the inclusion of global or local prior knowledge with bayesian methods.

BACKGROUND: Diffusion-weighted MRI (DWI) has potential to noninvasively characterize breast cancer lesions; models such as intravoxel incoherent motion (IVIM) provide pseudodiffusion parameters that reflect tissue perfusion, but are dependent on the details of acquisition and analysis strategy. PURPOSE: To examine the effect of fitting algorithms, including conventional least-squares (LSQ) and segmented (SEG) methods as well as Bayesian methods with global shrinkage (BSP) and local spatial (FBM) priors, on the power of IVIM parameters to differentiate benign and malignant breast lesions. STUDY TYPE: Prospective patient study. SUBJECTS: 61 patients with confirmed breast lesions. FIELD STRENGTH/SEQUENCE: DWI (bipolar SE-EPI, 13 b values) was included in a clinical MR protocol including T2 -weighted and dynamic contrast-enhanced MRI on a 3T scanner. ASSESSMENT: The IVIM model was fitted voxelwise in lesion regions of interest (ROIs), and derived parameters were compared across methods within benign and malignant subgroups (correlation, coefficients of variation). Area under receiver operator characteristic curves (ROC AUCs) were calculated to determine discriminatory power of parameter combinations from all fitting methods. STATISTICAL TESTS: Kruskal-Wallis, Mann-Whitney, Pearson correlation. RESULTS: All methods provided useful IVIM parameters; D was well-correlated across all methods (r > 0.8), with a wider range for f and D* (0.3-0.7). Fitting methods gave detectable differences in parameters, but all showed increased f and decreased D in malign lesions. D was the most discriminatory single parameter, with LSQ performing least well (AUC 0.83). In general, ROC AUCs were maximized by the inclusion of pseudodiffusion parameters, and by the use of Bayesian methods incorporating prior information (maximum AUC of 0.92 for BSP). DATA CONCLUSION: DWI performs well at classifying breast lesions, but careful consideration of analysis procedure can improve performance. D is the most discriminatory single parameter, but including pseudodiffusion parameters (f and D*) increases ROC AUC. Bayesian methods outperformed conventional least-squares and segmented fitting methods for breast lesion classification. LEVEL OF EVIDENCE: 3 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;50:1478-1488.

Support vector machine for breast cancer classification using diffusion-weighted MRI histogram features: Preliminary study.

BACKGROUND: Diffusion-weighted MRI (DWI) is currently one of the fastest developing MRI-based techniques in oncology. Histogram properties from model fitting of DWI are useful features for differentiation of lesions, and classification can potentially be improved by machine learning. PURPOSE: To evaluate classification of malignant and benign tumors and breast cancer subtypes using support vector machine (SVM). STUDY TYPE: Prospective. SUBJECTS: Fifty-one patients with benign (n = 23) and malignant (n = 28) breast tumors (26 ER+, whereof six were HER2+). FIELD STRENGTH/SEQUENCE: Patients were imaged with DW-MRI (3T) using twice refocused spin-echo echo-planar imaging with echo time / repetition time (TR/TE) = 9000/86 msec, 90 × 90 matrix size, 2 × 2 mm in-plane resolution, 2.5 mm slice thickness, and 13 b-values. ASSESSMENT: Apparent diffusion coefficient (ADC), relative enhanced diffusivity (RED), and the intravoxel incoherent motion (IVIM) parameters diffusivity (D), pseudo-diffusivity (D*), and perfusion fraction (f) were calculated. The histogram properties (median, mean, standard deviation, skewness, kurtosis) were used as features in SVM (10-fold cross-validation) for differentiation of lesions and subtyping. STATISTICAL TESTS: Accuracies of the SVM classifications were calculated to find the combination of features with highest prediction accuracy. Mann-Whitney tests were performed for univariate comparisons. RESULTS: For benign versus malignant tumors, univariate analysis found 11 histogram properties to be significant differentiators. Using SVM, the highest accuracy (0.96) was achieved from a single feature (mean of RED), or from three feature combinations of IVIM or ADC. Combining features from all models gave perfect classification. No single feature predicted HER2 status of ER + tumors (univariate or SVM), although high accuracy (0.90) was achieved with SVM combining several features. Importantly, these features had to include higher-order statistics (kurtosis and skewness), indicating the importance to account for heterogeneity. DATA CONCLUSION: Our findings suggest that SVM, using features from a combination of diffusion models, improves prediction accuracy for differentiation of benign versus malignant breast tumors, and may further assist in subtyping of breast cancer. LEVEL OF EVIDENCE: 3 Technical Efficacy: Stage 3 J. Magn. Reson. Imaging 2018;47:1205-1216.

Relative enhanced diffusivity: noise sensitivity, protocol optimization, and the relation to intravoxel incoherent motion.

OBJECTIVE: To explore the relationship between relative enhanced diffusivity (RED) and intravoxel incoherent motion (IVIM), as well as the impact of noise and the choice of intermediate diffusion weighting (b value) on the RED parameter. MATERIALS AND METHODS: A mathematical derivation was performed to cast RED in terms of the IVIM parameters. Noise analysis and b value optimization was conducted by using Monte Carlo calculations to generate diffusion-weighted imaging data appropriate to breast and liver tissue at three different signal-to-noise ratios. RESULTS: RED was shown to be approximately linearly proportional to the IVIM parameter f, inversely proportional to D and to follow an inverse exponential decay with respect to D*. The choice of intermediate b value was shown to be important in minimizing the impact of noise on RED and in maximizing its discriminatory power. RED was shown to be essentially a reparameterization of the IVIM estimates for f and D obtained with three b values. CONCLUSION: RED imaging in the breast and liver should be performed with intermediate b values of 100 and 50 s/mm2, respectively. Future clinical studies involving RED should also estimate the IVIM parameters f and D using three b values for comparison.

Segmentation of pre- and posttreatment diffuse glioma tissue subregions including resection cavities.

Background: Evaluating longitudinal changes in gliomas is a time-intensive process with significant interrater variability. Automated segmentation could reduce interrater variability and increase workflow efficiency for assessment of treatment response. We sought to evaluate whether neural networks would be comparable to expert assessment of pre- and posttreatment diffuse gliomas tissue subregions including resection cavities. Methods: A retrospective cohort of 647 MRIs of patients with diffuse gliomas (average 55.1 years; 29%/36%/34% female/male/unknown; 396 pretreatment and 251 posttreatment, median 237 days post-surgery) from 7 publicly available repositories in The Cancer Imaging Archive were split into training (536) and test/generalization (111) samples. T1, T1-post-contrast, T2, and FLAIR images were used as inputs into a 3D nnU-Net to predict 3 tumor subregions and resection cavities. We evaluated the performance of networks trained on pretreatment training cases (Pre-Rx network), posttreatment training cases (Post-Rx network), and both pre- and posttreatment cases (Combined networks). Results: Segmentation performance was as good as or better than interrater reliability with median dice scores for main tumor subregions ranging from 0.82 to 0.94 and strong correlations between manually segmented and predicted total lesion volumes (0.94 < R 2 values < 0.98). The Combined network performed similarly to the Pre-Rx network on pretreatment cases and the Post-Rx network on posttreatment cases with fewer false positive resection cavities (7% vs 59%). Conclusions: Neural networks that accurately segment pre- and posttreatment diffuse gliomas have the potential to improve response assessment in clinical trials and reduce provider burden and errors in measurement.

Discrimination of Breast Cancer from Healthy Breast Tissue Using a Three-component Diffusion-weighted MRI Model.

PURPOSE: Diffusion-weighted MRI (DW-MRI) is a contrast-free modality that has demonstrated ability to discriminate between predefined benign and malignant breast lesions. However, how well DW-MRI discriminates cancer from all other breast tissue voxels in a clinical setting is unknown. Here we explore the voxelwise ability to distinguish cancer from healthy breast tissue using signal contributions from the newly developed three-component multi-b-value DW-MRI model. EXPERIMENTAL DESIGN: Patients with pathology-proven breast cancer from two datasets (n = 81 and n = 25) underwent multi-b-value DW-MRI. The three-component signal contributions C 1 and C 2 and their product, C 1 C 2, and signal fractions F 1, F 2, and F 1 F 2 were compared with the image defined on maximum b-value (DWI max), conventional apparent diffusion coefficient (ADC), and apparent diffusion kurtosis (K app). The ability to discriminate between cancer and healthy breast tissue was assessed by the false-positive rate given a sensitivity of 80% (FPR80) and ROC AUC. RESULTS: Mean FPR80 for both datasets was 0.016 [95% confidence interval (CI), 0.008-0.024] for C 1 C 2, 0.136 (95% CI, 0.092-0.180) for C 1, 0.068 (95% CI, 0.049-0.087) for C 2, 0.462 (95% CI, 0.425-0.499) for F 1 F 2, 0.832 (95% CI, 0.797-0.868) for F 1, 0.176 (95% CI, 0.150-0.203) for F 2, 0.159 (95% CI, 0.114-0.204) for DWI max, 0.731 (95% CI, 0.692-0.770) for ADC, and 0.684 (95% CI, 0.660-0.709) for K app. Mean ROC AUC for C 1 C 2 was 0.984 (95% CI, 0.977-0.991). CONCLUSIONS: The C 1 C 2 parameter of the three-component model yields a clinically useful discrimination between cancer and healthy breast tissue, superior to other DW-MRI methods and obliviating predefining lesions. This novel DW-MRI method may serve as noncontrast alternative to standard-of-care dynamic contrast-enhanced MRI.

Multi-modal characterization of vasculature and nanoparticle accumulation in five tumor xenograft models.

Preclinical research has demonstrated that nanoparticles and macromolecules can accumulate in solid tumors due to the enhanced permeability and retention effect. However, drug loaded nanoparticles often fail to show increased efficacy in clinical trials. A better understanding of how tumor heterogeneity affects nanoparticle accumulation could help elucidate this discrepancy and help in patient selection for nanomedicine therapy. Here we studied five human tumor models with varying morphology and evaluated the accumulation of 100 nm polystyrene nanoparticles. Each tumor model was characterized in vivo using micro-computed tomography, contrast-enhanced ultrasound and diffusion-weighted and dynamic contrast-enhanced magnetic resonance imaging. Ex vivo, the tumors were sectioned for both fluorescence microscopy and histology. Nanoparticle uptake and distribution in the tumors were generally heterogeneous. Density of functional blood vessels measured by fluorescence microscopy correlated significantly (p = 0.0056) with nanoparticle accumulation and interestingly, inflow of microbubbles measured with ultrasound also showed a moderate but significant (p = 0.041) correlation with nanoparticle accumulation indicating that both amount of vessels and vessel morphology and perfusion predict nanoparticle accumulation. This indicates that blood vessel characterization using contrast-enhanced ultrasound imaging or other methods could be valuable for patient stratification for treatment with nanomedicines.