Machine learning-optimized targeted detection of alternative splicing.

RNA sequencing (RNA-seq) is widely adopted for transcriptome analysis but has inherent biases that hinder the comprehensive detection and quantification of alternative splicing. To address this, we present an efficient targeted RNA-seq method that greatly enriches for splicing-informative junction-spanning reads. Local splicing variation sequencing (LSV-seq) utilizes multiplexed reverse transcription from highly scalable pools of primers anchored near splicing events of interest.

Machine learning-optimized targeted detection of alternative splicing.

RNA-sequencing (RNA-seq) is widely adopted for transcriptome analysis but has inherent biases which hinder the comprehensive detection and quantification of alternative splicing. To address this, we present an efficient targeted RNA-seq method that greatly enriches for splicing-informative junction-spanning reads. Local Splicing Variation sequencing (LSV-seq) utilizes multiplexed reverse transcription from highly scalable pools of primers anchored near splicing events of interest.

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: RNA-Seq simulation through high fidelity 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.

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.

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.

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).

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.

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 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.

Identification and Characterization of Oncogenic 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 , encoding a guanine nucleotide exchange factor, as significantly mutated in lung adenocarcinomas lacking canonical oncogenic RTK/Ras/Raf pathway mutations.

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.

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.

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 interdependency. 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 dependencies. MAGOHB is the top gene dependency in cells with hemizygous MAGOH deletion, a pervasive genetic event that frequently occurs due to chromosome 1p loss.

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 () has been shown to promote cancer cell proliferation and tumorigenesis and to be genomically amplified in cancer cells. We recently reported that the gene is also subject to other types of somatic coding and noncoding genomic alterations in diverse cancer types.

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.

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.

Identification of focally amplified lineage-specific super-enhancers in human epithelial cancers.

Whole-genome analysis approaches are identifying recurrent cancer-associated somatic alterations in noncoding DNA regions. We combined somatic copy number analysis of 12 tumor types with tissue-specific epigenetic profiling to identify significant regions of focal amplification harboring super-enhancers. Copy number gains of noncoding regions harboring super-enhancers near KLF5, USP12, PARD6B and MYC are associated with overexpression of these cancer-related genes.

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.

miR-17-92 explains MYC oncogene addiction.

MYC regulates tumorigenesis by coordinating the expression of thousands of genes. We found that MYC appears to regulate the decisions between cell survival versus death and self-renewal versus senescence through the microRNA miR-17-92 cluster. Addiction to the MYC oncogene may therefore in fact be an addiction to miR-17-92.

Activation of Cre recombinase alone can induce complete tumor regression.

The Cre/loxP system is a powerful tool for generating conditional genomic recombination and is often used to examine the mechanistic role of specific genes in tumorigenesis. However, Cre toxicity due to its non-specific endonuclease activity has been a concern. Here, we report that tamoxifen-mediated Cre activation in vivo induced the regression of primary lymphomas in p53-/- mice.

MYC through miR-17-92 suppresses specific target genes to maintain survival, autonomous proliferation, and a neoplastic state.

The MYC oncogene regulates gene expression through multiple mechanisms, and its overexpression culminates in tumorigenesis. MYC inactivation reverses turmorigenesis through the loss of distinguishing features of cancer, including autonomous proliferation and survival. Here we report that MYC via miR-17-92 maintains a neoplastic state through the suppression of chromatin regulatory genes Sin3b, Hbp1, Suv420h1, and Btg1, as well as the apoptosis regulator Bim.

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.