Common variants in breast cancer risk loci predispose to distinct tumor subtypes.

BACKGROUND: Genome-wide association studies (GWAS) have identified multiple common breast cancer susceptibility variants. Many of these variants have differential associations by estrogen receptor (ER) status, but how these variants relate with other tumor features and intrinsic molecular subtypes is unclear. METHODS: Among 106,571 invasive breast cancer cases and 95,762 controls of European ancestry with data on 173 breast cancer variants identified in previous GWAS, we used novel two-stage polytomous logistic regression models to evaluate variants in relation to multiple tumor features (ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and grade) adjusting for each other, and to intrinsic-like subtypes. RESULTS: Eighty-five of 173 variants were associated with at least one tumor feature (false discovery rate < 5%), most commonly ER and grade, followed by PR and HER2. Models for intrinsic-like subtypes found nearly all of these variants (83 of 85) associated at p < 0.05 with risk for at least one luminal-like subtype, and approximately half (41 of 85) of the variants were associated with risk of at least one non-luminal subtype, including 32 variants associated with triple-negative (TN) disease. Ten variants were associated with risk of all subtypes in different magnitude. Five variants were associated with risk of luminal A-like and TN subtypes in opposite directions. CONCLUSION: This report demonstrates a high level of complexity in the etiology heterogeneity of breast cancer susceptibility variants and can inform investigations of subtype-specific risk prediction.

Mendelian randomisation study of smoking exposure in relation to breast cancer risk.

BACKGROUND: Despite a modest association between tobacco smoking and breast cancer risk reported by recent epidemiological studies, it is still equivocal whether smoking is causally related to breast cancer risk. METHODS: We applied Mendelian randomisation (MR) to evaluate a potential causal effect of cigarette smoking on breast cancer risk. Both individual-level data as well as summary statistics for 164 single-nucleotide polymorphisms (SNPs) reported in genome-wide association studies of lifetime smoking index (LSI) or cigarette per day (CPD) were used to obtain MR effect estimates. Data from 108,420 invasive breast cancer cases and 87,681 controls were used for the LSI analysis and for the CPD analysis conducted among ever-smokers from 26,147 cancer cases and 26,072 controls. Sensitivity analyses were conducted to address pleiotropy. RESULTS: Genetically predicted LSI was associated with increased breast cancer risk (OR 1.18 per SD, 95% CI: 1.07-1.30, P = 0.11 × 10-2), but there was no evidence of association for genetically predicted CPD (OR 1.02, 95% CI: 0.78-1.19, P = 0.85). The sensitivity analyses yielded similar results and showed no strong evidence of pleiotropic effect. CONCLUSION: Our MR study provides supportive evidence for a potential causal association with breast cancer risk for lifetime smoking exposure but not cigarettes per day among smokers.

Assessment of stromal tumor infiltrating lymphocytes and immunohistochemical features in invasive micropapillary breast carcinoma with long-term outcomes.

PURPOSE: We studied the long-term outcomes of invasive micropapillary carcinoma (IMPCs) of the breast in relation to stromal tumor infiltrating lymphocytes (sTILs), prognostic biomarkers and clinicopathological features. METHODS: Stage I-III IMPCs treated with upfront surgery at our institution (January 2000 and December 2016) were included. Central pathology review was performed and sTILs (including zonal distribution and hot spot analysis) and tumor-associated plasma cells (TAPC) were evaluated. Expression of P53, BCL2, FOXP3, and WT1, which are variably linked to breast cancer prognosis, was measured by immunohistochemistry using tissue microarrays. Time-to-event endpoints were distant recurrence free interval (DRFI) and breast cancer-specific survival (BCSS). RESULTS: We included 111 patients of whom 59% were pure IMPCs. Standard clinicopathological features were comparable between pure and non-pure IMPCs. Overall, the mean sTILs level was 20% with higher proportion of sTILs present at the invasive front. There were no significant differences between pure- and non-pure IMPCs in sTILs levels, nor in the spatial distribution of the hot spot regions or in the distribution of TAPC. Higher sTILs correlated with worse DRFI (HR = 1.55; p = 0.0172) and BCSS (HR = 2.10; p < 0.001). CONCLUSIONS: Clinicopathological features, geographical distribution of sTILs and TAPC are similar between pure and non-pure IMPCs. Despite a high proportion of grade 3 tumors and lymph node involvement, we observed a low rate of distant recurrences and breast cancer-related death in this cohort of stage I-III IMPCs treated with primary surgery. Caution in interpretation of the observed prognostic correlations is required given the very low number of events, warranting validation in other cohorts.

Differences in polygenic score distributions in European ancestry populations: implications for breast cancer risk prediction.

The 313-variant polygenic risk score (PRS313) provides a promising tool for breast cancer risk prediction. However, evaluation of the PRS313 across different European populations which could influence risk estimation has not been performed. Here, we explored the distribution of PRS313 across European populations using genotype data from 94,072 females without breast cancer, of European-ancestry from 21 countries participating in the Breast Cancer Association Consortium (BCAC) and 225,105 female participants from the UK Biobank. The mean PRS313 differed markedly across European countries, being highest in south-eastern Europe and lowest in north-western Europe. Using the overall European PRS313 distribution to categorise individuals leads to overestimation and underestimation of risk in some individuals from south-eastern and north-western countries, respectively. Adjustment for principal components explained most of the observed heterogeneity in mean PRS. Country-specific PRS distributions may be used to calibrate risk categories in individuals from different countries.

Genome-wide association study of treatment-related toxicity two years following radiotherapy for breast cancer.

BACKGROUND AND PURPOSE: Up to a quarter of breast cancer patients treated by surgery and radiotherapy experience clinically significant toxicity. If patients at high risk of adverse effects could be identified at diagnosis, their treatment could be tailored accordingly. This study was designed to identify common single nucleotide polymorphisms (SNPs) associated with toxicity two years following whole breast radiotherapy. MATERIALS AND METHODS: A genome-wide association study (GWAS) was performed in 1,640 breast cancer patients with complete SNP, clinical, treatment and toxicity data, recruited across 18 European and US centres into the prospective REQUITE cohort study. Toxicity data (CTCAE v4.0) were collected at baseline, end of radiotherapy, and annual follow-up. A total of 7,097,340 SNPs were tested for association with the residuals of toxicity endpoints, adjusted for clinical, treatment co-variates and population substructure. RESULTS: Quantile-quantile plots showed more associations with toxicity above the p < 5 × 10-5 level than expected by chance. Eight SNPs reached genome-wide significance. Nipple retraction grade ≥ 2 was associated with the rs188287402 variant (p = 2.80 × 10-8), breast oedema grade ≥ 2 with rs12657177 (p = 1.12 × 10-10), rs75912034 (p = 1.12 × 10-10), rs145328458 (p = 1.06 × 10-9) and rs61966612 (p = 1.23 × 10-9), induration grade ≥ 2 with rs77311050 (p = 2.54 × 10-8) and rs34063419 (p = 1.21 × 10-8), and arm lymphoedema grade ≥ 1 with rs643644 (p = 3.54 × 10-8). Heritability estimates across significant endpoints ranged from 25% to 39%. Our study did not replicate previously reported SNPs associated with breast radiation toxicity at the pre-specified significance level. CONCLUSIONS: This GWAS for long-term breast radiation toxicity provides further evidence for significant association of common SNPs with distinct toxicity endpoints.

Gene-Environment Interactions Relevant to Estrogen and Risk of Breast Cancer: Can Gene-Environment Interactions Be Detected Only among Candidate SNPs from Genome-Wide Association Studies?

In this study we aim to examine gene-environment interactions (GxEs) between genes involved with estrogen metabolism and environmental factors related to estrogen exposure. GxE analyses were conducted with 1970 Korean breast cancer cases and 2052 controls in the case-control study, the Seoul Breast Cancer Study (SEBCS). A total of 11,555 SNPs from the 137 candidate genes were included in the GxE analyses with eight established environmental factors. A replication test was conducted by using an independent population from the Breast Cancer Association Consortium (BCAC), with 62,485 Europeans and 9047 Asians. The GxE tests were performed by using two-step methods in GxEScan software. Two interactions were found in the SEBCS. The first interaction was shown between rs13035764 of NCOA1 and age at menarche in the GE|2df model (p-2df = 1.2 × 10-3). The age at menarche before 14 years old was associated with the high risk of breast cancer, and the risk was higher when subjects had homozygous minor allele G. The second GxE was shown between rs851998 near ESR1 and height in the GE|2df model (p-2df = 1.1 × 10-4). Height taller than 160 cm was associated with a high risk of breast cancer, and the risk increased when the minor allele was added. The findings were not replicated in the BCAC. These results would suggest specificity in Koreans for breast cancer risk.