Mitogen-activated protein kinase (MAPK) pathways form the backbone of signal transduction in the mammalian cell. Here we applied a systematic experimental and computational approach to map 2,269 interactions between human MAPK-related proteins and other cellular machinery and to assemble these data into functional modules. Multiple lines of evidence including conservation with yeast supported a core network of 641 interactions. Using small interfering RNA knockdowns, we observed that approximately one-third of MAPK-interacting proteins modulated MAPK-mediated signaling. We uncovered the Na-H exchanger NHE1 as a potential MAPK scaffold, found links between HSP90 chaperones and MAPK pathways and identified MUC12 as the human analog to the yeast signaling mucin Msb2. This study makes available a large resource of MAPK interactions and clone libraries, and it illustrates a methodology for probing signaling networks based on functional refinement of experimentally derived protein-interaction maps.

Influenza A virus is an RNA virus that encodes up to 11 proteins and this small coding capacity demands that the virus use the host cellular machinery for many aspects of its life cycle. Knowledge of these host cell requirements not only informs us of the molecular pathways exploited by the virus but also provides further targets that could be pursued for antiviral drug development. Here we use an integrative systems approach, based on genome-wide RNA interference screening, to identify 295 cellular cofactors required for early-stage influenza virus replication. Within this group, those involved in kinase-regulated signalling, ubiquitination and phosphatase activity are the most highly enriched, and 181 factors assemble into a highly significant host-pathogen interaction network. Moreover, 219 of the 295 factors were confirmed to be required for efficient wild-type influenza virus growth, and further analysis of a subset of genes showed 23 factors necessary for viral entry, including members of the vacuolar ATPase (vATPase) and COPI-protein families, fibroblast growth factor receptor (FGFR) proteins, and glycogen synthase kinase 3 (GSK3)-beta. Furthermore, 10 proteins were confirmed to be involved in post-entry steps of influenza virus replication. These include nuclear import components, proteases, and the calcium/calmodulin-dependent protein kinase (CaM kinase) IIbeta (CAMK2B). Notably, growth of swine-origin H1N1 influenza virus is also dependent on the identified host factors, and we show that small molecule inhibitors of several factors, including vATPase and CAMK2B, antagonize influenza virus replication.

Tumor cells have an altered metabolic phenotype characterized by increased glycolysis and diminished oxidative phosphorylation. Despite the suspected importance of glycolysis in tumorigenesis, the temporal relationship between oncogene signaling, in vivo tumor formation, and glycolytic pathway activity is poorly understood. Moreover, how glycolytic pathways are altered as tumors regress remains unknown. Here, we use a switchable model of Myc-driven liver cancer, along with hyperpolarized (13)C-pyruvate magnetic resonance spectroscopic imaging (MRSI) to visualize glycolysis in de novo tumor formation and regression. LDHA abundance and activity in tumors is tightly correlated to in vivo pyruvate conversion to lactate and is rapidly inhibited as tumors begin to regress, as are numerous glycolysis pathway genes. Conversion of pyruvate to alanine predominates in precancerous tissues prior to observable morphologic or histological changes. These results demonstrate that metabolic changes precede tumor formation and regression and are directly linked to the activity of a single oncogene.

Although oncogene-targeted therapy often elicits profound initial tumor responses in patients, responses are generally incomplete because some tumor cells survive initial therapy as residual disease that enables eventual acquired resistance. The mechanisms underlying tumor cell adaptation and survival during initial therapy are incompletely understood. Here, through the study of EGFR mutant lung adenocarcinoma, we show that NF-κB signaling is rapidly engaged upon initial EGFR inhibitor treatment to promote tumor cell survival and residual disease. EGFR oncogene inhibition induced an EGFR-TRAF2-RIP1-IKK complex that stimulated an NF-κB-mediated transcriptional survival program. The direct NF-κB inhibitor PBS-1086 suppressed this adaptive survival program and increased the magnitude and duration of initial EGFR inhibitor response in multiple NSCLC models, including a patient-derived xenograft. These findings unveil NF-κB activation as a critical adaptive survival mechanism engaged by EGFR oncogene inhibition and provide rationale for EGFR and NF-κB co-inhibition to eliminate residual disease and enhance patient responses.

Embryonic stem cells (ESCs) have adopted an accelerated cell-cycle program with shortened gap phases and precocious expression of cell-cycle regulatory proteins, including cyclins and cyclin-dependent kinases (CDKs). We examined the effect of CDK inhibition on the pathways regulating proliferation and survival of ESCs. We found that inhibiting cyclin-dependent kinase 1 (CDK1) leads to activation of the DNA damage response, nuclear p53 stabilization, activation of a subset of p53 target genes including NOXA, and negative regulation of the anti-apoptotic protein MCL1 in human and mouse ESCs, but not differentiated cells. We demonstrate that MCL1 is highly expressed in ESCs and loss of MCL1 leads to ESC death. Finally, we show that clinically relevant CDK1 inhibitors prevent formation of ESC-derived tumors and induce necrosis in established ESC-derived tumors. Our data demonstrate that ES cells are uniquely sensitive to CDK1 inhibition via a p53/NOXA/MCL1 pathway.

Most eukaryotic cells express small regulatory RNAs. The purpose of one class, the somatic endogenous siRNAs (endo-siRNAs), remains unclear. Here, we show that the endo-siRNA pathway promotes odor adaptation in C. elegans AWC olfactory neurons. In adaptation, the nuclear Argonaute NRDE-3, which acts in AWC, is loaded with siRNAs targeting odr-1, a gene whose downregulation is required for adaptation. Concomitant with increased odr-1 siRNA in AWC, we observe increased binding of the HP1 homolog HPL-2 at the odr-1 locus in AWC and reduced odr-1 mRNA in adapted animals. Phosphorylation of HPL-2, an in vitro substrate of the EGL-4 kinase that promotes adaption, is necessary and sufficient for behavioral adaptation. Thus, environmental stimulation amplifies an endo-siRNA negative feedback loop to dynamically repress cognate gene expression and shape behavior. This class of siRNA may act broadly as a rheostat allowing prolonged stimulation to dampen gene expression and promote cellular memory formation. PAPERFLICK:

Hepatocellular carcinoma (HCC) is a cancer with multiple aetiologies and widespread prevalence. Largely refractory to current treatments, HCC is the fourth leading cause of cancer-related deaths worldwide. MicroRNAs (miRNAs) are important regulators in HCCs. We aimed to identify tumour suppressor miRNAs during tumour regression in a conditional c-MYC-driven mouse model (LT2/MYC) of HCC, and to evaluate their therapeutic potential for HCC treatment.

Cancer cells utilize multiple nutrient scavenging mechanisms to support growth and survival in nutrient-poor, hypoxic tumor microenvironments. Among these mechanisms, macropinocytosis has emerged as an important pathway of extracellular nutrient acquisition in cancer cells, particularly in tumors with activated RAS signaling, such as pancreatic cancer. However, the absence of a clinically available inhibitor, as well as the gap of knowledge in macropinocytosis regulation, remain a hurdle for its use for cancer therapy. Here, we use the Informer set library to identify novel regulators of macropinocytosis-dependent growth in pancreatic cancer cells. Understanding how these regulators function will allow us to provide novel opportunities for therapeutic intervention.

Activation of the PI3K-mTOR pathway is central to breast cancer pathogenesis including resistance to many targeted therapies. The mTOR kinase forms two distinct complexes, mTORC1 and mTORC2, and understanding which is required for the survival of malignant cells has been limited by tools to selectively and completely impair either subcomplex. To address this, we used RMC-6272, a bi-steric molecule with a rapamycin-like moiety linked to an mTOR active-site inhibitor that displays >25-fold selectivity for mTORC1 over mTORC2 substrates. Complete suppression of mTORC1 by RMC-6272 causes apoptosis in ER+/HER2- breast cancer cell lines, particularly in those that harbor mutations in PIK3CA or PTEN, due to inhibition of the rapamycin resistant, mTORC1 substrate 4EBP1 and reduction of the pro-survival protein MCL1. RMC-6272 reduced translation of ribosomal mRNAs, MYC target genes, and components of the CDK4/6 pathway, suggesting enhanced impairment of oncogenic pathways compared to the partial mTORC1 inhibitor everolimus. RMC-6272 maintained efficacy in hormone therapy-resistant acquired cell lines and patient-derived xenografts (PDX), showed increased efficacy in CDK4/6 inhibitor treated acquired resistant cell lines versus their parental counterparts, and was efficacious in a PDX from a patient experiencing resistance to CDK4/6 inhibition. Bi-steric mTORC1-selective inhibition may be effective in overcoming multiple forms of therapy-resistance in ER+ breast cancers.

The mitotic spindle is the bipolar, microtubule-based structure that segregates chromosomes at each cell division. Aberrant spindles are frequently observed in cancer cells, but how oncogenic transformation affects spindle mechanics and function, particularly in the mechanical context of solid tumors, remains poorly understood. Here, we constitutively overexpress the oncogene cyclin D1 in human MCF10A cells to probe its effects on spindle architecture and response to compressive force. We find that cyclin D1 overexpression increases the incidence of spindles with extra poles, centrioles, and chromosomes. However, it also protects spindle poles from fracturing under compressive force, a deleterious outcome linked to multipolar cell divisions. Our findings suggest that cyclin D1 overexpression may adapt cells to increased compressive stress, contributing to its prevalence in cancers such as breast cancer by allowing continued proliferation in mechanically challenging environments.

Although targeted therapies often elicit profound initial patient responses, these effects are transient due to residual disease leading to acquired resistance. How tumors transition between drug responsiveness, tolerance and resistance, especially in the absence of preexisting subclones, remains unclear. In epidermal growth factor receptor (EGFR)-mutant lung adenocarcinoma cells, we demonstrate that residual disease and acquired resistance in response to EGFR inhibitors requires Aurora kinase A (AURKA) activity. Nongenetic resistance through the activation of AURKA by its coactivator TPX2 emerges in response to chronic EGFR inhibition where it mitigates drug-induced apoptosis. Aurora kinase inhibitors suppress this adaptive survival program, increasing the magnitude and duration of EGFR inhibitor response in preclinical models. Treatment-induced activation of AURKA is associated with resistance to EGFR inhibitors in vitro, in vivo and in most individuals with EGFR-mutant lung adenocarcinoma. These findings delineate a molecular path whereby drug resistance emerges from drug-tolerant cells and unveils a synthetic lethal strategy for enhancing responses to EGFR inhibitors by suppressing AURKA-driven residual disease and acquired resistance.

Expression of the oncogenic transcription factor MYC is disproportionately elevated in triple-negative breast cancer (TNBC), as compared to estrogen receptor-, progesterone receptor- or human epidermal growth factor 2 receptor-positive (RP) breast cancer. We and others have shown that MYC alters metabolism during tumorigenesis. However, the role of MYC in TNBC metabolism remains mostly unexplored. We hypothesized that MYC-dependent metabolic dysregulation is essential for the growth of MYC-overexpressing TNBC cells and may identify new therapeutic targets for this clinically challenging subset of breast cancer. Using a targeted metabolomics approach, we identified fatty acid oxidation (FAO) intermediates as being dramatically upregulated in a MYC-driven model of TNBC. We also identified a lipid metabolism gene signature in patients with TNBC that were identified from The Cancer Genome Atlas database and from multiple other clinical data sets, implicating FAO as a dysregulated pathway that is critical for TNBC cell metabolism. We found that pharmacologic inhibition of FAO catastrophically decreased energy metabolism in MYC-overexpressing TNBC cells and blocked tumor growth in a MYC-driven transgenic TNBC model and in a MYC-overexpressing TNBC patient-derived xenograft. These findings demonstrate that MYC-overexpressing TNBC shows an increased bioenergetic reliance on FAO and identify the inhibition of FAO as a potential therapeutic strategy for this subset of breast cancer.

Genome-wide association studies have identified over 70 single-nucleotide polymorphisms (SNPs) associated with breast cancer. A subset of these SNPs are associated with quantitative expression of nearby genes, but the functional effects of the majority remain unknown. We hypothesized that some risk SNPs may regulate alternative splicing. Using RNA-sequencing data from breast tumors and germline genotypes from The Cancer Genome Atlas, we tested the association between each risk SNP genotype and exon-, exon-exon junction- or transcript-specific expression of nearby genes. Six SNPs were associated with differential transcript expression of seven nearby genes at FDR < 0.05 (BABAM1, DCLRE1B/PHTF1, PEX14, RAD51L1, SRGAP2D and STXBP4). We next developed a Bayesian approach to evaluate, for each SNP, the overlap between the signal of association with breast cancer and the signal of association with alternative splicing. At one locus (SRGAP2D), this method eliminated the possibility that the breast cancer risk and the alternate splicing event were due to the same causal SNP. Lastly, at two loci, we identified the likely causal SNP for the alternative splicing event, and at one, functionally validated the effect of that SNP on alternative splicing using a minigene reporter assay. Our results suggest that the regulation of differential transcript isoform expression is the functional mechanism of some breast cancer risk SNPs and that we can use these associations to identify causal SNPs, target genes and the specific transcripts that may mediate breast cancer risk.

Acute myeloid leukemia (AML) is the most common acute leukemia in adults and is often resistant to conventional therapies. The MYC oncogene is commonly overexpressed in AML but has remained an elusive target. We aimed to examine the consequences of targeting MYC both directly and indirectly in AML overexpressing MYC/Myc due to trisomy 8/15 (human/mouse), FLT3-ITD mutation, or gene amplification.

Histologically normal tissue adjacent to the tumor (NAT) is commonly used as a control in cancer studies. However, little is known about the transcriptomic profile of NAT, how it is influenced by the tumor, and how the profile compares with non-tumor-bearing tissues. Here, we integrate data from the Genotype-Tissue Expression project and The Cancer Genome Atlas to comprehensively analyze the transcriptomes of healthy, NAT, and tumor tissues in 6506 samples across eight tissues and corresponding tumor types. Our analysis shows that NAT presents a unique intermediate state between healthy and tumor. Differential gene expression and protein-protein interaction analyses reveal altered pathways shared among NATs across tissue types. We characterize a set of 18 genes that are specifically activated in NATs. By applying pathway and tissue composition analyses, we suggest a pan-cancer mechanism of pro-inflammatory signals from the tumor stimulates an inflammatory response in the adjacent endothelium.

Breast cancer arises from breast epithelial cells that acquire genetic alterations leading to subsequent loss of tissue homeostasis. Several distinct epithelial subpopulations have been proposed, but complete understanding of the spectrum of heterogeneity and differentiation hierarchy in the human breast remains elusive. Here, we use single-cell mRNA sequencing (scRNAseq) to profile the transcriptomes of 25,790 primary human breast epithelial cells isolated from reduction mammoplasties of seven individuals. Unbiased clustering analysis reveals the existence of three distinct epithelial cell populations, one basal and two luminal cell types, which we identify as secretory L1- and hormone-responsive L2-type cells. Pseudotemporal reconstruction of differentiation trajectories produces one continuous lineage hierarchy that closely connects the basal lineage to the two differentiated luminal branches. Our comprehensive cell atlas provides insights into the cellular blueprint of the human breast epithelium and will form the foundation to understand how the system goes awry during breast cancer.

Estrogen, progesterone, and HER2 receptor-negative triple-negative breast cancers encompass the most clinically challenging subtype for which targeted therapeutics are lacking. We find that triple-negative tumors exhibit elevated MYC expression, as well as altered expression of MYC regulatory genes, resulting in increased activity of the MYC pathway. In primary breast tumors, MYC signaling did not predict response to neoadjuvant chemotherapy but was associated with poor prognosis. We exploit the increased MYC expression found in triple-negative breast cancers by using a synthetic-lethal approach dependent on cyclin-dependent kinase (CDK) inhibition. CDK inhibition effectively induced tumor regression in triple-negative tumor xenografts. The proapoptotic BCL-2 family member BIM is up-regulated after CDK inhibition and contributes to this synthetic-lethal mechanism. These results indicate that aggressive breast tumors with elevated MYC are uniquely sensitive to CDK inhibitors.

We used a quantitative, high-density genetic interaction map, or E-MAP (Epistatic MiniArray Profile), to interrogate the relationships within and between RNA-processing pathways. Due to their complexity and the essential roles of many of the components, these pathways have been difficult to functionally dissect. Here, we report the results for 107,155 individual interactions involving 552 mutations, 166 of which are hypomorphic alleles of essential genes. Our data enabled the discovery of links between components of the mRNA export and splicing machineries and Sem1/Dss1, a component of the 19S proteasome. In particular, we demonstrate that Sem1 has a proteasome-independent role in mRNA export as a functional component of the Sac3-Thp1 complex. Sem1 also interacts with Csn12, a component of the COP9 signalosome. Finally, we show that Csn12 plays a role in pre-mRNA splicing, which is independent of other signalosome components. Thus, Sem1 is involved in three separate and functionally distinct complexes.

Inferences about protein function are often made based on sequence homology to other gene products of known activities. This approach is valuable for small families of conserved proteins but can be difficult to apply to large superfamilies of proteins with diverse function. In this study we looked at sequence homology between members of the DJ-1/ThiJ/PfpI superfamily, which includes a human protein of unclear function, DJ-1, associated with inherited Parkinson's disease.

Herpesviruses constitute a family of large DNA viruses widely spread in vertebrates and causing a variety of different diseases. They possess dsDNA genomes ranging from 120 to 240 kbp encoding between 70 to 170 open reading frames. We previously reported the protein interaction networks of two herpesviruses, varicella-zoster virus (VZV) and Kaposi's sarcoma-associated herpesvirus (KSHV). In this study, we systematically tested three additional herpesvirus species, herpes simplex virus 1 (HSV-1), murine cytomegalovirus and Epstein-Barr virus, for protein interactions in order to be able to perform a comparative analysis of all three herpesvirus subfamilies. We identified 735 interactions by genome-wide yeast-two-hybrid screens (Y2H), and, together with the interactomes of VZV and KSHV, included a total of 1,007 intraviral protein interactions in the analysis. Whereas a large number of interactions have not been reported previously, we were able to identify a core set of highly conserved protein interactions, like the interaction between HSV-1 UL33 with the nuclear egress proteins UL31/UL34. Interactions were conserved between orthologous proteins despite generally low sequence similarity, suggesting that function may be more conserved than sequence. By combining interactomes of different species we were able to systematically address the low coverage of the Y2H system and to extract biologically relevant interactions which were not evident from single species.