View/Download PDF
Original Research Article
2019
:9;
43
doi:
10.25259/JCIS_70_2019
PDF

Mammographic Breast Density Assessed with Fully Automated Method and its Risk for Breast Cancer

Department of Medical Imaging Technology, Manipal College of Health Professions, Manipal, Karnataka, India.
School of Information Sciences, Manipal Institute of Technology, Manipal, Karnataka, India.
Department of Radiodiagnosis, Kasturba Medical College and Hospital, Manipal Academy of Higher Education, Manipal, Karnataka, India.
Corresponding author: Dr. Pendem Saikiran, Department of Medical Imaging Technology, Manipal College of Health Professions, Manipal Academy of Higher Education, Manipal - 576 104, Karnataka, India. kiran30188@gmail.com
Licence
This is an open-access article distributed under the terms of the Creative Commons Attribution-Non Commercial-Share Alike 4.0 License, which allows others to remix, tweak, and build upon the work non-commercially, as long as the author is credited and the new creations are licensed under the identical terms.
How to cite this article: Saikiran P, Ramzan R, Nandish S, Kamineni PD, Priyanka, John AM. Mammographic breast density assessed with fully automated method and its risk for breast cancer. J Clin Imaging Sci 2019;9:43.

Abstract

Objectives:

We evaluated the association between breast cancer and breast density (BD) measured using fully automated software. We also evaluated the performance of cancer risk models such as only clinical risk factors, density related measures, and both clinical risk factors and density-related measures for determining cancer risk.

Materials and Methods:

This is a retrospective case–control study. The data were collected from August 2015 to December 2018. Two hundred fifty women with breast cancer and 400 control subjects were included in this study. We evaluated the BD qualitatively using breast imaging-reporting and data system density and quantitatively using 3D slicer. We also collected clinical factors such as age, familial history of breast cancer, menopausal status, number of births, body mass index, and hormonal replacement therapy use. We calculated the odds ratio (OR) for BD to determine the risk of breast cancer. We performed receiver operating characteristic (ROC) curve to assess the performance of cancer risk models.

Results:

The OR for the percentage BD for second, third, and fourth quartiles was 1.632 (95% confidence intervals [CI]: 1.102–2.416), 2.756 (95% CI: 1.704–4.458), and 3.163 (95% CI: 1.356–5.61). The area under ROC curve for clinical risk factors only, mammographic density measures, combined mammographic, and clinical risk factors was 0.578 (95% CI: 0.45, 0.64), 0.684 (95% CI: 0.58, 0.75), and 0.724 (95% CI: 0.64, 0.80), respectively.

Conclusion:

Mammographic BD was found to be positively associated with breast cancer. The density related measures combined clinical risk factors, and density model had good discriminatory power in identifying the cancer risk.

Keywords

Breast density
Breast imaging-reporting and data system
3D slicer
Body mass index
Breast cancer

INTRODUCTION

Breast screening is widely recommended as women are at risk of breast cancer. There is number of risk factors associated with increased risk of breast cancer such as age, family history, menopausal status, and hormone replacement therapy (HRT). There are many cancer risk models used in practice such as the Gail, Claus, BRCAPRO, BOADICEA, and Rosner-Colditz models. They have not assessed the mammographic breast density (BD). The fibroglandular tissue in the breast on mammogram represents the mammographic BD and found to be independent risk factor of breast cancer. The estimated odds ratio (OR) ranges from 2.9 to 6.0. Females having extremely dense breasts are at higher risk of cancer than those with fatty breasts.[1-9]

There are qualitative and quantitative methods for assessing the mammographic BD. The qualitative method includes American College of Radiology’s breast imaging-reporting and data system (BI-RADS), which has four qualitative BD categories: Almost entirely fatty, scattered areas of fibroglandular density, heterogeneously dense, and extremely dense. BI-RADS is commonly used in clinical practice for reporting mammographic BD. However, there was too much variability in assigning grades by the radiologists; hence, there arise need for quantitative methods. Several quantitative methods have been proposed by different studies. The quantitative methods for measuring BD include interactive thresholding, semiautomatic, fully automatic, area based, and volumetric based methods. To the best of our knowledge, studies have evaluated the use of volumetric versus area density measures in breast cancer risk assessment who have primarily focused on determining which individual measure is a better predictor of risk and reported mixed.[10-13] Hence, aim of the study is to determine the association between breast cancer and BD in a case–control setting after adjusting for clinical risk factors and measured using fully automated 3D slicer software. We also evaluated the performance of cancer risk models such as only clinical risk factors, density related measures, and both clinical risk factors and density related measures.

MATERIALS AND METHODS

This is a retrospective age-matched case–control study. The study was approved by the Institutional Ethical Committee (IEC: 267/2018). A total of 250 women who underwent screening mammography at our hospital between August 2015 and December 2018 who were diagnosed with breast cancer were included as cases. A total of 400 women without history of breast cancer who underwent screening at our hospital between August 2015 and December 2018 were included as controls. The women who underwent two or more screening mammograms negative for cancer were included as controls. For all the women X-ray, mammography and ultrasonography (USG) breast was performed. All the mammograms were obtained using GE Senographe DMR plus Mammography (GE Healthcare Systems, Milwaukee, Wisconsin, United States). The USG breast was performed using Philips HD 11XE ultrasound (Philips Healthcare systems, Amsterdam, Netherlands).

BD assessment

The quantitative assessment of BD was performed using 3D slicer version 3.6.3 (Harvard University, Boston, UK). To avoid the effect of presence of cancer on the density measurement, density was calculated on the contralateral breast which is noncancerous. We calculated BD on craniocaudal view of the noncancerous breast for the patients and left craniocaudal view for the matched control subjects. The software automatically identifies the pixel with higher intensities [Figure 1a and b]. The density was measured in percentage. The percentage area density was calculated as the dense area divided by the total breast area.

Figure 1:: (a) Left craniocaudal mammogram of a 41-year-old woman with a negative screening exam (b) example of density segmentation using 3D slicer software tool showing automatic identification of mammographic density.

The qualitative assessment of BD was performed using BI- RADS (5th edition, American College of Radiology) Visual method. Two radiologists who have experience of more than 10 years have performed the qualitative assessment.

Clinical risk factors

Data regarding clinical risk factors were obtained from the medical records. Risk factors which are included in the study were: Age (years), number of births, family history of breast cancer (yes/no), HRT use at the time of screen (yes/no), body mass index (BMI), and menopausal status (premenopausal/ postmenopausal).

Statistical analysis

The data were analyzed using Statistical Package for the Social Sciences (SPSS Inc., Chicago, IL, USA) version 18.0. Descriptive analysis was performed for clinical risk factors and percentage mammographic density and expressed as mean and standard deviations or frequencies and proportions as appropriate. Logistic regression was performed to examine the effects of BI-RADS density on invasive breast cancer risk. Associations were expressed in OR and 95% confidence intervals (CI). The receiver operating characteristic (ROC) curve was used to evaluate the predictive performance of three risk models, such as risk factor model, mammographic density breast model, and combined risk factor, and mammographic density model. Pearson’s correlation was done to determine the correlation between qualitative and quantitative analysis. The interobserver variability about BD according to BI-RADS category was assessed using kappa value.

RESULTS

A total of 250 women with invasive cancer and 400 control subjects were included in the study. The results of the descriptive analysis of clinical risk factors and mammographic densities are mentioned in Table 1. The cases were found to have a larger proportion of patients who have familial history of breast cancer, who use HRT, postmenopausal, had BMI >30 kg/m2. The cases were found to have higher number of heterogeneously or extremely dense breasts. There was good correlation between the quantitative and qualitative measurement; Pearson’s correlation coefficient (r) was 0.79. The kappa value was 0.68, which shows good agreement in estimating BI-RADS category between two readers.

Table 1:: The characteristics of cases and control subjects.
Patient characteristics Cases n=250 Controls n=400
Age (years) mean (SD) 50.15 (10.5) 50.15 (10.5)
Number of births n (%)
0 4 (1.6) 12 (3)
1 72 (28.8) 136 (34)
2 116 (46.4) 162 (40.5)
3 46 (18.4) 66 (16.5)
4 12 (4.8) 24 (6)
Familial history of breast cancer n (%)
Yes 83 (33.2) 94 (23.5)
HRT use n (%)
Yes 29 (11.6) 27 (6.75)
BMI n (%)
<25 65 (26) 148 (37)
25–29 89 (35.6) 118 (29.5)
30+ 96 (38.4) 134 (33.5)
Menopausal status
Pre 50 (20) 287 (21.7)
Post 200 (80) 113 (78)
Percentage area density 23.04 (11.29) 17.86 (12.09)
BIRADS density (%)
BIRADS A–Predominantly fatty 32 (12.8) 109 (27.25)
BIRADS B–Scattered fibroglandular 100 (40) 160 (40)
BIRADS C–Heterogeneously dense 92 (36.8) 107 (26.75)
BIRADS D–Extremely dense 26 (10.4) 24 (6)
3D slicer area density n (%)
<15.54 28 (11.2) 93 (23.25)
15.54–30.19 85 (34) 152 (38)
30.19–45.28 88 (35.2) 101 (25.25)
>45.28 49 (19.6) 54 (13.5)

HRT: Hormone replacement therapy, BMI: Body mass index, BI-RADS: Breast imaging-reporting and data system

Association of BD and risk of cancer measured with 3D slicer and BI-RADS

The association of BD with cancer measured with quantitative measurement is given in Table 2. The association of BD with cancer was more in heterogeneous (OR 2.75, 95% CI 1.704–4.458) and extreme dense breasts (OR 3.16, 95% CI 1.356–5.610) compared to fatty and scattered heterogeneous (OR 1.63, 95% CI 1.102–2.416). The association of BD with cancer measured with BI-RADS category is given in Table 3. Clinical assessment with BI-RADS allowed slightly better discrimination of patients from controls.

Table 2:: The breast cancer risk associated with quantitative assessment of mammographic density.
Percentage density quartiles Cases (n) Controls (n) Fully adjusted odds ratio P-value Confidence interval
<15.54% (reference) 28 93 - - -
15.54–30.19 85 152 1.63 <0.015 1.102–2.416
30.19–45.28 88 101 2.75 <0.001 1.704–4.458
>45.28 49 54 3.16 <0.013 1.356–5.610
Table 3:: Breast cancer risk associated with BI-RADS and 3D slicer categories.
Density category BI-RADS (OR, CI) 3D slicer (OR, CI)
Entirely fatty Reference Reference
Scattered fibroglandular 1.76 (1.309–2.735) 1.63 (1.102–2.416)
Heterogeneously dense 3.00 (1.916–4.352) 2.75 (1.704–4.458)
Extremely dense 3.67 (1.586–5.924) 3.16 (1.356–5.610)
P-value <0.001 <0.05

BI-RADS: Breast imaging-reporting and data system, OR: Odds ratio, CI: Confidence interval

Discriminatory ability of the models

Figure 2 shows the ROC curves for the risk factor model, density model, combined risk factor, and density model. The ROC curve for the first model of association between breast cancer risk and clinical risk factors had area under curve of 0.578 (95% CI: 0.45, 0.64). The ROC curve for the second model of association between breast cancer risk and BD had area under curve of 0.684 (95% CI: 0.58, 0.75). The ROC curve for third model of association between breast cancer risk and combined clinical factors and BD had area under curve of 0.724 (95% CI: 0.64, 0.80). All three models have good significance in differentiating cases and controls and suggesting the risk of cancer more in patients compared to controls. However, the density related model performed significantly better than the first model (P < 0.001). The third model performed better than the first and second models (P < 0.001).

Figure 2:: Comparison of receiver operating characteristic curves for the three models.

DISCUSSION

In our age-matched case–control study, we assessed area BD using open source software 3D slicer; we found that scattered fibro glandular, heterogeneously dense, extremely dense breasts were found to be at risk of cancer. The finding of our study is concordant with the literature that increased levels of BD are associated with increased risk of breast cancer.[1-9] It is established that mammographic density reduces the diagnosis due to masking of the underlying pathology.[14,15] In our study, we noticed in 53 cases the breast lesions were found to be missing on X-ray mammography but detected in USG. Out of 53 cases, 45 cases belong to BI-RADS B, C, and D categories, which is evident of masking effect of BD in detection of breast lesions.

We found that mammographic density model, combined density, and clinical risk factor model improve discriminatory accuracy especially in women with dense breasts. The quantitative method of BD assessment such as area and volume based may improve the discrimination of women with dense breasts who are at risk of breast cancer, as qualitative assessment showed poor agreement among the readers.[16,17]

Our study noticed a good correlation (r = 0.79) between qualitative and quantitative assessment which is supported by literature. There are several studies which assessed the correlation between qualitative and quantitative and found good agreement between the two methods.[18-20]

In our study, the OR (1.63 CI 1.102–2.416, 2.75 CI 1.704– 4.458, 3.16 CI 1.356–5.601) was slightly higher compared to the studies who assessed BD using fully automated methods such as Libra (OR 1.66 CI 1.23–2.24, 2.17 CI 1.37–3.44, 2.71 CI 1.5–4.91) and Image J (OR 0.64 CI 0.26–1.59, 0.84 CI 0.53–1.36, 2.14 CI 1.44–3.16).[21,22]

Our study found an inverse relationship between patients age and mammographic density consistent with results of other studies.[23-26] We found higher density breasts in the menopausal age of 40–50 years. Checka et al.[27] found dense breasts in extreme age of 70–79 and >80. However, our study did not find dense breasts in an extreme age.

First, the strength of our study was the combined use of clinical risk factors and BD in assessing the risk of breast cancer in patients. The second strength of our study was the use of screening mammograms and images of noncancerous breast to avoid the influence of the presence of cancer on the density assessments.

There are few limitations of our study. First, our study did not include all the clinical risk factors of breast cancer. Second, sample size was less, a larger sample of may have revealed further relationships, and further studies should focus on recruiting larger sample sizes.

CONCLUSION

We found that mammographic density measured by a fully automated method using 3D slicer was a significant factor contributing to risk of cancer. Besides, its significance in breast cancer quantitative mammographic density assessment is not currently considered in screening programs. Our study findings recommend inclusion of quantitative BD measurements done using fully automated software in diagnostic radiology to improve the discriminatory accuracy of risk of cancer.

Declaration of patient consent

Patient’s consent not obtained as patients identity is not disclosed or compromised.

Financial support and sponsorship

Nil.

Conflicts of interest

There are no conflicts of interest.

References

  1. , , , , , , et al. Mammographic density and the risk and detection of breast cancer. N Engl J Med. 2007;356:227-36
    [CrossRef] [PubMed] [Google Scholar]
  2. , , , , , , et al. Pooled analysis of prospective cohort studies on height, weight, and breast cancer risk. Am J Epidemiol. 2000;152:514-27
    [CrossRef] [PubMed] [Google Scholar]
  3. , , , , , , et al. Mammographic breast density as a general marker of breast cancer risk. Cancer Epidemiol Biomarkers Prev. 2007;16:43-9
    [CrossRef] [Google Scholar]
  4. , . Breast cancer risk-assessment models. Breast Cancer Res. 2007;9:213
    [CrossRef] [PubMed] [Google Scholar]
  5. , , , , , , et al. Utility of relative and absolute measures of mammographic density vs clinical risk factors in evaluating breast cancer risk at time of screening mammography. Br J Radiol. 2016;89:1059
    [CrossRef] [PubMed] [Google Scholar]
  6. , , , , , , et al. Mammographic density phenotypes and risk of breast cancer: A meta-analysis. J Natl Cancer Inst. 2014;106:dju078
    [CrossRef] [PubMed] [Google Scholar]
  7. . Dense breast tissue as an important risk factor for breast cancer and implications for early detection. Breast Care (Basel). 2009;4:89-92
    [CrossRef] [PubMed] [Google Scholar]
  8. , . Mammographic density, hormone therapy, and risk of breast cancer. Cancer Epidemiol Biomarkers Prev. 2006;15:1750
    [CrossRef] [PubMed] [Google Scholar]
  9. , , , , , , et al. Mammographic features and breast cancer risk: Effects with time, age, and menopause status. J Natl Cancer Inst. 1995;87:1622-9
    [CrossRef] [PubMed] [Google Scholar]
  10. American College of Radiology Breast Imaging Reporting and Data System Atlas (BI-RADS® Atlas) In: Vol. Reston: American College of Radiology. 5
  11. , , . Imaging breast density: Established and emerging modalities. Transl Oncol. 2015;8:435-45
    [CrossRef] [PubMed] [Google Scholar]
  12. , , , , , , et al. Automated assessment of the composition of breast tissue revealed on tissue-thickness-corrected mammography. AJR Am J Roentgenol. 2003;180:257-62
    [CrossRef] [PubMed] [Google Scholar]
  13. , , , , , , et al. Automated volumetric breast density estimation: A comparison with visual assessment. Clin Radiol. 2013;68:690-5
    [CrossRef] [PubMed] [Google Scholar]
  14. , . Breast density and parenchymal patterns as markers of breast cancer risk: A meta-analysis. Cancer Epidemiol Biomarkers Prev. 2006;15:1159-69
    [CrossRef] [PubMed] [Google Scholar]
  15. , , , . Mammographic breast density and risk of breast cancer: Masking bias or causality? Eur J Epidemiol. 1998;14:315-20
    [Google Scholar]
  16. , , , , , , et al. Combining quantitative and qualitative breast density measures to assess breast cancer risk. Breast Cancer Res. 2017;19:97
    [CrossRef] [PubMed] [Google Scholar]
  17. , , , , , , et al. Variation in mammographic breast density assessments among radiologists in clinical practice: A multicenter observational study. Ann Intern Med. 2016;165:457-64
    [CrossRef] [PubMed] [Google Scholar]
  18. , , , , . Comparison of visual assessment of breast density in BI-RADS 4th and 5th editions with automated volumetric measurement. AJR Am J Roentgenol. 2017;209:703-8
    [CrossRef] [PubMed] [Google Scholar]
  19. , , , , . Qualitative versus quantitative mammographic breast density assessment: Applications for the US and abroad. Diagnostics (Basel). 2017;7:E30
    [CrossRef] [PubMed] [Google Scholar]
  20. , , , , , , et al. Breast cancer risk and mammographic density assessed with semiautomated and fully automated methods and BI-RADS. Radiology. 2017;282:348-55
    [CrossRef] [PubMed] [Google Scholar]
  21. , , , , , , et al. Comparison of fully and semi-automated area-based methods for measuring mammographic density and predicting breast cancer risk. Br J Cancer. 2014;110:1908-16
    [CrossRef] [PubMed] [Google Scholar]
  22. , , , , . Preliminary evaluation of the publicly available laboratory for breast radiodensity assessment (LIBRA) software tool: Comparison of fully automated area and volumetric density measures in a case-control study with digital mammography. Breast Cancer Res. 2015;17:117
    [CrossRef] [PubMed] [Google Scholar]
  23. , , , , , , et al. Age-specific trends in mammographic density: The Minnesota breast cancer family study. Am J Epidemiol. 2008;167:1027-36
    [CrossRef] [PubMed] [Google Scholar]
  24. , , , , , , et al. Mammographic density and ageing: A collaborative pooled analysis of cross-sectional data from 22 countries worldwide. PLoS Med. 2017;14:e1002335
    [CrossRef] [PubMed] [Google Scholar]
  25. , , , , , , et al. Association between mammographic density and age-related lobular involution of the breast. J Clin Oncol. 2010;28:2207-12
    [CrossRef] [PubMed] [Google Scholar]
  26. , , , , , , et al. Characteristics of mammographic breast density and associated factors for Chinese women: Results from an automated measurement. J Oncol. 2019;2019:1-9
    [CrossRef] [PubMed] [Google Scholar]
  27. , , , , . The relationship of mammographic density and age: Implications for breast cancer screening. AJR Am J Roentgenol. 2012;198:W292-5
    [CrossRef] [PubMed] [Google Scholar]
Show Sections