Isoform and pathway-specific regulation of post-transcriptional RNA processing in human cells.

Steady-state levels of RNA transcripts are controlled by their rates of synthesis and degradation. Here we used nascent RNA Bru-seq and BruChase-seq to profile RNA dynamics across 16 human cell lines as part of ENCODE4 Deeply Profiled Cell Lines collection. We show that RNA turnover dynamics differ widely between transcripts of different genes and between different classes of RNA. Gene set enrichment analysis (GSEA) revealed that transcripts encoding proteins belonging to the same pathway often show similar turnover dynamics. Furthermore, transcript isoforms show distinct dynamics suggesting that RNA turnover is important in regulating mRNA isoform choice. Finally, splicing across newly made transcripts appears to be cooperative with either all or none type splicing. These data sets generated as part of ENCODE4 illustrate the intricate and coordinated regulation of RNA dynamics in controlling gene expression to allow for the precise coordination of cellular functions.

Characterizing nascent transcription patterns of PROMPTs, eRNAs, and readthrough transcripts in the ENCODE4 deeply profiled cell lines.

Arising as co-products of canonical gene expression, transcription-associated lincRNAs, such as promoter upstream transcripts (PROMPTs), enhancer RNAs (eRNAs), and readthrough (RT) transcripts, are often regarded as byproducts of transcription, although they may be important for the expression of nearby genes. We identified regions of nascent expression of these lincRNA in 16 human cell lines using Bru-seq techniques, and found distinctly regulated patterns of PROMPT, eRNA, and RT transcription using the diverse biochemical approaches in the ENCODE4 deeply profiled cell lines collection. Transcription of these lincRNAs was influenced by sequence-specific features and the local or 3D chromatin landscape. However, these sequence and chromatin features do not describe the full spectrum of lincRNA expression variability we identify, highlighting the complexity of their regulation. This may suggest that transcription-associated lincRNAs are not merely byproducts, but rather that the transcript itself, or the act of its transcription, is important for genomic function.

KMT2D links TGF-ß signaling to noncanonical activin pathway and regulates pancreatic cancer cell plasticity.

Although KMT2D, also known as MLL2, is known to play an essential role in development, differentiation, and tumor suppression, its role in pancreatic cancer development is not well understood. Here, we discovered a novel signaling axis mediated by KMT2D, which links TGF-ß to the activin A pathway. We found that TGF-ß upregulates a microRNA, miR-147b, which in turn leads to post-transcriptional silencing of KMT2D. Loss of KMT2D induces the expression and secretion of activin A, which activates a noncanonical p38 MAPK-mediated pathway to modulate cancer cell plasticity, promote a mesenchymal phenotype, and enhance tumor invasion and metastasis in mice. We observed a decreased KMT2D expression in human primary and metastatic pancreatic cancer. Furthermore, inhibition or knockdown of activin A reversed the protumoral role of KMT2D loss. These findings support a tumor-suppressive role of KMT2D in pancreatic cancer and identify miR-147b and activin A as novel therapeutic targets.

GRHL2-controlled gene expression networks in luminal breast cancer.

Grainyhead like 2 (GRHL2) is an essential transcription factor for development and function of epithelial tissues. It has dual roles in cancer by supporting tumor growth while suppressing epithelial to mesenchymal transitions (EMT). GRHL2 cooperates with androgen and estrogen receptors (ER) to regulate gene expression. We explore genome wide GRHL2 binding sites conserved in three ER⍺/GRHL2 positive luminal breast cancer cell lines by ChIP-Seq. Interaction with the ER⍺/FOXA1/GATA3 complex is observed, however, only for a minor fraction of conserved GRHL2 peaks. We determine genome wide transcriptional dynamics in response to loss of GRHL2 by nascent RNA Bru-seq using an MCF7 conditional knockout model. Integration of ChIP- and Bru-seq pinpoints candidate direct GRHL2 target genes in luminal breast cancer. Multiple connections between GRHL2 and proliferation are uncovered, including transcriptional activation of ETS and E2F transcription factors. Among EMT-related genes, direct regulation of CLDN4 is corroborated but several targets identified in other cells (including CDH1 and ZEB1) are ruled out by both ChIP- and Bru-seq as being directly controlled by GRHL2 in luminal breast cancer cells. Gene clusters correlating positively (including known GRHL2 targets such as ErbB3, CLDN4/7) or negatively (including TGFB1 and TGFBR2) with GRHL2 in the MCF7 knockout model, display similar correlation with GRHL2 in ER positive as well as ER negative breast cancer patients. Altogether, this study uncovers gene sets regulated directly or indirectly by GRHL2 in luminal breast cancer, identifies novel GRHL2-regulated genes, and points to distinct GRHL2 regulation of EMT in luminal breast cancer cells. Video Abstract.

miRNA analysis reveals novel dysregulated pathways in amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease. Its complex pathogenesis and phenotypic heterogeneity hinder therapeutic development and early diagnosis. Altered RNA metabolism is a recurrent pathophysiologic theme, including distinct microRNA (miRNA) profiles in ALS tissues. We profiled miRNAs in accessible biosamples, including skin fibroblasts and whole blood and compared them in age- and sex-matched healthy controls versus ALS participants with and without repeat expansions to chromosome 9 open reading frame 72 (C9orf72; C9-ALS and nonC9-ALS), the most frequent ALS mutation. We identified unique and shared profiles of differential miRNA (DmiRNA) levels in each C9-ALS and nonC9-ALS tissues versus controls. Fibroblast DmiRNAs were validated by quantitative real-time PCR and their target mRNAs by 5-bromouridine and 5-bromouridine-chase sequencing. We also performed pathway analysis to infer biological meaning, revealing anticipated, tissue-specific pathways and pathways previously linked to ALS, as well as novel pathways that could inform future research directions. Overall, we report a comprehensive study of a miRNA profile dataset from C9-ALS and nonC9-ALS participants across two accessible biosamples, providing evidence of dysregulated miRNAs in ALS and possible targets of interest. Distinct miRNA patterns in accessible tissues may also be leveraged to distinguish ALS participants from healthy controls for earlier diagnosis. Future directions may look at potential correlations of miRNA profiles with clinical parameters.

Divergent Molecular and Cellular Responses to Low and High-Dose Ionizing Radiation.

Cancer risk after ionizing radiation (IR) is assumed to be linear with the dose; however, for low doses, definite evidence is lacking. Here, using temporal multi-omic systems analyses after a low (LD; 0.1 Gy) or a high (HD; 1 Gy) dose of X-rays, we show that, although the DNA damage response (DDR) displayed dose proportionality, many other molecular and cellular responses did not. Phosphoproteomics uncovered a novel mode of phospho-signaling via S12-PPP1R7, and large-scale dephosphorylation events that regulate mitotic exit control in undamaged cells and the G2/M checkpoint upon IR in a dose-dependent manner. The phosphoproteomics of irradiated DNA double-strand breaks (DSBs) repair-deficient cells unveiled extended phospho-signaling duration in either a dose-dependent (DDR signaling) or independent (mTOR-ERK-MAPK signaling) manner without affecting signal magnitude. Nascent transcriptomics revealed the transcriptional activation of genes involved in NRF2-regulated antioxidant defense, redox-sensitive ERK-MAPK signaling, glycolysis and mitochondrial function after LD, suggesting a prominent role for reactive oxygen species (ROS) in molecular and cellular responses to LD exposure, whereas DDR genes were prominently activated after HD. However, how and to what extent the observed dose-dependent differences in molecular and cellular responses may impact cancer development remain unclear, as the induction of chromosomal damage was found to be dose-proportional (10-200 mGy).

CDK12 regulates co-transcriptional splicing and RNA turnover in human cells.

The cyclin-dependent kinase CDK12 has garnered interest as a cancer therapeutic target as DNA damage response genes are particularly suppressed by loss of CDK12 activity. In this study, we assessed the acute effects of CDK12 inhibition on transcription and RNA processing using nascent RNA Bru-seq and BruChase-seq. Acute transcriptional changes were overall small after CDK12 inhibition but over 600 genes showed intragenic premature termination, including DNA repair and cell cycle genes. Furthermore, many genes showed reduced transcriptional readthrough past the end of genes in the absence of CDK12 activity. RNA turnover was dramatically affected by CDK12 inhibition and importantly, caused increased degradation of many transcripts from DNA damage response genes. We also show that co-transcriptional splicing was suppressed by CDK12 inhibition. Taken together, these studies reveal the roles of CDK12 in regulating transcription elongation, transcription termination, co-transcriptional splicing, and RNA turnover.

Myotubularin-related phosphatase 5 is a critical determinant of autophagy in neurons.

Autophagy is a conserved, multi-step process of capturing proteolytic cargo in autophagosomes for lysosome degradation. The capacity to remove toxic proteins that accumulate in neurodegenerative disorders attests to the disease-modifying potential of the autophagy pathway. However, neurons respond only marginally to conventional methods for inducing autophagy, limiting efforts to develop therapeutic autophagy modulators for neurodegenerative diseases. The determinants underlying poor autophagy induction in neurons and the degree to which neurons and other cell types are differentially sensitive to autophagy stimuli are incompletely defined. Accordingly, we sampled nascent transcript synthesis and stabilities in fibroblasts, induced pluripotent stem cells (iPSCs), and iPSC-derived neurons (iNeurons), thereby uncovering a neuron-specific stability of transcripts encoding myotubularin-related phosphatase 5 (MTMR5). MTMR5 is an autophagy suppressor that acts with its binding partner, MTMR2, to dephosphorylate phosphoinositides critical for autophagy initiation and autophagosome maturation. We found that MTMR5 is necessary and sufficient to suppress autophagy in iNeurons and undifferentiated iPSCs. Using optical pulse labeling to visualize the turnover of endogenously encoded proteins in live cells, we observed that knockdown of MTMR5 or MTMR2, but not the unrelated phosphatase MTMR9, significantly enhances neuronal degradation of TDP-43, an autophagy substrate implicated in several neurodegenerative diseases. Our findings thus establish a regulatory mechanism of autophagy intrinsic to neurons and targetable for clearing disease-related proteins in a cell-type-specific manner. In so doing, our results not only unravel novel aspects of neuronal biology and proteostasis but also elucidate a strategy for modulating neuronal autophagy that could be of high therapeutic potential for multiple neurodegenerative diseases.

KDM6A Regulates Cell Plasticity and Pancreatic Cancer Progression by Noncanonical Activin Pathway.

BACKGROUND & AIMS: Inactivating mutations of KDM6A, a histone demethylase, were frequently found in pancreatic ductal adenocarcinoma (PDAC). We investigated the role of KDM6A (lysine demethylase 6A) in PDAC development. METHODS: We performed a pancreatic tissue microarray analysis of KDM6A protein levels. We used human PDAC cell lines for KDM6A knockout and knockdown experiments. We performed bromouridine sequencing analysis to elucidate the effects of KDM6A loss on global transcription. We performed studies with Ptf1aCre; LSL-KrasG12D; Trp53R172H/+; Kdm6afl/fl or fl/Y, Ptf1aCre; Kdm6afl/fl or fl/Y, and orthotopic xenograft mice to investigate the impacts of Kdm6a deficiency on pancreatic tumorigenesis and pancreatitis. RESULTS: Loss of KDM6A was associated with metastasis in PDAC patients. Bromouridine sequencing analysis showed up-regulation of the epithelial-mesenchymal transition pathway in PDAC cells deficient in KDM6A. Loss of KDM6A promoted mesenchymal morphology, migration, and invasion in PDAC cells in vitro. Mechanistically, activin A and subsequent p38 activation likely mediated the role of KDM6A loss. Inhibiting either activin A or p38 reversed the effect. Pancreas-specific Kdm6a-knockout mice pancreata showed accelerated PDAC progression, developed a more aggressive undifferentiated type of PDAC, and increased metastases in the background of Kras and p53 mutations. Kdm6a-deficient pancreata in a pancreatitis model had a delayed recovery with increased PDAC precursor lesions compared with wild-type pancreata. CONCLUSIONS: Loss of KDM6A accelerates PDAC progression and metastasis, most likely by a noncanonical p38-dependent activin A pathway. KDM6A also promotes pancreatic tissue recovery from pancreatitis. Activin A might be used as a therapeutic target for KDM6A-deficient PDACs.

Dissecting regulatory pathways for transcription recovery following DNA damage reveals a non-canonical function of the histone chaperone HIRA.

Transcription restart after a genotoxic challenge is a fundamental yet poorly understood process. Here, we dissect the interplay between transcription and chromatin restoration after DNA damage by focusing on the human histone chaperone complex HIRA, which is required for transcription recovery post UV. We demonstrate that HIRA is recruited to UV-damaged chromatin via the ubiquitin-dependent segregase VCP to deposit new H3.3 histones. However, this local activity of HIRA is dispensable for transcription recovery. Instead, we reveal a genome-wide function of HIRA in transcription restart that is independent of new H3.3 and not restricted to UV-damaged loci. HIRA coordinates with ASF1B to control transcription restart by two independent pathways: by stabilising the associated subunit UBN2 and by reducing the expression of the transcription repressor ATF3. Thus, HIRA primes UV-damaged chromatin for transcription restart at least in part by relieving transcription inhibition rather than by depositing new H3.3 as an activating bookmark.

A CSB-PAF1C axis restores processive transcription elongation after DNA damage repair.

Bulky DNA lesions in transcribed strands block RNA polymerase II (RNAPII) elongation and induce a genome-wide transcriptional arrest. The transcription-coupled repair (TCR) pathway efficiently removes transcription-blocking DNA lesions, but how transcription is restored in the genome following DNA repair remains unresolved. Here, we find that the TCR-specific CSB protein loads the PAF1 complex (PAF1C) onto RNAPII in promoter-proximal regions in response to DNA damage. Although dispensable for TCR-mediated repair, PAF1C is essential for transcription recovery after UV irradiation. We find that PAF1C promotes RNAPII pause release in promoter-proximal regions and subsequently acts as a processivity factor that stimulates transcription elongation throughout genes. Our findings expose the molecular basis for a non-canonical PAF1C-dependent pathway that restores transcription throughout the human genome after genotoxic stress.

Publisher Correction: Regulation of cellular sterol homeostasis by the oxygen responsive noncoding RNA lincNORS.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Regulation of cellular sterol homeostasis by the oxygen responsive noncoding RNA lincNORS.

We hereby provide the initial portrait of lincNORS, a spliced lincRNA generated by the MIR193BHG locus, entirely distinct from the previously described miR-193b-365a tandem. While inducible by low O2 in a variety of cells and associated with hypoxia in vivo, our studies show that lincNORS is subject to multiple regulatory inputs, including estrogen signals. Biochemically, this lincRNA fine-tunes cellular sterol/steroid biosynthesis by repressing the expression of multiple pathway components. Mechanistically, the function of lincNORS requires the presence of RALY, an RNA-binding protein recently found to be implicated in cholesterol homeostasis. We also noticed the proximity between this locus and naturally occurring genetic variations highly significant for sterol/steroid-related phenotypes, in particular the age of sexual maturation. An integrative analysis of these variants provided a more formal link between these phenotypes and lincNORS, further strengthening the case for its biological relevance.

Discovery of Mitochondrial Transcription Inhibitors Active in Pancreatic Cancer Cells.

Mitochondrial dysfunction is a hallmark of cancer cells and targeting cancer mitochondria has emerged as a promising anti-cancer therapy. Previously, we repurposed chlorambucil by conjugating it to a mitochondrial targeting triphenylphosphonium (TPP) group to design Mito-Chlor, a novel agent that acts on mitochondria DNA (mtDNA). Herein, we show that Mito-Chlor, but not chlorambucil, inhibits the nascent transcription of mtDNA. Clustering analysis of transcriptomic profile of our Bru-seq database led to the identification of another mitochondrial transcription inhibitor SQD1, which inhibits the proliferation of MIA PaCa-2 cells with an IC50 of 1.3 µM. Interestingly, Mito-Chlor reduces expression of mitochondrial proteins, interferes with mitochondria membrane potential, and impairs oxidative phosphorylation while SQD1 does not. Both compounds increased cellular and mitochondrial reactive oxygen species and stimulated similar signaling pathways in response to oxidative stress. As mitochondrial transcription inhibitors and redox modulators, SQD1 and Mito-Chlor are promising for the treatment of pancreatic cancer by blocking mitochondrial function.

Principles of mRNA control by human PUM proteins elucidated from multimodal experiments and integrative data analysis.

The human PUF-family proteins, PUM1 and PUM2, posttranscriptionally regulate gene expression by binding to a PUM recognition element (PRE) in the 3'-UTR of target mRNAs. Hundreds of PUM1/2 targets have been identified from changes in steady-state RNA levels; however, prior studies could not differentiate between the contributions of changes in transcription and RNA decay rates. We applied metabolic labeling to measure changes in RNA turnover in response to depletion of PUM1/2, showing that human PUM proteins regulate expression almost exclusively by changing RNA stability. We also applied an in vitro selection workflow to precisely identify the binding preferences of PUM1 and PUM2. By integrating our results with prior knowledge, we developed a "rulebook" of key contextual features that differentiate functional versus nonfunctional PREs, allowing us to train machine learning models that accurately predict the functional regulation of RNA targets by the human PUM proteins.

Characterization of novel primary miRNA transcription units in human cells using Bru-seq nascent RNA sequencing.

MicroRNAs (miRNAs) are key contributors to gene regulatory networks. Because miRNAs are processed from RNA polymerase II transcripts, insight into miRNA regulation requires a comprehensive understanding of the regulation of primary miRNA transcripts. We used Bru-seq nascent RNA sequencing and hidden Markov model segmentation to map primary miRNA transcription units (TUs) across 32 human cell lines, allowing us to describe TUs encompassing 1443 miRNAs from miRBase and 438 from MirGeneDB. We identified TUs for 61 miRNAs with an unknown CAGE TSS signal for MirGeneDB miRNAs. Many primary transcripts containing miRNA sequences failed to generate mature miRNAs, suggesting that miRNA biosynthesis is under both transcriptional and post-transcriptional control. In addition to constitutive and cell-type specific TU expression regulated by differential promoter usage, miRNA synthesis can be regulated by transcription past polyadenylation sites (transcriptional read through) and promoter divergent transcription (PROMPTs). We identified 197 miRNA TUs with novel promoters, 97 with transcriptional read-throughs and 3 miRNA TUs that resemble PROMPTs in at least one cell line. The miRNA TU annotation data resource described here reveals a greater complexity in miRNA regulation than previously known and provides a framework for identifying cell-type specific differences in miRNA transcription in cancer and cell transition states.

Identifying cis Elements for Spatiotemporal Control of Mammalian DNA Replication.

The temporal order of DNA replication (replication timing [RT]) is highly coupled with genome architecture, but cis-elements regulating either remain elusive. We created a series of CRISPR-mediated deletions and inversions of a pluripotency-associated topologically associating domain (TAD) in mouse ESCs. CTCF-associated domain boundaries were dispensable for RT. CTCF protein depletion weakened most TAD boundaries but had no effect on RT or A/B compartmentalization genome-wide. By contrast, deletion of three intra-TAD CTCF-independent 3D contact sites caused a domain-wide early-to-late RT shift, an A-to-B compartment switch, weakening of TAD architecture, and loss of transcription. The dispensability of TAD boundaries and the necessity of these "early replication control elements" (ERCEs) was validated by deletions and inversions at additional domains. Our results demonstrate that discrete cis-regulatory elements orchestrate domain-wide RT, A/B compartmentalization, TAD architecture, and transcription, revealing fundamental principles linking genome structure and function.

Analysis of the androgen receptor-regulated lncRNA landscape identifies a role for ARLNC1 in prostate cancer progression.

The androgen receptor (AR) plays a critical role in the development of the normal prostate as well as prostate cancer. Using an integrative transcriptomic analysis of prostate cancer cell lines and tissues, we identified ARLNC1 (AR-regulated long noncoding RNA 1) as an important long noncoding RNA that is strongly associated with AR signaling in prostate cancer progression. Not only was ARLNC1 induced by the AR protein, but ARLNC1 stabilized the AR transcript via RNA-RNA interaction. ARLNC1 knockdown suppressed AR expression, global AR signaling and prostate cancer growth in vitro and in vivo. Taken together, these data support a role for ARLNC1 in maintaining a positive feedback loop that potentiates AR signaling during prostate cancer progression and identify ARLNC1 as a novel therapeutic target.

Transcriptional and post-transcriptional regulation of the ionizing radiation response by ATM and p53.

In response to ionizing radiation (IR), cells activate a DNA damage response (DDR) pathway to re-program gene expression. Previous studies using total cellular RNA analyses have shown that the stress kinase ATM and the transcription factor p53 are integral components required for induction of IR-induced gene expression. These studies did not distinguish between changes in RNA synthesis and RNA turnover and did not address the role of enhancer elements in DDR-mediated transcriptional regulation. To determine the contribution of synthesis and degradation of RNA and monitor the activity of enhancer elements following exposure to IR, we used the recently developed Bru-seq, BruChase-seq and BruUV-seq techniques. Our results show that ATM and p53 regulate both RNA synthesis and stability as well as enhancer element activity following exposure to IR. Importantly, many genes in the p53-signaling pathway were coordinately up-regulated by both increased synthesis and RNA stability while down-regulated genes were suppressed either by reduced synthesis or stability. Our study is the first of its kind that independently assessed the effects of ionizing radiation on transcription and post-transcriptional regulation in normal human cells.

Gene length as a biological timer to establish temporal transcriptional regulation.

Transcriptional timing is inherently influenced by gene length, thus providing a mechanism for temporal regulation of gene expression. While gene size has been shown to be important for the expression timing of specific genes during early development, whether it plays a role in the timing of other global gene expression programs has not been extensively explored. Here, we investigate the role of gene length during the early transcriptional response of human fibroblasts to serum stimulation. Using the nascent sequencing techniques Bru-seq and BruUV-seq, we identified immediate genome-wide transcriptional changes following serum stimulation that were linked to rapid activation of enhancer elements. We identified 873 significantly induced and 209 significantly repressed genes. Variations in gene size allowed for a large group of genes to be simultaneously activated but produce full-length RNAs at different times. The median length of the group of serum-induced genes was significantly larger than the median length of all expressed genes, housekeeping genes, and serum-repressed genes. These gene length relationships were also observed in corresponding mouse orthologs, suggesting that relative gene size is evolutionarily conserved. The sizes of transcription factor and microRNA genes immediately induced after serum stimulation varied dramatically, setting up a cascade mechanism for temporal expression arising from a single activation event. The retention and expansion of large intronic sequences during evolution have likely played important roles in fine-tuning the temporal expression of target genes in various cellular response programs.