- Research article
- Open Access
- Open Peer Review
Breast density and polymorphisms in genes coding for CYP1A2 and COMT: the Multiethnic Cohort
BMC Cancervolume 7, Article number: 30 (2007)
Mammographic density is a strong predictor of breast cancer risk and is increased by hormone replacement therapy (HRT). Some associations with genetic polymorphisms in enzymes involved in estrogen metabolism have been described. This cross-sectional analysis examined the relation between mammographic density and the CYP1A2*1F and COMT Val 58 Met polymorphisms among 332 breast cancer cases and 254 controls in the Hawaii component of the Multiethnic Cohort.
Mammographic density, before diagnosis in cases, was quantified by using a validated computer-assisted method. Blood samples were genotyped by standard PCR/RFLP methods. Adjusted mean percent density was calculated by genotype using mixed models with the unstructured covariance option.
A positive association between the C allele in the CYP1A2*1F gene and percent density, but not the dense area, was suggested (p = 0.11). The relation was limited to controls (p = 0.045), postmenopausal women not using HRT (p = 0.08), and normal weight subjects (p = 0.046). We did not observe any relation between the COMT Val 58 Met polymorphism and breast density.
The lack of an association between the CYP1A2 genotype and the size of the dense areas suggests an effect on the non-dense, i.e., fatty breast tissue. The discrepancies among studies may be due to differential susceptibility; changes in enzyme activity as a result of the CYP1A2*1F polymorphism may influence breast tissue differently depending on hormonal status. Larger studies with the ability to look at interactions would be useful to elucidate the influence of genetic variation in CYP1A2 and COMT on the risk of developing breast cancer.
Mammographic density has been shown to be independently associated with breast cancer risk . Although a link between mammographic density and circulating estrogen levels was only reported in one out of four reports [2–5], postmenopausal hormone replacement therapy (HRT) appears to increase breast density in observational and experimental studies [6, 7]. In addition, associations with genetic polymorphisms in enzymes involved in estrogen metabolism have been explored [8–14]. A previous study in Hawaii with predominately premenopausal women showed lower mammographic densities for women with the C allele in the CYP1A2*1F gene and the Met allele in the COMT gene . The CC genotype for CYP1A2*1F was also significantly associated with lower serum estradiol levels during the luteal phase . The particular polymorphisms are thought to lead to lower enzyme activity and they had been associated with breast cancer risk in some reports [16, 17]. Also, higher CYP1A2 activity as assessed by urinary caffeine metabolites was associated with higher mammographic densities among postmenopausal women with high malondialdehyde (MDA) levels, an indicator of lipid peroxidation . In a case-control study nested within the Multiethnic Cohort (MEC), postmenopausal women with at least one C allele in the CYP1A2*1F gene had a lower risk of breast cancer than women with the common alleles , but a case-control study in Shanghai with women aged 25–64 years did not observe this association . As for the Met allele in the COMT gene, higher breast density was found for postmenopausal HRT users in one report , lower breast density in women not using HRT in another study , lower breast density in pre- but not postmenopausal women , higher increase in breast density in response to estrogens during an intervention , and no association among postmenopausal women . To clarify the conflicting results related to these two polymorphisms, we linked the aforementioned case-control study  with a mammographic case-control study also nested within the MEC that had collected multiple mammograms over time  and examined the association between breast density and the CYP1A2*1F and the COMT Val 58 Met polymorphisms among the overlapping subjects.
The subjects for this analysis were originally recruited for two separate case-control studies [18, 20] nested within the MEC study . Both studies were approved by the Human Subject Committee at the University of Hawaii. All participants completed an informed consent form. Controls were randomly selected from the cohort and frequency matched by age and ethnicity for both studies. The genetic polymorphism case-control study included 1,339 breast cancer cases and 1,370 controls from Hawaii and Los Angeles . Incident breast cancer cases since 1995 were identified through the rapid reporting system of the Hawaii Tumor Registry and the Los Angeles County Cancer Surveillance Program. Cases and controls who agreed to participate in the study (74% for cases and 66% for controls) donated a blood sample. For the mammographic density case-control study , all incident cases diagnosed with breast cancer in Hawaii between cohort entry and December 2000 were eligible. Of them, 51% agreed to be in the study and for 44% of the eligible subjects a mammogram was located. The final sample size was 607 cases and 667 controls. After linking the breast density study with the genetic polymorphism study, 575 subjects had mammographic and genetic information available.
To assess breast density, the cranio-caudal views for both sides of several mammographic examinations were retrieved . All mammograms for cases were performed before treatment for breast cancer was initiated and only five mammograms were taken at the time of diagnosis. The films were scanned after removing personal identifiers and quantified with a computer-assisted method by one reader . The total and the dense areas of the breast were estimated, percent density was calculated as the ratio of the dense to the total area of the breast, and the values for the right and the left breast were averaged. For a subset of mammograms read in duplicate, the intraclass correlation coefficient was 0.96 for the size of the dense areas and 0.996 for the total breast area, resulting in a value of 0.974 for percent density.
DNA was extracted from white blood cells and genotyped by PCR/RFLP methods for the COMT Val 58 Met and CYP1A2*1F (A-154C) polymorphisms as described previously .
All statistical analyses were performed using SAS version 9 (Cary, NC) . In order to estimate adjusted mean percent density, mixed models with the unstructured covariance option were applied while adjusting for age at mammogram, BMI, case-control status, HRT use and type (estrogen vs. combined therapy), ethnicity, parity, age at menarche, age at first live birth, family history of breast cancer, and menopausal status, where appropriate. A square root transformation was applied to percent density because the variables were not normally distributed. We repeated some models using the log-transformed size of the dense areas as the dependent variable. Mixed models were used to account for the multiple density measures per person. The availability of repeated density assessments increases the statistical power of the analysis because it decreases the intraindividual variation in density measures . Age, menopausal status, and HRT use were included as time dependent variables; all other variables were only assessed once. Genotypes were coded according to the number of variant alleles as none, one, and two to test for a codominant model; coding for a recessive and a dominant model was also performed. Adjusted mean percent densities were calculated according to the number of variant alleles among all women and stratified by case status, ethnicity, menopausal status, and HRT use that were hypothesized a priori as possible modifiers. We also tested for interaction effects of genotype with these variables. Logistic regression was used to estimate the risk of breast cancer related to the two genetic polymorphisms.
Due to the study design, nearly half of the 575 subjects were breast cancer cases (Table 1). On average, three mammograms were included per woman, but the number of scanned films was higher for cases than for controls. At the time of the first mammogram, the mean age was 57.1 ± 9.3 years and 67% of subjects were postmenopausal. By the time of the last mammogram, on average 4.7 years later, 80% of women were postmenopausal. The median time between first mammogram and breast cancer diagnosis was 6.6 years. Japanese women constituted the largest ethnic group, followed by Caucasians and Native Hawaiians. Approximately two-thirds of women had used HRT at the time of at least one mammogram, but overall 44% of mammograms were performed without HRT use. A higher proportion of controls than cases reported ERT use, whereas combined HRT use was more common among cases. As expected, mean percent density was significantly higher among cases than controls; 36.7% for cases vs. 30.6% for controls.
Adjusted percent densities were significantly higher with an increasing number of variant alleles in the CYP1A2*1F polymorphism among all women; the difference between the AA and CC genotypes was 5.1% (p for trend = 0.11) (Table 2). The respective p-values for the recessive and the dominant model were 0.23 and 0.18 with density differences of 4.1% and 2.2% between groups (data not shown). Although we detected no significant interaction effects between the CYP1A2*1F polymorphism and case status, menopausal status, and BMI (p = 0.13, 0.38, and 0.37, respectively), we noted a number of differences after stratification by these variables. The association between the CYP1A2*1F allele and percent density was stronger among controls than cases. The difference between subjects with AA and CC genotype was 8.9% (p = 0.045), while there was no significant difference in breast density by genotype among cases (p = 0.86). Given the small number of mammograms among premenopausal women, the 9.6% difference between the two homozygous groups was not significant (p = 0.38). Among postmenopausal women, the 9.1% difference for mammograms taken without HRT use was close to significance (p = 0.08), whereas the difference among HRT users was only 2.5% (p = 0.48). Similarly, the positive association was limited to women with a BMI of less than 25 kg/m2 (9.4%, p = 0.046), whereas percent density did not vary by genotype among overweight and obese women. Stratification by ethnicity showed non-significant higher mean densities for Japanese CC carriers but no relation among Caucasians. The mean densities for the Japanese AA, AC, and CC carriers were 32.9%, 31.7%, and 39.4% (p = 0.44); the respective values for Caucasians were 26.0%, 27.9%, and 25.1% (p = 0.60).
Genotypes for COMT Val 58 Met were not associated with percent density in the entire population or in any subgroup (Table 2). The density difference between homozygous genotypes was less than 1% in the entire population. The respective p-values for the recessive and the dominant model were 0.83 and 0.70. No interaction effect of COMT Val 58 Met with case status, menopausal status, or BMI was suggested (p = 0.71, 0.90, and 0.20, respectively).
When we used the size of the dense areas as dependent variable, we found no difference in mean densities by genotype for CYP1A2 and COMT (p = 0.85 and 0.98).
Compared to women with the AA genotype, the risk estimates for breast cancer associated with one or two C alleles in the CYP1A2 gene were 0.86 (95% CI = 0.60–1.23) and 0.62 (95% CI = 0.30–1.30), respectively. For the COMT Val 58 Met polymorphism, the OR was 0.93 (95% CI = 0.64–1.35) for one Met allele and 0.61 (95% CI = 0.37–1.02) for two Met alleles as compared to the women with two Val alleles.
In this analysis of longitudinal mammographic data, we observed a weak association between percent density and the number of C alleles in the CYP1A2*1F gene, but the relation was limited to controls, postmenopausal women not using HRT, and normal weight women. The fact that the size of the dense areas was not associated with the CYP1A2 genotype suggests that the relation may be due to an effect on the non-dense, i.e., fatty breast tissue. We did not observe any relation between percent density and the Met polymorphism in the COMT gene.
In comparison, a Canadian study reported mean percent densities of 14% for the lowest and 31% for the highest quartile of CYP1A2 activity (p for trend = 0.01) among postmenopausal women not using HRT . The current findings are inconsistent with studies conducted in Hawaii that described lower percent densities  and estrogen levels  in premenopausal carriers of the C allele, as well as a lower breast cancer risk among postmenopausal women with at least one C allele for CYP1A2*1F . However, the estimated risks of 0.86 (95% CI = 0.60–1.23) for heterozygotes and 0.62 (95% CI = 0.30–1.30) for women with the CC genotype in the current study are very similar to the risk estimates in the larger case-control study (OR = 1.0, 0.9 and 0.7 for the AA, AC, and CC genotypes; p for trend = 0.03) that also included Latina and African-American women in Los Angeles. The findings with regard to the COMT Val 58 Met polymorphism are also inconsistent. Studies in Hawaii  and in Canada  suggested lower breast densities for premenopausal carriers of the Met allele. In contrast, one recent study found no relation among postmenopausal women  and two reports described associations for subgroups of postmenopausal women with the Met allele: higher breast density for HRT users  and lower breast density in women not using HRT .
These discrepant results for both genotypes may be due to chance given the limited sample size, but other explanations are also plausible. The small change in hormone levels due to a variant allele may not be observable, especially in light of the weak evidence for a relation between breast density and circulating estrogen levels [2–5]. Our observation that the relation between the CYP1A2 genotype and percent density was probably due to an inverse association with the non-dense, i.e., fatty breast tissue, may offer a plausible explanation for the discrepancy with the findings on breast cancer risk . As a speculation, we propose that one or two C alleles may be protective against breast cancer because the CYP1A2 genotype might influence the amount of adipose tissue in the breast and elsewhere. As apparent from the sex differences in body fat distribution, steroid hormones influence fat tissue, most likely through binding to sex steroid hormone receptors that regulate the production of adipocytokines in adipose tissue [25, 26]. Therefore, differences in circulating estrogen metabolite patterns due to altered enzyme activity might modify the amount of adipose tissue in the breast. Alternatively, the causality of the association between genotype and amount of non-dense tissue may be reversed. Because there is evidence that the expression level of CYP1A2 is upregulated by estrogens , postmenopausal women with more body fat and higher endogenous estrogen production may have higher CYP1A2 activity. The genetic polymorphism may only affect postmenopausal women with low body fat, low endogenous estrogens, and relatively low CYP1A2 activity, whereas the higher estrogen levels and CYP1A2 activity in women who use HRT , will develop breast cancer , or are overweight  may outweigh any changes in hormone patterns due to the CYP1A2 variant alleles. This idea of differential susceptibility would explain why the effects of the CYP1A2 genotype were only apparent in certain subgroups and why our previous findings in premenopausal women were so different [10, 15]. In addition, a recent experimental study showed that, among all CYP enzymes tested, CYP1A2 has the highest 2-hydroxylation (2-OH) activity on estrone, especially at low concentrations . The positive association between the CYP1A2 polymorphism and breast density among controls, women not taking HRT, and normal weight women could be due to the higher levels of 2-OH estrone, a metabolite with low estrogenic activity, among common allele carriers than among women with the variant allele resulting in greater mammographic densities among the latter.
Strengths of this study were our prospective data from a cohort study that minimized recall bias and the retrieval of longitudinal mammographic data. There are a number of serious limitations in our study. For instance, the response rates and proportions of ethnic groups by case status were not optimal [18, 20]. Selection bias due to the availability of mammograms may have led to a higher proportion of women at high risk for breast cancer than in the general population. Although we were able to show a strong association between breast cancer diagnosis and breast density in the mammographic density study , the scanned mammograms may not have covered the critical time of a woman's life. Unmeasured confounders, such as lipid peroxidation and intake of antioxidants [11, 32], may be partially responsible for the contradictory findings. Although CYP1A2 activity contributes to estrogen metabolism, polymorphisms in other genes should be considered by themselves and in combination. It is not quite clear yet how much the CYP1A2*1F allele influences actual enzyme activity. When caffeine metabolites were measured, enzyme activity was higher among carriers of the CYP1A2*1F allele in one , but not in another report . Given the stronger association of percent density with progesterone than estrogens [3, 7], polymorphisms related to progestin metabolism may provide additional insights.
This study among a subset of MEC members suggests higher percent densities for women carrying at least one C allele in the CYP1A2 gene than for those with two A alleles but no association between breast density and the Met allele in the COMT gene. It appears that the CYP1A2 genotype affects the non-dense area of the breast, but a possible mechanism of action remains to be defined. Given the weak association between percent density and estrogens, the alteration in enzyme activity due to the CYP1A2*1F polymorphism may be too small to influence breast density, in particular among women who are exposed to relatively high hormone levels, such as HRT users, overweight postmenopausal women, and women who will develop breast cancer. Until larger studies with the ability to look at gene-gene and gene-environment interactions are completed, it will remain unclear to what extent the variant alleles for CYP1A2 and COMT affect a woman's lifetime risk of developing breast cancer.
Boyd NF, Lockwood GA, Byng JW, Tritchler DL, Yaffe MJ: Mammographic densities and breast cancer risk. Cancer Epidemiol Biomarkers Prev. 1998, 7: 1133-1144.
Boyd NF, Stone J, Martin LJ, Jong R, Fishell E, Yaffe M, Hammond G, Minkin S: The association of breast mitogens with mammographic densities. Br J Cancer. 2002, 87: 876-882. 10.1038/sj.bjc.6600537.
Noh JJ, Maskarinec G, Pagano I, Cheung LW, Stanczyk FZ: Mammographic densities and circulating hormones: a cross-sectional study in premenopausal women. Breast. 2006, 15: 20-28. 10.1016/j.breast.2005.04.014.
Tamimi RM, Hankinson SE, Colditz GA, Byrne C: Endogenous sex hormone levels and mammographic density among postmenopausal women. Cancer Epidemiol Biomarkers Prev. 2005, 14: 2641-2647. 10.1158/1055-9965.EPI-05-0558.
Greendale GA, Palla SL, Ursin G, Laughlin GA, Crandall C, Pike MC, Reboussin BA: The association of endogenous sex steroids and sex steroid binding proteins with mammographic density: results from the Postmenopausal Estrogen/Progestin Interventions Mammographic Density Study. Am J Epidemiol. 2005, 162: 826-834. 10.1093/aje/kwi286.
Laya MB, Gallagher JC, Schreiman JS, Larson EB, Watson P, Weinstein L: Effect of postmenopausal hormonal replacement therapy on mammographic density and parenchymal pattern. Radiology. 1995, 196: 433-437.
Greendale GA, Reboussin BA, Slone S, Wasilauskas C, Pike MC, Ursin G: Postmenopausal hormone therapy and change in mammographic density. J Natl Cancer Inst. 2003, 95: 30-37.
Haiman CA, Bernstein L, Berg D, Ingles SA, Salane M, Ursin G: Genetic determinants of mammographic density. Breast Cancer Res. 2002, 4: R5-10.1186/bcr434.
Haiman CA, Hankinson SE, De Vivo I, Guillemette C, Ishibe N, Hunter DJ, Byrne C: Polymorphisms in steroid hormone pathway genes and mammographic density. Breast Cancer Res Treat. 2003, 77: 27-36. 10.1023/A:1021112121782.
Maskarinec G, Lurie G, Williams AE, Le Marchand L: An investigation of mammographic density and gene variants in healthy women. Int J Cancer. 2004, 112: 683-10.1002/ijc.20448.
Hong CC, Tang BK, Rao V, Agarwal S, Martin L, Tritchler D, Yaffe M, Boyd NF: Cytochrome P450 1A2 (CYP1A2) activity, mammographic density, and oxidative stress: a cross-sectional study. Breast Cancer Res. 2004, 6: R338-R351. 10.1186/bcr797.
Hong CC, Thompson HJ, Jiang C, Hammond GL, Tritchler D, Yaffe M, Boyd NF: Val158Met Polymorphism in catechol-O-methyltransferase gene associated with risk factors for breast cancer. Cancer Epidemiol Biomarkers Prev. 2003, 12: 838-847.
Warren R, Skinner J, Sala E, Denton E, Dowsett M, Folkerd E, Healey CS, Dunning A, Doody D, Ponder B, Luben RN, Day NE, Easton D: Associations among mammographic density, circulating sex hormones, and polymorphisms in sex hormone metabolism genes in postmenopausal women. Cancer Epidemiol Biomarkers Prev. 2006, 15: 1502-1508. 10.1158/1055-9965.EPI-05-0828.
Lord SJ, Mack WJ, Van Den BD, Pike MC, Ingles SA, Haiman CA, Wang W, Parisky YR, Hodis HN, Ursin G: Polymorphisms in genes involved in estrogen and progesterone metabolism and mammographic density changes in women randomized to postmenopausal hormone therapy: results from a pilot study. Breast Cancer Res. 2005, 7: R336-R344. 10.1186/bcr999.
Lurie G, Maskarinec G, Kaaks R, Stanczyk FZ, Le ML: Association of genetic polymorphisms with serum estrogens measured multiple times during a 2-year period in premenopausal women. Cancer Epidemiol Biomarkers Prev. 2005, 14: 1521-1527. 10.1158/1055-9965.EPI-04-0746.
Dunning AM, Healey CS, Pharoah PD, Teare MD, Ponder BA, Easton DF: A systematic review of genetic polymorphisms and breast cancer risk. Cancer Epidemiol Biomarkers Prev. 1999, 8: 843-854.
Coughlin SS, Piper M: Genetic polymorphisms and risk of breast cancer. Cancer Epidemiol Biomarkers Prev. 1999, 8: 1023-1032.
Le Marchand L, Donlon T, Kolonel LN, Henderson BE, Wilkens LR: Estrogen metabolism-related genes and breast cancer risk: the multiethnic cohort study. Cancer Epidemiol Biomarkers Prev. 2005, 14: 1998-2003. 10.1158/1055-9965.EPI-05-0076.
Long JR, Egan KM, Dunning L, Shu XO, Cai Q, Cai H, Dai Q, Holtzman J, Gao YT, Zheng W: Population-based case-control study of AhR (aryl hydrocarbon receptor) and CYP1A2 polymorphisms and breast cancer risk. Pharmacogenet Genomics. 2006, 16: 237-243.
Maskarinec G, Pagano I, Lurie G, Wilkens LR, Kolonel LN: Mammographic Density and Breast Cancer Risk: The Multiethnic Cohort. Am J Epidemiol. 2005, 162: 743-752. 10.1093/aje/kwi270.
Kolonel LN, Henderson BE, Hankin JH, Nomura AMY, Wilkens LR, Pike MC, Stram DO, Monroe KR, Earle ME, Nagamine FS: A multiethnic cohort in Hawaii and Los Angeles: baseline characteristics. Am J Epidemiol. 2000, 151: 346-357.
Byng JW, Boyd NF, Fishell E, Jong RA, Yaffe MJ: The quantitative analysis of mammographic densities. Phys Med Biol. 1994, 39: 1629-1638. 10.1088/0031-9155/39/10/008.
Inc. SASI: SAS OnlineDoc 9.1.2. 2004, Cary, NC, SAS Institute Inc.
Singer JD, Willett JB: Applied longitudinal data analysis: Modeling change and event occurrence. 2003, Oxford, Oxford University Press
Mayes JS, Watson GH: Direct effects of sex steroid hormones on adipose tissues and obesity. Obes Rev. 2004, 5: 197-216. 10.1111/j.1467-789X.2004.00152.x.
Laughlin GA, Barrett-Connor E, May S: Sex-specific association of the androgen to oestrogen ratio with adipocytokine levels in older adults: the Rancho Bernardo Study. Clin Endocrinol (Oxf). 2006, 65: 506-513. 10.1111/j.1365-2265.2006.02624.x.
Tsuchiya Y, Nakajima M, Yokoi T: Cytochrome P450-mediated metabolism of estrogens and its regulation in human. Cancer Lett. 2005, 227: 115-124. 10.1016/j.canlet.2004.10.007.
Lee SA, Ross RK, Pike MC: An overview of menopausal oestrogen-progestin hormone therapy and breast cancer risk. Br J Cancer. 2005
The Endogenous Hormones and Breast Cancer Collaborative Group: Endogenous sex hormones and breast cancer in postmenopausal women: reanalysis of nine prospective studies. J Natl Cancer Inst. 2002, 94: 606-616.
Key TJ, Appleby PN, Reeves GK, Roddam A, Dorgan JF, Longcope C, Stanczyk FZ, Stephenson HE, Falk RT, Miller R, Schatzkin A, Allen DS, Fentiman IS, Key TJ, Wang DY, Dowsett M, Thomas HV, Hankinson SE, Toniolo P, Akhmedkhanov A, Koenig K, Shore RE, Zeleniuch-Jacquotte A, Berrino F, Muti P, Micheli A, Krogh V, Sieri S, Pala V, Venturelli E, Secreto G, Barrett-Connor E, Laughlin GA, Kabuto M, Akiba S, Stevens RG, Neriishi K, Land CE, Cauley JA, Kuller LH, Cummings SR, Helzlsouer KJ, Alberg AJ, Bush TL, Comstock GW, Gordon GB, Miller SR, Longcope C: Body mass index, serum sex hormones, and breast cancer risk in postmenopausal women. J Natl Cancer Inst. 2003, 95: 1218-1226.
Cribb AE, Knight MJ, Dryer D, Guernsey J, Hender K, Tesch M, Saleh TM: Role of polymorphic human cytochrome P450 enzymes in estrone oxidation. Cancer Epidemiol Biomarkers Prev. 2006, 15: 551-558. 10.1158/1055-9965.EPI-05-0801.
Wiseman A: Oestrogen-receptors (ER) are likely to be promiscuous: wider role for oestrogens and mimics. Med Hypotheses. 2005, 65: 760-765. 10.1016/j.mehy.2005.04.019.
Sachse C, Brockmoller J, Bauer S, Roots I: Functional significance of a C-->A polymorphism in intron 1 of the cytochrome P450 CYP1A2 gene tested with caffeine. Br J Clin Pharmacol. 1999, 47: 445-449. 10.1046/j.1365-2125.1999.00898.x.
Takata K, Saruwatari J, Nakada N, Nakagawa M, Fukuda K, Tanaka F, Takenaka S, Mihara S, Marubayashi T, Nakagawa K: Phenotype-genotype analysis of CYP1A2 in Japanese patients receiving oral theophylline therapy. Eur J Clin Pharmacol. 2006, 62: 23-28. 10.1007/s00228-005-0057-z.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/7/30/prepub
The genetic polymorphism study was supported by an award from the U.S. Army Medical Research Materiel Command under DAMD (17-00-1-0283) to Dr. Loic Le Marchand and by a grant from the National Cancer Institute to Dr. Brian Henderson (R01 CA63464). The mammographic case-control study was funded by a grant from the National Cancer Institute (R01 CA 85265). The Multiethnic Cohort (MEC) has been supported by USPHS (National Cancer Insitute) Grant R37 CA 54281 (PI: Dr. L.N. Kolonel). We are grateful to Drs. Laurence N. Kolonel and Brian H. Henderson for accesss to the MEC study. We would also like to thank Drs. Lynne Wilkens and Ian Pagano for their advice on the statistical analysis.
The author(s) declare that they have no competing interests.
YT performed the statistical analysis and drafted the manuscript. GM designed the study and helped to perform the statistical analysis and draft the manuscript. LM designed the study, supervised the genetic analysis and critically commented on the drafted manuscript. GM and LM obtained funding for the previous projects. All authors read and approved the final manuscript.