Dynamic imaging of nascent RNA reveals general principles of transcription dynamics and stochastic splice site selection.

The activities of RNA polymerase and the spliceosome are responsible for the heterogeneity in the abundance and isoform composition of mRNA in human cells. However, the dynamics of these megadalton enzymatic complexes working in concert on endogenous genes have not been described. Here, we establish a quasi-genome-scale platform for observing synthesis and processing kinetics of single nascent RNA molecules in real time. We find that all observed genes show transcriptional bursting. We also observe large kinetic variation in intron removal for single introns in single cells, which is inconsistent with deterministic splice site selection. Transcriptome-wide footprinting of the U2AF complex, nascent RNA profiling, long-read sequencing, and lariat sequencing further reveal widespread stochastic recursive splicing within introns. We propose and validate a unified theoretical model to explain the general features of transcription and pervasive stochastic splice site selection.

The p53-induced RNA-binding protein ZMAT3 is a splicing regulator that inhibits the splicing of oncogenic CD44 variants in colorectal carcinoma.

p53 is an intensely studied tumor-suppressive transcription factor. Recent studies suggest that the RNA-binding protein (RBP) ZMAT3 is important in mediating the tumor-suppressive effects of p53. Here, we globally identify ZMAT3-regulated RNAs and their binding sites at nucleotide resolution in intact colorectal cancer (CRC) cells. ZMAT3 binds to thousands of mRNA precursors, mainly at intronic uridine-rich sequences and affects their splicing. The strongest alternatively spliced ZMAT3 target was CD44, a cell adhesion gene and stem cell marker that controls tumorigenesis. Silencing ZMAT3 increased inclusion of CD44 variant exons, resulting in significant up-regulation of oncogenic CD44 isoforms (CD44v) and increased CRC cell growth that was rescued by concurrent knockdown of CD44v Silencing p53 phenocopied the loss of ZMAT3 with respect to CD44 alternative splicing, suggesting that ZMAT3-mediated regulation of CD44 splicing is vital for p53 function. Collectively, our findings uncover a p53-ZMAT3-CD44 axis in growth suppression in CRC cells.

The splicing factor U2AF1 contributes to cancer progression through a noncanonical role in translation regulation.

Somatic mutations in the genes encoding components of the spliceosome occur frequently in human neoplasms, including myeloid dysplasias and leukemias, and less often in solid tumors. One of the affected factors, U2AF1, is involved in splice site selection, and the most common change, S34F, alters a conserved nucleic acid-binding domain, recognition of the 3' splice site, and alternative splicing of many mRNAs. However, the role that this mutation plays in oncogenesis is still unknown. Here, we uncovered a noncanonical function of U2AF1, showing that it directly binds mature mRNA in the cytoplasm and negatively regulates mRNA translation. This splicing-independent role of U2AF1 is altered by the S34F mutation, and polysome profiling indicates that the mutation affects translation of hundreds of mRNA. One functional consequence is increased synthesis of the secreted chemokine interleukin 8, which contributes to metastasis, inflammation, and cancer progression in mice and humans.

Argonaute-miRNA Complexes Silence Target mRNAs in the Nucleus of Mammalian Stem Cells.

In mammals, gene silencing by the RNA-induced silencing complex (RISC) is a well-understood cytoplasmic posttranscriptional gene regulatory mechanism. Here, we show that embryonic stem cells (ESCs) contain high levels of nuclear AGO proteins and that in ESCs nuclear AGO protein activity allows for the onset of differentiation. In the nucleus, AGO proteins interact with core RISC components, including the TNRC6 proteins and the CCR4-NOT deadenylase complex. In contrast to cytoplasmic miRNA-mediated gene silencing that mainly operates on cis-acting elements in mRNA 3' untranslated (UTR) sequences, in the nucleus AGO binding in the coding sequence and potentially introns also contributed to post-transcriptional gene silencing. Thus, nuclear localization of AGO proteins in specific cell types leads to a previously unappreciated expansion of the miRNA-regulated transcriptome.

A Muscle-Specific Enhancer RNA Mediates Cohesin Recruitment and Regulates Transcription In trans.

The enhancer regions of the myogenic master regulator MyoD give rise to at least two enhancer RNAs. Core enhancer eRNA (CEeRNA) regulates transcription of the adjacent MyoD gene, whereas DRReRNA affects expression of Myogenin in trans. We found that DRReRNA is recruited at the Myogenin locus, where it colocalizes with Myogenin nascent transcripts. DRReRNA associates with the cohesin complex, and this association correlates with its transactivating properties. Despite being expressed in undifferentiated cells, cohesin is not loaded on Myogenin until the cells start expressing DRReRNA, which is then required for cohesin chromatin recruitment and maintenance. Functionally, depletion of either cohesin or DRReRNA reduces chromatin accessibility, prevents Myogenin activation, and hinders muscle cell differentiation. Thus, DRReRNA ensures spatially appropriate cohesin loading in trans to regulate gene expression.

A non-radioactive, improved PAR-CLIP and small RNA cDNA library preparation protocol.

Crosslinking and immunoprecipitation (CLIP) methods are powerful techniques to interrogate direct protein-RNA interactions and dissect posttranscriptional gene regulatory networks. One widely used CLIP variant is photoactivatable ribonucleoside enhanced CLIP (PAR-CLIP) that involves in vivo labeling of nascent RNAs with the photoreactive nucleosides 4-thiouridine (4SU) or 6-thioguanosine (6SG), which can efficiently crosslink to interacting proteins using UVA and UVB light. Crosslinking of 4SU or 6SG to interacting amino acids changes their base-pairing properties and results in characteristic mutations in cDNA libraries prepared for high-throughput sequencing, which can be computationally exploited to remove abundant background from non-crosslinked sequences and help pinpoint RNA binding protein binding sites at nucleotide resolution on a transcriptome-wide scale. Here we present a streamlined protocol for fluorescence-based PAR-CLIP (fPAR-CLIP) that eliminates the need to use radioactivity. It is based on direct ligation of a fluorescently labeled adapter to the 3'end of crosslinked RNA on immobilized ribonucleoproteins, followed by isolation of the adapter-ligated RNA and efficient conversion into cDNA without the previously needed size fractionation on denaturing polyacrylamide gels. These improvements cut the experimentation by half to 2 days and increases sensitivity by 10-100-fold.

Proximity-CLIP provides a snapshot of protein-occupied RNA elements in subcellular compartments.

Methods for the systematic study of subcellular RNA localization are limited, and their development has lagged behind that of proteomic tools. We combined APEX2-mediated proximity biotinylation of proteins with photoactivatable ribonucleoside-enhanced crosslinking to simultaneously profile the proteome and the transcriptome bound by RNA-binding proteins in any given subcellular compartment. Our approach is fractionation independent and allows study of the localization of RNA processing intermediates, as well as the identification of regulatory RNA cis-acting elements occupied by proteins, in a cellular-compartment-specific manner. We used our method, Proximity-CLIP, to profile RNA and protein in the nucleus, in the cytoplasm, and at cell-cell interfaces. Among other insights, we observed frequent transcriptional readthrough continuing for several kilobases downstream of the canonical cleavage and polyadenylation site and a differential RBP occupancy pattern for mRNAs in the nucleus and cytoplasm. We observed that mRNAs localized to cell-cell interfaces often encoded regulatory proteins and contained protein-occupied CUG sequence elements in their 3' untranslated region.

Pharmacoepigenomics circuits induced by a novel retinoid-polyamine conjugate in human immortalized keratinocytes.

Retinoids are widely used in diseases spanning from dermatological lesions to cancer, but exhibit severe adverse effects. A novel all-trans-Retinoic Acid (atRA)-spermine conjugate (termed RASP) has shown previously optimal in vitro and in vivo anti-inflammatory and anticancer efficacy, with undetectable teratogenic and toxic side-effects. To get insights, we treated HaCaT cells which resemble human epidermis with IC50 concentration of RASP and analyzed their miRNA expression profile. Gene ontology analysis of their predicted targets indicated dynamic networks involved in cell proliferation, signal transduction and apoptosis. Furthermore, DNA microarrays analysis verified that RASP affects the expression of the same categories of genes. A protein-protein interaction map produced using the most significant common genes, revealed hub genes of nodal functions. We conclude that RASP is a synthetic retinoid derivative with improved properties, which possess the beneficial effects of retinoids without exhibiting side-effects and with potential beneficial effects against skin diseases including skin cancer.

Mammalian PNLDC1 is a novel poly(A) specific exonuclease with discrete expression during early development.

PNLDC1 is a homologue of poly(A) specific ribonuclease (PARN), a known deadenylase with additional role in processing of non-coding RNAs. Both enzymes were reported recently to participate in piRNA biogenesis in silkworm and C. elegans, respectively. To get insights on the role of mammalian PNLDC1, we characterized the human and mouse enzymes. PNLDC1 shows limited conservation compared to PARN and represents an evolutionary related but distinct group of enzymes. It is expressed specifically in mouse embryonic stem cells, human and mouse testes and during early mouse embryo development, while it fades during differentiation. Its expression in differentiated cells, is suppressed through methylation of its promoter by the de novo methyltransferase DNMT3B. Both enzymes are localized mainly in the ER and exhibit in vitro specificity restricted solely to 3' RNA or DNA polyadenylates. Knockdown of Pnldc1 in mESCs and subsequent NGS analysis showed that although the expression of the remaining deadenylases remains unaffected, it affects genes involved mainly in reprogramming, cell cycle and translational regulation. Mammalian PNLDC1 is a novel deadenylase expressed specifically in cell types which share regulatory mechanisms required for multipotency maintenance. Moreover, it could be involved both in posttranscriptional regulation through deadenylation and genome surveillance during early development.

Dysregulation of Lipid Metabolism in Mkp-1 Deficient Mice during Gram-Negative Sepsis.

Mitogen-activated protein kinase phosphatase (Mkp)-1 exerts its anti-inflammatory activities during Gram-negative sepsis by deactivating p38 and c-Jun N-terminal kinase (JNK). We have previously shown that Mkp-1+/+ mice, but not Mkp-1-/- mice, exhibit hypertriglyceridemia during severe sepsis. However, the regulation of hepatic lipid stores and the underlying mechanism of lipid dysregulation during sepsis remains an enigma. To understand the molecular mechanism underlying the sepsis-associated metabolic changes and the role of Mkp-1 in the process, we infected Mkp-1+/+ and Mkp-1-/- mice with Escherichia coli i.v., and assessed the effects of Mkp-1 deficiency on tissue lipid contents. We also examined the global gene expression profile in the livers via RNA-seq. We found that in the absence of E. coli infection, Mkp-1 deficiency decreased liver triglyceride levels. Upon E. coli infection, Mkp-1+/+ mice, but not Mkp-1-/- mice, developed hepatocyte ballooning and increased lipid deposition in the livers. E. coli infection caused profound changes in the gene expression profile of a large number of proteins that regulate lipid metabolism in wildtype mice, while these changes were substantially disrupted in Mkp-1-/- mice. Interestingly, in Mkp-1+/+ mice E. coli infection resulted in downregulation of genes that facilitate fatty acid synthesis but upregulation of Cd36 and Dgat2, whose protein products mediate fatty acid uptake and triglyceride synthesis, respectively. Taken together, our studies indicate that sepsis leads to a substantial change in triglyceride metabolic gene expression programs and Mkp-1 plays an important role in this process.

Epigenetic repression of antiviral genes by SARS-CoV-2 NSP1.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) evades the innate immune machinery through multiple viral proteins, including nonstructural protein 1 (NSP1). While NSP1 is known to suppress translation of host mRNAs, the mechanisms underlying its immune evasion properties remain elusive. By integrating RNA-seq, ribosome footprinting, and ChIP-seq in A549 cells we found that NSP1 predominantly represses transcription of immune-related genes by favoring Histone 3 Lysine 9 dimethylation (H3K9me2). G9a/GLP H3K9 methyltransferase inhibitor UNC0638 restored expression of antiviral genes and restricted SARS-CoV-2 replication. Our multi-omics study unravels an epigenetic mechanism underlying host immune evasion by SARS-CoV-2 NSP1. Elucidating the factors involved in this phenomenon, may have implications for understanding and treating viral infections and other immunomodulatory diseases.

Structural basis of MALAT1 RNA maturation and mascRNA biogenesis.

The metastasis-associated lung adenocarcinoma transcript 1 (MALAT1) long noncoding RNA (lncRNA) has key roles in regulating transcription, splicing, tumorigenesis, etc. Its maturation and stabilization require precise processing by RNase P, which simultaneously initiates the biogenesis of a 3' cytoplasmic MALAT1-associated small cytoplasmic RNA (mascRNA). mascRNA was proposed to fold into a transfer RNA (tRNA)-like secondary structure but lacks eight conserved linking residues required by the canonical tRNA fold. Here we report crystal structures of human mascRNA before and after processing, which reveal an ultracompact, quasi-tRNA-like structure. Despite lacking all linker residues, mascRNA faithfully recreates the characteristic 'elbow' feature of tRNAs to recruit RNase P and ElaC homolog protein 2 (ELAC2) for processing, which exhibit distinct substrate specificities. Rotation and repositioning of the D-stem and anticodon regions preclude mascRNA from aminoacylation, avoiding interference with translation. Therefore, a class of metazoan lncRNA loci uses a previously unrecognized, unusually streamlined quasi-tRNA architecture to recruit select tRNA-processing enzymes while excluding others to drive bespoke RNA biogenesis, processing and maturation.

Integrative transcriptomic and proteomic profiling of the effects of cell confluency on gene expression.

In this study we examine the impact of cell confluency on gene expression. We focused on Argonaute (AGO) protein dynamics and associated gene and protein expression in HEK293, A375, and SHSY5Y cell lines. As a consequence of cell confluency, AGO2 protein translocates into the nucleus. Therefore, we generated transcriptomic data using RNA sequencing to compare gene expression in subconfluent versus confluent cells, which highlighted significant alterations in gene regulation patterns directly corresponding to changes in cell density. Our study also encompasses miRNA profiling data obtained through small RNA sequencing, revealing miRNA expressional changes dependent on cellular confluency, as well as cellular localization. Finally, we derived proteomic data from mass spectrometry analyses following AGO1-4 immunoprecipitation, providing a comprehensive view of AGO interactome in both nuclear and cytoplasmic compartments under varying confluency. These datasets offer a detailed exploration of the cellular and molecular dynamics, influenced by cell confluency, presenting a valuable resource for further research in cellular biology, particularly in understanding the basic mechanisms of cell density in cancer cells.

Keratin 19 binds and regulates cytoplasmic HNRNPK mRNA targets in triple-negative breast cancer.

BACKGROUND: Heterogeneous nuclear ribonucleoprotein K (HNRNPK) regulates pre-mRNA processing and long non-coding RNA localization in the nucleus. It was previously shown that shuttling of HNRNPK to the cytoplasm promotes cell proliferation and cancer metastasis. However, the mechanism of HNRNPK cytoplasmic localization, its cytoplasmic RNA ligands, and impact on post-transcriptional gene regulation remain uncharacterized. RESULTS: Here we show that the intermediate filament protein Keratin 19 (K19) directly interacts with HNRNPK and sequesters it in the cytoplasm. Correspondingly, in K19 knockout breast cancer cells, HNRNPK does not localize in the cytoplasm, resulting in reduced cell proliferation. We comprehensively mapped HNRNPK binding sites on mRNAs and showed that, in the cytoplasm, K19-mediated HNRNPK-retention increases the abundance of target mRNAs bound to the 3' untranslated region (3'UTR) at the expected cytidine-rich (C-rich) sequence elements. Furthermore, these mRNAs protected by HNRNPK in the cytoplasm are typically involved in cancer progression and include the p53 signaling pathway that is dysregulated upon HNRNPK knockdown (HNRNPK KD) or K19 knockout (KRT19 KO). CONCLUSIONS: This study identifies how a cytoskeletal protein can directly regulate gene expression by controlling the subcellular localization of RNA-binding proteins to support pathways involved in cancer progression.

Matrin3 regulates mitotic spindle dynamics by controlling alternative splicing of CDC14B.

Matrin3 is an RNA-binding protein that regulates diverse RNA-related processes, including mRNA splicing. Although Matrin3 has been intensively studied in neurodegenerative diseases, its function in cancer remains unclear. Here, we report Matrin3-mediated regulation of mitotic spindle dynamics in colorectal cancer (CRC) cells. We comprehensively identified RNAs bound and regulated by Matrin3 in CRC cells and focused on CDC14B, one of the top Matrin3 targets. Matrin3 knockdown results in increased inclusion of an exon containing a premature termination codon in the CDC14B transcript and simultaneous down-regulation of the standard CDC14B transcript. Knockdown of CDC14B phenocopies the defects in mitotic spindle dynamics upon Matrin3 knockdown, and the elongated and misoriented mitotic spindle observed upon Matrin3 knockdown are rescued upon overexpression of CDC14B, suggesting that CDC14B is a key downstream effector of Matrin3. Collectively, these data reveal a role for the Matrin3/CDC14B axis in control of mitotic spindle dynamics.

Cytokine-enhanced cytolytic activity of exosomes from NK Cells.

Natural killer (NK) cells play key roles in immune surveillance against tumors and viral infection. NK cells distinguish abnormal cells from healthy cells by cell-cell interaction with cell surface proteins and then attack target cells via multiple mechanisms. In addition, extracellular vesicles (EVs) derived from NK cells (NK-EVs), including exosomes, possess cytotoxic capacity against tumor cells, but their characteristics and regulation by cytokines remain unknown. Here, we report that EVs derived from human NK-92 cells stimulated with IL-15 + IL-21 show enhanced cytotoxic capacity against tumor cells. Major cytolytic granules, granzyme B and granzyme H, are enriched by IL-15 + IL-21 stimulation in NK-EVs; however, knockout experiments reveal those cytolytic granules are independent of enhanced cytotoxic capacity. To find out the key molecules, mass spectrometry analyses were performed with different cytokine conditions, no cytokine, IL-15, IL-21, or IL-15 + IL-21. We then found that CD226 (DNAM-1) on NK-EVs is enriched by IL-15 + IL-21 stimulation and that blocking antibodies against CD226 reduced the cytolytic activity of NK-EVs. We also show NK-EVs are taken up by target cells via macropinocytosis. Collectively, our findings elucidate the novel properties of NK-EVs and the mechanism of their incorporation into target cells.

Proximity-CLIP and Expedited Non-Radioactive Library Preparation of Small RNA Footprints for Next-Generation Sequencing.

During the course of their life cycle, most RNAs move between several cellular environments where they associate with different RNA binding proteins (RBPs). Reciprocally, a significant portion of RBPs reside in more than a single cellular compartment, where they can interact with discrete RNAs and even exert distinct biological roles. Proximity-CLIP combines proximity biotinylation of proteins with photoactivatable ribonucleoside-enhanced protein-RNA crosslinking to simultaneously profile the proteome, including RBPs and the RBP-bound transcriptome, in any given subcellular compartment. Here we provide a detailed experimental protocol for Proximity-CLIP along with a simplified non-radioactive, small-RNA cDNA library preparation protocol. Published 2020 U.S. Government. Basic Protocol 1: Cell culture, 4SU labeling, proximity biotinylation, and crosslinking Basic Protocol 2: Cell extraction, streptavidin affinity purification, and on-beads trypsinization Basic Protocol 3: RNA footprints cDNA library preparation Support Protocol: Preparation of RNA-seq libraries from intact RNA.

Multiple capsid protein binding sites mediate selective packaging of the alphavirus genomic RNA.

The alphavirus capsid protein (Cp) selectively packages genomic RNA (gRNA) into the viral nucleocapsid to produce infectious virus. Using photoactivatable ribonucleoside crosslinking and an innovative biotinylated Cp retrieval method, here we comprehensively define binding sites for Semliki Forest virus (SFV) Cp on the gRNA. While data in infected cells demonstrate Cp binding to the proposed genome packaging signal (PS), mutagenesis experiments show that PS is not required for production of infectious SFV or Chikungunya virus. Instead, we identify multiple Cp binding sites that are enriched on gRNA-specific regions and promote infectious SFV production and gRNA packaging. Comparisons of binding sites in cytoplasmic vs. viral nucleocapsids demonstrate that budding causes discrete changes in Cp-gRNA interactions. Notably, Cp's top binding site is maintained throughout virus assembly, and specifically binds and assembles with Cp into core-like particles in vitro. Together our data suggest a model for selective alphavirus genome recognition and assembly.

Contribution of miRNAs, tRNAs and tRFs to Aberrant Signaling and Translation Deregulation in Lung Cancer.

Transcriptomics profiles of miRNAs, tRNAs or tRFs are used as biomarkers, after separate examination of several cancer cell lines, blood samples or biopsies. However, the possible contribution of all three profiles on oncogenic signaling and translation as a net regulatory effect, is under investigation. The present analysis of miRNAs and tRFs from lung cancer biopsies indicated putative targets, which belong to gene networks involved in cell proliferation, transcription and translation regulation. In addition, we observed differential expression of specific tRNAs along with several tRNA-related genes with possible involvement in carcinogenesis. Transfection of lung adenocarcinoma cells with two identified tRFs and subsequent NGS analysis indicated gene targets that mediate signaling and translation regulation. Broader analysis of all major signaling and translation factors in several biopsy specimens revealed a crosstalk between the PI3K/AKT and MAPK pathways and downstream activation of eIF4E and eEF2. Subsequent polysome profile analysis and 48S pre-initiation reconstitution experiments showed increased global translation rates and indicated that aberrant expression patterns of translation initiation factors could contribute to elevated protein synthesis. Overall, our results outline the modulatory effects that possibly correlate the expression of important regulatory non-coding RNAs with aberrant signaling and translation deregulation in lung cancer.

Competitive inhibition of human poly(A)-specific ribonuclease (PARN) by synthetic fluoro-pyranosyl nucleosides.

Poly(A)-specific ribonuclease (PARN) is a cap-interacting deadenylase that mediates, together with other exonucleases, the eukaryotic mRNA turnover and thus is actively involved in the regulation of gene expression. Aminoglycosides and natural nucleotides are the only reported modulators of human PARN activity, so far. In the present study, we show that synthetic nucleoside analogues bearing a fluoro-glucopyranosyl sugar moiety and benzoyl-modified cytosine or adenine as a base can effectively inhibit human PARN. Such nucleoside analogues exhibited substantial inhibitory effects, when tested against various cancer cell lines, as has been previously reported. Kinetic analysis showed that the inhibition of PARN is competitive and could not be released by altering Mg(II) concentration. Moreover, substitution of the 2', 4', or 6'-OH of the sugar moiety with acetyl and/or trityl groups was crucial for inhibitory efficacy. To understand how the nucleosides fit into the active site of PARN, we performed molecular docking experiments followed by molecular dynamics simulations. The in silico analysis showed that these compounds can efficiently dock into the active site of PARN. Our results support the idea that the sugar moiety mediates the stabilization of the nucleoside into the active site through interactions with catalytic amino acid residues. Taken together, our in vitro and in silico data suggest that human PARN is among the molecular targets of these compounds and could act therapeutically by lowering the mRNA turnover rate, thus explaining their known in vivo inhibitory effect at the molecular level.