Breast cancer genomes from CHEK2 c.1100delC mutation carriers lack somatic TP53 mutations and display a unique structural variant size distribution profile.

BACKGROUND: CHEK2 c.1100delC was the first moderate-risk breast cancer (BC) susceptibility allele discovered. Despite several genomic, transcriptomic and functional studies, however, it is still unclear how exactly CHEK2 c.1100delC promotes tumorigenesis. Since the mutational landscape of a tumor reflects the processes that have operated on its development, the aim of this study was to uncover the somatic genomic landscape of CHEK2-associated BC. METHODS: We sequenced primary BC (pBC) and normal genomes of 20 CHEK2 c.1100delC mutation carriers as well as their pBC transcriptomes. Including pre-existing cohorts, we exhaustively compared CHEK2 pBC genomes to those from BRCA1/2 mutation carriers, those that displayed homologous recombination deficiency (HRD) and ER- and ER+ pBCs, totaling to 574 pBC genomes. Findings were validated in 517 metastatic BC genomes subdivided into the same subgroups. Transcriptome data from 168 ER+ pBCs were used to derive a TP53-mutant gene expression signature and perform cluster analysis with CHEK2 BC transcriptomes. Finally, clinical outcome of CHEK2 c.1100delC carriers was compared with BC patients displaying somatic TP53 mutations in two well-described retrospective cohorts totaling to 942 independent pBC cases. RESULTS: BC genomes from CHEK2 mutation carriers were most similar to ER+ BC genomes and least similar to those of BRCA1/2 mutation carriers in terms of tumor mutational burden as well as mutational signatures. Moreover, CHEK2 BC genomes did not show any evidence of HRD. Somatic TP53 mutation frequency and the size distribution of structural variants (SVs), however, were different compared to ER+ BC. Interestingly, BC genomes with bi-allelic CHEK2 inactivation lacked somatic TP53 mutations and transcriptomic analysis indicated a shared biology with TP53 mutant BC. Moreover, CHEK2 BC genomes had an increased frequency of > 1 Mb deletions, inversions and tandem duplications with peaks at specific sizes. The high chromothripsis frequency among CHEK2 BC genomes appeared, however, not associated with this unique SV size distribution profile. CONCLUSIONS: CHEK2 BC genomes are most similar to ER+ BC genomes, but display unique features that may further unravel CHEK2-driven tumorigenesis. Increased insight into this mechanism could explain the shorter survival of CHEK2 mutation carriers that is likely driven by intrinsic tumor aggressiveness rather than endocrine resistance.

The circular RNome of primary breast cancer.

Circular RNAs (circRNAs) are a class of RNAs that is under increasing scrutiny, although their functional roles are debated. We analyzed RNA-seq data of 348 primary breast cancers and developed a method to identify circRNAs that does not rely on unmapped reads or known splice junctions. We identified 95,843 circRNAs, of which 20,441 were found recurrently. Of the circRNAs that match exon boundaries of the same gene, 668 showed a poor or even negative (R < 0.2) correlation with the expression level of the linear gene. In silico analysis showed only a minority (8.5%) of circRNAs could be explained by known splicing events. Both these observations suggest that specific regulatory processes for circRNAs exist. We confirmed the presence of circRNAs of CNOT2, CREBBP, and RERE in an independent pool of primary breast cancers. We identified circRNA profiles associated with subgroups of breast cancers and with biological and clinical features, such as amount of tumor lymphocytic infiltrate and proliferation index. siRNA-mediated knockdown of circCNOT2 was shown to significantly reduce viability of the breast cancer cell lines MCF-7 and BT-474, further underlining the biological relevance of circRNAs. Furthermore, we found that circular, and not linear, CNOT2 levels are predictive for progression-free survival time to aromatase inhibitor (AI) therapy in advanced breast cancer patients, and found that circCNOT2 is detectable in cell-free RNA from plasma. We showed that circRNAs are abundantly present, show characteristics of being specifically regulated, are associated with clinical and biological properties, and thus are relevant in breast cancer.

Functional Homologous Recombination (HR) Screening Shows the Majority of BRCA1/2-Mutant Breast and Ovarian Cancer Cell Lines Are HR-Proficient.

Tumors with a pathogenic BRCA1/2 mutation are homologous recombination (HR)-deficient (HRD) and consequently sensitive to platinum-based chemotherapy and Poly-[ADP-Ribose]-Polymerase inhibitors (PARPi). We hypothesized that functional HR status better reflects real-time HR status than BRCA1/2 mutation status. Therefore, we determined the functional HR status of 53 breast cancer (BC) and 38 ovarian cancer (OC) cell lines by measuring the formation of RAD51 foci after irradiation. Discrepancies between functional HR and BRCA1/2 mutation status were investigated using exome sequencing, methylation and gene expression data from 50 HR-related genes. A pathogenic BRCA1/2 mutation was found in 10/53 (18.9%) of BC and 7/38 (18.4%) of OC cell lines. Among BRCA1/2-mutant cell lines, 14/17 (82.4%) were HR-proficient (HRP), while 1/74 (1.4%) wild-type cell lines was HRD. For most (80%) cell lines, we explained the discrepancy between functional HR and BRCA1/2 mutation status. Importantly, 12/14 (85.7%) BRCA1/2-mutant HRP cell lines were explained by mechanisms directly acting on BRCA1/2. Finally, functional HR status was strongly associated with COSMIC single base substitution signature 3, but not BRCA1/2 mutation status. Thus, the majority of BRCA1/2-mutant cell lines do not represent a suitable model for HRD. Moreover, exclusively determining BRCA1/2 mutation status may not suffice for platinum-based chemotherapy or PARPi patient selection.

miRNA expression profiling of 51 human breast cancer cell lines reveals subtype and driver mutation-specific miRNAs.

INTRODUCTION: Breast cancer is a genetically and phenotypically complex disease. To understand the role of miRNAs in this molecular complexity, we performed miRNA expression analysis in a cohort of molecularly well-characterized human breast cancer cell lines to identify miRNAs associated with the most common molecular subtypes and the most frequent genetic aberrations. METHODS: Using a microarray carrying LNA™ modified oligonucleotide capture probes), expression levels of 725 human miRNAs were measured in 51 breast cancer cell lines. Differential miRNA expression was explored by unsupervised cluster analysis and was then associated with the molecular subtypes and genetic aberrations commonly present in breast cancer. RESULTS: Unsupervised cluster analysis using the most variably expressed miRNAs divided the 51 breast cancer cell lines into a major and a minor cluster predominantly mirroring the luminal and basal intrinsic subdivision of breast cancer cell lines. One hundred and thirteen miRNAs were differentially expressed between these two main clusters. Forty miRNAs were differentially expressed between basal-like and normal-like/claudin-low cell lines. Within the luminal-group, 39 miRNAs were associated with ERBB2 overexpression and 24 with E-cadherin gene mutations, which are frequent in this subtype of breast cancer cell lines. In contrast, 31 miRNAs were associated with E-cadherin promoter hypermethylation, which, contrary to E-cadherin mutation, is exclusively observed in breast cancer cell lines that are not of luminal origin. Thirty miRNAs were associated with p16INK4 status while only a few miRNAs were associated with BRCA1, PIK3CA/PTEN and TP53 mutation status. Twelve miRNAs were associated with DNA copy number variation of the respective locus. CONCLUSION: Luminal-basal and epithelial-mesenchymal associated miRNAs determine the subdivision of miRNA transcriptome of breast cancer cell lines. Specific sets of miRNAs were associated with ERBB2 overexpression, p16INK4a or E-cadherin mutation or E-cadherin methylation status, which implies that these miRNAs may contribute to the driver role of these genetic aberrations. Additionally, miRNAs, which are located in a genomic region showing recurrent genetic aberrations, may themselves play a driver role in breast carcinogenesis or contribute to a driver gene in their vicinity. In short, our study provides detailed molecular miRNA portraits of breast cancer cell lines, which can be exploited for functional studies of clinically important miRNAs.

Expression and Localization of Ferritin-Heavy Chain Predicts Recurrence for Breast Cancer Patients with a BRCA1/2 Mutation.

The ferritin-heavy chain (FTH1) is the catalytic subunit of the ferroxidase ferritin, which prevents oxidative DNA damage via intracellular iron storage. FTH1 was shown to be a prognostic marker for triple-negative breast cancer (BC) patients and associated with an enrichment of CD8+ effector T cells. However, whether the expression and localization of FTH1 are also associated with clinical outcome in other BC subtypes is unknown. Here, we investigated the association of FTH1 with time to survival in BCs from 222 BRCA1/2 mutation carriers by immunohistochemistry on tissue microarrays. In addition, for 51 of these patients, the association between FTH1 and specific subsets of T cells was evaluated on whole slides using automatic scoring algorithms. We revealed that nuclear FTH1 (nFTH1) expression, in multivariable analyses, was associated with a shorter disease-free (HR = 2.71, 95% CI = 1.49-4.92, p = 0.001) and metastasis-free survival (HR = 3.54, 95% CI = 1.45-8.66, p = 0.006) in patients carrying a BRCA1/2 mutation. However, we found no relation between cytoplasmic FTH1 expression and survival of BRCA1/2 mutation carriers. Moreover, we did not detect an association between FTH1 expression and the amount of CD45+ (p = 0.13), CD8+ (p = 0.18), CD4+ (p = 0.20) or FOXP3+ cells (p = 0.17). Consequently, the mechanism underlying the worse recurrence-free survival of nFTH1 expression in BRCA1/2 mutation carriers needs further investigation.

A pipeline for copy number profiling of single circulating tumour cells to assess intrapatient tumour heterogeneity.

Intrapatient tumour heterogeneity is likely a major determinant of clinical outcome in cancer patients. To assess heterogeneity in a minimally invasive manner, methods to perform single circulating tumour cell (CTC) genomics at high resolution are necessary. However, due to the rarity of CTCs, development of such methods is challenging. Here, we developed a modular single CTC analysis pipeline to assess intrapatient heterogeneity by copy number (CN) profiling. To optimize this pipeline, spike-in experiments using MCF-7 breast cancer cells were performed. The VyCAP puncher system was used to isolate single cells. The quality of whole genome amplification (WGA) products generated by REPLI-g and Ampli1™ methods, as well as the results from the Illumina Truseq and the Ampli1™ LowPass library preparation techniques, was compared. Moreover, a bioinformatic pipeline was designed to generate CN profiles from single CTCs. The optimal combination of Ampli1™ WGA and Illumina Truseq library preparation was successfully validated on patient-derived CTCs. In conclusion, we developed a novel modular pipeline to isolate single CTCs and subsequently generate detailed patient-derived CN profiles that allow assessment of intrapatient heterogeneity in future studies.

Metabotropic glutamate receptor 1 is associated with unfavorable prognosis in ER-negative and triple-negative breast cancer.

New therapies are an urgent medical need in all breast cancer subgroups. Metabotropic glutamate receptor 1 (mGluR1) is suggested as a potential new molecular target. We examined the prevalence mGluR1 expression in different clinically relevant breast cancer subgroups and determined its association with prognosis. In this retrospective cohort, 394 consecutive primary breast cancer tissues were incorporated into a tissue microarray and immunohistochemically stained for mGluR1. The prevalence of mGluR1 protein expression in different breast cancer subgroups was evaluated and correlated with metastasis-free survival (MFS) and overall survival (OS). In total, 56% (n = 219) breast cancer tissues had mGluR1 expression. In estrogen receptor (ER)-negative tumors, 31% (n = 18/58) had mGluR1 expression that was significantly associated with MFS (HR 5.00, 95% CI 1.03-24.35, p = 0.046) in multivariate analysis, independently from other prognostic factors. Of the 44 triple-negative breast cancer (TNBC), 25% (n = 11) expressed mGluR1. mGluR1 expression in TNBC was significantly associated with shorter MFS (HR 8.60, 95% CI 1.06-20.39, p = 0.044) and with poor OS (HR 16.07, 95% CI 1.16-223.10, p = 0.039). In conclusion, mGluR1 is frequently expressed in breast cancer. In ER-negative breast cancer and in TNBC mGluR1 protein expression is an unfavorable prognostic marker. This study provides rationale to explore mGluR1 as a novel target for breast cancer treatment, especially for the more aggressive TNBC.

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.

MicroRNAs as possible indicators of drug sensitivity in breast cancer cell lines.

MicroRNAs (miRNAs) regulate gene expression post-transcriptionally. In this way they might influence whether a cell is sensitive or resistant to a certain drug. So far, only a limited number of relatively small scale studies comprising few cell lines and/or drugs have been performed. To obtain a broader view on miRNAs and their association with drug response, we investigated the expression levels of 411 miRNAs in relation to drug sensitivity in 36 breast cancer cell lines. For this purpose IC50 values of a drug screen involving 34 drugs were associated with miRNA expression data of the same breast cancer cell lines. Since molecular subtype of the breast cancer cell lines is considered a confounding factor in drug association studies, multivariate analysis taking subtype into account was performed on significant miRNA-drug associations which retained 13 associations. These associations consisted of 11 different miRNAs and eight different drugs (among which Paclitaxel, Docetaxel and Veliparib). The taxanes, Paclitaxel and Docetaxel, were the only drugs having miRNAs in common: hsa-miR-187-5p and hsa-miR-106a-3p indicative of drug resistance while Paclitaxel sensitivity alone associated with hsa-miR-556-5p. Tivantinib was associated with hsa-let-7d-5p and hsa-miR-18a-5p for sensitivity and hsa-miR-637 for resistance. Drug sensitivity was associated with hsa-let-7a-5p for Bortezomib, hsa-miR-135a-3p for JNJ-707 and hsa-miR-185-3p for Panobinostat. Drug resistance was associated with hsa-miR-182-5p for Veliparib and hsa-miR-629-5p for Tipifarnib. Pathway analysis for significant miRNAs was performed to reveal biological roles, aiding to find a potential mechanistic link for the observed associations with drug response. By doing so hsa-miR-187-5p was linked to the cell cycle G2-M checkpoint in line with this checkpoint being the target of taxanes. In conclusion, our study shows that miRNAs could potentially serve as biomarkers for intrinsic drug resistance and that pathway analyses can provide additional information in this context.

GATA3 mRNA expression, but not mutation, associates with longer progression-free survival in ER-positive breast cancer patients treated with first-line tamoxifen for recurrent disease.

In breast cancer, GATA3 mutations have been associated with a favorable prognosis and the response to neoadjuvant aromatase inhibitor treatment. Therefore, we investigated whether GATA3 mutations predict the outcome of tamoxifen treatment in the advanced setting. In a retrospective study consisting of 235 hormone-naive patients with ER-positive breast cancer who received tamoxifen as first-line treatment for recurrent disease, GATA3 mutations (in 14.0% of patients) did not significantly associate with either the overall response rate (ORR) or with the length of progression-free survival (PFS) after the start of tamoxifen therapy. Interestingly, among 148 patients for whom both mutation and mRNA expression data were available, GATA3 mutations associated with an increased expression of GATA3. However, only 23.7% of GATA3 high tumors had a mutation. Evaluation of the clinical significance of GATA3 mRNA revealed that it was associated with prolonged PFS, but not with the ORR, also in multivariate analysis. Thus, GATA3 mRNA expression, but not GATA3 mutation, is an independent predictor of prolonged PFS in ER-positive breast cancer patients who received first-line tamoxifen for recurrent disease. Besides GATA3 mutation, other mechanisms must exist that underlie increased GATA3 levels.

Recurrent HOXB13 mutations in the Dutch population do not associate with increased breast cancer risk.

The HOXB13 p.G84E mutation has been firmly established as a prostate cancer susceptibility allele. Although HOXB13 also plays a role in breast tumor progression, the association of HOXB13 p.G84E with breast cancer risk is less evident. Therefore, we comprehensively interrogated the entire HOXB13 coding sequence for mutations in 1,250 non-BRCA1/2 familial breast cancer cases and 800 controls. We identified two predicted deleterious missense mutations, p.G84E and p.R217C, that were recurrent among breast cancer cases and further evaluated their association with breast cancer risk in a larger study. Taken together, 4,520 familial non-BRCA1/2 breast cancer cases and 3,127 controls were genotyped including the cases and controls of the whole gene screen. The concordance rate for the genotyping assays compared with Sanger sequencing was 100%. The prostate cancer risk allele p.G84E was identified in 18 (0.56%) of 3,187 cases and 16 (0.70%) of 2,300 controls (OR = 0.81, 95% CI = 0.41-1.59, P = 0.54). Additionally, p.R217C was identified in 10 (0.31%) of 3,208 cases and 2 (0.087%) of 2,288 controls (OR = 3.57, 95% CI = 0.76-33.57, P = 0.14). These results imply that none of the recurrent HOXB13 mutations in the Dutch population are associated with breast cancer risk, although it may be worthwhile to evaluate p.R217C in a larger study.

Breast cancer genome and transcriptome integration implicates specific mutational signatures with immune cell infiltration.

A recent comprehensive whole genome analysis of a large breast cancer cohort was used to link known and novel drivers and substitution signatures to the transcriptome of 266 cases. Here, we validate that subtype-specific aberrations show concordant expression changes for, for example, TP53, PIK3CA, PTEN, CCND1 and CDH1. We find that CCND3 expression levels do not correlate with amplification, while increased GATA3 expression in mutant GATA3 cancers suggests GATA3 is an oncogene. In luminal cases the total number of substitutions, irrespective of type, associates with cell cycle gene expression and adverse outcome, whereas the number of mutations of signatures 3 and 13 associates with immune-response specific gene expression, increased numbers of tumour-infiltrating lymphocytes and better outcome. Thus, while earlier reports imply that the sheer number of somatic aberrations could trigger an immune-response, our data suggests that substitutions of a particular type are more effective in doing so than others.

Understanding drugs in breast cancer through drug sensitivity screening.

With substantial numbers of breast tumors showing or acquiring treatment resistance, it is of utmost importance to develop new agents for the treatment of the disease, to know their effectiveness against breast cancer and to understand their relationships with other drugs to best assign the right drug to the right patient. To achieve this goal drug screenings on breast cancer cell lines are a promising approach. In this study a large-scale drug screening of 37 compounds was performed on a panel of 42 breast cancer cell lines representing the main breast cancer subtypes. Clustering, correlation and pathway analyses were used for data analysis. We found that compounds with a related mechanism of action had correlated IC50 values and thus grouped together when the cell lines were hierarchically clustered based on IC50 values. In total we found six clusters of drugs of which five consisted of drugs with related mode of action and one cluster with two drugs not previously connected. In total, 25 correlated and four anti-correlated drug sensitivities were revealed of which only one drug, Sirolimus, showed significantly lower IC50 values in the luminal/ERBB2 breast cancer subtype. We found expected interactions but also discovered new relationships between drugs which might have implications for cancer treatment regimens.

CD49f-based selection of circulating tumor cells (CTCs) improves detection across breast cancer subtypes.

Circulating tumor cells (CTCs) can be enumerated using CellSearch, but not all breast cancer subtypes, specifically those with epithelial-mesenchymal transition (EMT) characteristics, sufficiently express the enrichment (EpCAM) and selection (CK8/18/19) markers used in this method. While CD146 can detect EpCAM-negative CTCs, we here evaluated the value of various cytokeratins and CD49f to detect CK8/18/19-negative CTCs. The tested cytokeratins provided no substantial benefit, but adding CD49f to CK8/18/19 as a selection marker resulted in improved recovery of normal-like cell lines. Combined staining of CK8/18/19 and CD49f after CD146/EpCAM enrichment is likely to further improve CTC detection in breast cancer.