Over one million women in the United States receive biopsy diagnoses of benign breast disease (BBD) each year, which confer a 1.5–4.0-fold increase in breast cancer risk. Studies in the general population suggest that nonsteroidal anti-inflammatory agents (NSAID) lower breast cancer risk; however, associations among women with BBD are unknown. We assessed whether NSAID use among women diagnosed with BBD is associated with lower breast cancer risk. Participants included 3,080 women (mean age = 50.3 ± 13.5 years) in the Mayo BBD surgical biopsy cohort diagnosed between January 1, 1992 and December 31, 2001 who completed breast cancer risk factor questionnaires that assessed NSAID use, and whose biopsies underwent detailed pathology review, masked to outcome. Women were followed from date of BBD biopsy to breast cancer diagnosis (main outcome) or censoring (death, prophylactic mastectomy, reduction mammoplasty, lobular carcinoma in situ or last contact). Median follow-up time was 16.4 ± 6.0 years. Incident breast cancer was diagnosed among 312 women over a median follow-up of 9.9 years. Regular non-aspirin NSAID use was associated with lower breast cancer risk [HR = 0.63; 95% confidence interval (CI) = 0.46–0.85; P = 0.002] with trends of lower risk (highest tertiles of use vs. nonuse) for greater number of years used [HR = 0.55; 95% CI = 0.31–0.97; Ptrend = 0.003), days used per month (HR = 0.51; 95% CI = 0.33–0.80; Ptrend = 0.001) and lifetime number of doses taken (HR = 0.53; 95% CI = 0.31–0.89; Ptrend = 0.003). We conclude that nonaspirin NSAID use is associated with statistically significant lower breast cancer risk after a BBD biopsy, including a dose–response effect, suggesting a potential role for NSAIDs in breast cancer prevention among patients with BBD.

Over one million women are diagnosed with benign breast disease (BBD) in the United States annually, which defines a large and expanding pool of patients at increased breast cancer risk (1). Implementing effective breast cancer risk reduction strategies among these women could substantially reduce breast cancer incidence and mortality.

BBD includes morphologically and biologically diverse lesions, several of which may be present concurrently in a single biopsy. To facilitate clinical management, BBD biopsies are classified as nonproliferative, proliferative without atypia or atypical hyperplasia, denoting progressively increasing levels of breast cancer risk. Over 90% of BBD biopsies do not show atypia, and despite conferring only a 1.5- to 2.0-fold increase in breast cancer risk, nonatypical BBD represents the most frequent diagnosis preceding breast cancer (1, 2). Atypical hyperplasia increases breast cancer risk approximately 4-fold, but comprises less than 10% of BBD diagnoses (1, 2). BBD categories define risk among groups of women, but individual risk varies greatly within BBD categories (3), and we hypothesize that several putative carcinogenic mechanisms underlie risk of developing BBD and its progression to breast cancer.

Increased circulating estrogen levels are associated with proliferative BBD, delayed age-related involution of normal lobules from which BBD arises, and increased breast cancer risk (4–6). Endocrine-based prevention lowers risk of estrogen receptor (ER)-positive breast cancer by approximately 50% (7); however, putative nonhormonal carcinogenic mechanisms, including chronic inflammation and immune responses, are not directly targeted by these agents. Sterile inflammation and aging (“inflammaging”) and other processes such as obesity can contribute to a proposed “inflammogenesis” model of breast cancer (8–10). In this model, cyclooxygenases catalyze prostaglandin production, which drives inflammation and fibrosis, and activates multiple downstream carcinogenic mechanisms. In preclinical breast cancer models, activation of inducible COX-2 increases procarcinogenic inflammation; treatment with COX-2 inhibitors, such as nonsteroidal anti-inflammatory agents (NSAID), reduces inflammation and inhibits breast cancer development (11–15). In human studies, inflammatory mediators have been strongly implicated in breast carcinogenesis (16, 17), and elevated COX-2 expression in atypical hyperplasia has been associated with increased breast cancer risk (18). Furthermore, inflammation increases estrogen production via activation of aromatase and regular NSAID use is associated with lower circulating estrogen levels (19). Nonetheless, meta-analyses of studies of average risk women in the general population provide modest evidence that NSAIDs lower breast cancer risk, with significant heterogeneity in results among individual studies and multiple limitations, as noted in these reviews (20–26). However, a notable limitation in nearly all studies of NSAID use and breast cancer risk is a lack of information regarding whether women had procarcinogenic inflammation in their breast tissues when medications were used.

Cyst formation and chronic inflammation are integral components of BBD, as reflected in the historical term “chronic cystic mastitis”. Cysts lined by apocrine cells frequently contain fluid with high levels of androgens, inflammatory mediators and enzymes with both prostaglandin synthesizing (COX-2) and degrading functions (15-hydroxy prostaglandin dehydrogenase; refs. 27–29). Cysts are often associated with duct ectasia, consisting of dilated structures containing fluid and macrophages and encircled by mononuclear cells and fibrosis. An analysis of 1,467 patients with BBD, including 91 women who developed incident breast cancer, found that aspirin use was associated with significantly reduced breast cancer risk [OR = 0.46; 95% confidence interval (CI) = 0.22–0.98], whereas nonaspirin NSAIDs did not show significant protection (OR = 0.82; 95% CI = 0.41–1.62; ref. 30). These data suggest that NSAID use reduces breast cancer risk among patients with BBD, but the small sample size and lack of information on medication doses, follow-up time, and detailed BBD features limit conclusions. Accordingly, we analyzed associations between NSAID use and breast cancer risk in the Mayo Clinic BBD cohort study, which includes detailed breast cancer risk factor data with pathology review and long-term follow-up (2, 31).

Study sample

The Mayo BBD Cohort currently includes 13,455 women, ages 18–85 years, who underwent BBD biopsies between 1967 and 2001 at Mayo Clinic (Rochester, MN; refs. 2, 31). Women who had not been diagnosed with invasive or in situ breast cancer prior to or 6 months following BBD biopsy, had not undergone mastectomy or breast reduction and had not used prevention therapy were eligible. This analysis includes 3,080 (69%) of 4,482 potential participants biopsied between January 1, 1992 and December 31, 2001 who completed the study questionnaire that included self-reported NSAID use; prior to 1992 NSAID use was not queried. Mean age at biopsy (SD) was 50.3 (13.5) years with mean follow-up of 16.4 (6.0) years. Incident breast cancers (including ductal carcinoma in situ) were ascertained from study questionnaires, tumor registry, and review of medical records. The study protocol, including patient contact, written consent in the context of completing questionnaires, and patient follow-up, was approved by the Mayo Clinic Institutional Review Board, and in accord with the Common Rule.

Exposure information

Demographic and clinical characteristics and established breast cancer risk factors were collected using study questionnaires and comprehensive medical record review (2). Family history of breast cancer was categorized as strong, weak, or negative. A strong family history was defined as the patient having (1) at least one first-degree relative with breast cancer diagnosed before age 50 years or (2) two relatives with breast cancer at any age, with at least one being a first-degree relative. Breast cancer risk factors were collected using a baseline study questionnaire and comprehensive medical record review. Participants were asked, “Have you ever regularly take any nonsteroidal anti-inflammatory drugs (see names below), either over-the-counter or by prescription? (exclude occasional use of less than once per month). Specific drug names were listed, followed by fillable tables, in which the name of each drug was listed over three column headers “Average days per month used,” “On days used, average number of pills taken,” and “Total number of years taken.” Each of the three metrics of use was associated with “two-column bubble forms” for entering numerals from 00 to 99.”

Study coordinators classified agents into the following categories: ibuprofen and other prostaglandin synthesis inhibitors; naproxen, celecoxib, rofecoxib and other COX inhibitors; aspirin; other NSAIDs; and non-NSAIDs. Non-NSAID drugs were excluded from analyses.

Histologic examination

Diagnostic BBD biopsies were reviewed in detail and specific proliferative (e.g., ductal hyperplasia) and nonproliferative lesions (e.g., cysts, apocrine metaplasia, duct ectasia, columnar cell change) were recorded using a standardized form, masked to follow-up (DWV; refs. 2, 31, 32). An overall BBD classification based on the highest risk category was also assigned: nonproliferative, proliferative without atypia, or atypical hyperplasia. Level of lobular involution (LI), a feature inversely related to breast cancer risk, was assessed as: noninvoluted (0%); partial involution (1%–75%); or complete involution (>75%; ref. 32).

Statistical analysis

Associations for different non-aspirin NSAID classes were similar (see below) and these data were combined. We examined ever use overall and by number of years used, average number of days used per month, and total number of doses. Total dose was calculated by multiplying together average number of days used per month, average number of pills taken on days used and total number of years taken. Metrics of NSAID use were stratified as never used, and among users, in tertiles. We compared frequencies of key epidemiological features among NSAID users and nonusers using χ2 tests of significance. Among women with breast cancer, we compare NSAID use by breast cancer subtype using χ2tests.

Women were followed from date of BBD biopsy to either date of breast cancer diagnosis or censoring events, including death, prophylactic mastectomy, reduction mammoplasty, diagnosis of lobular carcinoma in situ or last contact. Patients with lobular carcinoma in situ (n = 5) were censored because of routine use of mastectomy during the period of the study. We compared breast cancer risk across use of NSAIDs using Kaplan–Meier curves (Supplementary Fig. S1) and then in Cox proportional hazards regression models with estimated HRs and 95% CIs. We performed unadjusted analyses and multivariate Cox regression models, adjusted for covariates, including: age at biopsy, BBD classification, level of LI, body mass index (BMI; kg/m2) and number of years between biopsy and questionnaire. We assessed dose–response effects related to frequency of NSAID use with breast cancer risk using tests for trend, modeling the ordered frequency of use variable as a one degree-of-freedom term in the Cox regression models. Analyses were carried out overall and by subgroups defined by demographic and clinical variables. We confirmed proportional hazards assumptions by fitting time-by-NSAID use interactions. We ran additional Cox regression analyses to assess, in turn, the potential confounding effects of Breast Imaging, Reporting and Data System (BI-RADs) mammographic density (33) and use of menopausal hormone replacement. We fit Fine and Gray proportional subdistribution hazards models to determine the potential impact of treating death as a competing risk rather than a censoring variable (34).

To account for the potential for NSAID use recall bias, we ran Cox regression sensitivity analyses that excluded follow-up prior to date of questionnaire, such that individuals with censoring or breast cancer events that preceded date of questionnaire were excluded from these analyses. We ran sensitivity analyses for specific drug classes, adjusted for use of other drug classes. All statistical tests were two-sided, and all analyses were carried out using SAS Studio, Release 3.7 (SAS Institute, Inc.).

NSAID use among women with BBD

Of the 3,080 (69%) women who provided NSAID use data, 312 women were diagnosed with incident breast cancer (230 invasive) over a median of 9.9 (SD, 5.9) years of follow-up. Breast cancer risk was not significantly different among women who provided NSAID data and those who did not and were excluded from the analysis (Cox model P = 0.46). Approximately 38% of women reported some level of regular NSAID use with minimal variation by year of BBD biopsy and similar levels of use for specific drug classes (Table 1). Increased total lifetime number of NSAID doses was higher among women 45 years of age and older at biopsy (P < 0.001), obese women (P < 0.001), and women with proliferative BBD (either with or without atypia) versus nonproliferative BBD (P = 0.02). Positive family history was marginally associated with increased NSAID use. Year of biopsy, level of LI, and age at first birth were unrelated to NSAID use.

Table 1.

Distributions of demographic and clinical variables by lifetime NSAID doses.

Lifetime NSAID Doses, n (%)
Missing01–150151–450451+Total
Attribute(n = 35)(n = 1892)(n = 411)(n = 372)(n = 370)(n = 3045)P
Year of biopsy       0.87 
 ≤1993 402 (64.6%) 78 (12.5%) 65 (10.5%) 77 (12.4%) 622 (20.4%)  
 1994–95 359 (60.7%) 80 (13.5%) 82 (13.9%) 70 (11.8%) 591 (19.4%)  
 1996–97 380 (60.9%) 88 (14.1%) 79 (12.7%) 77 (12.3%) 624 (20.5%)  
 1998–99 405 (63.8%) 79 (12.4%) 76 (12.0%) 75 (11.8%) 635 (20.9%)  
 2000–01 346 (60.4%) 86 (15.0%) 70 (12.2%) 71 (12.4%) 573 (18.8%)  
Age at biopsy       <0.001 
 <45 14 710 (65.7%) 138 (12.8%) 111 (10.3%) 122 (11.3%) 1081 (35.5%)  
 45–55 11 502 (56.0%) 145 (16.2%) 130 (14.5%) 119 (13.3%) 896 (29.4%)  
 >55 10 680 (63.7%) 128 (12.0%) 131 (12.3%) 129 (12.1%) 1068 (35.1%)  
Extent of lobular involution       0.75 
 None 339 (63.4%) 74 (13.8%) 63 (11.8%) 59 (11.0%) 535 (20.5%)  
 Partial 16 706 (61.2%) 149 (12.9%) 153 (13.3%) 146 (12.7%) 1154 (44.3%)  
 Complete 561 (61.2%) 129 (14.1%) 106 (11.6%) 120 (13.1%) 916 (35.2%)  
 Missing 286 59 50 45 440  
BBD       0.02 
 NP 21 1091 (64.0%) 216 (12.7%) 192 (11.3%) 207 (12.1%) 1706 (56.4%)  
 PDWA 640 (60.1%) 147 (13.8%) 153 (14.4%) 125 (11.7%) 1065 (35.2%)  
 AH 148 (57.8%) 47 (18.4%) 25 (9.8%) 36 (14.1%) 256 (8.5%)  
 Missing 13 18  
Family history of breast cancer       0.06 
 None 19 981 (64.7%) 189 (12.5%) 187 (12.3%) 160 (10.5%) 1517 (49.9%)  
 Weak 596 (59.5%) 145 (14.5%) 122 (12.2%) 138 (13.8%) 1001 (32.9%)  
 Strong 312 (59.7%) 77 (14.7%) 63 (12.0%) 71 (13.6%) 523 (17.2%)  
 Missing  
BMI at biopsy       <0.001 
 0–21 294 (69.2%) 52 (12.2%) 48 (11.3%) 31 (7.3%) 425 (14.0%)  
 22–25 547 (63.4%) 116 (13.4%) 102 (11.8%) 98 (11.4%) 863 (28.3%)  
 26–29 392 (59.7%) 92 (14.0%) 87 (13.2%) 86 (13.1%) 657 (21.6%)  
 ≥30 11 472 (56.6%) 117 (14.0%) 118 (14.1%) 127 (15.2%) 834 (27.4%)  
 Missing 187 (70.3%) 34 (12.8%) 17 (6.4%) 28 (10.5%) 266 (8.7%)  
Age first live birth/No. Children       0.35 
 <21, 1+ 392 (58.6%) 100 (14.9%) 79 (11.8%) 98 (14.6%) 669 (22.1%)  
 ≥21, 3+ 540 (64.5%) 102 (12.2%) 103 (12.3%) 92 (11.0%) 837 (27.6%)  
 ≥21, 1–2 15 698 (62.9%) 145 (13.1%) 137 (12.3%) 130 (11.7%) 1110 (36.7%)  
 Nulliparous 249 (60.4%) 62 (15.0%) 51 (12.4%) 50 (12.1%) 412 (13.6%)  
 Missing 13 17  
Number of years between biopsy and questionnaire       0.56 
 ≤10 19 743 (61.0%) 174 (14.3%) 158 (13.0%) 143 (11.7%) 1218 (40.0%)  
 11–15 14 1055 (62.8%) 214 (12.7%) 201 (12.0%) 211 (12.6%) 1681 (55.2%)  
 ≥16 94 (64.4%) 23 (15.8%) 13 (8.9%) 16 (11.0%) 146 (4.8%)  
Lifetime NSAID Doses, n (%)
Missing01–150151–450451+Total
Attribute(n = 35)(n = 1892)(n = 411)(n = 372)(n = 370)(n = 3045)P
Year of biopsy       0.87 
 ≤1993 402 (64.6%) 78 (12.5%) 65 (10.5%) 77 (12.4%) 622 (20.4%)  
 1994–95 359 (60.7%) 80 (13.5%) 82 (13.9%) 70 (11.8%) 591 (19.4%)  
 1996–97 380 (60.9%) 88 (14.1%) 79 (12.7%) 77 (12.3%) 624 (20.5%)  
 1998–99 405 (63.8%) 79 (12.4%) 76 (12.0%) 75 (11.8%) 635 (20.9%)  
 2000–01 346 (60.4%) 86 (15.0%) 70 (12.2%) 71 (12.4%) 573 (18.8%)  
Age at biopsy       <0.001 
 <45 14 710 (65.7%) 138 (12.8%) 111 (10.3%) 122 (11.3%) 1081 (35.5%)  
 45–55 11 502 (56.0%) 145 (16.2%) 130 (14.5%) 119 (13.3%) 896 (29.4%)  
 >55 10 680 (63.7%) 128 (12.0%) 131 (12.3%) 129 (12.1%) 1068 (35.1%)  
Extent of lobular involution       0.75 
 None 339 (63.4%) 74 (13.8%) 63 (11.8%) 59 (11.0%) 535 (20.5%)  
 Partial 16 706 (61.2%) 149 (12.9%) 153 (13.3%) 146 (12.7%) 1154 (44.3%)  
 Complete 561 (61.2%) 129 (14.1%) 106 (11.6%) 120 (13.1%) 916 (35.2%)  
 Missing 286 59 50 45 440  
BBD       0.02 
 NP 21 1091 (64.0%) 216 (12.7%) 192 (11.3%) 207 (12.1%) 1706 (56.4%)  
 PDWA 640 (60.1%) 147 (13.8%) 153 (14.4%) 125 (11.7%) 1065 (35.2%)  
 AH 148 (57.8%) 47 (18.4%) 25 (9.8%) 36 (14.1%) 256 (8.5%)  
 Missing 13 18  
Family history of breast cancer       0.06 
 None 19 981 (64.7%) 189 (12.5%) 187 (12.3%) 160 (10.5%) 1517 (49.9%)  
 Weak 596 (59.5%) 145 (14.5%) 122 (12.2%) 138 (13.8%) 1001 (32.9%)  
 Strong 312 (59.7%) 77 (14.7%) 63 (12.0%) 71 (13.6%) 523 (17.2%)  
 Missing  
BMI at biopsy       <0.001 
 0–21 294 (69.2%) 52 (12.2%) 48 (11.3%) 31 (7.3%) 425 (14.0%)  
 22–25 547 (63.4%) 116 (13.4%) 102 (11.8%) 98 (11.4%) 863 (28.3%)  
 26–29 392 (59.7%) 92 (14.0%) 87 (13.2%) 86 (13.1%) 657 (21.6%)  
 ≥30 11 472 (56.6%) 117 (14.0%) 118 (14.1%) 127 (15.2%) 834 (27.4%)  
 Missing 187 (70.3%) 34 (12.8%) 17 (6.4%) 28 (10.5%) 266 (8.7%)  
Age first live birth/No. Children       0.35 
 <21, 1+ 392 (58.6%) 100 (14.9%) 79 (11.8%) 98 (14.6%) 669 (22.1%)  
 ≥21, 3+ 540 (64.5%) 102 (12.2%) 103 (12.3%) 92 (11.0%) 837 (27.6%)  
 ≥21, 1–2 15 698 (62.9%) 145 (13.1%) 137 (12.3%) 130 (11.7%) 1110 (36.7%)  
 Nulliparous 249 (60.4%) 62 (15.0%) 51 (12.4%) 50 (12.1%) 412 (13.6%)  
 Missing 13 17  
Number of years between biopsy and questionnaire       0.56 
 ≤10 19 743 (61.0%) 174 (14.3%) 158 (13.0%) 143 (11.7%) 1218 (40.0%)  
 11–15 14 1055 (62.8%) 214 (12.7%) 201 (12.0%) 211 (12.6%) 1681 (55.2%)  
 ≥16 94 (64.4%) 23 (15.8%) 13 (8.9%) 16 (11.0%) 146 (4.8%)  

Note: Values presented as number (percent). P values calculated using χ2 tests. Strong family history defined as at least one first-degree relative diagnosed with breast cancer prior to age 50 years or two or more relatives diagnosed with breast cancer; other reports of a positive family history were defined as weak. Individuals with missing NSAID values provided information on ever NSAID use but not dosage information.

Abbreviations: Histologic impression: NP, nonproliferative disease; PDWA, proliferative disease without atypia; AH, atypical hyperplasia; BMI, body mass index.

Regular NSAID use and breast cancer risk

Regular NSAID use overall was associated with lower breast cancer risk (HR = 0.76; 95% CI = 0.59–0.99; P = 0.039) with trends suggesting lower risk with increased use (comparing highest tertile for each category of use versus non-use for each metric). Specifically, for number of years used: HR = 0.76; 95% CI = 0.52–1.10; Ptrend = 0.058, for days used per month: HR = 0.49; 95% CI = 0.18–1.33; Ptrend = 0.015 and for total number of doses taken: HR = 0.65, 95% CI = 0.43–0.98; Ptrend = 0.019 (Supplementary Table S1). Lower risk was associated with regular use of nonaspirin NSAIDs (but not aspirin), yielding HR = 0.63 (95% CI = 0.46–0.85; P = 0.002; Table 2). Compared with nonuse, all metrics of non-aspirin NSAID use showed dose–response effects, including years of use [HR = 0.55; 95% CI = 0.31–0.97; Ptrend = 0.003], days used per month (HR = 0.51; 95% CI = 0.33–0.80; Ptrend = 0.001) and total lifetime doses (HR = 0.53; 95% CI = 0.31–0.89; Ptrend = 0.003). Similar patterns for ever use and increased use were found for Ibuprofen and naproxen and their respective related agents, whereas regular aspirin use was not significantly related to breast cancer risk (Supplementary Table S1).

Table 2.

Associations of regular NSAID use with risk of breast cancer among women diagnosed with BBD.

No. eventsPerson-yearsHR (95% CI)aPaHR (95% CI)bPb
Nonaspirin NSAIDs (ibuprofen and other prostaglandin synthesis inhibitors, and naproxen, celecoxib, rofecoxib and other COX inhibitors) 
 Ever used    0.002  0.002 
  No 241 35345 1.00 (ref)  1.00 (ref)  
  Yes 65 14510 0.65 (0.50–0.86)  0.63 (0.46–0.85)  
 Number of years used    0.003c  0.003c 
  0 241 35345 1.00 (ref)  1.00 (ref)  
  1–3 25 4887 0.74 (0.49–1.12)  0.70 (0.45–1.09)  
  4–10 22 5455 0.59 (0.38–0.91)  0.60 (0.37–0.97)  
  11+ 18 4168 0.63 (0.39–1.01)  0.55 (0.31–0.97)  
 Days used per month    <0.001c  0.001c 
  0 241 35345 1.00 (ref)  1.00 (ref)  
  1–9 24 4149 0.84 (0.55–1.28)  0.80 (0.49–1.32)  
  10–29 18 4052 0.65 (0.40–1.04)  0.67 (0.40–1.13)  
  30+ 23 6309 0.53 (0.35–0.81)  0.51 (0.33–0.80)  
 Total number of doses    <0.001c  0.003c 
  0 241 35345 1.00 (ref)  1.00 (ref)  
  1–100 27 4770 0.83 (0.55–1.23)  0.69 (0.44–1.08)  
  101–360 21 4928 0.62 (0.40–0.97)  0.66 (0.41–1.06)  
  361+ 17 4811 0.51 (0.31–0.84)  0.53 (0.31–0.89)  
Aspirin 
 Ever used    0.564  0.426 
  No 250 41363 1.00 (ref)  1.00 (ref)  
  Yes 56 8492 1.09 (0.82–1.46)  0.88 (0.64–1.21)  
Number of years used    0.668c  0.390c 
  0 250 41363 1.00 (ref)  1.00 (ref)  
  1–5 24 3167 1.25 (0.82–1.90)  1.04 (0.66–1.65)  
  6–14 12 2561 0.77 (0.43–1.38)  0.61 (0.32–1.16)  
  15+ 20 2764 1.20 (0.76–1.89)  0.94 (0.57–1.56)  
 Days used per month    0.671c  0.468c 
  0 250 41363 1.00 (ref)  1.00 (ref)  
  1–29 15 2040 1.22 (0.72–2.05)  0.84 (0.46–1.55)  
  30+ 41 6451 1.05 (0.75–1.46)  0.89 (0.62–1.27)  
 Total number of doses    0.995c  0.293c 
  0 250 41363 1.00 (ref)  1.00 (ref)  
  1–150 23 2986 1.27 (0.83–1.95)  0.95 (0.59–1.55)  
  151–360 18 2581 1.15 (0.71–1.86)  0.99 (0.59–1.65)  
  361+ 15 2924 0.85 (0.50–1.42)  0.70 (0.40–1.24)  
No. eventsPerson-yearsHR (95% CI)aPaHR (95% CI)bPb
Nonaspirin NSAIDs (ibuprofen and other prostaglandin synthesis inhibitors, and naproxen, celecoxib, rofecoxib and other COX inhibitors) 
 Ever used    0.002  0.002 
  No 241 35345 1.00 (ref)  1.00 (ref)  
  Yes 65 14510 0.65 (0.50–0.86)  0.63 (0.46–0.85)  
 Number of years used    0.003c  0.003c 
  0 241 35345 1.00 (ref)  1.00 (ref)  
  1–3 25 4887 0.74 (0.49–1.12)  0.70 (0.45–1.09)  
  4–10 22 5455 0.59 (0.38–0.91)  0.60 (0.37–0.97)  
  11+ 18 4168 0.63 (0.39–1.01)  0.55 (0.31–0.97)  
 Days used per month    <0.001c  0.001c 
  0 241 35345 1.00 (ref)  1.00 (ref)  
  1–9 24 4149 0.84 (0.55–1.28)  0.80 (0.49–1.32)  
  10–29 18 4052 0.65 (0.40–1.04)  0.67 (0.40–1.13)  
  30+ 23 6309 0.53 (0.35–0.81)  0.51 (0.33–0.80)  
 Total number of doses    <0.001c  0.003c 
  0 241 35345 1.00 (ref)  1.00 (ref)  
  1–100 27 4770 0.83 (0.55–1.23)  0.69 (0.44–1.08)  
  101–360 21 4928 0.62 (0.40–0.97)  0.66 (0.41–1.06)  
  361+ 17 4811 0.51 (0.31–0.84)  0.53 (0.31–0.89)  
Aspirin 
 Ever used    0.564  0.426 
  No 250 41363 1.00 (ref)  1.00 (ref)  
  Yes 56 8492 1.09 (0.82–1.46)  0.88 (0.64–1.21)  
Number of years used    0.668c  0.390c 
  0 250 41363 1.00 (ref)  1.00 (ref)  
  1–5 24 3167 1.25 (0.82–1.90)  1.04 (0.66–1.65)  
  6–14 12 2561 0.77 (0.43–1.38)  0.61 (0.32–1.16)  
  15+ 20 2764 1.20 (0.76–1.89)  0.94 (0.57–1.56)  
 Days used per month    0.671c  0.468c 
  0 250 41363 1.00 (ref)  1.00 (ref)  
  1–29 15 2040 1.22 (0.72–2.05)  0.84 (0.46–1.55)  
  30+ 41 6451 1.05 (0.75–1.46)  0.89 (0.62–1.27)  
 Total number of doses    0.995c  0.293c 
  0 250 41363 1.00 (ref)  1.00 (ref)  
  1–150 23 2986 1.27 (0.83–1.95)  0.95 (0.59–1.55)  
  151–360 18 2581 1.15 (0.71–1.86)  0.99 (0.59–1.65)  
  361+ 15 2924 0.85 (0.50–1.42)  0.70 (0.40–1.24)  

Note: Category-specific numbers of events and person-years may not sum to overall numbers due to missing values for some variables.

aUnadjusted Cox proportional hazards analysis.

bCox proportional hazard analyses adjusted for age at biopsy, histologic impression, extent of lobular involution, BMI and number of years between biopsy and completion of questionnaire.

cTest for trend.

Nonaspirin NSAID use and breast cancer risk by selected patient and biopsy characteristics

Reduced breast cancer risk with regular non-aspirin NSAID use was found among women with nonproliferative BBD (HR = 0.52; 95% CI = 0.31–0.87; P = 0.013); associations for proliferative BBD with or without atypia were similar in direction, but failed to reach significance (Table 3). Similarly, significantly decreased risk with regular nonaspirin NSAID use was found among women ≤50 years of age (HR = 0.56; 95% CI = 0.36–0.87; P = 0.010), parous women (HR = 0.64; 95% CI = 0.46–0.89; P = 0.007), and obese women (HR = 0.45; 95% CI = 0.0.25–0.80; P = 0.007), whereas risks were nonsignificantly reduced among women who were >50 years of age, nulliparous, or nonobese, albeit in some analyses that included few events. Lower breast cancer risk was associated with several features included in the spectrum of nonproliferative BBD (cysts, apocrine metaplasia, duct ectasia, and columnar cell change) and were markedly lower for women with ≥3 such lesions (HR = 0.33; 95% CI = 0.17–0.65; P = 0.001). Among patients with BBD who later developed breast cancer, risks for NSAID use did not vary by estrogen receptor, progesterone receptor, or nodal status of the breast cancer (Supplementary Table S2).

Table 3.

Associations of regular nonaspirin NSAID use with BC risk in relation to demographic and clinical variables.

Stratification variableEver usePerson yearsNo. censoredNo. eventsHR (95% CI)aPa
Nonproliferative disease No 20,397 1130 111 1.00 (ref) 0.013 
 Yes 8,263 441 24 0.52 (0.31–0.87)  
Proliferative disease No 14,736 828 128 1.00 (ref) 0.083 
 Yes 6,178 324 41 0.72 (0.49–1.04)  
Proliferative disease without atypia No 12,158 686 86 1.00 (ref) 0.126 
 Yes 5,039 266 27 0.69 (0.43–1.11)  
Atypical hyperplasia No 2,578 142 42 1.00 (ref) 0.425 
 Yes 1,139 58 14 0.76 (0.39–1.49)  
Age ≤ 50 No 19,048 1008 106 1.00 (ref) 0.010 
 Yes 9,029 471 33 0.56 (0.36–0.87)  
Age > 50 No 16,297 961 135 1.00 (ref) 0.135 
 Yes 5,481 299 32 0.73 (0.48–1.10)  
Parous No 30,883 1714 206 1.00 (ref) 0.007 
 Yes 12,444 656 57 0.64 (0.46–0.89)  
Nulliparous No 4,463 255 35 1.00 (ref) 0.288 
 Yes 2,066 114 0.63 (0.26–1.49)  
BMI < 30 No 23,179 1283 162 1.00 (ref) 0.120 
 Yes 8,664 453 47 0.75 (0.53–1.08)  
BMI ≥ 30 No 8,960 486 75 1.00 (ref) 0.007 
 Yes 4,766 254 17 0.45 (0.25–0.80)  
Time between biopsy and survey < 12 years No 17,020 1028 124 1.00 (ref) 0.141 
 Yes 7,247 420 40 0.74 (0.49–1.11)  
Time between biopsy and survey ≥ 12 years No 18,325 941 117 1.00 (ref) 0.013 
 Yes 7,263 350 25 0.56 (0.35–0.89)  
Cyst – Present No 13,810 749 114 1.00 (ref) 0.008 
 Yes 6,168 325 29 0.56 (0.36–0.86)  
Cyst – Absent No 21,428 1215 125 1.00 (ref) 0.146 
 Yes 8,287 441 36 0.73 (0.48–1.12)  
Duct ectasia No 3,748 204 32 1.00 (ref) 0.038 
 Present Yes 1,601 84 0.36 (0.13–0.95)  
Duct ectasia No 31,466 1758 207 1.00 (ref) 0.028 
 Absent Yes 12,820 680 60 0.70 (0.51–0.96)  
Apocrine metaplasia No 10,477 567 80 1.00 (ref) 0.002 
 Present Yes 4,428 235 16 0.41 (0.23–0.73)  
Apocrine metaplasia No 24,761 1397 159 1.00 (ref) 0.165 
 Absent Yes 10,027 531 49 0.78 (0.54–1.11)  
Columnar cell alteration No 13,272 735 115 1.00 (ref) 0.040 
 Present Yes 6,033 312 37 0.66 (0.45–0.98)  
Columnar cell alteration No 21,966 1229 124 1.00 (ref) 0.028 
 Absent Yes 8,422 454 28 0.58 (0.36–0.94)  
No inflammatory lesionsb No 15,430 868 83 1.00 (ref) 0.263 
 Yes 5,895 315 23 0.73 (0.43–1.26)  
At least one inflammatory lesionb No 19,809 1096 156 1.00 (ref) 0.006 
 Yes 8,542 450 42 0.60 (0.42–0.87)  
Number inflammatory lesions = 1b No 7,045 405 45 1.00 (ref) 0.250 
 Yes 2,862 151 13 0.67 (0.34–1.33)  
Number inflammatory lesions = 2b No 5,661 310 52 1.00 (ref) 0.959 
 Yes 2,427 129 18 1.02 (0.57–1.81)  
Number inflammatory lesions = 3 or moreb No 7,103 381 59 1.00 (ref) 0.001 
 Yes 3,252 170 11 0.33 (0.17–0.65)  
Stratification variableEver usePerson yearsNo. censoredNo. eventsHR (95% CI)aPa
Nonproliferative disease No 20,397 1130 111 1.00 (ref) 0.013 
 Yes 8,263 441 24 0.52 (0.31–0.87)  
Proliferative disease No 14,736 828 128 1.00 (ref) 0.083 
 Yes 6,178 324 41 0.72 (0.49–1.04)  
Proliferative disease without atypia No 12,158 686 86 1.00 (ref) 0.126 
 Yes 5,039 266 27 0.69 (0.43–1.11)  
Atypical hyperplasia No 2,578 142 42 1.00 (ref) 0.425 
 Yes 1,139 58 14 0.76 (0.39–1.49)  
Age ≤ 50 No 19,048 1008 106 1.00 (ref) 0.010 
 Yes 9,029 471 33 0.56 (0.36–0.87)  
Age > 50 No 16,297 961 135 1.00 (ref) 0.135 
 Yes 5,481 299 32 0.73 (0.48–1.10)  
Parous No 30,883 1714 206 1.00 (ref) 0.007 
 Yes 12,444 656 57 0.64 (0.46–0.89)  
Nulliparous No 4,463 255 35 1.00 (ref) 0.288 
 Yes 2,066 114 0.63 (0.26–1.49)  
BMI < 30 No 23,179 1283 162 1.00 (ref) 0.120 
 Yes 8,664 453 47 0.75 (0.53–1.08)  
BMI ≥ 30 No 8,960 486 75 1.00 (ref) 0.007 
 Yes 4,766 254 17 0.45 (0.25–0.80)  
Time between biopsy and survey < 12 years No 17,020 1028 124 1.00 (ref) 0.141 
 Yes 7,247 420 40 0.74 (0.49–1.11)  
Time between biopsy and survey ≥ 12 years No 18,325 941 117 1.00 (ref) 0.013 
 Yes 7,263 350 25 0.56 (0.35–0.89)  
Cyst – Present No 13,810 749 114 1.00 (ref) 0.008 
 Yes 6,168 325 29 0.56 (0.36–0.86)  
Cyst – Absent No 21,428 1215 125 1.00 (ref) 0.146 
 Yes 8,287 441 36 0.73 (0.48–1.12)  
Duct ectasia No 3,748 204 32 1.00 (ref) 0.038 
 Present Yes 1,601 84 0.36 (0.13–0.95)  
Duct ectasia No 31,466 1758 207 1.00 (ref) 0.028 
 Absent Yes 12,820 680 60 0.70 (0.51–0.96)  
Apocrine metaplasia No 10,477 567 80 1.00 (ref) 0.002 
 Present Yes 4,428 235 16 0.41 (0.23–0.73)  
Apocrine metaplasia No 24,761 1397 159 1.00 (ref) 0.165 
 Absent Yes 10,027 531 49 0.78 (0.54–1.11)  
Columnar cell alteration No 13,272 735 115 1.00 (ref) 0.040 
 Present Yes 6,033 312 37 0.66 (0.45–0.98)  
Columnar cell alteration No 21,966 1229 124 1.00 (ref) 0.028 
 Absent Yes 8,422 454 28 0.58 (0.36–0.94)  
No inflammatory lesionsb No 15,430 868 83 1.00 (ref) 0.263 
 Yes 5,895 315 23 0.73 (0.43–1.26)  
At least one inflammatory lesionb No 19,809 1096 156 1.00 (ref) 0.006 
 Yes 8,542 450 42 0.60 (0.42–0.87)  
Number inflammatory lesions = 1b No 7,045 405 45 1.00 (ref) 0.250 
 Yes 2,862 151 13 0.67 (0.34–1.33)  
Number inflammatory lesions = 2b No 5,661 310 52 1.00 (ref) 0.959 
 Yes 2,427 129 18 1.02 (0.57–1.81)  
Number inflammatory lesions = 3 or moreb No 7,103 381 59 1.00 (ref) 0.001 
 Yes 3,252 170 11 0.33 (0.17–0.65)  

Note: Category-specific numbers of events and person-years may not sum to overall numbers due to missing values for some variables.

Abbreviation: BMI, body mass index.

aCox proportional hazard analyses adjusted for age, histologic impression, extent of lobular involution, BMI and number of years between biopsy and completion of questionnaire.

bNonproliferative BBD lesions include: cyst, duct ectasia, apocrine metaplasia, and columnar cell alteration.

Sensitivity analyses

Associations of ever versus never regular non-aspirin NSAID use with risk of breast cancer were similar when subsetting to proliferative BBD and atypical hyperplasia (see Table 3) and in analyses adjusted additionally for breast density (HR = 0.62; 95% CI, 0.46–0.84) and menopausal hormone use (HR = 0.61; 95% CI, 0.45–0.83). Accounting for death as a competing risk did not substantively change results (HR = 0.65; 95% CI, 0.48–0.88). Analyses that excluded follow-up prior to date of completion of the questionnaire (which effectively removed 197 participants who were diagnosed with breast cancer prior to completing the questionnaire and 604 participants right-censored prior to completion) yielded similar results (HR 0.57; 95% CI, 0.35–0.93). See Supplementary Table S3 for full results. Associations persisted when breast cancer risk related to a specific agent was adjusted for use of other agents.

Well-tolerated targeted breast cancer prevention strategies could offer great benefits for individual women diagnosed with BBD and reduce the burden of breast cancer in the population. We report that women in the Mayo BBD cohort who regularly used nonaspirin NSAIDs experienced significantly lower breast cancer risk, with suggestively greater effects among subsets of patients.

Regular use of non-aspirin NSAIDs was associated with 37% lower breast cancer risk in our analysis, with greater effects in relation to increased years or frequency of use and greater number of lifetime doses consumed. Sensitivity analyses that excluded women who were diagnosed with breast cancer prior to reporting NSAID use or for whom data were provided by next-of-kin yielded similar results, limiting concerns of bias. Adjusting for mammographic density, use of menopausal hormone therapy or use of multiple different NSAIDs did not alter these conclusions. Risk associations did not vary by estrogen receptor or nodal status of breast cancers.

In a prior report, with limited dosing information, aspirin use among women with BBD showed a significant reduction in breast cancer risk, whereas nonaspirin NSAID use was associated with nonsignificantly lower risk (30). We did not identify significant risk associations for aspirin use; however, 54.1% of aspirin users in our analysis were 55 years of age or older compared with 21.4% of nonaspirin users, and we did not collect data on dosage of the tablets consumed; thus, if risk reduction with aspirin use is greatest among younger women or only realized with high dose formulations, a protective effect may not have been found. Consistent with this view, risk associations were greatest among nonaspirin NSAID users who were less than 50 years of age. Several features of the premenopausal period are associated with inflammation, including menses, which is associated with increased circulating inflammatory cytokines (35), and postpartum involution, in which breast tissue undergoes dramatic remodeling. In preclinical models, postpartum involution drives the development of aggressive breast cancer and this effect can be mitigated with NSAIDs (14). In humans, dysregulated post-partum inflammation may contribute to transiently elevated breast cancer risk after childbirth (36). Effects of nonaspirin NSAID use on breast cancer risk were also greater among obese versus nonobese women. Obesity is associated with low-grade chronic inflammation, which may manifest in adipose rich tissues, such as the breast, and increase breast cancer risk through hormonal and nonhormonal pathways (37). Mechanistic studies show that sera from obese patients induce dramatic increases in macrophage production of COX-2, which can be linked to increased aromatase expression in preadipocytes, and can be inhibited with the NSAID celecoxib (38).

Meta-analyses among average risk women have generally found that aspirin and non-aspirin NSAID use are associated with lower breast cancer risk; however, results are variable (20–25) and magnitudes of effects are generally modest. Misclassification of drug use is an important concern in observational epidemiologic studies evaluating breast risk; however, random misclassification would generally result in under-estimation of true relationships and over-reporting among cases (recall bias) would likely produce a similar effect. Furthermore, NSAID use in large epidemiologic studies cannot be temporally related to inflammation in the breast, which would be anticipated to dilute protective effects. Previously, we reported that immune cells are increased in BBD biopsies versus normal tissues donated for research, and that specific immune responses may predict breast cancer risk among women with BBD (18, 39, 40), suggesting that targeting inflammation may have potential utility for prevention among women with BBD.

Our data suggest that lower breast cancer risk with NSAID use was more evident for non-proliferative BBD, which is the most common category of BBD diagnosis. Nonproliferative BBD frequently accompanies proliferative BBD or atypical hyperplasia, although pathology reports may only report the highest risk lesion. Apocrine cysts are among the most frequent manifestations of nonproliferative BBD, and while cysts are not breast cancer precursors, cyst fluid is rich in hormones and inflammatory mediators, which may create a microenvironment conducive to breast cancer development (27–29, 41). Furthermore, cysts may form and collapse, resulting in inflammation and scarring, even when no longer recognizable microscopically. Women with ≥3 nonproliferative histopathologic BBD lesions in BBD biopsies, a potential marker of diffuse bilateral disease, showed markedly lower breast cancer risk with nonaspirin NSAID use (HR = 0.33; 95% CI = 0.17–0.65; P = 0.001).

Our study leverages strengths of the Mayo BBD cohort, which includes collection of breast cancer risk factors, tissues and follow-up. Our ability to demonstrate a dose–response relationship between non-aspirin NSAID use and lower breast cancer risk among women with BBD is an additional strength. However, this population is comprised almost exclusively of white women, which limits generalizability. Furthermore, medication use was self-reported at a single time point and did not distinguish use of low-dose aspirin from higher doses, which may have limited our ability to identify a significant association for use of this drug. In addition, stratified analyses assessing NSAID associations in certain subsets of women resulted in modest sample sizes and decreased statistical power. We did not collect information about use of statins or bisphosphonates, which are frequently used at older ages, and could affect breast cancer risk. Results of multivariable analyses adjusted for potential confounders, produced only modest effects on risk estimates. Although anti-inflammatory agents are used sporadically for acute illnesses or injuries, our analysis focused on regular use and found strongest effects for higher cumulative doses.

In summary, our analysis of women with BBD demonstrates that regular use of nonaspirin NSAIDs is associated with an approximately 37% reduction in breast cancer risk. Given that BBD is associated with increased breast cancer risk, and frequently includes inflammatory lesions in breast tissues, our findings suggest that NSAIDs may represent an important prevention strategy, and therefore merits further clinical and mechanistic investigation. We propose that clinical prevention trials evaluating the efficacy of NSAIDs for breast cancer prevention among women with BBD that incorporate correlative studies to examine interactions between inflammation and estrogen pathways in breast tissues are warranted. A recently completed phase II trial of celecoxib 400 mg twice daily for 6 months, demonstrated good tolerance with favorable modulation of selected biomarkers (42).

M.H. Frost reports grants from NCI during the conduct of the study. C.M. Vachon reports grants from NIH/NCI during the conduct of the study; grants from Grail "STRIVE mammography study" outside the submitted work. No potential conflicts of interest were disclosed by the other authors.

M.E. Sherman: Conceptualization, resources, formal analysis, validation, methodology, writing-original draft, project administration, writing-review and editing. R.A. Vierkant: Conceptualization, formal analysis, supervision, writing-original draft, project administration, writing-review and editing. S. Kaggal: Formal analysis. T.L. Hoskin: Formal analysis, writing-review and editing. M.H. Frost: Data curation, supervision, project administration, writing-review and editing. L. Denison: Writing-review and editing. D.W. Visscher: Formal analysis, methodology. J.M. Carter: Writing-original draft, writing-review and editing. S.J. Winham: Formal analysis, supervision, writing-original draft, project administration, writing-review and editing. M.R. Jensen: Formal analysis. D.C. Radisky: Formal analysis, funding acquisition, writing-original draft, writing-review and editing. C.M. Vachon: Formal analysis, writing-original draft, writing-review and editing. A.C. Degnim: Conceptualization, resources, supervision, funding acquisition, writing-original draft, project administration, writing-review and editing.

This work was supported by the NCI of the NIH (CA229811 to M.E. Sherman and A.C. Degnim, and CA1187112 to A.C. Degnim).

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Morrow
M
,
Schnitt
SJ
,
Norton
L
. 
Current management of lesions associated with an increased risk of breast cancer
.
Nat Rev Clin Oncol
2015
;
12
:
227
38
.
2.
Hartmann
LC
,
Sellers
TA
,
Frost
MH
,
Lingle
WL
,
Degnim
AC
,
Ghosh
K
, et al
Benign breast disease and the risk of breast cancer
.
N Engl J Med
2005
;
353
:
229
37
.
3.
Pankratz
VS
,
Degnim
AC
,
Frank
RD
,
Frost
MH
,
Visscher
DW
,
Vierkant
RA
, et al
Model for individualized prediction of breast cancer risk after a benign breast biopsy
.
J Clin Oncol
2015
;
33
:
923
9
.
4.
Khodr
ZG
,
Sherman
ME
,
Pfeiffer
RM
,
Gierach
GL
,
Brinton
LA
,
Falk
RT
, et al
Circulating sex hormones and terminal duct lobular unit involution of the normal breast
.
Cancer Epidemiol Biomarkers Prev
2014
;
23
:
2765
73
.
5.
Schairer
C
,
Hill
D
,
Sturgeon
SR
,
Fears
T
,
Mies
C
,
Ziegler
RG
, et al
Serum concentrations of estrogens, sex hormone binding globulin, and androgens and risk of breast hyperplasia in postmenopausal women
.
Cancer Epidemiol Biomarkers Prev
2005
;
14
:
1660
5
.
6.
Key
T
,
Appleby
P
,
Barnes
I
,
Reeves
G
,
Endogenous
H
,
Breast Cancer Collaborative Group
. 
Endogenous sex hormones and breast cancer in postmenopausal women: reanalysis of nine prospective studies
.
J Natl Cancer Inst
2002
;
94
:
606
16
.
7.
Nelson
HD
,
Fu
R
,
Zakher
B
,
Pappas
M
,
McDonagh
M
. 
Medication use for the risk reduction of primary breast cancer in women: updated evidence report and systematic review for the US Preventive Services Task Force
.
JAMA
2019
;
322
:
868
86
.
8.
Harris
RE
,
Casto
BC
,
Harris
ZM
. 
Cyclooxygenase-2 and the inflammogenesis of breast cancer
.
World J Clin Oncol
2014
;
5
:
677
92
.
9.
Kennedy
BM
,
Harris
RE
. 
Cyclooxygenase and lipoxygenase gene expression in the inflammogenesis of breast cancer
.
Inflammopharmacology
2018
May 7.
doi:10.1007/s10787-018-0489-6. Epub ahead of print. PMID: 29736687.
10.
Fresques
T
,
Zirbes
A
,
Shalabi
S
,
Samson
S
,
Preto
S
,
Stampfer
MR
, et al
Breast tissue biology expands the possibilities for prevention of age-related breast cancers
.
Front Cell Dev Biol
2019
;
7
:
174
.
11.
Wallace
TR
,
Tarullo
SE
,
Crump
LS
,
Lyons
TR
. 
Studies of postpartum mammary gland involution reveal novel pro-metastatic mechanisms
.
J Cancer Metastasis Treat
2019
;
5
:
9
.
12.
Jindal
S
,
Gao
D
,
Bell
P
,
Albrektsen
G
,
Edgerton
SM
,
Ambrosone
CB
, et al
Postpartum breast involution reveals regression of secretory lobules mediated by tissue-remodeling
.
Breast Cancer Res
2014
;
16
:
R31
.
13.
Martinson
HA
,
Jindal
S
,
Durand-Rougely
C
,
Borges
VF
,
Schedin
P
. 
Wound healing-like immune program facilitates postpartum mammary gland involution and tumor progression
.
Int J Cancer
2015
;
136
:
1803
13
.
14.
Lyons
TR
,
O'Brien
J
,
Borges
VF
,
Conklin
MW
,
Keely
PJ
,
Eliceiri
KW
, et al
Postpartum mammary gland involution drives progression of ductal carcinoma in situ through collagen and COX-2
.
Nat Med
2011
;
17
:
1109
15
.
15.
Schedin
P
,
O'Brien
J
,
Rudolph
M
,
Stein
T
,
Borges
V
. 
Microenvironment of the involuting mammary gland mediates mammary cancer progression
.
J Mammary Gland Biol Neoplasia
2007
;
12
:
71
82
.
16.
Suman
S
,
Sharma
PK
,
Rai
G
,
Mishra
S
,
Arora
D
,
Gupta
P
, et al
Current perspectives of molecular pathways involved in chronic inflammation-mediated breast cancer
.
Biochem Biophys Res Commun
2016
;
472
:
401
9
.
17.
King
J
,
Mir
H
,
Singh
S
. 
Association of cytokines and chemokines in pathogenesis of breast cancer
.
Prog Mol Biol Transl Sci
2017
;
151
:
113
36
.
18.
Visscher
DW
,
Pankratz
VS
,
Santisteban
M
,
Reynolds
C
,
Ristimaki
A
,
Vierkant
RA
, et al
Association between cyclooxygenase-2 expression in atypical hyperplasia and risk of breast cancer
.
J Natl Cancer Inst
2008
;
100
:
421
7
.
19.
Gates
MA
,
Tworoger
SS
,
Eliassen
AH
,
Missmer
SA
,
Hankinson
SE
. 
Analgesic use and sex steroid hormone concentrations in postmenopausal women
.
Cancer Epidemiol Biomarkers Prev
2010
;
19
:
1033
41
.
20.
Qiao
Y
,
Yang
T
,
Gan
Y
,
Li
W
,
Wang
C
,
Gong
Y
, et al
Associations between aspirin use and the risk of cancers: a meta-analysis of observational studies
.
BMC Cancer
2018
;
18
:
288
.
21.
Lu
L
,
Shi
L
,
Zeng
J
,
Wen
Z
. 
Aspirin as a potential modality for the chemoprevention of breast cancer: a dose-response meta-analysis of cohort studies from 857,831 participants
.
Oncotarget
2017
;
8
:
40389
401
.
22.
Zhong
S
,
Chen
L
,
Zhang
X
,
Yu
D
,
Tang
J
,
Zhao
J
. 
Aspirin use and risk of breast cancer: systematic review and meta-analysis of observational studies
.
Cancer Epidemiol Biomarkers Prev
2015
;
24
:
1645
55
.
23.
Algra
AM
,
Rothwell
PM
. 
Effects of regular aspirin on long-term cancer incidence and metastasis: a systematic comparison of evidence from observational studies versus randomised trials
.
Lancet Oncol
2012
;
13
:
518
27
.
24.
Zhao
YS
,
Zhu
S
,
Li
XW
,
Wang
F
,
Hu
FL
,
Li
DD
, et al
Association between NSAIDs use and breast cancer risk: a systematic review and meta-analysis
.
Breast Cancer Res Treat
2009
;
117
:
141
50
.
25.
Takkouche
B
,
Regueira-Mendez
C
,
Etminan
M
. 
Breast cancer and use of nonsteroidal anti-inflammatory drugs: a meta-analysis
.
J Natl Cancer Inst
2008
;
100
:
1439
47
.
26.
de Pedro
M
,
Baeza
S
,
Escudero
MT
,
Dierssen-Sotos
T
,
Gomez-Acebo
I
,
Pollan
M
, et al
Effect of COX-2 inhibitors and other non-steroidal inflammatory drugs on breast cancer risk: a meta-analysis
.
Breast Cancer Res Treat
2015
;
149
:
525
36
.
27.
Celis
JE
,
Gromova
I
,
Gromov
P
,
Moreira
JM
,
Cabezon
T
,
Friis
E
, et al
Molecular pathology of breast apocrine carcinomas: a protein expression signature specific for benign apocrine metaplasia
.
FEBS Lett
2006
;
580
:
2935
44
.
28.
Gromov
P
,
Espinoza
JA
,
Gromova
I
. 
Molecular and diagnostic features of apocrine breast lesions
.
Expert Rev Mol Diagn
2015
;
15
:
1011
22
.
29.
Wells
CA
,
El-Ayat
GA
. 
Non-operative breast pathology: apocrine lesions
.
J Clin Pathol
2007
;
60
:
1313
20
.
30.
Gallicchio
L
,
McSorley
MA
,
Newschaffer
CJ
,
Thuita
LW
,
Huang
HY
,
Hoffman
SC
, et al
Nonsteroidal antiinflammatory drugs, cyclooxygenase polymorphisms, and the risk of developing breast carcinoma among women with benign breast disease
.
Cancer
2006
;
106
:
1443
52
.
31.
Visscher
DW
,
Frost
MH
,
Hartmann
LC
,
Frank
RD
,
Vierkant
RA
,
McCullough
AE
, et al
Clinicopathologic features of breast cancers that develop in women with previous benign breast disease
.
Cancer
2016
;
122
:
378
85
.
32.
Milanese
TR
,
Hartmann
LC
,
Sellers
TA
,
Frost
MH
,
Vierkant
RA
,
Maloney
SD
, et al
Age-related lobular involution and risk of breast cancer
.
J Natl Cancer Inst
2006
;
98
:
1600
7
.
33.
Sickles
E
,
D'Orsi
C
,
Bassett
L
, et al
ACR BI-RADS Atlas® Mammography, ACR BI-RADS Atlas®, Breast Imaging Reporting and Data System
.
Reston, VA
:
American College of Radiology
; 
2003
.
34.
Fine
JP
,
Gray
RJ
. 
A proportional hazards model for the subdistribution of a competing risk
.
J Am Stat Assoc
1999
;
94
:
496
509
.
35.
Whitcomb
BW
,
Mumford
SL
,
Perkins
NJ
,
Wactawski-Wende
J
,
Bertone-Johnson
ER
,
Lynch
KE
, et al
Urinary cytokine and chemokine profiles across the menstrual cycle in healthy reproductive-aged women
.
Fertil Steril
2014
;
101
:
1383
91
.
36.
Nichols
HB
,
Schoemaker
MJ
,
Cai
J
,
Xu
J
,
Wright
LB
,
Brook
MN
, et al
Breast cancer risk after recent childbirth: a pooled analysis of 15 prospective studies
.
Ann Intern Med
2019
;
170
:
22
30
.
37.
Agurs-Collins
T
,
Ross
SA
,
Dunn
BK
. 
The many faces of obesity and its influence on breast cancer risk
.
Front Oncol
2019
;
9
:
765
.
38.
Bowers
LW
,
Maximo
IX
,
Brenner
AJ
,
Beeram
M
,
Hursting
SD
,
Price
RS
, et al
NSAID use reduces breast cancer recurrence in overweight and obese women: role of prostaglandin-aromatase interactions
.
Cancer Res
2014
;
74
:
4446
57
.
39.
Degnim
AC
,
Hoskin
TL
,
Arshad
M
,
Frost
MH
,
Winham
SJ
,
Brahmbhatt
RA
, et al
Alterations in the immune cell composition in premalignant breast tissue that precede breast cancer development
.
Clin Cancer Res
2017
;
23
:
3945
52
.
40.
Kerekes
D
,
Visscher
DW
,
Hoskin
TL
,
Radisky
DC
,
Brahmbhatt
RD
,
Pena
A
, et al
CD56+ immune cell infiltration and MICA are decreased in breast lobules with fibrocystic changes
.
Breast Cancer Res Treat
2018
;
167
:
649
58
.
41.
Madej
P
,
Franik
G
,
Kurpas
P
,
Owczarek
A
,
Chudek
J
,
Olszanecka-Glinianowicz
M
. 
Evaluation of adipokines, inflammatory markers, and sex hormones in simple and complex breast cysts' fluid
.
Dis Markers
2016
;
2016
:
5174929
.
42.
Bayraktar
S
,
Baghaki
S
,
Wu
J
,
Liu
DD
,
Gutierrez-Barrera
AM
,
Bevers
TB
, et al
Biomarker modulation study of celecoxib for chemoprevention in women at increased risk for breast cancer: a phase II pilot study
.
Cancer Prev Res
2020
;
13
:
795
802
.