Integrative clustering reveals a novel split in the luminal A subtype of breast cancer with impact on outcome

Background: Breast cancer is a heterogeneous disease at the clinical and molecular level. In this study we integrate classifications extracted from five different molecular levels in order to identify integrated subtypes. Methods: Tumor tissue from 425 patients with primary breast cancer from the Oslo2 study was cut and blended, and divided into fractions for DNA, RNA and protein isolation and metabolomics, allowing the acquisition of representative and comparable molecular data. Patients were stratified into groups based on their tumor characteristics from five different molecular levels, using various clustering methods. Finally, all previously identified and newly determined subgroups were combined in a multilevel classification using a "cluster-of-clusters" approach with consensus clustering. Results: Based on DNA copy number data, tumors were categorized into three groups according to the complex arm aberration index. mRNA expression profiles divided tumors into five molecular subgroups according to PAM50 subtyping, and clustering based on microRNA expression revealed four subgroups. Reverse-phase protein array data divided tumors into five subgroups. Hierarchical clustering of tumor metabolic profiles revealed three clusters. Combining DNA copy number and mRNA expression classified tumors into seven clusters based on pathway activity levels, and tumors were classified into ten subtypes using integrative clustering. The final consensus clustering that incorporated all aforementioned subtypes revealed six major groups. Five corresponded well with the mRNA subtypes, while a sixth group resulted from a split of the luminal A subtype; these tumors belonged to distinct microRNA clusters. Gain-of-function studies using MCF-7 cells showed that microRNAs differentially expressed between the luminal A clusters were important for cancer cell survival. These microRNAs were used to validate the split in luminal A tumors in four independent breast cancer cohorts. In two cohorts the microRNAs divided tumors into subgroups with significantly different outcomes, and in another a trend was observed. Conclusions: The six integrated subtypes identified confirm the heterogeneity of breast cancer and show that finer subdivisions of subtypes are evident. Increasing knowledge of the heterogeneity of the luminal A subtype may add pivotal information to guide therapeutic choices, evidently bringing us closer to improved treatment for this largest subgroup of breast cancer.Peer reviewe

The impact of coding germline variants on contralateral breast cancer risk and survival

Evidence linking coding germline variants in breast cancer (BC)-susceptibility genes other than BRCA1, BRCA2, and CHEK2 with contralateral breast cancer (CBC) risk and breast cancer-specific survival (BCSS) is scarce. The aim of this study was to assess the association of protein-truncating variants (PTVs) and rare missense variants (MSVs) in nine known (ATM, BARD1, BRCA1, BRCA2, CHEK2, PALB2, RAD51C, RAD51D, and TP53) and 25 suspected BC-susceptibility genes with CBC risk and BCSS. Hazard ratios (HRs) and 95% confidence intervals (CIs) were estimated with Cox regression models. Analyses included 34,401 women of European ancestry diagnosed with BC, including 676 CBCs and 3,449 BC deaths; the median follow-up was 10.9 years. Subtype analyses were based on estrogen receptor (ER) status of the first BC. Combined PTVs and pathogenic/likely pathogenic MSVs in BRCA1, BRCA2, and TP53 and PTVs in CHEK2 and PALB2 were associated with increased CBC risk [HRs (95% CIs): 2.88 (1.70–4.87), 2.31 (1.39–3.85), 8.29 (2.53–27.21), 2.25 (1.55–3.27), and 2.67 (1.33–5.35), respectively]. The strongest evidence of association with BCSS was for PTVs and pathogenic/likely pathogenic MSVs in BRCA2 (ER-positive BC) and TP53 and PTVs in CHEK2 [HRs (95% CIs): 1.53 (1.13–2.07), 2.08 (0.95–4.57), and 1.39 (1.13–1.72), respectively, after adjusting for tumor characteristics and treatment]. HRs were essentially unchanged when censoring for CBC, suggesting that these associations are not completely explained by increased CBC risk, tumor characteristics, or treatment. There was limited evidence of associations of PTVs and/or rare MSVs with CBC risk or BCSS for the 25 suspected BC genes. The CBC findings are relevant to treatment decisions, follow-up, and screening after BC diagnosis.</p

Association of the CHEK2 c.1100delC variant, radiotherapy, and systemic treatment with contralateral breast cancer risk and breast cancer-specific survival

Background: Breast cancer (BC) patients with a germline CHEK2 c.1100delC variant have an increased risk of contralateral BC (CBC) and worse BC-specific survival (BCSS) compared to non-carriers.Aim: To assessed the associations of CHEK2 c.1100delC, radiotherapy, and systemic treatment with CBC risk and BCSS.Methods: Analyses were based on 82,701 women diagnosed with a first primary invasive BC including 963 CHEK2 c.1100delC carriers; median follow-up was 9.1 years. Differential associations with treatment by CHEK2 c.1100delC status were tested by including interaction terms in a multivariable Cox regression model. A multi-state model was used for further insight into the relation between CHEK2 c.1100delC status, treatment, CBC risk and death. Results: There was no evidence for differential associations of therapy with CBC risk by CHEK2 c.1100delC status. The strongest association with reduced CBC risk was observed for the combination of chemotherapy and endocrine therapy [HR (95% CI): 0.66 (0.55-0.78)]. No association was observed with radiotherapy.Results from the multi-state model showed shorter BCSS for CHEK2 c.1100delC carriers versus non-carriers also after accounting for CBC occurrence [HR (95% CI): 1.30 (1.09-1.56)].Conclusion: Systemic therapy was associated with reduced CBC risk irrespective of CHEK2 c.1100delC status. Moreover, CHEK2 c.1100delC carriers had shorter BCSS, which appears not to be fully explained by their CBC risk.Peer reviewe

Rare germline copy number variants (CNVs) and breast cancer risk.

Funder: CIHRGermline copy number variants (CNVs) are pervasive in the human genome but potential disease associations with rare CNVs have not been comprehensively assessed in large datasets. We analysed rare CNVs in genes and non-coding regions for 86,788 breast cancer cases and 76,122 controls of European ancestry with genome-wide array data. Gene burden tests detected the strongest association for deletions in BRCA1 (P = 3.7E-18). Nine other genes were associated with a p-value < 0.01 including known susceptibility genes CHEK2 (P = 0.0008), ATM (P = 0.002) and BRCA2 (P = 0.008). Outside the known genes we detected associations with p-values < 0.001 for either overall or subtype-specific breast cancer at nine deletion regions and four duplication regions. Three of the deletion regions were in established common susceptibility loci. To the best of our knowledge, this is the first genome-wide analysis of rare CNVs in a large breast cancer case-control dataset. We detected associations with exonic deletions in established breast cancer susceptibility genes. We also detected suggestive associations with non-coding CNVs in known and novel loci with large effects sizes. Larger sample sizes will be required to reach robust levels of statistical significance

Germline HOXB13 mutations p.G84E and p.R217C do not confer an increased breast cancer risk

In breast cancer, high levels of homeobox protein Hox-B13 (HOXB13) have been associated with disease progression of ER-positive breast cancer patients and resistance to tamoxifen treatment. Since HOXB13 p.G84E is a prostate cancer risk allele, we evaluated the association between HOXB13 germline mutations and breast cancer risk in a previous study consisting of 3,270 familial non-BRCA1/2 breast cancer cases and 2,327 controls from the Netherlands. Although both recurrent HOXB13 mutations p.G84E and p.R217C were not associated with breast cancer risk, the risk estimation for p.R217C was not very precise. To provide more conclusive evidence regarding the role of HOXB13 in breast cancer susceptibility, we here evaluated the association between HOXB13 mutations and increased breast cancer risk within 81 studies of the international Breast Cancer Association Consortium containing 68,521 invasive breast cancer patients and 54,865 controls. Both HOXB13 p.G84E and p.R217C did not associate with the development of breast cancer in European women, neither in the overall analysis (OR = 1.035, 95% CI = 0.859-1.246, P = 0.718 and OR = 0.798, 95% CI = 0.482-1.322, P = 0.381 respectively), nor in specific high-risk subgroups or breast cancer subtypes. Thus, although involved in breast cancer progression, HOXB13 is not a material breast cancer susceptibility gene.Peer reviewe

rs2735383, located at a microRNA binding site in the 3 ' UTR of NBS1, is not associated with breast cancer risk

NBS1, also known as NBN, plays an important role in maintaining genomic stability. Interestingly, rs2735383 G > C, located in a microRNA binding site in the 3'-untranslated region (UTR) of NBS1, was shown to be associated with increased susceptibility to lung and colorectal cancer. However, the relation between rs2735383 and susceptibility to breast cancer is not yet clear. Therefore, we genotyped rs2735383 in 1,170 familial non-BRCA1/2 breast cancer cases and 1,077 controls using PCR-based restriction fragment length polymorphism (RFLP-PCR) analysis, but found no association between rs2735383CC and breast cancer risk (OR = 1.214, 95% CI = 0.936-1.574, P = 0.144). Because we could not exclude a small effect size due to a limited sample size, we further analyzed imputed rs2735383 genotypes (r(2) > 0.999) of 47,640 breast cancer cases and 46,656 controls from the Breast Cancer Association Consortium (BCAC). However, rs2735383CC was not associated with overall breast cancer risk in European (OR = 1.014, 95% CI = 0.969-1.060, P = 0.556) nor in Asian women (OR = 0.998, 95% CI = 0.905-1.100, P = 0.961). Subgroup analyses by age, age at menarche, age at menopause, menopausal status, number of pregnancies, breast feeding, family history and receptor status also did not reveal a significant association. This study therefore does not support the involvement of the genotype at NBS1 rs2735383 in breast cancer susceptibility.Peer reviewe

Germline HOXB13 mutations p.G84E and p.R217C do not confer an increased breast cancer risk

Abstract: In breast cancer, high levels of homeobox protein Hox-B13 (HOXB13) have been associated with disease progression of ER-positive breast cancer patients and resistance to tamoxifen treatment. Since HOXB13 p.G84E is a prostate cancer risk allele, we evaluated the association between HOXB13 germline mutations and breast cancer risk in a previous study consisting of 3,270 familial non-BRCA1/2 breast cancer cases and 2,327 controls from the Netherlands. Although both recurrent HOXB13 mutations p.G84E and p.R217C were not associated with breast cancer risk, the risk estimation for p.R217C was not very precise. To provide more conclusive evidence regarding the role of HOXB13 in breast cancer susceptibility, we here evaluated the association between HOXB13 mutations and increased breast cancer risk within 81 studies of the international Breast Cancer Association Consortium containing 68,521 invasive breast cancer patients and 54,865 controls. Both HOXB13 p.G84E and p.R217C did not associate with the development of breast cancer in European women, neither in the overall analysis (OR = 1.035, 95% CI = 0.859–1.246, P = 0.718 and OR = 0.798, 95% CI = 0.482–1.322, P = 0.381 respectively), nor in specific high-risk subgroups or breast cancer subtypes. Thus, although involved in breast cancer progression, HOXB13 is not a material breast cancer susceptibility gene

CYP3A7*1C allele: linking premenopausal oestrone and progesterone levels with risk of hormone receptor-positive breast cancers

Funder: Breast Cancer Now (BCN); doi: https://doi.org/10.13039/100009794Funder: Cancer Research UK (CRUK); doi: https://doi.org/10.13039/501100000289Funder: RCUK | Medical Research Council (MRC); doi: https://doi.org/10.13039/501100000265Funder: U.S. Department of Health & Human Services | National Institutes of Health (NIH)Funder: Wellcome Trust (Wellcome); doi: https://doi.org/10.13039/100004440Funder: EC | EC Seventh Framework Programm | FP7 Ideas: European Research Council (FP7-IDEAS-ERC - Specific Programme: "Ideas" Implementing the Seventh Framework Programme of the European Community for Research, Technological Development and Demonstration Activities (2007 to 2013)); doi: https://doi.org/10.13039/100011199; Grant(s): HEALTH-F2-2009-223175, HEALTH-F2-2009-223175Funder: Genome Canada (Génome Canada); doi: https://doi.org/10.13039/100008762Funder: Gouvernement du Canada | Canadian Institutes of Health Research (Instituts de Recherche en Santé du Canada); doi: https://doi.org/10.13039/501100000024Funder: Quebec Breast cancer Foundation Genome QuebecFunder: U.S. Department of Health & Human Services | NIH | U.S. National Library of Medicine (NLM); doi: https://doi.org/10.13039/100000092Funder: EC | EC Seventh Framework Programm | FP7 Ideas: European Research Council (FP7-IDEAS-ERC - Specific Programme: "Ideas" Implementing the Seventh Framework Programme of the European Community for Research, Technological Development and Demonstration Activities (2007 to 2013))Funder: European Union’s Horizon 2020Funder: Deutsche Krebshilfe (German Cancer Aid); doi: https://doi.org/10.13039/501100005972Funder: BCAST - European Union’s Horizon 2020Funder: Breast Cancer Now; doi: https://doi.org/10.13039/501100007913Abstract: Background: Epidemiological studies provide strong evidence for a role of endogenous sex hormones in the aetiology of breast cancer. The aim of this analysis was to identify genetic variants that are associated with urinary sex-hormone levels and breast cancer risk. Methods: We carried out a genome-wide association study of urinary oestrone-3-glucuronide and pregnanediol-3-glucuronide levels in 560 premenopausal women, with additional analysis of progesterone levels in 298 premenopausal women. To test for the association with breast cancer risk, we carried out follow-up genotyping in 90,916 cases and 89,893 controls from the Breast Cancer Association Consortium. All women were of European ancestry. Results: For pregnanediol-3-glucuronide, there were no genome-wide significant associations; for oestrone-3-glucuronide, we identified a single peak mapping to the CYP3A locus, annotated by rs45446698. The minor rs45446698-C allele was associated with lower oestrone-3-glucuronide (−49.2%, 95% CI −56.1% to −41.1%, P = 3.1 × 10–18); in follow-up analyses, rs45446698-C was also associated with lower progesterone (−26.7%, 95% CI −39.4% to −11.6%, P = 0.001) and reduced risk of oestrogen and progesterone receptor-positive breast cancer (OR = 0.86, 95% CI 0.82–0.91, P = 6.9 × 10–8). Conclusions: The CYP3A7*1C allele is associated with reduced risk of hormone receptor-positive breast cancer possibly mediated via an effect on the metabolism of endogenous sex hormones in premenopausal women

Tumor expression, plasma levels and genetic polymorphisms of the coagulation inhibitor TFPI are associated with clinicopathological parameters and survival in breast cancer, in contrast to the coagulation initiator TF

Introduction Hypercoagulability in malignancy increases the risk of thrombosis, but is also involved in cancer progression. Experimental studies suggest that tissue factor (TF) and tissue factor pathway inhibitor (TFPI) are involved in cancer biology as a tumor- promoter and suppressor, respectively, but the clinical significance is less clear. Here, we aimed to investigate the clinical relevance of TF and TFPI genetic and phenotypic diversity in breast cancer. Methods The relationship between tumor messenger RNA (mRNA) expression and plasma levels of TF and TFPI (α and β), tagging single nucleotide polymorphisms (tagSNPs) in F3 (TF) (n = 6) and TFPI (n = 18), and clinicopathological characteristics and molecular tumor subtypes were explored in 152 treatment naive breast cancer patients. The effect of tumor expressed TF and TFPIα and TFPIβ on survival was investigated in a merged breast cancer dataset of 1881 patients. Results Progesterone receptor negative patients had higher mRNA expression of total TFPI (α + β) (P = 0.021) and TFPIβ (P = 0.014) in tumors. TF mRNA expression was decreased in grade 3 tumors (P = 0.003). In plasma, total TFPI levels were decreased in patients with larger tumors (P = 0.013). SNP haplotypes of TFPI, but not TF, were associated with specific clinicopathological characteristics like tumor size (odds ratio (OR) 3.14, P = 0.004), triple negativity (OR 2.4, P = 0.004), lymph node spread (OR 3.34, P = 0.006), and basal-like (OR 2.3, P = 0.011) and luminal B (OR 3.5, P = 0.005) molecular tumor subtypes. Increased expression levels of TFPIα and TFPIβ in breast tumors were associated with better outcome in all tumor subtypes combined (P = 0.007 and P = 0.005) and in multiple subgroups, including lymph node positive subjects (P = 0.006 and P = 0.034). Conclusions This study indicates that genetic and phenotypic variation of both TFPIα and TFPIβ, more than TF, are markers of cancer progression. Together with the previously demonstrated tumor suppressor effects of TFPI, the beneficial effect of tumor expressed TFPI on survival, renders TFPI as a potential anticancer agent, and the clinical significance of TFPI in cancer deserves further investigation