Interpretable deep learning regression for breast density estimation on MRI

Breast density, which is the ratio between fibroglandular tissue (FGT) and total breast volume, can be assessed qualitatively by radiologists and quantitatively by computer algorithms. These algorithms often rely on segmentation of breast and FGT volume. In this study, we propose a method to directly assess breast density on MRI, and provide interpretations of these assessments. We assessed breast density in 506 patients with breast cancer using a regression convolutional neural network (CNN). The input for the CNN were slices of breast MRI of 128 × 128 voxels, and the output was a continuous density value between 0 (fatty breast) and 1 (dense breast). We used 350 patients to train the CNN, 75 for validation, and 81 for independent testing. We investigated why the CNN came to its predicted density using Deep SHapley Additive exPlanations (SHAP). The density predicted by the CNN on the testing set was significantly correlated with the ground truth densities (N = 81 patients, Spearman's ρ = 0:86, P < 0:001). When inspecting what the CNN based its predictions on, we found that voxels in FGT commonly had positive SHAP-values, voxels in fatty tissue commonly had negative SHAP-values, and voxels in non-breast tissue commonly had SHAP-values near zero. This means that the prediction of density is based on the structures we expect it to be based on, namely FGT and fatty tissue. To conclude, we presented an interpretable deep learning regression method for breast density estimation on MRI with promising results.

[1]  Berkman Sahiner,et al.  Correlation between mammographic density and volumetric fibroglandular tissue estimated on breast MR images. , 2004, Medical physics.

[2]  Florentin Wörgötter,et al.  A deep learning framework for efficient analysis of breast volume and fibroglandular tissue using MR data with strong artifacts , 2019, International Journal of Computer Assisted Radiology and Surgery.

[3]  Shandong Wu,et al.  Automated fibroglandular tissue segmentation and volumetric density estimation in breast MRI using an atlas-aided fuzzy C-means method. , 2013, Medical physics.

[4]  Scott Lundberg,et al.  A Unified Approach to Interpreting Model Predictions , 2017, NIPS.

[5]  Jimmy Ba,et al.  Adam: A Method for Stochastic Optimization , 2014, ICLR.

[6]  J. Wolfe Breast patterns as an index of risk for developing breast cancer. , 1976, AJR. American journal of roentgenology.

[7]  V. McCormack,et al.  Breast Density and Parenchymal Patterns as Markers of Breast Cancer Risk: A Meta-analysis , 2006, Cancer Epidemiology Biomarkers & Prevention.

[8]  M. A. van den Bosch,et al.  Prediction Model For Extensive Ductal Carcinoma In Situ Around Early-Stage Invasive Breast Cancer , 2016, Investigative radiology.

[9]  Nico Karssemeijer,et al.  Using deep learning to segment breast and fibroglandular tissue in MRI volumes , 2017, Medical physics.

[10]  E. Morris,et al.  Contralateral parenchymal enhancement on dynamic contrast-enhanced MRI reproduces as a biomarker of survival in ER-positive/HER2-negative breast cancer patients , 2018, European Radiology.

[11]  Avanti Shrikumar,et al.  Learning Important Features Through Propagating Activation Differences , 2017, ICML.

[12]  Max A. Viergever,et al.  Automatic Slice Identification in 3D Medical Images with a ConvNet Regressor , 2016, LABELS/DLMIA@MICCAI.

[13]  M. Viergever,et al.  Complementary Value of Contralateral Parenchymal Enhancement on DCE-MRI to Prognostic Models and Molecular Assays in High-risk ER+/HER2− Breast Cancer , 2017, Clinical Cancer Research.

[14]  N. Boyd,et al.  Mammographic density and the risk and detection of breast cancer. , 2007, The New England journal of medicine.

[15]  L. Liberman,et al.  Breast imaging reporting and data system (BI-RADS). , 2002, Radiologic clinics of North America.

[16]  Samuel J. Magny,et al.  Breast Imaging Reporting and Data System , 2020, Definitions.

[17]  Brian B. Avants,et al.  N4ITK: Improved N3 Bias Correction , 2010, IEEE Transactions on Medical Imaging.

[18]  M. Viergever,et al.  Are contralateral parenchymal enhancement on dynamic contrast-enhanced MRI and genomic ER-pathway activity in ER-positive/HER2-negative breast cancer related? , 2019, European journal of radiology.

[19]  Min-Ying Su,et al.  Development of a quantitative method for analysis of breast density based on three-dimensional breast MRI. , 2008, Medical physics.

[20]  Ivan Dmitriev,et al.  Association between Parenchymal Enhancement of the Contralateral Breast in Dynamic Contrast-enhanced MR Imaging and Outcome of Patients with Unilateral Invasive Breast Cancer. , 2015, Radiology.