Introduction

Breast cancer is strongly related to various hormonal factors. We know that different hormones play a key role in the pathogenesis of invasive breast cancer in both, premenopausal and postmenopausal women [14], likely representing a direct effect of circulating hormones [57].

The patterns of risk associated with reproductive history suggest that a prolonged exposure to ovarian hormones increases breast cancer occurrence [8]. It has been well established that early age at menarche, late age at menopause, null parity and late age at first full term pregnancy (FFTP) are important risk factors for breast cancer [911]. These characteristics may have long lasting effects on hormone metabolism, although data regarding reproductive characteristics and hormone levels are still inconclusive [1216]

In an effort to measure the lifetime exposure to endogenous hormones, the concept of the “lifetime cumulative number of menstrual cycles” was introduced as an index that could reflect total ovarian exposure [17, 18]. The “lifetime cumulative number of cycles” summarizes information such as menarche, parity, menopause, regularity pattern, breast feeding and oral contraceptive (OC) use. In previous studies we, but also others, related the cumulative numbers of menstrual cycles to breast cancer risk [17, 18]. Ovarian function, and therefore, “lifetime number of menstrual cycles” is defined by age at menopause, and age at menopause is related to postmenopausal estrogen levels [19]; with this rationale we studied the association between the “lifetime number of menstrual cycles” and hormone levels after the menopause. The hypothesis fits in, where a common factor at hypothalamic level may affect the speed by which oocytes are being used, increasing postmenopausal gonoadothopin levels and maybe, stimulating ovarian testosterone release from the remaining theca cells, as well as non ovarian estrogens from adrenal origin.

The aim of this study was to evaluate the relation between the “lifetime cumulative number of menstrual cycles” and the reproductive characteristics that contribute to this estimate on one side, and plasma concentrations of sex hormones on the other side, in order to elucidate the mechanism by which reproductive risk factors could influence hormone-related cancers.

Study design and methods

The Prospect-EPIC study is one of the two Dutch cohorts participating in the European Prospective Investigation into Cancer and Nutrition; a multi-center prospective study carried out in 23 centers from 10 European countries [20].

A detailed description of the study population and data collection of the Prospect-EPIC study has been published elsewhere [21]. Briefly, 17,357 women aged 49–70, residing in Utrecht and its surroundings were recruited through a regional, population-based, breast cancer screening program between 1993 and 1997. Base-line characteristics were collected on the basis of two self-administered questionnaires (regarding lifestyle and nutrition). The lifestyle questionnaire contains information on demographic characteristics, family history, past and current morbidity, reproductive history (gynecologic and obstetric), smoking habits, and physical activity. In addition, anthropometric measurements (height, weight, and waist and hip circumference) were taken as well as a 30 ml blood sample. Within 24 h, samples of 4 ml serum, 9 ml citrate plasma, 2 ml white blood cells, and 2 ml red blood cells were fractionated into 0.5 ml aliquots and stored at −196°C. Blood samples were successfully drawn from 97.5% of the participants [21]. All women signed a written informed consent and the study was approved by the Institutional Review Board of the University Medical Center Utrecht.

Of all Prospect-EPIC participants a 10% random sample (n = 1736) was taken and hormone measurements were done. For the present study, participants who used hormone replacement therapy (HRT) or OC in the previous twelve months were excluded (n = 311), as well as those women with history any cancer other than non-melanoma cancer (n = 25). From the remaining 1,400 participants only those who underwent natural menopause were eligible (n = 860) because only in them we had reliable information about ovarian activity not masked by the absence of menstrual cycles due to hysterectomy. Natural menopause was defined as the complete cessation of menstrual bleeding for 12 months preceding enrollment due to natural causes (according to self-report).

We computed the “lifetime cumulative number of menstrual cycles” of the women in this study. For parous women, the number of menstrual cycles before a FFTP was also estimated as we previously described [17]. Briefly, the “lifetime cumulative number of menstrual cycles” was computed from age at first menstrual period to age at menopause. Ages at menarche, at menopause and at FFTP were recorded as continuous variables. Age at regularization of cycles was computed according to the time between menarche and the establishment of regular periods (immediately after menarche, after 1, 2, 3, 5, or 10 years, after FFTP or never regular). The regularity pattern was defined categorically as 24, 26, 28, 30 and 32 day cycles, or always irregular. Information about OC use included: age at start, number of years of use and age at last use. Data about number of and age at each pregnancy, stillbirth, miscarriage and abortion, if any, were available, as well as the number of weeks of breastfeeding for each pregnancy.

We took all the previously mentioned variables into account when computing the “lifetime cumulative number of menstrual cycles”, considering that each pregnancy, stillbirth, miscarriage or abortion, resulted in an absence of cycles for a period of 36, 28 and 12 weeks respectively. We assumed that lactation results in a 6-week absence of cycles [22]. Menstrual cycles due to OC use were computed as 28 days duration (see Appendix). Considering that irregular cycles are a combination of long and short menstrual periods, the number of irregular cycles was computed according to each participant’s regularity pattern; the total number of irregular years was recorded in order to make further adjustments. The number of menstrual cycles was impossible to compute in the participants who reported to be always irregular or never having regular periods, therefore we added the “irregular” category as such in the analysis (See appendix).

Bounded levels of estrogens and androgens were measured, levels of estrone (E1), estradiol (E2), androstenedione, testosterone, sex hormone binding globulin (SHBG) and dehydroepiandrostenedione sulfate (DHEAS) were measured in plasma using commercially available double-antibody radioimmunoassay (RIA) kits (Diagnostic System Laboratories Inc., Webster, TX, USA). The following kits were used: E1: DSL-8700, E2: DSL-39100, androstenedione: DSL-4200, TST: DSL-4100 SHBG: DSL-6300 and DHEAS: DSL-2700. The intra-assay coefficients of variation were 5.6%, 3.9%, 4.3%, 7.7%, 3.0% and 5.2% respectively. The inter-assay coefficients of variation were 11.1%, 4.1%, 6.3%, 8.1%, 4.0% and 5.3% respectively. All hormone assays were performed by the laboratory of the Hormones and Cancer group at the International Agency for Research on Caner (IARC). Although technically SHBG is not a hormone, for reasons of convenience it will be referred as such. We calculated the free-androgen index (FAI) dividing testosterone by SHBG and multiplying by 100, the SHGB/E2 index was calculated dividing SHBG by E2.

Data analysis

Means with their standard deviation (SD) values, or frequencies of baseline characteristics were calculated. Concentrations of all hormones were logarithmically transformed to produce approximately normal distributions. To make the results easily interpretable, we transformed the mean and its 95% confidence interval (CI) back to their original scale, resulting in geometric means and 95% CIs. Geometric mean levels and 95% CIs were calculated for each hormone (SHBG, E1, E2, androstenedione, testosterone, DHEAS), SHBG/estradiol index (SHBG/E2) and free androgen index (FAI).

Using analysis of covariance models, geometric mean hormone levels were calculated for the different categories of cumulative number of menstrual cycles (lifetime and before a FFTP). The same analysis was also used for categories of age at menarche, age at menopause, parity, number of children and age of FFTP. In order to build the adjusted model we classified confounders in three categories as follows: Group A contained characteristics clearly related to hormone levels such as age, years since menopause and body mass index (BMI); Group B contained factors possibly related to hormone levels like OC and HRT use, smoking habits, physical activity and alcohol consumption; group C contained characteristics unlikely related to hormone levels such as marital status, education level and family history of breast cancer. We first adjusted for all confounders of type A, B, and C, then we built models by deleting B and C types of characteristics, without affecting validity but increasing precision. In the final model we included confounders in group A and B, therefore we adjusted for age at intake as a continuous variable, number of years from menopause to inclusion in the study, BMI at examination, use of HRT ever (Y/N), use of OC ever (Y/N), smoking habits, alcohol consumption and physical activity. Information on current alcohol use was obtained from the food frequency questionnaire, total alcohol consumption in grams per day was calculated from questions in the number of glasses of wine, beer, fortified wines, and liquor/spirits per day, week, month or year. One glass of any alcoholic beverage was assumed to contain 10g of alcohol. With the daily physical activity data from the general questionnaire, the Voorrips score was calculated (and was entered in the model as a continuous variable), this score has been developed and validated especially for women in this age group [23]. For analysis of the “lifetime cumulative number of menstrual cycles” and the “number of cycles before a FFTP” variables, we also adjusted for more than one year of irregular cycles (Y/N). Given that age at menopause can be determined using age at intake in combination of years since menopause; when analyzing age at menopause, the covariate years since menopause was not entered in the model.

To test the significance of the models’ linearity, median values in each category for every variable were computed accordingly. Subsequently these values were included as a continuous variable in a linear regression model to test for linear trends. Because there is no clear natural hierarchy, the “always irregular” category was not included when testing for trend for the “lifetime number of menstrual cycles” and the “number of menstrual cycles before a FFTP” variables. Additional mutual adjustments for reproductive characteristics did not show any difference in the observed trends. The Statistical Package for Social Sciences for WINDOWS, version 14.0 (Chicago, IL), was used for all statistical analysis, and all tests were two sided.

Results

General and reproductive characteristics of the Prospect-EPIC women that were included in this study are shown on Table  1. On average, participants were 60.2 years old at enrollment; had a BMI of 26.2 kg/m2 and the majority of them (69.4%) had at least low technical-secondary education. 113 (13.1 %) were nulliparous, and for the parous, the mean age for a FFTP was 25.5 years old. Mean age at menarche and menopause were 13.6 and 49.6 years old respectively, the time from menopause to inclusion, was on average 10.6 years. In 786 participants (91.4%) the lifetime cumulative number of cycles was calculated, as we had incomplete data regarding reproductive history in 74 women. The median number of lifetime menstrual cycles was 451.3, corresponding to 34.7 years of menstrual activity if considering an average of 13 cycles per year. From the 747 parous participants, 64 had incomplete data to compute the number of menstrual cycles before a FFT, in the remaining 683 participants (91.4%), the median number of menstrual cycles before a FFTP was 146.7 (approximately 11.3 years of menstrual activity). These results are consistent with the previous data on such estimates for a larger sample of the Prospect-EPIC participants as we previously published [17]. Table 2 shows the geometric mean hormone levels for all participants with 95% CIs of SHBG, E1, E2, androstenedione, testosterone, DHEAS, SHBG/E2 and FAI.

Table 1 General and reproductive characteristics of the study participants (n = 860)
Table 2 Geometric mean serum sex hormone and SHBG levels

Table 3 shows the unadjusted geometric mean hormone levels and 95% CI for each hormone according to some general and reproductive characteristics. When testing for trend participants in the lower quartile of “lifetime cumulative of menstrual cycles” had higher levels of SHBG, androstenedione, testosterone and DHEAS than the participants with the highest quartile (P trend = 0.015, 0.003, 0.052 and 0.035 respectively), and lower levels of SHBG/E2 (P trend = 0.031).

Table 3 Hormones unadjusted means (95%CI) according to patient characteristics

In table 4 adjusted geometric means for each hormone according to the reproductive characteristics are shown. Participants in the lowest quartile of cumulative number of menstrual cycles tended to have higher SHBG (P trend = 0.013) and lower SHBG/E2 (P trend = 0.050) levels. Androstenedione and DHEAS levels increased with increasing age at menarche (P trend = 0.027 and 0.022 respectively), a similar tendency was observed for E2, although the P for trend did not reach statistical significance (P trend = .058). Women with more children, tended to have higher SBHG levels (P trend = 0.030), number of children was also related to lower SHBG/E2 (p trend = 0.010) levels. No association was found between number of cycles before a FFTP, parity, age at FFTP and any of the hormones. Age at menopause related inversely with androstenedione and testosterone; women with older age at menopause tend to have lower levels of such hormones (p trend = 0.024 and 0.015 respectively).

Table 4 Hormones adjusted means (95%CI) for serum sex hormones and SHBG levels according to reproductive characteristics

Discussion

The results of the present study suggest that the “lifetime cumulative number of cycles” is related to SHBG and that some reproductive characteristics are also associated to postmenopausal serum sex hormone levels. We noted significant associations after adjusting for important confounders, the “lifetime cumulative number of menstrual cycles” related inversely with levels of SHBG and positively with SHBG/E2. Age at menarche directly relates to both DHEAS and androstenedione. Number of children had a significant relation with levels of SHBG and SHBG/E2 and age at menopause was indirectly associated to androstenedione and testosterone.

Several studies have investigated the relation between reproductive risk factors and endogenous sex hormone levels in postmenopausal women, but with conflicting results [6, 13, 15, 2427]. The patterns of risk associated with reproductive history suggest that prolonged exposure to ovarian hormones increases breast cancer occurrence [8, 28, 29]. The cumulative number of menstrual cycles has been used as a proxy for ovarian activity, women in the highest quartile of the “lifetime cumulative number of menstrual cycles” estimate, have an increased breast cancer risk when compared with women in the lowest quartile [17, 18]. A significant inverse relation between number of menstrual cycles and SHBG was observed in our study. It has been previously described that SHBG decrease the proportion of E2 that is able to leave the circulation and enter the cells, and an inverse relation between concentration of SHBG and breast cancer has been observed in epidemiological studies [1, 10]. Our findings suggest therefore, that postmenopausal free levels of estrogens or androgens may be related to the “lifetime cumulative number of menstrual cycles” reflecting lifetime exposure to ovarian hormones.

Early age at menarche and late age at menopause have long being recognized as risk factors for breast cancer. A 1-year delay in the onset of menarche is associated with a 5% reduction in risk for developing breast cancer in later life [30]. We found that age at menarche has a direct relation with serum levels of DHEAS and androstenedione measured after the menopause. Androgens have been associated with breast cancer risk, women with breast cancer tend to have higher levels of DHEAS and androstenedione [2, 4, 31, 32]. It has been proposed that in postmenopausal women DHEAS might act as a moderate estrogen agonist, through competitive binding of its metabolite 5-androstene-3B,17B-diol [33]. It is also well known that the peripheral conversion of adrenal androstenedione is the principal source of E1 in postmenopausal women [34], therefore we would expect an inverse association between both DHEAS and androstenedione and age at menarche. Madigan et al. showed, after adjusting for BMI and years since menopause, that young age at menarche, was associated with higher levels of androstenedione, but also of E1 and E2 [15]. Our results are in conflict with such findings, however several studies have found no association between age at menarche and SHBG, estrogens or androgens [14, 16, 35, 36].

Childbearing and a young age at FFTP protect against breast cancer and additional pregnancies further reduce the risk [7, 10, 37, 38]. We found that parous women tend to have lower levels of E1, however this difference was not statistically significant (P = .072). Lamar et al. [14] found that nulliparous women had higher testosterone levels, a finding that could contribute to the increased risk of breast cancer in nulliparous women. Our findings also showed that the number of live births is directly related to SHBG, and indirectly to SHBG/E2 and FAI (P for trend for FAI did not reach statistical significance, P = 0.064). This observation is in conflict with the known association between low SBGH levels and breast cancer risk, and the protective effect that the number of pregnancies has on breast cancer [28, 29, 38, 39]. It is important to mention that a woman’s genetic and environmental profile may influence both, her ability to have children and her postmenopausal sex hormone concentrations, making these observations difficult to interpret. It has been hypothesized that the overall effect of parity is due to a protective effect of a young age at FFTP [40], another proposed mechanism is through lactation, as breastfeeding induces final differentiation of the terminal duct epithelium making it relatively insensitive to hormonal stimulation [41, 42]. Despite these observations we found no relation between age at FFTP and any of the hormones. No relation between parity, number of pregnancies or age at FFTP and serum hormone levels was reported by several studies. [12, 1416, 43, 44].

We observed an indirect association between age at menopause and both, androstenedione and testosterone. Several studies have failed in finding an association between age at menopause and sex hormones [12, 14, 16], however a direct association with menopause and androstenedione has been reported [15], as well as a relation between age at menopause and E2 in a group of 173 postmenopausal overweight women [13].

To appreciate the findings of our study, some strengths and limitations need to be mentioned. The main advantage is our large sample size, to the best of our knowledge this is the largest study assessing the relationship between reproductive risk factors and serum sex hormone levels in naturally menopausal women and the first one evaluating the relationship between sex hormones and an index that can be used as a proxy of ovarian exposure. Another advantage is that we attempted to adjust for all relevant covariates, in order to obtain fully adjusted estimates. The participants in this cross-sectional study originated from a random sample of a large cohort, and blood samples where drawn from 97.5% of the 17,357 participants, favoring the generalization of our results. Given the prospective nature of the Prospect-EPIC study, and because participants were recruited through an existing population-based breast cancer screening program, we included only disease-free participants at baseline blood collection.

A limitation of this study is that only a single blood sample was assayed per participant, however it has been shown that a single measure provides substantial information in postmenopausal women for E1, E2 and SHBG [45]. Serial measurements of serum hormones in postmenopausal women suggest that a single measurement can reliably categorize women, at least over the short term. Intraclass correlation coefficients of hormone measurements in serum collected over a 3-year period were reported to be between 0.68 and 0.78 for the estrogens, 0.66 and 0.88 for the androgens, and 0.92 for SHBG [46]. We did not have information regarding premenopausal sex hormone concentrations in the participants of this study; it remains unknown whether some underlying hormonal profile could have influenced both, reproductive capabilities and postmenopausal sex hormone concentrations. It is also possible that premenopausal hormone levels are more related to reproductive factors than postmenopausal levels, and that is as a proxy for premenopausal levels that postmenopausal levels relate to reproductive characteristics. We should consider that the data used regarding reproductive characteristics are prone to recall errors. However, it is unlikely that recall errors are related to endogenous hormone levels themselves, if any misclassification occurred, it would have been non-differential.

In summary, in this group of naturally postmenopausal women, the “lifetime cumulative number menstrual cycles” was related only to SHBG, suggesting that free levels of estrogens or androgens may be related to this index, reflecting lifetime exposure to ovarian hormones. Our findings show that some reproductive characteristics are related to postmenopausal serum sex hormone levels, however we do not find clear evidence to explain the mechanism by which reproductive risk factors influence hormone-related cancers in naturally postmenopausal women. It seems that reproductive and menstrual characteristics do not affect postmenopausal breast cancer risk via a long-term effect on sex hormone levels. The association of reproductive risk factors and breast cancer could be a result of prolonged exposure to ovarian hormones, and not only the product of a direct, long-lasting effect of sex hormone serum concentrations. Further research is needed in order to unravel the biological mechanisms through which reproductive characteristics relate to breast cancer risk.