SLAM-seq defines direct gene-regulatory functions of the BRD4-MYC axis

PIN1 Provides Dynamic Control of MYC in Response to Extrinsic Signals

PIN1 is a phosphorylation-directed member of the peptidyl-prolyl cis/trans isomerase (PPIase) family that facilitates conformational changes in phosphorylated targets such as c-MYC (MYC). Following signaling events that mediate phosphorylation of MYC at Serine 62, PIN1 establishes structurally distinct pools of MYC through its trans-cis and cis-trans isomerization activity at Proline 63. Through these isomerization steps, PIN1 functionally regulates MYC's stability, the molecular timing of its DNA binding and transcriptional activity, and its subnuclear localization. Recently, our group showed that Serine 62 phosphorylated MYC can associate with the inner basket of the nuclear pore (NP) in a PIN1-dependent manner. The poised euchromatin at the NP basket enables rapid cellular response to environmental signals and cell stress, and PIN1-mediated trafficking of MYC calibrates this response. In this perspective, we describe the molecular aspects of PIN1 target recognition and PIN1's function in the context of its temporal and spatial regulation of MYC.

Metabolic Profiling Reveals a Dependency of Human Metastatic Breast Cancer on Mitochondrial Serine and One-Carbon Unit Metabolism

Breast cancer is the most common cancer among American women and a major cause of mortality. To identify metabolic pathways as potential targets to treat metastatic breast cancer, we performed metabolomics profiling on the breast cancer cell line MDA-MB-231 and its tissue-tropic metastatic subclones. Here, we report that these subclones with increased metastatic potential display an altered metabolic profile compared with the parental population. In particular, the mitochondrial serine and one-carbon (1C) unit pathway is upregulated in metastatic subclones. Mechanistically, the mitochondrial serine and 1C unit pathway drives the faster proliferation of subclones through enhanced de novo purine biosynthesis. Inhibition of the first rate-limiting enzyme of the mitochondrial serine and 1C unit pathway, serine hydroxymethyltransferase (SHMT2), potently suppresses proliferation of metastatic subclones in culture and impairs growth of lung metastatic subclones at both primary and metastatic sites in mice. Some human breast cancers exhibit a significant association between the expression of genes in the mitochondrial serine and 1C unit pathway with disease outcome and higher expression of SHMT2 in metastatic tumor tissue compared with primary tumors. In addition to breast cancer, a few other cancer types, such as adrenocortical carcinoma and kidney chromophobe cell carcinoma, also display increased SHMT2 expression during disease progression. Together, these results suggest that mitochondrial serine and 1C unit metabolism plays an important role in promoting cancer progression, particularly in late-stage cancer. IMPLICATIONS: This study identifies mitochondrial serine and 1C unit metabolism as an important pathway during the progression of a subset of human breast cancers.

Inferring Causal Gene Regulatory Networks from Coupled Single-Cell Expression Dynamics Using Scribe

Here, we present Scribe (https://github.com/aristoteleo/Scribe-py), a toolkit for detecting and visualizing causal regulatory interactions between genes and explore the potential for single-cell experiments to power network reconstruction. Scribe employs restricted directed information to determine causality by estimating the strength of information transferred from a potential regulator to its downstream target. We apply Scribe and other leading approaches for causal network reconstruction to several types of single-cell measurements and show that there is a dramatic drop in performance for "pseudotime"-ordered single-cell data compared with true time-series data. We demonstrate that performing causal inference requires temporal coupling between measurements. We show that methods such as "RNA velocity" restore some degree of coupling through an analysis of chromaffin cell fate commitment. These analyses highlight a shortcoming in experimental and computational methods for analyzing gene regulation at single-cell resolution and suggest ways of overcoming it.

Targeting MYC Proteins for Tumor Therapy

Targeting the function of MYC oncoproteins holds the promise of achieving conceptually new and effective anticancer therapies that can be applied to a broad range of tumors. The nature of the target however—a broadly, possibly universally acting transcription factor that has no enzymatic activity and is largely unstructured unless complexed with partner proteins—has so far defied the development of clinically applicable MYC-directed therapies. At the same time, lingering questions about exactly which functions of MYC proteins account for their pervasive oncogenic role in human tumors and need to be targeted have prevented the development of effective therapies using surrogate targets that act in critical MYC-dependent pathways. In this review, we therefore argue that rigorous testing of critical oncogenic functions and protein/protein interactions and new chemical approaches to target them are necessary to successfully eradicate MYC-driven tumors.

Thioguanosine Conversion Enables mRNA-Lifetime Evaluation by RNA Sequencing Using Double Metabolic Labeling (TUC-seq DUAL)

Temporal information about cellular RNA populations is essential to understand the functional roles of RNA. We have developed the hydrazine/NH₄ Cl/OsO₄ -based conversion of 6-thioguanosine (6sG) into A', where A' constitutes a 6-hydrazino purine derivative. A' retains the Watson-Crick base-pair mode and is efficiently decoded as adenosine in primer extension assays and in RNA sequencing. Because 6sG is applicable to metabolic labeling of freshly synthesized RNA and because the conversion chemistry is fully compatible with the conversion of the frequently used metabolic label 4-thiouridine (4sU) into C, the combination of both modified nucleosides in dual-labeling setups enables high accuracy measurements of RNA decay. This approach, termed TUC-seq DUAL, uses the two modified nucleosides in subsequent pulses and their simultaneous detection, enabling mRNA-lifetime evaluation with unprecedented precision.

A Protocol for Transcriptome-Wide Inference of RNA Metabolic Rates in Mouse Embryonic Stem Cells

The relative ease of mouse Embryonic Stem Cells (mESCs) culture and the potential of these cells to differentiate into any of the three primary germ layers: ectoderm, endoderm and mesoderm (pluripotency), makes them an ideal and frequently used ex vivo system to dissect how gene expression changes impact cell state and differentiation. These efforts are further supported by the large number of constitutive and inducible mESC mutants established with the aim of assessing the contributions of different pathways and genes to cell homeostasis and gene regulation. Gene product abundance is controlled by the modulation of the rates of RNA synthesis, processing, and degradation. The ability to determine the relative contribution of these different RNA metabolic rates to gene expression control using standard RNA-sequencing approaches, which only capture steady state abundance of transcripts, is limited. In contrast, metabolic labeling of RNA with 4-thiouridine (4sU) coupled with RNA-sequencing, allows simultaneous and reproducible inference of transcriptome wide synthesis, processing, and degradation rates. Here we describe, a detailed protocol for 4sU metabolic labeling in mESCs that requires short 4sU labeling times at low concentration and minimally impacts cellular homeostasis. This approach presents a versatile method for in-depth characterization of the gene regulatory strategies governing gene steady state abundance in mESC.

Target gene-independent functions of MYC oncoproteins

Oncoproteins of the MYC family are major drivers of human tumorigenesis. Since a large body of evidence indicates that MYC proteins are transcription factors, studying their function has focused on the biology of their target genes. Detailed studies of MYC-dependent changes in RNA levels have provided contrasting models of the oncogenic activity of MYC proteins through either enhancing or repressing the expression of specific target genes, or as global amplifiers of transcription. In this Review, we first summarize the biochemistry of MYC proteins and what is known (or is unclear) about the MYC target genes. We then discuss recent progress in defining the interactomes of MYC and MYCN and how this information affects central concepts of MYC biology, focusing on mechanisms by which MYC proteins modulate transcription. MYC proteins promote transcription termination upon stalling of RNA polymerase II, and we propose that this mechanism enhances the stress resilience of basal transcription. Furthermore, MYC proteins coordinate transcription elongation with DNA replication and cell cycle progression. Finally, we argue that the mechanism by which MYC proteins regulate the transcription machinery is likely to promote tumorigenesis independently of global or relative changes in the expression of their target genes.

Comparative structure-function analysis of bromodomain and extraterminal motif (BET) proteins in a gene-complementation system

The widely expressed bromodomain and extraterminal motif (BET) proteins bromodomain-containing protein 2 (BRD2), BRD3, and BRD4 are multifunctional transcriptional regulators that bind acetylated chromatin via their conserved tandem bromodomains. Small molecules that target BET bromodomains are being tested for various diseases but typically do not discern between BET family members. Genomic distributions and protein partners of BET proteins have been described, but the basis for differences in BET protein function within a given lineage remains unclear. By establishing a gene knockout-rescue system in a Brd2-null erythroblast cell line, here we compared a series of mutant and chimeric BET proteins for their ability to modulate cell growth, differentiation, and gene expression. We found that the BET N-terminal halves bearing the bromodomains convey marked differences in protein stability but do not account for specificity in BET protein function. Instead, when BET proteins were expressed at comparable levels, their specificity was largely determined by the C-terminal half. Remarkably, a chimeric BET protein comprising the N-terminal half of the structurally similar short BRD4 isoform (BRD4S) and the C-terminal half of BRD2 functioned similarly to intact BRD2. We traced part of the BRD2-specific activity to a previously uncharacterized short segment predicted to harbor a coiled-coil (CC) domain. Deleting the CC segment impaired BRD2's ability to restore growth and differentiation, and the CC region functioned in conjunction with the adjacent ET domain to impart BRD2-like activity onto BRD4S. In summary, our results identify distinct BET protein domains that regulate protein turnover and biological activities.

Lisa: inferring transcriptional regulators through integrative modeling of public chromatin accessibility and ChIP-seq data

We developed Lisa (http://lisa.cistrome.org/) to predict the transcriptional regulators (TRs) of differentially expressed or co-expressed gene sets. Based on the input gene sets, Lisa first uses histone mark ChIP-seq and chromatin accessibility profiles to construct a chromatin model related to the regulation of these genes. Using TR ChIP-seq peaks or imputed TR binding sites, Lisa probes the chromatin models using in silico deletion to find the most relevant TRs. Applied to gene sets derived from targeted TF perturbation experiments, Lisa boosted the performance of imputed TR cistromes and outperformed alternative methods in identifying the perturbed TRs.

Determining mRNA Stability by Metabolic RNA Labeling and Chemical Nucleoside Conversion

The varying rates at which mRNAs decay are tightly coordinated with transcriptional changes to shape gene expression during development and disease. But currently available RNA sequencing approaches lack the temporal information to determine the relative contribution of RNA biogenesis, processing and turnover to the establishment of steady-state gene expression profiles.Here, we describe a protocol that combines metabolic RNA labeling with chemical nucleoside conversion by thiol-linked alkylation of 4-thiouridine to determine RNA stability in cultured cells (SLAMseq). When coupled to cost-effective mRNA 3' end sequencing approaches, SLAMseq determines the half-life of polyadenylated transcripts in a global and transcript-specific manner using untargeted or targeted cDNA library preparation protocols.We provide a step-by-step instruction for time-resolved mRNA 3' end sequencing, which augments traditional RNA-seq approaches to acquire the temporal resolution necessary to study the molecular principles that control gene expression.

MYC in Germinal Center-derived lymphomas: Mechanisms and therapeutic opportunities

The rearrangement of immunoglobulin loci during the germinal center reaction is associated with an increased risk of chromosomal translocations that activate oncogenes such as MYC, BCL2 or BCL6, thus contributing to the development of B-cell lymphomas. MYC and BCL2 activation are initiating events in Burkitt's (BL) and Follicular Lymphoma (FL), respectively, but can occur at later stages in other subtypes such as Diffuse Large-B Cell Lymphoma (DLBCL). MYC can also be activated during the progression of FL to the transformed stage. Thus, either DLBCL or FL can give rise to aggressive double-hit lymphomas (DHL) with concurrent activation of MYC and BCL2. Research over the last three decades has improved our understanding of the functions of these oncogenes and the basis for their cooperative action in lymphomagenesis. MYC, in particular, is a transcription factor that contributes to cell activation, growth and proliferation, while concomitantly sensitizing cells to apoptosis, the latter being blocked by BCL2. Here, we review our current knowledge about the role of MYC in germinal center B-cells and lymphomas, discuss MYC-induced dependencies that can sensitize cancer cells to select pharmacological inhibitors, and illustrate their therapeutic potential in aggressive lymphomas-and in particular in DHL, in combination with BCL2 inhibitors.

Lessons from eRNAs: understanding transcriptional regulation through the lens of nascent RNAs

Nascent transcription assays, such as global run-on sequencing (GRO-seq) and precision run-on sequencing (PRO-seq), have uncovered a myriad of unstable RNAs being actively produced from numerous sites genome-wide. These transcripts provide a more complete and immediate picture of the impact of regulatory events. Transcription factors recruit RNA polymerase II, effectively initiating the process of transcription; repressors inhibit polymerase recruitment. Efficiency of recruitment is dictated by sequence elements in and around the RNA polymerase loading zone. A combination of sequence elements and RNA binding proteins subsequently influence the ultimate stability of the resulting transcript. Some of these transcripts are capable of providing feedback on the process, influencing subsequent transcription. By monitoring RNA polymerase activity, nascent assays provide insights into every step of the regulated process of transcription.

BET Inhibition Induces HEXIM1- and RAD51-Dependent Conflicts between Transcription and Replication

BET bromodomain proteins are required for oncogenic transcription activities, and BET inhibitors have been rapidly advanced into clinical trials. Understanding the effects of BET inhibition on processes such as DNA replication will be important for future clinical applications. Here, we show that BET inhibition, and specifically inhibition of BRD4, causes replication stress through a rapid overall increase in RNA synthesis. We provide evidence that BET inhibition acts by releasing P-TEFb from its inhibitor HEXIM1, promoting interference between transcription and replication. Unusually, these transcription-replication conflicts do not activate the ATM/ATR-dependent DNA damage response but recruit the homologous recombination factor RAD51. Both HEXIM1 and RAD51 promote BET inhibitor-induced fork slowing but also prevent a DNA damage response. Our data suggest that BET inhibitors slow replication through concerted action of transcription and recombination machineries and shed light on the importance of replication stress in the action of this class of experimental cancer drugs.

The Leukemogenic TCF3-HLF Complex Rewires Enhancers Driving Cellular Identity and Self-Renewal Conferring EP300 Vulnerability

The chimeric transcription factor TCF3-HLF defines an incurable acute lymphoblastic leukemia subtype. Here we decipher the regulome of endogenous TCF3-HLF and dissect its essential transcriptional components and targets by functional genomics. We demonstrate that TCF3-HLF recruits HLF binding sites at hematopoietic stem cell/myeloid lineage associated (super-) enhancers to drive lineage identity and self-renewal. Among direct targets, hijacking an HLF binding site in a MYC enhancer cluster by TCF3-HLF activates a conserved MYC-driven transformation program crucial for leukemia propagation in vivo. TCF3-HLF pioneers the cooperation with ERG and recruits histone acetyltransferase p300 (EP300), conferring susceptibility to EP300 inhibition. Our study provides a framework for targeting driving transcriptional dependencies in this fatal leukemia.

Nascent RNA analyses: tracking transcription and its regulation

The programmes that direct an organism's development and maintenance are encoded in its genome. Decoding of this information begins with regulated transcription of genomic DNA into RNA. Although transcription and its control can be tracked indirectly by measuring stable RNAs, it is only by directly measuring nascent RNAs that the immediate regulatory changes in response to developmental, environmental, disease and metabolic signals are revealed. Multiple complementary methods have been developed to quantitatively track nascent transcription genome-wide at nucleotide resolution, all of which have contributed novel insights into the mechanisms of gene regulation and transcription-coupled RNA processing. Here we critically evaluate the array of strategies used for investigating nascent transcription and discuss the recent conceptual advances they have provided.

Combinational therapeutic targeting of BRD4 and CDK7 synergistically induces anticancer effects in head and neck squamous cell carcinoma

The bromodomain and extra-terminal domain protein BRD4 has been recognized as a key oncogenic driver and a druggable target against cancer. However, these BRD4 inhibitors as monotherapy were moderate in efficacy in preclinical models. Here we utilized a small-scale drug synergy screen that combined the BRD4 inhibitor (JQ1) with 8 epigenetic or transcriptional targeted chemicals and identified THZ1 (a CDK7 inhibitor) acting synergistically with JQ1 against head neck squamous cell carcinoma (HNSCC). Combinational JQ1 and THZ1 treatment impaired cell proliferation, induced apoptosis and senescence, which were largely recapitulated by dual BRD4 and CDK7 knockdown. Combinational treatment inhibited tumor growth and progression in 4NQO-induced HNSCC and xenograft animal models. RNA-sequencing analyses identified hundreds of differentially expressed genes modulated by JQ1 and THZ1, which were significantly enriched in categories including cell cycle and apoptosis. Mechanistically, combinational treatment reduced H3K27ac enrichment in the super-enhancer region of YAP1, which inactivated its transcription and in turn induced anti-proliferative and pro-apoptotic effects. Combined BRD4 and CDK7 upregulation associated with worst prognosis in HNSCC patients. Collectively, our findings reveal a novel therapeutic strategy of pharmacological inhibitions of BRD4 and CDK7 against HNSCC.

Non-canonical transcriptional consequences of BET inhibition in cancer

Inhibition of the bromo and extra-terminal domain (BET) protein family in preclinical studies has demonstrated that BET proteins are critical for cancer progression and important therapeutic targets. Downregulation of the MYC oncogene, CDK6, BCL2 and FOSL1 are just a few examples of the effects of BET inhibitors that can lead to cell cycle arrest and apoptosis in cancer cells. However, BET inhibitors have had little success in the clinic as a single agent, and there are an increasing number of reports of resistance to BET inhibition emerging after sustained treatment of cancer cells in vitro. Here we summarize the non-canonical consequences of BET inhibition in cancer, and discuss how these may both lead to resistance and inform rational combinations that could greatly enhance the clinical application of these inhibitors.

An efficient auxin-inducible degron system with low basal degradation in human cells

Auxin-inducible degron technology allows rapid and controlled protein depletion. However, basal degradation without auxin and inefficient auxin-inducible depletion have limited its utility. We have identified a potent auxin-inducible degron system composed of auxin receptor F-box protein AtAFB2 and short degron miniIAA7. The system showed minimal basal degradation and enabled rapid auxin-inducible depletion of endogenous human transmembrane, cytoplasmic and nuclear proteins in 1 h with robust functional phenotypes.

An early Myc-dependent transcriptional program orchestrates cell growth during B-cell activation

Upon activation, lymphocytes exit quiescence and undergo substantial increases in cell size, accompanied by activation of energy-producing and anabolic pathways, widespread chromatin decompaction, and elevated transcriptional activity. These changes depend upon prior induction of the Myc transcription factor, but how Myc controls them remains unclear. We addressed this issue by profiling the response to LPS stimulation in wild-type and c-myc-deleted primary mouse B-cells. Myc is rapidly induced, becomes detectable on virtually all active promoters and enhancers, but has no direct impact on global transcriptional activity. Instead, Myc contributes to the swift up- and down-regulation of several hundred genes, including many known regulators of the aforementioned cellular processes. Myc-activated promoters are enriched for E-box consensus motifs, bind Myc at the highest levels, and show enhanced RNA Polymerase II recruitment, the opposite being true at down-regulated loci. Remarkably, the Myc-dependent signature identified in activated B-cells is also enriched in Myc-driven B-cell lymphomas: hence, besides modulation of new cancer-specific programs, the oncogenic action of Myc may largely rely on sustained deregulation of its normal physiological targets.

Tagger-A Swiss army knife for multiomics to dissect cell type-specific mechanisms of gene expression in mice

A deep understanding of how regulation of the multiple levels of gene expression in mammalian tissues give rise to complex phenotypes has been impeded by cellular diversity. A handful of techniques were developed to tag-select nucleic acids of interest in specific cell types, thereby enabling their capture. We expanded this strategy by developing the Tagger knock-in mouse line bearing a quad-cistronic transgene combining enrichment tools for nuclei, nascent RNA, translating mRNA, and mature microRNA (miRNA). We demonstrate that Tagger can capture the desired nucleic acids, enabling multiple omics approaches to be applied to specific cell types in vivo using a single transgenic mouse line.

This Methods and Resources paper describes Tagger, a knock-in mouse line bearing a quad-cistronic transgene that enables the capture of translating mRNAs, mature miRNAs, pulse-labeled total RNA, and the nucleus, all from specific cells of complex tissues.

Chemo-enzymatic treatment of RNA to facilitate analyses

Labeling RNA is a recurring problem to make RNA compatible with state-of-the-art methodology and comes in many flavors. Considering only cellular applications, the spectrum still ranges from site-specific labeling of individual transcripts, for example, for live-cell imaging of mRNA trafficking, to metabolic labeling in combination with next generation sequencing to capture dynamic aspects of RNA metabolism on a transcriptome-wide scale. Combining the specificity of RNA-modifying enzymes with non-natural substrates has emerged as a valuable strategy to modify RNA site- or sequence-specifically with functional groups suitable for subsequent bioorthogonal reactions and thus label RNA with reporter moieties such as affinity or fluorescent tags. In this review article, we will cover chemo-enzymatic approaches (a) for in vitro labeling of RNA for application in cells, (b) for treatment of total RNA, and (c) for metabolic labeling of RNA. This article is categorized under: RNA Processing < RNA Editing and Modification RNA Methods < RNA Analyses in vitro and In Silico RNA Methods < RNA Analyses in Cells.

BORIS promotes chromatin regulatory interactions in treatment-resistant cancer cells

The CCCTC-binding factor (CTCF), which anchors DNA loops that organize the genome into structural domains, has a central role in gene control by facilitating or constraining interactions between genes and their regulatory elements^1,2. In cancer cells, the disruption of CTCF binding at specific loci by somatic mutation^3,4 or DNA hypermethylation⁵ results in the loss of loop anchors and consequent activation of oncogenes. By contrast, the germ-cell-specific paralogue of CTCF, BORIS (brother of the regulator of imprinted sites, also known as CTCFL)⁶, is overexpressed in several cancers^7-9, but its contributions to the malignant phenotype remain unclear. Here we show that aberrant upregulation of BORIS promotes chromatin interactions in ALK-mutated, MYCN-amplified neuroblastoma¹⁰ cells that develop resistance to ALK inhibition. These cells are reprogrammed to a distinct phenotypic state during the acquisition of resistance, a process defined by the initial loss of MYCN expression followed by subsequent overexpression of BORIS and a concomitant switch in cellular dependence from MYCN to BORIS. The resultant BORIS-regulated alterations in chromatin looping lead to the formation of super-enhancers that drive the ectopic expression of a subset of proneural transcription factors that ultimately define the resistance phenotype. These results identify a previously unrecognized role of BORIS-to promote regulatory chromatin interactions that support specific cancer phenotypes.

Sequencing cell-type-specific transcriptomes with SLAM-ITseq

Analysis of cell-type-specific transcriptomes is vital for understanding the biology of tissues and organs in the context of multicellular organisms. In this Protocol Extension, we combine a previously developed cell-type-specific metabolic RNA labeling method (thiouracil (TU) tagging) and a pipeline to detect the labeled transcripts by a novel RNA sequencing (RNA-seq) method, SLAMseq (thiol (SH)-linked alkylation for the metabolic sequencing of RNA). By injecting a uracil analog, 4-thiouracil, into transgenic mice that express cell-type-specific uracil phosphoribosyltransferase (UPRT), an enzyme required for 4-thiouracil incorporation into newly synthesized RNA, only cells expressing UPRT synthesize thiol-containing RNA. Total RNA isolated from a tissue of interest is then sequenced with SLAMseq, which introduces thymine to cytosine (T>C) conversions at the sites of the incorporated 4-thiouracil. The resulting sequencing reads are then mapped with the T>C-aware alignment software, SLAM-DUNK, which allows mapping of reads containing T>C mismatches. The number of T>C conversions per transcript is further analyzed to identify which transcripts are synthesized in the UPRT-expressing cells. Thus, our method, SLAM-ITseq (SLAMseq in tissue), enables cell-specific transcriptomics without laborious FACS-based cell sorting or biochemical isolation of the labeled transcripts used in TU tagging. In the murine tissues we assessed previously, this method identified ~5,000 genes that are expressed in a cell type of interest from the total RNA pool from the tissue. Any laboratory with access to a high-throughput sequencer and high-power computing can adapt this protocol with ease, and the entire pipeline can be completed in <5 d.

Acute depletion of CTCF directly affects MYC regulation through loss of enhancer-promoter looping

Numerous pieces of evidence support the complex, 3D spatial organization of the genome dictates gene expression. CTCF is essential to define topologically associated domain boundaries and to facilitate the formation of insulated chromatin loop structures. To understand CTCF's direct role in global transcriptional regulation, we integrated the miniAID-mClover3 cassette to the endogenous CTCF locus in a human pediatric B-ALL cell line, SEM, and an immortal erythroid precursor cell line, HUDEP-2, to allow for acute depletion of CTCF protein by the auxin-inducible degron system. In SEM cells, CTCF loss notably disrupted intra-TAD loops and TAD integrity in concurrence with a reduction in CTCF-binding affinity, while showing no perturbation to nuclear compartment integrity. Strikingly, the overall effect of CTCF's loss on transcription was minimal. Whole transcriptome analysis showed hundreds of genes differentially expressed in CTCF-depleted cells, among which MYC and a number of MYC target genes were specifically downregulated. Mechanically, acute depletion of CTCF disrupted the direct interaction between the MYC promoter and its distal enhancer cluster residing ∼1.8 Mb downstream. Notably, MYC expression was not profoundly affected upon CTCF loss in HUDEP-2 cells suggesting that CTCF could play a B-ALL cell line specific role in maintaining MYC expression.

NASC-seq monitors RNA synthesis in single cells

Sequencing of newly synthesised RNA can monitor transcriptional dynamics with great sensitivity and high temporal resolution, but is currently restricted to populations of cells. Here, we develop new transcriptome alkylation-dependent single-cell RNA sequencing (NASC-seq), to monitor newly synthesised and pre-existing RNA simultaneously in single cells. We validate the method on pre-labelled RNA, and by demonstrating that more newly synthesised RNA was detected for genes with known high mRNA turnover. Monitoring RNA synthesis during Jurkat T-cell activation with NASC-seq reveals both rapidly up- and down-regulated genes, and that induced genes are almost exclusively detected as newly transcribed. Moreover, the newly synthesised and pre-existing transcriptomes after T-cell activation are distinct, confirming that NASC-seq simultaneously measures gene expression corresponding to two time points in single cells. Altogether, NASC-seq enables precise temporal monitoring of RNA synthesis at single-cell resolution during homoeostasis, perturbation responses and cellular differentiation.

BRD4 Promotes Gastric Cancer Progression and Metastasis through Acetylation-Dependent Stabilization of Snail

Cancer metastasis, a leading cause of death in patients, is associated with aberrant expression of epigenetic modifiers, yet it remains poorly defined how epigenetic readers drive metastatic growth and whether epigenetic readers are targetable to control metastasis. Here, we report that bromodomain-containing protein 4 (BRD4), a histone acetylation reader and emerging anticancer therapeutic target, promotes progression and metastasis of gastric cancer. The abundance of BRD4 in human gastric cancer tissues correlated with shortened metastasis-free gastric cancer patient survival. Consistently, BRD4 maintained invasiveness of cancer cells in vitro and their dissemination at distal organs in vivo. Surprisingly, BRD4 function in this context was independent of its putative transcriptional targets such as MYC or BCL2, but rather through stabilization of Snail at posttranslational levels. In an acetylation-dependent manner, BRD4 recognized acetylated lysine 146 (K146) and K187 on Snail to prevent Snail recognition by its E3 ubiquitin ligases FBXL14 and β-Trcp1, thereby inhibiting Snail polyubiquitination and proteasomal degradation. Accordingly, genome-wide transcriptome analyses identified that BRD4 and Snail regulate a partially shared metastatic gene signature in gastric cancer cells. These findings reveal a noncanonical posttranscriptional regulatory function of BRD4 in maintaining cancer growth and dissemination, with immediate translational implications for treating gastric metastatic malignancies with clinically available bromodomain inhibitors. SIGNIFICANCE: These findings reveal a novel posttranscriptional regulatory function of the epigenetic reader BRD4 in cancer metastasis via stabilizing Snail, with immediate translational implication for treating metastatic malignancies with clinically available bromodomain inhibitors. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/79/19/4869/F1.large.jpg.

BRD4 bimodal binding at promoters and drug-induced displacement at Pol II pause sites associates with I-BET sensitivity

Background

Deregulated transcription is a major driver of diseases such as cancer. Bromodomain and extra-terminal (BET) proteins (BRD2, BRD3, BRD4 and BRDT) are chromatin readers essential for maintaining proper gene transcription by specifically binding acetylated lysine residues. Targeted displacement of BET proteins from chromatin, using BET inhibitors (I-BETs), is a promising therapy, especially for acute myeloid leukemia (AML), and evaluation of resistance mechanisms is necessary to optimize the clinical efficacy of these drugs.

Results

To uncover mechanisms of intrinsic I-BET resistance, we quantified chromatin binding and displacement for BRD2, BRD3 and BRD4 after dose response treatment with I-BET151, in sensitive and resistant in vitro models of leukemia, and mapped BET proteins/I-BET interactions genome wide using antibody- and compound-affinity capture methods followed by deep sequencing. The genome-wide map of BET proteins sensitivity to I-BET revealed a bimodal pattern of binding flanking transcription start sites (TSSs), in which drug-mediated displacement from chromatin primarily affects BRD4 downstream of the TSS and prolongs the pausing of RNA Pol II. Correlation of BRD4 binding and drug-mediated displacement at RNA Pol II pause sites with gene expression revealed a differential behavior of sensitive and resistant tumor cells to I-BET and identified a BRD4 signature at promoters of sensitive coding and non-coding genes.

Conclusions

We provide evidence that I-BET-induced shift of Pol II pausing at promoters via displacement of BRD4 is a determinant of intrinsic I-BET sensitivity. This finding may guide pharmacological treatment to enhance the clinical utility of such targeted therapies in AML and potentially other BET proteins-driven diseases.

Electronic supplementary material

The online version of this article (10.1186/s13072-019-0286-5) contains supplementary material, which is available to authorized users.

Pharmacological Modulation of Transcriptional Coregulators in Cancer

Upon binding of transcription factors to cis-regulatory DNA sequences, transcriptional coregulators are required for the activation or suppression of chromatin-dependent transcriptional signaling. These coregulators are frequently implicated in oncogenesis via causal roles in dysregulated, malignant transcriptional control and represent one of the fastest-growing target classes in small-molecule drug discovery. However, challenges in targeting coregulators include identifying evidence of cancer-specific genetic dependency, matching the pharmacologically addressable protein fold to a functional role in disease pathology, and achieving the necessary selectivity to exploit a given genetic dependency. We discuss here how recent trends in cancer pharmacology have confronted these challenges, positioning coregulators as tractable targets in the development of new cancer therapies.

Quantification of experimentally induced nucleotide conversions in high-throughput sequencing datasets

BACKGROUND

Methods to read out naturally occurring or experimentally introduced nucleic acid modifications are emerging as powerful tools to study dynamic cellular processes. The recovery, quantification and interpretation of such events in high-throughput sequencing datasets demands specialized bioinformatics approaches.

RESULTS

Here, we present Digital Unmasking of Nucleotide conversions in K-mers (DUNK), a data analysis pipeline enabling the quantification of nucleotide conversions in high-throughput sequencing datasets. We demonstrate using experimentally generated and simulated datasets that DUNK allows constant mapping rates irrespective of nucleotide-conversion rates, promotes the recovery of multimapping reads and employs Single Nucleotide Polymorphism (SNP) masking to uncouple true SNPs from nucleotide conversions to facilitate a robust and sensitive quantification of nucleotide-conversions. As a first application, we implement this strategy as SLAM-DUNK for the analysis of SLAMseq profiles, in which 4-thiouridine-labeled transcripts are detected based on T > C conversions. SLAM-DUNK provides both raw counts of nucleotide-conversion containing reads as well as a base-content and read coverage normalized approach for estimating the fractions of labeled transcripts as readout.

CONCLUSION

Beyond providing a readily accessible tool for analyzing SLAMseq and related time-resolved RNA sequencing methods (TimeLapse-seq, TUC-seq), DUNK establishes a broadly applicable strategy for quantifying nucleotide conversions.

Transcriptional Reprogramming and Novel Therapeutic Approaches for Targeting Prostate Cancer Stem Cells

Prostate cancer is the most common malignancy in men and the second cause of cancer-related deaths in western countries. Despite the progress in the treatment of localized prostate cancer, there is still lack of effective therapies for the advanced forms of the disease. Most patients with advanced prostate cancer become resistant to androgen deprivation therapy (ADT), which remains the main therapeutic option in this setting, and progress to lethal metastatic castration-resistant prostate cancer (mCRPC). Current therapies for prostate cancer preferentially target proliferating, partially differentiated, and AR-dependent cancer cells that constitute the bulk of the tumor mass. However, the subpopulation of tumor-initiating or tumor-propagating stem-like cancer cells is virtually resistant to the standard treatments causing tumor relapse at the primary or metastatic sites. Understanding the pathways controlling the establishment, expansion and maintenance of the cancer stem cell (CSC) subpopulation is an important step toward the development of more effective treatment for prostate cancer, which might enable ablation or exhaustion of CSCs and prevent treatment resistance and disease recurrence. In this review, we focus on the impact of transcriptional regulators on phenotypic reprogramming of prostate CSCs and provide examples supporting the possibility of inhibiting maintenance and expansion of the CSC pool in human prostate cancer along with the currently available methodological approaches. Transcription factors are key elements for instructing specific transcriptional programs and inducing CSC-associated phenotypic changes implicated in disease progression and treatment resistance. Recent studies have shown that interfering with these processes causes exhaustion of CSCs with loss of self-renewal and tumorigenic capability in prostate cancer models. Targeting key transcriptional regulators in prostate CSCs is a valid therapeutic strategy waiting to be tested in clinical trials.

Functional-genetic approaches to understanding drug response and resistance

Drug development remains a slow and expensive process, while the effective use of established therapeutics is widely hampered by our limited understanding of response and resistance mechanisms. Functional-genetic tools such as CRISPR/Cas9, advanced RNAi methods, and targeted protein degradation, together with other emerging technologies such as time-resolved and single-cell transcriptomics, fundamentally change the way we can search for candidate therapeutic targets and evaluate them before drug development. In addition, for already available therapeutics these tools open vast opportunities for probing response mechanisms and predictive biomarkers, and thereby guide the development of personalized therapies. Here, we review promising applications and remaining limitations of recently established functional-genetic tools for high-throughput screening and the in-depth analysis of candidate targets and established drugs.

Identification and characterization of cancer vulnerabilities via targeted protein degradation

Target(ed) protein degradation (TPD) is a novel paradigm in drug discovery and a promising therapeutic strategy. TPD is based on small-molecules that catalyze the degradation of proteins by re-directing the ubiquitination activity of ubiquitin E3 ligases. Its unique molecular pharmacology enables robust, selective and fast elimination of proteins in cellular assays and in vivo. In addition to possible clinical applications, TPD is also emerging as an attractive alternative to traditional pharmacologic or genetic perturbation strategies. Directly acting degraders, as well as chemical-genetics derivatives offer unique opportunities in the pre-clinical identification, characterization and mechanistic validation of therapeutic targets.

MYC Recruits SPT5 to RNA Polymerase II to Promote Processive Transcription Elongation

The MYC oncoprotein binds to promoter-proximal regions of virtually all transcribed genes and enhances RNA polymerase II (Pol II) function, but its precise mode of action is poorly understood. Using mass spectrometry of both MYC and Pol II complexes, we show here that MYC controls the assembly of Pol II with a small set of transcription elongation factors that includes SPT5, a subunit of the elongation factor DSIF. MYC directly binds SPT5, recruits SPT5 to promoters, and enables the CDK7-dependent transfer of SPT5 onto Pol II. Consistent with known functions of SPT5, MYC is required for fast and processive transcription elongation. Intriguingly, the high levels of MYC that are expressed in tumors sequester SPT5 into non-functional complexes, thereby decreasing the expression of growth-suppressive genes. Altogether, these results argue that MYC controls the productive assembly of processive Pol II elongation complexes and provide insight into how oncogenic levels of MYC permit uncontrolled cellular growth.

•

MYC enhances productive transcription by defining the protein composition of Pol II

•

MYC directly binds SPT5 and hands it over to Pol II in a CDK7-dependent manner

•

Transfer of SPT5 increases speed and processivity of Pol II

•

MYC’s effects on Pol II function shape its tumor-specific gene expression profile

Recruitment of BRCA1 limits MYCN-driven accumulation of stalled RNA polymerase

MYC is an oncogenic transcription factor that binds globally to active promoters and promotes transcriptional elongation by RNA polymerase II (RNAPII)^1,2. Deregulated expression of the paralogous protein MYCN drives the development of neuronal and neuroendocrine tumours and is often associated with a particularly poor prognosis³. Here we show that, similar to MYC, activation of MYCN in human neuroblastoma cells induces escape of RNAPII from promoters. If the release of RNAPII from transcriptional pause sites (pause release) fails, MYCN recruits BRCA1 to promoter-proximal regions. Recruitment of BRCA1 prevents MYCN-dependent accumulation of stalled RNAPII and enhances transcriptional activation by MYCN. Mechanistically, BRCA1 stabilizes mRNA decapping complexes and enables MYCN to suppress R-loop formation in promoter-proximal regions. Recruitment of BRCA1 requires the ubiquitin-specific protease USP11, which binds specifically to MYCN when MYCN is dephosphorylated at Thr58. USP11, BRCA1 and MYCN stabilize each other on chromatin, preventing proteasomal turnover of MYCN. Because BRCA1 is highly expressed in neuronal progenitor cells during early development⁴ and MYC is less efficient than MYCN in recruiting BRCA1, our findings indicate that a cell-lineage-specific stress response enables MYCN-driven tumours to cope with deregulated RNAPII function.

Recruitment of BRCA1 limits MYCN-driven accumulation of stalled RNA polymerase

MYC is an oncogenic transcription factor that binds globally to active promoters and promotes transcriptional elongation by RNA polymerase II (RNAPII)^1,2. Deregulated expression of the paralogous protein MYCN drives the development of neuronal and neuroendocrine tumours and is often associated with a particularly poor prognosis³. Here we show that, similar to MYC, activation of MYCN in human neuroblastoma cells induces escape of RNAPII from promoters. If the release of RNAPII from transcriptional pause sites (pause release) fails, MYCN recruits BRCA1 to promoter-proximal regions. Recruitment of BRCA1 prevents MYCN-dependent accumulation of stalled RNAPII and enhances transcriptional activation by MYCN. Mechanistically, BRCA1 stabilizes mRNA decapping complexes and enables MYCN to suppress R-loop formation in promoter-proximal regions. Recruitment of BRCA1 requires the ubiquitin-specific protease USP11, which binds specifically to MYCN when MYCN is dephosphorylated at Thr58. USP11, BRCA1 and MYCN stabilize each other on chromatin, preventing proteasomal turnover of MYCN. Because BRCA1 is highly expressed in neuronal progenitor cells during early development⁴ and MYC is less efficient than MYCN in recruiting BRCA1, our findings indicate that a cell-lineage-specific stress response enables MYCN-driven tumours to cope with deregulated RNAPII function.

Therapeutic Targeting of BRD4 in Head Neck Squamous Cell Carcinoma

The bromodomain and extraterminal family members are epigenetic readers and transcriptional coactivators which are critically involved in various biological processes including tumorigenesis. BRD4 has been increasingly appreciated as a key oncogene and promising anticancer target. Here, we sought to characterize the expression of BRD4 and its tumorigenic roles as well as therapeutic targeting in HNSCC.

Methods: Expression of BRD4 mRNA and protein was determined by bioinformatics interrogation of publically available databases, primary HNSCC samples and 4NQO-induced HNSCC animal model. The tumorigenic roles of BRD4 in HNSCC were evaluated by genetic and pharmacological approach in vitro and in vivo. Therapeutic efficiency of BRD4 targeting by JQ1 was assessed in three preclinical models including xenograft model, 4NQO-induced model and patients-derived xenograft model. Gene candidates responsible for therapeutic effects of JQ1 were identified by transcriptional profiling in HNSCC cells after JQ1 exposure.

Results: Significant upregulation of BRD4 was found in primary HNSCC samples and 4NQO-induced HNSCC model. Its overexpression associated with aggressive clinicopathological features and inferior overall and disease-free survival. BRD4 depletion by genetic silencing or pharmacological inhibition impaired cell proliferation, migration and invasion and reduced tumor growth and metastasis in vivo. Transcriptional profiling of HNSCC cells following JQ1 exposure identified hundreds of genes which might mediated its antitumor effects and enriched in cancer-relevant pathways. A novel prognostic risk score derived from JQ1-regulated genes was developed to stratify patients into subgroups with favorable or inferior prognosis.

Conclusions: Our findings reveal that BRD4 serves as a novel and critical mediator underlying tumorigenesis and a robust prognostic biomarker in HNSCC. Therapeutic targeting of BRD4 represents a potent and promising strategy against HNSCC.

Long noncoding RNA expression profiling in cancer: Challenges and opportunities

In recent years, technological advances in transcriptome profiling revealed that the repertoire of human RNA molecules is more diverse and extended than originally thought. This diversity and complexity mainly derive from a large ensemble of noncoding RNAs. Because of their key roles in cellular processes important for normal development and physiology, disruption of noncoding RNA expression is intrinsically linked to human disease, including cancer. Therefore, studying the noncoding portion of the transcriptome offers the prospect of identifying novel therapeutic and diagnostic targets. Although evidence of the relevance of noncoding RNAs in cancer is accumulating, we still face many challenges when it comes to accurately profiling their expression levels. Some of these challenges are inherent to the technologies employed, whereas others are associated with characteristics of the noncoding RNAs themselves. In this review, we discuss the challenges related to long noncoding RNA expression profiling, highlight how cancer long noncoding RNAs provide new opportunities for cancer diagnosis and treatment, and reflect on future developments.

Synergistic activity of BET inhibitor MK-8628 and PLK inhibitor Volasertib in preclinical models of medulloblastoma

Medulloblastoma is the most prevalent central nervous system tumor in children. Targeted treatment approaches for patients with high-risk medulloblastoma are needed as current treatment regimens are not curative in many cases and cause significant therapy-related morbidity. Medulloblastoma harboring MYC amplification have the most aggressive clinical course and worst outcome. Targeting the BET protein BRD4 has significant anti-tumor effects in preclinical models of MYC-amplified medulloblastoma, however, in most cases these are not curative. We here assessed the therapeutic efficacy of the orally bioavailable BRD4 inhibitor, MK-8628, in preclinical models of medulloblastoma. MK-8628 showed therapeutic efficacy against in vitro and in vivo models of MYC-amplified medulloblastoma by inducing apoptotic cell death and cell cycle arrest. Gene expression analysis of cells treated with MK-8628 showed that anti-tumor effects were accompanied by significant repression of MYC transcription as well as disruption of MYC-regulated transcriptional programs. Additionally, we found that targeting of MYC protein stability through pharmacological PLK1 inhibition showed synergistic anti-medulloblastoma effects when combined with MK-8628 treatment. Thus, MK-8628 is effective against preclinical high-risk medulloblastoma models and its effects can be enhanced through simultaneous targeting of PLK1.

Interactome Rewiring Following Pharmacological Targeting of BET Bromodomains

Targeting bromodomains (BRDs) of the bromo-and-extra-terminal (BET) family offers opportunities for therapeutic intervention in cancer and other diseases. Here, we profile the interactomes of BRD2, BRD3, BRD4, and BRDT following treatment with the pan-BET BRD inhibitor JQ1, revealing broad rewiring of the interaction landscape, with three distinct classes of behavior for the 603 unique interactors identified. A group of proteins associate in a JQ1-sensitive manner with BET BRDs through canonical and new binding modes, while two classes of extra-terminal (ET)-domain binding motifs mediate acetylation-independent interactions. Last, we identify an unexpected increase in several interactions following JQ1 treatment that define negative functions for BRD3 in the regulation of rRNA synthesis and potentially RNAPII-dependent gene expression that result in decreased cell proliferation. Together, our data highlight the contributions of BET protein modules to their interactomes allowing for a better understanding of pharmacological rewiring in response to JQ1.

Global role for coactivator complexes in RNA polymerase II transcription

SAGA and TFIID are related transcription complexes, which were proposed to alternatively deliver TBP at different promoter classes. Recent genome-wide studies in yeast revealed that both complexes are required for the transcription of a vast majority of genes by RNA polymerase II raising new questions about the role of coactivators.

Contributions of regulated transcription and mRNA decay to the dynamics of gene expression

Organisms have acquired sophisticated regulatory networks that control gene expression in response to cellular perturbations. Understanding of the mechanisms underlying the coordinated changes in gene expression in response to external and internal stimuli is a fundamental issue in biology. Recent advances in high-throughput technologies have enabled the measurement of diverse biological information, including gene expression levels, kinetics of gene expression, and interactions among gene expression regulatory molecules. By coupling these technologies with quantitative modeling, we can now uncover the biological roles and mechanisms of gene regulation at the system level. This review consists of two parts. First, we focus on the methods using uridine analogs that measure synthesis and decay rates of RNAs, which demonstrate how cells dynamically change the regulation of gene expression in response to both internal and external cues. Second, we discuss the underlying mechanisms of these changes in kinetics, including the functions of transcription factors and RNA-binding proteins. Overall, this review will help to clarify a system-level view of gene expression programs in cells. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Turnover and Surveillance > Regulation of RNA Stability RNA Methods > RNA Analyses in vitro and In Silico.

Eukaryotic core promoters and the functional basis of transcription initiation

RNA polymerase II (Pol II) core promoters are specialized DNA sequences at transcription start sites of protein-coding and non-coding genes that support the assembly of the transcription machinery and transcription initiation. They enable the highly regulated transcription of genes by selectively integrating regulatory cues from distal enhancers and their associated regulatory proteins. In this Review, we discuss the defining properties of gene core promoters, including their sequence features, chromatin architecture and transcription initiation patterns. We provide an overview of molecular mechanisms underlying the function and regulation of core promoters and their emerging functional diversity, which defines distinct transcription programmes. On the basis of the established properties of gene core promoters, we discuss transcription start sites within enhancers and integrate recent results obtained from dedicated functional assays to propose a functional model of transcription initiation. This model can explain the nature and function of transcription initiation at gene starts and at enhancers and can explain the different roles of core promoters, of Pol II and its associated factors and of the activating cues provided by enhancers and the transcription factors and cofactors they recruit.

SLAM-ITseq: sequencing cell type-specific transcriptomes without cell sorting

Cell type-specific transcriptome analysis is an essential tool for understanding biological processes in which diverse types of cells are involved. Although cell isolation methods such as fluorescence-activated cell sorting (FACS) in combination with transcriptome analysis have widely been used so far, their time-consuming and harsh procedures limit their applications. Here, we report a novel in vivo metabolic RNA sequencing method, SLAM-ITseq, which metabolically labels RNA with 4-thiouracil in a specific cell type in vivo followed by detection through an RNA-seq-based method that specifically distinguishes the thiolated uridine by base conversion. This method has successfully identified the cell type-specific transcriptome in three different tissues: endothelial cells in brain, epithelial cells in intestine and adipocytes in white adipose tissue. As this method does not require isolation of cells or RNA prior to the transcriptomic analysis, SLAM-ITseq provides an easy yet accurate snapshot of the transcriptional state in vivo.

MYC-induced metabolic stress and tumorigenesis

The MYC oncogene is commonly altered across human cancers. Distinct from the normal MYC proto-oncogene, which is under tight transcriptional, translational, and post-translational control, deregulated oncogenic MYC drives imbalanced, non-linear amplification of transcription that results in oncogenic 'stress.' The term 'stress' had been a euphemism for our lack of mechanistic understanding, but synthesis of many studies over the past decade provides a more coherent picture of oncogenic MYC driving metastable cellular states, particularly altered metabolism, that activate and depend on cellular stress response pathways to allow for continued growth and survival. Both deregulated metabolism and these stress response pathways represent vulnerabilities that can be exploited therapeutically.