- Research article
- Open Access
Risk factors for breast cancer development by tumor characteristics among women with benign breast disease
Breast Cancer Research volume 23, Article number: 34 (2021)
Among women diagnosed with invasive breast cancer, 30% have a prior diagnosis of benign breast disease (BBD). Thus, it is important to identify factors among BBD patients that elevate invasive cancer risk. In the general population, risk factors differ in their associations by clinical pathologic features; however, whether women with BBD show etiologic heterogeneity in the types of breast cancers they develop remains unknown.
Using a nested case-control study of BBD and breast cancer risk conducted in a community healthcare plan (Kaiser Permanente Northwest), we assessed relationships of histologic features in BBD biopsies and patient characteristics with subsequent breast cancer risk and tested for heterogeneity of associations by estrogen receptor (ER) status, tumor grade, and size. The study included 514 invasive breast cancer cases (median follow-up of 9 years post-BBD diagnosis) and 514 matched controls, diagnosed with proliferative or non-proliferative BBD between 1971 and 2006, with follow-up through mid-2015. Odds ratios (ORs) and 95% confidence intervals (CIs) were obtained using multivariable polytomous logistic regression models.
Breast cancers were predominantly ER-positive (86%), well or moderately differentiated (73%), small (74% < 20 mm), and stage I/II (91%). Compared to patients with non-proliferative BBD, proliferative BBD with atypia conferred increased risk for ER-positive cancer (OR = 5.48, 95% CI = 2.14–14.01) with only one ER-negative case, P-heterogeneity = 0.45. The presence of columnar cell lesions (CCLs) at BBD diagnosis was associated with a 1.5-fold increase in the risk of both ER-positive and ER-negative tumors, with a 2-fold increase (95% CI = 1.21–3.58) observed among postmenopausal women (56%), independent of proliferative BBD status with and without atypia. We did not identify statistically significant differences in risk factor associations by tumor grade or size.
Most tumors that developed after a BBD diagnosis in this cohort were highly treatable low-stage ER-positive tumors. CCL in BBD biopsies may be associated with moderately increased risk, independent of BBD histology, and irrespective of ER status.
It is well established that benign breast diseases (BBDs) increase the risk of breast cancer in women . In the USA, biopsy-confirmed BBD diagnoses are 4 times more common than invasive breast cancers, affecting ~ 1 million women annually. BBDs include a range of pathologies with approximately 30% comprising proliferative BBD lesions (PDWA, proliferative disease, i.e., without atypia) and 65% encompassing non-proliferative benign lesions [2, 3]. Proliferative BBD with atypia (atypical hyperplasia) represents 3–4% of BBD diagnoses, has been associated with a 4–5-fold increased breast cancer (BC) risk [3, 4], and is considered as an indication for possible chemoprevention .
Since women diagnosed with BBD comprise a large proportion of future breast cancer cases (~ 30%), it is important to identify factors associated with subsequent BC . Having undergone a clinically indicated breast biopsy, women with BBD provide the opportunity for evaluation of whether histopathologic features increase BC risk independently of other patient characteristics. Providing accurate invasive BC risk estimates for BBD patients could improve clinical management of this large group of women, including determination of whether watchful waiting, surgery, or chemoprevention should be offered . In a model that predicted overall BC risk for BBD patients, key risk factors were parity/age at first birth and family history of breast cancer, as well as histologic features, including columnar cell lesions (CCLs), radial scars, sclerosing adenosis, and lobular involution . Within the Kaiser BBD study, previous analyses examined associations of lifestyle, reproductive, and pathologic characteristics that are associated with subsequent risk of breast cancer [7, 8]; however, whether risks associated with these factors are similar or vary by the characteristics of the tumors diagnosed among BBD patients is not well understood and has not been previously assessed.
Breast cancer is a heterogeneous disease, with data strongly supporting different associations of risk factors such as parity and genetic susceptibility loci by tumor subtypes defined by hormone receptor status and other clinical/pathologic features [9,10,11,12]. Whether similar etiologic heterogeneity exists for breast cancer arising in BBD patients is unknown, and to date, very few cohort studies have been able to assess this hypothesis [7, 13, 14]. We therefore assessed patient and BBD histopathological features associated with BC risk and evaluated whether relationships varied by tumor characteristics.
We utilized data from a previously designed case-control study of BBD and breast cancer risk , nested within a BBD cohort from Kaiser Permanente Northwest (KPNW), which provides medical care for approximately ½ million members located in Southwest Washington and Northwest Oregon. The source population was comprised of women who had biopsy-confirmed BBD diagnosed between August 3, 1971, and December 31, 2006 (N = 15,395; age range 21–85 years); one control had a BBD diagnosis in 2012. Invasive breast cancer diagnoses were obtained through record linkage with the KPNW Tumor Registry, as previously described . Briefly, cases were defined as women diagnosed with invasive breast cancer at least 1 year after BBD diagnosis and with no prior history of in situ lesions. The follow-up rate for the KPNW Tumor Registry since its inception in 1960 has been excellent, accounting for 98% of patients (living or dead) even if they are no longer health plan members. Controls were women with a biopsy for BBD who were alive but had not developed breast cancer during the same follow-up calendar period as that for the corresponding cases. For each case, a control case was randomly selected using risk-set sampling with replacement, matched at age at diagnosis of BBD (± 1 year) and, implicitly, given the risk-set sampling, on duration of membership of KPNW . In addition, each control was required not to have undergone a mastectomy before the date of diagnosis of breast cancer for her matched case. If a selected control did not have breast tissue for evaluation or had no risk factor information, a replacement control was selected. There were 514 cases and 514 matched controls available for these analyses.
We obtained information on the following variables from the KPNW Tumor Registry: date of invasive breast cancer diagnosis, date of death, tumor histology and behavior using ICD-O coding, grade, AJCC clinical staging variable, tumor size, and lymph node status. Immunohistochemistry (IHC) data on key markers were obtained for ER (which has been measured on cases since the mid-1970s), PR (first measured in 1983), and HER2 status (first measured in 1988). The percentage of tumors analyzed for ER status increased from 0% in 1980 to 18% in 1989, 54% in 1994, and 85% in 2006; HER2 status was available on 43% of tumors until 2001 and on 68% from 2002 onward.
Risk factor information
Risk factor data were captured using routine medical records, which include information on clinic visits, prescriptions, operations, and laboratory testing. These data were linked using the KPNW unique health record numbers to identify each individual member. From these records, data were available for various exposures. We focused on established or suspected breast cancer risk factors, specifically family history of breast cancer, age at menarche, parity, age at first birth, body mass index (BMI) at benign biopsy, menopausal status, and menopausal hormone use. KPNW has a well-characterized mammography screening program, and by 1993, more than 75% of women over 45 had had a screening mammogram . For this reason, we also report characteristics stratified by this date (Supplemental Table 1).
BBD histology was assessed according to the Page classification criteria , as follows: proliferative disease with atypia (ADH), if atypical hyperplasia (either ductal or lobular) was present; proliferative disease without atypia, if epithelial hyperplasia without atypia (either moderate or florid) OR fibroadenoma (either complex-no atypia or complex-atypia) OR sclerosing adenosis OR radial scar OR papilloma was present; non-proliferative, if non-proliferative lesion (either cysts, fibrosis, or apocrine metaplasia) OR mild epithelial hyperplasia without atypia OR fibroadenoma was present. Pathologist assessment of biopsies was blinded to case-control status as previously described . CCL pathology and terminal duct lobular unit (TDLU) involution status were additionally ascertained on digitized BBD H&E images using the Aperio Scanscope system. Furthermore, as involution status is not typically reported in clinical pathology, we obtained more detailed semi-quantitative measures of TDLU involution previously reported to be associated with breast cancer risk among BBD patients (see Supplemental Methods).
We used multiple imputation to impute missing data for risk factors. The largest amount of missingness (21%) was seen for the age at menarche variable (see Supplemental Methods and Appendix A). We did not impute missing data on MHT use, as the percentage of missingness was too large (45%) to allow for stable imputation. All variables that were correlated or associated with the outcome variables or that were potentially related to the missingness of other imputed variables were included in the imputation models (see Supplemental Methods) [17, 18]. Variables were imputed as continuous and later categorized for further analyses. Multiple imputations were performed using IVEware (0.3 version, https://www.src.isr.umich.edu/software/iveware-documentation/iveware-with-sas/). Details of imputation steps are delineated in a flow chart (Supplementary Figure 1 and Appendix A in Supplemental Material). All calculations presented in this paper were conducted for 5 imputed datasets separately; estimates from the different imputed datasets were combined and variances were computed using Rubin’s formula as implemented in SAS 9.4, PROC MIANALYZE.
Descriptive statistics of demographic and tumor characteristics by calendar year of BBD diagnosis and age at breast cancer diagnosis were assessed using chi-squared or Fisher’s exact tests. Conditional logistic regression and unconditional logistic models adjusted for the matching factors yielded similar estimates (Supplemental Table 2). We therefore present odds ratios (ORs) and 95% confidence intervals (CIs) for demographic, reproductive, or tissue factors (explanatory variables) for overall breast cancer risk from unconditional regression models. The unconditional logistic regression models included matching factors, as well as continuous age at BBD diagnosis and follow-up period from BBD diagnosis to breast cancer diagnosis; other key risk factors included in the models were family history of breast cancer in 1st-degree relatives, history of bilateral oophorectomy, and parity. For the main analysis to determine associations by breast cancer subtype (comparing ER, PR, and HER2 status or clinicopathologically defined subtypes), we used polytomous logistic regression models adjusted for matching factors and using the same variables as for the overall BC model. Heterogeneity between factors was assessed using polytomous logistic regression analyses restricted to cases (case-only analyses) with the tumor characteristics (ER, tumor size, and grade) as the outcome variable. Models were also stratified by menopausal status. A P ≤ 0.05 was considered statistically significant and all tests were two-sided. All analyses were performed using SAS V9.4.
Tumor characteristics by patient characteristics, age at breast cancer diagnosis, and BBD calendar year at diagnosis
Characteristics of the BBD patients are detailed in Supplemental Tables 1 and 2. The median age of BBD diagnosis was 51.5 years and over 60% of cases were diagnosed between 1980 and 1999. As mammography screening became more common after 1993, we observed an increased frequency of proliferative disease without atypia (38.5 vs 27.6%) and with atypia (7.3 vs 4.0%, Supplemental Table 1) subsequent to 1993. We also observed increased detection of prevalent breast cancer and a stage shift, with a 12% increase in diagnosed stage I tumors after 1993 (Supplemental Table 1). The median follow-up between BBD and breast cancer diagnosis was 9.0 years (IQR = 4.4, 15.8 years), with a median age at breast cancer diagnosis of 62.7 years. Descriptive characteristics of the BBD patients subsequently diagnosed with breast cancer, overall, and stratified by age at breast cancer diagnosis, are presented in Table 1. Most cases were older than 50 years (87%) at diagnosis. Breast cancers were mostly diagnosed from 1996 to 2005, and 55% of cases were diagnosed within 10 years of their initial BBD diagnosis.
Tumors diagnosed in this population were predominantly small (71.6% ≤ 20 mm), well or moderately differentiated (73.6%), ER-positive (85.9%), PR-positive (71.2%), HER2-negative (79.4%), lymph node-negative (73.6%), and of ductal histology (83.5%). The invasive cases were overwhelmingly of low stage, with < 10% of cases being diagnosed as stage III or IV. As expected, ER status differed significantly by age at breast cancer diagnosis, with a higher proportion (24.5%) of ER-negative breast cancers diagnosed among women ≤ 50 compared to those older than 50 (12.8%). We also assessed whether there were differences in tumor characteristics by calendar period before and during/after 1993 . Of the tumor characteristics evaluated, HER2 status showed a statistically significant difference by BBD diagnosis before and during/after 1993, with a higher proportion of HER2-negative cases diagnosed after 1993 compared to prior years (84.2 vs 73.9%, Supplemental Table 1). After 1993, significantly more of the tumors occurred among women ≥ 50 years of age, were of smaller size (10–20 mm; 44.4% during/after 1993 vs 34.8% before 1993), and had no or mild involution (48.93% during/after 1993 vs 36.56% before 1993, Supplemental Table 1). Histologic grade data were not routinely reported prior to 1993.
Breast cancer risk factors among women with BBD
Association results for all cases combined for established risk factors, including BBD characteristics, are shown in Supplemental Table 2. We found that younger age at first full-term birth and history of bilateral oophorectomy were inversely associated with breast cancer risk, whereas positive family history of breast cancer in a 1st-degree relative, increasing severity of BBD histology, and presence of CCL at BBD diagnosis were associated with increased breast cancer risk (Supplemental Table 2). A representative image of CCL is shown in Fig. 1. CCL with atypia, also known as flat epithelial atypia, is a more severe lesion that has been suggested to be associated with increased risk of breast cancer ; however, we were unable to assess this association in the present study as only 2 controls and 1 case had flat epithelial atypia. Lobular involution, which has been proposed in other BBD patient populations as a key risk factor for subsequent breast cancer [6, 19], was weakly inversely associated with breast cancer risk in our population [complete vs no involution, OR (95% CI) = 0.89 (0.65, 1.24)]. Neither radial scar nor sclerosing adenosis conferred a significant risk among those with proliferative BBD disease (data not shown).
Breast cancer risk by ER status among women with BBD
Risk associations by ER status for patient characteristics and histologic features of BBD are presented in Table 2. Most tumors were ER-positive, which limited the power to detect heterogeneity; as a result, patterns of association observed for overall invasive breast cancer were generally consistent with those observed for ER-positive breast cancer risk. Having an age at first birth < 30 years was associated with reduced risk of ER-positive (OR = 0.69, 95% CI = 0.49–0.98), but not ER-negative (OR = 1.08, 95% CI = 0.51–2.30) breast cancer (P-heterogeneity = 0.24). Compared with patients with non-proliferative BBD, those with proliferative BBD with atypia had a greater than fivefold increased risk for ER-positive disease (OR = 5.48, 95% CI = 2.14–14.01). There was only one ER-negative case; hence, too few cases to provide a reliable estimate in this subgroup. After accounting for BBD histology, the presence of CCLs at BBD diagnosis was associated with a 1.5-fold increased risk for both ER-positive (95% CI = 1.03–2.29) and ER-negative (95% CI = 0.73–3.07) tumors (P-heterogeneity = 0.94).
Breast cancer risk by grade and tumor size among women with BBD
While breast cancer risk associations for patient and histologic characteristics were generally consistent by tumor grade, we found suggestive evidence that higher levels of involution were inversely associated with reduced risk among well-differentiated tumors (complete vs no involution, OR (95% CI) = 0.51 (0.29, 0.90), Supplemental Table 3); this association was not evident among moderately or poorly differentiated tumors, P-het = 0.054. Associations for patient and histologic characteristics by tumor size (< 20 mm indicating small and > 20 mm larger tumors) are presented in Supplemental Table 4. We did not observe any significant heterogeneity by tumor size, although associations of the severity of BBD histology (P < 0.0001) and presence of CCLs (P = 0.038) with increased breast cancer risk were somewhat stronger among patients with smaller tumors.
Breast cancer risk by menopausal status at BBD diagnosis and before and after 1993
We performed additional analyses stratified by menopausal status using factors that have been significantly associated with breast cancer risk, in our study population (Table 3). The association of CCLs with elevated breast cancer risk was most apparent for postmenopausal women (OR = 2.08, 95% CI = 1.21, 3.58), P-het = 0.09.
We did not find evidence of heterogeneity in risk factor associations before and after 1993 when 75% of women over 45 had had a screening mammogram (Supplemental Table 5), suggesting associations are robust.
There are few BBD cohort studies with comprehensive follow-up and detailed characteristics on subsequently diagnosed invasive tumors. Among a large, well-characterized cohort of patients diagnosed with BBD, we expanded upon previous analyses in this population [7, 8] that evaluated associations of well-established breast cancer risk factors and BBD features with breast cancer risk to determine herein if there exists possible etiologic heterogeneity using tumor characteristic data obtained from the long-standing high-quality Kaiser tumor cancer registry . Determining risk associations by tumor subtypes has been increasingly recognized as an important area of research [20, 21]. We comprehensively analyzed histopathologic features of the BBD biopsy and clinicopathologic characteristics of the breast tumors that developed subsequent to BBD diagnosis. We found that most breast cancers diagnosed among BBD patients within this general community healthcare plan were the low-stage, ER-positive tumors that tend to be highly responsive to treatments. Compared with patients with non-proliferative BBD, those with proliferative BBD with atypia had an over fivefold increased risk of ER-positive breast cancer. Our analyses provided limited evidence for heterogeneity in risk factor associations that we evaluated by tumor characteristics. Importantly, histology and presence of CCLs at the time of BBD were independently associated with subsequent breast cancer risk irrespective of ER status, or tumor size or grade. These data provide further support for CCLs as a breast cancer risk factor, especially for postmenopausal BBD patients.
Patient characteristics and breast cancer risk
We extended prior findings by evaluating etiologic heterogeneity for risk factors thought to be most relevant for subsequent risk among BBD patients. Consistent with multiple studies of women diagnosed with sporadic breast cancer [9,10,11,12] and with the Mayo Clinic BBD cohort , we observed that parity/age at first birth was significantly associated with future breast cancer risk. Also consistent with limited data from other BBD cohorts , we found that history of bilateral oophorectomy was associated with reduced breast cancer risk, whereas a positive family history of breast cancer in a 1st-degree relative tended to be associated with increased breast cancer risk. We were able to further evaluate potential etiologic heterogeneity in risk associations, because our BBD study is nested within a single, large, well-defined population with access to health care and with a long-standing tumor tissue registry that has collected data since the 1960s. Although the number of ER-negative tumors was limited, we confirmed patterns of association by ER status that have also been observed among women diagnosed with sporadic breast cancers [9,10,11,12]. For example, early age at first birth showed an inverse association with ER-positive tumors suggesting that reproductive risk factor associations among BBD patients may be similar to those observed in the general population.
It has long been established that ADH is a high-risk precursor lesion, conferring a 4–5-fold increase in the risk of breast cancer development . Indeed, our analysis of BBD patients showed that compared to non-proliferative lesions, a diagnosis of ADH was associated with over 5-fold increased breast cancer risk; however, this association was limited to the risk of ER-positive breast cancer, with little or no risk observed for ER-negative disease. Limited data from cohorts have reported on the relation of ADH by hormone receptor subtype. Consistent with our results, a previous population-based case-control study, CASH, found that history of any benign breast disease was associated only with increased risk of ER-positive luminal A tumors (OR = 1.89, 95% CI 1.43–2.50) . These findings support the hypothesis that benign lesions are more likely to be hormone receptor positive with less genomic instability [23,24,25,26]. Moreover, the finding that ADH might be more relevant for ER-positive breast cancer risk is also consistent with hormonal chemoprevention trials which show a significant reduction in risk in women diagnosed with ADH [5, 27].
Proliferative BBD without atypia
Proliferative disease without atypia is a conglomerate of multiple different pathologies. Radial scars are proliferative lesions that visually appear similar to tumors on mammograms. Pathologically, they are associated with epithelial elements and/or other proliferative lesions such as sclerosing adenosis. Recent analyses in two Swedish cohorts recruited through mammography screening programs from 2001 to 2013 with over 75,000 subjects did not find any significant heterogeneity when evaluating the associations of non-proliferative or proliferative BBD lesions with risk of molecularly defined subtypes of breast cancer; however, a limitation in this study was that there was no separation of proliferative lesions with or without atypia . Our analysis using the Dupont and Page classification of BBD histology showed that PDWA was associated with increased risk of both ER+ and ER− disease. This is consistent with the hypothesis that these lesions serve as precursor lesions in the natural history of breast cancer . The Nurses’ Health Study and the Mayo BBD study showed radial scars to be one of the PDWA lesions associated with breast cancer risk, with an almost 2-fold increased risk [6, 29]. The Nurses’ Health Study cohorts showed an independent association after adjusting for BBD histology. After restricting analyses to BBD patients with proliferative disease, we did not find a significant increase in risk associated with sclerosing adenosis or radial scar. Reasons for these discrepancies could be differences in the study populations as the calendar periods for all three of the cohorts overlap; the Nurses BBD cohort were women who reported a first diagnosis of BBD between 1976 and 1998; the Mayo BBD cohort is a hospital-based referral center that may see more high-risk women but with a similar calendar period to our study with BBD diagnoses occurring between 1967 and 2001. As our study was embedded in a healthcare organization where the women are actively followed and had access to mammography screening, breast cancers detected in our study may be found earlier than in other cohorts which makes direct comparisons difficult.
Columnar cell lesions
While its generally accepted that ADH, lobular neoplasia, and ductal carcinoma in situ (DCIS) are precursor lesions for breast cancer [30, 31], evidence is accumulating that CCLs may also be a precursor, although conveying lower risk particularly among populations with access to mammography screening [32,33,34]. Our data support CCLs as a common risk factor for both ER-positive and ER-negative breast cancers with relative risk estimates of around 1.5 after accounting for BBD histology, for both tumor types. Molecular analyses of CCLs suggest that some alterations occur early and mimic those observed in coincident established precursor lesions such as DCIS [3, 8, 23, 33, 35,36,37,38,39,40,41,42]. There are three other studies that have evaluated CCLs in BBD cohorts, Nashville, Nurses, and Mayo, and our data are consistent with the findings of all three [23, 35, 36, 40]. Only the Nurses BBD cohort obtained tumor characteristics on cases and also did not find any significant differences by ER status or grade, which is consistent with our data. Current management of CCL, according to the Mayo Clinic review, suggests that these patients should be managed with annual clinical breast exams and mammography . Our data suggest CCL might be a risk factor for both ER+ and ER− disease, which has not been observed previously and requires further investigation.
Reduced lobular involution has been previously shown to be a significant factor associated with elevated breast cancer risk among women with BBD [19, 43, 44]. In the present study, we found that involution was weakly inversely related with breast cancer risk for ER+ disease. Further, higher levels of lobular involution were inversely and significantly associated with risk for well-differentiated tumors, a finding which was not observed for moderately or poorly differentiated tumors. As there were relatively few hormone-negative cases in this population, we had limited power to test for differences by ER status. A limitation of this study is the absence of data on menopausal hormone therapy (MHT) use after BBD diagnoses, as MHT uptake was rising at the same time that mammography screening was increasingly being adopted . Given previous studies showing that recent MHT use reduces involution levels among current but not former users , we hypothesize that MHT use post-BBD diagnosis may be an unmeasured negative confounder attenuating involution associations towards the null. Overall, it is unclear whether screening practices, unmeasured MHT use, or/and both might confound observed associations with involution. Additional contemporary studies in other populations with managed health care might help clarify the relationship of involution with future breast cancer risk.
This is a unique cohort of women enrolled in a general community healthcare plan, providing the population access to screening and preventative services that may not be typical for other subsets of the US population. It is not currently known whether women with BBD are followed and screened more closely compared to women without BBD diagnoses. Current data from a study of over 42,000 screened women in Spain found that women with a previous benign breast disease diagnosis had a higher cumulative risk of screen-detected cancer and interval cancers, consistent with data supporting BBD as a risk factor for breast cancer regardless of the mode of detection . Whether women with BBD in our cohort are more likely than the general population to participate in screening is not known but could be the subject of future research.
A limitation of our study is that risk estimates are based on a BBD patient population diagnosed on excisional biopsies during a calendar period spanning the adoption of widespread mammography screening, which became more commonplace in KPNW around 1993; thus, associations may not be reflective of BBD diagnosed in more recent years. Since 1995, advances in breast imaging technologies have resulted in a shift in diagnostic biopsy procedures to core needle biopsies, which currently comprise about 80% of biopsies in the USA . Based on data from the US Breast Cancer Surveillance Consortium, risks associated with high-risk ADH lesions on excisional biopsy were lower when diagnosed via core biopsy (6.7% vs 5.0%), perhaps reflecting the size of the ADH focus . Another limitation as noted earlier was the absence of risk factor data after BBD diagnosis and in particular on MHT use, which has been noted to be prevalent at KPNW at this time  and could have biased the results, especially with respect to the findings for involution. The small number of patients with ER-negative breast cancers is another limitation. Our study had multiple strengths: it was a nested case-control study embedded within a large, well-characterized BBD cohort with lengthy follow-up and access to archival BBD tissues and well-established detailed tumor registry data, the latter aspect allowing us to provide one of the most detailed analyses to date of tumor characteristics of breast cancers diagnosed among the vast majority of patients. Moreover, because we studied women in a healthcare management organization, we could also evaluate temporal changes in access to mammography screening, data that are limited in other cohorts.
Summary and conclusions
Within this BBD cohort, the largest to date with longitudinal follow-up, we provide breast cancer risk factor associations within an HMO patient population . Our data provide further evidence that PDWA is associated with both ER-positive and ER-negative disease. Further, we show CCLs are associated with moderate increases in breast cancer risk, independent of BBD histology, and irrespective of ER status, in agreement with previous studies. Given the predominance of low-stage ER-positive tumors that developed among this cohort, our findings suggest that invasive cancers that develop subsequent to a BBD diagnosis are likely highly treatable with low mortality. Histologic evaluation of BBD biopsies is a promising avenue for the identification of new risk factors for different molecular subtypes of breast cancer and could inform the natural history of disease; however, given complex relationships between screening and diagnosis, along with secular changes in risk factor prevalence, contemporary prospective studies are needed to clarify the relationships of factors that may influence progression of BBD to cancer.
Availability of data and materials
The datasets generated or analyzed during the current study are not publicly available due to data privacy of patients. The authors will make the data available upon reasonable request and with permission of the Kaiser Permanente Center for Health Research in Portland, Oregon.
Dupont WD, Page DL. Risk factors for breast cancer in women with proliferative breast disease. N Engl J Med. 1985;312(3):146–51.
Allison KH, Abraham LA, Weaver DL, Tosteson AN, Nelson HD, Onega T, Geller BM, Kerlikowske K, Carney PA, Ichikawa LE, Buist DS, Elmore JG. Trends in breast biopsy pathology diagnoses among women undergoing mammography in the United States: a report from the Breast Cancer Surveillance Consortium. Cancer. 2015;121(9):1369–78. https://doi.org/10.1002/cncr.29199.
Neal L, Sandhu NP, Hieken TJ, Glazebrook KN, Mac Bride MB, Dilaveri CA, Wahner-Roedler DL, Ghosh K, Visscher DW. Diagnosis and management of benign, atypical, and indeterminate breast lesions detected on core needle biopsy. Mayo Clin Proc. 2014;89(4):536–47.
Hartmann LC, Degnim AC, Santen RJ, Dupont WD, Ghosh K. Atypical hyperplasia of the breast--risk assessment and management options. N Engl J Med. 2015;372(1):78–89.
Cuzick J, Sestak I, Thorat MA. Impact of preventive therapy on the risk of breast cancer among women with benign breast disease. Breast. 2015;24(Suppl 2):S51–5.
Pankratz VS, Degnim AC, Frank RD, Frost MH, Visscher DW, Vierkant RA, Hieken TJ, Ghosh K, Tarabishy Y, Vachon CM, Radisky DC, Hartmann LC. Model for individualized prediction of breast cancer risk after a benign breast biopsy. J Clin Oncol. 2015;33(8):923–9.
Arthur R, Wang Y, Ye K, Glass AG, Ginsberg M, Loudig O, Rohan T. Association between lifestyle, menstrual/reproductive history, and histological factors and risk of breast cancer in women biopsied for benign breast disease. Breast Cancer Res Treat. 2017;165(3):623–31.
Kabat GC, Jones JG, Olson N, Negassa A, Duggan C, Ginsberg M, Kandel RA, Glass AG, Rohan TE. A multi-center prospective cohort study of benign breast disease and risk of subsequent breast cancer. Cancer Causes Control. 2010;21(6):821–8.
Palmer JR, Boggs DA, Wise LA, Ambrosone CB, Adams-Campbell LL, Rosenberg L. Parity and lactation in relation to estrogen receptor negative breast cancer in African American women. Cancer Epidemiol Biomark Prev. 2011;20(9):1883–91.
Palmer JR, Viscidi E, Troester MA, et al. Parity, lactation, and breast cancer subtypes in African American women: results from the AMBER Consortium. J Natl Cancer Inst. 2014;106(10):dju237. https://doi.org/10.1093/jnci/dju237.
Sighoko D, Ogundiran T, Ademola A, Adebamowo C, Chen L, Odedina S, Anetor I, Ndom P, Gakwaya A, Ojengbede O, Huo D, Olopade OI. Breast cancer risk after full-term pregnancies among African women from Nigeria, Cameroon, and Uganda. Cancer. 2015;121(13):2237–43. https://doi.org/10.1002/cncr.29305.
Yang XR, Chang-Claude J, Goode EL, Couch FJ, Nevanlinna H, Milne RL, Gaudet M, Schmidt MK, Broeks A, Cox A, Fasching PA, Hein R, Spurdle AB, Blows F, Driver K, Flesch-Janys D, Heinz J, Sinn P, Vrieling A, Heikkinen T, Aittomaki K, Heikkila P, Blomqvist C, Lissowska J, Peplonska B, Chanock S, Figueroa J, Brinton L, Hall P, Czene K, Humphreys K, Darabi H, Liu J, Van ‘t Veer LJ, van Leeuwen FE, Andrulis IL, Glendon G, Knight JA, Mulligan AM, O’Malley FP, Weerasooriya N, John EM, Beckmann MW, Hartmann A, Weihbrecht SB, Wachter DL, Jud SM, Loehberg CR, Baglietto L, English DR, Giles GG, McLean CA, Severi G, Lambrechts D, Vandorpe T, Weltens C, Paridaens R, Smeets A, Neven P, Wildiers H, Wang X, Olson JE, Cafourek V, Fredericksen Z, Kosel M, Vachon C, Cramp HE, Connley D, Cross SS, Balasubramanian SP, Reed MW, Dork T, Bremer M, Meyer A, Karstens JH, Ay A, Park-Simon TW, Hillemanns P, Arias Perez JI, Menendez Rodriguez P, Zamora P, Benitez J, Ko YD, Fischer HP, Hamann U, Pesch B, Bruning T, Justenhoven C, Brauch H, Eccles DM, Tapper WJ, Gerty SM, Sawyer EJ, Tomlinson IP, Jones A, Kerin M, Miller N, McInerney N, Anton-Culver H, Ziogas A, et al. Associations of breast cancer risk factors with tumor subtypes: a pooled analysis from the Breast Cancer Association Consortium studies. J Natl Cancer Inst. 2011;103(3):250–63.
Hartmann LC, Sellers TA, Frost MH, Lingle WL, Degnim AC, Ghosh K, Vierkant RA, Maloney SD, Pankratz VS, Hillman DW, Suman VJ, Johnson J, Blake C, Tlsty T, Vachon CM, Melton LJ 3rd, Visscher DW. Benign breast disease and the risk of breast cancer. N Engl J Med. 2005;353(3):229–37.
Tamimi RM, Rosner B, Colditz GA. Evaluation of a breast cancer risk prediction model expanded to include category of prior benign breast disease lesion. Cancer. 2010;116(21):4944–53.
Rothman KGS, Lash T. Modern epidemiology. Philadelphia: Lippincott Williams & Wilkins; 2008.
Glass AG, Lacey JV Jr, Carreon JD, Hoover RN. Breast cancer incidence, 1980-2006: combined roles of menopausal hormone therapy, screening mammography, and estrogen receptor status. J Natl Cancer Inst. 2007;99(15):1152–61.
Ali AM, Dawson SJ, Blows FM, Provenzano E, Ellis IO, Baglietto L, Huntsman D, Caldas C, Pharoah PD. Comparison of methods for handling missing data on immunohistochemical markers in survival analysis of breast cancer. Br J Cancer. 2011;104(4):693–9.
Moons KG, Donders RA, Stijnen T, Harrell FE Jr. Using the outcome for imputation of missing predictor values was preferred. J Clin Epidemiol. 2006;59(10):1092–101.
Baer HJ, Collins LC, Connolly JL, Colditz GA, Schnitt SJ, Tamimi RM. Lobule type and subsequent breast cancer risk: results from the Nurses’ Health Studies. Cancer. 2009;115(7):1404–11.
Lee A, Mavaddat N, Wilcox AN, Cunningham AP, Carver T, Hartley S, Babb de Villiers C, Izquierdo A, Simard J, Schmidt MK, Walter FM, Chatterjee N, Garcia-Closas M, Tischkowitz M, Pharoah P, Easton DF, Antoniou AC. BOADICEA: a comprehensive breast cancer risk prediction model incorporating genetic and nongenetic risk factors. Genet Med. 2019;21(8):1708–18.
Mavaddat N, Rebbeck TR, Lakhani SR, Easton DF, Antoniou AC. Incorporating tumour pathology information into breast cancer risk prediction algorithms. Breast Cancer Res. 2010;12(3):R28.
Gaudet MM, Press MF, Haile RW, Lynch CF, Glaser SL, Schildkraut J, Gammon MD, Douglas Thompson W, Bernstein JL. Risk factors by molecular subtypes of breast cancer across a population-based study of women 56 years or younger. Breast Cancer Res Treat. 2011;130(2):587–97.
Boulos FI, Dupont WD, Simpson JF, Schuyler PA, Sanders ME, Freudenthal ME, Page DL. Histologic associations and long-term cancer risk in columnar cell lesions of the breast: a retrospective cohort and a nested case-control study. Cancer. 2008;113(9):2415–21.
Fitzgibbons PL, Henson DE, Hutter RV. Benign breast changes and the risk for subsequent breast cancer: an update of the 1985 consensus statement. Cancer Committee of the College of American Pathologists. Arch Pathol Lab Med. 1998;122(12):1053–5.
Lopez-Garcia MA, Geyer FC, Natrajan R, Kreike B, Mackay A, Grigoriadis A, Reis-Filho JS, Weigelt B. Transcriptomic analysis of tubular carcinomas of the breast reveals similarities and differences with molecular subtype-matched ductal and lobular carcinomas. J Pathol. 2010;222(1):64–75.
Page DL, Dupont WD, Rogers LW, Rados MS. Atypical hyperplastic lesions of the female breast. A long-term follow-up study. Cancer. 1985;55(11):2698–708.
Cuzick J, Sestak I, Cawthorn S, Hamed H, Holli K, Howell A, Forbes JF, Investigators I-I. Tamoxifen for prevention of breast cancer: extended long-term follow-up of the IBIS-I breast cancer prevention trial. Lancet Oncol. 2015;16(1):67–75.
Holm J, Eriksson L, Ploner A, Eriksson M, Rantalainen M, Li J, Hall P, Czene K. Assessment of breast cancer risk factors reveals subtype heterogeneity. Cancer Res. 2017;77(13):3708–17.
Aroner SA, Collins LC, Connolly JL, Colditz GA, Schnitt SJ, Rosner BA, Hankinson SE, Tamimi RM. Radial scars and subsequent breast cancer risk: results from the Nurses’ Health Studies. Breast Cancer Res Treat. 2013;139(1):277–85.
Allred DC, Mohsin SK, Fuqua SA. Histological and biological evolution of human premalignant breast disease. Endocr Relat Cancer. 2001;8(1):47–61.
van Diest PJ. Ductal carcinoma in situ in breast carcinogenesis. J Pathol. 1999;187(4):383–4.
Lee S, Medina D, Tsimelzon A, Mohsin SK, Mao S, Wu Y, Allred DC. Alterations of gene expression in the development of early hyperplastic precursors of breast cancer. Am J Pathol. 2007;171(1):252–62.
Simpson PT, Gale T, Reis-Filho JS, Jones C, Parry S, Sloane JP, Hanby A, Pinder SE, Lee AH, Humphreys S, Ellis IO, Lakhani SR. Columnar cell lesions of the breast: the missing link in breast cancer progression? A morphological and molecular analysis. Am J Surg Pathol. 2005;29(6):734–46.
Sinn HP. Breast cancer precursors: lessons learned from molecular genetics. J Mol Med (Berl). 2009;87(2):113–5.
Aroner SA, Collins LC, Schnitt SJ, Connolly JL, Colditz GA, Tamimi RM. Columnar cell lesions and subsequent breast cancer risk: a nested case-control study. Breast Cancer Res. 2010;12(4):R61.
Boulos FI, Dupont WD, Schuyler PA, Sanders ME, Page DL, Fedda FA, Simpson JF. Clinicopathologic characteristics of carcinomas that develop after a biopsy containing columnar cell lesions: evidence against a precursor role. Cancer. 2012;118(9):2372–7.
Collins LC. Precursor lesions of the low-grade breast neoplasia pathway. Surg Pathol Clin. 2018;11(1):177–97.
Jung MM, Colditz GA, Collins LC, Schnitt SJ, Connolly JL, Tamimi RM. Lifetime physical activity and the incidence of proliferative benign breast disease. Cancer Causes Control. 2011;22(9):1297–305.
Meares AL, Frank RD, Degnim AC, Vierkant RA, Frost MH, Hartmann LC, Winham SJ, Visscher DW. Mucocele-like lesions of the breast: a clinical outcome and histologic analysis of 102 cases. Hum Pathol. 2016;49:33–8.
Said SM, Visscher DW, Nassar A, Frank RD, Vierkant RA, Frost MH, Ghosh K, Radisky DC, Hartmann LC, Degnim AC. Flat epithelial atypia and risk of breast cancer: a Mayo cohort study. Cancer. 2015;121(10):1548–55.
Verschuur-Maes AH, Moelans CB, de Bruin PC, van Diest PJ. Analysis of gene copy number alterations by multiplex ligation-dependent probe amplification in columnar cell lesions of the breast. Cell Oncol (Dordr). 2014;37(2):147–54.
Verschuur-Maes AH, van Gils CH, van den Bosch MA, De Bruin PC, van Diest PJ. Digital mammography: more microcalcifications, more columnar cell lesions without atypia. Modern Pathol. 2011;24(9):1191–7.
Figueroa JD, Pfeiffer RM, Brinton LA, Palakal M, Degnim AC, Radisky D, Hartmann LC, Frost M, Stallings-Mann ML, DP, Visscher D, Sherman ME. Standardized measures of lobular involution and subsequent breast cancer risk among women with benign breast disease. Breast Cancer Res Treat. 2016;159(1):163–72.
Milanese TR, Hartmann LC, Sellers TA, Frost MH, Vierkant RA, Maloney SD, Pankratz VS, Degnim AC, Vachon CM, Reynolds CA, Thompson RA, Melton LJ 3rd, Goode EL, Visscher DW. Age-related lobular involution and risk of breast cancer. J Natl Cancer Inst. 2006;98(22):1600–7.
Figueroa JD, Pfeiffer RM, Patel DA, Linville L, Brinton LA, Gierach GL, Yang XR, Papathomas D, Visscher D, Mies C, Degnim AC, Anderson WF, Hewitt S, Khodr ZG, Clare SE, Storniolo AM, Sherman ME. Terminal duct lobular unit involution of the normal breast: implications for breast cancer etiology. J Natl Cancer Inst. 2014;106(10).
Roman M, Quintana MJ, Ferrer J, Sala M, Castells X. Cumulative risk of breast cancer screening outcomes according to the presence of previous benign breast disease and family history of breast cancer: supporting personalised screening. Br J Cancer. 2017;116(11):1480–5.
Menes TS, Kerlikowske K, Lange J, Jaffer S, Rosenberg R, Miglioretti DL. Subsequent breast cancer risk following diagnosis of atypical ductal hyperplasia on needle biopsy. JAMA Oncol. 2017;3(1):36–41.
We thank Pat Madigan (NCI) for data extraction and quality control checks.
Thomas Rohan is supported by the Breast Cancer Research Foundation (BCRF-20-140). Tumor characteristic data were supported by the Intramural Research Program of the National Cancer Institute, USA. Jonine Figueroa is supported in part by Wellcome Trust grant 207800/Z/17/Z.
Ethics approval and consent to participate
The study used routinely collected data and archival BBD tissue no longer needed for treatment or clinical care. Hence, no verbal or written consent was obtained. All data were anonymized, and the study was approved by the Committee on Clinical Investigations of the Albert Einstein College of Medicine, the Kaiser Permanente Northwest Biospecimen Review Committee, and the National Institutes of Health (NIH) Office of Human Subjects Research (OHSR).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplemental methods. Supplementary Figure 1. Steps of multiple imputation (for more details see Appendix A). Supplementary Table 1. Characteristics of BBD histology and breast cancers among cases by calendar year of BBD and breast cancer diagnoses (N = 514). Supplemental Table 2. Comparison of associations between select patient characteristics and histologic features with breast cancer risk in conditional and unconditional logistic regression models using multiple imputation dataset (N = 1028). Supplemental Table 3. Associations between select patient characteristics and histologic features with breast cancer risk by tumor grade (N = 922). Supplemental Table 4. Associations between select patient characteristics and histologic features with breast cancer risk by tumor size (N = 1013). Supplementary Table 5. Associations between demographic and histologic features with breast cancer risk by BBD calendar year before vs. after 1993 (N = 1028). Appendix A. Distribution of numeric variables before (top) and after (bottom) multiple imputations.
About this article
Cite this article
Figueroa, J.D., Gierach, G.L., Duggan, M.A. et al. Risk factors for breast cancer development by tumor characteristics among women with benign breast disease. Breast Cancer Res 23, 34 (2021). https://doi.org/10.1186/s13058-021-01410-1
- Benign breast disease
- Estrogen receptor
- Columnar cell lesions
- Breast cancer
- Tumor characteristics