Genetic Variants at Chromosomes 2q35, 5p12, 6q25.1, 10q26.13, and 16q12.1 Influence the Risk of Breast Cancer in Men

Summary

Male breast cancer accounts for approximately 1% of all breast cancer. To date, risk factors for male breast cancer are poorly defined, but certain risk factors and genetic features appear common to both male and female breast cancer. Genome-wide association studies (GWAS) have recently identified common single nucleotide polymorphisms (SNPs) that influence female breast cancer risk; 12 of these have been independently replicated. To examine if these variants contribute to male breast cancer risk, we genotyped 433 male breast cancer cases and 1,569 controls. Five SNPs showed a statistically significant association with male breast cancer: rs13387042 (2q35) (odds ratio (OR)  = 1.30, p = 7.98×10⁻⁴), rs10941679 (5p12) (OR = 1.26, p = 0.007), rs9383938 (6q25.1) (OR = 1.39, p = 0.004), rs2981579 (FGFR2) (OR = 1.18, p = 0.03), and rs3803662 (TOX3) (OR = 1.48, p = 4.04×10⁻⁶). Comparing the ORs for male breast cancer with the published ORs for female breast cancer, three SNPs—rs13387042 (2q35), rs3803662 (TOX3), and rs6504950 (COX11)—showed significant differences in ORs (p<0.05) between sexes. Breast cancer is a heterogeneous disease; the relative risks associated with loci identified to date show subtype and, based on these data, gender specificity. Additional studies of well-defined patient subgroups could provide further insight into the biological basis of breast cancer development.

Introduction

Breast cancer does not exclusively affect females. Around 300 men in the UK and 1,900 men in the US are diagnosed with the disease each year [1]. The average age at incidence of male breast cancer is somewhat different to that seen for female breast cancer, with the disease typically affecting men 5–10 years later than women. Perhaps because male breast cancer is not common, few risk factors have been demonstrated to influence disease risk, but tentative associations with obesity, lack of exercise, excess alcohol consumption, gynaecomastia, past benign breast disease, past liver disease, infertility, diabetes and exposure to ionising radiation have been suggested [2], [3].

Investigation of susceptibility genes for male breast cancer has been limited. It has however been shown that approximately 10% of men with breast cancer carry BRCA2 mutations, while mutations in BRCA1 are exceedingly rare [4]. The relative risk of breast cancer in men associated with BRCA2 mutations is high [5]. Recently the CHEK2 1100delC variant has been found to give a 10-fold risk of male breast cancer independent of BRCA1 or BRCA2 [6]. Mutations in these genes are rare in the general population and it is likely that much of the genetic contribution to female breast cancer risk can be attributed to the co-inheritance of multiple low risk common variants [7]. Recent genome-wide association studies (GWAS) have shown associations between single nucleotide polymorphisms (SNPs) mapping to a dozen or more loci and female breast cancer risk in European populations, each conferring odds ratios (ORs) of 1.04–1.43 [8]–[14]. To explore the possibility that the same risk variants influence male breast cancer risk, we conducted a case-control study of male breast cancer, genotyping 12 SNPs annotating the loci that have the strongest and most consistent associations with female breast cancer.

Materials and Methods

457 cases of male breast cancer were recruited in a population-based case-control study of the genetic, environmental and behavioral causes of male breast cancer being conducted in England and Wales. Potential cases were all men resident in these countries aged 18–79 with newly diagnosed breast cancer since January 1^st, 2005, identified through notifications by treatment centres and systematic regular listings of cases from regional cancer registries. 98% of cases for whom registry data has been received have been histologically confirmed. The median age at diagnosis of cases was 65.5 years (interquartile range: 59–72).

A total of 1608 unmatched controls were available for genotyping; 535 men were ascertained through our ongoing breast cancer studies and a further 1073 were healthy male and female individuals from the UK Genetic Lung Cancer Predisposition Study (GELCAPS) [15]. The decision to include a second control set was made a priori, with the aim of increasing statistical power. We saw no evidence for an effect of control group on the overall effect estimate for each SNP. Collection of blood samples from all subjects was undertaken with informed consent and relevant ethical review committee approval.

DNA was extracted from venous blood samples using conventional methodologies and quantified by Picogreen (Invitrogen, Carlsbad CA). SNPs were chosen for analysis on the basis of validated associations with female breast cancer from recent GWAS [8]–[14]. Genotyping of rs11249433, rs13387042, rs4973768, rs10941679, rs16886165, rs9383938, rs13281615, rs865686, rs2981579, rs3817198, rs3803662 and rs6504950 was performed by allele-specific PCR using KASPar chemistry (Kbioscience, Hertfordshire, UK). Each DNA plate contained 5% sample duplication to assess genotyping concordance between duplicate pairs. We attempted to genotype 2119 samples (including duplicates, n = 54) and excluded samples (n = 49; 11 cases, 34 controls and four members of a duplicate pair) in which no-calls were observed for two or more SNPs. Genotyping QC statistics were therefore computed on 2070 samples (Figure S1). Final locus and sample completion rates were >99.9%. The mean genotype concordance between duplicate pairs was 99.8%. We excluded a further 18 subjects due to self-reported non-European ancestry (13 cases and 5 controls). No SNP genotypes showed significant deviation from the proportions expected under Hardy-Weinberg equilibrium in controls (Table S1).

ORs and 95% confidence intervals (CI) were calculated using unconditional logistic regression. The odds ratio for each SNP was determined by fitting multiplicative and unconstrained genetic models. P-values were computed from likelihood ratio test statistics. Case-only unconditional logistic regression was used to test the significance of association with age at diagnosis. Deviation of genotype proportions from Hardy-Weinberg equilibrium was assessed in controls using an exact test [16]. To compare formally the ORs in males with the equivalent published ORs for female disease, we assumed both sets of ORs were log-normally distributed. Then under the null hypothesis that the OR in males is equal to the OR in females, the difference between the estimated log ORs is normally distributed with mean zero and variance equal to the sum of the squared standard errors of the two estimates. From this we obtained a χ² statistic for each comparison (1 degree of freedom [d.f.]) and from the sum of the χ² statistics a global test for all comparisons (12 d.f.). Statistical analyses were performed using the Genotype Libraries and Utilities (GLU) package (http://code.google.com/p/glu-genetics) and R [17].

Results/Discussion

433 male breast cancer cases and 1569 controls were successfully genotyped according to our predefined QC criteria. The majority of cases were diagnosed with invasive breast cancer (n = 399 (92%)) while a further 31 (7%) were ductal carcinoma in situ. Three cases (<1%) were of unknown histology. Table 1 shows the OR for male breast cancer associated with each of the 12 SNPs previously reported to be associated with female breast cancer risk. For five SNPs, rs13387042 (2q35), rs10941679 (5p12), rs9383938 (6q25.1), rs2981579 (FGFR2) and rs3803662 (TOX3), the risk allele for female breast cancer was associated with increased risk of male breast cancer (p<0.05). Two SNPs, rs13387042 (2q35) and rs3803662 (TOX3), remained significant below the Bonferroni adjusted threshold for independent tests of p<4.12×10⁻³.

Comparing OR_male estimates with those for female breast cancer (OR_female) there were two SNPs, rs13387042 (2q35) and rs3803882 (TOX3) for which the OR_male was significantly higher than the OR_female, albeit not after adjusting for multiple testing (rs13387042, OR_male:OR_female p = 0.03; rs3803882, OR_male:OR_female p = 0.04; Table 2). rs3803662 (TOX3) showed the strongest association with male breast cancer (OR_male  = 1.48; 95% CI 1.26–1.75, p = 4.04×10⁻⁶) with an excess relative risk that was more than twice the female estimate (OR_female  = 1.20; 95% CI 1.16–1.24) [9]. Similarly, the excess risk conferred by rs13387042 (2q35) in males (OR_male  = 1.30; 95% CI 1.11–1.51, p = 7.98×10⁻⁴) was more than double that observed in females (OR_female  = 1.12; 95% CI 1.09–1.15) [11].

For one SNP (rs6504950, COX11) the OR_male was in the opposite direction to that reported for female breast cancer (OR_male  = 0.90; 95% CI 0.76–1.06, OR_female  = 1.05; 95% CI 1.03–1.07) [8], [9] and was inconsistent with the female estimate (OR_male:OR_female p = 0.04; Table 2). For the other nine SNPs that we tested the OR_male estimates were consistent with the OR_female estimates. Comparing the combined estimates of all 12 SNPs, however, there was nominal evidence that the male ORs differed from the female ORs (p = 0.03; Table 2).

The frequency of female breast tumors that are estrogen receptor (ER) positive varies, particularly according to menopausal status at diagnosis [18]. Based on a sample of almost 3,000 patients the proportion is typically between 64% and 79% [18]. In contrast, male breast tumors, tend to be overwhelmingly ER-positive (>90%) [19]. In the current study estrogen receptor status was known for 251 male breast cancer cases, 246 (98%) of which had ER-positive tumors. For nine of the 12 SNPs that we genotyped, OR_female estimates stratified according to ER status have been reported for Caucasian populations (Tables S2a and S2b). In females, the OR for ER-positive disease is stronger than the OR for ER-negative disease for all nine of these loci and this difference is significant for all but two of them (rs16886165 (MAP3KI) and rs3817198 (LSP1) [8], [11], [20]. Given the predominance of ER-positive tumors in male disease we also compared the OR_male with the OR_female for ER-positive disease (Table S2a). There was nominally significant evidence overall that the male ORs differed from those for ER-positive female disease (p = 0.05). We also tested for a difference between the OR_male estimates for these nine SNPs and the OR_female estimates for ER-negative disease (Table S2b); there was stronger evidence of a difference (p = 0.01).

Finally, we assessed the relationship between genotype and age at onset of male breast cancer (Table S3) for each of the 12 loci. There was no evidence for a trend with age at diagnosis.

We have shown, for the first time that common genetic variants influence susceptibility to male breast cancer. Furthermore we have demonstrated that for at least a subset of known susceptibility loci the risk allele for female breast cancer is also associated with increased risk of disease in males. To our knowledge these 433 male breast cancer cases represent the largest single series to date; despite this, we lacked power to detect modest relative risks for all but the most common variants. For example we had only 40% power to detect an OR of 1.15 for a variant with a minor allele frequency of 30% at a significance level of 5%. The lack of a statistically significant association with male breast cancer risk for seven of the 12 SNPs that we tested may, therefore, simply reflect a lack of power.

Notably, for two of the three SNPs for which the OR_male was inconsistent with the OR_female (rs13387042 (2q35) and rs3803882 (TOX3)) the association in males was stronger than that in females. While the OR_male estimates were slightly closer to the ORs for ER-positive disease in females, it is noticeable that these are the two SNPs that show the largest effects on ER-negative disease risk in females. Although the significance of this observation, if any, is not yet clear, our data on male breast cancer alongside the published associations with female breast cancer, clearly implicate the 2q35 and 16q12.1 loci in the aetiology of breast cancer, irrespective of gender and tumor pathology.

Given that the majority of female breast cancer risk loci identified to date demonstrate a degree of specificity for ER-positive or ER-negative disease [8], [11], [20], [21] it seems likely that subtype specific GWAS will lead to the identification of additional risk loci. Our analyses suggest that GWAS of male breast cancer may also lead to the identification of novel breast cancer risk loci in males and that these should provide further insight into the biological basis of male and female breast cancer development.

Supporting Information

Attachment 1

Attachment 2

Attachment 3

Attachment 4