Obesity-associated changes in molecular biology of primary breast cancer.

Obesity is associated with an increased risk of developing breast cancer (BC) and worse prognosis in BC patients, yet its impact on BC biology remains understudied in humans. This study investigates how the biology of untreated primary BC differs according to patients' body mass index (BMI) using data from >2,000 patients. We identify several genomic alterations that are differentially prevalent in overweight or obese patients compared to lean patients. We report evidence supporting an ageing accelerating effect of obesity at the genetic level. We show that BMI-associated differences in bulk transcriptomic profile are subtle, while single cell profiling allows detection of more pronounced changes in different cell compartments. These analyses further reveal an elevated and unresolved inflammation of the BC tumor microenvironment associated with obesity, with distinct characteristics contingent on the estrogen receptor status. Collectively, our analyses imply that obesity is associated with an inflammaging-like phenotype. We conclude that patient adiposity may play a significant role in the heterogeneity of BC and should be considered for BC treatment tailoring.

Cis-regulatory mutations associate with transcriptional and post-transcriptional deregulation of gene regulatory programs in cancers.

Most cancer alterations occur in the noncoding portion of the human genome, where regulatory regions control gene expression. The discovery of noncoding mutations altering the cells' regulatory programs has been limited to few examples with high recurrence or high functional impact. Here, we show that transcription factor binding sites (TFBSs) have similar mutation loads to those in protein-coding exons. By combining cancer somatic mutations in TFBSs and expression data for protein-coding and miRNA genes, we evaluate the combined effects of transcriptional and post-transcriptional alterations on the regulatory programs in cancers. The analysis of seven TCGA cohorts culminates with the identification of protein-coding and miRNA genes linked to mutations at TFBSs that are associated with a cascading trans-effect deregulation on the cells' regulatory programs. Our analyses of cis-regulatory mutations associated with miRNAs recurrently predict 12 mature miRNAs (derived from 7 precursors) associated with the deregulation of their target gene networks. The predictions are enriched for cancer-associated protein-coding and miRNA genes and highlight cis-regulatory mutations associated with the dysregulation of key pathways associated with carcinogenesis. By combining transcriptional and post-transcriptional regulation of gene expression, our method predicts cis-regulatory mutations related to the dysregulation of key gene regulatory networks in cancer patients.

DNA copy number motifs are strong and independent predictors of survival in breast cancer.

Somatic copy number alterations are a frequent sign of genome instability in cancer. A precise characterization of the genome architecture would reveal underlying instability mechanisms and provide an instrument for outcome prediction and treatment guidance. Here we show that the local spatial behavior of copy number profiles conveys important information about this architecture. Six filters were defined to characterize regional traits in copy number profiles, and the resulting Copy Aberration Regional Mapping Analysis (CARMA) algorithm was applied to tumors in four breast cancer cohorts (n = 2919). The derived motifs represent a layer of information that complements established molecular classifications of breast cancer. A score reflecting presence or absence of motifs provided a highly significant independent prognostic predictor. Results were consistent between cohorts. The nonsite-specific occurrence of the detected patterns suggests that CARMA captures underlying replication and repair defects and could have a future potential in treatment stratification.

Author Correction: Landscape of somatic mutations in 560 breast cancer whole-genome sequences.

In the Methods section of this Article, 'greater than' should have been 'less than' in the sentence 'Putative regions of clustered rearrangements were identified as having an average inter-rearrangement distance that was at least 10 times greater than the whole-genome average for the individual sample. '. The Article has not been corrected.

The Longitudinal Transcriptional Response to Neoadjuvant Chemotherapy with and without Bevacizumab in Breast Cancer.

Purpose: Chemotherapy-induced alterations to gene expression are due to transcriptional reprogramming of tumor cells or subclonal adaptations to treatment. The effect on whole-transcriptome mRNA expression was investigated in a randomized phase II clinical trial to assess the effect of neoadjuvant chemotherapy with the addition of bevacizumab.Experimental Design: Tumor biopsies and whole-transcriptome mRNA profiles were obtained at three fixed time points with 66 patients in each arm. Altogether, 358 specimens from 132 patients were available, representing the transcriptional state before treatment start, at 12 weeks and after treatment (25 weeks). Pathologic complete response (pCR) in breast and axillary nodes was the primary endpoint.Results: pCR was observed in 15 patients (23%) receiving bevacizumab and chemotherapy and 8 patients (12%) receiving only chemotherapy. In the estrogen receptor-positive patients, 11 of 54 (20%) treated with bevacizumab and chemotherapy achieved pCR, while only 3 of 57 (5%) treated with chemotherapy reached pCR. In patients with estrogen receptor-positive tumors treated with combination therapy, an elevated immune activity was associated with good response. Proliferation was reduced after treatment in both treatment arms and most pronounced in the combination therapy arm, where the reduction in proliferation accelerated during treatment. Transcriptional alterations during therapy were subtype specific, and the effect of adding bevacizumab was most evident for luminal-B tumors.Conclusions: Clinical response and gene expression response differed between patients receiving combination therapy and chemotherapy alone. The results may guide identification of patients likely to benefit from antiangiogenic therapy. Clin Cancer Res; 23(16); 4662-70. ©2017 AACR.

Recommended Guidelines for Validation, Quality Control, and Reporting of TP53 Variants in Clinical Practice.

Accurate assessment of TP53 gene status in sporadic tumors and in the germline of individuals at high risk of cancer due to Li-Fraumeni Syndrome (LFS) has important clinical implications for diagnosis, surveillance, and therapy. Genomic data from more than 20,000 cancer genomes provide a wealth of information on cancer gene alterations and have confirmed TP53 as the most commonly mutated gene in human cancer. Analysis of a database of 70,000 TP53 variants reveals that the two newly discovered exons of the gene, exons 9ß and 9γ, generated by alternative splicing, are the targets of inactivating mutation events in breast, liver, and head and neck tumors. Furthermore, germline rearrange-ments in intron 1 of TP53 are associated with LFS and are frequently observed in sporadic osteosarcoma. In this context of constantly growing genomic data, we discuss how screening strategies must be improved when assessing TP53 status in clinical samples. Finally, we discuss how TP53 alterations should be described by using accurate nomenclature to avoid confusion in scientific and clinical reports. Cancer Res; 77(6); 1250-60. ©2017 AACR.

Assessment of TP53 Polymorphisms and MDM2 SNP309 in Premenopausal Breast Cancer Risk.

Germline polymorphic variants in cancer predisposition genes such as TP53 have been shown to impact the risk of premenopausal cancer. Accordingly, the aim of this study was to assess the spectrum of polymorphisms in TP53 and its negative regulatory gene, MDM2 (SNP309:T>G) in patients with premenopausal breast cancer. Our findings in a cohort of 40 female patients demonstrate no significant correlation between the studied polymorphisms and risk of premenopausal breast cancer. Although one polymorphism is found in high frequency in this cohort (rs1800372:A>G, 9.0%), it was not associated with the risk of developing cancer before the age of 35 years in an extended cohort of 1,420 breast cancer cases. Functional studies of the rs1800372:A>G polymorphic allele reveal that it does not affect p53 transactivation function. Further study of variants or mutations in other cancer susceptibility genes is warranted to refine our understanding of the germline contribution to premenopausal breast cancer susceptibility.

TP53 Mutations in Breast and Ovarian Cancer.

Breast and ovarian cancers are the second and fifth leading causes of cancer deaths among women. Both breast and ovarian cancers are highly heterogeneous and are presented with diverse morphology, natural history, and response to therapy. In recent years, international efforts have led to extensive molecular characterization of both breast and ovarian tumors and identified biologically and clinically relevant subtypes of the diseases based on these molecular features. The role of TP53 in tumor initiation and progression is context dependent, and abrogation of the TP53 pathway seems to be essential for the development of basal-like breast cancers and high-grade serous ovarian cancers. These subtypes of breast and ovarian cancer show several genomic similarities including high frequency of TP53 mutation, which seems to be an early, initiating, and driving alteration in these cancer subtypes.

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.

Direct Transcriptional Consequences of Somatic Mutation in Breast Cancer.

Disordered transcriptomes of cancer encompass direct effects of somatic mutation on transcription, coordinated secondary pathway alterations, and increased transcriptional noise. To catalog the rules governing how somatic mutation exerts direct transcriptional effects, we developed an exhaustive pipeline for analyzing RNA sequencing data, which we integrated with whole genomes from 23 breast cancers. Using X-inactivation analyses, we found that cancer cells are more transcriptionally active than intermixed stromal cells. This is especially true in estrogen receptor (ER)-negative tumors. Overall, 59% of substitutions were expressed. Nonsense mutations showed lower expression levels than expected, with patterns characteristic of nonsense-mediated decay. 14% of 4,234 rearrangements caused transcriptional abnormalities, including exon skips, exon reusage, fusions, and premature polyadenylation. We found productive, stable transcription from sense-to-antisense gene fusions and gene-to-intergenic rearrangements, suggesting that these mutation classes drive more transcriptional disruption than previously suspected. Systematic integration of transcriptome with genome data reveals the rules by which transcriptional machinery interprets somatic mutation.

Landscape of somatic mutations in 560 breast cancer whole-genome sequences.

We analysed whole-genome sequences of 560 breast cancers to advance understanding of the driver mutations conferring clonal advantage and the mutational processes generating somatic mutations. We found that 93 protein-coding cancer genes carried probable driver mutations. Some non-coding regions exhibited high mutation frequencies, but most have distinctive structural features probably causing elevated mutation rates and do not contain driver mutations. Mutational signature analysis was extended to genome rearrangements and revealed twelve base substitution and six rearrangement signatures. Three rearrangement signatures, characterized by tandem duplications or deletions, appear associated with defective homologous-recombination-based DNA repair: one with deficient BRCA1 function, another with deficient BRCA1 or BRCA2 function, the cause of the third is unknown. This analysis of all classes of somatic mutation across exons, introns and intergenic regions highlights the repertoire of cancer genes and mutational processes operating, and progresses towards a comprehensive account of the somatic genetic basis of breast cancer.

The Sub-Cellular Localization of WRAP53 Has Prognostic Impact in Breast Cancer.

WRAP53 protein controls intracellular trafficking of DNA repair proteins, the telomerase enzyme, and splicing factors. Functional loss of the protein has been linked to carcinogenesis, premature aging and neurodegeneration. The aim of this study was to investigate the prognostic significance of WRAP53 protein expression in breast cancer. A tissue microarray was constructed from primary breast tumors and immunostained by a polyclonal WRAP53 antibody to assess the protein expression pattern. Two different patient cohorts with long term follow-up were studied; a test- and a validation set of 154 and 668 breast tumor samples respectively. Breast cancer patients with tumor cells lacking the expression of WRAP53 in the nucleus had a significantly poorer outcome compared to patients with tumor cells expressing this protein in the nuclei (HR = 1.95, 95%CI = 1.09-3.51, p = 0.025). Nuclear localization of WRAP53 was further shown to be an independent marker of prognosis in multivariate analysis (HR = 2.57, 95%CI = 1.27-5.19, p = 0.008), and also significantly associated with better outcome in patients with TP53 mutation. Here we show that the sub-cellular localization of the WRAP53 protein has a significant impact on breast cancer survival, and thus has a potential as a clinical marker in diagnostics and treatment.

Lymphocyte Invasion in IC10/Basal-Like Breast Tumors Is Associated with Wild-Type TP53.

UNLABELLED: Lymphocytic infiltration is associated with better prognosis in several epithelial malignancies including breast cancer. The tumor suppressor TP53 is mutated in approximately 30% of breast adenocarcinomas, with varying frequency across molecular subtypes. In this study of 1,420 breast tumors, we tested for interaction between TP53 mutation status and tumor subtype determined by PAM50 and integrative cluster analysis. In integrative cluster 10 (IC10)/basal-like breast cancer, we identify an association between lymphocytic infiltration, determined by an expression score, and retention of wild-type TP53. The expression-derived score agreed with the degree of lymphocytic infiltration assessed by pathologic review, and application of the Nanodissect algorithm was suggestive of this infiltration being primarily of cytotoxic T lymphocytes (CTL). Elevated expression of this CTL signature was associated with longer survival in IC10/Basal-like tumors. These findings identify a new link between the TP53 pathway and the adaptive immune response in estrogen receptor (ER)-negative breast tumors, suggesting a connection between TP53 inactivation and failure of tumor immunosurveillance. IMPLICATIONS: The association of lymphocytic invasion of ER-negative breast tumors with the retention of wild-type TP53 implies a novel protective connection between TP53 function and tumor immunosurveillance.

TP53 mutation spectrum in breast cancer is subtype specific and has distinct prognostic relevance.

PURPOSE: In breast cancer, the TP53 gene is frequently mutated and the mutations have been associated with poor prognosis. The prognostic impact of the different types of TP53 mutations across the different molecular subtypes is still poorly understood. Here, we characterize the spectrum and prognostic significance of TP53 mutations with respect to the PAM50 subtypes and integrative clusters (IC). EXPERIMENTAL DESIGN: TP53 mutation status was obtained for 1,420 tumor samples from the METABRIC cohort by sequencing all coding exons using the Sanger method. RESULTS: TP53 mutations were found in 28.3% of the tumors, conferring a worse overall and breast cancer-specific survival [HR = 2.03; 95% confidence interval (CI), 1.65-2.48, P < 0.001], and were also found to be an independent marker of poor prognosis in estrogen receptor-positive cases (HR = 1.86; 95% CI, 1.39-2.49, P < 0.001). The mutation spectrum of TP53 varied between the breast cancer subtypes, and individual alterations showed subtype-specific association. TP53 mutations were associated with increased mortality in patients with luminal B, HER2-enriched, and normal-like tumors, but not in patients with luminal A and basal-like tumors. Similar observations were made in ICs, where mutation associated with poorer outcome in IC1, IC4, and IC5. The combined effect of TP53 mutation, TP53 LOH, and MDM2 amplification on mortality was additive. CONCLUSION: This study reveals that TP53 mutations have different clinical relevance in molecular subtypes of breast cancer, and suggests diverse roles for TP53 in the biology underlying breast cancer development.

AMG 900, pan-Aurora kinase inhibitor, preferentially inhibits the proliferation of breast cancer cell lines with dysfunctional p53.

Aurora kinases play important roles in cell division and are frequently overexpressed in human cancer. AMG 900 is a novel pan-Aurora kinase inhibitor currently being tested in Phase I clinical trials. We aimed to evaluate the in vitro activity of AMG 900 in a panel of 44 human breast cancer and immortalized cell lines and identify predictors of response. AMG 900 inhibited proliferation at low nanomolar concentrations in all cell lines tested. Response was further classified based on the induction of lethality. 25 cell lines were classified as highly sensitive (lethality at 10 nM of AMG 900 >10 %), 19 cell lines as less sensitive to AMG 900 (lethality at 10 nM of AMG 900 <10 %). Traditional molecular subtypes of breast cancer did not predict for this differential response. There was a weak association between AURKA amplification and response to AMG 900 (response ratio = 2.53, p = 0.09). mRNA expression levels of AURKA, AURKB, and AURKC and baseline protein levels of Aurora kinases A and B did not significantly associate with response. Cell lines with TP53 loss of function mutations (RR = 1.86, p = 0.004) and low baseline p21 protein levels (RR = 2.28, p = 0.0004) were far more likely to be classified as highly sensitive to AMG 900. AMG 900 induced p53 and p21 protein expression in cell lines with wt TP53. AMG 900 caused the accumulation of cells with >4 N DNA content in a majority of cell lines independently of sensitivity and p53 status. AMG 900 induced more pronounced apoptosis in highly sensitive p53-dysfunctional cell lines. We have found that AMG 900 is highly active in breast cancer cell lines and that TP53 loss of function mutations as well as low baseline expression of p21 protein predict strongly for increased sensitivity to this compound in vitro.

AZGP1 and SPDEF mRNA expression differentiates breast carcinoma from ovarian serous carcinoma.

The ANPEP, AZGP1, and SPDEF genes were previously found to be overexpressed in breast compared to ovarian carcinoma effusions. The present study validated this finding in a larger cohort consisting of both primary and metastatic tumors. ANPEP, AZGP1, and SPDEF mRNA expression was investigated in 83 breast carcinomas (57 primary carcinomas and 26 effusions) and 40 ovarian carcinomas (20 primary carcinomas and 20 effusions) using qPCR. ANPEP protein expression was immunohistochemically analyzed in 53 breast carcinoma effusions and patient-matched primary carcinomas (n = 25) and lymph node metastases (n = 16). mRNA and protein levels were studied for association with tumor type and anatomic site, and for clinical role in breast carcinoma. AZGP1 and SPDEF mRNA was overexpressed in breast compared to ovarian carcinoma (both p < 0.001). AZGP1 mRNA was overexpressed in primary breast carcinoma compared to effusions (p < 0.001), with opposite findings for ANPEP (p = 0.044). AZGP1 mRNA expression correlated with positive ER status (p = 0.032) and grade 1 histology (p = 0.011), whereas SPDEF mRNA levels were associated with positive ER (p = 0.002) and PR (p = 0.013) status and tamoxifen treatment (p = 0.004). ANPEP protein expression was higher in breast carcinoma effusions compared to primary tumors and lymph node metastases (both p = 0.001). ANPEP, AZGP1, and SPDEF levels were unrelated to disease-free or overall survival. This is the first study documenting ANPEP, AZGP1, and SPDEF expression in breast carcinoma effusions. AZGP1 and SPDEF may be novel molecular markers for the differentiation of breast from ovarian carcinoma. ANPEP may be involved in breast carcinoma progression in view of its overexpression in effusions compared to solid specimens.

The landscape of cancer genes and mutational processes in breast cancer.

All cancers carry somatic mutations in their genomes. A subset, known as driver mutations, confer clonal selective advantage on cancer cells and are causally implicated in oncogenesis, and the remainder are passenger mutations. The driver mutations and mutational processes operative in breast cancer have not yet been comprehensively explored. Here we examine the genomes of 100 tumours for somatic copy number changes and mutations in the coding exons of protein-coding genes. The number of somatic mutations varied markedly between individual tumours. We found strong correlations between mutation number, age at which cancer was diagnosed and cancer histological grade, and observed multiple mutational signatures, including one present in about ten per cent of tumours characterized by numerous mutations of cytosine at TpC dinucleotides. Driver mutations were identified in several new cancer genes including AKT2, ARID1B, CASP8, CDKN1B, MAP3K1, MAP3K13, NCOR1, SMARCD1 and TBX3. Among the 100 tumours, we found driver mutations in at least 40 cancer genes and 73 different combinations of mutated cancer genes. The results highlight the substantial genetic diversity underlying this common disease.

The life history of 21 breast cancers.

Cancer evolves dynamically as clonal expansions supersede one another driven by shifting selective pressures, mutational processes, and disrupted cancer genes. These processes mark the genome, such that a cancer's life history is encrypted in the somatic mutations present. We developed algorithms to decipher this narrative and applied them to 21 breast cancers. Mutational processes evolve across a cancer's lifespan, with many emerging late but contributing extensive genetic variation. Subclonal diversification is prominent, and most mutations are found in just a fraction of tumor cells. Every tumor has a dominant subclonal lineage, representing more than 50% of tumor cells. Minimal expansion of these subclones occurs until many hundreds to thousands of mutations have accumulated, implying the existence of long-lived, quiescent cell lineages capable of substantial proliferation upon acquisition of enabling genomic changes. Expansion of the dominant subclone to an appreciable mass may therefore represent the final rate-limiting step in a breast cancer's development, triggering diagnosis.

Mutational processes molding the genomes of 21 breast cancers.

All cancers carry somatic mutations. The patterns of mutation in cancer genomes reflect the DNA damage and repair processes to which cancer cells and their precursors have been exposed. To explore these mechanisms further, we generated catalogs of somatic mutation from 21 breast cancers and applied mathematical methods to extract mutational signatures of the underlying processes. Multiple distinct single- and double-nucleotide substitution signatures were discernible. Cancers with BRCA1 or BRCA2 mutations exhibited a characteristic combination of substitution mutation signatures and a distinctive profile of deletions. Complex relationships between somatic mutation prevalence and transcription were detected. A remarkable phenomenon of localized hypermutation, termed "kataegis," was observed. Regions of kataegis differed between cancers but usually colocalized with somatic rearrangements. Base substitutions in these regions were almost exclusively of cytosine at TpC dinucleotides. The mechanisms underlying most of these mutational signatures are unknown. However, a role for the APOBEC family of cytidine deaminases is proposed.

The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups.

The elucidation of breast cancer subgroups and their molecular drivers requires integrated views of the genome and transcriptome from representative numbers of patients. We present an integrated analysis of copy number and gene expression in a discovery and validation set of 997 and 995 primary breast tumours, respectively, with long-term clinical follow-up. Inherited variants (copy number variants and single nucleotide polymorphisms) and acquired somatic copy number aberrations (CNAs) were associated with expression in ~40% of genes, with the landscape dominated by cis- and trans-acting CNAs. By delineating expression outlier genes driven in cis by CNAs, we identified putative cancer genes, including deletions in PPP2R2A, MTAP and MAP2K4. Unsupervised analysis of paired DNA­RNA profiles revealed novel subgroups with distinct clinical outcomes, which reproduced in the validation cohort. These include a high-risk, oestrogen-receptor-positive 11q13/14 cis-acting subgroup and a favourable prognosis subgroup devoid of CNAs. Trans-acting aberration hotspots were found to modulate subgroup-specific gene networks, including a TCR deletion-mediated adaptive immune response in the 'CNA-devoid' subgroup and a basal-specific chromosome 5 deletion-associated mitotic network. Our results provide a novel molecular stratification of the breast cancer population, derived from the impact of somatic CNAs on the transcriptome.