Breast density and breast cancer-specific survival by detection mode

Breast density is known to affect breast cancer risk and screening sensitivity, but it may also be associated with breast cancer survival. The interpretation of results from previous studies on breast density and survival is complicated by the association between detection mode and survival. Here, we studied the effect of breast density on breast cancer-specific survival for different detection modes (screen-detected, interval ≤ 24 or > 24 months, non-participant). Data from the Nijmegen (Dutch) breast cancer screening programme were used. Women diagnosed with invasive breast cancer between 1975 and 2011 were included. Breast density was assessed visually, based on a dichotomized Wolfe scale: ‘fatty breasts’ (≤25%) and ‘dense breasts’ (> 25%). Cox proportional hazard regression was used to obtain hazard ratios (HR). We identified 2742 eligible women, with a breast pattern available for 2233 women. A diagnosis of interval cancer (HR 2.06, 95% CI 1.62–2.61) led to a significantly increased risk of breast cancer death compared with screen-detected cancer. No significant cause-specific survival difference between women with dense and fatty breasts was observed (HR 0.94, 95% CI 0.77–1.15). The hazard was only higher for women with dense breasts among interval cancers ≤24 m (HR 1.07, 95% CI 0.74–1.56). The hazard appeared to be lower for women with dense breasts than for women with fatty breasts among screen-detected (HR 0.77, 95% CI 0.53–1.11) and interval cancers > 24 m (HR 0.80, 95% CI 0.53–1.20). None of the effects were statistically significant. Detection mode is strongly associated with breast cancer death. No clear association is apparent between breast density and breast cancer death, regardless of detection mode.


Background
Fibroglandular breast tissue is radiodense and appears as white areas on a mammogram. High breast density, i.e., a high proportion of dense tissue, has been associated with an increased breast cancer risk. The relative risk is estimated at 4 to 6 when women with more than 75% dense tissue are compared with women with less than 25% dense tissue [1]. In addition, dense breast tissue masks tumours on mammograms: tumour detection is decreased when breast density is high [2]. This results in a lower programme sensitivity in breast cancer screening in women with dense breasts.
It is unclear to what degree breast density is associated with breast cancer survival. Two different hypotheses may be suggested. High breast density delays cancer detection due to the masking effect [2]. Tumours occurring in dense breasts may thus on average have progressed to a higher stage, with worse prognostic characteristics, at diagnosis. This would affect breast cancer survival. Some studies have indeed observed an association between high breast density and indicators of tumour progression, although results are inconclusive [3][4][5][6][7]. Examples include larger tumour size and lymph node involvement.
The second hypothesis is based on an association with tumour aggressiveness. The etiological pathway through which breast density increases breast cancer risk is unknown. Breast tumours often originate from epithelial tissue. A higher number of epithelial cells and a potentially increased proliferation in dense tissue may thus explain the association, although an important role for the stromal tissue has also been suggested [8,9]. Whether the natural history of these tumours differs from tumours arising in fatty breasts has, however, not been determined. Some breast cancer risk factors (e.g., bodymass index [BMI] and several reproductive factors) are known to affect the risk of tumour subtypes differently, which may translate into a different disease prognosis and breast cancer survival [10][11][12]. The evidence on breast density and its possible involvement in pathways that lead to more aggressive tumours is still ambiguous [3][4][5][6][7][13][14][15][16][17][18]. A higher tumour aggressiveness, as indicated by tumour markers including hormone receptor status, also results in higher stage at diagnosis and poorer breast cancer survival.
The difficulty with survival analyses is that they may be affected by lead time bias when the study is conducted in a screening setting. Lead time is the time by which a breast cancer diagnosis is advanced due to early detection [19]. If survival time is measured as the time since diagnosis, the advanced diagnosis would increase survival time, even without an actual effect on breast cancer survival.
Since the distribution of screen-detected and interval tumours is expected to be different in women with high breast density compared with women with low breast density, this may affect survival estimates as well. To gain more insight into these different effects, the association between density and breast cancer survival could be studied separately for each detection mode, in addition to studying overall breast cancer survival.
Previous studies have addressed the effect of breast density on breast cancer survival or breast cancer mortality, but the results have been inconclusive [20][21][22][23][24][25][26][27][28]. The studies by Chiu et al. and Olsson et al. appeared to suggest an increased mortality in dense breasts, whereas other studies found no significant difference or even observed a potential survival benefit. Detection mode (i.e., screen-detected cancer, interval cancer, or cancer in non-participants) was not taken into account in most of these studies, apart from Olsson [21] and Eriksson [24], making it difficult to compare results and explain the observed differences.
In this study, we addressed the effect of breast density pattern on breast cancer-specific survival. Breast cancerspecific survival was studied for different modes of detection, in order to remove or reduce the potential role of lead and length time bias in our associations.

Population
The study was based on data from the Nijmegen (the Netherlands) breast cancer screening population. A pilot screening programme started here in 1975, with women 35 years or older being invited [29]. In 1989, a centralized nationwide programme was introduced in the Netherlands. This programme was aimed at women 50 to 69 years old. The screening age in Nijmegen was adapted accordingly. Since 1998, Dutch women aged 70 to 74 years also receive a screening invitation. The current programme invites women 50 to 74 years to mammographic screening every 2 years. Every mammogram is read independently by two certified screening radiologists. Women included in our study received at least one screening invitation before their diagnosis. All women consented to the use of their anonymous data for scientific research.

Breast cancer patients
Women diagnosed with invasive breast cancer between 1975 and 2011, both via screening and outside the screening programme, were included in our study. Women were excluded if it was unclear if the tumour truly had an invasive histology (n = 37). The diagnosis date was defined as the date at which clinical mammography was performed (n = 2655). When this date was unavailable, the date was set at 1 month prior to surgery (n = 77) or the date was obtained from the pathology report (n = 10). Women were excluded from the study when these dates were missing (n = 11). Information on vital status and cause of death were obtained via linkage with municipality records.
Women were followed from their diagnosis date until date of death, date of moving out of the area, or the end of study follow-up (December 31st, 2011). Follow-up was not always complete: it was unclear if some women were still alive, but no report of relocation or death was available (n = 652). The end date was then set at the year after the last screening invitation for a sensitivity analysis. We studied the effect of this censoring on the estimates.
Mode of detection was categorised into screendetected, interval, and non-participant cancer. Screendetected cancer was defined as all diagnoses after a positive screening exam. Interval cancer was further divided into tumours occurring within 24 months (≤24 m) and more than 24 months (> 24 m) after a negative screening exam. The category 'non-participant cancer' consists of cancer diagnoses in women who never participated in breast cancer screening.

Breast density
Screening is based on two-view mammography at first examination: mediolateral oblique (MLO) and craniocaudal (CC). In the time period that we studied, subsequent exams consisted of one view (MLO), with the CC view being used on indication. Digital mammography was introduced in the Nijmegen screening programme in 2007. Density was, however, only assessed on filmscreen mammograms.
Breast density was assessed visually by a trained research assistant. These assessments were based on all available mammography views. A quantified version of the Wolfe breast pattern scale was used, which consisted of the following categories: (1) < 5%, (2) 5-25%, (3) 26-75%, and (4) > 75% breast density. This scale was dichotomized for all analyses into 'fatty breasts' (≤25%) and 'dense breasts' (> 25%). A previous study has shown that there is a strong correlation between the Wolfe and the BI-RADS scale [30]. Breast density measurements were based on the last screening round before diagnosis (n = 1515). If this was not available, breast density was based on the diagnostic mammogram (n = 612) or the round preceding the last screening round (n = 106). Women were excluded from the analyses on breast density when no breast pattern was available in our dataset (n = 509).

Statistical analyses
Cox proportional hazard regression was used to assess the effect of breast density on breast cancer-specific survival. We thus report cause-specific survival estimates, rather than e..g. overall survival or disease-free survival. Breast density was included as a determinant in the analyses. This resulted in hazard ratios (HR) and corresponding 95% confidence intervals (CI). The time scale was from time of diagnosis to date of breast cancer death. Analyses were stratified by mode of detection. All analyses were corrected for age at diagnosis.
Kaplan-Meier curves were examined and timedependent covariates were included in the regression model to test the proportional hazards assumptions. All tests were two-sided, with a p-value < 0.05 being considered statistically significant. The analyses were performed with SPSS 22.0.

Results
We identified 2742 eligible women with an invasive tumour, with a breast pattern available for 2233 women. The age at diagnosis ranged from 38 to 97 years. Most women had a fatty breast pattern (n = 1456, 65.2%). There were 536 breast cancer deaths in this population. Breast density was known for 479 breast cancer deaths, of whom 173 (36.1%) had dense breasts. There were 149 breast cancer deaths among the 1257 screen-detected cases, 129 among the 557 interval cases ≤24 m, 154 among the 604 interval cases > 24 m, and 104 among the 324 non-participant cases.
Breast cancer-specific survival by detection mode The age-adjusted effect of detection mode on breast cancer-specific survival is presented in Table 1 Breast cancer-specific survival by breast density As shown in Table 2, there was no significant difference in breast cancer-specific survival between dense and fatty breasts when all detection modes were combined (age-adjusted HR 0.94, 95% CI 0.77-1.15). The hazard was only slightly higher for women with dense breasts among interval tumours ≤24 m (age-adjusted HR 1.07, 95% CI 0.74-1.56). The hazard appeared to be lower in women with dense breasts than fatty breasts among screen-detected (age-adjusted HR 0.77, 95% CI 0.53-1.11) and interval tumours > 24 m (age-adjusted HR 0.80, 95% CI 0.53-1.20). When screen-detected and interval tumours ≤24 m were combined, there appeared to be no difference in hazard (age-adjusted HR 1.01, 95% CI 0.78-1.31). None of the effects were, however, statistically significant. Adjustment of the overall effect for age and mode of detection resulted in a HR of 0.84 (95% CI 0.68-1.03), possibly suggesting an inverse association between high breast density and breast cancer-specific survival.
Additional adjustment for year of diagnosis did not change the estimates (data not shown). Furthermore, sensitivity analyses showed that censoring the follow-up after last invitation in women with apparent loss to follow-up had no effect on the estimates (data not shown).

Discussion
There was no clear overall association between breast density and breast cancer-specific survival. A high breast density might lead to missed detection and thus a later stage at diagnosis, when tumours are harder to treat (i.e. a true worsening of survival). Since screen-detected tumours are relatively less frequent in women with dense breasts due to missed detection, however, survival estimates are less affected by lead and length time bias in women with high breast density. These biases make survival associated with low breast density appear better, thus further increasing the difference in breast cancerspecific survival. By adjusting for mode of detection, we attempted to remove lead and length time bias from our association. This adjustment, however, also removed that part of the true association between breast density and breast cancer-specific survival due to missed detection's influence on later stage at diagnosis and subsequent worse survival. The possible inverse association between breast density and breast cancer-specific survival after adjustment would thus have to be attributed to other factors.
Breast cancer-specific survival in interval cancer is always expected to be worse than in screen-detected cancer [31]. The difference in breast cancer-specific survival (HR 2.06, 95% CI 1.62-2.61) may partly be explained by lead or length time bias, resulting in an artificial increase in survival time due to the advanced diagnosis in screendetected cancer or the detection of slow-growing tumours. The remaining breast cancer-specific survival difference would be explained by a better prognosis as a result of early detection. Detection mode was previously shown to give important prognostic information, even after adjustment for well-known prognostic characteristics [32].
The distribution of interval and screen-detected tumours is expected to differ between women with high and low breast density. As a result of the masking effect of dense tissue and subsequent missed detection, the risk of interval tumours is higher [20,[33][34][35][36]. Tumours are thus detected at a later stage, when they are harder to treat. Newer mammography systems are expected to decrease the masking effect. Although mammographic technology has greatly evolved between 1975 and 2011, adjustment for year of diagnosis did not appear to alter the estimates (results not shown). This suggests that the association, or rather lack of association, is consistent over time.
After adjusting the overall estimates for detection mode, women with dense breasts appeared to have a survival advantage compared to women with fatty breasts. This potential inverse association between breast density would have to be explained by other factors than missed detection, e.g. an inverse association with tumour aggressiveness, but we have to be cautious with the interpretation of these results. There are several studies on the association between breast density and prognostic    tumour characteristics [3][4][5][6][7][13][14][15][16][17][18]. Results on associations with tumour markers of aggressiveness, e.g. hormone receptor status, appear to be inconclusive. Some studies observed associations with tumour markers of progression, e.g. large tumour size and lymph node involvement, which may be a consequence of late detection. Most studies were not able to differentiate between screen-detected and interval tumours, thus limiting conclusions that may be drawn on the contribution of late detection to these associations. Several previous studies have addressed the association between breast density and breast cancer mortality or breast cancer-specific survival, with large variation in design, outcomes, and conclusions (Table 3) [20][21][22][23][24][25][26][27][28]. Olsson et al. was the only one to stratify by mode of detection [21]. Dense breasts appeared to be associated with an increased risk of breast cancer death in symptomatic and non-symptomatic tumours. Gierach et al. adjusted for mode of detection in the fully adjusted model, but they did not present the stratified results [25]. They found no significant difference in breast cancer-specific survival between dense and fatty breasts. This was true for most other studies as well. Olsen et al. found a higher mortality for dense breasts, but case fatality appeared to be lower [20]. Six studies presented estimates that were adjusted for indicators of tumour progression and/or type of treatment. These adjustments are likely to remove (a large part of) the negative effect of missed detection on breast cancer-specific survival as well. Assuming that there is no biological association with tumour aggressiveness, the adjusted estimates are expected to show no effect. Olsson et al. and Chiu et al., however, found the opposite: the already positive hazard ratios further increased after these adjustments [21,27]. Maskarinec et al. and Zhang et al. only presented adjusted estimates, which suggested no effect [26,28]. This was in line with the findings from Gierach and colleagues [25].
Stratification by detection mode removes, or at least reduces, the effects of missed detection, lead time bias, and length time bias. In the resulting groups (i.e. screendetected and interval cancer), there are still other effects that need to be considered. A substantial number of interval tumours in the first year after screening in women with high density may actually be slow-growing tumours that were masked during screening (i.e. false-negative results). The first-year interval tumours in women with a low breast density, on the other hand, could be more aggressive fast-growing tumours. Delayed detection due to high breast density could also still have some influence on observed breast cancer-specific survival in symptomatic cases. Furthermore, excessive opportunistic screening provides a challenge in differentiating between symptomatic and non-symptomatic tumours. Opportunistic screening refers to screening activities that take place outside the organised screening programme. For example, asymptomatic women may ask their GP for a referral to the hospital rather than or even in addition to attending the regular screening programme. These women are classified as non-participants, despite participating in 'opportunistic screening'. Previous research has shown that opportunistic screening occurs frequently [37], and this may have also affected our estimates for the non-participants.
Our study had some limitations. We used a qualitative method to measure breast density, which may have led to non-differential misclassification. Breast pattern was also missing in a relatively large group. This was mainly the result of the introduction of digital mammography. Digital mammography was introduced in the clinic several years before it was introduced in the Dutch screening programme. Breast density pattern was therefore more often missing for 'interval tumours > 24 m' and 'non-participant' tumours than for 'screen-detected' or 'interval ≤ 24 m' tumours. It seems unlikely, however, that the missing data is associated with breast density status or breast cancer death. Another limitation is the lack of potential confounder data. We only have accurate and complete information on age at diagnosis and year of diagnosis. We did, however, stratify by mode of detection, which can also be seen as a proxy method for dealing with lead and length time bias. Finally, a certain degree of informative censoring may have potentially affected our findings.

Conclusions
Detection mode undoubtedly has a strong association with breast cancer death, but there is no clear association between breast density and breast cancer death. This is regardless of detection mode. Breast density is known to affect screening performance and screening decreases the risk of breast cancer death. Further research is still, however, needed on the interplay between breast density, detection mode, and breast cancerspecific survival. Disentangling the different components that influence the observed association between breast density and breast cancer survival will provide new insights into the potential role of breast density in both tumour biology and early detection.