Mutations in E3 ubiquitin ligase Parkin have been linked to familial Parkinson's disease. Accumulating evidence suggests that Parkin is a tumor suppressor, but the underlying mechanism is poorly understood. Here we show that Parkin is an E3 ubiquitin ligase for hypoxia-inducible factor 1α (HIF-1α). Parkin interacts with HIF-1α and promotes HIF-1α degradation through ubiquitination, which in turn inhibits metastasis of breast cancer cells. Parkin downregulation in breast cancer cells promotes metastasis, which can be inhibited by targeting HIF-1α with RNA interference or the small-molecule inhibitor YC-1. We further identify lysine 477 (K477) of HIF-1α as a major ubiquitination site for Parkin. K477R HIF-1α mutation and specific cancer-associated Parkin mutations largely abolish the functions of Parkin to ubiquitinate HIF-1α and inhibit cancer metastasis. Importantly, Parkin expression is inversely correlated with HIF-1α expression and metastasis in breast cancer. Our results reveal an important mechanism for Parkin in tumor suppression and HIF-1α regulation.

Despite toxic side effects and limited durable response, the current standard-of-care treatment for high grade serous ovarian cancer (HGSOC) remains platinum/taxane-based chemotherapy. Given that the overall prognosis has not improved drastically over the past several decades, there is a critical need to understand the underlying mechanisms that lead to tumour development and progression.

Dysregulation of PI3K/Akt signaling is a dominant feature in basal-like or triple-negative breast cancers (TNBC). However, the mechanisms regulating this pathway are largely unknown in this subset of aggressive tumors. Here we demonstrate that the transcription factor SOX4 is a key regulator of PI3K signaling in TNBC. Genomic and proteomic analyses coupled with mechanistic studies identified TGFBR2 as a direct transcriptional target of SOX4 and demonstrated that TGFBR2 is required to mediate SOX4-dependent PI3K signaling. We further report that SOX4 and the SWI/SNF ATPase SMARCA4, which are uniformly overexpressed in basal-like tumors, form a previously unreported complex that is required to maintain an open chromatin conformation at the TGFBR2 regulatory regions in order to mediate TGFBR2 expression and PI3K signaling. Collectively, our findings delineate the mechanism by which SOX4 and SMARCA4 cooperatively regulate PI3K/Akt signaling and suggest that this complex may play an essential role in TNBC genesis and/or progression.

Inhibition of the HER2/ERBB2 receptor is a keystone to treating HER2-positive malignancies, particularly breast cancer, but a significant fraction of HER2-positive (HER2+) breast cancers recur or fail to respond. Anti-HER2 monoclonal antibodies, like trastuzumab or pertuzumab, and ATP active site inhibitors like lapatinib, commonly lack durability because of adaptive changes in the tumor leading to resistance. HER2+ cell line responses to inhibition with lapatinib were analyzed by RNAseq and ChIPseq to characterize transcriptional and epigenetic changes. Motif analysis of lapatinib-responsive genomic regions implicated the pioneer transcription factor FOXA1 as a mediator of adaptive responses. Lapatinib in combination with FOXA1 depletion led to dysregulation of enhancers, impaired adaptive upregulation of HER3, and decreased proliferation. HER2-directed therapy using clinically relevant drugs (trastuzumab with or without lapatinib or pertuzumab) in a 7-day clinical trial designed to examine early pharmacodynamic response to antibody-based anti-HER2 therapy showed reduced FOXA1 expression was coincident with decreased HER2 and HER3 levels, decreased proliferation gene signatures, and increased immune gene signatures. This highlights the importance of the immune response to anti-HER2 antibodies and suggests that inhibiting FOXA1-mediated adaptive responses in combination with HER2 targeting is a potential therapeutic strategy.

KRAS is one of the most frequently mutated proto-oncogenes in pancreatic ductal adenocarcinoma (PDAC) and aberrantly activated in triple-negative breast cancer (TNBC). A profound role of microRNAs (miRNAs) in the pathogenesis of human cancer is being uncovered, including in cancer therapy. Using in silico prediction algorithms, we identified miR-873 as a potential regulator of KRAS, and we investigated its role in PDAC and TNBC. We found that reduced miR-873 expression is associated with shorter patient survival in both cancers. miR-873 expression is significantly repressed in PDAC and TNBC cell lines and inversely correlated with KRAS levels. We demonstrate that miR-873 directly bound to the 3' UTR of KRAS mRNA and suppressed its expression. Notably, restoring miR-873 expression induced apoptosis; recapitulated the effects of KRAS inhibition on cell proliferation, colony formation, and invasion; and suppressed the activity of ERK and PI3K/AKT, while overexpression of KRAS rescued the effects mediated by miR-873. Moreover, in vivo delivery of miR-873 nanoparticles inhibited KRAS expression and tumor growth in PDAC and TNBC tumor models. In conclusion, we provide the first evidence that miR-873 acts as a tumor suppressor by targeting KRAS and that miR-873-based gene therapy may be a therapeutic strategy in PDAC and TNBC.

Triple-negative breast cancer (TNBC), which is largely synonymous with the basal-like molecular subtype, is the 5th leading cause of cancer deaths for women in the United States. The overall prognosis for TNBC patients remains poor given that few treatment options exist; including targeted therapies (not FDA approved), and multi-agent chemotherapy as standard-of-care treatment. TNBC like other complex diseases is governed by the perturbations of the complex interaction networks thereby elucidating the underlying molecular mechanisms of this disease in the context of network principles, which have the potential to identify targets for drug development. Here, we present an integrated "omics" approach based on the use of transcriptome and interactome data to identify dynamic/active protein-protein interaction networks (PPINs) in TNBC patients. We have identified three highly connected modules, EED, DHX9, and AURKA, which are extremely activated in TNBC tumors compared to both normal tissues and other breast cancer subtypes. Based on the functional analyses, we propose that these modules are potential drivers of proliferation and, as such, should be considered candidate molecular targets for drug development or drug repositioning in TNBC. Consistent with this argument, we repurposed steroids, anti-inflammatory agents, anti-infective agents, cardiovascular agents for patients with basal-like breast cancer. Finally, we have performed essential metabolite analysis on personalized genome-scale metabolic models and found that metabolites such as sphingosine-1-phosphate and cholesterol-sulfate have utmost importance in TNBC tumor growth.

Mounting evidence supports a role for the immune system in breast cancer outcomes. The ability to distinguish highly immunogenic tumors susceptible to anti-tumor immunity from weakly immunogenic or inherently immune-resistant tumors would guide development of therapeutic strategies in breast cancer. Genomic, transcriptomic and clinical data from The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) breast cancer cohorts were used to examine statistical associations between tumor mutational burden (TMB) and the survival of patients whose tumors were assigned to previously-described prognostic immune subclasses reflecting favorable, weak or poor immune-infiltrate dispositions (FID, WID or PID, respectively). Tumor immune subclasses were associated with survival in patients with high TMB (TMB-Hi, P < 0.001) but not in those with low TMB (TMB-Lo, P = 0.44). This statistical relationship was confirmed in the METABRIC cohort (TMB-Hi, P = 0.047; TMB-Lo, P = 0.39), and also found to hold true in the more-indolent Luminal A tumor subtype (TMB-Hi, P = 0.011; TMB-Lo, P = 0.91). In TMB-Hi tumors, the FID subclass was associated with prolonged survival independent of tumor stage, molecular subtype, age and treatment. Copy number analysis revealed the reproducible, preferential amplification of chromosome 1q immune-regulatory genes in the PID immune subclass. These findings demonstrate a previously unappreciated role for TMB as a determinant of immune-mediated survival of breast cancer patients and identify candidate immune-regulatory mechanisms associated with immunologically cold tumors. Immune subtyping of breast cancers may offer opportunities for therapeutic stratification.

Chromosome 11q13-14 amplification is a defining feature of high-risk hormone receptor-positive (HR+) breast cancer; however, the mechanism(s) by which this amplicon contributes to breast tumorigenesis remains unclear. In the current study, proteogenomic analyses of >3000 breast tumors from the TCGA, METABRIC and CPTAC studies demonstrated that carnitine palmitoyltransferase 1A (CPT1A), which is localized to this amplicon, is overexpressed at the mRNA and protein level in aggressive luminal tumors, strongly associated with indicators of tumor proliferation and a predictor of poor prognosis. In vitro genetic studies demonstrated that CPT1A is required for and can promote luminal breast cancer proliferation, survival, as well as colony and mammosphere formation. Since CPT1A is the rate-limiting enzyme during fatty acid oxidation (FAO), our data indicate that FAO may be essential for these tumors. Pharmacologic inhibition of FAO prevented in vitro and in vivo tumor growth and cell proliferation as well as promoted apoptosis in luminal breast cancer cells and orthotopic xenograft tumor models. Collectively, our data establish an oncogenic role for CPT1A and FAO in HR+ luminal tumors and provide preclinical evidence and rationale supporting further investigation of FAO as a potential therapeutic opportunity for the treatment of HR+ breast cancer.

Circulating miRNAs are known to play important roles in intercellular communication. However, the effects of exosomal miRNAs on cells are not fully understood.

Breast cancer heterogeneity occurs as a consequence of the dysregulation of numerous oncogenic pathways as well as many non-genetic factors, including tumor microenvironmental stresses such as hypoxia, lactic acidosis, and glucose deprivation. Although the importance of these non-genetic factors is well recognized, it is not clear how to integrate these factors within the genetic framework of cancer as the next logical step in understanding tumor heterogeneity.

Recent studies have emphasized the importance of pathway-specific interpretations for understanding the functional relevance of gene alterations in human cancers. Although signaling activities are often conceptualized as linear events, in reality, they reflect the activity of complex functional networks assembled from modules that each respond to input signals. To acquire a deeper understanding of this network structure, we developed an approach to deconstruct pathways into modules represented by gene expression signatures. Our studies confirm that they represent units of underlying biological activity linked to known biochemical pathway structures. Importantly, we show that these signaling modules provide tools to dissect the complexity of oncogenic states that define disease outcomes as well as response to pathway-specific therapeutics. We propose that this model of pathway structure constitutes a framework to study the processes by which information propogates through cellular networks and to elucidate the relationships of fundamental modules to cellular and clinical phenotypes.

Proliferative and invasive breast tumors evolve heterogeneously in individual patients, posing significant challenges in identifying new druggable targets for precision, effective therapy. Here we present a functional multi-omics method, interaction-Correlated Multi-omic Aberration Patterning (iC-MAP), which dissects intra-tumor heterogeneity and identifies in situ the oncogenic consequences of multi-omics aberrations that drive proliferative and invasive tumors. First, we perform chromatin activity-based chemoproteomics (ChaC) experiments on breast cancer (BC) patient tissues to identify genetic/transcriptomic alterations that manifest as oncogenically active proteins. ChaC employs a biotinylated small molecule probe that specifically binds to the oncogenically active histone methyltransferase G9a, enabling sorting/enrichment of a G9a-interacting protein complex that represents the predominant BC subtype in a tissue. Second, using patient transcriptomic/genomic data, we retrospectively identified some G9a interactor-encoding genes that showed individualized iC-MAP. Our iC-MAP findings represent both new diagnostic/prognostic markers to identify patient subsets with incurable metastatic disease and targets to create individualized therapeutic strategies.

Somatic mutations have been extensively characterized in breast cancer, but the effects of these genetic alterations on the proteomic landscape remain poorly understood. Here we describe quantitative mass-spectrometry-based proteomic and phosphoproteomic analyses of 105 genomically annotated breast cancers, of which 77 provided high-quality data. Integrated analyses provided insights into the somatic cancer genome including the consequences of chromosomal loss, such as the 5q deletion characteristic of basal-like breast cancer. Interrogation of the 5q trans-effects against the Library of Integrated Network-based Cellular Signatures, connected loss of CETN3 and SKP1 to elevated expression of epidermal growth factor receptor (EGFR), and SKP1 loss also to increased SRC tyrosine kinase. Global proteomic data confirmed a stromal-enriched group of proteins in addition to basal and luminal clusters, and pathway analysis of the phosphoproteome identified a G-protein-coupled receptor cluster that was not readily identified at the mRNA level. In addition to ERBB2, other amplicon-associated highly phosphorylated kinases were identified, including CDK12, PAK1, PTK2, RIPK2 and TLK2. We demonstrate that proteogenomic analysis of breast cancer elucidates the functional consequences of somatic mutations, narrows candidate nominations for driver genes within large deletions and amplified regions, and identifies therapeutic targets.

Elucidating the molecular drivers of human breast cancers requires a strategy that is capable of integrating multiple forms of data and an ability to interpret the functional consequences of a given genetic aberration. Here we present an integrated genomic strategy based on the use of gene expression signatures of oncogenic pathway activity (n = 52) as a framework to analyze DNA copy number alterations in combination with data from a genome-wide RNA-mediated interference screen. We identify specific DNA amplifications and essential genes within these amplicons representing key genetic drivers, including known and new regulators of oncogenesis. The genes identified include eight that are essential for cell proliferation (FGD5, METTL6, CPT1A, DTX3, MRPS23, EIF2S2, EIF6 and SLC2A10) and are uniquely amplified in patients with highly proliferative luminal breast tumors, a clinical subset of patients for which few therapeutic options are effective. This general strategy has the potential to identify therapeutic targets within amplicons through an integrated use of genomic data sets.

Tumor cells rely on increased glycolytic capacity to promote cell growth and progression. While glycolysis is known to be upregulated in the majority of triple negative (TNBC) or basal-like subtype breast cancers, the mechanism remains unclear. Here, we used integrative genomic analyses to identify a subset of basal-like tumors characterized by increased expression of the oncogenic transcription factor SOX4 and its co-factor the SWI/SNF ATPase SMARCA4. These tumors are defined by unique gene expression programs that correspond with increased tumor proliferation and activation of key metabolic pathways, including glycolysis. Mechanistically, we demonstrate that the SOX4-SMARCA4 complex mediates glycolysis through direct transcriptional regulation of Hexokinase 2 (HK2) and that aberrant HK2 expression and altered glycolytic capacity are required to mediate SOX4-SMARCA4-dependent cell growth. Collectively, we have defined the SOX4-SMARCA4-HK2 signaling axis in basal-like breast tumors and established that this axis promotes metabolic reprogramming which is required to maintain tumor cell growth.

BRAF-driven colorectal cancer is among the poorest prognosis subtypes of colon cancer. Previous studies suggest that BRAF-mutant serrated cancers frequently exhibit Microsatellite Instability (MSI) and elevated levels of WNT signaling. The loss of tumor-suppressor Smad4 in oncogenic BRAF-V600E mouse models promotes rapid serrated tumor development and progression, and SMAD4 mutations co-occur in human patient tumors with BRAF-V600E mutations. This study assesses the role of SMAD4 in early-stage serrated tumorigenesis. SMAD4 loss promotes microsatellite stable (MSS) serrated tumors in an oncogenic BRAF-V600E context, providing a model for MSS serrated cancers. Inactivation of Msh2 in these mice accelerated tumor formation, and whole-exome sequencing of both MSS and MSI serrated tumors derived from these mouse models revealed that all serrated tumors developed oncogenic WNT mutations, predominantly in the WNT-effector gene Ctnnb1 (β-catenin). Mouse models mimicking the oncogenic β-catenin mutation show that the combination of three oncogenic mutations (Ctnnb1, Braf, and Smad4) are critical to drive rapid serrated dysplasia formation. Re-analysis of human tumor data reveals BRAF-V600E mutations co-occur with oncogenic mutations in both WNT and SMAD4/TGFβ pathways. These findings identify SMAD4 as a critical factor in early-stage serrated cancers and helps broaden the knowledge of this rare but aggressive subset of colorectal cancer.