BEERS2: RNA-Seq simulation through high fidelity in silico modeling.

Simulation of RNA-seq reads is critical in the assessment, comparison, benchmarking and development of bioinformatics tools. Yet the field of RNA-seq simulators has progressed little in the last decade. To address this need we have developed BEERS2, which combines a flexible and highly configurable design with detailed simulation of the entire library preparation and sequencing pipeline. BEERS2 takes input transcripts (typically fully length messenger RNA transcripts with polyA tails) from either customizable input or from CAMPAREE simulated RNA samples. It produces realistic reads of these transcripts as FASTQ, SAM or BAM formats with the SAM or BAM formats containing the true alignment to the reference genome. It also produces true transcript-level quantification values. BEERS2 combines a flexible and highly configurable design with detailed simulation of the entire library preparation and sequencing pipeline and is designed to include the effects of polyA selection and RiboZero for ribosomal depletion, hexamer priming sequence biases, GC-content biases in polymerase chain reaction (PCR) amplification, barcode read errors and errors during PCR amplification. These characteristics combine to make BEERS2 the most complete simulation of RNA-seq to date. Finally, we demonstrate the use of BEERS2 by measuring the effect of several settings on the popular Salmon pseudoalignment algorithm.

Identification of cancer-cytotoxic modulators of PDE3A by predictive chemogenomics.

High cancer death rates indicate the need for new anticancer therapeutic agents. Approaches to discovering new cancer drugs include target-based drug discovery and phenotypic screening. Here, we identified phosphodiesterase 3A modulators as cell-selective cancer cytotoxic compounds through phenotypic compound library screening and target deconvolution by predictive chemogenomics. We found that sensitivity to 6-(4-(diethylamino)-3-nitrophenyl)-5-methyl-4,5-dihydropyridazin-3(2H)-one, or DNMDP, across 766 cancer cell lines correlates with expression of the gene PDE3A, encoding phosphodiesterase 3A. Like DNMDP, a subset of known PDE3A inhibitors kill selected cancer cells, whereas others do not. Furthermore, PDE3A depletion leads to DNMDP resistance. We demonstrated that DNMDP binding to PDE3A promotes an interaction between PDE3A and Schlafen 12 (SLFN12), suggestive of a neomorphic activity. Coexpression of SLFN12 with PDE3A correlates with DNMDP sensitivity, whereas depletion of SLFN12 results in decreased DNMDP sensitivity. Our results implicate PDE3A modulators as candidate cancer therapeutic agents and demonstrate the power of predictive chemogenomics in small-molecule discovery.

The interaction between Myc and Miz1 is required to antagonize TGFbeta-dependent autocrine signaling during lymphoma formation and maintenance.

The Myc protein suppresses the transcription of several cyclin-dependent kinase inhibitors (CKIs) via binding to Miz1; whether this interaction is important for Myc's ability to induce or maintain tumorigenesis is not known. Here we show that the oncogenic potential of a point mutant of Myc (MycV394D) that is selectively deficient in binding to Miz1 is greatly attenuated. Binding of Myc to Miz1 is continuously required to repress CKI expression and inhibit accumulation of trimethylated histone H3 at Lys 9 (H3K9triMe), a hallmark of cellular senescence, in T-cell lymphomas. Lymphomas that arise express high amounts of transforming growth factor beta-2 (TGFbeta-2) and TGFbeta-3. Upon Myc suppression, TGFbeta signaling is required to induce CKI expression and cellular senescence and suppress tumor recurrence. Binding of Myc to Miz1 is required to antagonize growth suppression and induction of senescence by TGFbeta. We demonstrate that, since lymphomas express high levels of TGFbeta, they are poised to elicit an autocrine program of senescence upon Myc inactivation, demonstrating that TGFbeta is a key factor that establishes oncogene addiction of T-cell lymphomas.

Addiction to multiple oncogenes can be exploited to prevent the emergence of therapeutic resistance.

Many cancers exhibit sensitivity to the inhibition of a single genetic lesion, a property that has been successfully exploited with oncogene-targeted therapeutics. However, inhibition of single oncogenes often fails to result in sustained tumor regression due to the emergence of therapy-resistant cells. Here, we report that MYC-driven lymphomas frequently acquire activating mutations in ß-catenin, including a previously unreported mutation in a splice acceptor site. Tumors with these genetic lesions are highly dependent on ß-catenin for their survival and the suppression of ß-catenin resulted in marked apoptosis causally related to a decrease in Bcl-xL expression. Using a novel inducible inhibitor of ß-catenin, we illustrate that, although MYC withdrawal or ß-catenin inhibition alone results in initial tumor regression, most tumors ultimately recurred, mimicking the clinical response to single-agent targeted therapy. Importantly, the simultaneous combined inhibition of both MYC and ß-catenin promoted more rapid tumor regression and successfully prevented tumor recurrence. Hence, we demonstrated that MYC-induced tumors are addicted to mutant ß-catenin, and the combined inactivation of MYC and ß-catenin induces sustained tumor regression. Our results provide a proof of principle that targeting multiple oncogene addicted pathways can prevent therapeutic resistance.

Lymphomas that recur after MYC suppression continue to exhibit oncogene addiction.

The suppression of oncogenic levels of MYC is sufficient to induce sustained tumor regression associated with proliferative arrest, differentiation, cellular senescence, and/or apoptosis, a phenomenon known as oncogene addiction. However, after prolonged inactivation of MYC in a conditional transgenic mouse model of Eµ-tTA/tetO-MYC T-cell acute lymphoblastic leukemia, some of the tumors recur, recapitulating what is frequently observed in human tumors in response to targeted therapies. Here we report that these recurring lymphomas express either transgenic or endogenous Myc, albeit in many cases at levels below those in the original tumor, suggesting that tumors continue to be addicted to MYC. Many of the recurring lymphomas (76%) harbored mutations in the tetracycline transactivator, resulting in expression of the MYC transgene even in the presence of doxycycline. Some of the remaining recurring tumors expressed high levels of endogenous Myc, which was associated with a genomic rearrangement of the endogenous Myc locus or activation of Notch1. By gene expression profiling, we confirmed that the primary and recurring tumors have highly similar transcriptomes. Importantly, shRNA-mediated suppression of the high levels of MYC in recurring tumors elicited both suppression of proliferation and increased apoptosis, confirming that these tumors remain oncogene addicted. These results suggest that tumors induced by MYC remain addicted to overexpression of this oncogene.

BEERS2: RNA-Seq simulation through high fidelity in silico modeling.

Simulation of RNA-seq reads is critical in the assessment, comparison, benchmarking, and development of bioinformatics tools. Yet the field of RNA-seq simulators has progressed little in the last decade. To address this need we have developed BEERS2, which combines a flexible and highly configurable design with detailed simulation of the entire library preparation and sequencing pipeline. BEERS2 takes input transcripts (typically fully-length mRNA transcripts with polyA tails) from either customizable input or from CAMPAREE simulated RNA samples. It produces realistic reads of these transcripts as FASTQ, SAM, or BAM formats with the SAM or BAM formats containing the true alignment to the reference genome. It also produces true transcript-level quantification values. BEERS2 combines a flexible and highly configurable design with detailed simulation of the entire library preparation and sequencing pipeline and is designed to include the effects of polyA selection and RiboZero for ribosomal depletion, hexamer priming sequence biases, GC-content biases in PCR amplification, barcode read errors, and errors during PCR amplification. These characteristics combine to make BEERS2 the most complete simulation of RNA-seq to date. Finally, we demonstrate the use of BEERS2 by measuring the effect of several settings on the popular Salmon pseudoalignment algorithm.

A genome-scale CRISPR screen reveals PRMT1 as a critical regulator of androgen receptor signaling in prostate cancer.

Androgen receptor (AR) signaling is the central driver of prostate cancer across disease states. While androgen deprivation therapy (ADT) is effective in the initial treatment of prostate cancer, resistance to ADT or to next-generation androgen pathway inhibitors invariably arises, most commonly through the re-activation of the AR axis. Thus, orthogonal approaches to inhibit AR signaling in advanced prostate cancer are essential. Here, via genome-scale CRISPR-Cas9 screening, we identify protein arginine methyltransferase 1 (PRMT1) as a critical mediator of AR expression and signaling. PRMT1 regulates the recruitment of AR to genomic target sites and the inhibition of PRMT1 impairs AR binding at lineage-specific enhancers, leading to decreased expression of key oncogenes, including AR itself. In addition, AR-driven prostate cancer cells are uniquely susceptible to combined AR and PRMT1 inhibition. Our findings implicate PRMT1 as a key regulator of AR output and provide a preclinical framework for co-targeting of AR and PRMT1 in advanced prostate cancer.

Sequential translocation of an Escherchia coli two-partner secretion pathway exoprotein across the inner and outer membranes.

In Gram-negative bacteria, a variety of high molecular weight 'exoproteins' are translocated across the outer membrane (OM) via the two-partner secretion (TPS) pathway by interacting with a dedicated transporter. It is unclear, however, whether the translocation of exoproteins across the OM is coupled to their translocation across the inner membrane (IM). To address this question, we separated the production of an Escherichia coli O157:H7 exoprotein (OtpA) and its transporter (OtpB) temporally by placing otpA and otpB under the control of distinct regulatable promoters. We found that when both full-length and truncated forms of OtpA were expressed prior to OtpB, a significant fraction of the exoprotein was secreted. The results indicate that OtpA can be translocated into the periplasm and briefly remain secretion-competent. Furthermore, by engineering cysteine residues into OtpA and using disulphide bond formation as a reporter of periplasmic localization, we obtained additional evidence that the C-terminus of OtpA enters the periplasm before the N-terminus is translocated across the OM even when OtpA and OtpB are expressed simultaneously. Taken together, our results demonstrate that the translocation of a TPS exoprotein across the OM can occur independently from its translocation across the IM.

RNA-binding proteins of COSMIC importance in cancer.

Herculean efforts by the Wellcome Sanger Institute, the National Cancer Institute, and the National Human Genome Research Institute to sequence thousands of tumors representing all major cancer types have yielded more than 700 genes that contribute to neoplastic growth when mutated, amplified, or deleted. While some of these genes (now included in the COSMIC Cancer Gene Census) encode proteins previously identified in hypothesis-driven experiments (oncogenic transcription factors, protein kinases, etc.), additional classes of cancer drivers have emerged, perhaps none more surprisingly than RNA-binding proteins (RBPs). Over 40 RBPs responsible for virtually all aspects of RNA metabolism, from synthesis to degradation, are recurrently mutated in cancer, and just over a dozen are considered major cancer drivers. This Review investigates whether and how their RNA-binding activities pertain to their oncogenic functions. Focusing on several well-characterized steps in RNA metabolism, we demonstrate that for virtually all cancer-driving RBPs, RNA processing activities are either abolished (the loss-of-function phenotype) or carried out with low fidelity (the LoFi phenotype). Conceptually, this suggests that in normal cells, RBPs act as gatekeepers maintaining proper RNA metabolism and the "balanced" proteome. From the practical standpoint, at least some LoFi phenotypes create therapeutic vulnerabilities, which are beginning to be exploited in the clinic.

A predominant enhancer co-amplified with the SOX2 oncogene is necessary and sufficient for its expression in squamous cancer.

Amplification and overexpression of the SOX2 oncogene represent a hallmark of squamous cancers originating from diverse tissue types. Here, we find that squamous cancers selectively amplify a 3' noncoding region together with SOX2, which harbors squamous cancer-specific chromatin accessible regions. We identify a single enhancer e1 that predominantly drives SOX2 expression. Repression of e1 in SOX2-high cells causes collapse of the surrounding enhancers, remarkable reduction in SOX2 expression, and a global transcriptional change reminiscent of SOX2 knockout. The e1 enhancer is driven by a combination of transcription factors including SOX2 itself and the AP-1 complex, which facilitates recruitment of the co-activator BRD4. CRISPR-mediated activation of e1 in SOX2-low cells is sufficient to rebuild the e1-SOX2 loop and activate SOX2 expression. Our study shows that squamous cancers selectively amplify a predominant enhancer to drive SOX2 overexpression, uncovering functional links among enhancer activation, chromatin looping, and lineage-specific copy number amplifications of oncogenes.

Chromatin Looping Shapes KLF5-Dependent Transcriptional Programs in Human Epithelial Cancers.

Activation of transcription factors is a key driver event in cancer. We and others have recently reported that the Krüppel-like transcription factor KLF5 is activated in multiple epithelial cancer types including squamous cancer and gastrointestinal adenocarcinoma, yet the functional consequences and the underlying mechanisms of this activation remain largely unknown. Here we demonstrate that activation of KLF5 results in strongly selective KLF5 dependency for these cancer types. KLF5 bound lineage-specific regulatory elements and activated gene expression programs essential to cancer cells. HiChIP analysis revealed that multiple distal KLF5 binding events cluster and synergize to activate individual target genes. Immunoprecipitation-mass spectrometry assays showed that KLF5 interacts with other transcription factors such as TP63 and YAP1, as well as the CBP/EP300 acetyltransferase complex. Furthermore, KLF5 guided the CBP/EP300 complex to increase acetylation of H3K27, which in turn enhanced recruitment of the bromodomain protein BRD4 to chromatin. The 3D chromatin architecture aggregated KLF5-dependent BRD4 binding to activate polymerase II elongation at KLF5 target genes, which conferred a transcriptional vulnerability to proteolysis-targeting chimera-induced degradation of BRD4. Our study demonstrates that KLF5 plays an essential role in multiple epithelial cancers by activating cancer-related genes through 3D chromatin loops, providing an evidence-based rationale for targeting the KLF5 pathway. SIGNIFICANCE: An integrative 3D genomics methodology delineates mechanisms underlying the function of KLF5 in multiple epithelial cancers and suggests potential strategies to target cancers with aberrantly activated KLF5.

Multi-Omics Analysis Identifies MGA as a Negative Regulator of the MYC Pathway in Lung Adenocarcinoma.

Genomic analysis of lung adenocarcinomas has revealed that the MGA gene, which encodes a heterodimeric partner of the MYC-interacting protein MAX, is significantly mutated or deleted in lung adenocarcinomas. Most of the mutations are loss of function for MGA, suggesting that MGA may act as a tumor suppressor. Here, we characterize both the molecular and cellular role of MGA in lung adenocarcinomas and illustrate its functional relevance in the MYC pathway. Although MGA and MYC interact with the same binding partner, MAX, and recognize the same E-box DNA motif, we show that the molecular function of MGA appears to be antagonistic to that of MYC. Using mass spectrometry-based affinity proteomics, we demonstrate that MGA interacts with a noncanonical PCGF6-PRC1 complex containing MAX and E2F6 that is involved in gene repression, while MYC is not part of this MGA complex, in agreement with previous studies describing the interactomes of E2F6 and PCGF6. Chromatin immunoprecipitation-sequencing and RNA sequencing assays show that MGA binds to and represses genes that are bound and activated by MYC. In addition, we show that, as opposed to the MYC oncoprotein, MGA acts as a negative regulator for cancer cell proliferation. Our study defines a novel MYC/MAX/MGA pathway, in which MYC and MGA play opposite roles in protein interaction, transcriptional regulation, and cellular proliferation. IMPLICATIONS: This study expands the range of key cancer-associated genes whose dysregulation is functionally equivalent to MYC activation and places MYC within a linear pathway analogous to cell-cycle or receptor tyrosine kinase/RAS/RAF pathways in lung adenocarcinomas.

Identification and Characterization of Oncogenic SOS1 Mutations in Lung Adenocarcinoma.

Lung adenocarcinomas are characterized by mutations in the receptor tyrosine kinase (RTK)/Ras/Raf pathway, with up to 75% of cases containing mutations in known driver genes. However, the driver alterations in the remaining cases are yet to be determined. Recent exome sequencing analysis has identified SOS1, encoding a guanine nucleotide exchange factor, as significantly mutated in lung adenocarcinomas lacking canonical oncogenic RTK/Ras/Raf pathway mutations. Here, we demonstrate that ectopic expression of lung adenocarcinoma-derived mutants of SOS1 induces anchorage-independent cell growth in vitro and tumor formation in vivo. Biochemical experiments suggest that these mutations lead to overactivation of the Ras pathway, which can be suppressed by mutations that disrupt either the Ras-GEF or putative Rac-GEF activity of SOS1. Transcriptional profiling reveals that the expression of mutant SOS1 leads to the upregulation of MYC target genes and genes associated with Ras transformation. Furthermore, we demonstrate that an AML cancer cell line harboring a lung adenocarcinoma-associated mutant SOS1 is dependent on SOS1 for survival and is also sensitive to MEK inhibition. Our work provides experimental evidence for the role of SOS1 as an oncogene and suggests a possible therapeutic strategy to target SOS1-mutated cancers. IMPLICATIONS: This study demonstrates that SOS1 mutations found in lung adenocarcinoma are oncogenic and that MEK inhibition may be a therapeutic avenue for the treatment of SOS1-mutant cancers.

Somatic Superenhancer Duplications and Hotspot Mutations Lead to Oncogenic Activation of the KLF5 Transcription Factor.

The Krüppel-like family of transcription factors plays critical roles in human development and is associated with cancer pathogenesis. Krüppel-like factor 5 gene (KLF5) has been shown to promote cancer cell proliferation and tumorigenesis and to be genomically amplified in cancer cells. We recently reported that the KLF5 gene is also subject to other types of somatic coding and noncoding genomic alterations in diverse cancer types. Here, we show that these alterations activate KLF5 by three distinct mechanisms: (i) Focal amplification of superenhancers activates KLF5 expression in squamous cell carcinomas; (ii) Missense mutations disrupt KLF5-FBXW7 interactions to increase KLF5 protein stability in colorectal cancer; (iii) Cancer type-specific hotspot mutations within a zinc-finger DNA binding domain of KLF5 change its DNA binding specificity and reshape cellular transcription. Utilizing data from CRISPR/Cas9 gene knockout screening, we reveal that cancer cells with KLF5 overexpression are dependent on KLF5 for their proliferation, suggesting KLF5 as a putative therapeutic target.Significance: Our observations, together with previous studies that identified oncogenic properties of KLF5, establish the importance of KLF5 activation in human cancers, delineate the varied genomic mechanisms underlying this occurrence, and nominate KLF5 as a putative target for therapeutic intervention in cancer. Cancer Discov; 8(1); 108-25. ©2017 AACR.This article is highlighted in the In This Issue feature, p. 1.

Identification of ADAR1 adenosine deaminase dependency in a subset of cancer cells.

Systematic exploration of cancer cell vulnerabilities can inform the development of novel cancer therapeutics. Here, through analysis of genome-scale loss-of-function datasets, we identify adenosine deaminase acting on RNA (ADAR or ADAR1) as an essential gene for the survival of a subset of cancer cell lines. ADAR1-dependent cell lines display increased expression of interferon-stimulated genes. Activation of type I interferon signaling in the context of ADAR1 deficiency can induce cell lethality in non-ADAR1-dependent cell lines. ADAR deletion causes activation of the double-stranded RNA sensor, protein kinase R (PKR). Disruption of PKR signaling, through inactivation of PKR or overexpression of either a wildtype or catalytically inactive mutant version of the p150 isoform of ADAR1, partially rescues cell lethality after ADAR1 loss, suggesting that both catalytic and non-enzymatic functions of ADAR1 may contribute to preventing PKR-mediated cell lethality. Together, these data nominate ADAR1 as a potential therapeutic target in a subset of cancers.

Genome-scale analysis identifies paralog lethality as a vulnerability of chromosome 1p loss in cancer.

Functional redundancy shared by paralog genes may afford protection against genetic perturbations, but it can also result in genetic vulnerabilities due to mutual interdependency1-5. Here, we surveyed genome-scale short hairpin RNA and CRISPR screening data on hundreds of cancer cell lines and identified MAGOH and MAGOHB, core members of the splicing-dependent exon junction complex, as top-ranked paralog dependencies6-8. MAGOHB is the top gene dependency in cells with hemizygous MAGOH deletion, a pervasive genetic event that frequently occurs due to chromosome 1p loss. Inhibition of MAGOHB in a MAGOH-deleted context compromises viability by globally perturbing alternative splicing and RNA surveillance. Dependency on IPO13, an importin-ß receptor that mediates nuclear import of the MAGOH/B-Y14 heterodimer9, is highly correlated with dependency on both MAGOH and MAGOHB. Both MAGOHB and IPO13 represent dependencies in murine xenografts with hemizygous MAGOH deletion. Our results identify MAGOH and MAGOHB as reciprocal paralog dependencies across cancer types and suggest a rationale for targeting the MAGOHB-IPO13 axis in cancers with chromosome 1p deletion.

An alternative splicing switch in FLNB promotes the mesenchymal cell state in human breast cancer.

Alternative splicing of mRNA precursors represents a key gene expression regulatory step and permits the generation of distinct protein products with diverse functions. In a genome-scale expression screen for inducers of the epithelial-to-mesenchymal transition (EMT), we found a striking enrichment of RNA-binding proteins. We validated that QKI and RBFOX1 were necessary and sufficient to induce an intermediate mesenchymal cell state and increased tumorigenicity. Using RNA-seq and eCLIP analysis, we found that QKI and RBFOX1 coordinately regulated the splicing and function of the actin-binding protein FLNB, which plays a causal role in the regulation of EMT. Specifically, the skipping of FLNB exon 30 induced EMT by releasing the FOXC1 transcription factor. Moreover, skipping of FLNB exon 30 is strongly associated with EMT gene signatures in basal-like breast cancer patient samples. These observations identify a specific dysregulation of splicing, which regulates tumor cell plasticity and is frequently observed in human cancer.

Copy-number and gene dependency analysis reveals partial copy loss of wild-type SF3B1 as a novel cancer vulnerability.

Genomic instability is a hallmark of human cancer, and results in widespread somatic copy number alterations. We used a genome-scale shRNA viability screen in human cancer cell lines to systematically identify genes that are essential in the context of particular copy-number alterations (copy-number associated gene dependencies). The most enriched class of copy-number associated gene dependencies was CYCLOPS (Copy-number alterations Yielding Cancer Liabilities Owing to Partial losS) genes, and spliceosome components were the most prevalent. One of these, the pre-mRNA splicing factor SF3B1, is also frequently mutated in cancer. We validated SF3B1 as a CYCLOPS gene and found that human cancer cells harboring partial SF3B1 copy-loss lack a reservoir of SF3b complex that protects cells with normal SF3B1 copy number from cell death upon partial SF3B1 suppression. These data provide a catalog of copy-number associated gene dependencies and identify partial copy-loss of wild-type SF3B1 as a novel, non-driver cancer gene dependency.

BIM mediates oncogene inactivation-induced apoptosis in multiple transgenic mouse models of acute lymphoblastic leukemia.

Oncogene inactivation in both clinical targeted therapies and conditional transgenic mouse cancer models can induce significant tumor regression associated with the robust induction of apoptosis. Here we report that in MYC-, RAS-, and BCR-ABL-induced acute lymphoblastic leukemia (ALL), apoptosis upon oncogene inactivation is mediated by the same pro-apoptotic protein, BIM. The induction of BIMin the MYC- and RAS-driven leukemia is mediated by the downregulation of miR-17-92. Overexpression of miR-17-92 blocked the induction of apoptosis upon oncogene inactivation in the MYC and RAS-driven but not in the BCR-ABL-driven ALL leukemia. Hence, our results provide novel insight into the mechanism of apoptosis upon oncogene inactivation and suggest that induction of BIM-mediated apoptosis may be an important therapeutic approach for ALL.