- Research article
- Open Access
Diffusion tensor imaging for characterizing tumor microstructure and improving diagnostic performance on breast MRI: a prospective observational study
Breast Cancer Research volume 21, Article number: 102 (2019)
Diffusion-weighted imaging (DWI) can increase breast MRI diagnostic specificity due to the tendency of malignancies to restrict diffusion. Diffusion tensor imaging (DTI) provides further information over conventional DWI regarding diffusion directionality and anisotropy. Our study evaluates DTI features of suspicious breast lesions detected on MRI to determine the added diagnostic value of DTI for breast imaging.
With IRB approval, we prospectively enrolled patients over a 3-year period who had suspicious (BI-RADS category 4 or 5) MRI-detected breast lesions with histopathological results. Patients underwent multiparametric 3 T MRI with dynamic contrast-enhanced (DCE) and DTI sequences. Clinical factors (age, menopausal status, breast density, clinical indication, background parenchymal enhancement) and DCE-MRI lesion parameters (size, type, presence of washout, BI-RADS category) were recorded prospectively by interpreting radiologists. DTI parameters (apparent diffusion coefficient [ADC], fractional anisotropy [FA], axial diffusivity [λ1], radial diffusivity [(λ2 + λ3)/2], and empirical difference [λ1 − λ3]) were measured retrospectively. Generalized estimating equations (GEE) and least absolute shrinkage and selection operator (LASSO) methods were used for univariate and multivariate logistic regression, respectively. Diagnostic performance was internally validated using the area under the curve (AUC) with bootstrap adjustment.
The study included 238 suspicious breast lesions (95 malignant, 143 benign) in 194 women. In univariate analysis, lower ADC, axial diffusivity, and radial diffusivity were associated with malignancy (OR = 0.37–0.42 per 1-SD increase, p < 0.001 for each), as was higher FA (OR = 1.45, p = 0.007). In multivariate analysis, LASSO selected only ADC (OR = 0.41) as a predictor for a DTI-only model, while both ADC (OR = 0.41) and FA (OR = 0.88) were selected for a model combining clinical and imaging parameters. Post-hoc analysis revealed varying association of FA with malignancy depending on the lesion type. The combined model (AUC = 0.81) had a significantly better performance than Clinical/DCE-MRI-only (AUC = 0.76, p < 0.001) and DTI-only (AUC = 0.75, p = 0.002) models.
DTI significantly improves diagnostic performance in multivariate modeling. ADC is the most important diffusion parameter for distinguishing benign and malignant breast lesions, while anisotropy measures may help further characterize tumor microstructure and microenvironment.
Dynamic contrast-enhanced (DCE) breast MRI is established to be the most sensitive tool for the detection of breast cancer [1,2,3,4,5]. As a result, its utility in both screening and diagnostic setting has increased rapidly over the past two decades. Although breast MRI specificity and positive predictive value have improved since its inception, DCE MRI continues to result in many false positives and unnecessary biopsies [2, 6, 7]. In fact, recent studies have demonstrated that as few as one in five biopsy recommendations based on DCE MRI yield malignancy [8,9,10]. These unnecessary biopsies can result in increased health care costs, patient anxiety, and delays in breast cancer treatment. Accordingly, identification of MRI methods that complement DCE techniques and improve breast MRI specificity without lowering its sensitivity is an important area of active research.
Diffusion-weighted imaging (DWI) has emerged as an adjunct to DCE-MRI that can improve the detection and characterization of breast cancer [11,12,13]. DWI interrogates the in vivo mobility of water molecules, which in turn may provide information on microstructural characteristics of tissue, including cell density and presence of macromolecules and cell membranes. Numerous prior studies have shown that breast cancers feature impeded diffusion and appear as the areas of hyperintensity on DWI with correspondingly low apparent diffusion coefficient (ADC) compared to normal fibroglandular tissue . The best explored application of DWI for breast cancer is decreasing the false-positive rate and increasing the diagnostic specificity when used in addition to conventional DCE-MRI [15, 16]. Although DWI is prone to susceptibility and field inhomogeneity artifacts, diffusion sequences can be acquired quickly and can help to identify lesions warranting biopsy .
Diffusion tensor imaging (DTI) is an extension of conventional DWI that interrogates water motion in six or more directions to characterize diffusion directionality (anisotropy) in addition to ADC . Water diffusion within biological tissue is often anisotropic due to directionally dependent restriction imposed by microstructural architecture. It is hypothesized that normal mammary ducts allow water to diffuse more freely in a direction parallel to the walls of the ducts whereas proliferating neoplastic cells reduce diffusion anisotropy by blocking ducts [18, 19]. Although there is a general consensus that malignant lesions demonstrate reduced diffusion on DWI in comparison with most benign and normal fibroglandular tissues, there are conflicting results regarding the added diagnostic utility of DTI parameters such as fractional anisotropy (FA), a measure of diffusion directionality in which higher FA indicates more anisotropic diffusion oriented along a single direction and lower FA indicates more equal diffusion in all directions [18, 20,21,22,23,24]. Some studies have reported lower FA in benign lesions compared with malignant lesions, attributed to the differences in microscopic composition and organization [20, 23,24,25], while others found no such differences in FA [18, 21, 22].
The purpose of our study was to evaluate the DTI features of suspicious breast lesions detected on 3 T MRI and to determine whether DTI can statistically improve diagnostic performance over conventional assessment.
This prospective study was approved by our Institutional Review Board and was compliant with the Health Insurance Portability and Accountability Act (HIPAA). All patients provided informed consent allowing us to review the MRI images, medical records, and pathology results. Enrolled patients were 18 years or older and underwent 3 T breast MRI, including DCE and DTI sequences, from October 2010 to December 2013. DTI sequences were appended to the standard clinical MRI examination; therefore, enrolled patients did not undergo additional MRI examinations. Patients with MRI-detected lesions characterized as Breast Imaging Reporting and Data System (BI-RADS) category 4 or 5 who underwent core needle biopsy (CNB) and/or surgical excision were eligible for the study. Clinical indications for breast MRI included high-risk screening, extent of disease evaluation, and problem solving. In subjects with known existing breast cancer, the eligible lesion must have been distinct from the previously biopsy-proven cancer. Subjects receiving neoadjuvant chemotherapy less than 6 months prior to MRI were excluded. Patients unable or unwilling to provide informed consent or undergo the entire MRI examination were also excluded.
Breast MRIs were performed using a 3 T Philips Achieva Tx MRI scanner (Philips Healthcare, Best, The Netherlands) with a dedicated 16-channel bilateral breast coil (Mammo-Trak, Philips Healthcare, Best, The Netherlands). MRI sequences were obtained in the axial orientation, and each MRI exam included DWI, T2-weighted fast spin-echo, T1-weighted non-fat-suppressed, and T1-weighted fat-suppressed DCE-MRI sequences with one precontrast and three post-contrast acquisitions. DCE-MRI was acquired with T1-weighted fat-suppressed 3D fast gradient echo (eTHRIVE) sequences with parallel imaging technique (sensitivity encoding; SENSE). The following imaging parameters were utilized: repetition time (TR)/echo time (TE), 5.9/3 ms; flip angle, 10°; matrix size, 440 × 660; field of view (FOV), 22 × 33 cm; in-plane voxel size, 0.5 mm; and slice thickness, 1.3 mm. Post-contrast sequences were acquired with k-space centered at 120, 300, and 480 s after contrast injection. The contrast agent administered was 0.1 mmol/kg body weight gadoteridol (ProHance, Bracco Diagnostics, Milan, Italy). The DCE scan time was 2 min and 57 s per acquisition.
In order to minimize disruption of the clinical portion of breast MRI examinations, diffusion tensor imaging was performed immediately following DCE imaging using a 2D diffusion-weighted single-shot spin-echo-prepared echo-planar imaging (EPI) sequence with parallel imaging and fat suppression (spectral attenuated inversion recovery (SPAIR)) with the following parameters: SENSE reduction factor, 3; averages, 2; TR/TE, 5336/61 ms; matrix, 240 × 240; FOV, 36 × 36 cm; in-plane voxel size, 1.5 mm; slice thickness, 5 mm; and gap, 0. Diffusion gradients were applied in six directions with b values of 0, 100, and 800 s/mm2. The total diffusion imaging acquisition time was 3 min and 28 s.
Clinical MRI interpretation
Clinical interpretations for all MRI studies were performed prospectively by fellowship-trained radiologists specializing in breast imaging. Lesions were assessed using American College of Radiology (ACR) BI-RADS breast MRI lexicon , and lesion kinetic features were measured using computer-assisted diagnosis (CAD) software (CADstream v. 5.2.7, Merge Healthcare, Chicago, IL). For DCE kinetics, enhancement curve types of persistent, plateau, or washout were categorized for each voxel by evaluating the change in signal intensity from the initial (at 120 s) to the final (at 480 s) post-contrast scan, with washout defined as > 10% decrease in signal intensity, persistent > 10% increase, and plateau < 10% change. Recorded lesion characteristics included the number of lesions per patient, lesion type (focus, mass, non-mass enhancement [NME]), size, DCE kinetic pattern of worst curve type (defined as most suspicious, with washout > plateau > persistent ), final BI-RADS assessment and recommendation.
This information was entered into our clinical MRI database. Because the clinical evaluation of DCE-MRI images was performed prospectively, radiologists were blinded to the result of histopathology at the time of interpretation. DTI data was not reviewed by radiologists for BI-RADS assessment. Histopathology results from CNB and/or excision biopsy were later extracted from the clinical record for the purposes of this study.
DTI analysis was performed offline by trained research scientists who were blinded to the lesion pathology outcomes. Diffusion tensor images were first spatially registered using a commercially available 3D affine transformation algorithm (Diffusion Registration tool, Philips Healthcare, Best, The Netherlands), with b = 0 s/mm2 images as reference, to minimize the artifacts due to motion and eddy current-based image distortion . Voxel-based DTI parametric maps were then calculated and analyzed using in-house custom semi-automated software developed in ImageJ (National Institutes of Health, Bethesda, MD) based on standard methods . With DTI, the MRI signal obtained with diffusion weighting is reduced in intensity proportional to the water mobility and is described by S = S0e–bD, where S0 is the signal intensity without diffusion weighting, S is the signal intensity with diffusion weighting, D is the diffusion tensor, and b is the applied diffusion sensitization . A diffusion tensor matrix is derived for each image voxel and diagonalized to obtain the diffusion tensor eigenvalues λ1, λ2, and λ3 and scalar values describing the magnitude or rate of diffusion along each of the three principal axes (from largest to smallest) of the diffusion tensor ellipsoid (in mm2/s). From those, a number of other rotationally invariant DTI parameters can be calculated. ADC (also known as mean diffusivity, MD, or averaged diffusivity, Dav), which describes the degree of mobility or hindrance of water molecules, is given by:
(an alternative way to calculate ADC from the diffusion tensors vs. conventional DWI that does not account for anisotropy). Mean axial diffusivity is defined as λ1. Mean radial diffusivity is defined as the average of λ2 and λ3.
FA is a unitless measure of the degree of directionality of diffusion, ranging from 0 (completely isotropic) to 1 (completely anisotropic), given by:
Additionally, the empirical parameter λ1 − λ3 is a non-normalized measure related to both diffusion mobility and directionality, which has been reported in previous studies to be significantly lower in malignant lesions compared to normal fibroglandular tissue and benign breast lesions [19, 22]. For qualitative interpretation, combined diffusion-weighted images were also calculated as the geometric average of unidirectional b = 800 s/mm2 images.
A region of interest (ROI) corresponding to the BI-RADS 4 or 5 lesion on DCE-MRI was manually defined on DTI under the supervision of an experienced radiologist. This ROI was selected on the combined DWI at the central slice of the lesion to include any hyperintensity and to avoid obvious areas of cyst, necrosis, or fat by referring to T1- and T2-weighted images. A semi-automated thresholding tool enabled exclusion from the ROI any voxels with very low DWI signal intensity corresponding to adipose or normal fibroglandular tissue . The ROI was propagated to all DTI parametric maps, and the mean voxel value was calculated to characterize the lesion on each map. Propagating lesion ROIs directly from DCE-MRI was not possible due to the common spatial distortions inherent to echo-planar imaging-based DTI datasets (caused by B0 field inhomogeneities). Examples of breast DTI parametric maps are shown for several benign and malignant lesions in Figs. 1, 2, 3, and 4.
An analysis was performed at the lesion level, and the primary outcome was histopathological determination of malignant vs. benign. Univariate associations between malignancy and clinical, DCE-MRI, and DTI parameters were explored using generalized estimating equations (GEE)-based logistic regression to account for any correlation between multiple lesions from the same woman. The area under the receiver operating characteristic (ROC) curve (AUC) and odds ratio (OR) were calculated for each parameter. Any variable which appeared to be right-skewed based on visual inspection of the histogram was log-transformed prior to inclusion in the model. ORs for continuous variables were scaled to correspond to the difference per 1-SD increase in the variable. Spearman’s rank correlation was used to examine the pairwise relationships between continuous variables.
Multivariate logistic regression models were developed using the least absolute shrinkage and selection operator (LASSO), which is a machine learning technique that simultaneously performs variable selection and parameter regularization to limit overfitting . The LASSO penalty/regularization parameter was selected to minimize the model deviance estimated using leave-one-patient-out cross-validation. Three primary models were generated based on (1) Clinical/DCE-MRI parameters only, (2) DTI parameters only, and (3) the combination of both sets of parameters. Due to substantial collinearity among the DTI parameters, the primary DTI-based models were developed considering ADC and FA only. Sensitivity analyses to determine how model performance was impacted by changes in specific input variables were also performed where ADC was replaced with axial and radial diffusivity and where FA was replaced with λ1 − λ3. Model performance for discriminating between malignant and benign lesions was summarized using the AUC after optimism adjustment using the bootstrap to account for training and testing using the same dataset . The bootstrap was also used to compare AUC estimates from different models. Bootstrap resampling was performed by patient rather than by lesion to account for the non-independence of lesions from the same patient.
An exploratory subgroup analysis of DTI parameters within the lesion groups defined by size (< 1 cm and ≥ 1 cm) and type (masses and non-masses) was conducted to assess whether the relationships between the DTI parameters and malignancy differed between any of these groups. This was done using GEE-based logistic regression models with interaction terms corresponding to the comparisons of interest, specifically size × ADC, type × ADC, size × FA, and type × FA. Lastly, a post-hoc multivariate LASSO model was developed, which included all Clinical/DCE-MRI parameters, DTI parameters, and the same interactions between size, type, and DTI parameters. All statistical calculations were performed using the statistical computing language R (version 3.1.1; R Foundation for Statistical Computing, Vienna, Austria). Throughout, two-sided tests were used, with statistical significance defined as p < 0.05.
Patient cohort and lesion characteristics
Two hundred sixty-six women with 354 MRI-detected BI-RADS category 4 or 5 lesions were enrolled in our study prior to CNB and/or surgical excision. Of those, DTI was not performed in the MRI examinations for 23 women (24 lesions) and had technical or image post-processing issues in 13 women (16 lesions). Another 31 women (69 lesions) were then excluded: 42 lesions lacked adequate reference standard to determine the pathologic outcome, and 27 lesions were not evaluable because poor image quality or small lesion size precluded DTI measures. Lastly, 5 women with 7 lesions were excluded for incomplete clinical data: 2 women (2 lesions) in whom BPE was not evaluable due to prior mastectomy and 3 women (5 lesions) who did not have a recent mammogram available to assess breast density. After all exclusions, there were 238 lesions in 194 women included in the analysis.
Patient and lesion characteristics are summarized in Table 1. The median patient age was 51 years (range 23 to 83 years). The clinical indication for breast MRI was screening or problem solving in 68 (35.1%) patients and evaluating the extent of disease for newly diagnosed breast cancer in 128 (64.9%) patients. The majority (80.9%) of patients had a single BI-RADS 4 or 5 lesion detected on MRI while 15.5% had 2 lesions and 3.6% had 3 lesions. The median lesion diameter was 1.1 cm (range 0.4 to 11.4 cm), and most lesions were masses (n = 135, 56.7%), were classified as BI-RADS 4 (n = 218, 91.6%), and demonstrated washout on delayed phase DCE images (n = 198, 83.2%). Ninety-five of 238 (39.9%) lesions were malignant and 143 (60.1%) were benign. Lesion subtype information was available in 90 of 95 malignancies, which is summarized in Table 1.
Performance of DCE-MRI and DTI parameters in discriminating malignant and benign lesions
In univariate analysis, all clinical and DCE-MRI parameters evaluated were significantly associated with malignancy except for mass vs. non-mass lesion type (p = 0.98), as summarized in Table 2. From DTI, lower ADC was significantly associated with malignancy (OR = 0.37 per 1-SD increase, p < 0.001), as were axial (OR = 0.42, p < 0.001) and radial (OR = 0.40, p < 0.001) diffusivity. Higher FA was significantly associated with malignancy (OR = 1.45, p = 0.007), though λ1 − λ3 was not (OR = 0.96, p = 0.77).
There was substantial collinearity among the DTI parameters, at least in part due to their definitions which are all functions of the three diffusion eigenvalues. The pairwise correlations of ADC with axial (r = 0.81) and radial (r = 0.95) diffusivity were high, as was the correlation between FA and λ1 − λ3 (r = 0.89). FA was correlated with ADC (r = − 0.51), primarily through its correlation with radial diffusivity (r = − 0.73) rather than the axial component (r = 0.04). By contrast, λ1 − λ3 was only weakly correlated with ADC (r = − 0.07) because λ1 − λ3 was positively correlated with axial diffusivity (r = 0.48) but negatively correlated with radial diffusivity (r = − 0.36).
The multivariate LASSO models are summarized in Table 3. The LASSO selected all of the Clinical/DCE-MRI parameters except post-menopausal status and selected only ADC (OR = 0.41 per 1-SD increase) for the DTI-only model. In the combined Clinical/DCE-MRI+DTI model, the LASSO selected both ADC (OR = 0.41) and FA (OR = 0.88). The Clinical/DCE-MRI model (AUC = 0.76) and DTI-only model (AUC = 0.75) had similar discrimination performance (ΔAUC = − 0.01, 95% CI − 0.10 to 0.06, p = 0.54). The Clinical/DCE-MRI+DTI model (AUC = 0.81) had significantly better performance than both the Clinical/DCE-MRI model (ΔAUC = 0.05, 95% CI 0.02 to 0.10, p < 0.001) and the DTI-only model (ΔAUC = 0.07, 95% CI 0.03 to 0.13, p = 0.002). Corresponding ROC curves are shown in Fig. 5. Since the OR for FA was relatively modest at 0.88, the incremental impact of FA was further examined. The performance of the Clinical/DCE-MRI+DTI model was compared with a Clinical/DCE-MRI+ADC model (including the same parameters but without FA), and the improvement in AUC from adding FA was not significant overall (ΔAUC = 0.00, 95% CI − 0.01 to 0.02, p = 0.81). The performance results from all of these models were very similar when the models were re-fit using axial and radial diffusivity instead of ADC and using λ1 − λ3 instead of FA.
Subgroup analysis of DTI parameters in small vs. large and mass vs. non-mass lesions
An exploratory subgroup analysis was conducted of small (< 1 cm, n = 98) vs. large (≥ 1 cm, n = 140) lesions and mass (n = 135) vs. non-mass (NME n = 99; foci n = 4) lesions to assess how the DTI parameters were related to malignancy within these lesion subgroups. In univariate modeling, the predictive values of ADC and FA for malignancy were not significantly different between small and large lesions (p = 0.23 and p = 0.43, respectively) but were different between masses and non-masses (p = 0.002 and p = 0.006, respectively), Fig. 6. Specifically, ADC was more strongly predictive of malignancy in masses (OR = 0.19, p < 0.001) than non-masses (OR = 0.64, p = 0.060), while FA was significantly predictive only for masses (OR = 2.02, p < 0.001, Table 4) in univariate analysis.
Since both lesion type and size were related—non-masses (82/103; 80%) were more likely to be large (> 1 cm) than masses (58/135; 43%, p < 0.0001)—and ADC and FA were moderately correlated (r = − 0.51), they were included together in a multivariate model to assess their independent associations with malignancy in mass and non-mass groups (Table 4). The predictive value of FA remained significantly different between masses and non-masses in the multivariate model (OR 1.62 vs. 0.54, p = 0.009). Interestingly, the OR values implied opposite independent associations of FA with malignancy in the two groups, with malignant masses having higher FA values than benign masses but malignant non-masses having lower FA than benign non-masses. By contrast, the ORs for ADC were no longer significantly different between masses and non-masses in the multivariate model (OR 0.23 vs. 0.42, p = 0.20), with lower ADC values associated with malignancy in both lesion type groups.
We further explored why the FA association with malignancy varied between masses and non-masses. Of the malignant lesions, 51 of 54 (94.4%) masses were invasive cancers compared to 22 of 36 (61.1%) non-masses (p < 0.001). FA was significantly lower in DCIS lesions compared to invasive lesions (0.22 vs. 0.30, p = 0.026). On the other hand, FA was significantly higher in benign non-masses than in benign masses (0.26 vs. 0.21, p = 0.040). Among benign lesions, masses were most commonly fibroadenoma (25/83 [30.1%]), while non-masses were most commonly fibrocystic change (15/60 [25.0%]) and usual ductal hyperplasia (10/60 [16.7%]).
To account for differential associations with malignancy depending on the lesion subgroup, an additional LASSO model combining Clinical/DCE-MRI parameters, DTI parameters, and size- and type-DTI interactions was then fit to assess the impact of incorporating the subgroup analysis results into the modeling (Table 5). This model allowed the ORs for ADC and FA to depend on whether the lesion was large or small and whether it was a mass or not. The resulting model with size- and type-DTI interactions (AUC = 0.85, 95% CI 0.82 to 0.90) demonstrated a substantial improvement in discrimination performance compared to the basic Clinical/DCE-MRI model (ΔAUC = 0.09, 95% CI 0.04 to 0.13, p < 0.001) and a modest though statistically significant improvement compared to the Clinical/DCE-MRI+DTI model without interactions shown in Table 3 (ΔAUC = 0.03, 95% CI 0.01 to 0.07, p = 0.018, Fig. 7).
The results of our prospective study of women with suspicious 3 T DCE-MRI-detected BI-RADS 4 and 5 lesions who underwent CNB and/or surgical excision showed that the addition of DTI to conventional breast MRI assessments may improve the ability to distinguish between benign and malignant lesions. Overall, malignancies exhibited lower ADC, lower axial (λ1) and radial ([λ2 + λ3]/2) diffusivity, and higher FA on DTI than benign lesions. Further analysis demonstrated that FA association with malignancy differed for masses and non-masses, and among the diffusion parameters, ADC contributed most significantly to the overall diagnostic performance. Machine learning-based LASSO modeling incorporating multiple clinical, conventional DCE and DTI parameters achieved high diagnostic performance in differentiating malignant from benign lesions (AUC = 0.85).
Prior studies have demonstrated that lesion size and worst delayed phase kinetics on DCE-MRI are significant predictors of malignancy [32,33,34], and the results of our study further validate the significance of these parameters, with malignant lesions being larger and more likely to demonstrate washout on delayed phase in both univariate and multivariate analyses. As expected, BI-RADS category 5 lesions were more likely to be malignant than BI-RADS 4 lesions. Lesion type (mass vs. non-mass) was not found to be a significant parameter. Other patient-level characteristics associated with higher odds for lesion malignancy were older age, post-menopausal status, lower BPE, and lower breast density. These findings agree with established evidence of increasing risk of breast cancer with age , as well as reduced mammographic breast density and BPE that go along with increasing age . Furthermore, higher BPE (moderate or marked levels) has been shown to be associated with higher abnormal interpretation (BI-RADS 0,3,4,5) and biopsy rates and lower specificity (i.e., higher false-positive rate) .
The latest (fifth) edition of the ACR BI-RADS lexicon contains descriptors for lesion morphology and contrast kinetics but not for diffusion characteristics , which reflects the current clinical practice of DWI being used only in select imaging centers as an adjunctive technique to conventional DCE-MRI. Results from multiple single-center studies, including prior work at our institution, have demonstrated the value of quantitative ADC measures in discriminating malignant from benign lesions [15, 37]. Results of a recent multicenter trial further confirmed that by implementing an ADC cutoff, DWI has potential to reduce the rate of unnecessary biopsies prompted by conventional DCE-MRI . The pathologic basis for impeded diffusion and lower ADC values in malignant lesions has been proposed to arise in part from the higher cellularity and more restricted extracellular environment of breast cancers compared to benign lesions, supported by imaging pathologic comparisons in some studies [11, 23]. In our study, univariate analysis identified ADC as a significant parameter, and LASSO modeling selected ADC for both the DTI-only and DCE+DTI models. Our findings add to a growing body of evidence supporting the use of ADC as an independent biomarker for differentiating between malignant and benign lesions.
In addition to the magnitude of diffusion provided by ADC, other DTI parameters provide information on the direction and anisotropy of diffusion, potentially allowing for further characterization of the underlying breast tissue. Although there is a general consensus in the literature supporting impeded diffusion (reflected by low ADC) as a feature of breast malignancies, there are conflicting results regarding the added utility of other DTI parameters in differentiating between malignant and benign lesions. FA is the most studied DTI parameter aside from ADC, and while some studies have reported higher FA in malignant lesions compared with benign ones [20, 23,24,25], others have found no significant difference [18, 21, 22]. The univariate analysis in our study showed mean FA to be significantly higher in malignant lesions compared to benign lesions. Mean axial (λ1) and radial ([λ2 + λ3]/2) diffusivity were also significantly lower in malignancies. On the other hand, the empirical parameter λ1 − λ3, a proposed alternate measure of diffusion directionality that has been reported to be lower in malignant lesions compared to normal breast tissue and benign lesions [19, 22], did not significantly distinguish malignant and benign lesions in our study.
Post-hoc subanalyses demonstrated interesting differences in FA association with malignancy based on the lesion type of mass or non-mass. Higher FA was associated with malignancy for masses, while lower FA was associated with malignancy for non-mass lesions in multivariate modeling, which has not been previously reported. Our findings suggest this opposite association of FA with malignancy likely relates to biologic differences of the typical benign pathologies represented within each lesion type. In our study, benign masses demonstrated lower FA than benign non-masses. Within masses, the most common false positives on MRI are fibroadenomas, where FA has been reported to be lower than for malignancies in multiple prior studies in addition to ours [20, 23,24,25] and attributed to their characteristic myxoid extracellular matrix and fibrous stroma with absent or compressed tubular structures . On the other hand, within non-masses, common benign pathologies such as ductal hyperplasia or pseudoangiomatous stromal hyperplasia (PASH) may grow diffusely while maintaining some of the normal ductal architecture and native FA levels of the intervening fibroglandular tissue. Further investigation with more detailed pathological assessment is needed to better understand the microstructural and microenvironmental characteristics of breast lesions influencing diffusion anisotropy measures and how this information may be used to improve diagnostic accuracy. Regardless, our findings suggest DTI anisotropy metrics must be considered in the context with lesion type for diagnostic purposes.
The results of our multivariate LASSO analysis showed that a model incorporating all clinical factors (except menopausal status), DCE-MRI parameters, and DTI parameters (ADC and FA) and accounting for interactions achieved the best diagnostic performance in differentiating malignant from benign lesions, represented the highest AUC (0.85). Comparison of several multivariate models further demonstrated that among the DTI parameters, ADC contributed most significantly to the overall diagnostic performance and suggested FA added incremental value only after accounting for the interactions with lesion type. ADC maps can be acquired using the standard DWI technique with three orthogonal diffusion-sensitizing gradients, while multiple additional gradient directions must be applied to acquire sufficient data to reconstruct diffusion tensors in DTI. Our results suggest that using only the less time-consuming standard DWI acquisition in conjunction with conventional DCE-MRI may be sufficient for improving the diagnostic performance.
Along with demonstrating potential clinical utility, the study also identified that further technical developments are needed to address breast DTI image quality issues. Lesion evaluability on DTI was limited in part by technical issues inherent to the single-shot echo-planar imaging (EPI) technique. EPI is widely used for diffusion imaging but suffers from limited spatial resolution, spatial distortion, and frequent artifacts . These issues are further magnified for breast imaging due to the particular challenges of off-isocenter imaging, air-tissue interfaces, and significant fat content in the breast . To mitigate these potential sources of error, exams found on visual assessment to have significant imaging artifacts were excluded from the study, but these issues may still have contributed to reduced lesion evaluability in the remaining exams, with 27/272 (10%) lesions with adequate DTI quality and reference standard still deemed non-evaluable. The results of recent multicenter breast DWI trials also identified reliable image quality to be a challenge, with 13 to 29% of cases being excluded for DWI technical issues [16, 39, 40]. However, a range of emerging technical advancements in DWI acquisition strategies hold potential to improve image quality [41,42,43,44].
There are several strengths of our study. To our knowledge, our study contains the largest prospective cohort of patients to date who have undergone breast DTI for detecting breast malignancy. After excluding lesions without direct pathologic sampling and with other missing clinical data, 238 lesions included in our study were available for statistical analysis. Our study featured a prospective design, and interpreting radiologists and research associates conducting quantitative DTI analyses were blinded to the histopathologic results. All breast MRIs in our study were performed at high-field strength (3 T), which has been shown in prior studies to have improved diagnostic performance compared to 1.5 T due to superior signal-to-noise ratio (SNR), which can in turn enable increased spatial resolution [45,46,47]. To increase the robustness of our statistical analysis, we used an advanced machine learning-based LASSO modeling technique, which simultaneously performs variable selection and parameter regularization to limit overfitting, along with bootstrap optimism adjustment to account for training and testing using the same dataset.
There are also several important limitations to our study. First, although our study contains the largest cohort of breast DTI patients to date, it was performed at a single institution, which may limit the generalizability of our results. DTI was performed with minimal six gradient directions and with 1.5 × 1.5 × 5 mm spatial resolution. It is possible additional directions and/or higher spatial resolution could better elucidate directionality in breast lesions; however, this would require additional acquisition time and may incur additional motion and other artifacts. In order to minimize disruption of the clinical portion of breast MRI examinations and in case patients were not able to tolerate the full examination, DTI was performed after the completion of the standard MRI protocol including DCE-MRI (approximately 10 min after contrast injection). Although we have previously shown no significant effect on breast tumor ADC measures using our imaging protocol , which was also verified across multiple independent breast DWI studies , it may be preferable to acquire DTI before contrast injection to avoid any possible confounding effects on other DTI parameters. Furthermore, ROIs were manually defined on DTI maps for each exam after comparing DCE-MRI and DTI. Propagating lesion ROIs directly from DCE images was not reliable due to the spatial distortions common to EPI-based breast DWI and DTI datasets (caused by field inhomogeneities). Manual ROI definition is prone to operator dependence and sampling error, especially for irregularly shaped masses or NME. To address this limitation, we used a semi-automated ROI tool to avoid non-tumor voxels, which has been shown in a previous study to improve the inter-reader reproducibility of breast lesion ADC values without introducing bias vs. manual ROI measures . Alternate ROI approaches such as sampling a small hotspot region could provide different performance results for distinguishing benign from malignant lesions.
Evaluating a combination of clinical, DCE-MRI, and diffusion parameters may improve the ability to distinguish between benign and malignant lesions on breast MRI, thereby decreasing false-positive diagnoses and avoiding unnecessary biopsies. ADC was the most important diffusion parameter for distinguishing benign and malignant breast lesions, supporting continued use of standard DWI sequences for feasible clinical implementation. However, our results suggest DTI may enable further biologic characterization relating to variations in tumor microstructure and microenvironment, which warrants further investigation.
Availability of data and materials
The datasets generated and analyzed during the current study are not publicly available due to identifiable patient information but are available from the corresponding author on reasonable request.
Dynamic contrast enhanced
Apparent diffusion coefficient
Diffusion tensor imaging
Health Insurance Portability and Accountability Act
Breast Imaging Reporting and Data System
Core needle biopsy
Field of view
American College of Radiology
Region of interest
Generalized estimating equations
Receiver operating characteristic
Area under the receiver operating characteristic curve
Least absolute shrinkage and selection operator
Riedl CC, Luft N, Bernhart C, Weber M, Bernathova M, Tea MK, Rudas M, Singer CF, Helbich TH. Triple-modality screening trial for familial breast cancer underlines the importance of magnetic resonance imaging and questions the role of mammography and ultrasound regardless of patient mutation status, age, and breast density. J Clin Oncol. 2015;33(10):1128–35.
Lehman CD, Isaacs C, Schnall MD, Pisano ED, Ascher SM, Weatherall PT, Bluemke DA, Bowen DJ, Marcom PK, Armstrong DK, et al. Cancer yield of mammography, MR, and US in high-risk women: prospective multi-institution breast cancer screening study. Radiology. 2007;244(2):381–8.
Kuhl C, Weigel S, Schrading S, Arand B, Bieling H, Konig R, Tombach B, Leutner C, Rieber-Brambs A, Nordhoff D, et al. Prospective multicenter cohort study to refine management recommendations for women at elevated familial risk of breast cancer: the EVA trial. J Clin Oncol. 2010;28(9):1450–7.
Sardanelli F, Podo F, D’Agnolo G, Verdecchia A, Santaquilani M, Musumeci R, Trecate G, Manoukian S, Morassut S, de Giacomi C, et al. Multicenter comparative multimodality surveillance of women at genetic-familial high risk for breast cancer (HIBCRIT study): interim results. Radiology. 2007;242(3):698–715.
Berg WA, Zhang Z, Lehrer D, Jong RA, Pisano ED, Barr RG, Bohm-Velez M, Mahoney MC, Evans WP 3rd, Larsen LH, et al. Detection of breast cancer with addition of annual screening ultrasound or a single screening MRI to mammography in women with elevated breast cancer risk. Jama. 2012;307(13):1394–404.
Elmore JG, Armstrong K, Lehman CD, Fletcher SW. Screening for breast cancer. Jama. 2005;293(10):1245–56.
Peters NH, Borel Rinkes IH, Zuithoff NP, Mali WP, Moons KG, Peeters PH. Meta-analysis of MR imaging in the diagnosis of breast lesions. Radiology. 2008;246(1):116–24.
Lee JM, Ichikawa L, Valencia E, Miglioretti DL, Wernli K, Buist DSM, Kerlikowske K, Henderson LM, Sprague BL, Onega T, et al. Performance benchmarks for screening breast MR imaging in community practice. Radiology. 2017;285(1):44–52.
Niell BL, Gavenonis SC, Motazedi T, Chubiz JC, Halpern EP, Rafferty EA, Lee JM. Auditing a breast MRI practice: performance measures for screening and diagnostic breast MRI. J Am Coll Radiol. 2014;11(9):883–9.
Strigel RM, Rollenhagen J, Burnside ES, Elezaby M, Fowler AM, Kelcz F, Salkowski L, DeMartini WB. Screening breast MRI outcomes in routine clinical practice: comparison to BI-RADS benchmarks. Acad Radiol. 2017;24(4):411–7.
Guo Y, Cai YQ, Cai ZL, Gao YG, An NY, Ma L, Mahankali S, Gao JH. Differentiation of clinically benign and malignant breast lesions using diffusion-weighted imaging. J Magn Reson Imaging. 2002;16(2):172–8.
Kinoshita T, Yashiro N, Ihara N, Funatu H, Fukuma E, Narita M. Diffusion-weighted half-Fourier single-shot turbo spin echo imaging in breast tumors: differentiation of invasive ductal carcinoma from fibroadenoma. J Comput Assist Tomogr. 2002;26(6):1042–6.
Sinha S, Lucas-Quesada FA, Sinha U, DeBruhl N, Bassett LW. In vivo diffusion-weighted MRI of the breast: potential for lesion characterization. J Magn Reson Imaging. 2002;15(6):693–704.
Partridge SC, Nissan N, Rahbar H, Kitsch AE, Sigmund EE. Diffusion-weighted breast MRI: clinical applications and emerging techniques. J Magn Reson Imaging. 2017;45(2):337–55.
Zhang L, Tang M, Min Z, Lu J, Lei X, Zhang X. Accuracy of combined dynamic contrast-enhanced magnetic resonance imaging and diffusion-weighted imaging for breast cancer detection: a meta-analysis. Acta Radiol. 2016;57(6):651–60.
Rahbar H, Zhang Z, Chenevert TL, Romanoff J, Kitsch A, Hanna LG, Harvey SM, Moy L, DeMartini WB, Dogan BE, Yang W, Wang L, Joe BN, Oh KY, Neal CH, McDonald ES, Schnall M, Lehman C, Comstock C, Partridge SC. Utility of Diffusion Weighted Imaging to Decrease Unnecessary Biopsies Prompted by Breast MRI: A Trial of the ECOG-ACRIN Cancer Research Group (A6702). Clin Cancer Res. 2019;25(6):1756-65. https://0-doi-org.brum.beds.ac.uk/10.1158/1078-0432.CCR-18-2967. Epub 2019 Jan 15. PMID: 30647080.
Le Bihan D, Mangin JF, Poupon C, Clark CA, Pappata S, Molko N, Chabriat H. Diffusion tensor imaging: concepts and applications. J Magn Reson Imaging. 2001;13(4):534–46.
Partridge SC, Ziadloo A, Murthy R, White SW, Peacock S, Eby PR, DeMartini WB, Lehman CD. Diffusion tensor MRI: preliminary anisotropy measures and mapping of breast tumors. J Magn Reson Imaging. 2010;31(2):339–47.
Furman-Haran E, Grobgeld D, Nissan N, Shapiro-Feinberg M, Degani H. Can diffusion tensor anisotropy indices assist in breast cancer detection? J Magn Reson Imaging. 2016;44(6):1624–32.
Baltzer PA, Schafer A, Dietzel M, Grassel D, Gajda M, Camara O, Kaiser WA. Diffusion tensor magnetic resonance imaging of the breast: a pilot study. Eur Radiol. 2011;21(1):1–10.
Cakir O, Arslan A, Inan N, Anik Y, Sarisoy T, Gumustas S, Akansel G. Comparison of the diagnostic performances of diffusion parameters in diffusion weighted imaging and diffusion tensor imaging of breast lesions. Eur J Radiol. 2013;82(12):e801–6.
Eyal E, Shapiro-Feinberg M, Furman-Haran E, Grobgeld D, Golan T, Itzchak Y, Catane R, Papa M, Degani H. Parametric diffusion tensor imaging of the breast. Investig Radiol. 2012;47(5):284–91.
Jiang R, Ma Z, Dong H, Sun S, Zeng X, Li X. Diffusion tensor imaging of breast lesions: evaluation of apparent diffusion coefficient and fractional anisotropy and tissue cellularity. Br J Radiol. 2016;89(1064):20160076.
Teruel JR, Goa PE, Sjobakk TE, Ostlie A, Fjosne HE, Bathen TF. Diffusion weighted imaging for the differentiation of breast tumors: from apparent diffusion coefficient to high order diffusion tensor imaging. J Magn Reson Imaging. 2016;43(5):1111–21.
Tsougos I, Svolos P, Kousi E, Athanassiou E, Theodorou K, Arvanitis D, Fezoulidis I, Vassiou K. The contribution of diffusion tensor imaging and magnetic resonance spectroscopy for the differentiation of breast lesions at 3T. Acta Radiol. 2014;55(1):14–23.
Morris EA, Comstock CE, Lee CH, et al. ACR BI-RADS® magnetic resonance imaging. In: ACR BI-RADS® Atlas, Breast Imaging Reporting and Data System. Reston: American College of Radiology; 2013.
Arlinghaus LR, Welch EB, Chakravarthy AB, Xu L, Farley JS, Abramson VG, Grau AM, Kelley MC, Mayer IA, Means-Powell JA, et al. Motion correction in diffusion-weighted MRI of the breast at 3T. J Magn Reson Imaging. 2011;33(5):1063–70.
Basser PJ, Pierpaoli C. Microstructural and physiological features of tissues elucidated by quantitative-diffusion-tensor MRI. J Magn Reson B. 1996;111(3):209–19.
Rahbar H, Kurland BF, Olson ML, Kitsch AE, Scheel JR, Chai X, Usoro J, Lehman CD, Partridge SC. Diffusion-weighted breast magnetic resonance imaging: a semiautomated voxel selection technique improves interreader reproducibility of apparent diffusion coefficient measurements. J Comput Assist Tomogr. 2016;40(3):428–35.
Tibshirani R. Regression shrinkage and selection via the Lasso. J R Stat Soc Ser B. 1996;58(1):267–88.
Steyerberg EW, Harrell FE Jr, Borsboom GJ, Eijkemans MJ, Vergouwe Y, Habbema JD. Internal validation of predictive models: efficiency of some procedures for logistic regression analysis. J Clin Epidemiol. 2001;54(8):774–81.
Wang LC, DeMartini WB, Partridge SC, Peacock S, Lehman CD. MRI-detected suspicious breast lesions: predictive values of kinetic features measured by computer-aided evaluation. AJR Am J Roentgenol. 2009;193(3):826–31.
Liberman L, Mason G, Morris EA, Dershaw DD. Does size matter? Positive predictive value of MRI-detected breast lesions as a function of lesion size. AJR Am J Roentgenol. 2006;186(2):426–30.
Demartini WB, Kurland BF, Gutierrez RL, Blackmore CC, Peacock S, Lehman CD. Probability of malignancy for lesions detected on breast MRI: a predictive model incorporating BI-RADS imaging features and patient characteristics. Eur Radiol. 2011;21(8):1609–17.
King V, Gu Y, Kaplan JB, Brooks JD, Pike MC, Morris EA. Impact of menopausal status on background parenchymal enhancement and fibroglandular tissue on breast MRI. Eur Radiol. 2012;22(12):2641–7.
Ray KM, Kerlikowske K, Lobach IV, Hofmann MB, Greenwood HI, Arasu VA, Hylton NM, Joe BN. Effect of background parenchymal enhancement on breast MR imaging interpretive performance in community-based practices. Radiology. 2018;286(3):822–9.
Chen X, Li WL, Zhang YL, Wu Q, Guo YM, Bai ZL. Meta-analysis of quantitative diffusion-weighted MR imaging in the differential diagnosis of breast lesions. BMC Cancer. 2010;10:693.
Le Bihan D, Poupon C, Amadon A, Lethimonnier F. Artifacts and pitfalls in diffusion MRI. J Magn Reson Imaging. 2006;24(3):478–88.
Newitt DC, Zhang Z, Gibbs JE, Partridge SC, Chenevert TL, Rosen MA, Bolan PJ, Marques HS, Aliu S, Li W, et al. Test-retest repeatability and reproducibility of ADC measures by breast DWI: results from the ACRIN 6698 trial. J Magn Reson Imaging. 2019;49(6):1617–28.
Partridge SC, Zhang Z, Newitt DC, Gibbs JE, Chenevert TL, Rosen MA, Bolan PJ, Marques HS, Romanoff J, Cimino L, et al. Diffusion-weighted MRI findings predict pathologic response in neoadjuvant treatment of breast cancer: the ACRIN 6698 Multicenter Trial. Radiology. 2018;289(3):618–27.
Barentsz MW, Taviani V, Chang JM, Ikeda DM, Miyake KK, Banerjee S, van den Bosch MA, Hargreaves BA, Daniel BL. Assessment of tumor morphology on diffusion-weighted (DWI) breast MRI: diagnostic value of reduced field of view DWI. J Magn Reson Imaging. 2015;42(6):1656–65.
Hancu I, Lee SK, Hulsey K, Lenkinski R, Holland D, Sperl JI, Tan ET. Distortion correction in diffusion-weighted imaging of the breast: performance assessment of prospective, retrospective, and combined (prospective + retrospective) approaches. Magn Reson Med. 2017;78(1):247–53.
Teruel JR, Fjøsne HE, Østlie A, Holland D, Dale AM, Bathen TF, Goa PE. Magn Reson Med. 2015;74(4):1138-44. https://0-doi-org.brum.beds.ac.uk/10.1002/mrm.25489. Epub 2014 Oct 16. PMID: 25323982.
Bogner W, Pinker-Domenig K, Bickel H, Chmelik M, Weber M, Helbich TH, Trattnig S, Gruber S. Readout-segmented echo-planar imaging improves the diagnostic performance of diffusion-weighted MR breast examinations at 3.0 T. Radiology. 2012;263(1):64–76.
Uematsu T, Kasami M, Yuen S, Igarashi T, Nasu H. Comparison of 3- and 1.5-T dynamic breast MRI for visualization of spiculated masses previously identified using mammography. AJR Am J Roentgenol. 2012;198(6):W611–7.
Rahbar H, DeMartini WB, Lee AY, Partridge SC, Peacock S, Lehman CD. Accuracy of 3 T versus 1.5 T breast MRI for pre-operative assessment of extent of disease in newly diagnosed DCIS. Eur J Radiol. 2015;84(4):611–6.
Butler RS, Chen C, Vashi R, Hooley RJ, Philpotts LE. 3.0 tesla vs 1.5 tesla breast magnetic resonance imaging in newly diagnosed breast cancer patients. World J Radiol. 2013;5(8):285–94.
Nguyen VT, Rahbar H, Olson ML, Liu CL, Lehman CD, Partridge SC. Diffusion-weighted imaging: effects of intravascular contrast agents on apparent diffusion coefficient measures of breast malignancies at 3 tesla. J Magn Reson Imaging. 2015;42(3):788–800.
Dorrius MD, Dijkstra H, Oudkerk M, Sijens PE. Effect of b value and pre-admission of contrast on diagnostic accuracy of 1.5-T breast DWI: a systematic review and meta-analysis. Eur Radiol. 2014;24(11):2835–47.
The authors gratefully acknowledge the patients who participated in the study and Debosmita Biswas, M.S., for the assistance with illustrations.
Funding for the study was provided by the National Cancer Institute (NCI) through grants R01 CA151326 and R01 CA207290.
Ethics approval and consent to participate
Institutional IRB approval was obtained (FHCRC IRO#7206).
Consent for publication
The authors disclose research grants and/or non-financial research support from GE Healthcare (SCP, HR, DSH), Philips Healthcare (SCP, HR, DSH), Toshiba America Medical Systems (DSH), and Siemens Medical Solutions USA (DSH) outside the submitted work. The other authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Luo, J., Hippe, D.S., Rahbar, H. et al. Diffusion tensor imaging for characterizing tumor microstructure and improving diagnostic performance on breast MRI: a prospective observational study. Breast Cancer Res 21, 102 (2019). https://0-doi-org.brum.beds.ac.uk/10.1186/s13058-019-1183-3
- Diffusion tensor imaging (DTI)
- Dynamic contrast-enhanced (DCE) MRI
- Breast MRI
- Apparent diffusion coefficient (ADC)
- Fractional anisotropy (FA)
- Suspicious lesions
- False positives