- Open Access
Genetic variants in the MRPS30 region and postmenopausal breast cancer risk
Genome Medicinevolume 3, Article number: 42 (2011)
The Erratum to this article has been published in Genome Medicine 2012 4:19
Genome-wide association studies have identified several genomic regions that are associated with breast cancer risk, but these provide an explanation for only a small fraction of familial breast cancer aggregation. Genotype by environment interactions may contribute further to such explanation, and may help to refine the genomic regions of interest.
We examined genotypes for 4,988 SNPs, selected from recent genome-wide studies, and four randomized hormonal and dietary interventions among 2,166 women who developed invasive breast cancer during the intervention phase of the Women's Health Initiative (WHI) clinical trial (1993 to 2005), and one-to-one matched controls. These SNPs derive from 3,224 genomic regions having pairwise squared correlation (r2) between adjacent regions less than 0.2. Breast cancer and SNP associations were identified using a test statistic that combined evidence of overall association with evidence for SNPs by intervention interaction.
The combined 'main effect' and interaction test led to a focus on two genomic regions, the fibroblast growth factor receptor two (FGFR2) and the mitochondrial ribosomal protein S30 (MRPS30) regions. The ranking of SNPs by significance level, based on this combined test, was rather different from that based on the main effect alone, and drew attention to the vicinities of rs3750817 in FGFR2 and rs7705343 in MRPS30. Specifically, rs7705343 was included with several FGFR2 SNPs in a group of SNPs having an estimated false discovery rate < 0.05. In further analyses, there were suggestions (nominal P < 0.05) that hormonal and dietary intervention hazard ratios varied with the number of minor alleles of rs7705343.
Genotype by environment interaction information may help to define genomic regions relevant to disease risk. Combined main effect and intervention interaction analyses raise novel hypotheses concerning the MRPS30 genomic region and the effects of hormonal and dietary exposures on postmenopausal breast cancer risk.
Genome-wide association studies have identified a substantial number of common genetic variants that are associated with risk, for each of several diseases. However, most such associations are weak and account for only a small fraction of familial disease aggregation . In the case of breast cancer, seven reproducible genetic susceptibility alleles were estimated to explain about 5% of heritability . Studies of low frequency genetic variants, gene-gene interactions, genotype by environment interaction, and shared environment have been suggested  as means to identify the 'missing heritability' for complex diseases, along with more thorough study of variants within genomic regions of interest.
Closely related to this is the role of genetic variants in model discrimination and disease risk prediction. A recent multiple-cohort analysis of ten common genetic variants that reliably associate with breast cancer concluded that 'the level of predicted breast cancer risk among most women changed little' when these SNPs were added to existing risk assessment models . In response, an accompanying editorial  pointed out that cellular networks within which the SNPs operate may associate more strongly with risk than do tagging SNPs alone, that gene-gene and gene-environment interactions are 'likely to be profoundly important', and that associations with breast cancer subtypes may be more impressive.
A challenge to pursuing the gene-environment concept is the typical difficulty in assessing key environmental exposures. For example, given the well-established association between obesity and post-menopausal breast cancer risk, one might expect that total energy consumption and other dietary factors may influence breast cancer risk, possibly in a manner that depends on genetic factors that relate to hormone metabolism, growth factors, or inflammation. However, dietary data are attended by random and systematic assessment biases that may seriously attenuate and distort estimated associations .
Randomized controlled intervention trials can provide highly desirable settings for the incorporation of genotype by environment interactions into genetic association analyses. First, the intervention group assignment is known with precision, and secondly, this assignment is statistically independent of underlying genotype by virtue of randomization. This latter feature also allows highly efficient case-only test statistics [6–8] to be used for genotype by intervention interaction testing.
The Women's Health Initiative (WHI) randomized controlled trial included four randomized and controlled comparisons among postmenopausal women in a partial factorial design [9, 10]. Specifically, it comprised a postmenopausal hormone therapy component that involved two non-overlapping trials: estrogen versus placebo (E-alone trial) among women who were post-hysterectomy, and estrogen plus progestin versus placebo (E+P trial) among women with a uterus; a low-fat dietary modification (DM) versus usual diet component, and a calcium and vitamin D (CaD) versus placebo supplementation component.
An elevation of breast cancer risk triggered the early stopping of the E+P trial in 2002 [11, 12]. In the E-alone trial, which was stopped early in 2004 primarily due to an elevation of stroke risk , there was a surprising suggestion of a reduction in breast cancer risk in the intervention group, as well as apparent interactions of the E-alone hazard ratio with several other breast cancer risk factors . The DM trial continued to its planned termination in 2005. While overall it provided non-significant evidence of a breast cancer reduction over its 8.1-year average follow-up period, the breast cancer hazard ratio was significantly lower in the quartile of women who had a comparatively high fat content in their diet at baseline . These women made a larger dietary change if assigned to the low-fat diet intervention. The CaD trial did not yield evidence of an effect on breast cancer risk .
We studied 4,988 SNPs in relation to breast cancer incidence and clinical trials intervention effects during the intervention phase of the WHI clinical trial. Nearly all of these SNPs were selected as the top-ranked SNPs according to significance level for association with breast cancer in the NCI Cancer Genetic Markers of Susceptibility (C-GEMS) genome-wide association study , while the remaining 244 were selected based on published data from the Breast Cancer Association Consortium genome-wide association study . These SNPs were scattered throughout the genome. In fact, they arise from 3,224 distinct loci when a squared pairwise correlation (r2) between adjacent regions of less than 0.2 is used to define new loci. We ranked SNPs according to a null hypothesis test that combined evidence of overall breast cancer association with evidence of interaction with one or more of the randomized clinical trial intervention assignments.
Materials and methods
Study design and population
Enrollees in WHI trials were postmenopausal women aged 50 to 79 years who met component-specific eligibility criteria . Women were randomized to a hormone therapy component, or a DM component, or both. At the one-year anniversary from enrollment, participating women could be further randomized into a CaD supplementation component. A total of 68,132 women were enrolled into the trials between 1993 and 1998, among which there were 10,739 in E-alone, 16,608 in E+P, 48,835 in DM, and 36,282 in CaD components. Details about distributions of demographic variables and breast cancer risk factors in the study cohort were published previously . For the DM trial we chose to focus interaction testing on the subset of 12,208 women having baseline percentage of energy from fat in the upper quartile, and we denote the DM intervention in this sub-cohort by DMQ.
Case and control selection
All 2,242 invasive breast cancer cases that developed between randomization and the end of the trial intervention phase (31 March 2005) were considered for inclusion, among which a total of 2,166 (96.6%) cases had adequate quantity and quality of DNA. This leads to analyses based on 247 cases for E-alone, 471 cases for E+P, 428 cases for DMQ, 1,049 cases for CaD (cases arising after CaD randomization only), and corresponding controls that were one-to-one matched to cases on baseline age, self-reported ethnicity, participation in each trial component, years since randomization, and baseline hysterectomy status.
Genotyping and data cleaning methods at Perlegen Sciences (Mountain View, CA, USA) have been described . The average call rate for these SNPs was 99.8%, and the average concordance rate for 157 blind duplicate samples was also 99.8%.
Principal component analysis was used to characterize population structure and to identify genotyping artifacts. The top 20 principal components did not associate with common sources of experimental variability (for example, date of sample processing or hybridization performance for either chip design). The first ten principal components were found to account for 86% of the total SNP genotype variation, while the first four principal components provided good separation among the major self-reported 'ethnicities' (white, black, Hispanic, Asian/Pacific Islander, northern versus southern European ancestry).
A five-component test statistic was used for each SNP to test association with breast cancer. The first 'main effect' component arose as score test from a standard logistic regression of case (1) versus control (0) status on number of minor SNP alleles and potential confounding factors. The logistic regression model included the (log transformed) Gail 5-year breast cancer risk score , previous hormone use (indicators for < 5, 5 to 10, and ≥10 years for each of estrogen and estrogen plus progestin), and (log transformed) body mass index. Also included are variables used for matching controls to cases in control selection. In addition, eigenvectors from the first ten principal components from correlation analysis of the genotype data were included to adjust for population stratification . The other four test statistic components were case-only tests for dependence of intervention odds ratios on SNP genotype for each of E-alone, E+P, DMQ, and CaD. These statistics arise as score tests in logistic regression of active (1) versus placebo or usual diet (0) randomization assignment on the number of minor SNP alleles with logistic regression location parameter offset by log q/(1 - q), where q is the fraction of women assigned to active intervention for the pertinent clinical trial component. The main effect test statistic is asymptotically independent of each of the case-only test statistics , and the interaction tests for E-alone and E+P are independent since they are based on non-overlapping sets of women. A 'sandwich' variance estimator was used to allow for possible correlations among the other pairs of case-only test statistics. A chi-square test with five degrees of freedom was then used to test SNP association with breast cancer, for each of the SNPs. Further details about this joint test procedure are included here as Additional file 1.
SNPs of interest in these association tests were subsequently examined for evidence of main effect and interaction effects separately. The latter once again employed case-only analyses, and for descriptive purposes, intervention odds ratios were estimated separately at zero, one, and two minor SNP alleles. A likelihood ratio test with two degrees of freedom assessed SNP by intervention interaction in these analyses.
The potential of SNP by clinical trial interactions to contribute to the ability to discriminate between breast cancer cases and controls was evaluated by estimating areas under the receiver operating characteristic curves (AUC), and associated confidence intervals.
Some further analyses were carried out with breast cancers classified according to either the estrogen receptor status or the progesterone receptor status of the breast tumor. All significance levels (P-values) are two-sided.
This research conforms to the Helsinki Declaration and pertinent legislation, and has been approved by the Institutional Review Board of the Fred Hutchinson Cancer Research Center. All women included in this report provided informed consent that permitted their biospecimens and data to be used in the present research project.
Simultaneous tests of main effect and interaction with clinical trial interventions
Table 1 presents the top 20 SNPs ranked by P-value of the combined test of main effect and interaction. Among the 4,988 SNPs evaluated, six SNPs have the joint test P-value less than 10-6 and a false discovery rate (FDR) less than 0.0005, all in the FGFR2 (fibroblast growth factor receptor 2) region in chromosome region 10q16. Immediately following are several SNPs from the MRPS30 (mitochondrial ribosomal protein S30) region in chromosome region 5p12. Of these SNPs, rs7705343 is included in the set of SNPs having FDR < 0.05, while close-by SNP rs13159598 is also among SNPs having FDR < 0.10.
Table 1 also shows P-values and rankings for these SNPs under the main effect association test alone. While P-values for FGFR2 SNPs tend to be somewhat diluted by the inclusion of the interaction information in the test statistic, the ordering of these SNPs is rather different under the two-testing procedures. For example, SNP rs3750817, which is in a somewhat separate linkage disequilibrium bin from tagging SNP rs2981582 , has a comparatively higher ranking with the combined test. We have previously reported suggestive evidence of interaction of rs3750817 with E-alone and E+P , and DMQ .
SNPs in the MRPS30 region of chromosome 5p12 have a higher ranking overall with the combined versus the main effect test. Moreover, the ordering of SNPs within this region is considerably altered by the inclusion of the interaction information. These analyses point to the genomic region in proximity of rs7705343 as relevant to breast cancer risk. Figure 1 shows squared pairwise correlations (r2) among SNPs in the MRPS30 region of chromosome 5p12. The combined test rankings tend to decrease as one moves from rs7705343 to the tagging SNP rs4415084 at the opposite end of this genomic region of approximately 230 kb.
Table 2 shows P-values individually for the five components of the combined test, for the eight SNPs in the MRPS30 region. Most of the association information derives from the main effect test, but the intervention interaction tests have rather different P-values across these SNPs, with rs7705343 having nominally significant (P < 0.05) interactions with each of E-alone, DMQ, and CaD, while interactions in relation to rs4415084 are not significant for any of the interventions.
Table 3 shows estimated intervention odds ratios and 95% confidence intervals as a function of the number of minor alleles of rs7705343 for each of the four interventions. The GG genotype is associated with lower intervention ORs for each of E-alone, DMQ, and CaD. Additional file 2 provides corresponding information with breast cancers classified according to estrogen receptor or progesterone receptor positivity. No clear variations by tumor receptor status were suggested, through statistical power for detecting moderate variations with tumor type is limited.
The majority (86%) of the case-control samples are from European-ancestry populations. In Additional files 3 and 4 we provide P-values for interaction between trial components and SNPs in the MRPS30 region, and the estimated intervention odds ratios and 95% confidence intervals as a function of the number of minor alleles of rs7705343 among women of European ancestry specifically. The patterns that we observe are quite similar to the overall patterns.
We also examined the joint associations of these FGFR2 and MRPS30 SNPs with hormonal and dietary intervention effects, using case-only analysis. Based on logistic regression applied to cases in DMQ, where the indicator for active treatment is regressed on genotypes of rs3750817 and rs7705343 together, both SNPs showed nominally significant interactions. The P-values for rs3750817 and rs7705343 were 0.0059 and 0.037. When E-alone was similarly considered, rs3750817 and rs7705343 had P-values of 0.053 and 0.043 in the joint interaction model.
The AUC was calculated from logistic regression analyses that included clinical trial randomization assignments for each of the four interventions and potential confounding factors. This gave an AUC (95% confidence interval) of 0.594 (0.578, 0.611). When main effect indicator variables were added for one and two minor alleles of rs3750817 and rs7705343, the AUC increased to 0.610 (0.594, 0.627). When SNP by intervention interaction indicator variables were also included, the AUC increased further to 0.621 (0.604, 0.637). A bootstrap test of significance for the genotype by intervention terms gave a nominal P-value of 0.007.
We evaluated the association between 4,988 SNPs and invasive breast cancer incidence in the WHI clinical trial through the use of a statistic that combines SNP main effect information with SNP by intervention interaction information for each of four randomized interventions. This view of the data provided a clear focus on two genomic regions, the FGFR2 region of chromosome 10 q, which has a very strong main effect along with suggestive evidence for interaction, and the MRPS30 region of chromosome 5 p, which shows evidence of a comparatively smaller main effect and suggestive evidence for interaction. The inclusion of the clinical trial interventions in this testing procedure leads to interest in subregions containing FGFR2 SNP rs3750817 and MRPS30 SNP rs7705343 that are some distance from their associated tagging SNPs, possibly suggesting more than one regulatory element in these non-coding genomic regions.
We have previously [9, 10] discussed these data in relation to FGFR2. The eight MRPS30 SNPs considered here fall in a linkage disequilibrium region of approximately 230 kb from downstream of fibroblast growth factor 10 (FGF10) to downstream of MRPS30, with a minimum squared correlation among SNPs of 0.80 (Figure 1). FGF10/FGFR2 signaling [26–29] could be relevant to these associations, though there is a recombination hotspot between the FGF10 gene and the 5p12 SNPs studied here.
Our analyses suggest that interactions of these two SNPs with WHI clinical trial interventions lead to a detectable increase in the ability to distinguish breast cancer cases from controls. Note, however, that AUC values in this context may be optimistic in view of our procedure for identifying SNPs of interest. Moreover, since the interactions identified in the study have yet to be confirmed by replication studies, the increase in AUC detected here is of exploratory nature as well. Also note that AUCs estimated here tend to be somewhat low due to age matching in the case-control sample.
When our combined test is separated into its constituents, one observes nominally significant evidence of interaction of MRPS30 SNP rs7705343 with three of the four WHI interventions. Given the manner in which we ranked SNPs, these analyses (Tables 2 and 3) should be regarded as exploratory and such interactions will need to be confirmed separately. Unfortunately, other clinical trial data are not available for this purpose, and confirmation in observational study settings will involve the challenge of reliable ascertainment of the relevant hormonal or dietary exposures, and will need to be carried out in a case-control rather than case-only model. Hence, quite large numbers of cases and controls will be needed, as may be accessible through cohort consortia.
It is interesting to see a significant interaction of rs7705343 with E-alone with the estimated intervention OR below 1.0 for the GG genotype, and an insignificant interaction of rs7705343 with E+P with the estimated intervention OR greater than 1 for the GG genotype. Few interactions with study subject characteristics have been suggested for E+P , with FGFR2 SNP rs3750817 as a possible exception . In contrast, interactions with several subject characteristics have been identified for E-alone, including family history of breast cancer, benign breast disease , and again FGFR2 SNP rs3750817 . A possible explanation is that the progestin in E+P tends to overwhelm the minor variations in hormone therapy hazard ratios that would otherwise occur, giving rise to a strong and fairly uniform risk elevation.
Study strengths include its nesting within the randomized controlled WHI clinical trial, implying randomization assignments that are known and that are statistically independent of genotype and the related ability to use case-only analyses for intervention testing. Other strengths of the study include the use of pre-diagnostic blood specimens, collected and stored according to a standardized protocol, and quality-controlled SNP genotyping.
A limitation of the study is that the average age at enrollment was 63 years in the WHI controlled trials, with many women well past menopause at enrollment. We have reported, in combined clinical trials and observational studies analyses, higher breast cancer hazard ratios for E+P and E-alone among women who first use these preparations soon after the menopause, compared to those using them later [30, 31]. Hence, the magnitude of the odds ratios shown here may be lower than would apply to typical hormone therapy users.
Simultaneous consideration of overall association and intervention interaction point to genomic regions in the vicinity of FGFR2 and MRPS30 genes as relevant to breast cancer risk among postmenopausal women. Moreover, subregions that were not otherwise the focus of interest, in the vicinity of SNPs rs3750817 and rs7705343, were identified as worthy of further study by virtue of suggestive interactions with hormonal and dietary interventions. These analyses represent an early step in assessing the role of genotype by 'environment' interactions to help explain familial breast cancer patterns, or as a contributor to risk discrimination.
area under the receiver operating characteristic curve
- CaD trial:
calcium and vitamin D versus placebo supplementation component
- DM trial:
low-fat dietary modification versus usual diet component
low-fat dietary modification trial in the subset of women having baseline percentage of energy from fat in the upper quartile
- E-alone trial:
estrogen versus placebo
- E+P trial:
estrogen plus progestin versus placebo
false discovery rate
- FGF10 :
fibroblast growth factor 10
- FGFR2 :
fibroblast growth factor receptor 2
- MRPS30 :
mitochondrial ribosomal protein S30
single nucleotide polymorphism
Women's Health Initiative.
Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, Hunter DJ, McCarthy MI, Ramos EM, Cardon LR, Chakravarti A, Cho JH, Guttmacher AE, Kong A, Kruglyak L, Mardis E, Rotimi CN, Slatkin M, Valle D, Whittemore AS, Boehnke M, Clark AG, Eichler EE, Gibson G, Haines JL, Mackay TFC, McCarroll SA, Visscher PM: Finding the missing heritability of complex diseases. Nature. 2009, 461: 747-753. 10.1038/nature08494.
Pharoah PDP, Antoniou AC, Easton DF, Ponder BAJ: Polygenes, risk prediction, and targeted prevention of breast cancer. N Engl J Med. 2008, 358: 2796-2803. 10.1056/NEJMsa0708739.
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.
Devilee P, Rookus MA: A tiny step closer to personalized risk prediction for breast cancer. N Engl J Med. 2010, 362: 1043-1045. 10.1056/NEJMe0912474.
Prentice RL, Shaw PA, Bingham SA, Beresford SA, Caan B, Neuhouser ML, Patterson RE, Stefanick ML, Satterfield S, Thomson CA, Snetselaar L, Thomas A, Tinker LF: Biomarker-calibrated energy and protein consumption and increased cancer risk among postmenopausal women. Am J Epidemiol. 2009, 169: 977-989. 10.1093/aje/kwp008.
Self SG, Longton G, Kopecky KJ, Liang KY: On estimating HLA/disease association with application to a study of aplastic anemia. Biometrics. 1991, 47: 53-61. 10.2307/2532495.
Piegorsch WW, Weinberg CR, Taylor JA: Non-hierarchical logistic models and case-only designs for assessing susceptibility in population-based case-control studies. Stat Med. 1994, 13: 153-162. 10.1002/sim.4780130206.
Vittinghoff E, Bauer DC: Case-only analysis of treatment-covariate iterations in clinical trials. Biometrics. 2006, 62: 769-776. 10.1111/j.1541-0420.2006.00511.x.
The Women's Health Initiative Study Group: Design of the Women's Health Initiative clinical trial and observational study. Control Clin Trials. 1998, 19: 61-109.
Prentice RL, Anderson GL: The Women's Health Initiative: lessons learned. Ann Rev Public Health. 2007, 29: 131-150.
Rossouw JE, Anderson GL, Prentice RL, LaCroix AZ, Kooperberg C, Stefanick ML, Jackson RD, Beresford SA, Howard BV, Johnson KC, Kotchen JM, Ockene J, Writing Group for the Women's Health Initiative Investigators: Risks and benefits of estrogen plus progestin in healthy postmenopausal women: principal results from the Women's Health Initiative randomized controlled trial. JAMA. 2002, 288: 321-333. 10.1001/jama.288.3.321.
Chlebowski RT, Hendrix SL, Langer RD, Stefanick ML, Gass M, Lane D, Rodabough RJ, Gilligan MA, Cyr MG, Thomson CA, Khandekar J, Petrovitch H, McTiernan A, WHI Investigators: Influence of estrogen plus progestin on breast cancer and mammography in healthy postmenopausal women: the Women's Health Initiative randomized trial. JAMA. 2003, 289: 3243-3253. 10.1001/jama.289.24.3243.
Anderson GL, Limacher M, Assaf AR, Bassford T, Beresford SA, Black H, Bonds D, Brunner R, Brzyski R, Caan B, Chlebowski R, Curb D, Gass M, Hays J, Heiss G, Hendrix S, Howard BV, Hsia J, Hubbell A, Jackson R, Johnson KC, Judd H, Kotchen JM, Kuller L, LaCroix AZ, Lane D, Langer RD, Lasser N, Lewis CE, Manson J, Women's Health Initiative Steering Committee, et al: Effects of conjugated equine estrogen in postmenopausal women with hysterectomy: the Women's Health Initiative randomized controlled trial. JAMA. 2004, 291: 1701-1712.
Stefanick ML, Anderson GL, Margolis KL, Hendrix SL, Rodabough RJ, Paskett ED, Lane DS, Hubbell FA, Assaf AR, Sarto GE, Schenken RS, Yasmeen S, Lessin L, Chlebowski RT, WHI Investigators: Effects of conjugated equine estrogens on breast cancer and mammography screening in postmenopausal women with hysterectomy. JAMA. 2006, 295: 1647-1657. 10.1001/jama.295.14.1647.
Prentice RL, Caan B, Chlebowski RT, Patterson R, Kuller LH, Ockene JK, Margolis KL, Limacher MC, Manson JE, Parker LM, Paskett E, Phillips L, Robbins J, Rossouw JE, Sarto GE, Shikany JM, Stefanick ML, Thomson CA, Van Horn L, Vitolins MZ, Wactawski-Wende J, Wallace RB, Wassertheil-Smoller S, Whitlock E, Yano K, Adams-Campbell L, Anderson GL, Assaf AR, Beresford SA, Black HR, et al: Low-fat dietary pattern and risk of invasive breast cancer: the Women's Health Initiative randomized controlled dietary modification trial. JAMA. 2006, 295: 629-642. 10.1001/jama.295.6.629.
Chlebowski RT, Johnson KC, Kooperberg C, Pettinger M, Wactawski-Wende J, Rohan T, Rossouw J, Lane D, O'Sullivan MJ, Yasmeen S, Hiatt RA, Shikany JM, Vitolins M, Khandekar J, Hubbell FA, Women's Health Initiative Investigators: Calcium plus vitamin D supplementation and the risk of breast cancer. J Natl Cancer Inst. 2008, 100: 1581-1591. 10.1093/jnci/djn360.
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.
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, SEARCH collaborators, 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, et al: Genome-wide association study identifies novel breast cancer susceptibility loci. Nature. 2007, 447: 1087-1093. 10.1038/nature05887.
Hays JL, Hunt JR, Hubbell FA, Anderson GL, Limacher M, Allen C, Rossouw JE: The women's health initiative recruitment methods and results. Ann Epidmiol. 2002, 13 (9 Supp): S18-S77.
Saccone SS, Hinrich AL, Saccone NL, Chase GA, Konvicka K, Madden PA, Breslau N, Johnson EO, Hatsukami D, Pomerleau O, Swan GE, Goate AM, Rutter J, Bertelsen S, Fox L, Fugman D, Martin NG, Montgomery GW, Wang JC, Ballinger DG, Rice JP, Bierut LJ: Cholinergic nicotinic receptor genes implicated in a nicotine dependence association study targeting 348 candidate genes with 3713 SNPs. Hum Mol Genet. 2006, 16: 36-49. 10.1093/hmg/ddl438.
Gail MH, Constantino JP, Bryant J, Croyle R, Freedman L, Helzlsouer K, Vogel V: Weighing the risks and benefits of tamoxifen treatment for preventing breast cancer. J Natl Cancer Inst. 1999, 91: 1829-1846. 10.1093/jnci/91.21.1829.
Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D: Principal component analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006, 38: 904-909. 10.1038/ng1847.
Dai J, LeBlanc M, Kooperberg C, Prentice RL: On Two-stage Hypothesis Testing Procedures via Asymptotically Independent Statistics. 2010, University of Washington Biostat Working Paper Series 366, [http://www.bepress.com/uwbiostat/]
Prentice RL, Huang Y, Hinds DA, Peters U, Pettinger M, Cox DR, Beilharz E, Chlebowski RT, Rossouw JE, Caan B, Ballinger DG: Variation in the FGFR2 gene and the effects of postmenopausal hormone therapy on invasive breast cancer. Cancer Epidemiol Biomarkers Prev. 2009, 18: 3079-3085. 10.1158/1055-9965.EPI-09-0611.
Prentice RL, Huang Y, Hinds DA, Peters U, Cox DR, Beilharz E, Chlebowski RT, Rossouw JE, Caan B, Ballinger DG: Variation in the FGFR2 gene and the effects of a low-fat dietary pattern on invasive breast cancer. Cancer Epidemiol Biomarkers Prev. 2010, 19: 74-79. 10.1158/1055-9965.EPI-09-0663.
Theodorou V, Boer M, Weigelt B, Jonkers J, van der Valk M, Hilkens J: FGF10 is an oncogene activated by MMTV insertional mutagenesis in mouse mammary tumors and overexpressed in a subset of human breast carcinomas. Oncogene. 2004, 23: 6047-6055. 10.1038/sj.onc.1207816.
Katoh M: Cancer genomics and genetics of FGFR2 (Review). Int J Oncol. 2008, 33: 233-237.
Lü J, Izvolsky KI, Qian J, Cardoso WV: Identification of FGF10 targets in the embryonic lung epithelium during bud morphogenesis. J Biol Chem. 2005, 280: 4834-4841.
Nomura S, Yoshitomi H, Takano S, Shida T, Kobayashi S, Ohtsuka M, Kimura F, Shimizu H, Yoshidome H, Kato A, Miyazaki M: FGF10/FGFR2 signal induces cell migration and invasion in pancreatic cancer. Br J Cancer. 2008, 99: 305-313. 10.1038/sj.bjc.6604473.
Prentice RL, Chlebowski RT, Stefanick ML, Manson JE, Langer RD, Pettinger M, Hendrix SL, Hubbell FA, Kooperberg C, Kuller LH, Lane DS, McTiernan A, O'Sullivan MJ, Rossouw JE, Anderson GL: Conjugated equine estrogens and breast cancer risk in the Women's Health Initiative clinical trial and observational study. Am J Epidemiol. 2008, 167: 1407-1415. 10.1093/aje/kwn090.
Prentice RL, Chlebowski RT, Stefanick ML, Manson JE, Pettinger M, Hendrix SL, Hubbell FA, Kooperberg C, Kuller LH, Lane DS, McTiernan A, Jo O'Sullivan M, Rossouw JE, Anderson GL: Estrogen plus progestin therapy and breast cancer in recently postmenopausal women. Am J Epidemiol. 2008, 167: 1207-1216. 10.1093/aje/kwn044.
Decisions concerning study design, data collection and analysis, interpretation of the results, the preparation of the manuscript, or the decision to submit the manuscript for publication resided with committees composed of WHI investigators that included NHLBI representatives. Program Office: (National Heart, Lung, and Blood Institute, Bethesda, MD, USA) Jacques Rossouw, Shari Ludlam, Joan McGowan, Leslie Ford, and Nancy Geller. Clinical Coordinating Center: (Fred Hutchinson Cancer Research Center, Seattle, WA, USA) Ross Prentice, Garnet Anderson, Andrea LaCroix, Charles L Kooperberg; (Medical Research Labs, Highland Heights, KY, USA) Evan Stein; (University of California at San Francisco, San Francisco, CA, USA) Steven Cummings. Clinical Centers: (Albert Einstein College of Medicine, Bronx, NY, USA) Sylvia Wassertheil-Smoller; (Baylor College of Medicine, Houston, TX, USA) Haleh Sangi-Haghpeykar; (Brigham and Women's Hospital, Harvard Medical School, Boston, MA, USA) JoAnn E Manson; (Brown University, Providence, RI, USA) Charles B Eaton; (Emory University, Atlanta, GA, USA) Lawrence S Phillips; (Fred Hutchinson Cancer Research Center, Seattle, WA, USA) Shirley Beresford; (George Washington University Medical Center, Washington, DC, USA) Lisa Martin; (Los Angeles Biomedical Research Institute at Harbor- UCLA Medical Center, Torrance, CA, USA) Rowan Chlebowski; (Kaiser Permanente Center for Health Research, Portland, OR, USA) Erin LeBlanc; (Kaiser Permanente Division of Research, Oakland, CA, USA) Bette Caan; (Medical College of Wisconsin, Milwaukee, WI, USA) Jane Morley Kotchen; (MedStar Research Institute/Howard University, Washington, DC, USA) Barbara V Howard; (Northwestern University, Chicago/Evanston, IL, USA) Linda Van Horn; (Rush Medical Center, Chicago, IL, USA) Henry Black; (Stanford Prevention Research Center, Stanford, CA, USA). Marcia L Stefanick; (State University of New York at Stony Brook, Stony Brook, NY, USA) Dorothy Lane; (The Ohio State University, Columbus, OH, USA) Rebecca Jackson; (University of Alabama at Birmingham, Birmingham, AL, USA) Cora E Lewis; (University of Arizona, Tucson/Phoenix, AZ, USA) Cynthia A Thomson; (University at Buffalo, Buffalo, NY, USA) Jean Wactawski-Wende; (University of California at Davis, Sacramento, CA, USA) John Robbins; (University of California at Irvine, CA, USA) F Allan Hubbell; (University of California at Los Angeles, Los Angeles, CA, USA) Lauren Nathan; (University of California at San Diego, LaJolla/Chula Vista, CA, USA) Robert D Langer; (University of Cincinnati, Cincinnati, OH, USA) Margery Gass; (University of Florida, Gainesville/Jacksonville, FL, USA) Marian Limacher; (University of Hawaii, Honolulu, HI, USA) J David Curb; (University of Iowa, Iowa City/Davenport, IA, USA) Robert Wallace; (University of Massachusetts/Fallon Clinic, Worcester, MA, USA) Judith Ockene; (University of Medicine and Dentistry of New Jersey, Newark, NJ, USA) Norman Lasser; (University of Miami, Miami, FL, USA) Mary Jo O'Sullivan; (University of Minnesota, Minneapolis, MN, USA) Karen Margolis; (University of Nevada, Reno, NV) Robert Brunner; (University of North Carolina, Chapel Hill, NC, USA) Gerardo Heiss; (University of Pittsburgh, Pittsburgh, PA, USA) Lewis Kuller; (University of Tennessee Health Science Center, Memphis, TN, USA) Karen C Johnson; (University of Texas Health Science Center, San Antonio, TX, USA) Robert Brzyski; (University of Wisconsin, Madison, WI, USA) Gloria E Sarto; (Wake Forest University School of Medicine, Winston-Salem, NC, USA) Mara Vitolins; (Wayne State University School of Medicine/Hutzel Hospital, Detroit, MI, USA) Michael S Simon. Women's Health Initiative Memory Study: (Wake Forest University School of Medicine, Winston-Salem, NC, USA) Sally Shumaker. This work was supported by the National Heart, Lung, and Blood Institute, National Institutes of Health, US Department of Health and Human Services [contracts HHSN268200764314C, N01WH22110, 24152, 32100-2, 32105-6, 32108-9, 32111-13, 32115, 32118-19, 32122, 42107-26, 42129-32, and 44221]. Clinical Trials Registration: ClinicalTrials.gov identifier, NCT00000611. The work of Dr Prentice was partially supported by grants CA53996 and CA148065 from the National Cancer Institute.
RTC reports receiving consulting fees from AstraZeneca, Novartis, Pfizer, and Eli Lilly, lecture fees from AstraZeneca and Novartis, and grant support from Amgen. No other potential conflict of interest relevant to this article was reported.
All authors were involved in development and/or critical review and revision of the manuscript. Additionally, DB, DH, DC, and EB had primary responsibility for project genotyping; YH, DH and RP had primary responsibility for data analysis; RC, JR, AM, TR and RP had responsibility for clinical data; and DB, UP and RP had primary administrative responsibility for this research project.
An erratum to this article is available at http://dx.doi.org/10.1186/gm318.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.