Skip to main content

Genetic variation at CYP3A is associated with age at menarche and breast cancer risk: a case-control study

Abstract

Introduction

We have previously shown that a tag single nucleotide polymorphism (rs10235235), which maps to the CYP3A locus (7q22.1), was associated with a reduction in premenopausal urinary estrone glucuronide levels and a modest reduction in risk of breast cancer in women age ≤50 years.

Methods

We further investigated the association of rs10235235 with breast cancer risk in a large case control study of 47,346 cases and 47,570 controls from 52 studies participating in the Breast Cancer Association Consortium. Genotyping of rs10235235 was conducted using a custom Illumina Infinium array. Stratified analyses were conducted to determine whether this association was modified by age at diagnosis, ethnicity, age at menarche or tumor characteristics.

Results

We confirmed the association of rs10235235 with breast cancer risk for women of European ancestry but found no evidence that this association differed with age at diagnosis. Heterozygote and homozygote odds ratios (ORs) were OR = 0.98 (95% CI 0.94, 1.01; P = 0.2) and OR = 0.80 (95% CI 0.69, 0.93; P = 0.004), respectively (Ptrend = 0.02). There was no evidence of effect modification by tumor characteristics. rs10235235 was, however, associated with age at menarche in controls (Ptrend = 0.005) but not cases (Ptrend = 0.97). Consequently the association between rs10235235 and breast cancer risk differed according to age at menarche (Phet = 0.02); the rare allele of rs10235235 was associated with a reduction in breast cancer risk for women who had their menarche age ≥15 years (ORhet = 0.84, 95% CI 0.75, 0.94; ORhom = 0.81, 95% CI 0.51, 1.30; Ptrend = 0.002) but not for those who had their menarche age ≤11 years (ORhet = 1.06, 95% CI 0.95, 1.19, ORhom = 1.07, 95% CI 0.67, 1.72; Ptrend = 0.29).

Conclusions

To our knowledge rs10235235 is the first single nucleotide polymorphism to be associated with both breast cancer risk and age at menarche consistent with the well-documented association between later age at menarche and a reduction in breast cancer risk. These associations are likely mediated via an effect on circulating hormone levels.

Introduction

Family history is a well-established risk factor for breast cancer. First-degree relatives of women with breast cancer have an approximately twofold increased risk of developing the disease relative to the general population [1]. Twin studies are consistent with this familial clustering having, at least in part, a genetic origin [2, 3]. Mutations in high-risk susceptibility genes (mainly BRCA1 and BRCA2) explain most large multiple-case families, but account for only 15 to 20% of the excess familial risk [4]. Genome-wide association studies [5, 6] have identified more than 70 common variants that are associated with breast cancer susceptibility but they account for only another approximately 15% of the excess familial risk. The so-called ‘missing heritability’ may be explained by common variants with very small effects and/or by rarer variants with larger effects, neither of which can be identified by current genome-wide association studies. A statistically efficient alternative is to increase power by trying to identify variants associated with known quantitative phenotypic markers of susceptibility to breast cancer [7], and then to test them for association with breast cancer risk. This approach might also improve our understanding of the biological mechanisms involved in breast cancer pathogenesis.

Endogenous sex hormones are well-established risk factors for breast cancer in postmenopausal women [8]; the evidence in premenopausal women is less consistent, with some, but not all, studies suggesting an association between higher circulating levels of estrogens and increased breast cancer risk [917]. Genetic factors influence the levels of endogenous sex hormones [18] and therefore single nucleotide polymorphisms (SNPs) in genes regulating these hormonal pathways are good candidates for being breast cancer predisposition variants. We have previously studied 642 SNPs tagging 42 genes that might influence sex hormone levels in 729 healthy premenopausal women of European ancestry in relation to cyclic variations in oestrogen levels during the menstrual cycle. We found that the minor allele of rs10273424, which maps 50 kb 3′ to CYP3A5, was associated with a reduction of 22% (95% confidence interval (CI) = –28%, –15%; P = 10-9) in levels of urinary oestrone glucuronide, a metabolite that is highly correlated with serum oestradiol levels [19]. Analysis of 10,551 breast cancer cases and 17,535 controls of European ancestry demonstrated that the minor allele of rs10235235, a proxy for rs10273424 (r2 = 1.0), was also associated with a weak reduction in breast cancer risk but only in women aged 50 years or younger at diagnosis (odds ratio (OR) = 0.91, 95% CI = 0.83, 0.99; P = 0.03) [19].

The aim of the present study was to further investigate an association between rs10235235 and breast cancer risk using a much larger set of subjects – the Breast Cancer Association Consortium (BCAC) – comprising data from 49 additional studies, and to assess whether there was evidence of effect modification by age at diagnosis, ethnicity, age at menarche or tumour characteristics.

Materials and methods

Sample selection

Samples for the case–control analyses were drawn from 52 studies participating in the BCAC: 41 studies from populations of predominantly European ancestry, nine studies of Asian ancestry and two studies of African-American ancestry. The majority were population-based or hospital-based case–control studies, but some studies were nested in cohorts, selected samples by age, oversampled for cases with a family history or selected samples on the basis of tumour characteristics (Table S1 in Additional file 1). Studies provided ~2% of samples in duplicate for quality control purposes (see below). Study subjects were recruited on protocols approved by the Institutional Review Boards at each participating institution, and all subjects provided written informed consent (Additional file 2).

Genotyping and post-genotyping quality control

Genotyping for rs10235235 was carried out as part of a collaboration between the BCAC and three other consortia (the Collaborative Oncological Gene-environment Study (COGS)). Full details of SNP selection, array design, genotyping and post-genotyping quality control have been published [5]. Briefly, three categories of SNPs were chosen for inclusion in the array: SNPs selected on the basis of pooled genome-wide association study data; SNPs selected for the fine-mapping of published risk loci; and candidate SNPs selected on the basis of previous analyses or specific hypotheses. rs10235235 was a candidate SNP selected on the basis of our previous analyses [19].

For the COGS project overall, genotyping of 211,155 SNPs in 114,225 samples was conducted using a custom Illumina Infinium array (iCOGS; Illumina, San Diego, CA, USA) in four centres. Genotypes were called using Illumina’s proprietary GenCall algorithm. Standard quality control measures were applied across all SNPs and all samples genotyped as part of the COGS project. Samples were excluded for any of the following reasons: genotypically not female XX (XY, XXY or XO, n = 298); overall call rate <95% (n = 1,656); low or high heterozygosity (P < 10-6, separately for individuals of European, Asian and African-American ancestry, n = 670); individuals not concordant with previous genotyping within the BCAC (n = 702); individuals where genotypes for the duplicate sample appeared to be from a different individual (n = 42); cryptic duplicates within studies where the phenotypic data indicated that the individuals were different, or between studies where genotype data indicated samples were duplicates (n = 485); first-degree relatives (n = 1,981); phenotypic exclusions (n = 527); or concordant replicates (n = 2,629).

Ethnic outliers were identified by multidimensional scaling, combining the iCOGS array data with the three Hapmap2 populations, based on a subset of 37,000 uncorrelated markers that passed quality control (including ~1,000 selected as ancestry informative markers). Most studies were predominantly of a single ancestry (European or Asian), and women with >15% minority ancestry, based on the first two components, were excluded (n = 1,244). Two studies from Singapore (SGBCC) and Malaysia (MYBRCA; see Table S1 in Additional file 1 for all full study names) contained a substantial fraction of women of mixed European/Asian ancestry (probably of South Asian ancestry). For these studies, no exclusions for ethnic outliers were made, but principal components analysis (see below) was used to adjust for inflation in these studies. Similarly, for the two African-American studies (NBHS and SCCS), no exclusions for ethnic outliers were made.

Principal component analyses were carried out separately for the European, Asian and African-American subgroups, based on a subset of 37,000 uncorrelated SNPs. For the analyses of European subjects, we included the first six principal components as covariates, together with a seventh component derived specific to one study (LMBC) for which there was substantial inflation not accounted for by the components derived from the analysis of all studies. Addition of further principal components did not reduce inflation further. Two principal components were included for the studies conducted in Asian populations and two principal components were included for the African-American studies.

For the main analyses of rs10235235 and breast cancer risk, we excluded women from three studies (BBCS, BIGGS and UKBGS) that were genotyped in the hypothesis-generating study (n = 5,452) [19] and women with non-invasive cancers (ductal carcinoma in situ/lobular carcinoma in situ, n = 2,663) or cancers of uncertain status (n = 960)). After exclusions there were 47,346 invasive breast cancer case samples and 47,570 control samples from 49 studies (38 from populations of predominantly European ancestry, nine Asian and two African-American) used in the analysis (Tables S1 and S2 in Additional file 1). After quality control exclusions (above) the call rate for rs10235235 was 100% (one no call in 94,916 samples), and for the controls there was no evidence of deviation from Hardy–Weinberg equilibrium in any of the contributing studies (Table S2 in Additional file 1).

We did not test for an association between rs10235235 and age at menarche in our hypothesis-generating study [19]. Therefore, to maximise our power to detect an association, we included menarche data from BBCS cases (n = 2,508) and controls (n = 1,650) and from UKBGS cases (n = 3,388) and controls (n = 4,081) in this analysis. Age at menarche was not available for samples from BIGGS. Full details of genotyping of rs10235235 in BBCS and UKBGS samples have been published previously [19]. Briefly, genotyping was carried out using competitive allele-specific polymerase chain reaction KASPar chemistry (KBiosciences Ltd, Hoddesdon, Hertfordshire, UK). Call rates were 98.0% (BBCS) and 96.6% (UKBGS); there was no evidence for deviation from Hardy–Weinberg equilibrium (P = 0.29 (BBCS); P = 0.92 (UKBGS)), and the duplicate concordance based on a 1% (BBCS) and 5% (UKBGS) random sample of duplicates was 100% for both studies.

Statistical analysis

We estimated per-allele and genotypic log odds ratios (ORs) for the European, Asian and African-American subgroups separately using logistic regression, adjusted for principal components and study [5]. To test for departure from a multiplicative model we compared multiplicative and unconstrained models using a one degree of freedom likelihood ratio test. Heterogeneity in ORs between studies within each subgroup (European, Asian and African-American), and between subgroups, was assessed using the Cochrane Q statistic and quantified using the I2 measure [20].

Analyses stratified by oestrogen receptor status (+/–), progesterone receptor status (+/–), morphology (ductal or lobular), grade (1,2,3), lymph node involvement (+/–) or age at diagnosis (≤50 and >50 years) were restricted to studies of European ancestry due to the small number of studies of Asian and African-American ancestry. In addition, studies were excluded if they had selected cases on the basis of the stratifying variable, or had collected data on that variable for less than 5% of cases or less than 10 cases in total. Availability of data for each of the stratifying variables in each study is shown in Table S3 in Additional file 1. To assess the relationship between each of the stratifying variables and genotype, stratum-specific ORs were calculated using logistic regression. Cases in each stratum were compared with all control subjects, adjusted for study and principal components. Case-only logistic regression was used to test for heterogeneity between strata (binary stratifying variables) or across strata (stratifying variables with three or more strata). P values were estimated using likelihood ratio tests with one degree of freedom.

We assessed whether rs10235235 was associated with age at menarche in cases and controls separately. Studies that had not collected data on age at menarche in both cases and controls were excluded (Table S4 in Additional file 1). We used linear regression, adjusted for principal components and study, to estimate the relationship between age at menarche (years) and rs10235235 genotype (0, 1, 2 rare alleles) and logistic regression adjusted for principal components and study to estimate the association between age at menarche and breast cancer risk. To test for effect modification of an association between rs10235235 and breast cancer risk by age at menarche, we used logistic regression adjusted for principal components, study and age at menarche (grouped as ≤11, 12, 13, 14 and ≥15 years) with and without an interaction term(s). We considered four models: no interaction (zero interaction terms); assuming a linear interaction between genotype and menarche group (one interaction term); assuming a linear interaction between genotype and menarche group but allowing the linear term to differ between women who were heterozygous and those who were homozygous for the rare allele (two interaction terms); and one interaction term for each possible genotype/menarche group combination (eight interaction terms). Nested models were compared using likelihood ratio tests. All statistical analyses were performed using STATA version 11.0 (StataCorp, College Station, TX, USA). All P values reported are two-sided.

Results

The case–control analysis comprised genotype data for 47,346 invasive breast cancer cases and 47,569 controls from 49 studies, including 80,518 (84.8%) subjects of self-reported European ancestry, 12,419 (13.1%) of self-reported Asian ancestry and 1,978 (2.1%) of self-reported African-American ancestry. The mean (± standard deviation) age at diagnosis was 56.1 (± 11.6) years for European cases, 51.1 (± 10.5) years for Asian cases and 53.1 (± 10.7) years for African-American cases. There were ethnic differences in the estimated minor allele frequency (MAF) of rs10235235 (Q = 7317.1, two degrees of freedom; P for heterogeneity (Phet) = 0). The overall MAF for European control women was 0.089 (95% CI = 0.087, 0.091), but with strong evidence of between-study heterogeneity (Phet = 1 × 10-22) that was accounted for by the three Finnish studies (HEBCS, MAF = 0.15; KBCP, MAF = 0.21; and OBCS, MAF = 0.15; Phet = 0.01); no evidence of heterogeneity remained after taking account of these studies (MAF = 0.087 (95% CI = 0.085, 0.089); Phet = 0.23). Relative to Europeans, the overall MAF was higher for African-Americans (0.213, 95% CI = 0.195, 0.232; Phet = 0.26) but much lower for Asians (0.002; 95% CI = 0.001, 0.002), with strong evidence of between-study heterogeneity for the latter (Phet = 4 × 10-14).

The case–control analysis was consistent with a modest association between rs10235235 and breast cancer risk for women of European ancestry, with an estimated per-allele OR of 0.96 (95% CI = 0.93, 0.99; P for linear trend (Ptrend) = 0.02). Genotype-specific ORs were 0.98 (95% CI = 0.94, 1.01; P = 0.21) for AG versus AA (Figure 1A) and 0.80 (95% CI = 0.69, 0.93; P = 0.004) for GG versus AA (Figure 1B), with no evidence of between-study heterogeneity for either OR estimate (Phet = 0.44, I2 = 1.9% and Phet = 0.76, I2 = 0.0% for heterozygote and homozygote OR estimates respectively). There was, however, marginally significant evidence that the genotypic OR estimates departed from those expected under a multiplicative model with the inverse association of the GG genotype being more than the square of that of the AG genotype (test for deviation from multiplicative model, P = 0.04).

Figure 1
figure 1

Association of rs10235235 with breast cancer risk for women of European ancestry. Forest plots of the association of the rs10235235 AG (heterozygote) genotype (A) and GG (homozygote) genotype (B) with breast cancer risk for women of European ancestry. Horizontal lines, 95% confidence intervals (CIs); square boxes, study-specific fixed-effects estimates; diamond, combined, fixed-effects estimate of the odds ratio (OR) and 95% CI. Vertical line, null effect (OR = 1.0); dashed vertical line, estimated heterozygote OR (A) and estimated homozygote OR (B). Homozygote ORs for six studies (CTS, DEMOKRITOS, kConFab/AOCS, NBCS, NBHS and RPCI) could not be estimated because there were no GG homozygotes among cases or among controls in each of these studies (see Table S2 in Additional file 1).

Data for rs10235235 in women of Asian or African-American ancestry were more limited, with just two African-American studies (1,046 cases and 932 controls) and nine Asian studies (5,795 cases and 6,624 controls). In addition, this SNP was sufficiently rare in Asian populations (MAF = 0.002) that we were unable to estimate the heterozygote OR in two Asian studies (SEBCS, one carrier among 1,114 cases and no carriers among 1,129 controls; TWBCS, one carrier among 236 controls and no carriers among 774 cases; Table S2 in Additional file 1) and we could not estimate a homozygote OR for any Asian study (Table S2 in Additional file 1). There was no clear evidence that this SNP was associated with breast cancer risk for women of Asian ancestry (heterozygote OR = 1.06, 95% CI = 0.76, 1.49) or African-American ancestry (heterozygote and homozygote ORs were OR = 1.09, 95% CI = 0.90, 1.32 and OR = 0.94, 95% CI = 0.62, 1.42 respectively; Figure S1 in Additional file 1). This analysis, however, had low power to detect associations in non-Europeans and these OR estimates were not inconsistent with the magnitude of the observed OR estimates for European women (Phet = 0.51).

Stratifying cases by oestrogen receptor (Phet = 0.83) or progesterone receptor (Phet = 0.19) status, tumour grade (Phet = 0.63) or nodal involvement at diagnosis (Phet = 0.51) showed no evidence of effect modification (Table 1). There was some evidence of effect modification by morphology (Phet = 0.03). For ductal cancers we estimated a very modest reduction of risk for heterozygotes (ORhet = 0.98, 95% CI = 0.93, 1.02; P = 0.30) and a stronger, significant reduction for homozygotes (ORhom = 0.74, 95% CI = 0.61, 0.90; P = 0.003). For lobular cancers there was no such trend (ORhet = 1.07, 95% CI = 0.98, 1.17; P = 0.14 and ORhom = 0.91, 95% CI = 0.64, 1.27; P = 0.57).

Table 1 Association of rs10235235 with risk of breast cancer for women of European ancestry: stratified analysis

The SNP rs10235235 maps to a locus (CYP3A) that has been considered an a priori candidate for involvement in determining age at menopause and age at menarche [21, 22]. Stratifying cases by age at diagnosis (≤50 or >50 years) as a proxy for menopausal status at diagnosis showed no evidence of effect modification (Phet = 0.89; Table 2), and excluding cases who were diagnosed between age 46 and 55 as potentially perimenopausal did not alter this result (Phet = 0.28). Data on age at menarche were available for 21,736 cases and 22,686 controls (Table S4 in Additional file 1); to increase the power of the analysis we included additional data from BBCS and UKBGS (5,737 cases, 5,572 controls; Table S4 in Additional file 1) [19]. There was a 1.5% (95% CI = 0.5%, 2.7%; P = 0.004) reduction in breast cancer risk associated with each additional year’s increase in age at menarche. Mean age at menarche was positively associated with number of copies of the minor allele of rs10235235 for controls (Ptrend = 0.005; Table 3) but not for cases (Ptrend = 0.97; Table 3). Consequently, there was an inverse trend in the magnitude of the heterozygote and homozygote breast cancer ORs with mean age at menarche (Phet = 0.02; Table 4); being a carrier of one or two rare alleles of rs10235235 was associated with an estimated 16% (ORhet = 0.84, 95% CI = 0.75, 0.94; P = 0.003) or 19% (ORhom = 0.81, 95% CI = 0.51, 1.30; P = 0.39) (Ptrend = 0.002) reduction in breast cancer risk for women who had their menarche at ages ≥15 years but there was no evidence of reduction for those with a menarche at age ≤11 years (ORhet = 1.06, 95% CI = 0.95, 1.19; P = 0.30 and ORhom = 1.07, 95% CI = 0.67, 1.72; P = 0.78) (Ptrend = 0.29). There was no evidence that the inverse trend in the magnitude of ORs with mean age at menarche differed between heterozygous and homozygous carriers (P = 0.97) and no evidence that the trend was nonlinear (P = 0.70).

Table 2 rs10235235 and risk of breast cancer for women of European ancestry by age at diagnosis
Table 3 Association of rs10235235 with age at menarche for women of European ancestry by case-control status
Table 4 rs10235235 and risk of breast cancer for women of European ancestry by age at menarche

Discussion

This study of more than 47,000 breast cancer cases and 47,000 controls has confirmed that rs10235235, mapping to 7q22.1 (CYP3A), is associated with a reduction in breast cancer risk for women of European ancestry. Previously, our hypothesis-generating study of 10,000 breast cancer cases and 17,000 controls found a per-allele OR estimate of 0.96 (95% CI = 0.90, 1.02; P = 0.2), with marginally significant evidence of an inverse association for breast cancer diagnosed age 50 years or younger (OR = 0.91, 95% CI = 0.83, 0.99; P = 0.03) but no evidence of an association for breast cancer at later ages (OR = 1.01, 95% CI = 0.93, 1.10; P = 0.82) [19]. In this considerably larger study, we found a heterozygote OR estimate of 0.98 (95% CI = 0.94, 1.01; P = 0.21) and a homozygote OR estimate of 0.80 (95% CI = 0.69, 0.93; P = 0.004) with marginally significant evidence that the inverse association for homozygotes is greater than predicted by a multiplicative model (P = 0.04).

To our knowledge, rs10235235 is the first SNP to be associated with both breast cancer risk and age at menarche, consistent with the well-documented association between later age at menarche and a reduction in breast cancer risk [23]. Genome-wide association studies have identified more than 70 breast cancer risk variants [5, 6] and more than 30 variants associated with age at menarche [22], none of which map to the CYP3A locus. rs10235235 was originally identified on the basis of a highly significant association with hormone levels, accounting for 4.9% of the variation in premenopausal urinary oestrone glucuronide levels [19]. In this current analysis, rs10235235 accounted for only 0.01% of the variation across controls in age at menarche and we estimate that this SNP explains just 0.01% of the familial excess breast cancer risk. Our data thus illustrate the potential statistical efficiency of studies of intermediate phenotypes in the identification of rarer (MAF < 10%) risk alleles with modest associations. Our analysis shows some inconsistency with a recent genome-wide study of circulating oestradiol, testosterone and sex hormone-binding globulin in postmenopausal women [24]. In that study there was no genome-wide significant association observed with plasma oestradiol levels in either the primary analysis of approximately 1,600 postmenopausal women who were not taking postmenopausal hormones at blood draw or the secondary analysis that included approximately 900 current postmenopausal hormone users. Further studies will be needed to determine whether the lack of an association between CYP3A variants and postmenopausal plasma oestradiol levels reflects a difference in the menopausal status of the study subjects, the hormone/metabolite that was analysed or chance.

One possible explanation for the apparent effect modification of the rs10235235–breast cancer risk association by age at menarche is that this is a function of genotyping a marker SNP rather than the true causal variant. For example, if rs10235235 was perfectly correlated with a causal variant, SNP X, with a MAF substantially lower than that of rs10235235 (D′ ~ 1.0, r2 < 1.0), then there would be three types of chromosome in the population: type i, chromosomes carrying the common allele of rs10235235 and the common allele of SNP X; type ii, chromosomes carrying the rare allele of rs10235235 and the common allele of SNP X; and type iii, chromosomes carrying the rare allele of rs10235235 and the rare (protective) allele of SNP X. Only chromosomes carrying the rare allele of rs10235235 and the rare (protective) allele of SNP X (type iii) would be enriched in controls. Genotyping the marker (rs10235235) rather than the causal variant leads to misclassification. As the causal variant is associated with a protective effect on breast cancer risk, the proportion of chromosomes carrying both the rare allele of the causal variant and the marker (type iii) compared with the common allele of the causal variant and the rare allele of the marker (type ii) will be greater in controls than in cases such that the extent of misclassification will be greater for cases than controls. This will attenuate the association between genotype and age at menarche to a greater extent in cases than in controls creating an apparent effect modification. Fine mapping and functional studies will be required to identify the causal variant and to determine the true relationship between the causal variant, age at menarche and breast cancer risk.

Despite our original finding of a strong association between rs10235235 and hormone levels, we found no evidence that the association between this SNP and breast cancer risk differed by the hormone receptor status of the tumour, and nor did we find any evidence that the association differed by stage, grade or lymph node involvement. There was marginally significant evidence that the association between rs10235235 and breast cancer risk differed between ductal and lobular cancers (Phet = 0.03). Given the number of stratified analyses that we carried out (six stratifying variables) and given that there is no biological basis to support an interaction between rs10235235 and morphology, this is probably a chance observation.

In contrast to our earlier study [19], we found no evidence of an interaction with age at diagnosis when we stratified cases by age ≤/>50 years, either including or excluding cases diagnosed between age 46 and 55 years as potentially perimenopausal. We used age at diagnosis as a proxy for menopausal status at diagnosis because menopausal status at diagnosis is difficult to determine by questionnaire, especially given the use of hormone replacement therapies; while information on age at diagnosis was available for all but 1.4% (n = 554) of cases, information on age at natural menopause was missing for 65.6% (n = 26,552) of cases of European ancestry. Similarly, although rs10235235 is a plausible candidate for association with age at menopause, we did not test this due to the limited amount of data on age at natural menopause for controls of European ancestry (n = 11,294, 28.2%) and the difficulty in ascertaining whether treatment for breast cancer had influenced reported age at menopause for cases.

The strengths of our study include the large size of this combined analysis, and the availability of information on tumour characteristics for the majority of cases and on age at menarche for the majority of cases and controls. Limitations include low power of the study to examine an association between genotype and breast cancer risk for non-Europeans.

Conclusions

In summary, we have confirmed that rs10235235 is associated with breast cancer, have shown for the first time that rs10235235 is associated with age at menarche in controls and have suggested a potential mechanism for these associations. rs10235235, which maps to the CYP3A locus, probably tags a causal variant that affects expression of one or more CYP3A genes.

Abbreviations

BCAC:

Breast Cancer Association Consortium

CI:

confidence interval

COGS:

Collaborative Oncological Gene-environment Study

MAF:

minor allele frequency

OR:

odds ratio

P trend :

P value for linear trend

SNP:

single nucleotide polymorphism.

References

  1. Collaborative Group on Hormonal Factors in Breast Cancer: Familial breast cancer: collaborative reanalysis of individual data from 52 epidemiological studies including 58,209 women with breast cancer and 101,986 women without the disease. Lancet. 2001, 358: 1389-1399.

    Article  Google Scholar 

  2. Lichtenstein P, Holm NV, Verkasalo PK, Iliadou A, Kaprio J, Koskenvuo M, Pukkala E, Skytthe A, Hemminki K: Environmental and heritable factors in the causation of cancer–analyses of cohorts of twins from Sweden, Denmark, and Finland. N Engl J Med. 2000, 343: 78-85. 10.1056/NEJM200007133430201.

    Article  CAS  PubMed  Google Scholar 

  3. Peto J, Mack TM: High constant incidence in twins and other relatives of women with breast cancer. Nat Genet. 2000, 26: 411-414. 10.1038/82533.

    Article  CAS  PubMed  Google Scholar 

  4. Pharoah PD, Antoniou A, Bobrow M, Zimmern RL, Easton DF, Ponder BA: Polygenic susceptibility to breast cancer and implications for prevention. Nat Genet. 2002, 31: 33-36. 10.1038/ng853.

    Article  CAS  PubMed  Google Scholar 

  5. Michailidou K, Hall P, Gonzalez-Neira A, Ghoussaini M, Dennis J, Milne RL, Schmidt MK, Chang-Claude J, Bojesen SE, Bolla MK, Wang Q, Dicks E, Lee A, Turnbull C, Rahman N, Fletcher O, Peto J, Gibson L, Dos Santos Silva I, Nevanlinna H, Muranen TA, Aittomaki K, Blomqvist C, Czene K, Irwanto A, Liu J, Waisfisz Q, Meijers-Heijboer H, Adank M, van der Luijt RB, et al: Large-scale genotyping identifies 41 new loci associated with breast cancer risk. Nat Genet. 2013, 45: 353-361. 10.1038/ng.2563.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Garcia-Closas M, Couch FJ, Lindstrom S, Michailidou K, Schmidt MK, Brook MN, Orr N, Rhie SK, Riboli E, Feigelson HS, Le Marchand L, Buring JE, Eccles D, Miron P, Fasching PA, Brauch H, Chang-Claude J, Carpenter J, Godwin AK, Nevanlinna H, Giles GG, Cox A, Hopper JL, Bolla MK, Wang Q, Dennis J, Dicks E, Howat WJ, Schoof N, Bojesen SE, et al: Genome-wide association studies identify four ER negative-specific breast cancer risk loci. Nat Genet. 2013, 45: 392-398. 10.1038/ng.2561.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Dunning AM, Dowsett M, Healey CS, Tee L, Luben RN, Folkerd E, Novik KL, Kelemen L, Ogata S, Pharoah PD, Easton DF, Day NE, Ponder BA: Polymorphisms associated with circulating sex hormone levels in postmenopausal women. J Natl Cancer Inst. 2004, 96: 936-945. 10.1093/jnci/djh167.

    Article  CAS  PubMed  Google Scholar 

  8. Key T, Appleby P, Barnes I, Reeves G: Endogenous sex hormones and breast cancer in postmenopausal women: reanalysis of nine prospective studies. J Natl Cancer Inst. 2002, 94: 606-616. 10.1093/jnci/94.8.606.

    Article  CAS  PubMed  Google Scholar 

  9. Walker K, Bratton DJ, Frost C: Premenopausal endogenous oestrogen levels and breast cancer risk: a meta-analysis. Br J Cancer. 2011, 105: 1451-1457. 10.1038/bjc.2011.358.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Fortner RT, Eliassen AH, Spiegelman D, Willett WC, Barbieri RL, Hankinson SE: Premenopausal endogenous steroid hormones and breast cancer risk: results from the Nurses’ Health Study II. Breast Cancer Res. 2013, 15: R19-10.1186/bcr3394.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Helzlsouer KJ, Alberg AJ, Bush TL, Longcope C, Gordon GB, Comstock GW: A prospective study of endogenous hormones and breast cancer. Cancer Detect Prev. 1994, 18: 79-85.

    CAS  PubMed  Google Scholar 

  12. Wysowski DK, Comstock GW, Helsing KJ, Lau HL: Sex hormone levels in serum in relation to the development of breast cancer. Am J Epidemiol. 1987, 125: 791-799.

    CAS  PubMed  Google Scholar 

  13. Rosenberg CR, Pasternack BS, Shore RE, Koenig KL, Toniolo PG: Premenopausal estradiol levels and the risk of breast cancer: a new method of controlling for day of the menstrual cycle. Am J Epidemiol. 1994, 140: 518-525.

    CAS  PubMed  Google Scholar 

  14. Thomas HV, Key TJ, Allen DS, Moore JW, Dowsett M, Fentiman IS, Wang DY: A prospective study of endogenous serum hormone concentrations and breast cancer risk in premenopausal women on the island of Guernsey. Br J Cancer. 1997, 75: 1075-1079. 10.1038/bjc.1997.183.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Kabuto M, Akiba S, Stevens RG, Neriishi K, Land CE: A prospective study of estradiol and breast cancer in Japanese women. Cancer Epidemiol Biomarkers Prev. 2000, 9: 575-579.

    CAS  PubMed  Google Scholar 

  16. Kaaks R, Berrino F, Key T, Rinaldi S, Dossus L, Biessy C, Secreto G, Amiano P, Bingham S, Boeing H, de Mesquita HB B, Chang-Claude J, Clavel-Chapelon F, Fournier A, van Gils CH, Gonzalez CA, Gurrea AB, Critselis E, Khaw KT, Krogh V, Lahmann PH, Nagel G, Olsen A, Onland-Moret NC, Overvad K, Palli D, Panico S, Peeters P, Quiros JR, Roddam A, et al: Serum sex steroids in premenopausal women and breast cancer risk within the European Prospective Investigation into Cancer and Nutrition (EPIC). J Natl Cancer Inst. 2005, 97: 755-765. 10.1093/jnci/dji132.

    Article  CAS  PubMed  Google Scholar 

  17. Eliassen AH, Missmer SA, Tworoger SS, Spiegelman D, Barbieri RL, Dowsett M, Hankinson SE: Endogenous steroid hormone concentrations and risk of breast cancer among premenopausal women. J Natl Cancer Inst. 2006, 98: 1406-1415. 10.1093/jnci/djj376.

    Article  CAS  PubMed  Google Scholar 

  18. Stone J, Folkerd E, Doody D, Schroen C, Treloar SA, Giles GG, Pike MC, English DR, Southey MC, Hopper JL, Dowsett M: Familial correlations in postmenopausal serum concentrations of sex steroid hormones and other mitogens: a twins and sisters study. J Clin Endocrinol Metab. 2009, 94: 4793-4800. 10.1210/jc.2009-0713.

    Article  CAS  PubMed  Google Scholar 

  19. Johnson N, Walker K, Gibson LJ, Orr N, Folkerd E, Haynes B, Palles C, Coupland B, Schoemaker M, Jones M, Broderick P, Sawyer E, Kerin M, Tomlinson IP, Zvelebil M, Chilcott-Burns S, Tomczyk K, Simpson G, Williamson J, Hillier SG, Ross G, Houlston RS, Swerdlow A, Ashworth A, Dowsett M, Peto J, Dos Santos Silva I, Fletcher O: CYP3A variation, premenopausal estrone levels, and breast cancer risk. J Natl Cancer Inst. 2012, 104: 657-669. 10.1093/jnci/djs156.

    Article  CAS  PubMed  Google Scholar 

  20. Higgins JP, Thompson SG, Deeks JJ, Altman DG: Measuring inconsistency in meta-analyses. BMJ. 2003, 327: 557-560. 10.1136/bmj.327.7414.557.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Dvornyk V, Waqar ul H: Genetics of age at menarche: a systematic review. Hum Reprod Update. 2012, 18: 198-210. 10.1093/humupd/dmr050.

    Article  CAS  PubMed  Google Scholar 

  22. He C, Murabito JM: Genome-wide association studies of age at menarche and age at natural menopause. Mol Cell Endocrinol. 2014, 382: 767-779. 10.1016/j.mce.2012.05.003.

    Article  CAS  PubMed  Google Scholar 

  23. Colditz GA, Baer HJ, Tamimi RM: Breast cancer. Cancer Epidemiology and Prevention. Edited by: Schottenfeld D, Fraumeni JF. 2006, New York: Oxford University Press Inc, 995-1011. 3

    Chapter  Google Scholar 

  24. Prescott J, Thompson DJ, Kraft P, Chanock SJ, Audley T, Brown J, Leyland J, Folkerd E, Doody D, Hankinson SE, Hunter DJ, Jacobs KB, Dowsett M, Cox DG, Easton DF, De Vivo I: Genome-wide association study of circulating estradiol, testosterone, and sex hormone-binding globulin in postmenopausal women. PLoS One. 2012, 7: e37815-10.1371/journal.pone.0037815.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors thank all of the individuals who took part in these studies and all of the researchers, clinicians, technicians and administrative staff who have enabled this work to be carried out.

ABCFS would like to thank Maggie Angelakos, Judi Maskiell and Gillian Dite. ABCS would like to thank Ellen van der Schoot and Sanquin Amsterdam. The ACP study wishes to thank the participants in the Thai Breast Cancer study. Special thanks also go to the Thai Ministry of Public Health (MOPH) doctors and nurses who helped with the data collection process. The study would like to thank Dr Prat Boonyawongviroj, the former Permanent Secretary of MOPH and Dr Pornthep Siriwanarungsan, the Department Director-General of Disease Control who have supported the study throughout. BBCS would like to thank Eileen Williams, Elaine Ryder-Mills and Kara Sargus. BIGGS would like to thank Niall McInerney, Gabrielle Colleran, Andrew Rowan and Angela Jones. CNIO-BCS would like to thank Charo Alonso, Tais Moreno, Guillermo Pita, Primitiva Menendez and Anna González-Neira. The authors would like to acknowledge the contribution of the staff of the Génome Québec-genotyping unit under the supervision of Dr Sylvie LaBoissière, as well as Frédérick Robidoux from the McGill University and Génome Québec Innovation Centre. ESTHER would like to thank Hartwig Ziegler, Sonja Wolf and Volker Hermann. GC-HBOC would like to thank Bernd Frank. HEBCS would like to thank Dr Sofia Khan, Dr Kirsimari Aaltonen and Dr Karl von Smitten, and research nurses Irja Erkkilä and Virpi Palola. KBCP would like to thank Eija Myöhänen and Helena Kemiläinen. kConFab/AOCS would like to thank Heather Thorne, Eveline Niedermayr, the AOCS Management Group (D Bowtell, G Chenevix-Trench, A deFazio, D Gertig, A Green, P Webb) and the ACS Management Group (A Green, P Parsons, N Hayward, P Webb, D Whiteman). LAABC thanks all of the study participants and the entire data collection team, especially Annie Fung and June Yashiki. LMBC would like to thank Gilian Peuteman, Dominiek Smeets, Thomas Van Brussel and Kathleen Corthouts.

MARIE would like to thank Tracy Slanger, Elke Mutschelknauss, Ramona Salazar, S Behrens, R Birr, W Busch, U Eilber, B Kaspereit, N Knese and K Smit. MBCSG would like to thank Siranoush Manokian, Bernard Peissel and Daniela Zaffaroni of the Fondazione Istituto Nazionale dei Tumori, Milan, Monica Barile of the Istituto Europeo di Oncologia, Milan and Loris Bernard and personnel of the Cogentech Cancer Genetic Test Laboratory, Milan, Italy. MTLGEBCS would like to thank Martine Tranchant (Cancer Genomics Laboratory, CRCHUQ), Marie-France Valois, Annie Turgeon and Lea Heguy (McGill University Health Center, Royal Victoria Hospital; McGill University) for DNA extraction, sample management and skillful technical assistance. JS is Chairholder of the Canada Research Chair in Oncogenetics. MYBRCA would like to thank Phuah Sze Yee, Peter Kang, Kang In Nee, Kavitta Sivanandan, Shivaani Mariapun, Yoon Sook-Yee, Daphne Lee, Teh Yew Ching and Nur Aishah Mohd Taib for DNA Extraction and patient recruitment. NBHS thanks study participants and research staff for their contributions and commitment to the study. OBCS would like to thank Meeri Otsukka and Kari Mononen. OFBCR would like to thank Teresa Selander and Nayana Weerasooriya. ORIGO thanks E Krol-Warmerdam and J Blom for patient accrual, administering questionnaires and managing clinical information. The LUMC survival data were retrieved from the Leiden hospital-based cancer registry system (ONCDOC) with the help of Dr J Molenaar. PBCS would like to thank Louise Brinton, Mark Sherman, Stephen Chanock, Neonila Szeszenia-Dabrowska, Beata Peplonska, Witold Zatonski, Pei Chao and Michael Stagner. pKARMA would like to thank The Swedish Medical Research Counsel. RBCS would like to thank Petra Bos, Jannet Blom, Ellen Crepin, Elisabeth Huijskens, Annette Heemskerk and the Erasmus MC Family Cancer Clinic. SASBAC would like to thank The Swedish Medical Research Counsel. SBCGS thanks study participants and research staff for their contributions and commitment to the study. SBCS would like to thank Sue Higham, Helen Cramp and Dan Connley. SEARCH would like to thank The SEARCH and EPIC teams. SGBCC would like to thank the participants and research coordinator Kimberley Chua. SKKDKFZS are grateful to all of the patients for their participation and thank the physicians and other hospital staff, scientists, research assistants and study staff who contributed to the patient recruitment, data collection and sample preparation. UKBGS thanks Breakthrough Breast Cancer and the Institute of Cancer Research for support and funding of the Breakthrough Generations Study, and the study participants, study staff, and the doctors, nurses and other healthcare providers and health information sources who have contributed to the study.

Consortia members

The GENICA network: Dr Margarete Fischer-Bosch-Institute of Clinical Pharmacology, Stuttgart, and University of Tübingen, Germany (Christina Justenhoven, Hiltrud Brauch); Department of Internal Medicine, Evangelische Kliniken Bonn gGmbH, Johanniter Krankenhaus, Bonn, Germany (Yon-Dschun Ko, Christian Baisch); Institute of Pathology, University of Bonn, Germany (Hans-Peter Fischer); Molecular Genetics of Breast Cancer, Deutsches Krebsforschungszentrum (DKFZ), Heidelberg, Germany (Ute Hamann); Institute for Prevention and Occupational Medicine of the German Social Accident Insurance (IPA), Bochum, Germany (Thomas Bruening, Beate Pesch, Sylvia Rabstein, Anne Spickenheuer); and Institute for Occupational Medicine and Maritime Medicine, University Medical Center Hamburg-Eppendorf, Germany (Volker Harth).

kConFab Investigators: David Amor, Lesley Andrews, Yoland Antill, Shane Armitage, Rosemary Balleine, Agnes Bankier, Patti Bastick, John Beilby, Barbara Bennett, Ian Bennett, Anneke Blackburn, Michael Bogwitz, Meagan Brennan, Melissa Brown, Michael Buckley, Matthew Burgess, Jo Burke, Phyllis Butow, Ian Campbell, Alice Christian, Georgia Chenevix-Trench, Christine Clarke, Alison Colley, Dick Cotton, Bronwyn Culling, Margaret Cummings, Sarah-Jane Dawson, Anna DeFazio, Martin Delatycki, Rebecca Dickson, Alexander Dobrovic, Tracy Dudding, Ted Edkins, Stacey Edwards, Gelareh Farshid, Susan Fawcett, Georgina Fenton, Michael Field, James Flanagan, Peter Fong, John Forbes, Stephen Fox, Juliet French, Clara Gaff, Mac Gardner, Mike Gattas, Graham Giles, Grantley Gill, Jack Goldblatt, Sian Greening, Scott Grist, Eric Haan, Marion Harris, Stewart Hart, Nick Hayward, Sue Healey, Louise Heiniger, John Hopper, Clare Hunt, Paul James, Mark Jenkins, Rick Kefford, Alexa Kidd, Belinda Kiely, Judy Kirk, James Kollias, Jessica Koehler, Serguei Kovalenko, Sunil Lakhani, Jennifer Leary, Geoff Lindeman, Lara Lipton, Liz Lobb, Graham Mann, Deborah Marsh, Bettina Meiser, Roger Milne, Gillian Mitchell, Shona O’Connell, Nick Pachter, Briony Patterson, Lester Peters, Kelly Phillips, Melanie Price, Lynne Purser, Tony Reeve, Edwina Rickard, Bridget Robinson, Barney Rudzki, Elizabeth Salisbury, Christobel Saunders, Joe Sambrook, Jodi Saunus, Robyn Sayer, Clare Scott, Elizabeth Scott, Rodney Scott, Adrienne Sexton, Raghwa Sharma, Andrew Shelling, Peter Simpson, Melissa Southey, Amanda Spurdle, Graeme Suthers, Pamela Sykes, Jessica Taylor, Ella Thompson, Heather Thorne, Sharron Townshend, Alison Trainer, Kathy Tucker, Janet Tyler, Jane Visvader, Logan Walker, Paul Waring, Robin Ward, Bev Warner, Rachael Williams, Ingrid Winship, Mary Ann Young (Peter MacCallum Cancer Center, Melbourne, Australia).

The Australian Ovarian Cancer Study Group: David D Bowtell, Adele C Green, Georgia Chenevix-Trench, Anna deFazio, Dorota Gertig, Penelope M Webb (Peter MacCallum Cancer Center, Melbourne, Australia).

Financial support

Part of this work was supported by the European Community’s Seventh Framework Programme under grant agreement number 223175 (grant number HEALTH-F2-2009-223175) (COGS). This work was partly supported by the Canadian Institutes of Health Research for the ‘CIHR Team in Familial Risks of Breast Cancer’ program (JS, DFE), and the Ministry of Economic Development, Innovation and Export Trade of Quebec – grant number PSR-SIIRI-701 (JS, DFE, PH).

The ABCFS and OFBCR work was supported by the United States National Cancer Institute, National Institutes of Health (NIH) under RFA-CA-06-503 and through cooperative agreements with members of the Breast Cancer Family Registry (BCFR) and Principal Investigators, including Cancer Care Ontario (U01 CA69467), Northern California Cancer Center (U01 CA69417) and University of Melbourne (U01 CA69638). Samples from the NC-BCFR were processed and distributed by the Coriell Institute for Medical Research. The content of this manuscript does not necessarily reflect the views or policies of the National Cancer Institute or any of the collaborating centers in the BCFR, nor does mention of trade names, commercial products, or organizations imply endorsement by the US Government or the BCFR. ABCFS was also supported by the National Health and Medical Research Council of Australia, the New South Wales Cancer Council, the Victorian Health Promotion Foundation (Australia) and the Victorian Breast Cancer Research Consortium. JLH is a National Health and Medical Research Council (NHMRC) Australia Fellow and a Victorian Breast Cancer Research Consortium Group Leader. MCS is a NHMRC Senior Research Fellow and a Victorian Breast Cancer Research Consortium Group Leader. The ABCS study was supported by the Dutch Cancer Society (grants NKI 2001-2423 and 2007-3839) and the Dutch National Genomics Initiative. The ACP study is funded by the Breast Cancer Research Trust, UK. The work of the BBCC was partly funded by ELAN-Fond of the University Hospital of Erlangen. BBCS is funded by Cancer Research UK and Breakthrough Breast Cancer, and acknowledges NHS funding to the NIHR Biomedical Research Centre and the National Cancer Research Network. BCAC is funded by CR-UK (C1287/A10118 and C1287/A12014). Meetings of the BCAC have been funded by the European Union COST programme (BM0606). DFE is a Principal Research Fellow of CR-UK. ES (BIGGS) is supported by NIHR Comprehensive Biomedical Research Centre, Guy’s & St. Thomas’ NHS Foundation Trust in partnership with King’s College London, UK. IT is supported by the Oxford Biomedical Research Centre. The BSUCH study was supported by the Dietmar-Hopp Foundation, the Helmholtz Society and the German Cancer Research Center (DKFZ). CGPS was supported by the Chief Physician Johan Boserup and Lise Boserup Fund, the Danish Medical Research Council and Herlev Hospital. CNIO-BCS was supported by the Genome Spain Foundation, the Red Temática de Investigación Cooperativa en Cáncer and grants from the Asociación Española Contra el Cáncer and the Fondo de Investigación Sanitario (PI081583 and PI081120). CTS was supported by the California Breast Cancer Act of 1993, the NIH (grants R01 CA77398 and the Lon V Smith Foundation (LVS39420)) and the California Breast Cancer Research Fund (contract 97-10500). Collection of cancer incidence data used in this study was supported by the California Department of Public Health as part of the statewide cancer reporting program mandated by California Health and Safety Code Section 103885. The ESTHER study was supported by a grant from the Baden Württemberg Ministry of Science, Research and Arts. Additional cases were recruited in the context of the VERDI study, which was supported by a grant from the German Cancer Aid (Deutsche Krebshilfe).

GC-HBOC was supported by Deutsche Krebshilfe (107054), the Dietmar-Hopp Foundation, the Helmholtz Society and the German Cancer Research Centre (DKFZ). GENICA was funded by the Federal Ministry of Education and Research (BMBF) Germany grants 01KW9975/5, 01KW9976/8, 01KW9977/0 and 01KW0114, the Robert Bosch Foundation, Stuttgart, Deutsches Krebsforschungszentrum (DKFZ), Heidelberg, Institute for Prevention and Occupational Medicine of the German Social Accident Insurance (IPA), Bochum, as well as the Department of Internal Medicine, Evangelische Kliniken Bonn gGmbH, Johanniter Krankenhaus, Bonn, Germany. HEBCS was supported by the Academy of Finland (132473), Helsinki University Central Hospital Research Fund, the Sigrid Juselius Foundation, the Finnish Cancer Society and the Nordic Cancer Union. HERPACC was supported by a Grant-in-Aid for Scientific Research on Priority Areas and on Innovative Area from the Ministry of Education, Science, Sports, Culture and Technology of Japan and by a Grant-in-Aid for the Third Term Comprehensive 10-Year Strategy for Cancer Control from Ministry Health, Labour and Welfare of Japan. HMBCS was supported by short-term fellowships from the German Academic Exchange Program (NVB) and the Friends of Hannover Medical School (NVB). KBCP was financially supported by the special Government Funding (EVO) of Kuopio University Hospital grants, Cancer Fund of North Savo, the Finnish Cancer Organizations, the Academy of Finland and by the strategic funding of the University of Eastern Finland. kConFab is supported by grants from the National Breast Cancer Foundation, the NHMRC, the Queensland Cancer Fund, the Cancer Councils of New South Wales, Victoria, Tasmania and South Australia and the Cancer Foundation of Western Australia. The kConFab Clinical Follow Up Study was funded by the NHMRC (145684, 288704, 454508). Financial support for the AOCS was provided by the United States Army Medical Research and Materiel Command (DAMD17-01-1-0729), the Cancer Council of Tasmania and Cancer Foundation of Western Australia and the NHMRC (199600). GC-T and P Webb are supported by the NHMRC. LAABC is supported by grants (1RB-0287, 3PB-0102, 5PB-0018, 10PB-0098) from the California Breast Cancer Research Program. Incident breast cancer cases were collected by the USC Cancer Surveillance Program (CSP), which is supported under subcontract by the California Department of Health. CSP is also part of the National Cancer Institute’s Division of Cancer Prevention and Control Surveillance, Epidemiology, and End Results Program, under contract number N01CN25403. LMBC is supported by the ‘Stichting tegen Kanker’ (232-2008 and 196-2010). DL is supported by the KULPFV/10/016-SymBioSysII.

The MARIE study was supported by the Deutsche Krebshilfe e.V. (70-2892-BR I), the Hamburg Cancer Society, the German Cancer Research Center and the genotype work in part by the Federal Ministry of Education and Research (BMBF) Germany (01KH0402). MBCSG was funded by grants from Italian Association for Cancer Research (AIRC, IG 8713), and by Italian citizens who allocated the 5 × 1,000 share of their tax payment in support of the Fondazione IRCCS Istituto Nazionale dei Tumori, according to Italian laws (INT-Institutional strategic projects ‘5x1000’). MCBCS was supported by the NIH grants CA116167 and CA128978, an NIH Specialized Program of Research Excellence (SPORE) in Breast Cancer (CA116201), the Breast Cancer Research Foundation, and a generous gift from the David F and Margaret T Grohne Family Foundation and the Ting Tsung and Wei Fong Chao Foundation. MCCS cohort recruitment was funded by VicHealth and Cancer Council Victoria. MCCS was further supported by Australian NHMRC grants 209057, 251553 and 504711 and by infrastructure provided by Cancer Council Victoria. MEC was support by NIH grants CA63464, CA54281, CA098758 and CA132839. The work of MTLGEBCS was supported by the Quebec Breast Cancer Foundation, the Canadian Institutes of Health Research for the ‘CIHR Team in Familial Risks of Breast Cancer’ program (grant number CRN-87521) and the Ministry of Economic Development, Innovation and Export Trade (grant number PSR-SIIRI-701). MYBRCA is funded by research grants from the Malaysian Ministry of Science, Technology and Innovation (MOSTI), Malaysian Ministry of Higher Education (UM.C/HlR/MOHE/06) and Cancer Research Initiatives Foundation (CARIF). Additional controls were recruited by the Singapore Eye Research Institute, which was supported by a grant from the Biomedical Research Council, Singapore (BMRC08/1/35/19/550) and the National Medical Research Council, Singapore (NMRC/CG/SERI/2010). NBCS was supported by grants from the Norwegian Research council, 155218/V40, 175240/S10 to A-LB-D, FUGE-NFR 181600/V11 to VNK and a Swizz Bridge Award to A-LB-. NBHS was supported by NIH grant R01CA100374. Biological sample preparation was conducted the Survey and Biospecimen Shared Resource, which is supported by P30 CA68485. OBCS was supported by research grants from the Finnish Cancer Foundation, the Academy of Finland, the University of Oulu, and the Oulu University Hospital. The ORIGO study was supported by the Dutch Cancer Society (RUL 1997-1505) and the Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL CP16). PBCS was funded by Intramural Research Funds of the National Cancer Institute, Department of Health and Human Services, USA. The pKARMA study was supported by Märit and Hans Rausings Initiative Against Breast Cancer. RBCS was funded by the Dutch Cancer Society (DDHK 2004-3124, DDHK 2009-4318).

The SASBAC study was supported by funding from the Agency for Science, Technology and Research of Singapore (A*STAR), the US NIH and the Susan G Komen Breast Cancer Foundation. SBCGS was supported primarily by NIH grants R01CA64277, R01CA148667, and R37CA70867. Biological sample preparation was conducted the Survey and Biospecimen Shared Resource, which is supported by P30 CA68485. SBCS was supported by Yorkshire Cancer Research S295, S299, S305PA. SCCS is supported by a grant from the National Institutes of Health (R01 CA092447). Data on SCCS cancer cases used in this publication were provided by the Alabama Statewide Cancer Registry; Kentucky Cancer Registry, Lexington, KY; Tennessee Department of Health, Office of Cancer Surveillance; Florida Cancer Data System; North Carolina Central Cancer Registry, North Carolina Division of Public Health; Georgia Comprehensive Cancer Registry; Louisiana Tumor Registry; Mississippi Cancer Registry; South Carolina Central Cancer Registry; Virginia Department of Health, Virginia Cancer Registry; and Arkansas Department of Health, Cancer Registry, Little Rock. The Arkansas Central Cancer Registry is fully funded by a grant from National Program of Cancer Registries, Centers for Disease Control and Prevention (CDC). Data on SCCS cancer cases from Mississippi were collected by the Mississippi Cancer Registry, which participates in the National Program of Cancer Registries of the CDC. The contents of this publication are solely the responsibility of the authors and do not necessarily represent the official views of the CDC or the Mississippi Cancer Registry. SEARCH is funded by programme grants from Cancer Research UK (C490/A10124 and C8197/A10123) and NIH grant 5U01CA098216-07. SEBCS was supported by the Korea Health 21 R&D Project (AO30001), Ministry of Health and Welfare, Republic of Korea. SGBCC is funded by the National Medical Research Council start-up Grant and Centre Grant (NMRC/CG/NCIS /2010). Additional controls were recruited by the Singapore Consortium of Cohort Studies-Multi-ethnic cohort (SCCS-MEC), which was funded by the Biomedical Research Council, grant number: 05/1/21/19/425. SKKDKFZS is supported by the DKFZ, Heidelberg, Germany. KJ-B (SZBCS) is a fellow of International PhD program, Postgraduate School of Molecular Medicine, Warsaw Medical University, supported by the Polish Foundation of Science. TBCS was funded by The National Cancer Institute, Thailand. TNBCC was supported by an NIH Specialized Program of Research Excellence (SPORE) in Breast Cancer (CA116201), the Breast Cancer Research Foundation, a generous gift from the David F and Margaret T Grohne Family Foundation and the Ting Tsung and Wei Fong Chao Foundation; The Stefanie Spielman Breast Cancer Fund and the OSU Comprehensive Cancer Center; the European Union (European Social Fund – ESF) and Greek national funds through the Operational Program ‘Education and Lifelong Learning’ of the National Strategic Reference Framework (NSRF) – Research Funding Program of the General Secretariat for Research & Technology: ARISTEIA. TWBCS is supported by the Taiwan Biobank project of the Institute of Biomedical Sciences, Academia Sinica, Taiwan. UKBGS is funded by Breakthrough Breast Cancer and the Institute of Cancer Research (ICR). ICR acknowledges NHS funding to the NIHR Biomedical Research Centre.

Author information

Authors and Affiliations

Authors

Consortia

Corresponding author

Correspondence to Nichola Johnson.

Additional information

Competing interests

The authors state that they have no competing interests.

Authors’ contributions

OF, FD and NO performed the statistical analyses. OF, IdSS and NJ drafted the manuscript. NJ, FD, NO, LG, MEJ, MJS, EJF, BPH, MG-C, MDo, AA, AJS, JP, IdSS and OF comprised the writing group that was responsible for the interpretation of the results and for critically reviewing the manuscript. AC, AJ, AHW, AMa, BBu, C-YS, DL, ES, GC-T, HN, HBre, HBra, ILA, JC-C, J-YC, JLH, LBa, MKB, HMi, PAF, PR, RW, SEB, TD, MKS and UH also significantly contributed to the interpretation of the results. OF, IdSS, NJ, JP, LG, DFE, MKB and JW conceived of the original design of the study and participated in subject recruitment and in acquisition of data. JBen, AG-N, RM, DCT, DV, FB, CL, JD, JS and KMi carried out the genotyping and/or data analysis. FD, NO, MEJ, MJS, EJF, BPH, JLH, MCS, GSD, CA, MKS, AB, LJVV, FA, KMu, ALo, PAF, MWB, ABE, SPR, ES, IT, MK, NM, BBu, FMa, AS, CS, PG, TT, EC, FMe, SEB, BGN, HF, RMi, MPZ, JIAP, JBen, LBe, HA-C, AZ, CCD, HBre, HMü, VA, AKD, AMe, JH, CRB, RKS, HBra, CJ, Y-DK, The GENICA Network, HN, TAM, KA, CB, KMa, TD, NVB, NNA, ALi, AMa, VK, V-MK, JMH, GC-T, JBee, kConFab Investigators, Australian Ovarian Cancer Study Group, AHW, DVdB, C-CT, DL, DS, PN, HW, JC-C, AR, SN, DF-J, PR, PP, BBo, VP, FJC, JEO, XW, ZF, VSP, GGG, GS, LBa, CH, JS, MSG, FL, MDu, PS, ST, CHY, SYP, BKC, VNK, GGA, A-LB-D, WZ, RW, KP, AJ-V, MG, ILA, JAK, GG, AMM, PD, JF, SJC, JLis, MES, PH, NS, MHo, AH, RAO, MT-L, JLiu, AC, IWB, MWRR, SSC, WB, LBS, PDPP, AMD, MS, DK, D-YN, SKP, J-YC, MHa, HMi, WYL, AT, UH, AF, TR, HUU, AJ, JLu, KJ-B, KD, SSa, VG, PB, JM, SSl, AET, CV, DY, C-YS, J-CY, C-SH, M-FH, AG-N, DCT, DV, FB, CL, JD, KMi, MKB, JW, DFE, MG-C, MDo, AA and AJS made substantial contributions in recruiting subjects and acquiring data, and in critically reviewing the manuscript. All authors take responsibility for the work and read and approved the final version of the manuscript.

Isabel dos Santos Silva and Olivia Fletcher contributed equally to this work.

Electronic supplementary material

13058_2013_3558_MOESM1_ESM.zip

Additional file 1: Contains Table S1 presenting details of participating BCAC studies; Table S2 presenting rs10235235 genotypes for breast cancer cases and controls from 49 BCAC studies; Table S3 presenting availability of data on age at diagnosis, hormone receptor status, morphology, grade and nodal status for breast cancer cases from 38 European BCAC studies; Table S4 presenting availability of data on age at menarche for breast cancer cases and controls from 40 European BCAC studies; and Figure S1 showing association of the rs10235235-AG genotype with breast cancer risk for women of Asian and African-American ancestry.(ZIP 213 KB)

Additional file 2: Presents details of ethical committees that approved each study.(XLSX 12 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

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Johnson, N., Dudbridge, F., Orr, N. et al. Genetic variation at CYP3A is associated with age at menarche and breast cancer risk: a case-control study. Breast Cancer Res 16, R51 (2014). https://doi.org/10.1186/bcr3662

Download citation

  • Received:

  • Accepted:

  • Published:

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

Keywords