Atlas-based probabilistic fibroglandular tissue segmentation in breast MRI

A publicly available deep learning model and dataset for segmentation of breast, fibroglandular tissue, and vessels in breast MRI

Breast density, or the amount of fibroglandular tissue (FGT) relative to the overall breast volume, increases the risk of developing breast cancer. Although previous studies have utilized deep learning to assess breast density, the limited public availability of data and quantitative tools hinders the development of better assessment tools. Our objective was to (1) create and share a large dataset of pixel-wise annotations according to well-defined criteria, and (2) develop, evaluate, and share an automated segmentation method for breast, FGT, and blood vessels using convolutional neural networks. We used the Duke Breast Cancer MRI dataset to randomly select 100 MRI studies and manually annotated the breast, FGT, and blood vessels for each study. Model performance was evaluated using the Dice similarity coefficient (DSC). The model achieved DSC values of 0.92 for breast, 0.86 for FGT, and 0.65 for blood vessels on the test set. The correlation between our model's predicted breast density and the manually generated masks was 0.95. The correlation between the predicted breast density and qualitative radiologist assessment was 0.75. Our automated models can accurately segment breast, FGT, and blood vessels using pre-contrast breast MRI data. The data and the models were made publicly available.

Automatic and fast segmentation of breast region-of-interest (ROI) and density in MRIs

Accurate segmentation of the breast region of interest (BROI) and breast density (BD) is a significant challenge during the analysis of breast MR images. Most of the existing methods for breast segmentation are semi-automatic and limited in their ability to achieve accurate results. This is because of difficulties in removing landmarks from noisy magnetic resonance images (MRI) due to similar intensity levels and the close connection to BROI. This study proposes an innovative, fully automatic and fast segmentation approach to identify and remove landmarks such as the heart and pectoral muscles. The BROI segmentation is carried out with a framework consisting of three major steps. Firstly, we use adaptive wiener filtering and k-means clustering to minimize the influence of noises, preserve edges and remove unwanted artefacts. The second step systematically excludes the heart area by utilizing active contour based level sets where initial contour points are determined by the maximum entropy thresholding and convolution method. Finally, a pectoral muscle is removed by using morphological operations and local adaptive thresholding on MR images. Prior to the elimination of the pectoral muscle, the MR image is sub divided into three sections: left, right, and central based on the geometrical information. Subsequently, a BD segmentation is achieved with 4 level fuzzy c-means (FCM) thresholding on the denoised BROI segmentation. The proposed method is validated using the 1350 breast images from 15 female subjects. The pixel-based quantitative analysis showed excellent segmentation results when compared with manually drawn BROI and BD. Furthermore, the presented results in terms of evaluation matrices: Acc, Sp, AUC, MR, P, Se and DSC demonstrate the high quality of segmentations using the proposed method. The average computational time for the segmentation of BROI and BD is 1 minute and 50 seconds.

Breast MRI segmentation for density estimation: Do different methods give the same results and how much do differences matter?

PURPOSE

To compare two methods of automatic breast segmentation with each other and with manual segmentation in a large subject cohort. To discuss the factors involved in selecting the most appropriate algorithm for automatic segmentation and, in particular, to investigate the appropriateness of overlap measures (e.g., Dice and Jaccard coefficients) as the primary determinant in algorithm selection.

METHODS

Two methods of breast segmentation were applied to the task of calculating MRI breast density in 200 subjects drawn from the Avon Longitudinal Study of Parents and Children, a large cohort study with an MRI component. A semiautomated, bias-corrected, fuzzy C-means (BC-FCM) method was combined with morphological operations to segment the overall breast volume from in-phase Dixon images. The method makes use of novel, problem-specific insights. The resulting segmentation mask was then applied to the corresponding Dixon water and fat images, which were combined to give Dixon MRI density values. Contemporaneously acquired T1 - and T2 -weighted image datasets were analyzed using a novel and fully automated algorithm involving image filtering, landmark identification, and explicit location of the pectoral muscle boundary. Within the region found, fat-water discrimination was performed using an Expectation Maximization-Markov Random Field technique, yielding a second independent estimate of MRI density.

RESULTS

Images are presented for two individual women, demonstrating how the difficulty of the problem is highly subject-specific. Dice and Jaccard coefficients comparing the semiautomated BC-FCM method, operating on Dixon source data, with expert manual segmentation are presented. The corresponding results for the method based on T1 - and T2 -weighted data are slightly lower in the individual cases shown, but scatter plots and interclass correlations for the cohort as a whole show that both methods do an excellent job in segmenting and classifying breast tissue.

CONCLUSIONS

Epidemiological results demonstrate that both methods of automated segmentation are suitable for the chosen application and that it is important to consider a range of factors when choosing a segmentation algorithm, rather than focus narrowly on a single metric such as the Dice coefficient.

MRI Background Parenchymal Enhancement Is Not Associated with Breast Cancer

Background

Previously, a strong positive association between background parenchymal enhancement (BPE) at magnetic resonance imaging (MRI) and breast cancer was reported in high-risk populations. We sought to determine, whether this was also true for non-high-risk patients.

Methods

540 consecutive patients underwent breast MRI for assessment of breast findings (BI-RADS 0–5, non-high-risk screening (no familial history of breast cancer, no known genetic mutation, no prior chest irradiation, or previous breast cancer diagnosis)) and subsequent histological work-up. For this IRB-approved study, BPE and fibroglandular tissue FGT were retrospectively assessed by two experienced radiologists according to the BI-RADS lexicon. Pearson correlation coefficients were calculated to explore associations between BPE, FGT, age and final diagnosis of breast cancer. Subsequently, multivariate logistic regression analysis, considering covariate colinearities, was performed, using final diagnosis as the target variable and BPE, FGT and age as covariates.

Results

Age showed a moderate negative correlation with FGT (r = -0.43, p<0.001) and a weak negative correlation with BPE (r = -0.28, p<0.001). FGT and BPE correlated moderately (r = 0.35, p<0.001). Final diagnosis of breast cancer displayed very weak negative correlations with FGT (r = -0.09, p = 0.046) and BPE (r = -0.156, p<0.001) and weak positive correlation with age (r = 0.353, p<0.001). On multivariate logistic regression analysis, the only independent covariate for prediction of breast cancer was age (OR 1.032, p<0.001).

Conclusions

Based on our data, neither BPE nor FGT independently correlate with breast cancer risk in non-high-risk patients at MRI. Our model retained only age as an independent risk factor for breast cancer in this setting.

Imaging Breast Density: Established and Emerging Modalities

Mammographic density has been proven as an independent risk factor for breast cancer. Women with dense breast tissue visible on a mammogram have a much higher cancer risk than women with little density. A great research effort has been devoted to incorporate breast density into risk prediction models to better estimate each individual’s cancer risk. In recent years, the passage of breast density notification legislation in many states in USA requires that every mammography report should provide information regarding the patient’s breast density. Accurate definition and measurement of breast density are thus important, which may allow all the potential clinical applications of breast density to be implemented. Because the two-dimensional mammography-based measurement is subject to tissue overlapping and thus not able to provide volumetric information, there is an urgent need to develop reliable quantitative measurements of breast density. Various new imaging technologies are being developed. Among these new modalities, volumetric mammographic density methods and three-dimensional magnetic resonance imaging are the most well studied. Besides, emerging modalities, including different x-ray–based, optical imaging, and ultrasound-based methods, have also been investigated. All these modalities may either overcome some fundamental problems related to mammographic density or provide additional density and/or compositional information. The present review article aimed to summarize the current established and emerging imaging techniques for the measurement of breast density and the evidence of the clinical use of these density methods from the literature.

Dose-response effects of aerobic exercise on estrogen among women at high risk for breast cancer: a randomized controlled trial

Medical and surgical interventions for elevated breast cancer risk (e.g., BRCA1/2 mutation, family history) focus on reducing estrogen exposure. Women at elevated risk may be interested in less aggressive approaches to risk reduction. For example, exercise might reduce estrogen, yet has fewer serious side effects and less negative impact than surgery or hormonal medications. Randomized controlled trial. Increased risk defined by risk prediction models or BRCA mutation status. Eligibility: Age 18-50, eumenorrheic, non-smokers, and body mass index (BMI) between 21 and 50 kg/m(2). 139 were randomized. Treadmill exercise: 150 or 300 min/week, five menstrual cycles. Control group maintained exercise <75 min/week.

PRIMARY OUTCOME

Area under curve (AUC) for urinary estrogen. Secondary measures: urinary progesterone, quantitative digitized breast dynamic contrast-enhanced magnetic resonance imaging background parenchymal enhancement. Mean age 34 years, mean BMI 26.8 kg/m(2). A linear dose-response relationship was observed such that every 100 min of exercise is associated with 3.6 % lower follicular phase estrogen AUC (linear trend test, p = 0.03). No changes in luteal phase estrogen or progesterone levels. There was also a dose-response effect noted: for every 100 min of exercise, there was a 9.7 % decrease in background parenchymal enhancement as measured by imaging (linear trend test, p = 0.009). Linear dose-response effect observed to reduce follicular phase estrogen exposure measured via urine and hormone sensitive breast tissue as measured by imaging. Future research should explore maintenance of effects and extent to which findings are repeatable in lower risk women. Given the high benefit to risk ratio, clinicians can inform young women at increased risk that exercise may blunt estrogen exposure while considering whether to try other preventive therapies.

Automatic segmentation of breast MR images through a Markov random field statistical model

An algorithm dedicated to automatic segmentation of breast magnetic resonance images is presented in this paper. Our approach is based on a pipeline that includes a denoising step and statistical segmentation. The noise removal preprocessing relies on an anisotropic diffusion scheme, whereas the statistical segmentation is conducted through a Markov random field model. The continuous updating of all parameters governing the diffusion process enables automatic denoising, and the partial volume effect is also addressed during the labeling step. To assess the relevance, the Jaccard similarity coefficient was computed. Experiments were conducted on synthetic data and breast magnetic resonance images extracted from a high-risk population. The relevance of the approach for the dataset is highlighted, and we demonstrate accuracy superior to that of traditional clustering algorithms. The results emphasize the benefits of both denoising guided by input data and the inclusion of spatial dependency through a Markov random field. For example, the Jaccard coefficient for the clinical data was increased by 114%, 109%, and 140% with respect to a K-means algorithm and, respectively, for the adipose, glandular and muscle and skin components. Moreover, the agreement between the manual segmentations provided by an experienced radiologist and the automatic segmentations performed with this algorithm was good, with Jaccard coefficients equal to 0.769, 0.756, and 0.694 for the above-mentioned classes.

Automated fibroglandular tissue segmentation and volumetric density estimation in breast MRI using an atlas-aided fuzzy C-means method

PURPOSE

Breast magnetic resonance imaging (MRI) plays an important role in the clinical management of breast cancer. Studies suggest that the relative amount of fibroglandular (i.e., dense) tissue in the breast as quantified in MR images can be predictive of the risk for developing breast cancer, especially for high-risk women. Automated segmentation of the fibroglandular tissue and volumetric density estimation in breast MRI could therefore be useful for breast cancer risk assessment.

METHODS

In this work the authors develop and validate a fully automated segmentation algorithm, namely, an atlas-aided fuzzy C-means (FCM-Atlas) method, to estimate the volumetric amount of fibroglandular tissue in breast MRI. The FCM-Atlas is a 2D segmentation method working on a slice-by-slice basis. FCM clustering is first applied to the intensity space of each 2D MR slice to produce an initial voxelwise likelihood map of fibroglandular tissue. Then a prior learned fibroglandular tissue likelihood atlas is incorporated to refine the initial FCM likelihood map to achieve enhanced segmentation, from which the absolute volume of the fibroglandular tissue (|FGT|) and the relative amount (i.e., percentage) of the |FGT| relative to the whole breast volume (FGT%) are computed. The authors' method is evaluated by a representative dataset of 60 3D bilateral breast MRI scans (120 breasts) that span the full breast density range of the American College of Radiology Breast Imaging Reporting and Data System. The automated segmentation is compared to manual segmentation obtained by two experienced breast imaging radiologists. Segmentation performance is assessed by linear regression, Pearson's correlation coefficients, Student's paired t-test, and Dice's similarity coefficients (DSC).

RESULTS

The inter-reader correlation is 0.97 for FGT% and 0.95 for |FGT|. When compared to the average of the two readers' manual segmentation, the proposed FCM-Atlas method achieves a correlation of r = 0.92 for FGT% and r = 0.93 for |FGT|, and the automated segmentation is not statistically significantly different (p = 0.46 for FGT% and p = 0.55 for |FGT|). The bilateral correlation between left breasts and right breasts for the FGT% is 0.94, 0.92, and 0.95 for reader 1, reader 2, and the FCM-Atlas, respectively; likewise, for the |FGT|, it is 0.92, 0.92, and 0.93, respectively. For the spatial segmentation agreement, the automated algorithm achieves a DSC of 0.69 ± 0.1 when compared to reader 1 and 0.61 ± 0.1 for reader 2, respectively, while the DSC between the two readers' manual segmentation is 0.67 ± 0.15. Additional robustness analysis shows that the segmentation performance of the authors' method is stable both with respect to selecting different cases and to varying the number of cases needed to construct the prior probability atlas. The authors' results also show that the proposed FCM-Atlas method outperforms the commonly used two-cluster FCM-alone method. The authors' method runs at ∼5 min for each 3D bilateral MR scan (56 slices) for computing the FGT% and |FGT|, compared to ∼55 min needed for manual segmentation for the same purpose.

CONCLUSIONS

The authors' method achieves robust segmentation and can serve as an efficient tool for processing large clinical datasets for quantifying the fibroglandular tissue content in breast MRI. It holds a great potential to support clinical applications in the future including breast cancer risk assessment.

Deformable registration for quantifying longitudinal tumor changes during neoadjuvant chemotherapy

PURPOSE

To evaluate DRAMMS, an attribute-based deformable registration algorithm, compared to other intensity-based algorithms, for longitudinal breast MRI registration, and to show its applicability in quantifying tumor changes over the course of neoadjuvant chemotherapy.

METHODS

Breast magnetic resonance images from 14 women undergoing neoadjuvant chemotherapy were analyzed. The accuracy of DRAMMS versus five intensity-based deformable registration methods was evaluated based on 2,380 landmarks independently annotated by two experts, for the entire image volume, different image subregions, and patient subgroups. The registration method with the smallest landmark error was used to quantify tumor changes, by calculating the Jacobian determinant maps of the registration deformation.

RESULTS

DRAMMS had the smallest landmark errors (6.05 ± 4.86 mm), followed by the intensity-based methods CC-FFD (8.07 ± 3.86 mm), NMI-FFD (8.21 ± 3.81 mm), SSD-FFD (9.46 ± 4.55 mm), Demons (10.76 ± 6.01 mm), and Diffeomorphic Demons (10.82 ± 6.11 mm). Results show that registration accuracy also depends on tumor versus normal tissue regions and different patient subgroups.

CONCLUSIONS

The DRAMMS deformable registration method, driven by attribute-matching and mutual-saliency, can register longitudinal breast magnetic resonance images with a higher accuracy than several intensity-matching methods included in this article. As such, it could be valuable for more accurately quantifying heterogeneous tumor changes as a marker of response to treatment.