- Research article
- Open Access
Reproductive factors and risk of hormone receptor positive and negative breast cancer: a cohort study
BMC Cancer volume 13, Article number: 584 (2013)
The association of reproductive factors with hormone receptor (HR)-negative breast tumors remains uncertain.
Within the EPIC cohort, Cox proportional hazards models were used to describe the relationships of reproductive factors (menarcheal age, time between menarche and first pregnancy, parity, number of children, age at first and last pregnancies, time since last full-term childbirth, breastfeeding, age at menopause, ever having an abortion and use of oral contraceptives [OC]) with risk of ER-PR- (n = 998) and ER+PR+ (n = 3,567) breast tumors.
A later first full-term childbirth was associated with increased risk of ER+PR+ tumors but not with risk of ER-PR- tumors (≥35 vs. ≤19 years HR: 1.47 [95% CI 1.15-1.88] p trend < 0.001 for ER+PR+ tumors; ≥35 vs. ≤19 years HR: 0.93 [95% CI 0.53-1.65] p trend = 0.96 for ER-PR- tumors; P het = 0.03). The risk associations of menarcheal age, and time period between menarche and first full-term childbirth with ER-PR-tumors were in the similar direction with risk of ER+PR+ tumors (phet = 0.50), although weaker in magnitude and statistically only borderline significant. Other parity related factors such as ever a full-term birth, number of births, age- and time since last birth were associated only with ER+PR+ malignancies, however no statistical heterogeneity between breast cancer subtypes was observed. Breastfeeding and OC use were generally not associated with breast cancer subtype risk.
Our study provides possible evidence that age at menarche, and time between menarche and first full-term childbirth may be associated with the etiology of both HR-negative and HR-positive malignancies, although the associations with HR-negative breast cancer were only borderline significant.
Breast cancer is a complex and heterogeneous disease with a variety of histo-pathological and molecular subtypes with diverse clinical outcomes and relationships with established risk factors [1–3]. The major sub-classification of clinical breast tumors is based on the detection of estrogen (ER) and progesterone (PR) receptors and guides targeted therapies and provides important prognostic information . The presence or absence of hormone receptors, along with human epidermal growth factor-2 (HER2) also broadly correspond to more detailed molecular subclassification of breast tumors, as determined by microarray-based gene expression profiling coupled to hierarchical clustering analyses [5–7]. In addition to the clinical use of ER and PR, epidemiological data indicate that the association of reproductive history with breast cancer differs by the expression of ER and PR receptors .
Factors that influence the lifetime cumulative exposure to hormones during reproductive life, such as the age at menarche, age at first child-birth, time between age at menarche and first child birth, number of children, use of oral contraceptives (OC) and breastfeeding, have been suggested to be associated with risk of hormone receptor (HR)-positive malignancies (ER-positive or joint ER+PR+)[2, 8–11]. However, distinct risk factors for HR-negative (ER-negative, or joint ER-PR) cancer are debated [2, 3, 8, 12] and the etiologies of ER+PR+ and ER-PR- tumors remain unclear.
The incidence of HR-negative disease drops remarkably after menopause  suggesting that ovarian derived sex steroid hormones do have an impact on HR-negative tumors and recent studies are starting to show risk associations of reproductive factors with HR-negative malignancies [2, 3]. In fact, opposite risk associations between ER-PR- and ER+PR+ tumors have been observed with parity [11, 14], age at first pregnancy  and breastfeeding. Nonetheless, due to the rarity and heterogeneous nature of HR-negative breast tumors, epidemiological studies have been hindered by small sample sizes resulting in inconsistent risk associations between reproductive factors and HR-negative disease [8, 11, 15, 16].
The incidence of HR-negative disease drops remarkably after menopause  suggesting that ovarian derived sex steroid hormones do have an impact on HR-negative tumors and recent studies are starting to show risk associations of reproductive factors with HR-negative malignancies [2, 3]. In fact, opposite risk associations between ER-PR- and ER+PR+ tumors have been observed with parity [11, 14], age at first pregnancy  and breastfeeding . Nonetheless, due to the rarity and heterogeneous nature of HR-negative breast tumors, epidemiological studies have been hindered by small sample sizes resulting in inconsistent risk associations between reproductive factors and HR-negative disease [3, 8, 14, 15].
The European Prospective Investigation into Cancer and Nutrition (EPIC) is a multi-center prospective cohort study designed to investigate the relationships between diet, nutrition and metabolic factors and cancer, consisting of approximately 360,000 women and 150,000 men aged mostly between 25–70 years [16, 17]. All participants were enrolled between 1992 and 2000 and came from 23 regional and national research centers located in 10 western European countries: Denmark, France, Italy, Germany, Greece Norway, Spain, Sweden, The Netherlands and the United Kingdom. Extensive details about the standardized procedures for recruitment, measuring baseline anthropometry, questionnaires on current habitual diet, reproductive and menstrual history, exogenous hormone use [OC and hormone replacement therapy (HRT) use], medical history, lifetime smoking and alcohol consumption history, occupation, level of education and physical activity and biological sample collection at study centers are given elsewhere [16, 17]. All subjects gave written informed consent. The Internal Review Boards of the International Agency for Research on Cancer and the local ethics committees in participating countries approved the analyses based on EPIC participants.
Of the approximately 360,000 female participants in EPIC, women were excluded a priori if they had a history of cancer prior to recruitment or were missing a diagnosis or censoring date, thus leaving 345,153 participants. At the time of this analysis, three EPIC study centers, (Granada, Murcia and Malmo), did not provide any information on breast tumor hormone receptor status and therefore were excluded from this analysis (n = 26,091). Women were further excluded if they were missing questionnaire data (n = 526) or were missing data on age at menarche, age at menopause, age at first full-term birth, ever use of OCs, number of full-term births, age at last full-term childbirth and duration of breastfeeding (n = 7,439). This left a total of 311,097 women with 9,456 first primary invasive breast cancer cases from 10 countries for the present analysis.
Questionnaire data and classification of reproductive variables
The details of standardized procedures for collecting baseline information on the age at first and last menstruation, parity, breastfeeding, exogenous hormone use, and hysterectomy from the general lifestyle questionnaire has been previously reported [17, 18]. Briefly, in Greece, Italy, the Netherlands, Sweden and the United Kingdom, age at menarche was asked in years. In the other countries, age at menarche was asked in defined categories (≤8, 9, 10,…, 18, 19 or >19 years). The number of full-term pregnancies (defined as the sum of all live and stillborn children born) and spontaneous or induced abortions were also collected at baseline, together with the ages of the first three and last deliveries and the ages at first and last induced or spontaneous abortions and stillbirths. Except for Norway and the Swedish center Umeå, where information about multiple pregnancies was available, the number of pregnancies is overestimated as multiple pregnancies were counted as different pregnancies. The length of time between menarche and age at first pregnancy was estimated among women who had menarche between the ages of 8 and 20 years (time between menarche and first full-term birth = age at first full-term birth – age at menarche).
Women were considered postmenopausal at recruitment if they had had no menstrual cycles in the last 12 months, were older than 55 years (if the menstrual cycle history was missing), or had a bilateral oophorectomy. Women who were aged 46–55 years and had incomplete or were missing questionnaire data on menstrual history were classified with a peri/-or of unknown menopausal status. Women were deemed premenopausal if they reported regular menstrual cycles in the last 12 months or if they were younger than 46 years of age (if the menstrual cycle history was missing).
The details of standardized procedures for measuring height and weight at EPIC study centers has also been previously reported . In most countries, height, weight and waist and hip circumferences were measured to the nearest centimeter and kilogram, in light clothing, according to standardized protocols. In Norway, Umeå and a large proportion from France, subjects’ height and weight were measured and self-reported by the cohort participants themselves, following detailed instructions [17, 19]. For subjects that had neither self-reported nor measured weight or height data, the center-, age- and gender-specific average weight and height values were imputed for anthropometry variables used for adjustment purposes only. A sensitivity analysis that restricted the adjusted variables to those without imputation showed similar results to those presented (data not shown).
Prospective ascertainment of breast cancer cases and the coding of receptor status
In all countries (except for France, Germany and Greece) incident breast cancer cases were identified using record linkage with cancer and pathology registries. In France, Germany and Greece, cancer occurrence was prospectively ascertained through linkage with health insurance records and regular direct contact with participants and their next of kin, and all reported breast cancer cases were then systematically verified against clinical and pathological records. Mortality data were coded according to the 10th Revision of the International Statistical Classification of Diseases, Injuries, and Causes of Death (ICD-10), and cancer incidence data were coded according to the International Classification of Diseases for Oncology (ICD-O-2). Invasive (primary, malignant) breast cancer cases were classified as per the International Classification of Diseases for Oncology (Topography C50), second revision (ICD-O-2). Breast tumor receptor status was standardized across EPIC centers using the following criteria for a positive expression: ≥10% cells stained, any ‘plus-system’ description, ≥20 fmol/mg, an Allred score of ≥3, an IRS ≥2, or an H-score ≥10 .
Vital status was collected from regional or national mortality registries. The last updates of endpoint data for cancer incidence and vital status were between 2005 and 2010, depending on the center. Women were considered at risk from the time of recruitment until breast cancer diagnosis or censoring (age at death, loss to follow up, end of follow up, or diagnosis of other cancer) respectively. A total of 7,095 breast cancer cases had information on ER status (5,723 ER-positive, 1,372 ER-negative); of which, 5,843 had further information on PR status (3,567 ER+PR+, 1,078 ER+PR-, 200 ER-PR+, 998 ER-PR-).
Associations between reproductive factors and the risk of breast cancer subtype were evaluated using Cox proportional hazards models to estimate hazard ratios (HR) and 95% confidence intervals (CIs). Breast cancer subtypes were defined as jointly classified ER+PR+ or ER-PR- breast tumors. Results for ER-positive versus ER-negative (ignoring PR status); and PR-positive versus PR-negative (ignoring ER status) were generally similar to the jointly defined ERPR breast cancer subtypes and have been included in Additional file 1: Table S1. Results for breast tumors with discordant ER and PR status and unknown ER and/or PR status have been reported in Additional file 2: Table S2. All analyses were stratified by age at recruitment in one-year categories and by study center, to prevent violations of the proportional-hazard assumption. Trend tests across levels of exposure categories were performed on the continuous categorical variables entered as ordered, quantitative variables into the models.
Age at menarche was categorized as ≤12, 13–14 and ≥15 years and time between menarche and first full-term childbirth as <10 and ≥10 years. Parity related variables were divided into the following categories ever vs. never, number of full-term pregnancies (1, 2, ≥3), age at first full-term childbirth as ≤19, 20–24, 25-29,30-34, ≥35 years, age at last full-term childbirth since recruitment as ≤24, 25–29, 30-34, ≥35 years and time since last child birth as ≤20 and >20 years. Breastfeeding was categorized as ever versus never, and ≤1 month, 2-3, 4–6, 7–12, 13–17 and ≥18 months for total cumulative duration of breastfeeding. Dichotomized categories of ever vs. never having had a spontaneous or induced abortion, ever vs. never OC use, and current versus not currently using OCs (at baseline) also were analyzed. The duration of OC use was categorized into ≤1, 2–4, 5–9, and ≥10 years. Age at menopause was divided into the categories ≤48, 49–50, 51–54 and ≥55 years.
A basic model stratified by age and center and a multivariable model further adjusted for body mass index (BMI kg/m2, as a continuous variable), height (as a continuous variable), menopausal status at enrolment (premenopausal, peri-/unknown menopausal, postmenopausal [natural and surgical menopause], HRT use (premenopausal, ever use, never use and missing in postmenopausal women only), smoking status (current, former, never, missing), baseline alcohol consumption (non-consumers, 0.1–6 g/day, 6-12 g/day, 12-24 g/day, 24-60 g/day and greater than 60 g/day, missing), physical activity (Cambridge Index: active, moderately active, moderately inactive and inactive, missing ), education level (none, primary school, technical/professional school, secondary school, longer education including university degree, missing) were assessed. Missing values (generally <2%) were accounted for by creating an extra category in each covariate.
To avoid collinearity when studying the joint effect of the number of full-term pregnancies, age at first and last full-term childbirth and time since last childbirth, we used the approach described by Heuch et al.. In analyses including age at last full-term childbirth and time since last childbirth in an age adjusted model, the general age effect was represented by the age effect among nulliparous women. We assigned constant values for age at full-term childbirth and time since last full-term childbirth (corresponding to the reference categories) to nulliparous women, and indicator variables were introduced in the model to ensure that the risk estimates reflected effects in parous women only.
Differences in risk estimates of a given factor and across breast cancer subtypes were analyzed using the data augmentation method as described by Lunn and McNeil, using a likelihood ratio test to compare the model with and without interaction terms between the exposure of interest and breast cancer subtype . Women were considered at risk of a given breast cancer subtype until they were diagnosed with a different competing breast cancer subtype or were diagnosed with breast cancer and the receptor status information was missing. These women were censored at the time of occurrence of the competing breast cancer subtype . To assess whether breast cancer subtype reproductive risk factors changed across women after menopause, left and/or right side censoring was used to count person years within defined age periods <50 years, and ≥50 years. As no differences were observed between risk estimates of reproductive risk factors and breast cancer subtype risk across the age-bands we report results for all women combined. A sensitivity analysis restricting to cases with any indication of an ER and PR expression versus a complete absence of ER and PR expression (0% cells stained, a “-“ description (i.e. a negative/minus symbol description), 0 fmol/mg, an Allred score of 0, an IRS = 0, or an H-score = 0) was also performed. Heterogeneity in the risk associations between subgroups by age at diagnosis (<50 vs. ≥50 years), OC use, center, median BMI, age at first pregnancy and ever having breastfed were also examined using the Cochran’s Q statistic. A previous analysis on postmenopausal HRT use has been reported in the EPIC cohort, therefore HRT use as a predictor of breast cancer risk by HR status was not included in this analysis . All statistical analyses were performed using the SAS software package, version 9.2 (SAS Institute, Cary, NC).
A total cohort of 311,097 women was followed for a sum of 3,346,356 person years. Women were recruited into the EPIC study at the median age of 51.1 years (Table 1). At the time of recruitment, 46.5% of the women were postmenopausal and the median age of menopause was 50.0 years. For parous women, the median age at first full-term birth across the EPIC centers was 24.0 and the median time between menarche and first child birth was 11.0 years. The median age at last full-term pregnancy was 29.0 years, and a median time of 22.9 years had passed since the last childbirth. Of the women who had breastfed (83.8%), the median time of breastfeeding was 6.0 months for all pregnancies.
A statistically significant heterogeneity between the risk of ER-PR- (n = 998) and ER+PR+ (n = 3,567) tumors was observed for age at first full-term childbirth (Table 2). While Cox regression models showed no association with risk of ER-PR-malignancies, a first full-term childbirth after the age of 35 was associated with an increased risk of ER+PR+ tumors (≥35 vs. ≤19 years HR: 1.47 [95% CI 1.15-1.88] p trend < 0.001, P het = 0.03).
For ER-PR- breast tumors, similar risk associations to ER+PR+ breast tumors were observed with increasing menarcheal age (≥15 vs. ≤ 13 years ER-PR- HR: 0.84 [95% CI 0.69-1.03], P trend =0.17; ER+PR+ HR: 0.76 [95% CI 0.68-0.85], P trend <0.001; P het = 0.48), and among parous women, with longer time between menarche and first full-term childbirth (≥10 vs. < 10 years ER-PR- HR: 1.15 [95% CI 0.99-1.34], P trend =0.09; ER+PR+ HR: 1.22 [95% CI 1.12-1.33], P trend <0.001; P het = 0.52). Although the relative risk estimates for ER-PR- breast malignancies were in a similar direction to the ER+PR+ tumors and no statistical heterogeneity between the breast cancer subtypes was observed, it should be noted that risk associations for ER-PR- malignancies were weaker in magnitude and failed to reach statistical significance.
Ever a full-term childbirth, number of childbirths, age- and time since last full-term childbirth were generally not associated with the risk of ER-PR- malignancies but were associated with ER+PR+ tumors, however no statistical heterogeneity between the breast cancer subtypes was observed. Moreover, OC use, duration of OC use, breast feeding, ever having an abortion (spontaneous or induced) and age at menopause showed no significant associations with either ER-PR- or ER+PR+ breast cancer.
When all of the pregnancy related variables were examined together in the same model, risk associations for an increased number of full-term births, and later first and last full-term childbirth with ER-PR- tumors appeared to slightly strengthen; however they remained statistically not significant and had wider confidence intervals (Table 3). When analyses were restricted to postmenopausal women, the risk estimates for a later age at first full-term childbirth with ER+PR+ tumors were stronger (≥35 vs. ≤19 years HR: 1.64 [95% CI 1.05-2.56] p trend = 0.002) compared to the estimates for all the women combined (≥35 vs. ≤19 years HR: 1.30 [95% CI 0.95-1.79] p trend = 0.02).
Subgroup analyses showed no heterogeneity in most of the risk estimates of the reproductive variables for ER-PR- and ER+PR+ tumors (data not shown), with the exception of the duration of breastfeeding by age at first birth among parous women (Pheterogeneity = 0.002). Among women who had their first full-term childbirth before the age of 25, a longer duration of total breastfeeding showed an indication of a decreased risk association with ER+PR+ tumors (≥18vs. ≤ 1 month HR: 0.86 [95% CI 0.67-1.10], Ptrend = 0.01). In contrast, among women who had their first full-term childbirth after the age of 25, a longer cumulative duration of breastfeeding was associated with a particular increased risk of ER+PR+ tumors (≥18vs. ≤ 1 month HR: 1.50 [95% CI 1. 13–1.99], Ptrend = 0.005). For ER-PR- tumors, risk estimates in the similar direction to the ER+PR+ tumors were observed, however statistically not significant (first childbirth before the age of 25, ≥18vs. ≤ 1 month of breast feeding HR: 0.90 [95% CI 0.58-1.39], Ptrend = 0.89; first childbirth after the age of 25, ≥18vs. ≤ 1 month of breast feeding HR: 1.31 [95% CI 0.72-2.39], Ptrend = 0.41).
In the sensitivity analysis restricted to cases with any indication of a positive ER and PR expression versus a complete absence of ER and PR expression with reproductive factors showed similar patterns with risk of joint ER+PR+ and ER-PR- breast cancer subtypes (data not shown).
Our results showed a significantly heterogeneous risk association for age at first full-term pregnancy by receptor status of the tumor, where later first full-term childbirth was associated with increased risk of ER+PR+ tumors but not with risk of ER-PR- tumors. Menarcheal age and time period between menarche and first full-term childbirth showed suggestively similar risk associations with both ER-PR- and ER+PR+ tumors, however the risk estimates for ER-PR- tumors were generally weaker than for their ER+PR+ counterparts and only borderline significant. Although the heterogeneity between breast cancer subtypes was not statistically significant, other parity related factors (such as ever having a full-term childbirth, number of full-term childbirths, age- and time since last full-term childbirth) were associated only with ER+PR+ malignancies. Finally, the factors related to breastfeeding and OC use were generally not associated with HR-positive or HR-negative breast cancer risk.
Previous prospective studies investigating the association of reproductive factors with HR-positive breast cancer have shown relatively consistent inverse risk associations with increasing menarcheal age, ever having a full-term childbirth and particularly a full-term childbirth at an early age [3, 9, 14, 25]. However, consensus on the associations with HR-negative tumors has not been reached because previous prospective studies have lacked sufficient sample sizes [3, 9, 26] and because of the heterogeneous nature of HR-negative subtypes . A more recent study within the Women’s Health Initiative  showed that ever having a full-term childbirth was associated with an increased risk of triple-negative breast cancer (ER-, PR- and HER2-) (n = 307) and the positive association was strengthened with an increasing number of full-term births. We were unable to confirm the positive risk association with ever having a full-term childbirth and increasing number of full-term pregnancies with HR-negative tumors, however we did not have information on triple negative tumors.
We observed for both ER-PR- and ER+PR+ tumors, similar risk associations with increasing menarcheal age and a longer time between menarche and first full-term childbirth. A recent study within the EPIC cohort previously reported on the association of menarche with both ER-PR- and ER+PR+ tumors within the context of childhood growth and earlier sexual maturity . This study extends onto this study of menarche and growth and focuses on the time period between sexual maturity and first full-term pregnancy, thus illustrating the complex and entwined nature of endocrine-sensitive tumors with hormones during different life phases of growth, sexual maturity and reproduction. The inverse association of menarcheal age with increased breast cancer risk is thought to be resultant of a longer exposure to estrogens during a women’s reproductive life  but may also reflect early pubertal years characterized by more intensive and increased exposure to estrogen [27, 28]. Estrogens have been long established to have a late-stage growth promoting effect on estrogen sensitive tumors , however, evidence suggests that estrogens may also play an important role in earlier developmental stages of both HR-positive and -negative tumor types . Mammary stem cells have been shown to respond to sex steroid hormones without having a clear expression of an ER or PR . Further, the EPIC cohort, showed that pre-diagnostic levels of estrogens were associated with both HR-negative and HR-positive postmenopausal breast cancer . The longer time between menarche and a women’s first full-term childbirth would equate to a longer period of time with undifferentiated breast epithelial tissue and a shorter period of time that the breast is resistant to malignant transformation  and thus may have etiological importance in the formation of ER-PR- tumors as well.
We observed a significantly different risk association for a later age at first full-term childbirth with risk of HR-negative and HR-positive tumors, whereas the associations for parity related factors (such as ever having a full-term childbirth, age a last full-term childbirth and time since last full-term childbirth) appeared to aggregate around HR-positive tumors. The role of a pregnancy with the risk of breast cancer is thought to stem from two major avenues, firstly, hormonal changes before and after pregnancy and secondly , dramatic structural changes in the ductal system of the breast after pregnancy [26, 32]. A full-term childbirth is associated with a long term post-pregnancy reduction in levels of circulating hormones . Before a women’s first pregnancy the breast contains a high proportion of undifferentiated ducts and associated alveolar buds. Complete differentiation only occurs during pregnancy and lactation via complex morphological, physiological, and molecular changes . Terminally differentiated epithelial cells are more resistant to carcinogenic influences because of lower proliferation rates and longer DNA repair phases . The distinct inverse risk association for an earlier age at first full-term childbirth with ER+PR+ disease could be due to a shorter exposure to higher levels of ovarian estrogens and a shorter period of time of undifferentiated breast epithelial cells.
Recent prospective studies have reported reduced risk associations with breastfeeding with both ER-PR- and ER+PR+ breast cancer [3, 14]. In the current analysis, we also observed an inverse risk association for both ER-PR- and ER+PR+ malignancies with a longer cumulative duration of breastfeeding however, this was restricted to women who had an early full-term childbirth. In contrast to the recent studies, among women who had a later first full-term childbirth, an increased risk with ER-PR- and ER+PR+ breast tumors with a longer total duration of breastfeeding was observed. Breastfeeding is thought to protect a woman from developing breast cancer by increasing breast differentiation, postponing the return of the ovulatory menstrual cycle post-pregnancy, and/or changing the hormonal environment of the breast [9, 26, 32]. The inverse risk association of HR-negative and HR-positive tumors with breastfeeding coupled with a longer duration among young first time mothers could convey a similar protection. In addition, lactation at a younger age would also mean a shorter period of undifferentiated breast epithelial tissue .
OC use has been extensively studied by many epidemiological studies and most studies have found either no association or a moderate increased risk of overall breast cancer, particularly among very young women and recent OC users . More recent case–control studies investigating the risk of HR-defined breast cancer have started showing relationships of OC use with HR-negative breast cancer [33–37]. In the current study, we were unable to confirm significant risk associations of ever OC use and a longer duration of OC use with HR-negative tumors. There were only a small number of baseline OC users within this analytical cohort, and unfortunately within the EPIC cohort, information of dose or type of OC used, date of last use and information on changes to OC use after baseline are not available. Similarly, the estrogens and progestins in oral contraceptives differ in type and concentration  and this could be hiding a risk association.
Major strengths of this study are its prospective design and large number of incident cases with receptor information. The large case numbers allowed an in-depth analysis of reproductive-related relative risks, describing risk associations among women of predominantly premenopausal and postmenopausal age. The large case numbers also enabled us to examine subgroup effects such as age at first birth, age at diagnosis, BMI and breastfeeding. To our knowledge, this analysis uses the largest number of incident cases of ER-PR- malignancies, although future prospective studies with a greater number of ER-PR- cases are necessary to characterize the associations, which are of substantially smaller magnitude when compared to their HR-positive counterparts. Our study does have its limitations. The determination of ER and PR status in breast tumors has become a standard part of breast cancer diagnosis and is used to predict endocrine therapy response . While a number of studies have shown that the classification of the ER and PR in tumors is relatively robust [38, 39], the accuracy of classifying an ER or PR-negative tumor remains controversial [40, 41]. In the analysis that compared women with a complete absence of ER and PR expression to women with any indication of an ER and PR positive expression, ER-PR- risk estimates remained unchanged. Furthermore, proportions of ER and PR–negative tumors in the EPIC cohort are in line with previous reports [13, 42]. In addition, the inclusion of PR provides an indication of a functional estrogen pathway  and thus a joint ER-PR- may be more reflective of a true ER-negative tumor. At the time of this study, additional information on HER2 expression to determine breast cancer subtypes into more detailed molecular sub-classifications could not be completed because of insufficient information on HER2. However, as the routine assessment of HER2 is relatively more recent than ER and PR assessment, future cohort analyses will be able to include HER2.
In conclusion, our study provides evidence that later age at first full-term childbirth is associated with an increased risk of ER+PR+ tumors but not with ER-PR- tumors. Moreover, age at menarche and time between menarche and first full-term childbirth may be associated with the etiology of both HR-negative and HR-positive malignancies, although associations were only borderline significant for HR-negative tumors. Further studies with more incident cases of ER-PR- tumors are needed to provide more precise risk estimation for reproductive factors with HR-negative tumors.
European Prospective Investigation into Cancer and Nutrition
- 95% CI:
95% confidence interval
Human epidermal growth factor-2
Body mass index.
Weigelt B, Reis-Filho JS: Histological and molecular types of breast cancer: is there a unifying taxonomy?. Nat Rev Clin Oncol. 2009, 6 (12): 718-730. 10.1038/nrclinonc.2009.166.
Althuis MD, Fergenbaum JH, Garcia-Closas M, Brinton LA, Madigan MP, Sherman ME: Etiology of hormone receptor-defined breast cancer: a systematic review of the literature. Cancer Epidemiol Biomarkers Prev. 2004, 13 (10): 1558-1568.
Ma H, Bernstein L, Pike MC, Ursin G: Reproductive factors and breast cancer risk according to joint estrogen and progesterone receptor status: a meta-analysis of epidemiological studies. Breast Cancer Res. 2006, 8 (4): R43-10.1186/bcr1525.
Rakha EA, El-Sayed ME, Green AR, Paish EC, Powe DG, Gee J, Nicholson RI, Lee AH, Robertson JF, Ellis IO: Biologic and clinical characteristics of breast cancer with single hormone receptor positive phenotype. J Clin Oncol. 2007, 25 (30): 4772-4778. 10.1200/JCO.2007.12.2747.
Geyer FC, Reis-Filho JS: Microarray-base gene expression profiling as a clinical tool for breast cancer management: Are we there yet?. Int J Surg Pathol. 2008, 17 (4): 285-302.
Wirapati P, Sotiriou C, Kunkel S, Farmer P, Pradervand S, Haibe-Kains B, Desmedt C, Ignatiadis M, Sengstag T, Schutz F, Goldstein DR, Piccart M, Delorenzi M: Meta-analysis of gene expression profiles in breast cancer: toward a unified understanding of breast cancer subtyping and prognosis signatures. Breast Cancer Res. 2008, 10 (4): R65-10.1186/bcr2124.
Perou CM, Sorlie T, Eisen MB, de RM V, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA, Fluge O, Pergamenschikov A, Williams C, Zhu SX, Lonning PE, Borresen-Dale AL, Brown PO, Botstein D: Molecular portraits of human breast tumours. Nature. 2000, 406 (6797): 747-752. 10.1038/35021093.
Tamimi RM, Baer HJ, Marotti J, Galan M, Galaburda L, Fu Y, Deitz AC, Connolly JL, Schnitt SJ, Colditz GA, Collins LC: Comparison of molecular phenotypes of ductal carcinoma in situ and invasive breast cancer. Breast Cancer Res. 2008, 10 (4): R67-10.1186/bcr2128.
Lord SJ, Bernstein L, Johnson KA, Malone KE, McDonald JA, Marchbanks PA, Simon MS, Strom BL, Press MF, Folger SG, Burkman RT, Deapen D, Spirtas R, Ursin G: 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-1730. 10.1158/1055-9965.EPI-07-2824.
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, 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-263. 10.1093/jnci/djq526.
Millikan RC, Newman B, Tse CK, Moorman PG, Conway K, Dressler LG, Smith LV, Labbok MH, Geradts J, Bensen JT, Jackson S, Nyante S, Livasy C, Carey L, Earp HS, Perou CM: Epidemiology of basal-like breast cancer. Breast Cancer Res Treat. 2008, 109 (1): 123-139. 10.1007/s10549-007-9632-6.
Gierach GL, Burke A, Anderson WF: Epidemiology of triple negative breast cancers. Breast Dis. 2010, 32 (1–2): 5-24.
Yasui Y, Potter JD: The shape of age-incidence curves of female breast cancer by hormone-receptor status. Cancer Causes Control. 1999, 10 (5): 431-437. 10.1023/A:1008970121595.
Phipps AI, Chlebowski RT, Prentice R, McTiernan A, Wactawski-Wende J, Kuller LH, Adams-Campbell LL, Lane D, Stefanick ML, Vitolins M, Kabat GC, Rohan TE, Li CI: Reproductive history and oral contraceptive use in relation to risk of triple-negative breast cancer. J Natl Cancer Inst. 2011, 103 (6): 470-477. 10.1093/jnci/djr030.
Phipps AI, Malone KE, Porter PL, Daling JR, Li CI: Reproductive and hormonal risk factors for postmenopausal luminal, HER-2-overexpressing, and triple-negative breast cancer. Cancer. 2008, 113 (7): 1521-1526. 10.1002/cncr.23786.
Riboli E, Kaaks R: The EPIC project: rationale and study design. European Prospective Investigation into Cancer and nutrition. Int J Epidemiol. 1997, 26 (1): S6-S14.
Riboli E, Hunt KJ, Slimani N, Ferrari P, Norat T, Fahey M, Charrondiere UR, Hemon B, Casagrande C, Vignat J, Overvad K, Tjonneland A, Clavel-Chapelon F, Thiebaut A, Wahrendorf J, Boeing H, Trichopoulos D, Trichopoulou A, Vineis P, Palli D, Bueno-de-Mesquita HB, Peeters PH, Lund E, Engeset D, Gonzalez CA, Barricarte A, Berglund G, Hallmans G, Day NE, Key TJ, et al: European Prospective Investigation into Cancer and Nutrition (EPIC): study populations and data collection. Public Health Nutr. 2002, 5 (6B): 1113-1124. 10.1079/PHN2002394.
Onland-Moret NC, Peeters PH, van Gils CH, Clavel-Chapelon F, Key T, Tjonneland A, Trichopoulou A, Kaaks R, Manjer J, Panico S, Palli D, Tehard B, Stoikidou M, Bueno-de-Mesquita HB, Boeing H, Overvad K, Lenner P, Quiros JR, Chirlaque MD, Miller AB, Khaw KT, Riboli E: Age at menarche in relation to adult height: the EPIC study. Am J Epidemiol. 2005, 162 (7): 623-632. 10.1093/aje/kwi260.
Haftenberger M, Lahmann PH, Panico S, Gonzalez CA, Seidell JC, Boeing H, Giurdanella MC, Krogh V, Bueno-de-Mesquita HB, Peeters PH, Skeie G, Hjartaker A, Rodriguez M, Quiros JR, Berglund G, Janlert U, Khaw KT, Spencer EA, Overvad K, Tjonneland A, Clavel-Chapelon F, Tehard B, Miller AB, Klipstein-Grobusch K, Benetou V, Kiriazi G, Riboli E, Slimani N: Overweight, obesity and fat distribution in 50- to 64-year-old participants in the European Prospective Investigation into Cancer and Nutrition (EPIC). Public Health Nutr. 2002, 5 (6B): 1147-1162. 10.1079/PHN2002396.
James RE, Lukanova A, Dossus L, Becker S, Rinaldi S, Tjonneland A, Olsen A, Overvad K, Mesrine S, Engel P, Clavel-Chapelon F, Chang-Claude J, Vrieling A, Boeing H, Schutze M, Trichopoulou A, Lagiou P, Trichopoulos D, Palli D, Krogh V, Panico S, Tumino R, Sacerdote C, Rodriguez L, Buckland G, Sanchez MJ, Amiano P, Ardanaz E, Bueno-de-Mesquita B, Ros MM, et al: Postmenopausal serum sex steroids and risk of hormone receptor-positive and -negative breast cancer: a nested case–control study. Cancer Prev Res. 2011, 4 (10): 1626-1635. 10.1158/1940-6207.CAPR-11-0090.
Wareham NJ, Jakes RW, Rennie KL, Schuit J, Mitchell J, Hennings S, Day NE: Validity and repeatability of a simple index derived from the short physical activity questionnaire used in the European Prospective Investigation into Cancer and Nutrition (EPIC) study. Public Health Nutr. 2003, 6 (4): 407-413.
Heuch I, Albrektsen G, Kvale G: Modeling effects of age at and time since delivery on subsequent risk of cancer. Epidemiology. 1999, 10 (6): 739-746. 10.1097/00001648-199911000-00015.
Lunn M, McNeil D: Applying Cox regression to competing risks. Biometrics. 1995, 51 (2): 524-532. 10.2307/2532940.
Ritte R, Lukanova A, Berrino F, Dossus L, Tjonneland A, Olsen A, Overvad TF, Overvad K, Clavel-Chapelon F, Fournier A, Fagherazzi G, Rohrmann S, Teucher B, Boeing H, Aleksandrova K, Trichopoulou A, Lagiou P, Trichopoulos D, Palli D, Sieri S, Panico S, Tumino R, Vineis P, Quiros JR, Buckland G, Sanchez MJ, Amiano P, Chirlaque MD, Ardanaz E, Sund M, et al: Adiposity, hormone replacement therapy use and breast cancer risk by age and hormone receptor status: a large prospective cohort study. Breast Cancer Res. 2012, 14 (3): R76-10.1186/bcr3186.
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-1259. 10.1093/aje/kwp036.
Hankinson SE, Tamimi R, Hunter D: Breast cancer. Textbook of Cancer Epidemiology. Edited by: Adami H, Hunter D, Trichopoulos D. 2008, Oxford: Oxford University Press, 403-445. 2
Ritte R, Lukanova A, Tjonneland A, Olsen A, Overvad K, Mesrine S, Fagherazzi G, Dossus L, Teucher B, Steindorf K, Boeing H, Aleksandrova K, Trichopoulou A, Lagiou P, Trichopoulos D, Palli D, Grioni S, Mattiello A, Tumino R, Sacerdote C, Quiros JR, Buckland G, Molina-Montes E, Chirlaque MD, Ardanaz E, Amiano P, Bueno-de-Mesquita B, Van DF, Van Gils CH, Peeters PH, et al: Height, age at menarche and risk of hormone receptor-positive and -negative breast cancer: a cohort study. Int J Cancer J Int du cancer. 2013, 132 (11): 2619-2629. 10.1002/ijc.27913.
Apter D, Reinila M, Vihko R: Some endocrine characteristics of early menarche, a risk factor for breast cancer, are preserved into adulthood. Int J Cancer. 1989, 44 (5): 783-787. 10.1002/ijc.2910440506.
Dietel M, Lewis MA, Shapiro S: Hormone replacement therapy: pathobiological aspects of hormone-sensitive cancers in women relevant to epidemiological studies on HRT: a mini-review. Hum Reprod. 2005, 20 (8): 2052-2060. 10.1093/humrep/dei043.
Allred DC, Brown P, Medina D: The origins of estrogen receptor alpha-positive and estrogen receptor alpha-negative human breast cancer. Breast Cancer Res. 2004, 6 (6): 240-245. 10.1186/bcr938.
Asselin-Labat ML, Vaillant F, Sheridan JM, Pal B, Wu D, Simpson ER, Yasuda H, Smyth GK, Martin TJ, Lindeman GJ, Visvader JE: Control of mammary stem cell function by steroid hormone signalling. Nature. 2010, 465 (7299): 798-802. 10.1038/nature09027.
Russo IH, Russo J: Pregnancy-induced changes in breast cancer risk. J Mammary Gland Biol Neoplasia. 2011, 16 (3): 221-233. 10.1007/s10911-011-9228-y.
Ma H, Bernstein L, Ross RK, Ursin G: Hormone-related risk factors for breast cancer in women under age 50 years by estrogen and progesterone receptor status: results from a case–control and a case-case comparison. Breast Cancer Res. 2006, 8 (4): R39-10.1186/bcr1514.
Dolle JM, Daling JR, White E, Brinton LA, Doody DR, Porter PL, Malone KE: Risk factors for triple-negative breast cancer in women under the age of 45 years. Cancer Epidemiol Biomarkers Prev. 2009, 18 (4): 1157-1166. 10.1158/1055-9965.EPI-08-1005.
Althuis MD, Brogan DD, Coates RJ, Daling JR, Gammon MD, Malone KE, Schoenberg JB, Brinton LA: Breast cancers among very young premenopausal women (United States). Cancer Causes Control. 2003, 14 (2): 151-160. 10.1023/A:1023006000760.
Sweeney C, Giuliano AR, Baumgartner KB, Byers T, Herrick JS, Edwards SL, Slattery ML: Oral, injected and implanted contraceptives and breast cancer risk among U.S. Hispanic and non-Hispanic white women. Int J Cancer. 2007, 121 (11): 2517-2523. 10.1002/ijc.22970.
Rosenberg L, Boggs DA, Wise LA, Adams-Campbell LL, Palmer JR: Oral contraceptive use and estrogen/progesterone receptor-negative breast cancer among African American women. Cancer Epidemiol Biomarkers Prev. 2010, 19 (8): 2073-2079. 10.1158/1055-9965.EPI-10-0428.
Magne N, Toillon RA, Castadot P, Ramaioli A, Namer M: Different clinical impact of estradiol receptor determination according to the analytical method: a study on 1940 breast cancer patients over a period of 16 consecutive years. Breast Cancer Res Treat. 2006, 95 (2): 179-184. 10.1007/s10549-005-9065-z.
Chebil G, Bendahl PO, Idvall I, Ferno M: Comparison of immunohistochemical and biochemical assay of steroid receptors in primary breast cancer–clinical associations and reasons for discrepancies. Acta Oncol. 2003, 42 (7): 719-725. 10.1080/02841860310004724.
Welsh AW, Moeder CB, Kumar S, Gershkovich P, Alarid ET, Harigopal M, Haffty BG, Rimm DL: Standardization of estrogen receptor measurement in breast cancer suggests false-negative results are a function of threshold intensity rather than percentage of positive cells. J Clin Oncol. 2011, 29 (22): 2978-2984. 10.1200/JCO.2010.32.9706.
Hede K: Breast cancer testing scandal shines spotlight on black box of clinical laboratory testing. J Natl Cancer Inst. 2008, 100 (12): 836-837. 10.1093/jnci/djn200.
Chu KC, Anderson WF, Fritz A, Ries LA, Brawley OW: Frequency distributions of breast cancer characteristics classified by estrogen receptor and progesterone receptor status for eight racial/ethnic groups. Cancer. 2001, 92 (1): 37-45. 10.1002/1097-0142(20010701)92:1<37::AID-CNCR1289>3.0.CO;2-F.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/13/584/prepub
We would like to thank Sabine Rohrmann, Jutta Schmitt and Jutta Kneisel for their assistance during the collection of hormone receptor status data, and we thank all the EPIC cohort participants for their contributions to data collection at baseline recruitment and during follow-up. Finally, the comments from the anonymous reviewers are also greatly acknowledged.
This work was (partly) supported by a grant from the German Research Foundation, Graduiertenkolleg 793: Epidemiology of communicable and chronic noncommunicable diseases and their interrelationships. The coordination of EPIC is financially supported by the European Commission (DG-SANCO) and the International Agency for Research on Cancer. The national cohorts are supported by Danish Cancer Society (Denmark); Ligue contre le Cancer, Mutuelle Générale de l’Education Nationale, Institut National de la Santé et de la Recherche Médicale (France); Deutsche Krebshilfe, Deutsches Krebsforschungszentrum and Federal Ministry of Education and Research (Germany); the Hellenic Health Foundation and the Stavros Niarchos Foundation (Greece); Italian Association for Research on Cancer (AIRC) and National Research Council (Italy); Dutch Ministry of Public Health, Welfare and Sports (VWS), Netherlands Cancer Registry (NKR), LK Research Funds, Dutch Prevention Funds, Dutch ZON (Zorg Onderzoek Nederland), World Cancer Research Fund (WCRF), Statistics Netherlands (The Netherlands); ERC-2009-AdG 232997 and Nordforsk, (Norway); Health Research Fund (FIS), Regional Governments of Andalucía, Asturias, Basque Country, Murcia (no. 6236) and Navarra, ISCIII RTICC ’Red Temática de Investigación Cooperativa en Cáncer (R06/0020) (Spain); Swedish Cancer Society, Swedish Scientific Council and Regional Government of Skåne and Västerbotten (Sweden); Cancer Research UK, Medical Research Council, (United Kingdom).
The authors declare they have no competing interests.
RR, AL, and RK contributed to the conception of the current analysis and all authors were involved in the design and acquisition of data from the EPIC cohort. RR, AL and RK contributed to the analysis and all authors contributed to the interpretation of the data. RR, KT, AL and RK drafted the manuscript and all authors revised the final draft critically for important critical content. All authors have given final approval of the version to be published.
Electronic supplementary material
Additional file 1: Table S1: Reproductive factors and risk of ER-positive vs. ER-negative and PR-positive vs. PR-negative breast cancer in all women. (DOCX 49 KB)
Additional file 2: Table S2: Reproductive factors and risk of discordant breast cancer subtypes in all women. (DOCX 32 KB)
About this article
Cite this article
Ritte, R., Tikk, K., Lukanova, A. et al. Reproductive factors and risk of hormone receptor positive and negative breast cancer: a cohort study. BMC Cancer 13, 584 (2013). https://doi.org/10.1186/1471-2407-13-584
- Reproductive factors
- Risk factors
- Oral contraceptive
- Breast cancer