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

Reproductive aging-associated common genetic variants and the risk of breast cancer

Abstract

Introduction

A younger age at menarche and an older age at menopause are well established risk factors for breast cancer. Recent genome-wide association studies have identified several novel genetic loci associated with these two traits. However, the association between these loci and breast cancer risk is unknown.

Methods

In this study, we investigated 19 and 17 newly identified single nucleotide polymorphisms (SNPs) from the ReproGen Consortium that have been associated with age at menarche and age at natural menopause, respectively, and assessed their associations with breast cancer risk in 6 population-based studies among up to 3,683 breast cancer cases and 34,174 controls in white women of European ancestry. In addition, we used these SNPs to calculate genetic risk scores (GRSs) based on their associations with each trait.

Results

After adjusting for age and potential population stratification, two age at menarche associated SNPs (rs1079866 and rs7821178) and one age at natural menopause associated SNP (rs2517388) were associated with breast cancer risk (p values, 0.003, 0.009 and 0.023, respectively). The odds ratios for breast cancer corresponding to per-risk-allele were 1.14 (95% CI, 1.05 to 1.24), 1.08 (95% CI, 1.02 to 1.15) and 1.10 (95% CI, 1.01 to 1.20), respectively, and were in the direction predicted by their associations with age at menarche or age at natural menopause. These associations did not appear to be attenuated by further controlling for self-reported age at menarche, age at natural menopause, or known breast cancer susceptibility loci. Although we did not observe a statistically significant association between any GRS for reproductive aging and breast cancer risk, the 4th and 5th highest quintiles of the younger age at menarche GRS had odds ratios of 1.14 (95% CI, 1.01 to 1.28) and 1.13 (95% CI, 1.00 to 1.27), respectively, compared to the lowest quintile.

Conclusions

Our study suggests that three genetic variants, independent of their associations with age at menarche or age at natural menopause, were associated with breast cancer risk and may contribute modestly to breast cancer risk prediction; however, the combination of the 19 age at menarche or the 17 age at natural menopause associated SNPs did not appear to be useful for identifying a high risk subgroup for breast cancer.

Introduction

A younger age at menarche and an older age at menopause are well-established risk factors for the development of breast cancer [1]. In the general population, the risk of breast cancer decreases by 10% for each 2-year delay in menarche [2] but increases by 3% for each year that menopause is delayed [3]. These associations are consistent with the hypothesis that breast cancer risk is related to the extent of steroid hormone exposure during a woman's reproductive years, which drives breast mitotic activity and determines the probability of tumorigenic somatic events [4].

Recently, genome-wide association studies (GWAS) have identified several new common genetic loci associated with either age at menarche or age at natural menopause. Four independent GWAS of age at menarche have identified two novel loci at LIN28B and 9q31.2 [5–8], and two GWAS of age at natural menopause have identified four novel loci on chromosomes 5, 6, 19, and 20 [5, 9]. Most recently, the ReproGen Consortium, which consisted of these initial GWASs and many additional studies, has conducted expanded meta-analyses for age at menarche [10] and age at natural menopause [11] and reported more new loci identified for each trait. Given the well-established associations of age at menarche and age at natural menopause with breast cancer risk, we set out to assess whether these common genetic loci influence breast cancer risk and whether a genetic risk score (GRS) for these reproductive events might be useful for identifying a high-risk subgroup for breast cancer. Furthermore, since the reproductive risk factors have been observed to be differentially associated with breast cancer by tumor histological subtypes [12–16], we assessed these genetic associations by tumor histological subtypes defined by estrogen receptor (ER) status.

We therefore conducted a meta-analysis of six population-based studies to investigate the association between genetic loci associated with age at menarche or age at natural menopause and breast cancer risk. We assessed 19 and 17 single-nucleotide polymorphisms (SNPs) that have been previously reported to be linked to age at menarche [10] and age at natural menopause [11], respectively, among up to 3,683 breast cancer cases and 34,174 controls in women of European ancestry and evaluated whether these SNPs were differentially associated with breast cancer subtypes defined by ER status in two studies in which such data were available.

Materials and methods

Study population

The ReproGen Consortium was formed by more than 30 studies in the US and Europe to investigate the genetics of reproductive aging traits [10, 11]. Our analysis used data from six population-based studies from the ReproGen Consortium: the Nurses' Health Study (NHS), the Women's Genome Health Study (WGHS), the SardiNIA Breast Cancer Study (SardiNIA), the Rotterdam Study I and II (RSI+II), the Framingham Heart Study (FHS), and the Atherosclerosis Risk in Communities Study (ARIC). Each study had at least 200 breast cancer cases. Four studies were prospective cohort studies, one was a nested case-control study, and one was a case-control study. A description of the six studies is provided in Table 1, and more information is given in Additional file 1. Briefly, breast cancer cases occurring in defined populations during specific periods of time were identified by structured questionnaires, medical records, or linkage with a nationwide registry of cancer or death index or both. By the time we conducted this study, the majority of the women in these studies had passed through menopause. As most of the participants in these studies were European whites, we restricted analyses to women of European ancestry. We excluded subjects with missing information on age. Two studies (NHS and WGHS) provided information on the ER status of the breast tumors for a subset of the cases. This information was extracted from medical records. Each study was approved by the relevant local institutional review boards.

Table 1 List of participating studies and number of case and control subjects

Genotype data

We analyzed genotypes for 19 and 17 independent SNPs with reported associations with age at menarche and age at natural menopause, respectively, in the ReproGen Consortium, in which all SNPs achieved genome-wide significance in the meta-analysis of each trait (combined stage 1 and replication P value of less than 1 × 10-8) [10, 11]. None of these SNPs has been reported to be associated with breast cancer risk in previous GWAS and this is likely because of the very stringent P value threshold used to declare genome-wide significance (usually, P values were less than at least 1 × 10-7). As positive controls, 10 SNPs with consistently reported associations with breast cancer as shown in recent GWAS were included [17–19]. All 46 SNPs are listed in Table S1 of Additional file 2. Genotypes used in this analysis have been previously described [10, 11]. Complete genotype data from a total of up to 3,683 cases and 34,174 control subjects were available for analysis after the exclusions described in the 'Study population' section.

Breast cancer risk factors

The six studies from the ReproGen Consortium provided information on one or more of the following risk factors for breast cancer: age (continuous, at study entry or diagnosis), age at menarche (continuous, between 9 and 17 years), age at natural menopause (continuous, between 40 and 60 years), age at first live birth (less than 20, 20 to 24, 25 to 29 or no birth, at least 30 years), family history of breast cancer in first-degree relatives (yes/no), alcohol consumption (less than 5, 5 to 15, 15 to 30, at least 30 g/day), parity (0, 1 to 2, at least 3), menopausal hormone therapy (ever/never), oral contraceptive (OC) use (ever/never), and adult body mass index (BMI) (continuous).

Genetic risk score computation

The GRS was calculated on the basis of the 19 and 17 independent SNPs identified in previous studies as being associated with age at menarche and age at natural menopause, respectively [10, 11]. As a younger age at menarche and an older age at menopause are independently associated with an elevated breast cancer risk, we computed separate GRSs for a younger age at menarche and an older age at natural menopause. The risk allele was defined as an allele that was associated with a younger age at menarche or an older age at natural menopause. Two methods were used to determine the GRS: a simple count method (count GRS) and a weighted method (weighted GRS). We assumed an additive genetic model for each SNP, applying a linear weighting of 0, 1, or 2 to genotypes containing 0, 1, or 2 risk alleles, respectively. The count method assumes that each SNP contributes equally to the risk of breast cancer. The count GRS was calculated by simply summing the number of risk alleles of each SNP. For the weighted GRS, each SNP was weighted by β-coefficients obtained from the replication studies of recent meta-analyses of two traits [10, 11]. The weighted GRS was calculated by multiplying each β-coefficient by the number of corresponding risk alleles (0, 1, or 2) and then summing the products. To simplify interpretation and facilitate comparison with the count GRS, the weighted GRS was further divided by twice the sum of the β-coefficients and then multiplied by the total number of risk alleles. To provide a positive control and also to control for potential confounding by known breast cancer-associated genetic variants, a count GRS was computed on the basis of the 10 SNPs with consistently reported associations with breast cancer [19]: rs2981582, rs3803662, rs11249433, rs7716600, rs13387042, rs889312, rs13281615, rs999737, rs3817198, and rs1045485.

Statistical analysis

In each of the six studies, we performed logistic regression to evaluate the association with breast cancer for each of the 46 candidate SNPs, assuming an additive genetic model. Logistic regression was also used to analyze the association between GRS and breast cancer by including both GRSs for age at menarche and age at natural menopause in the model as the main effects. The GRSs were modeled as continuous variables or categorized into quintiles, and the cutoff points for quintiles were based on the WGHS population, which is the largest prospective cohort population among all participating studies. This approach was applied to each of the six participating studies. Odds ratios (ORs) and 95% confidence intervals (CIs) were estimated from logistic regression. To control for potential confounding by population stratification, we adjusted for the top principal components of genetic variation chosen for each study. We adjusted for age in the main model. To examine whether the genetic association of each of the candidate SNPs or GRSs with breast cancer is mediated through the onset of menarche or natural menopause, we then adjusted for self-reported age at menarche and age at natural menopause in the main model. Other conventional risk factors for breast cancer - including age at first live birth, family history of breast cancer in first-degree relatives, alcohol consumption, parity, menopausal hormone therapy, OC use, and adult BMI - were further included in the model to control for potential confounding in studies which had such data available. To examine whether these genetic associations differ by breast cancer subtypes, in each of the two studies that provided information on ER status, we then investigated the genetic association of each of the candidate SNPs or GRSs with breast cancer in subgroup analysis by ER histological status (positive or negative).

Forest plots were used to present study-specific ORs and 95% CIs. We then performed meta-analyses by using the fixed-effects model to estimate summary ORs from study-specific estimatesthat were weighted by the inverse of the variance of each study. As the meta-analyses restricted to prospective cohort studies or case-control studies yielded similar results, we present results from only the meta-analysis of all six participating studies. We also tested the heterogeneity of associations across studies as well as across different tumor subtypes by using the Q test [20].

All statistical analyses were performed by using SAS version 9.1 software (SAS Institute Inc., Cary, NC, USA). Power calculations were carried out by using Quanto (University of Southern California, Los Angeles, CA, USA). All P values were based on two-sided tests and were considered statistically significant if less than 0.05. Because SNPs were selected on the basis of an a priori hypothesis, adjustments for multiple comparison tests were not performed.

Results

The six participating studies contributed 3,683 breast cancer cases and 34,174 controls of self-reported white women of European ancestry (Table 1), all with available data on age and the 46 candidate SNPs, and at least one of the conventional risk factors considered. Of the 3,683 cases, about 52% were from the four prospective cohort studies (WGHS, RSI+II, FHS, and ARIC), about 30% were from the nested case-control study in NHS, and about 18% were from the population-based case-control study in SardiNIA. ER status was known for 2,087 cases in the NHS and the WGHS. On average, compared with the controls, the cases had a younger age at menarche and an older age at natural menopause. The expected associations with breast cancer were generally observed for the conventional risk factors across all of the studies (Table S2 of Additional file 3). The associations of the 46 candidate SNPs with age at menarche or age at natural menopause in the six studies were consistent with the original findings from the two meta-analyses [10, 11].

Table 2 shows the risk allele frequency and the corresponding per-risk-allele OR of breast cancer for each of the 46 candidate SNPs. The results are arranged in order of the strength of statistical significance (P value). The allele frequency for each SNP in the controls was similar to those reported for populations of European descent [21–23]. After adjusting for age and potential population stratification, we found that, among the 19 candidate SNPs for a younger age at menarche, two SNPs, rs1079866 and rs7821178, were significantly associated with breast cancer risk and had corresponding per-risk-allele ORs of 1.14 (96% CI = 1.05 to 1.24; P value = 0.003; P for heterogeneity = 0.37) and 1.08 (95% CI = 1.02 to 1.15; P value = 0.009; P for heterogeneity = 0.43), respectively. The SNP rs1079866 is located about 250 kb away from the INHBA gene on chromosome 7, whereas SNP rs7821178 is about 181 kb away from the PXMP3 gene (also known as PEX3) on chromosome 8. The strongest GWAS hit for age at menarche, rs7759938 at LIN28B on chromosome 6, was not found to be associated with breast cancer risk (P value = 0.60). Of the 17 candidate SNPs associated with an older age at natural menopause, one SNP, rs2517388, was significantly associated with breast cancer risk with a per-risk-allele OR of 1.10 (95% CI = 1.01 to 1.20; P value = 0.023; P for heterogeneity = 0.08). This SNP is an intronic SNP in the ASH2L gene on chromosome 8. The study-specific and summary ORs for the three associated SNPs are shown in Figure 1. Further adjustment for conventional risk factors - including age at menarche, age at natural menopause, age at first live birth, family history of breast cancer in first-degree relatives, alcohol consumption, parity, menopausal hormone therapy, OC use, and adult BMI - did not change the results substantially. For candidate loci for age at menarche and age at natural menopause, the findings did not differ materially when we further adjusted for known breast cancer-associated SNPs.

Table 2 Association of candidate single-nucleotide polymorphism loci and the risk of breast cancer
Figure 1
figure 1

Forest plots for the three candidate loci (rs1079866, rs7821178, and rs2517388) in association with breast cancer risk. Per-risk-allele odds ratios (ORs) and 95% confidence intervals (CIs) were obtained from unconditional logistic regression in each study, and age and potential population stratification were adjusted for. The size of the box is inversely proportional to the standard error of the log OR estimate. P values for heterogeneity across studies are 0.37, 0.43, and 0.08, respectively. ARIC, Atherosclerosis Risk in Communities Study; FHS, Framingham Heart Study; NHS, Nurses' Health Study; RSI+II, Rotterdam Study I, II; SardiNIA, SardiNIA Breast Cancer Study; WGHS, Women's Genome Health Study.

To evaluate the combined effect of candidate SNPs on breast cancer risk, we calculated a GRS for each trait by using either a count GRS or a weighted GRS approach. The mean values of count and weighted GRSs were 20.41 and 20.03, respectively, for age at menarche and 16.21 and 14.32, respectively, for age at natural menopause (Table 3). Based on the count GRS for a younger age at menarche, the OR for breast cancer associated with each point scored, corresponding to 1 risk allele, was 1.01 (95% CI = 1.00 to 1.03) after age and potential population stratification were adjusted for. ORs did not increase linearly across quintiles of GRS for age at menarche (P for trend = 0.06). Compared with women in the lowest quintile, women in the fourth and fifth quintiles had ORs for breast cancer of 1.14 (95% CI = 1.01 to 1.28) and 1.13 (95% CI = 1.00 to 1.27), respectively. Results were similar when analyses were performed by using weighted GRS. Overall, we did not observe statistically significant associations between breast cancer risk and age at natural menopause when either count or weighted GRS was used.

Table 3 Association between genetic risk score and risk of breast cancer

In secondary analyses, we then determined whether the associations of the 46 candidate SNPs with breast cancer vary across tumor subtypes defined by ER status in the NHS and the WGHS (Table 4). For the two SNPs (rs1079866 and rs7821178) that had reported associations with age at menarche and that were associated with overall breast cancer risk, we found no statistically significant evidence that the associations differed across subtypes (P for heterogeneity = 0.31 and 0.66, respectively), although rs1079866 appeared to have a stronger association with ER+ tumors (per-allele OR = 1.26; 95% CI = 1.12 to 1.41) than with ER- tumors (per-allele OR = 1.11; 95% CI = 0.89 to 1.38). Of note, one SNP that had a reported association with age at menarche, rs17188434, had a significantly stronger association with ER- tumors (per-allele OR = 1.51; 95% CI = 1.15 to 1.98) than with ER+ tumors (per-allele OR = 1.08; 95% CI = 0.92 to 1.26; P for heterogeneity = 0.035). Another SNP that had a reported association with age at menarche, rs17268785, was associated with a decreased risk of ER- tumors (per-allele OR = 0.83; 95% CI = 0.68 to 1.00) but an increased risk of ER+ tumors (per-allele OR = 1.07; 95% CI = 0.96 to 1.19; P for heterogeneity = 0.023). For the SNP that had a reported association with age at natural menopause and that was associated with overall breast cancer risk, we observed a stronger association with ER- tumors (per-allele OR = 1.32; 95% CI = 1.09 to 1.61) than ER+ tumors (per-allele OR = 1.11; 95% CI = 1.00 to 1.24); however, the test for heterogeneity was not statistically significant (P for heterogeneity = 0.12). When the count GRS for age at menarche or age at natural menopause was applied to ER+ and ER- breast cancer separately, the trend in the OR for ER+ tumors was very similar to that for overall breast cancer. The ER- tumor data suggested a somewhat different pattern, although the statistical power was limited for this subtype (Figure 2).

Table 4 Association of candidate single-nucleotide polymorphism loci and risk of breast cancer by estrogen receptor status in the Nurses' Health Study and Women's Genome Health Study
Figure 2
figure 2

The associations between groups defined by quintiles of genetic risk scores (GRSs) and risk of breast cancer by estrogen receptor (ER) status in the Nurses' Health Study and the Women's Genome Health Study. (a) Count GRS for age at menarche. (b) Count GRS for age at natural menopause. (c) Count GRS for breast cancer-associated SNPs. CI, confidence interval.

Of the 10 candidate SNPs with consistently reported associations with breast cancer risk, five SNPs (rs11249433, rs3803662, rs2981582, rs13387042, and rs999737) appeared to have a stronger association with ER+ tumor than ER- tumors, and rs3803662 reached statistical significance (P for heterogeneity = 0.008) with per-risk-allele ORs of 1.26 (95% CI = 1.16 to 1.37) and 0.98 (95% CI = 0.82 to 1.16) for ER+ and ER- tumors, respectively (Table 4). Two breast cancer candidate SNPs, rs1045485 and rs3817198, did not show statistically significant associations with overall risk. However, rs1045485 appeared to have a stronger association with ER- tumors (per-allele OR = 1.20; 95% CI = 0.95 to 1.52) than ER+ tumors (per-allele OR = 1.06; 95% CI = 0.95 to 1.19; P for heterogeneity = 0.35), and rs3817198 was associated with a decreased risk of ER- tumors (per-allele OR = 0.88; 95% CI = 0.74 to 1.04) but an increased risk of ER+ tumors (per-allele OR = 1.06; 95% CI = 0.98 to 1.15; P for heterogeneity = 0.047).

In these analyses, we further confirmed statistically significant associations with breast cancer risk for 8 of the 10 candidate SNPs that were identified previously in published GWAS of breast cancer (most P values were less than 0.001) (Table 2). We did not observe a statistically significant association for either LSP1-rs3817198 or CASP8-rs1045485 (both with P values of 0.13) in our study, although the direction of the associations was consistent with that of previous reports [21, 22]. We also calculated, as a positive control, a count GRS based on these 10 SNPs. We found that each score point increase, corresponding to one-risk-allele increase, was significantly associated with an OR of 1.13 (95% CI = 1.11 to 1.15) for breast cancer (Table 3). Compared with women in the lowest quintile, women in the highest quintile had an OR for breast cancer of 1.89 (95% CI = 1.67 to 2.14). For this GRS, the trend in log odds was significantly steeper for ER+ than for ER- tumors (P for heterogeneity < 0.001), and the OR across quintiles was no longer monotonic in ER- tumors (Figure 2).

Discussion

In this large meta-analysis of six population-based studies, we investigated whether 19 loci linked with age at menarche and 17 loci linked with age at natural menopause were associated with breast cancer risk among up to 3,683 breast cancer cases and 34,174 controls. We found that two SNPs with reported associations with age at menarche and one SNP with a reported association with age at natural menopause were significantly associated with breast cancer risk. However, no statistically significant associations were found for GRSs that combined all 19 or 17 loci associated with each trait, although the association for age-at-menarche GRS was marginally statistically significant. We confirmed most of the candidate loci for breast cancer which were identified in previous GWAS. Some of these associations appeared to differ by tumor subtypes defined by ER status.

In our analyses, most of the candidate SNPs, including the strongest GWAS hit for age at menarche or age at natural menopause, were not found to be associated with breast cancer risk. This is not necessarily surprising given that age at menarche and age at natural menopause are relatively weak risk factors [2, 3], and all candidate SNPs collectively explain only a small portion of the variation of each trait [10, 11]. However, two candidate SNPs for age at menarche, rs1079866 and rs7821178, and one candidate SNP for age at natural menopause, rs2517388, were found to be associated with breast cancer risk. These associations were not attenuated after we further adjusted for self-reported age at menarche and age at natural menopause, suggesting these three genetic loci were associated with breast cancer risk independently of their associations with age at menarche or age at natural menopause. It is possible that these genetic loci have pleiotropic effects on reproductive timing as well as other biological processes leading to breast cancer, and the observed associations might be due largely to other biological consequences of these risk variants that do not manifest themselves as changes in age at menarche or age at natural menopause. Alternatively, it is also possible that the relatively crude assignment of these reproductive events to a single chronological year is not sufficiently accurate to capture the biological effect of these processes on breast cancer risk and the genetic variants contribute independent information on the underlying biological risk. The three candidate SNPs also contributed to breast cancer risk independently of the known susceptibility loci for breast cancer, as further adjustment for breast cancer loci did not materially alter the results.

We found statistically significant evidence of association with breast cancer for eight of the 10 breast cancer susceptibility loci examined: FGFR2-rs2981582, TNRC9-rs3803662, 1p-rs11249433, 5p-rs7716600, 2q35-rs13387042, MAP3K1-rs889312, 8q24-rs13281615, and RAD51L1- rs999737. The direction and magnitude of these associations were consistent with those of previous reports [17, 18, 22–25]. We did not observe a statistically significant association for either LSP1-rs3817198 or CASP8-rs1045485. However, these two SNPs had relatively small reported effects that our study might not have been able to detect. When the 10 candidate SNPs were combined by using a polygenic risk score, the relative risk for women in the highest quintile was about twice that in the lowest quintile, and this is in accordance with other published results [19, 26]. In this study, none of the 10 breast cancer susceptibility loci was significantly associated with age at menarche or age at natural menopause, and this is in line with a previous report [27].

Given that most of the candidate loci for age at menarche and age at natural menopause were not associated with breast cancer risk, it is not surprising that there were no statistically significant associations for the polygenetic risk scores that combined all candidate loci for each trait. To conduct a post hoc and exploratory analysis, we created a polygenetic risk score by including only the three candidate loci associated with either age at menarche or age at natural menopause and with breast cancer risk and found that each risk allele increment was associated with an approximately 17% increased risk for breast cancer. Women with four or more risk alleles had an approximately 60% increased risk for breast cancer in comparison with those with two risk alleles or less. When we further combined the three associated SNPs with the 10 breast cancer susceptibility loci to create a polygenetic risk score, each risk allele increment was associated with an approximately 18% increased risk for breast cancer. For women with 14 or more risk alleles (the highest quintile), the risk for breast cancer increased threefold in comparison with those with 10 or less (the lowest quintile). Because the former group constitutes approximately 20% of the study population, the GRS that combines the three candidate SNPs for age at menarche and age at natural menopause and the identified breast cancer susceptibility loci might be useful for identifying a subgroup of women with a high genetic risk for breast cancer. Further research is needed to confirm this finding.

It has been hypothesized that the risk of ER+ breast cancer is positively associated with a woman's cumulative lifetime exposure to endogenous ovarian hormones [28]. A younger age at menarche [12, 15, 29] and an older age at menopause [30] have been observed to be more consistently associated with ER+ than ER- tumors. In this report, we found that candidate loci for age at menarche and age at natural menopause may also be differentially associated with tumor subtypes defined by ER status. Of the three candidate loci that were found to be associated with overall breast cancer risk, rs1079866 was more strongly associated with ER+ tumors, rs7821178 was equally associated with both, whereas rs2517388 was more strongly associated with ER- tumors, although differences were not statistically significant. Importantly, two candidate loci for age at menarche, rs17188434 and rs17268785, had significantly different associations with ER+ and ER- tumors. Whereas both SNPs were not significantly associated with overall and ER+ breast cancer, the former showed a statistically significant positive association with ER- tumors, whereas the latter showed a statistically significant inverse association with ER- tumor. These findings provide further support for the notion that ER+ and ER- tumors are the result of different etiologic pathways [31].

Although common genetic variants that influence the intermediate phenotypes or risk factors have been hypothesized to subsequently affect disease risk, few studies have assessed the association between these genetic variants and disease risk or, furthermore, whether these associations are mediated through the intermediate phenotypes. Chen and colleagues [32] investigated obesity-linked genetic variants in relation to breast cancer risk but found no statistically significant association. To our knowledge, ours is the first study to evaluate the associations of candidate loci for age at menarche and age at natural menopause with breast cancer risk. One of the strengths of our study is the relatively large combined sample size achieved through international collaboration. We had adequate statistical power (80%) to detect an OR of 1.12 for SNPs with a minor allele frequency (MAF) of 0.10 and an OR of 1.09 for SNPs with an MAF of 0.20. However, our analysis of ER+ tumors was less adequately powered, as the ER status was not available for all cases, and the study had limited statistical power for ER- tumors. One limitation in our study is the multiple comparisons that could lead to false-positive results. Although none of the candidate SNPs with a reported association with age at menarche or age at natural menopause survived Bonferroni correction in the test of breast cancer association, this correction is considered to be overly conservative given that the candidates were chosen on the basis of promising hypotheses. Another potential limitation of our study comes from differences in the study population and designs and methods of collecting risk factors and genetic marker data across studies. However, the findings were generally consistent across studies, arguing for the robustness of our results. Finally, as our analyses were restricted to women of European ancestry, results from this study may not be generalizable to other ethnic groups.

Conclusions

In summary, in this large analysis of the association of several novel candidate loci for age at menarche and age at natural menopause with breast cancer risk, we observed that three loci - two for age at menarche and one for age at natural menopause - were significantly associated with breast cancer risk independently of their associations with each trait and independently of known breast cancer susceptibility loci. These associations may differ by tumor subtypes defined by ER status. A combination of all 19 loci associated with age at menarche or 17 loci associated with age at natural menopause did not appear to be helpful for identifying a high-risk subgroup for breast cancer.

Abbreviations

ARIC:

Atherosclerosis Risk in Communities Study

BMI:

body mass index

CI:

confidence interval

ER:

estrogen receptor

FHS:

Framingham Heart Study

GRS:

genetic risk score

GWAS:

genome-wide association studies

MAF:

minor allele frequency

NHS:

Nurses' Health Study

OC:

oral contraceptive

OR:

odds ratio

RSI+II:

Rotterdam Study I and II

SardiNIA:

SardiNIA Breast Cancer Study

SNP:

single-nucleotide polymorphism

WGHS:

Women's Genome Health Study.

References

  1. Kelsey JL, Horn-Ross PL: Breast cancer: magnitude of the problem and descriptive epidemiology. Epidemiol Rev. 1993, 15: 7-16.

    CAS  PubMed  Google Scholar 

  2. Hsieh CC, Trichopoulos D, Katsouyanni K, Yuasa S: Age at menarche, age at menopause, height and obesity as risk factors for breast cancer: associations and interactions in an international case-control study. Int J Cancer. 1990, 46: 796-800. 10.1002/ijc.2910460508.

    Article  CAS  PubMed  Google Scholar 

  3. Breast cancer and hormone replacement therapy: collaborative reanalysis of data from 51 epidemiological studies of 52,705 women with breast cancer and 108,411 women without breast cancer. Collaborative Group on Hormonal Factors in Breast Cancer. Lancet. 1997, 350: 1047-1059.

  4. Pike MC, Spicer DV, Dahmoush L, Press MF: Estrogens, progestogens, normal breast cell proliferation, and breast cancer risk. Epidemiol Rev. 1993, 15: 17-35.

    CAS  PubMed  Google Scholar 

  5. He C, Kraft P, Chen C, Buring JE, Pare G, Hankinson SE, Chanock SJ, Ridker PM, Hunter DJ, Chasman DI: Genome-wide association studies identify loci associated with age at menarche and age at natural menopause. Nat Genet. 2009, 41: 724-728. 10.1038/ng.385.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Ong KK, Elks CE, Li S, Zhao JH, Luan J, Andersen LB, Bingham SA, Brage S, Smith GD, Ekelund U, Gillson CJ, Glaser B, Golding J, Hardy R, Khaw KT, Kuh D, Luben R, Marcus M, McGeehin MA, Ness AR, Northstone K, Ring SM, Rubin C, Sims MA, Song K, Strachan DP, Vollenweider P, Waeber G, Waterworth DM, Wong A, et al: Genetic variation in LIN28B is associated with the timing of puberty. Nat Genet. 2009, 41: 729-733. 10.1038/ng.382.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Perry JR, Stolk L, Franceschini N, Lunetta KL, Zhai G, McArdle PF, Smith AV, Aspelund T, Bandinelli S, Boerwinkle E, Cherkas L, Eiriksdottir G, Estrada K, Ferrucci L, Folsom AR, Garcia M, Gudnason V, Hofman A, Karasik D, Kiel DP, Launer LJ, van Meurs J, Nalls MA, Rivadeneira F, Shuldiner AR, Singleton A, Soranzo N, Tanaka T, Visser JA, Weedon MN, et al: Meta-analysis of genome-wide association data identifies two loci influencing age at menarche. Nat Genet. 2009, 41: 648-650. 10.1038/ng.386.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Sulem P, Gudbjartsson DF, Rafnar T, Holm H, Olafsdottir EJ, Olafsdottir GH, Jonsson T, Alexandersen P, Feenstra B, Boyd HA, Aben KK, Verbeek AL, Roeleveld N, Jonasdottir A, Styrkarsdottir U, Steinthorsdottir V, Karason A, Stacey SN, Gudmundsson J, Jakobsdottir M, Thorleifsson G, Hardarson G, Gulcher J, Kong A, Kiemeney LA, Melbye M, Christiansen C, Tryggvadottir L, Thorsteinsdottir U, Stefansson K: Genome-wide association study identifies sequence variants on 6q21 associated with age at menarche. Nat Genet. 2009, 41: 734-738. 10.1038/ng.383.

    Article  CAS  PubMed  Google Scholar 

  9. Stolk L, Zhai G, van Meurs JB, Verbiest MM, Visser JA, Estrada K, Rivadeneira F, Williams FM, Cherkas L, Deloukas P, Soranzo N, de Keyzer JJ, Pop VJ, Lips P, Lebrun CE, van der Schouw YT, Grobbee DE, Witteman J, Hofman A, Pols HA, Laven JS, Spector TD, Uitterlinden AG: Loci at chromosomes 13, 19 and 20 influence age at natural menopause. Nat Genet. 2009, 41: 645-647. 10.1038/ng.387.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Elks CE, Perry JR, Sulem P, Chasman DI, Franceschini N, He C, Lunetta KL, Visser JA, Byrne EM, Cousminer DL, Gudbjartsson DF, Esko T, Feenstra B, Hottenga JJ, Koller DL, Kutalik Z, Lin P, Mangino M, Marongiu M, McArdle PF, Smith AV, Stolk L, van Wingerden SH, Zhao JH, Albrecht E, Corre T, Ingelsson E, Hayward C, Magnusson PK, Smith EN, et al: Thirty new loci for age at menarche identified by a meta-analysis of genome-wide association studies. Nat Genet. 2010, 42: 1077-1085. 10.1038/ng.714.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Stolk L, Perry JR, Chasman DI, He C, Mangino M, Sulem P, Barbalic M, Broer L, Byrne EM, Ernst F, Esko T, Franceschini N, Gudbjartsson DF, Hottenga JJ, Kraft P, McArdle PF, Porcu E, Shin SY, Smith AV, van Wingerden S, Zhai G, Zhuang WV, Albrecht E, Alizadeh BZ, Aspelund T, Bandinelli S, Lauc LB, Beckmann JS, Boban M, Boerwinkle E, et al: Meta-analyses identify 13 loci associated with age at menopause and highlight DNA repair and immune pathways. Nat Genet. 2012, 44: 260-268. 10.1038/ng.1051.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. 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: 1558-1568.

    CAS  PubMed  Google Scholar 

  13. Kotsopoulos J, Chen WY, Gates MA, Tworoger SS, Hankinson SE, Rosner BA: Risk factors for ductal and lobular breast cancer: results from the nurses' health study. Breast Cancer Res. 2010, 12: R106-10.1186/bcr2790.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Li CI, Malone KE, Daling JR, Potter JD, Bernstein L, Marchbanks PA, Strom BL, Simon MS, Press MF, Ursin G, Burkman RT, Folger SG, Norman S, McDonald JA, Spirtas R: Timing of menarche and first full-term birth in relation to breast cancer risk. Am J Epidemiol. 2008, 167: 230-239.

    Article  PubMed  Google Scholar 

  15. 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: R43-10.1186/bcr1525.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Reeves GK, Pirie K, Green J, Bull D, Beral V: Reproductive factors and specific histological types of breast cancer: prospective study and meta-analysis. Br J Cancer. 2009, 100: 538-544. 10.1038/sj.bjc.6604853.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Thomas G, Jacobs KB, Kraft P, Yeager M, Wacholder S, Cox DG, Hankinson SE, Hutchinson A, Wang Z, Yu K, Chatterjee N, Garcia-Closas M, Gonzalez-Bosquet J, Prokunina-Olsson L, Orr N, Willett WC, Colditz GA, Ziegler RG, Berg CD, Buys SS, McCarty CA, Feigelson HS, Calle EE, Thun MJ, Diver R, Prentice R, Jackson R, Kooperberg C, Chlebowski R, Lissowska J, et al: A multistage genome-wide association study in breast cancer identifies two new risk alleles at 1p11.2 and 14q24.1 (RAD51L1). Nat Genet. 2009, 41: 579-584. 10.1038/ng.353.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Pharoah PD, Antoniou AC, Easton DF, Ponder BA: Polygenes, risk prediction, and targeted prevention of breast cancer. N Engl J Med. 2008, 358: 2796-2803. 10.1056/NEJMsa0708739.

    Article  CAS  PubMed  Google Scholar 

  19. Wacholder S, Hartge P, Prentice R, Garcia-Closas M, Feigelson HS, Diver WR, Thun MJ, Cox DG, Hankinson SE, Kraft P, Rosner B, Berg CD, Brinton LA, Lissowska J, Sherman ME, Chlebowski R, Kooperberg C, Jackson RD, Buckman DW, Hui P, Pfeiffer R, Jacobs KB, Thomas GD, Hoover RN, Gail MH, Chanock SJ, Hunter DJ: Performance of common genetic variants in breast-cancer risk models. N Engl J Med. 2010, 362: 986-993. 10.1056/NEJMoa0907727.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Higgins JP, Thompson SG: Quantifying heterogeneity in a meta-analysis. Stat Med. 2002, 21: 1539-1558. 10.1002/sim.1186.

    Article  PubMed  Google Scholar 

  21. Cox A, Dunning AM, Garcia-Closas M, Balasubramanian S, Reed MW, Pooley KA, Scollen S, Baynes C, Ponder BA, Chanock S, Lissowska J, Brinton L, Peplonska B, Southey MC, Hopper JL, McCredie MR, Giles GG, Fletcher O, Johnson N, dos Santos Silva I, Gibson L, Bojesen SE, Nordestgaard BG, Axelsson CK, Torres D, Hamann U, Justenhoven C, Brauch H, Chang-Claude J, Kropp S, et al: A common coding variant in CASP8 is associated with breast cancer risk. Nat Genet. 2007, 39: 352-358. 10.1038/ng1981.

    Article  CAS  PubMed  Google Scholar 

  22. Easton DF, Pooley KA, Dunning AM, Pharoah PD, Thompson D, Ballinger DG, Struewing JP, Morrison J, Field H, Luben R, Wareham N, Ahmed S, Healey CS, Bowman R, Meyer KB, Haiman CA, Kolonel LK, Henderson BE, Le Marchand L, Brennan P, Sangrajrang S, Gaborieau V, Odefrey F, Shen CY, Wu PE, Wang HC, Eccles D, Evans DG, Peto J, Fletcher O, et al: Genome-wide association study identifies novel breast cancer susceptibility loci. Nature. 2007, 447: 1087-1093. 10.1038/nature05887.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Stacey SN, Manolescu A, Sulem P, Rafnar T, Gudmundsson J, Gudjonsson SA, Masson G, Jakobsdottir M, Thorlacius S, Helgason A, Aben KK, Strobbe LJ, Albers-Akkers MT, Swinkels DW, Henderson BE, Kolonel LN, Le Marchand L, Millastre E, Andres R, Godino J, Garcia-Prats MD, Polo E, Tres A, Mouy M, Saemundsdottir J, Backman VM, Gudmundsson L, Kristjansson K, Bergthorsson JT, Kostic J, et al: Common variants on chromosomes 2q35 and 16q12 confer susceptibility to estrogen receptor-positive breast cancer. Nat Genet. 2007, 39: 865-869. 10.1038/ng2064.

    Article  CAS  PubMed  Google Scholar 

  24. Hunter DJ, Kraft P, Jacobs KB, Cox DG, Yeager M, Hankinson SE, Wacholder S, Wang Z, Welch R, Hutchinson A, Wang J, Yu K, Chatterjee N, Orr N, Willett WC, Colditz GA, Ziegler RG, Berg CD, Buys SS, McCarty CA, Feigelson HS, Calle EE, Thun MJ, Hayes RB, Tucker M, Gerhard DS, Fraumeni JF, Hoover RN, Thomas G, Chanock SJ: A genome-wide association study identifies alleles in FGFR2 associated with risk of sporadic postmenopausal breast cancer. Nat Genet. 2007, 39: 870-874. 10.1038/ng2075.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Stacey SN, Manolescu A, Sulem P, Thorlacius S, Gudjonsson SA, Jonsson GF, Jakobsdottir M, Bergthorsson JT, Gudmundsson J, Aben KK, Strobbe LJ, Swinkels DW, van Engelenburg KC, Henderson BE, Kolonel LN, Le Marchand L, Millastre E, Andres R, Saez B, Lambea J, Godino J, Polo E, Tres A, Picelli S, Rantala J, Margolin S, Jonsson T, Sigurdsson H, Jonsdottir T, Hrafnkelsson J, et al: Common variants on chromosome 5p12 confer susceptibility to estrogen receptor-positive breast cancer. Nat Genet. 2008, 40: 703-706. 10.1038/ng.131.

    Article  CAS  PubMed  Google Scholar 

  26. Reeves GK, Travis RC, Green J, Bull D, Tipper S, Baker K, Beral V, Peto R, Bell J, Zelenika D, Lathrop M: Incidence of breast cancer and its subtypes in relation to individual and multiple low-penetrance genetic susceptibility loci. JAMA. 2010, 304: 426-434. 10.1001/jama.2010.1042.

    Article  CAS  PubMed  Google Scholar 

  27. Travis RC, Reeves GK, Green J, Bull D, Tipper SJ, Baker K, Beral V, Peto R, Bell J, Zelenika D, Lathrop M: Gene-environment interactions in 7610 women with breast cancer: prospective evidence from the Million Women Study. Lancet. 2010, 375: 2143-2151. 10.1016/S0140-6736(10)60636-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Clavel-Chapelon F: Differential effects of reproductive factors on the risk of pre- and postmenopausal breast cancer. Results from a large cohort of French women. Br J Cancer. 2002, 86: 723-727. 10.1038/sj.bjc.6600124.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. 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: 250-263. 10.1093/jnci/djq526.

    Article  PubMed  Google Scholar 

  30. 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: 470-477. 10.1093/jnci/djr030.

    Article  PubMed  PubMed Central  Google Scholar 

  31. 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: 240-245. 10.1186/bcr938.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Chen F, Wilkens LR, Monroe KR, Stram DO, Kolonel LN, Henderson BE, Le Marchand L, Haiman CA: No association of risk variants for diabetes and obesity with breast cancer: the Multiethnic Cohort and PAGE studies. Cancer Epidemiol Biomarkers Prev. 2011, 20: 1039-1042. 10.1158/1055-9965.EPI-11-0135.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The NHS breast cancer GWAS scan was performed as part of the Cancer Genetic Markers of Susceptibility initiative of the National Cancer Institute (NCI). We particularly acknowledge the contributions of Robert Hoover, Amy Hutchinson, Kevin Jacobs, and Gilles Thomas. The present research is supported by CA 40356 and U01-CA98233 from the NCI. We acknowledge the study participants in the NHS for their contribution in making this study possible. The WGHS is supported by HL 043851 and HL69757 from the National Heart, Lung, and Blood Institute (NHLBI) and CA 047988 from the NCI, the Donald W Reynolds Foundation, and the Fondation Leducq, and collaborative scientific support and funding for genotyping were provided by Amgen (Thousand Oaks, CA, USA). We thank Joseph Miletich and Alex Parker and the technical staff at Amgen for their collaboration and scientific support in performing the genotyping for the WGHS. The SardiNIA Breast Cancer Study was funded by National Institutes of Health (NIH) contract NO1-AG-1-2109 from the National Institute of Aging (NIA) to the SardiNIA team (principal investigator: Manuela Uda). We thank Grazia Palomba, Antonio Cossu, Francesco Tanda, Mario Budroni, and the other members of the Sardinian Translational Oncology Group for their contribution in collecting Sardinian patients with breast cancer. The generation and management of GWAS genotype data for the Rotterdam Study are supported by The Netherlands Organisation for Scientific Research (NWO) Investments (175.010.2005.011 and 911-03-012). This study is funded by the Research Institute for Diseases in the Elderly (RIDE2) (014-93-015) and The Netherlands Genomics Initiative/NWO project 050-060-810. We thank Pascal Arp, Mila Jhamai, Marijn Verkerk, Lizbeth Herrera, and Marjolein Peters for their help in creating the GWAS database and Karol Estrada and Maksim V Struchalin for their support in creation and analysis of imputed data. The Rotterdam Study is funded by the Erasmus Medical Center and Erasmus University, Rotterdam; the Netherlands Organisation for Health Research and Development (ZonMw); the Research Institute for Diseases in the Elderly (RIDE); the Ministry of Education, Culture and Science; the Ministry for Health, Welfare and Sports; the European Commission (DG XII); and the Municipality of Rotterdam. The authors are grateful to the study participants, the staff from the Rotterdam Study, and the participating general practitioners and pharmacists. The authors acknowledge the essential role of the Cohorts for Heart and Aging Research in Genome Epidemiology (CHARGE) Consortium in the development and support of this manuscript. CHARGE members consist of The Netherlands' Rotterdam Study, the FHS, the Cardiovascular Health Study, the NHLBI's ARIC study, and the NIA's Iceland Age, Gene/Environment Susceptibility (AGES) study. The FHS phenotype-genotype analyses were supported in part by the NIA (Genetics of Reproductive Life Period and Health Outcomes, R21AG032598; JMM, KLL). The FHS of the NHLBI of the NIH and Boston University School of Medicine was supported by NHLBI's FHS contract N01-HC-25195 and its contract with Affymetrix, Inc. (Santa Clara, CA, USA) for genotyping services (contract N02-HL-6-4278). Analyses reflect intellectual input and resource development from the FHS investigators participating in the SNP Health Association Resource (SHARe) project. A portion of this research was conducted by using the Linux Cluster for Genetic Analysis (LinGA-II) funded by the Robert Dawson Evans Endowment of the Department of Medicine at Boston University School of Medicine and Boston Medical Center. The ARIC study is carried out as a collaborative study supported by NHLBI contracts N01-HC-55015, N01-HC-55016, N01-HC-55018, N01-HC-55019, N01-HC-55020, N01-HC-55021, N01-HC-55022, R01HL087641, R01HL59367, R01HL086694, and RC2 HL102419; National Human Genome Research Institute contract U01HG004402; and NIH contract HHSN268200625226C. The authors thank the staff and participants of the ARIC study for their important contributions. Infrastructure was partly supported by grant UL1RR025005, a component of the NIH, and NIH Roadmap for Medical Research.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Chunyan He.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

CH performed the meta-analysis, wrote the manuscript, helped to conceive and design the experiments, helped to perform the primary statistical analyses in each study, shared responsibility for the interpretation of results and critical revision of the manuscript, and contributed to a critical revision of the manuscript for important intellectual content. DIC helped to conceive and design the experiments, helped to perform the primary statistical analyses in each study, shared responsibility for the interpretation of results and critical revision of the manuscript, and contributed to a critical revision of the manuscript for important intellectual content. DJH helped to conceive and design the experiments, contributed to a critical revision of the manuscript for important intellectual content, and participated in the original design, subject recruitment, acquisition of data, biospecimen collection for the studies, and the genotyping and quality control of genotype and other data. JD, SJH, RR, and SS helped to perform the primary statistical analyses in each study and shared responsibility for the interpretation of results and critical revision of the manuscript. SJC, LC, EWD, JMM, PMR, and BHS contributed to a critical revision of the manuscript for important intellectual content and participated in the original design, subject recruitment, acquisition of data, biospecimen collection for the studies, and the genotyping and quality control of genotype and other data. JEB, LF-R, NF, SEH, AH, KLL, GP, EP, FR, LMR, GLS, LS, and AGU participated in the original design, subject recruitment, acquisition of data, biospecimen collection for the studies, and the genotyping and quality control of genotype and other data. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Supplementary methods for study population. (DOCX 51 KB)

13058_2011_2958_MOESM2_ESM.DOCX

Additional file 2: Table S1: Information on the 46 candidate SNP loci identified in previous genome-wide association studies for age at menarche, age at natural menopause and breast cancer. (DOCX 48 KB)

13058_2011_2958_MOESM3_ESM.DOCX

Additional file 3: Table S2: Characteristics of non-genetic risk factors for breast cancer in each participating study. (DOCX 40 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Reprints and permissions

About this article

Cite this article

He, C., Chasman, D.I., Dreyfus, J. et al. Reproductive aging-associated common genetic variants and the risk of breast cancer. Breast Cancer Res 14, R54 (2012). https://doi.org/10.1186/bcr3155

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/bcr3155

Keywords