Skip to main content
  • Research article
  • Open access
  • Published:

Parity, breastfeeding, and breast cancer risk by hormone receptor status and molecular phenotype: results from the Nurses’ Health Studies



Epidemiologic data suggest that parity increases risk of hormone receptor-negative breast cancer and that breastfeeding attenuates this association. Prospective data, particularly on the joint effects of higher parity and breastfeeding, are limited.


We investigated parity, breastfeeding, and breast cancer risk by hormone-receptor (estrogen (ER) and progesterone receptor (PR)) and molecular subtypes (luminal A, luminal B, HER2-enriched, and basal-like) in the Nurses’ Health Study (NHS; 1976–2012) and NHSII (1989–2013). A total of 12,452 (ER+ n = 8235; ER− n = 1978) breast cancers were diagnosed among 199,514 women. We used Cox proportional hazards models, adjusted for breast cancer risk factors, to calculate hazard ratios (HR) and 95% confidence intervals (CI).


Parous women had lower risk of ER+ breast cancer (vs. nulliparous, HR = 0.82 [0.77–0.88]); no association was observed for ER− disease (0.98 [0.84–1.13]; Phet = 0.03). Among parous women, breastfeeding was associated with lower risk of ER− (vs. never 0.82 [0.74–0.91]), but not ER+, disease (0.99 [0.94–1.05]; Phet < 0.001). Compared to nulliparous women, higher parity was inversely associated with luminal B breast cancer regardless of breastfeeding (≥ 3 children: ever breastfed, 0.78 [0.62–0.98]; never breastfed, 0.76 [0.58–1.00]) and luminal A disease only among women who had breastfed (≥ 3 children, 0.84 [0.71–0.99]). Basal-like breast cancer risk was suggestively higher among women with higher parity who never breastfed; associations were null among those who ever breastfed.


This study provides evidence that breastfeeding is inversely associated with hormone receptor-negative breast cancers, representing an accessible and cost-effective risk-reduction strategy for aggressive disease subtypes.


Data from recent prospective [1, 2], case-control [3,4,5,6], and case-only studies [7,8,9] suggest that parity increases risk of hormone receptor-negative breast cancer, while breastfeeding decreases risk (reviewed in [10]) or attenuates the detrimental effect of parity [3, 5, 11]. A general limitation of many prior prospective studies on these associations [1, 2, 12,13,14,15,16] is the inclusion of relatively few estrogen receptor (ER)-negative cases, given that only ~ 20% of tumors lack ER expression. Prospective data on the combination of higher parity and history of breastfeeding and ER− breast cancer risk are sparse [1, 11]. The African American Breast Cancer Epidemiology and Risk (AMBER) Consortium [11], including both case-control [4, 5] and cohort studies [1, 17], is the largest study to date (ER+, n = 2446; ER−, n = 1252). Within the consortium, parity was associated with increased risk of ER− breast cancer; this risk was attenuated among women who reported ever breastfeeding. Further, breastfeeding was only inversely associated with ER− disease; no association was observed for ER+ breast cancers.

To our knowledge, there are no prospective data on the combined effects of parity and breastfeeding by molecular phenotype (e.g., luminal subtypes, HER2-enriched, basal-like). Case-control [5, 7, 18] and case-only [9] studies suggest an increased risk of basal-like breast cancer only among parous women who did not breastfeed and an inverse association between breastfeeding and basal-like disease. In an earlier investigation of reproductive factors in the Nurses’ Health Studies (NHS, NHSII), we observed an inverse association between breastfeeding and basal-like breast cancer, but did not investigate the combined effect of parity and breastfeeding [13].

If breastfeeding is confirmed to reduce risk of hormone receptor-negative breast cancer, it represents an accessible and cost-effective risk-reduction strategy for this aggressive subtype. We hypothesized heterogeneous associations between parity and breastfeeding and ER+ and ER− breast cancer risk, with inverse associations hypothesized between parity and ER+ disease risk and between breastfeeding and ER− disease risk, and a positive association between parity and ER− breast cancer subtypes among women who did not breastfeed. Given intriguing, but inconsistent, prior results on the association between parity, breastfeeding, and breast cancer by subtype, and sparse prospective data on the combined effect of higher parity and breastfeeding, we examined these associations in the NHS and NHSII among almost 200,000 women.


Study population: the Nurses’ Health Study and Nurses’ Health Study II

The NHS was initiated in 1976 when 121,701 registered nurses, ages 30–55, completed and returned a mailed questionnaire [19]. The NHSII began in 1989 with 116,429 female registered nurses ages 25–42, using the same protocols. Participants have been followed biennially to update information on lifestyle and reproductive factors, and ascertain disease diagnoses. Cumulative follow-up rates are ≥ 94% in both cohorts.

Parity and breastfeeding assessment

Participants reported pregnancies and breastfeeding on the biennial questionnaire (pregnancies: NHS 1976–1984, 1996, and NHSII 1989–2009; breastfeeding NHS 1986 and NHSII 1993, 1997). Parity was investigated as nulliparous vs. parous and total number of pregnancies (1, 2, 3, 4+ children); breastfeeding was investigated as nulliparous/never vs. ever and total duration (< 6, 7–11, ≥ 12 months). Further, we cross-classified parity with breastfeeding (nulliparous, parous/ever breastfed, parous/never breastfed; and nulliparous, never breastfed and 1, 2, or ≥ 3 children, ever breastfed and 1, 2, or ≥ 3 children).

Covariate assessment

Age at each questionnaire was calculated using date of birth and questionnaire return date. Age at menarche and height were collected on baseline questionnaires. Other covariates included (date of collection) weight at age 18 (NHS 1980; NHSII 1989), height (NHS 1976; NHSII 1989), weight (biennially), oral contraceptive use (OC; NHS biennially until 1982; NHSII biennially), postmenopausal hormone therapy (HT; biennially), alcohol consumption (every 4 years, beginning in NHS 1980; NHSII 1991), age at first birth (NHS 1976 and biennially until 1982; NHS 1989 and biennially until 2011), menopausal status and age at menopause (biennially), diagnosis of benign breast disease (biennially), and family history of breast cancer (NHS 1982, and every 4 years beginning in 1988; NHSII 1989, and every 4 years beginning in 1997). Body mass index (BMI, kg/m2) at age 18 was calculated using self-reported weight and height.

Breast cancer case ascertainment

Participants reported disease status on the biennial questionnaires. Cases included in this analysis reported no prior cancer diagnosis (except non-melanoma skin) and were diagnosed with invasive breast cancer through June 2014 (NHS) or June 2013 (NHSII). A total of 12,452 (NHS n = 8807; NHSII n = 3645) incident breast cancer cases were identified; 82% had data on tumor ER status (ER+ n = 8235; ER− n = 1978). A study physician confirmed breast cancer cases and extracted invasive vs. in situ, hormone receptor status, and tumor characteristics from the medical record. Vital status was ascertained through June 2015, using the National Death Index. Of the 12,452 breast cancer cases identified, 11,369 (91%) were confirmed through medical record review; remaining cases were participant confirmed. Given the high confirmation rate by medical record for breast cancer in this cohort (99%), both medical record and participant-confirmed cases are included.

Tumor tissue collection

Details of tumor tissue collection have been detailed previously [20]. Tumor tissue samples were available from 4058 breast cancer cases (NHS, n = 3126; NHSII, n = 932). Tissue microarrays (TMAs) were created at the Dana Farber/Harvard Cancer Center (Boston, MA) with three 0.6-mm cores taken from each breast tumor included in the TMA block.

Immunohistochemical (IHC) analyses

IHC staining for ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2), epidermal growth factor receptor (EGFR), and cytokeratin (CK) 5/6 was performed using a Dako Autostainer (Dako Corporation, Carpinteria, CA). IHC methods have been described previously [21]. Immunostained TMA slides were evaluated by a study pathologist. Tumor sections with any nuclear staining for ER or PR were classified as ER+ or PR+; ER− or PR− had complete absence of staining. Tumors were considered HER2+ when more than 10% of cells showed moderate or strong staining. Tumors were classified as CK5/6+ or EGFR+ if there was evidence of cytoplasmic or membranous staining.

Classification of molecular subtypes

ER, PR, HER2, CK5/6, EGFR, and tumor grade were used to classify tumors into molecular subtypes; classification by IHC is a reasonable proxy for gene expression data [22,23,24]. Cases that were ER+ and/or PR+, HER2− and of histologic grade 1 and 2 were classified as luminal A cancers; cases that were either (a) ER+ and/or PR+ and HER2+ or (b) ER+ and/or PR+, HER2−, and of histologic grade 3 were classified as luminal B; cases that were ER−, PR−, and HER2+ were classified as HER2-enriched; and cases that were negative for ER, PR, and HER2 and positive for CK 5/6 and/or EGFR were categorized as basal-like. Cases that lacked expression of all five markers were considered “unclassified.” Of the invasive tumors on tissue microarrays, 3454 could be classified into the luminal A (n = 1903; 55%), luminal B (n = 934; 27%), HER2-enriched (n = 186; 5%), basal-like (n = 341; 10%), or unclassified (n = 90; 3%) subtypes. The few unclassified tumors were excluded from further analyses.

Statistical analysis

After exclusions for prior cancer (n = 4528) or missing data on parity, age at first birth, or breastfeeding (n = 34,088), 199,514 women remained in the analysis. We calculated person-years beginning at date of baseline questionnaire return and ending at the earliest of date of diagnosis of any cancer (except non-melanoma skin cancer), date of death, or end of follow-up (NHS: June 1, 2012; NHSII: June 1, 2013). We used multivariable Cox proportional hazards models to calculate hazard ratios (HR) and 95% confidence intervals (CI) for risk of breast cancer by hormone receptor subtype (ER+, ER−; ER+/PR+, ER+/PR−, ER−/PR−) and by molecular phenotype (luminal A, luminal B, HER2-enriched, basal-like). Models were stratified by age in months, follow-up year, and cohort. Covariates were updated whenever possible throughout follow-up. Covariates in the final models were height, BMI at age 18, weight change since age 18, age at menarche, age at natural menopause, menopausal status, HT use, alcohol consumption, total physical activity, family history of breast cancer, and history of benign breast disease. Models restricted to parous women were additionally adjusted for age at first birth; in sensitivity analyses, we adjusted for years between menarche and first birth and years between last birth and minimum of current age and age at menopause. In analyses restricted to parous women, parity and breastfeeding were mutually adjusted. We further evaluated associations for ER+ and ER− disease in subgroups defined by time between last pregnancy and diagnosis.

Heterogeneity (Phet) was assessed using a likelihood ratio test comparing models assuming the same association between the exposures and breast cancer subtypes to one assuming different associations for disease subtypes [25].

P values were considered statistically significant at < 0.05; all statistical tests were two-sided. Analyses were conducted in SAS 9.3 (Cary, NC).


A total of 93% of NHS participants and 67% of NHSII participants were parous at study baseline; among parous women, the majority reported ever breastfeeding (63% NHS, 79% NHSII). Median ages at recruitment for nulliparous women were 43.1 years in the NHS and 32.6 years in the NHSII; parous women were somewhat older (e.g., 2+ children, ever breastfed, median age 42.8 years in the NHS and 35.5 in the NHSII) (Table 1).

Table 1 Characteristics of the Nurses’ Health Study (NHS) and NHSII participants at baseline (NHS: 1976; NHSII: 1989)

The associations between parity, breastfeeding, and breast cancer risk differed by tumor ER status (Phet = 0.03). Parity was associated with lower risk of ER+ breast cancer (parous vs. nulliparous: HR = 0.82 [95% CI 0.77–0.88]) (Table 2), but not associated with ER− disease (HR = 0.98 [0.84–1.13]). In contrast, ever breastfeeding was only associated with reduced risk of ER− breast cancer (HR = 0.83 [0.75–0.92]), and not ER+ disease (HR = 0.99 [0.94–1.04]), relative to parous women who never breastfed and nulliparous women. Results were similar regardless of duration of breastfeeding (e.g., ≥ 12 months vs. never, ER− HR = 0.83 [0.73–0.94], Ptrend = 0.05; ER+ HR = 1.02 [0.96–1.09], Ptrend = 0.22). The association between parity and ER+ tumors did not differ by breastfeeding history (e.g., parous, never breastfed: HR = 0.83 [0.77–0.90]; parous, ever breastfed: HR = 0.82 [0.76–0.88]), and parity was not associated with ER− breast cancer risk in either breastfeeding history subgroup (parous, never breastfed: HR = 1.11 [0.94–1.31]; parous, ever breastfed: HR = 0.92 [0.79–1.08]).

Table 2 Multivariable-adjusteda hazard ratios (HRs) and 95% confidence intervals for breast cancer overall and by ER status in relation to parity and breastfeeding: Nurses’ Health Study (1976–2012) and Nurses’ Health Study II (1989–2013)

We next restricted the analyses to parous women and evaluated number of children, breastfeeding, and breast cancer risk (Table 2). Relative to women with parity of 1, only parity of 4 or higher was associated with lower ER+ breast cancer risk (≥ 4 children HR = 0.91 [0.82–1.00]). No association was observed in the other parity categories (e.g., 3 children HR = 1.03 [0.94–1.13]), or in any category for ER− breast cancer. The association between breastfeeding and breast cancer risk among parous women was similar to that observed in the analysis among all participants, including nulliparous women (e.g., ever vs. never, ER+ HR = 0.99 [0.94–1.05]; ER− HR = 0.82 [0.74–0.91]). Duration of breastfeeding was associated with a trend toward higher risk of ER+ breast cancer (Ptrend among ever breastfed = 0.02), but no significant associations were observed in categories of duration (e.g., ≥ 12 months vs. never, HR = 1.04 [0.97–1.11]). The associations between breastfeeding and ER− breast cancer risk were similar regardless of breastfeeding duration (e.g., ≥ 12 months vs. never, HR = 0.80 [0.70–0.92]; < 6 months vs. never, HR = 0.84 [0.75–0.95]; Ptrend among ever breastfed = 0.35). Patterns were similar for ER+/PR+ and ER−/PR− breast cancer (Additional file 1: Table S1). Statistical adjustment for years between menarche and age at first birth and years between age at last birth and minimum of current age and age at menopause did not impact the results (data not shown).

We next considered the joint association of parity and breastfeeding and breast cancer risk by molecular phenotypes (Table 3). For luminal A breast cancer, parity of 3 or higher, relative to nulliparity, was associated with significantly lower risk only among women who breastfed (≥ 3 children, ever breastfed: HR = 0.84 [0.71–0.99]; never breastfed: HR = 0.93 [0.77–1.12]). In contrast, for luminal B breast cancer, the association between parity and breast cancer risk was similar regardless of breastfeeding history (e.g., ≥ 3 children, ever breastfed: HR = 0.78 [0.62–0.98]; never breastfed: HR = 0.76 [0.58–1.00]). We observed the suggestion of increased risk of basal-like tumors among parous women who never breastfed (e.g., ≥ 3 children, never breastfed: HR = 1.43 [0.92–2.23]); this increased risk was not observed among those who ever breastfed (HR = 1.05 [0.70–1.57]). Risk of HER2-enriched disease was higher in parous, as compared to nulliparous, women, regardless of breastfeeding (parous/never breastfed, HR = 2.09 [1.02–4.30]; parous/ever breastfed, HR = 1.90 [0.95–3.80]).

Table 3 Multivariable-adjusteda hazard ratios (HRs) and 95% confidence intervals for molecular breast cancer subtypes in relation to parity and breastfeeding: Nurses’ Health Study (1976–2012) and Nurses’ Health Study II (1989–2013)

Overall, results were similar in the individual cohorts (NHS and NHSII). We observed significant heterogeneity between cohorts for ER+ disease (p = 0.02) in analyses considering the cross-classification of parity and breastfeeding. For ER+ disease, parity ≥ 3 among women who never breastfed was more strongly associated with breast cancer risk in the NHS (HR = 0.75 [0.67–0.84]) than in the NHSII (HR = 0.85 [0.67–1.08]); associations were similar in other categories (e.g., ≥3 children ever breastfed, NHS HR = 0.76 [0.69–0.85]; NHSII HR = 0.79 [0.69–0.89]).

Finally, we evaluated parity, breastfeeding, and breast cancer risk in strata of time between last pregnancy and diagnosis (< 10 vs. ≥ 10 years). Inverse associations between parity and ER+ disease and breastfeeding and ER− disease were observed only in women diagnosed 10+ years after last pregnancy (Additional file 1: Table S2). No associations were observed among participants diagnosed within 10 years of last pregnancy.


This large, prospective study provides additional evidence of lower risk of ER− breast cancer among women who have breastfed and, to our knowledge, is the first prospective investigation on the joint effect of parity and breastfeeding on breast cancer risk by molecular phenotype. While breastfeeding was inversely associated with ER− breast cancer, no association was observed for the ER+ subtype. In contrast, parity was associated with lower risk of ER+ breast cancer, but not associated with ER− disease. In analyses by molecular phenotype, higher parity was inversely associated with the luminal subtypes. Parity was associated with higher risk of basal-like disease among women who did not breastfeed and HER2-enriched breast cancer risk regardless of breastfeeding history. These results provide confirmatory evidence of heterogeneity in the associations between parity and breastfeeding and breast cancer subtypes and support a risk-reducing role of breastfeeding for ER- breast cancer.

The breast undergoes proliferation and differentiation during pregnancy, in preparation for lactation. Terminal differentiation of the terminal ductal lobular unit (TDLU) occurs in the final trimester; this is hypothesized to be an important mechanism linking full-term pregnancy to long-term reduced risk of breast cancer, shielding the breast tissue from carcinogenic transformation [26]. Recent research has identified parity-related gene expression signatures in normal breast tissue [27,28,29], providing further evidence for the long-term impact of pregnancy on the breast. A “parity signature” including differential expression in inflammatory and immune-related pathways may be retained in ER+, but not ER−, tumors [27] supporting a differential role for parity in the etiology of ER+ vs. ER− breast cancer. However, data are limited and the parity-related gene signatures are not entirely consistent across studies.

Direct data describing underlying physiologic mechanisms between breastfeeding and breast cancer subtypes are limited [30]. Post-pregnancy, and at cessation of lactation, the breast undergoes involution. Involution has features of wound healing [31, 32], which has established parallels to tumorigenic microenvironments, including involvement of cytokines, growth factors, and matrix metalloproteases (MMPs) (reviewed in [33]). Breastfeeding delays the process of involution, and it has been hypothesized that the more gradual involution associated with lactation (i.e., from exclusive breastfeeding through weaning) may reduce longer-term breast cancer risk. Further, the postpartum hormonal milieu differs in breastfeeding vs. non-breastfeeding women (e.g., breastfeeding is associated with lower estradiol and higher prolactin levels for up to 12 weeks postpartum and has been associated with a slower decline in postpartum osteoprotegerin levels) [34, 35]. This postpartum hormonal milieu may differentially impact risk in the context of involution in the immediate postpartum vs. after breastfeeding and may differentially impact development of ER+ vs. ER− disease. Finally, murine models show a lower percentage of ER+ luminal cells in the mammary glands of parous, relative to virgin, mice with no change in the number of basal-like progenitor cells, though these basal-like progenitor cells demonstrated characteristics suggesting reduced carcinogenic potential in parous mice [36]. Nonetheless, involution proximate to delivery (i.e., in the absence of breastfeeding) may potentiate basal-like progenitors toward a carcinogenic pathway.

While some prior research is suggestive of a positive association between parity and ER− breast cancer, data to date are not consistent (reviewed in [37]). The AMBER consortium is the largest study to date [11]. Parity was associated with 33% increased risk of ER− disease (n = 1252; parous vs. nulliparous, RR = 1.33 [1.11–1.59]). Higher risk was observed with higher parity among women who did not breastfeed (e.g., ≥ 4 vs. 1 birth, RR = 1.68 [1.15–2.44]) but not among women reporting ever breastfeeding (RR = 1.33 [0.91–1.95]). Ever parous, or number of children, were not significantly associated with ER− breast cancer in the current study, regardless of breastfeeding. Importantly, the AMBER consortium includes women of African-American descent, whereas the NHS/NHSII study population is predominantly non-Hispanic white. ER−/PR− (and triple negative) breast cancer is more common among African American women, as are higher parity [38] and lower rates of breastfeeding [39]. Consistent with most prior prospective studies, ever breastfeeding was significantly associated with lower risk of ER− disease, and not associated with ER+ disease, in the current study. These results are in line with a recent meta-analysis of prospective studies (ever vs. never, summary RRs: ER−/PR−, 0.84 [0.72–0.97]; ER+/PR+, 1.00 [0.90–1.10]) [10]. In analyses evaluating duration of breastfeeding, similar associations were observed between breastfeeding and ER− breast cancer regardless of duration, suggesting even shorter-term breastfeeding (cumulative < 6 months) may decrease ER− breast cancer risk. In contrast, a trend toward higher risk was observed between breastfeeding duration and ER+ breast cancer, though no significant associations were observed in any of the categories of duration. To our knowledge, this has not previously been reported. Previous prospective studies have not observed significant associations between breastfeeding duration and breast cancer risk by subtype [2, 12, 14,15,16]. Additional studies with detailed data on breastfeeding duration are warranted.

Few prior studies have examined the joint effects of parity and breastfeeding and breast cancer risk by molecular phenotype [5, 7, 9, 18]. Millikan et al., in the Carolina Breast Cancer Study [5], observed higher risk of basal-like breast cancer among women with high parity in the absence of breastfeeding (parity 3+, never breastfed: OR = 1.9 [1.1–3.3]); breastfeeding attenuated the increased risk (OR = 1.3 [0.7–2.3]). A case-only investigation reported similar associations [9]. In our prior analysis by molecular phenotype in the NHS/NHSII, we observed a positive association between parity and HER2-enriched disease, but no association with the other subtypes [13]. Breastfeeding was associated with a reduced risk of basal-like breast cancer (e.g., 7+ months vs. never, among parous women, HR = 0.65 [0.49–0.87]); however, the joint effects of parity and breastfeeding were not considered. In the current investigation, consistent with prior studies, parity was suggestively positively associated with basal-like breast cancer, but only among women with high parity and who never breastfed. The positive association of parity with HER2 disease did not differ by breastfeeding.

We provide the largest prospective study to date on parity, breastfeeding, and breast cancer by subtype, and add to the limited literature on these associations by molecular phenotype. While we expect that parity was accurately reported, there may be misclassification of breastfeeding duration. We would expect any misclassification to be nondifferential. Tumor tissue to define molecular phenotype was available for 33% of the cases, resulting in small case numbers in some subgroups.


Our results support a role for breastfeeding in reducing risk of hormone receptor-negative breast cancer, adding to the established benefits of breastfeeding. While further research is needed to describe the mechanistic pathway resulting in this decreased risk, our investigation provides additional compelling evidence that breastfeeding is a modifiable risk factor for the breast cancer subtypes with the fewest targeted therapies and the least favorable outcomes.



African American Breast Cancer Epidemiology and Risk (AMBER) Consortium


Body mass index


Confidence interval




Epidermal growth factor receptor


Estrogen receptor


Human epidermal growth factor receptor 2


Hazard ratio


Hormone therapy


Matrix metalloproteases


Nurses’ Health Study


Oral contraceptive


Progesterone receptor


Terminal ductal lobular unit


Tissue microarray


  1. 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 Biomarkers Prev. 2011;20(9):1883–91.

    Article  CAS  Google Scholar 

  2. Phipps AI, Chlebowski RT, Prentice R, McTiernan A, Wactawski-Wende J, Kuller LH, Adams-Campbell LL, Lane D, Stefanick ML, Vitolins M, et al. Reproductive history and oral contraceptive use in relation to risk of triple-negative breast cancer. J Natl Cancer Inst. 2011;103(6):470–7.

    Article  Google Scholar 

  3. Work ME, John EM, Andrulis IL, Knight JA, Liao Y, Mulligan AM, Southey MC, Giles GG, Dite GS, Apicella C, et al. Reproductive risk factors and oestrogen/progesterone receptor-negative breast cancer in the breast Cancer family registry. Br J Cancer. 2014;110(5):1367–77.

    Article  CAS  Google Scholar 

  4. Ambrosone CB, Zirpoli G, Ruszczyk M, Shankar J, Hong CC, McIlwain D, Roberts M, Yao S, McCann SE, Ciupak G, et al. Parity and breastfeeding among African-American women: differential effects on breast cancer risk by estrogen receptor status in the Women's circle of health study. Cancer Causes Control. 2014;25(2):259–65.

    Article  Google Scholar 

  5. Millikan RC, Newman B, Tse CK, Moorman PG, Conway K, Dressler LG, Smith LV, Labbok MH, Geradts J, Bensen JT, et al. Epidemiology of basal-like breast cancer. Breast Cancer Res Treat. 2008;109(1):123–39.

    Article  Google Scholar 

  6. Lord SJ, Bernstein L, Johnson KA, Malone KE, McDonald JA, Marchbanks PA, Simon MS, Strom BL, Press MF, Folger SG, et al. Breast cancer risk and hormone receptor status in older women by parity, age of first birth, and breastfeeding: a case-control study. Cancer Epidemiol Biomarkers Prev. 2008;17(7):1723–30.

    Article  CAS  Google Scholar 

  7. Redondo CM, Gago-Dominguez M, Ponte SM, Castelo ME, Jiang X, Garcia AA, Fernandez MP, Tome MA, Fraga M, Gude F, et al. Breast feeding, parity and breast cancer subtypes in a Spanish cohort. PLoS One. 2012;7(7):e40543.

    Article  CAS  Google Scholar 

  8. Shinde SS, Forman MR, Kuerer HM, Yan K, Peintinger F, Hunt KK, Hortobagyi GN, Pusztai L, Symmans WF. Higher parity and shorter breastfeeding duration: association with triple-negative phenotype of breast cancer. Cancer. 2010;116(21):4933–43.

    Article  Google Scholar 

  9. Kwan ML, Kushi LH, Weltzien E, Maring B, Kutner SE, Fulton RS, Lee MM, Ambrosone CB, Caan BJ. Epidemiology of breast cancer subtypes in two prospective cohort studies of breast cancer survivors. Breast Cancer Res. 2009;11(3):R31.

    Article  Google Scholar 

  10. Islami F, Liu Y, Jemal A, Zhou J, Weiderpass E, Colditz G, Boffetta P, Weiss M. Breastfeeding and breast cancer risk by receptor status--a systematic review and meta-analysis. Ann Oncol. 2015;26(12):2398–407.

    CAS  PubMed  PubMed Central  Google Scholar 

  11. Palmer JR, Viscidi E, Troester MA, Hong CC, Schedin P, Bethea TN, Bandera EV, Borges V, McKinnon C, Haiman CA, et al. Parity, lactation, and breast cancer subtypes in African American women: results from the AMBER consortium. J Natl Cancer Inst. 2014;106:10.

    Article  Google Scholar 

  12. Ritte R, Tikk K, Lukanova A, Tjonneland A, Olsen A, Overvad K, Dossus L, Fournier A, Clavel-Chapelon F, Grote V, et al. Reproductive factors and risk of hormone receptor positive and negative breast cancer: a cohort study. BMC Cancer. 2013;13:584.

    Article  Google Scholar 

  13. Sisti JS, Collins LC, Beck AH, Tamimi RM, Rosner BA, Eliassen AH. Reproductive risk factors in relation to molecular subtypes of breast cancer: results from the nurses’ health studies. Int J Cancer. 2016;138(10):2346–56.

    Article  CAS  Google Scholar 

  14. Ellingjord-Dale M, Vos L, Tretli S, Hofvind S, Dos-Santos-Silva I, Ursin G. Parity, hormones and breast cancer subtypes - results from a large nested case-control study in a national screening program. Breast Cancer Res. 2017;19(1):10.

    Article  Google Scholar 

  15. Horn J, Opdahl S, Engstrom MJ, Romundstad PR, Tretli S, Haugen OA, Bofin AM, Vatten LJ, Asvold BO. Reproductive history and the risk of molecular breast cancer subtypes in a prospective study of Norwegian women. Cancer Causes Control. 2014;25(7):881–9.

    Article  Google Scholar 

  16. Ma H, Henderson KD, Sullivan-Halley J, Duan L, Marshall SF, Ursin G, Horn-Ross PL, Largent J, Deapen DM, Lacey JV Jr, et al. Pregnancy-related factors and the risk of breast carcinoma in situ and invasive breast cancer among postmenopausal women in the California teachers study cohort. Breast Cancer Res. 2010;12(3):R35.

    Article  Google Scholar 

  17. Setiawan VW, Monroe KR, Wilkens LR, Kolonel LN, Pike MC, Henderson BE. Breast cancer risk factors defined by estrogen and progesterone receptor status: the multiethnic cohort study. Am J Epidemiol. 2009;169(10):1251–9.

    Article  Google Scholar 

  18. 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.

    Article  Google Scholar 

  19. Colditz GA, Hankinson SE. The Nurses’ health study: lifestyle and health among women. Nat Rev Cancer. 2005;5(5):388–96.

    Article  CAS  Google Scholar 

  20. Collins LC, Marotti JD, Baer HJ, Tamimi RM. Comparison of estrogen receptor results from pathology reports with results from central laboratory testing. J Natl Cancer Inst. 2008;100(3):218–21.

    Article  CAS  Google Scholar 

  21. Tamimi RM, Baer HJ, Marotti J, Galan M, Galaburda L, Fu Y, Deitz AC, Connolly JL, Schnitt SJ, Colditz GA, et al. Comparison of molecular phenotypes of ductal carcinoma in situ and invasive breast cancer. Breast Cancer Res. 2008;10(4):R67.

    Article  Google Scholar 

  22. Brenton JD, Carey LA, Ahmed AA, Caldas C. Molecular classification and molecular forecasting of breast cancer: ready for clinical application? J Clin Oncol. 2005;23(29):7350–60.

    Article  CAS  Google Scholar 

  23. Nielsen TO, Hsu FD, Jensen K, Cheang M, Karaca G, Hu Z, Hernandez-Boussard T, Livasy C, Cowan D, Dressler L, et al. Immunohistochemical and clinical characterization of the basal-like subtype of invasive breast carcinoma. Clin Cancer Res. 2004;10(16):5367–74.

    Article  Google Scholar 

  24. Abd El-Rehim DM, Ball G, Pinder SE, Rakha E, Paish C, Robertson JF, Macmillan D, Blamey RW, Ellis IO. High-throughput protein expression analysis using tissue microarray technology of a large well-characterised series identifies biologically distinct classes of breast cancer confirming recent cDNA expression analyses. Int J Cancer. 2005;116(3):340–50.

    Article  CAS  Google Scholar 

  25. Wang M, Spiegelman D, Kuchiba A, Lochhead P, Kim S, Chan AT, Poole EM, Tamimi R, Tworoger SS, Giovannucci E, et al. Statistical methods for studying disease subtype heterogeneity. Stat Med. 2016;35(5):782–800.

    Article  Google Scholar 

  26. Russo J, Mailo D, Hu YF, Balogh G, Sheriff F, Russo IH. Breast differentiation and its implication in cancer prevention. Clinical Cancer Res. 2005;11(2 Pt 2):931s–6s.

    CAS  Google Scholar 

  27. Rotunno M, Sun X, Figueroa J, Sherman ME, Garcia-Closas M, Meltzer P, Williams T, Schneider SS, Jerry DJ, Yang XR, et al. Parity-related molecular signatures and breast cancer subtypes by estrogen receptor status. Breast Cancer Res. 2014;16(4):R74.

    Article  Google Scholar 

  28. Asztalos S, Gann PH, Hayes MK, Nonn L, Beam CA, Dai Y, Wiley EL, Tonetti DA. Gene expression patterns in the human breast after pregnancy. Cancer Prev Res (Phila). 2010;3(3):301–11.

    Article  CAS  Google Scholar 

  29. Belitskaya-Levy I, Zeleniuch-Jacquotte A, Russo J, Russo IH, Bordas P, Ahman J, Afanasyeva Y, Johansson R, Lenner P, Li X, et al. Characterization of a genomic signature of pregnancy identified in the breast. Cancer Prev Res (Phila). 2011;4(9):1457–64.

    Article  Google Scholar 

  30. Faupel-Badger JM, Arcaro KF, Balkam JJ, Eliassen AH, Hassiotou F, Lebrilla CB, Michels KB, Palmer JR, Schedin P, Stuebe AM, et al. Postpartum remodeling, lactation, and breast cancer risk: summary of a National Cancer Institute-sponsored workshop. J Natl Cancer Inst. 2013;105(3):166–74.

    Article  Google Scholar 

  31. 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(8):1803–13.

    Article  CAS  Google Scholar 

  32. O'Brien J, Hansen K, Barkan D, Green J, Schedin P, O'Brien J, Hansen K, Barkan D, Green J, Schedin P. Non-steroidal anti-inflammatory drugs target the pro-tumorigenic extracellular matrix of the postpartum mammary gland. Int J Dev Biol. 2011;55(7–9):745–55.

    Article  Google Scholar 

  33. Arwert EN, Hoste E, Watt FM. Epithelial stem cells, wound healing and cancer. Nat Rev Cancer. 2012;12(3):170–80.

    Article  CAS  Google Scholar 

  34. Naylor KE, Rogers A, Fraser RB, Hall V, Eastell R, Blumsohn A. Serum osteoprotegerin as a determinant of bone metabolism in a longitudinal study of human pregnancy and lactation. J Clin Endocrinol Metab. 2003;88(11):5361–5.

    Article  CAS  Google Scholar 

  35. Yen SS, Jaffe RB, Barbieri RL. Reproductive endocrinology. 4th ed. Philadelphia: W.B. Saunders Company; 1999.

    Google Scholar 

  36. Meier-Abt F, Milani E, Roloff T, Brinkhaus H, Duss S, Meyer DS, Klebba I, Balwierz PJ, van Nimwegen E, Bentires-Alj M: Parity induces differentiation and reduces Wnt/ notch signaling ratio and proliferation potential of basal stem/progenitor cells isolated from mouse mammary epithelium. Breast Cancer Res 2013, 15(2):R36.

  37. Anderson KN, Schwab RB, Martinez ME. Reproductive risk factors and breast cancer subtypes: a review of the literature. Breast Cancer Res Treat. 2014;144(1):1–10.

    Article  Google Scholar 

  38. Statistics NCfH: National Vital Statistics Reports Volumne 64, number 1, January 15, 2015. 2015:1–68.

  39. McDowell MM, Wang CY, Kennedy-Stephenson J: Breastfeeding in the United States: findings from the national health and nutrition examination surveys, 1999-2006. NCHS Data Brief 2008(5):1–8.

Download references


We would like to thank the participants and staff of the Nurses’ Health Study for their valuable contributions as well as the following state cancer registries for their help: AL, AZ, AR, CA, CO, CT, DE, FL, GA, ID, IL, IN, IA, KY, LA, ME, MD, MA, MI, NE, NH, NJ, NY, NC, ND, OH, OK, OR, PA, RI, SC, TN, TX, VA, WA, WY. The authors assume full responsibility for analyses and interpretation of these data.


This work was funded by National Institute of Health grants UM1 CA186107, UM1 CA176726, and R01 CA050385. RT Fortner was supported in part by T32 CA009001.

Availability of data and materials

The datasets analysed during the current study are available by application by following the instructions here:

Author information

Authors and Affiliations



RTF, SEH, RMT, and AHE were responsible for the conception and design. LCC, SEH, RMT, and AHE were responsible for the acquisition of data. RTF, JS, BC, BR, and AHE were responsible for the statistical analyses and interpretation. RTF wrote the first draft of the manuscript; all authors revised the manuscript for important intellectual content. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Renée T. Fortner.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Institutional Review Board of the Brigham and Women’s Hospital (Boston, MA). Participants provided informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1:

Table S1. Multivariate-adjusteda hazard ratios (HRs) and 95% confidence intervals for breast cancer by ER/PR status in relation to parity and lactation: Nurses’ Health Study (1976–2012) and Nurses’ Health Study II (1989–2013). Table S2. Multivariate-adjusteda hazard ratios (HRs) and 95% confidence intervals for breast cancer by ER status in relation to parity and lactation, stratified by time since last birth: Nurses’ Health Study (1976–2012) and Nurses’ Health Study II (1989–2013) (DOCX 30 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Fortner, R.T., Sisti, J., Chai, B. et al. Parity, breastfeeding, and breast cancer risk by hormone receptor status and molecular phenotype: results from the Nurses’ Health Studies. Breast Cancer Res 21, 40 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: