High miR-30 Expression Associates with Improved Breast Cancer Patient Survival and Treatment Outcome.

Deregulated miRNA expression has been suggested in several stages of breast cancer pathogenesis. We have studied the miR-30 family, in particular miR-30d, in relation to breast cancer patient survival and treatment outcomes. With tumor specimens from 1238 breast cancer patients, we analyzed the association of miR-30d expression with tumor characteristics with the 5-year occurrence of breast cancer-specific death or distant metastasis (BDDM), and with 10-year breast cancer survival (BCS). We conducted a two-stage drug-screen to investigate the impact of miR-30 family members (miR-30a-30e) on sensitivity to doxorubicin and lapatinib in six breast cancer cell lines HCC1937, HCC1954, MDA-MB-361, MCF7, MDA-MB-436 and CAL-120, using drug sensitivity scores (DSS) to compare the miR-30 family mimics to their specific inhibitors. The study was complemented with Ingenuity Pathway Analysis (IPA) with the METABRIC data. We found that while high miR-30d expression is typical for aggressive tumors, it predicts better metastasis-free (pBDDM = 0.035, HR = 0.63, 95% CI = 0.4-0.9) and breast cancer-specific survival (pBCS = 0.018, HR = 0.61, 95% CI = 0.4-0.9), especially in HER2-positive (pBDDM = 0.0009), ER-negative (pBDDM = 0.003), p53-positive (pBDDM = 0.011), and highly proliferating (pBDDM = 0.0004) subgroups, and after adjuvant chemotherapy (pBDDM = 0.035). MiR-30d predicted survival independently of standard prognostic markers (pBDDM = 0.0004). In the drug-screening test, the miR-30 family sensitized the HER2-positive HCC1954 cell line to lapatinib (p < 10-2) and HCC1937, MDA-MB-361, MDA-MB-436 and CAL120 to doxorubicin (p < 10-4) with an opposite impact on MCF7. According to the pathway analysis, the miR-30 family has a suppressive effect on cell motility and metastasis in breast cancer. Our results suggest prognostic and predictive potential for the miR-30 family, which warrants further investigation.

Association of germline variation with the survival of women with BRCA1/2 pathogenic variants and breast cancer.

Germline genetic variation has been suggested to influence the survival of breast cancer patients independently of tumor pathology. We have studied survival associations of genetic variants in two etiologically unique groups of breast cancer patients, the carriers of germline pathogenic variants in BRCA1 or BRCA2 genes. We found that rs57025206 was significantly associated with the overall survival, predicting higher mortality of BRCA1 carrier patients with estrogen receptor-negative breast cancer, with a hazard ratio 4.37 (95% confidence interval 3.03-6.30, P = 3.1 × 10-9). Multivariable analysis adjusted for tumor characteristics suggested that rs57025206 was an independent survival marker. In addition, our exploratory analyses suggest that the associations between genetic variants and breast cancer patient survival may depend on tumor biological subgroup and clinical patient characteristics.

Meta-analysis of three genome-wide association studies identifies two loci that predict survival and treatment outcome in breast cancer.

The majority of breast cancers are driven by the female hormone oestrogen via oestrogen receptor (ER) alpha. ER-positive patients are commonly treated with adjuvant endocrine therapy, however, resistance is a common occurrence and aside from ER-status, no unequivocal predictive biomarkers are currently in clinical use. In this study, we aimed to identify constitutional genetic variants influencing breast cancer survival among ER-positive patients and specifically, among endocrine-treated patients. We conducted a meta-analysis of three genome-wide association studies comprising in total 3,136 patients with ER-positive breast cancer of which 2,751 had received adjuvant endocrine therapy. We identified a novel locus (rs992531 at 8p21.2) associated with reduced survival among the patients with ER-positive breast cancer (P = 3.77 × 10-8). Another locus (rs7701292 at 5q21.3) was associated with reduced survival among the endocrine-treated patients (P = 2.13 × 10-8). Interaction analysis indicated that the survival association of rs7701292 is treatment-specific and independent of conventional prognostic markers. In silico functional studies suggest plausible biological mechanisms for the observed survival associations and a functional link between the putative target genes of the rs992531 and rs7701292 (RHOBTB2 and RAB9P1, respectively). We further explored the genetic interaction between rs992531 and rs7701292 and found a significant, treatment-specific interactive effect on survival among ER-positive, endocrine-treated patients (hazard ratio = 6.97; 95% confidence interval, 1.79-27.08, Pinteraction= 0.036). This is the first study to identify a genetic interaction that specifically predicts treatment outcome. These findings may provide predictive biomarkers based on germ line genotype informing more personalized treatment selection.

Family history influences the tumor characteristics and prognosis of breast cancers developing during postmenopausal hormone therapy.

Long term use of postmenopausal hormone therapy (HT) has been reported to increase breast cancer risk. On the other hand, observational studies suggest that breast cancers diagnosed during HT may have a more favorable prognosis. While family history is a risk factor for breast cancer, and genetic factors also influence prognosis, the role of family history in combination with HT use has been little studied. We investigated the relationship between HT, family history, and prognosis in 584 (267 exposed) familial and 952 (460 exposed) non-familial breast cancer cases, using three survival end points: death from breast cancer (BCS), distant disease free survival (DDFS), and local recurrence free survival (LRFS). Among non-familial cases, HT was associated with better BCS (HR 0.63, 95% CI 0.41-0.94; p = 0.025), and DDFS (HR 0.58, 95% CI 0.40-0.85; p = 0.005), with a consistent but not statistically significant effect in LRFS. This effect was not seen in familial cases (HR > 1.0), and family history was found to interact with HT in BCS (p(interaction) = 0.0067) (BC-death) and DDFS (p(interaction) = 0.0070). There was phenotypic heterogeneity between HT-associated tumors in familial and non-familial cases, particularly on estrogen receptor (ER) status, although the interaction between HT and family history appears to be at least partially independent of these markers (p = 0.0370 after adjustment for standard prognostic factors). If confirmed by further studies, our results suggest that family history should be taken into consideration in clinical counseling before beginning a HT regimen.

Two missense mutations in KCNQ1 cause pituitary hormone deficiency and maternally inherited gingival fibromatosis.

Familial growth hormone deficiency provides an opportunity to identify new genetic causes of short stature. Here we combine linkage analysis with whole-genome resequencing in patients with growth hormone deficiency and maternally inherited gingival fibromatosis. We report that patients from three unrelated families harbor either of two missense mutations, c.347G>T p.(Arg116Leu) or c.1106C>T p.(Pro369Leu), in KCNQ1, a gene previously implicated in the long QT interval syndrome. Kcnq1 is expressed in hypothalamic GHRH neurons and pituitary somatotropes. Co-expressing KCNQ1 with the KCNE2 ß-subunit shows that both KCNQ1 mutants increase current levels in patch clamp analyses and are associated with reduced pituitary hormone secretion from AtT-20 cells. In conclusion, our results reveal a role for the KCNQ1 potassium channel in the regulation of human growth, and show that growth hormone deficiency associated with maternally inherited gingival fibromatosis is an allelic disorder with cardiac arrhythmia syndromes caused by KCNQ1 mutations.

TP53-based interaction analysis identifies cis-eQTL variants for TP53BP2, FBXO28, and FAM53A that associate with survival and treatment outcome in breast cancer.

TP53 overexpression is indicative of somatic TP53 mutations and associates with aggressive tumors and poor prognosis in breast cancer. We utilized a two-stage SNP association study to detect variants associated with breast cancer survival in a TP53-dependent manner. Initially, a genome-wide study (n = 575 cases) was conducted to discover candidate SNPs for genotyping and validation in the Breast Cancer Association Consortium (BCAC). The SNPs were then tested for interaction with tumor TP53 status (n = 4,610) and anthracycline treatment (n = 17,828). For SNPs interacting with anthracycline treatment, siRNA knockdown experiments were carried out to validate candidate genes.In the test for interaction between SNP genotype and TP53 status, we identified one locus, represented by rs10916264 (p(interaction) = 3.44 × 10-5; FDR-adjusted p = 0.0011) in estrogen receptor (ER) positive cases. The rs10916264 AA genotype associated with worse survival among cases with ER-positive, TP53-positive tumors (hazard ratio [HR] 2.36, 95% confidence interval [C.I] 1.45 - 3.82). This is a cis-eQTL locus for FBXO28 and TP53BP2; expression levels of these genes were associated with patient survival specifically in ER-positive, TP53-mutated tumors. Additionally, the SNP rs798755 was associated with survival in interaction with anthracycline treatment (p(interaction) = 9.57 × 10-5, FDR-adjusted p = 0.0130). RNAi-based depletion of a predicted regulatory target gene, FAM53A, indicated that this gene can modulate doxorubicin sensitivity in breast cancer cell lines.If confirmed in independent data sets, these results may be of clinical relevance in the development of prognostic and predictive marker panels for breast cancer.

Patient survival and tumor characteristics associated with CHEK2:p.I157T - findings from the Breast Cancer Association Consortium.

BACKGROUND: P.I157T is a CHEK2 missense mutation associated with a modest increase in breast cancer risk. Previously, another CHEK2 mutation, the protein truncating c.1100delC has been associated with poor prognosis of breast cancer patients. Here, we have investigated patient survival and characteristics of breast tumors of germ line p.I157T carriers. METHODS: We included in the analyses 26,801 European female breast cancer patients from 15 studies participating in the Breast Cancer Association Consortium. We analyzed the association between p.I157T and the clinico-pathological breast cancer characteristics by comparing the p.I157T carrier tumors to non-carrier and c.1100delC carrier tumors. Similarly, we investigated the p.I157T associated risk of early death, breast cancer-associated death, distant metastasis, locoregional relapse and second breast cancer using Cox proportional hazards models. Additionally, we explored the p.I157T-associated genomic gene expression profile using data from breast tumors of 183 Finnish female breast cancer patients (ten p.I157T carriers) (GEO: GSE24450). Differential gene expression analysis was performed using a moderated t test. Functional enrichment was investigated using the DAVID functional annotation tool and gene set enrichment analysis (GSEA). The tumors were classified into molecular subtypes according to the St Gallen 2013 criteria and the PAM50 gene expression signature. RESULTS: P.I157T was not associated with increased risk of early death, breast cancer-associated death or distant metastasis relapse, and there was a significant difference in prognosis associated with the two CHEK2 mutations, p.I157T and c.1100delC. Furthermore, p.I157T was associated with lobular histological type and clinico-pathological markers of good prognosis, such as ER and PR expression, low TP53 expression and low grade. Gene expression analysis suggested luminal A to be the most common subtype for p.I157T carriers and CDH1 (cadherin 1) target genes to be significantly enriched among genes, whose expression differed between p.I157T and non-carrier tumors. CONCLUSIONS: Our analyses suggest that there are fundamental differences in breast tumors of CHEK2:p.I157T and c.1100delC carriers. The poor prognosis associated with c.1100delC cannot be generalized to other CHEK2 mutations.

FANCM c.5101C>T mutation associates with breast cancer survival and treatment outcome.

Breast cancer (BC) is a heterogeneous disease, and different tumor characteristics and genetic variation may affect the clinical outcome. The FANCM c.5101C > T nonsense mutation in the Finnish population associates with increased risk of breast cancer, especially for triple-negative breast cancer patients. To investigate the association of the mutation with disease prognosis, we studied tumor phenotype, treatment outcome, and patient survival in 3,933 invasive breast cancer patients, including 101 FANCM c.5101C > T mutation carriers and 3,832 non-carriers. We also examined association of the mutation with nuclear immunohistochemical staining of DNA repair markers in 1,240 breast tumors. The FANCM c.5101C > T mutation associated with poor 10-year breast cancer-specific survival (hazard ratio (HR)=1.66, 95% confidence interval (CI) 1.09-2.52, p = 0.018), with a more pronounced survival effect among familial cases (HR = 2.93, 95% CI 1.5-5.76, p = 1.80 × 10-3 ). Poor disease outcome of the carriers was also found among the estrogen receptor (ER) positive subgroup of patients (HR = 1.8, 95% CI 1.09-2.98, p = 0.021). Reduced survival was seen especially among patients who had not received radiotherapy (HR = 3.43, 95% CI 1.6-7.34, p = 1.50 × 10-3 ) but not among radiotherapy treated patients (HR = 1.35, 95% CI 0.82-2.23, p = 0.237). Significant interaction was found between the mutation and radiotherapy (p = 0.040). Immunohistochemical analyses show that c.5101C > T carriers have reduced PAR-activity. Our results suggest that FANCM c.5101C > T nonsense mutation carriers have a reduced breast cancer survival but postoperative radiotherapy may diminish this survival disadvantage.

Polygenic risk score is associated with increased disease risk in 52 Finnish breast cancer families.

The risk of developing breast cancer is increased in women with family history of breast cancer and particularly in families with multiple cases of breast or ovarian cancer. Nevertheless, many women with a positive family history never develop the disease. Polygenic risk scores (PRSs) based on the risk effects of multiple common genetic variants have been proposed for individual risk assessment on a population level. We investigate the applicability of the PRS for risk prediction within breast cancer families. We studied the association between breast cancer risk and a PRS based on 75 common genetic variants in 52 Finnish breast cancer families including 427 genotyped women and pedigree information on ~4000 additional individuals by comparing the affected to healthy family members, as well as in a case-control dataset comprising 1272 healthy population controls and 1681 breast cancer cases with information on family history. Family structure was summarized using the BOADICEA risk prediction model. The PRS was associated with increased disease risk in women with family history of breast cancer as well as in women within the breast cancer families. The odds ratio (OR) for breast cancer within the family dataset was 1.55 [95 % CI 1.26-1.91] per unit increase in the PRS, similar to OR in unselected breast cancer cases of the case-control dataset (1.49 [1.38-1.62]). High PRS-values were informative for risk prediction in breast cancer families, whereas for the low PRS-categories the results were inconclusive. The PRS is informative in women with family history of breast cancer and should be incorporated within pedigree-based clinical risk assessment.

Heterogeneity of luminal breast cancer characterised by immunohistochemical expression of basal markers.

BACKGROUND: Luminal A breast cancer defined as hormone receptor positive and human epidermal growth factor receptor 2 (HER2) negative is known to be heterogeneous. Previous study showed that luminal A tumours with the expression of basal markers ((cytokeratin (CK) 5 or CK5/6) or epidermal growth factor receptor (EGFR)) were associated with poorer prognosis compared with those that stained negative for basal markers. Prompted by this study, we assessed whether tumour characteristics and risk factors differed by basal marker status within luminal A tumours. METHODS: We pooled 5040 luminal A cases defined by immunohistochemistry (4490 basal-negative ((CK5 (or CK5/6))- and EGFR-) and 550 basal-positive ((CK5 (or CK5/6+)) or EGFR+)) from eight studies participating in the Breast Cancer Association Consortium. Case-case comparison was performed using unconditional logistic regression. RESULTS: Tumour characteristics and risk factors did not vary significantly by the expression of basal markers, although results suggested that basal-positive luminal tumours tended to be smaller and node negative, and were more common in women with a positive family history and lower body mass index. CONCLUSIONS: Most established breast cancer risk factors were similar in basal-positive and basal-negative luminal A tumours. The non-significant but suggestive differences in tumour features and family history warrant further investigations.

Genetic variation in the immunosuppression pathway genes and breast cancer susceptibility: a pooled analysis of 42,510 cases and 40,577 controls from the Breast Cancer Association Consortium.

Immunosuppression plays a pivotal role in assisting tumors to evade immune destruction and promoting tumor development. We hypothesized that genetic variation in the immunosuppression pathway genes may be implicated in breast cancer tumorigenesis. We included 42,510 female breast cancer cases and 40,577 controls of European ancestry from 37 studies in the Breast Cancer Association Consortium (2015) with available genotype data for 3595 single nucleotide polymorphisms (SNPs) in 133 candidate genes. Associations between genotyped SNPs and overall breast cancer risk, and secondarily according to estrogen receptor (ER) status, were assessed using multiple logistic regression models. Gene-level associations were assessed based on principal component analysis. Gene expression analyses were conducted using RNA sequencing level 3 data from The Cancer Genome Atlas for 989 breast tumor samples and 113 matched normal tissue samples. SNP rs1905339 (A>G) in the STAT3 region was associated with an increased breast cancer risk (per allele odds ratio 1.05, 95 % confidence interval 1.03-1.08; p value = 1.4 × 10(-6)). The association did not differ significantly by ER status. On the gene level, in addition to TGFBR2 and CCND1, IL5 and GM-CSF showed the strongest associations with overall breast cancer risk (p value = 1.0 × 10(-3) and 7.0 × 10(-3), respectively). Furthermore, STAT3 and IL5 but not GM-CSF were differentially expressed between breast tumor tissue and normal tissue (p value = 2.5 × 10(-3), 4.5 × 10(-4) and 0.63, respectively). Our data provide evidence that the immunosuppression pathway genes STAT3, IL5, and GM-CSF may be novel susceptibility loci for breast cancer in women of European ancestry.

A polymorphism in the base excision repair gene PARP2 is associated with differential prognosis by chemotherapy among postmenopausal breast cancer patients.

BACKGROUND: Personalized therapy considering clinical and genetic patient characteristics will further improve breast cancer survival. Two widely used treatments, chemotherapy and radiotherapy, can induce oxidative DNA damage and, if not repaired, cell death. Since base excision repair (BER) activity is specific for oxidative DNA damage, we hypothesized that germline genetic variation in this pathway will affect breast cancer-specific survival depending on treatment. METHODS: We assessed in 1,408 postmenopausal breast cancer patients from the German MARIE study whether cancer specific survival after adjuvant chemotherapy, anthracycline chemotherapy, and radiotherapy is modulated by 127 Single Nucleotide Polymorphisms (SNPs) in 21 BER genes. For SNPs with interaction terms showing p<0.1 (likelihood ratio test) using multivariable Cox proportional hazard analyses, replication in 6,392 patients from nine studies of the Breast Cancer Association Consortium (BCAC) was performed. RESULTS: rs878156 in PARP2 showed a differential effect by chemotherapy (p=0.093) and was replicated in BCAC studies (p=0.009; combined analysis p=0.002). Compared to non-carriers, carriers of the variant G allele (minor allele frequency=0.07) showed better survival after chemotherapy (combined allelic hazard ratio (HR)=0.75, 95% 0.53-1.07) and poorer survival when not treated with chemotherapy (HR=1.42, 95% 1.08-1.85). A similar effect modification by rs878156 was observed for anthracycline-based chemotherapy in both MARIE and BCAC, with improved survival in carriers (combined allelic HR=0.73, 95% CI 0.40-1.32). None of the SNPs showed significant differential effects by radiotherapy. CONCLUSIONS: Our data suggest for the first time that a SNP in PARP2, rs878156, may together with other genetic variants modulate cancer specific survival in breast cancer patients depending on chemotherapy. These germline SNPs could contribute towards the design of predictive tests for breast cancer patients.

SNP-SNP interaction analysis of NF-κB signaling pathway on breast cancer survival.

In breast cancer, constitutive activation of NF-κB has been reported, however, the impact of genetic variation of the pathway on patient prognosis has been little studied. Furthermore, a combination of genetic variants, rather than single polymorphisms, may affect disease prognosis. Here, in an extensive dataset (n = 30,431) from the Breast Cancer Association Consortium, we investigated the association of 917 SNPs in 75 genes in the NF-κB pathway with breast cancer prognosis. We explored SNP-SNP interactions on survival using the likelihood-ratio test comparing multivariate Cox' regression models of SNP pairs without and with an interaction term. We found two interacting pairs associating with prognosis: patients simultaneously homozygous for the rare alleles of rs5996080 and rs7973914 had worse survival (HRinteraction 6.98, 95% CI=3.3-14.4, P=1.42E-07), and patients carrying at least one rare allele for rs17243893 and rs57890595 had better survival (HRinteraction 0.51, 95% CI=0.3-0.6, P = 2.19E-05). Based on in silico functional analyses and literature, we speculate that the rs5996080 and rs7973914 loci may affect the BAFFR and TNFR1/TNFR3 receptors and breast cancer survival, possibly by disturbing both the canonical and non-canonical NF-κB pathways or their dynamics, whereas, rs17243893-rs57890595 interaction on survival may be mediated through TRAF2-TRAIL-R4 interplay. These results warrant further validation and functional analyses.

Annexin A1 expression in a pooled breast cancer series: association with tumor subtypes and prognosis.

BACKGROUND: Annexin A1 (ANXA1) is a protein related with the carcinogenesis process and metastasis formation in many tumors. However, little is known about the prognostic value of ANXA1 in breast cancer. The purpose of this study is to evaluate the association between ANXA1 expression, BRCA1/2 germline carriership, specific tumor subtypes and survival in breast cancer patients. METHODS: Clinical-pathological information and follow-up data were collected from nine breast cancer studies from the Breast Cancer Association Consortium (BCAC) (n = 5,752) and from one study of familial breast cancer patients with BRCA1/2 mutations (n = 107). ANXA1 expression was scored based on the percentage of immunohistochemical staining in tumor cells. Survival analyses were performed using a multivariable Cox model. RESULTS: The frequency of ANXA1 positive tumors was higher in familial breast cancer patients with BRCA1/2 mutations than in BCAC patients, with 48.6 % versus 12.4 %, respectively; P <0.0001. ANXA1 was also highly expressed in BCAC tumors that were poorly differentiated, triple negative, EGFR-CK5/6 positive or had developed in patients at a young age. In the first 5 years of follow-up, patients with ANXA1 positive tumors had a worse breast cancer-specific survival (BCSS) than ANXA1 negative (HRadj = 1.35; 95 % CI = 1.05-1.73), but the association weakened after 10 years (HRadj = 1.13; 95 % CI = 0.91-1.40). ANXA1 was a significant independent predictor of survival in HER2+ patients (10-years BCSS: HRadj = 1.70; 95 % CI = 1.17-2.45). CONCLUSIONS: ANXA1 is overexpressed in familial breast cancer patients with BRCA1/2 mutations and correlated with poor prognosis features: triple negative and poorly differentiated tumors. ANXA1 might be a biomarker candidate for breast cancer survival prediction in high risk groups such as HER2+ cases.

Polymorphism at 19q13.41 Predicts Breast Cancer Survival Specifically after Endocrine Therapy.

PURPOSE: Although most patients with estrogen receptor (ER)-positive breast cancer benefit from endocrine therapies, a significant proportion do not. Our aim was to identify inherited genetic variations that might predict survival among patients receiving adjuvant endocrine therapies. EXPERIMENTAL DESIGN: We performed a meta-analysis of two genome-wide studies; Helsinki Breast Cancer Study, 805 patients, with 240 receiving endocrine therapy and Prospective study of Outcomes in Sporadic versus Hereditary breast cancer, 536 patients, with 155 endocrine therapy patients, evaluating 486,478 single-nucleotide polymorphisms (SNP). The top four associations from the endocrine treatment subgroup were further investigated in two independent datasets totaling 5,011 patients, with 3,485 receiving endocrine therapy. RESULTS: A meta-analysis identified a common SNP rs8113308, mapped to 19q13.41, associating with reduced survival among endocrine-treated patients [hazard ratio (HR), 1.69; 95% confidence interval (CI), 1.37-2.07; P = 6.34 × 10(-7)] and improved survival among ER-negative patients, with a similar trend in ER-positive cases not receiving endocrine therapy. In a multivariate analysis adjusted for conventional prognostic factors, we found a significant interaction between the rs8113308 and endocrine treatment, indicating a predictive, treatment-specific effect of the SNP rs8113308 on breast cancer survival, with the per-allele HR for interaction 2.16 (95% CI, 1.30-3.60; Pinteraction = 0.003) and HR = 7.77 (95% CI, 0.93-64.71) for the homozygous genotype carriers. A biologic rationale is suggested by in silico functional analyses. CONCLUSIONS: Our findings suggest carrying the rs8113308 rare allele may identify patients who will not benefit from adjuvant endocrine treatment.

Common germline polymorphisms associated with breast cancer-specific survival.

INTRODUCTION: Previous studies have identified common germline variants nominally associated with breast cancer survival. These associations have not been widely replicated in further studies. The purpose of this study was to evaluate the association of previously reported SNPs with breast cancer-specific survival using data from a pooled analysis of eight breast cancer survival genome-wide association studies (GWAS) from the Breast Cancer Association Consortium. METHODS: A literature review was conducted of all previously published associations between common germline variants and three survival outcomes: breast cancer-specific survival, overall survival and disease-free survival. All associations that reached the nominal significance level of P value <0.05 were included. Single nucleotide polymorphisms that had been previously reported as nominally associated with at least one survival outcome were evaluated in the pooled analysis of over 37,000 breast cancer cases for association with breast cancer-specific survival. Previous associations were evaluated using a one-sided test based on the reported direction of effect. RESULTS: Fifty-six variants from 45 previous publications were evaluated in the meta-analysis. Fifty-four of these were evaluated in the full set of 37,954 breast cancer cases with 2,900 events and the two additional variants were evaluated in a reduced sample size of 30,000 samples in order to ensure independence from the previously published studies. Five variants reached nominal significance (P <0.05) in the pooled GWAS data compared to 2.8 expected under the null hypothesis. Seven additional variants were associated (P <0.05) with ER-positive disease. CONCLUSIONS: Although no variants reached genome-wide significance (P <5 x 10(-8)), these results suggest that there is some evidence of association between candidate common germline variants and breast cancer prognosis. Larger studies from multinational collaborations are necessary to increase the power to detect associations, between common variants and prognosis, at more stringent significance levels.

The SNP rs6500843 in 16p13.3 is associated with survival specifically among chemotherapy-treated breast cancer patients.

We have utilized a two-stage study design to search for SNPs associated with the survival of breast cancer patients treated with adjuvant chemotherapy. Our initial GWS data set consisted of 805 Finnish breast cancer cases (360 treated with adjuvant chemotherapy). The top 39 SNPs from this stage were analyzed in three independent data sets: iCOGS (n=6720 chemotherapy-treated cases), SUCCESS-A (n=3596), and POSH (n=518). Two SNPs were successfully validated: rs6500843 (any chemotherapy; per-allele HR 1.16, 95% C.I. 1.08-1.26, p=0.0001, p(adjusted)=0.0091), and rs11155012 (anthracycline therapy; per-allele HR 1.21, 95% C.I. 1.08-1.35, p=0.0010, p(adjusted)=0.0270). The SNP rs6500843 was found to specifically interact with adjuvant chemotherapy, independently of standard prognostic markers (p(interaction)=0.0009), with the rs6500843-GG genotype corresponding to the highest hazard among chemotherapy-treated cases (HR 1.47, 95% C.I. 1.20-1.80). Upon trans-eQTL analysis of public microarray data, the rs6500843 locus was found to associate with the expression of a group of genes involved in cell cycle control, notably AURKA, the expression of which also exhibited differential prognostic value between chemotherapy-treated and untreated cases in our analysis of microarray data. Based on previously published information, we propose that the eQTL genes may be connected to the rs6500843 locus via a RBFOX1-FOXM1 -mediated regulatory pathway.

Identification of novel genetic markers of breast cancer survival.

BACKGROUND: Survival after a diagnosis of breast cancer varies considerably between patients, and some of this variation may be because of germline genetic variation. We aimed to identify genetic markers associated with breast cancer-specific survival. METHODS: We conducted a large meta-analysis of studies in populations of European ancestry, including 37954 patients with 2900 deaths from breast cancer. Each study had been genotyped for between 200000 and 900000 single nucleotide polymorphisms (SNPs) across the genome; genotypes for nine million common variants were imputed using a common reference panel from the 1000 Genomes Project. We also carried out subtype-specific analyses based on 6881 estrogen receptor (ER)-negative patients (920 events) and 23059 ER-positive patients (1333 events). All statistical tests were two-sided. RESULTS: We identified one new locus (rs2059614 at 11q24.2) associated with survival in ER-negative breast cancer cases (hazard ratio [HR] = 1.95, 95% confidence interval [CI] = 1.55 to 2.47, P = 1.91 x 10(-8)). Genotyping a subset of 2113 case patients, of which 300 were ER negative, provided supporting evidence for the quality of the imputation. The association in this set of case patients was stronger for the observed genotypes than for the imputed genotypes. A second locus (rs148760487 at 2q24.2) was associated at genome-wide statistical significance in initial analyses; the association was similar in ER-positive and ER-negative case patients. Here the results of genotyping suggested that the finding was less robust. CONCLUSIONS: This is currently the largest study investigating genetic variation associated with breast cancer survival. Our results have potential clinical implications, as they confirm that germline genotype can provide prognostic information in addition to standard tumor prognostic factors.

INPP4B and RAD50 have an interactive effect on survival after breast cancer.

Genes sharing similar genomic landscape have the potential to interactively orchestrate certain clinicopathological features of a disease. Deletion of the RAD50 gene is a common event particularly in basal-like breast cancer, and often occurs together with deletions of BRCA1, RB1, TP53, PTEN, and INPP4B. In this study, we investigate whether these co-deleted genes have interactive effects on survival in breast cancer. Using publicly available TCGA data, we employed Cox's proportional hazards models to test whether genomic deletions of these genes, or reduced protein or transcript levels associate with breast cancer patient survival in an interactive manner. Further validation was obtained at the transcriptional level by including 1,596 additional cases from 13 publicly available gene expression data sets from the KM-plotter database. Our results indicate that RAD50 and INPP4B associate interactively with breast cancer survival at the transcriptional, translational, and genomic levels in the TCGA data set (p (interaction) < 0.05). While neither of the genes was independently prognostic on its own, low INPP4B levels in combination with above median RAD50 abundance associated with increased hazard, both at the mRNA (HR 2.39, 95 % CI 1.20-4.76) and protein (HR 2.92, 95 % CI 1.42-6.00) levels, whereas concomitant deletion or low expression of both genes associated with unexpectedly improved survival. A similar pattern was observed in the KM-plotter data set (p (interaction) = 0.0067). We find that RAD50 and INPP4B expression levels have a synergistic influence on breast cancer survival, possibly through their effects on treatment response.

Performance of automated scoring of ER, PR, HER2, CK5/6 and EGFR in breast cancer tissue microarrays in the Breast Cancer Association Consortium.

Breast cancer risk factors and clinical outcomes vary by tumour marker expression. However, individual studies often lack the power required to assess these relationships, and large-scale analyses are limited by the need for high throughput, standardized scoring methods. To address these limitations, we assessed whether automated image analysis of immunohistochemically stained tissue microarrays can permit rapid, standardized scoring of tumour markers from multiple studies. Tissue microarray sections prepared in nine studies containing 20 263 cores from 8267 breast cancers stained for two nuclear (oestrogen receptor, progesterone receptor), two membranous (human epidermal growth factor receptor 2 and epidermal growth factor receptor) and one cytoplasmic (cytokeratin 5/6) marker were scanned as digital images. Automated algorithms were used to score markers in tumour cells using the Ariol system. We compared automated scores against visual reads, and their associations with breast cancer survival. Approximately 65-70% of tissue microarray cores were satisfactory for scoring. Among satisfactory cores, agreement between dichotomous automated and visual scores was highest for oestrogen receptor (Kappa = 0.76), followed by human epidermal growth factor receptor 2 (Kappa = 0.69) and progesterone receptor (Kappa = 0.67). Automated quantitative scores for these markers were associated with hazard ratios for breast cancer mortality in a dose-response manner. Considering visual scores of epidermal growth factor receptor or cytokeratin 5/6 as the reference, automated scoring achieved excellent negative predictive value (96-98%), but yielded many false positives (positive predictive value = 30-32%). For all markers, we observed substantial heterogeneity in automated scoring performance across tissue microarrays. Automated analysis is a potentially useful tool for large-scale, quantitative scoring of immunohistochemically stained tissue microarrays available in consortia. However, continued optimization, rigorous marker-specific quality control measures and standardization of tissue microarray designs, staining and scoring protocols is needed to enhance results.