m6A-driven SF3B1 translation control steers splicing to direct genome integrity and leukemogenesis.

SF3B1 is the most mutated splicing factor (SF) in myelodysplastic syndromes (MDSs), which are clonal hematopoietic disorders with variable risk of leukemic transformation. Although tumorigenic SF3B1 mutations have been extensively characterized, the role of "non-mutated" wild-type SF3B1 in cancer remains largely unresolved. Here, we identify a conserved epitranscriptomic program that steers SF3B1 levels to counteract leukemogenesis. Our analysis of human and murine pre-leukemic MDS cells reveals dynamic regulation of SF3B1 protein abundance, which affects MDS-to-leukemia progression in vivo. Mechanistically, ALKBH5-driven 5' UTR m6A demethylation fine-tunes SF3B1 translation directing splicing of central DNA repair and epigenetic regulators during transformation. This impacts genome stability and leukemia progression in vivo, supporting an integrative analysis in humans that SF3B1 molecular signatures may predict mutational variability and poor prognosis. These findings highlight a post-transcriptional gene expression nexus that unveils unanticipated SF3B1-dependent cancer vulnerabilities.

Probing the dynamic RNA structurome and its functions.

RNA is a key regulator of almost every cellular process, and the structures adopted by RNA molecules are thought to be central to their functions. The recent fast-paced evolution of high-throughput sequencing-based RNA structure mapping methods has enabled the rapid in vivo structural interrogation of entire cellular transcriptomes. Collectively, these studies are shedding new light on the long underestimated complexity of the structural organization of the transcriptome - the RNA structurome. Moreover, recent analyses are challenging the view that the RNA structurome is a static entity by revealing how RNA molecules establish intricate networks of alternative intramolecular and intermolecular interactions and that these ensembles of RNA structures are dynamically regulated to finely tune RNA functions in living cells. This new understanding of how RNA can shape cell phenotypes has important implications for the development of RNA-targeted therapeutic strategies.

Oncogenic translation directs spliceosome dynamics revealing an integral role for SF3A3 in breast cancer.

Splicing is a central RNA-based process commonly altered in human cancers; however, how spliceosomal components are co-opted during tumorigenesis remains poorly defined. Here we unravel the core splice factor SF3A3 at the nexus of a translation-based program that rewires splicing during malignant transformation. Upon MYC hyperactivation, SF3A3 levels are modulated translationally through an RNA stem-loop in an eIF3D-dependent manner. This ensures accurate splicing of mRNAs enriched for mitochondrial regulators. Altered SF3A3 translation leads to metabolic reprogramming and stem-like properties that fuel MYC tumorigenic potential in vivo. Our analysis reveals that SF3A3 protein levels predict molecular and phenotypic features of aggressive human breast cancers. These findings unveil a post-transcriptional interplay between splicing and translation that governs critical facets of MYC-driven oncogenesis.

Dnmt3L antagonizes DNA methylation at bivalent promoters and favors DNA methylation at gene bodies in ESCs.

The de novo DNA methyltransferase 3-like (Dnmt3L) is a catalytically inactive DNA methyltransferase that cooperates with Dnmt3a and Dnmt3b to methylate DNA. Dnmt3L is highly expressed in mouse embryonic stem cells (ESCs), but its function in these cells is unknown. Through genome-wide analysis of Dnmt3L knockdown in ESCs, we found that Dnmt3L is a positive regulator of methylation at the gene bodies of housekeeping genes and, more surprisingly, is also a negative regulator of methylation at promoters of bivalent genes. Dnmt3L is required for the differentiation of ESCs into primordial germ cells (PGCs) through the activation of the homeotic gene Rhox5. We demonstrate that Dnmt3L interacts with the Polycomb PRC2 complex in competition with the DNA methyltransferases Dnmt3a and Dnmt3b to maintain low methylation levels at the H3K27me3 regions. Thus, in ESCs, Dnmt3L counteracts the activity of de novo DNA methylases to maintain hypomethylation at promoters of bivalent developmental genes.

SHAPEwarp-web: sequence-agnostic search for structurally homologous RNA regions across databases of chemical probing data.

RNA molecules perform a variety of functions in cells, many of which rely on their secondary and tertiary structures. Chemical probing methods coupled with high-throughput sequencing have significantly accelerated the mapping of RNA structures, and increasingly large datasets of transcriptome-wide RNA chemical probing data are becoming available. Analogously to what has been done for decades in the protein world, this RNA structural information can be leveraged to aid the discovery of structural similarity to a known RNA (or RNA family), which, in turn, can inform about the function of transcripts. We have previously developed SHAPEwarp, a sequence-agnostic method for the search of structurally homologous RNA segments in a database of reactivity profiles derived from chemical probing experiments. In its original implementation, however, SHAPEwarp required substantial computational resources, even for moderately sized databases, as well as significant Linux command line know-how. To address these limitations, we introduce here SHAPEwarp-web, a user-friendly web interface to rapidly query large databases of RNA chemical probing data for structurally similar RNAs. Aside from featuring a completely rewritten core, which speeds up by orders of magnitude the search inside large databases, the web server hosts several high-quality chemical probing databases across multiple species. SHAPEwarp-web is available from https://shapewarp.incarnatolab.com.

Genome-scale deconvolution of RNA structure ensembles.

RNA structure heterogeneity is a major challenge when querying RNA structures with chemical probing. We introduce DRACO, an algorithm for the deconvolution of coexisting RNA conformations from mutational profiling experiments. Analysis of the SARS-CoV-2 genome using dimethyl sulfate mutational profiling with sequencing (DMS-MaPseq) and DRACO, identifies multiple regions that fold into two mutually exclusive conformations, including a conserved structural switch in the 3' untranslated region. This work may open the way to dissecting the heterogeneity of the RNA structurome.

The acetyltransferase p300 is recruited in trans to multiple enhancer sites by lncSmad7.

The histone acetyltransferase p300 (also known as KAT3B) is a general transcriptional coactivator that introduces the H3K27ac mark on enhancers triggering their activation and gene transcription. Genome-wide screenings demonstrated that a large fraction of long non-coding RNAs (lncRNAs) plays a role in cellular processes and organ development although the underlying molecular mechanisms remain largely unclear (1,2). We found 122 lncRNAs that interacts directly with p300. In depth analysis of one of these, lncSmad7, is required to maintain ESC self-renewal and it interacts to the C-terminal domain of p300. lncSmad7 also contains predicted RNA-DNA Hoogsteen forming base pairing. Combined Chromatin Isolation by RNA precipitation followed by sequencing (ChIRP-seq) together with CRISPR/Cas9 mutagenesis of the target sites demonstrate that lncSmad7 binds and recruits p300 to enhancers in trans, to trigger enhancer acetylation and transcriptional activation of its target genes. Thus, these results unveil a new mechanism by which p300 is recruited to the genome.

5-Ethynyluridine: A Bio-orthogonal Uridine Variant for mRNA-Based Therapies and Vaccines.

The identification of pseudo- and N1 -methylpseudo-uridine (Ψ and mΨ, respectively) as immunosilent uridine analogues has propelled the development of mRNA-based vaccines and therapeutics. Here, we have characterised another uridine analogue, 5-ethynyluridine (EU), which has an ethynyl moiety. We show that this uridine analogue does not cause immune activation in human macrophages, as it does not induce interleukin-6 secretion or expression of the inflammatory and antiviral genes MX1, PKR, and TAP2. Moreover, EU allows for prolonged expression, as shown with mRNA coding for yellow fluorescent protein (YFP). Side-by-side comparisons of EU with unmodified, Ψ, and mΨ revealed that EU-modified mRNA is expressed at lower levels, but confers similar stability and low immunogenicity to the other uridine analogues. Furthermore, structure analysis of modified mRNAs suggests that the observed phenotype is largely independent of RNA folding. Thus, EU is a potential candidate for RNA-based vaccines and therapeutics.

Intragenic DNA methylation prevents spurious transcription initiation.

In mammals, DNA methylation occurs mainly at CpG dinucleotides. Methylation of the promoter suppresses gene expression, but the functional role of gene-body DNA methylation in highly expressed genes has yet to be clarified. Here we show that, in mouse embryonic stem cells, Dnmt3b-dependent intragenic DNA methylation protects the gene body from spurious RNA polymerase II entry and cryptic transcription initiation. Using different genome-wide approaches, we demonstrate that this Dnmt3b function is dependent on its enzymatic activity and recruitment to the gene body by H3K36me3. Furthermore, the spurious transcripts can either be degraded by the RNA exosome complex or capped, polyadenylated, and delivered to the ribosome to produce aberrant proteins. Elongating RNA polymerase II therefore triggers an epigenetic crosstalk mechanism that involves SetD2, H3K36me3, Dnmt3b and DNA methylation to ensure the fidelity of gene transcription initiation, with implications for intragenic hypomethylation in cancer.

A novel SHAPE reagent enables the analysis of RNA structure in living cells with unprecedented accuracy.

Due to the mounting evidence that RNA structure plays a critical role in regulating almost any physiological as well as pathological process, being able to accurately define the folding of RNA molecules within living cells has become a crucial need. We introduce here 2-aminopyridine-3-carboxylic acid imidazolide (2A3), as a general probe for the interrogation of RNA structures in vivo. 2A3 shows moderate improvements with respect to the state-of-the-art selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) reagent NAI on naked RNA under in vitro conditions, but it significantly outperforms NAI when probing RNA structure in vivo, particularly in bacteria, underlining its increased ability to permeate biological membranes. When used as a restraint to drive RNA structure prediction, data derived by SHAPE-MaP with 2A3 yields more accurate predictions than NAI-derived data. Due to its extreme efficiency and accuracy, we can anticipate that 2A3 will rapidly take over conventional SHAPE reagents for probing RNA structures both in vitro and in vivo.

Hmga2 protein loss alters nuclear envelope and 3D chromatin structure.

BACKGROUND: The high-mobility group Hmga family of proteins are non-histone chromatin-interacting proteins which have been associated with a number of nuclear functions, including heterochromatin formation, replication, recombination, DNA repair, transcription, and formation of enhanceosomes. Due to its role based on dynamic interaction with chromatin, Hmga2 has a pathogenic role in diverse tumors and has been mainly studied in a cancer context; however, whether Hmga2 has similar physiological functions in normal cells remains less explored. Hmga2 was additionally shown to be required during the exit of embryonic stem cells (ESCs) from the ground state of pluripotency, to allow their transition into epiblast-like cells (EpiLCs), and here, we use that system to gain further understanding of normal Hmga2 function. RESULTS: We demonstrated that Hmga2 KO pluripotent stem cells fail to develop into EpiLCs. By using this experimental system, we studied the chromatin changes that take place upon the induction of EpiLCs and we observed that the loss of Hmga2 affects the histone mark H3K27me3, whose levels are higher in Hmga2 KO cells. Accordingly, a sustained expression of polycomb repressive complex 2 (PRC2), responsible for H3K27me3 deposition, was observed in KO cells. However, gene expression differences between differentiating wt vs Hmga2 KO cells did not show any significant enrichments of PRC2 targets. Similarly, endogenous Hmga2 association to chromatin in epiblast stem cells did not show any clear relationships with gene expression modification observed in Hmga2 KO. Hmga2 ChIP-seq confirmed that this protein preferentially binds to the chromatin regions associated with nuclear lamina. Starting from this observation, we demonstrated that nuclear lamina underwent severe alterations when Hmga2 KO or KD cells were induced to exit from the naïve state and this phenomenon is accompanied by a mislocalization of the heterochromatin mark H3K9me3 within the nucleus. As nuclear lamina (NL) is involved in the organization of 3D chromatin structure, we explored the possible effects of Hmga2 loss on this phenomenon. The analysis of Hi-C data in wt and Hmga2 KO cells allowed us to observe that inter-TAD (topologically associated domains) interactions in Hmga2 KO cells are different from those observed in wt cells. These differences clearly show a peculiar compartmentalization of inter-TAD interactions in chromatin regions associated or not to nuclear lamina. CONCLUSIONS: Overall, our results indicate that Hmga2 interacts with heterochromatic lamin-associated domains, and highlight a role for Hmga2 in the crosstalk between chromatin and nuclear lamina, affecting the establishment of inter-TAD interactions.

Multiplexed genome engineering by Cas12a and CRISPR arrays encoded on single transcripts.

The ability to modify multiple genetic elements simultaneously would help to elucidate and control the gene interactions and networks underlying complex cellular functions. However, current genome engineering technologies are limited in both the number and the type of perturbations that can be performed simultaneously. Here, we demonstrate that both Cas12a and a clustered regularly interspaced short palindromic repeat (CRISPR) array can be encoded in a single transcript by adding a stabilizer tertiary RNA structure. By leveraging this system, we illustrate constitutive, conditional, inducible, orthogonal and multiplexed genome engineering of endogenous targets using up to 25 individual CRISPR RNAs delivered on a single plasmid. Our method provides a powerful platform to investigate and orchestrate the sophisticated genetic programs underlying complex cell behaviors.

Genome-wide mapping of SARS-CoV-2 RNA structures identifies therapeutically-relevant elements.

SARS-CoV-2 is a betacoronavirus with a linear single-stranded, positive-sense RNA genome, whose outbreak caused the ongoing COVID-19 pandemic. The ability of coronaviruses to rapidly evolve, adapt, and cross species barriers makes the development of effective and durable therapeutic strategies a challenging and urgent need. As for other RNA viruses, genomic RNA structures are expected to play crucial roles in several steps of the coronavirus replication cycle. Despite this, only a handful of functionally-conserved coronavirus structural RNA elements have been identified to date. Here, we performed RNA structure probing to obtain single-base resolution secondary structure maps of the full SARS-CoV-2 coronavirus genome both in vitro and in living infected cells. Probing data recapitulate the previously described coronavirus RNA elements (5' UTR and s2m), and reveal new structures. Of these, ∼10.2% show significant covariation among SARS-CoV-2 and other coronaviruses, hinting at their functionally-conserved role. Secondary structure-restrained 3D modeling of these segments further allowed for the identification of putative druggable pockets. In addition, we identify a set of single-stranded segments in vivo, showing high sequence conservation, suitable for the development of antisense oligonucleotide therapeutics. Collectively, our work lays the foundation for the development of innovative RNA-targeted therapeutic strategies to fight SARS-related infections.

Molecular dynamics of targeting CD38 in multiple myeloma.

Multiple functions of CD38 need exploring to expand clinical application of anti-CD38 antibodies in multiple myeloma (MM). We investigated membrane dynamics of MM cells and subsequent events when CD38 is targeted by therapeutic antibodies. Human MM cells (BF01) were co-cultured in vitro with therapeutic antibody (or control immunoglobulin G) and analysed using gene expression profiling. Microvesicles from antibody-exposed cells were analysed for differential gene and microRNA (miRNA) expression, and for phenotypic characterisation. Exposure of BF01 cells to anti-CD38 antibody resulted in CD38 membrane redistribution, upregulation of metabolism-related genes and downregulation of genes involved in cell cycle processes. Microvesicles derived from antibody-exposed cells showed increased CD73 and CD39 expression, presence of programmed death-ligand 1 and significant up-/down-modulation of miRNAs. Microvesicles accumulated around immunoglobulin Fc receptor-positive (FcR+ ) cells. Upon internalisation, natural killer cells displayed significantly increased expression of genes related to activation and immune response, and downregulation of genes involved in the cell cycle. Cells may use microvesicles to transmit signals distally as part of a survival strategy. Microvesicles are equipped on their surface with enzymatic machinery leading to production of tolerogenic adenosine. Further, they are internalised in FcR+ cells with significant functional modifications. These observations have relevance for improving anti-CD38 therapeutic antibodies through targeting this mechanism and its sequelae.

Structure and regulation of coronavirus genomes: state-of-the-art and novel insights from SARS-CoV-2 studies.

Coronaviruses (CoV) are positive-sense single-stranded RNA viruses, harboring the largest viral RNA genomes known to date. Apart from the primary sequence encoding for all the viral proteins needed for the generation of new viral particles, certain regions of CoV genomes are known to fold into stable structures, controlling several aspects of CoV life cycle, from the regulation of the discontinuous transcription of subgenomic mRNAs, to the packaging of the genome into new virions. Here we review the current knowledge on CoV RNA structures, discussing it in light of the most recent discoveries made possible by analyses of the SARS-CoV-2 genome.

In vivo analysis of influenza A mRNA secondary structures identifies critical regulatory motifs.

The influenza A virus (IAV) is a continuous health threat to humans as well as animals due to its recurring epidemics and pandemics. The IAV genome is segmented and the eight negative-sense viral RNAs (vRNAs) are transcribed into positive sense complementary RNAs (cRNAs) and viral messenger RNAs (mRNAs) inside infected host cells. A role for the secondary structure of IAV mRNAs has been hypothesized and debated for many years, but knowledge on the structure mRNAs adopt in vivo is currently missing. Here we solve, for the first time, the in vivo secondary structure of IAV mRNAs in living infected cells. We demonstrate that, compared to the in vitro refolded structure, in vivo IAV mRNAs are less structured but exhibit specific locally stable elements. Moreover, we show that the targeted disruption of these high-confidence structured domains results in an extraordinary attenuation of IAV replicative capacity. Collectively, our data provide the first comprehensive map of the in vivo structural landscape of IAV mRNAs, hence providing the means for the development of new RNA-targeted antivirals.

RNA Framework: an all-in-one toolkit for the analysis of RNA structures and post-transcriptional modifications.

RNA is emerging as a key regulator of a plethora of biological processes. While its study has remained elusive for decades, the recent advent of high-throughput sequencing technologies provided the unique opportunity to develop novel techniques for the study of RNA structure and post-transcriptional modifications. Nonetheless, most of the required downstream bioinformatics analyses steps are not easily reproducible, thus making the application of these techniques a prerogative of few laboratories. Here we introduce RNA Framework, an all-in-one toolkit for the analysis of most NGS-based RNA structure probing and post-transcriptional modification mapping experiments. To prove the extreme versatility of RNA Framework, we applied it to both an in-house generated DMS-MaPseq dataset, and to a series of literature available experiments. Notably, when starting from publicly available datasets, our software easily allows replicating authors' findings. Collectively, RNA Framework provides the most complete and versatile toolkit to date for a rapid and streamlined analysis of the RNA epistructurome. RNA Framework is available for download at: http://www.rnaframework.com.

In vivo probing of nascent RNA structures reveals principles of cotranscriptional folding.

Defining the in vivo folding pathway of cellular RNAs is essential to understand how they reach their final native conformation. We here introduce a novel method, named Structural Probing of Elongating Transcripts (SPET-seq), that permits single-base resolution analysis of transcription intermediates' secondary structures on a transcriptome-wide scale, enabling base-resolution analysis of the RNA folding events. Our results suggest that cotranscriptional RNA folding in vivo is a mixture of cooperative folding events, in which local RNA secondary structure elements are formed as they get transcribed, and non-cooperative events, in which 5'-halves of long-range helices get sequestered into transient non-native interactions until their 3' counterparts have been transcribed. Together our work provides the first transcriptome-scale overview of RNA cotranscriptional folding in a living organism.

High-throughput single-base resolution mapping of RNA 2΄-O-methylated residues.

Functional characterization of the transcriptome requires tools for the systematic investigation of RNA post-transcriptional modifications. 2΄-O-methylation (2΄-OMe) of the ribose moiety is one of the most abundant post-transcriptional modifications of RNA, although its systematic analysis is difficult due to the lack of reliable high-throughput mapping methods. We describe here a novel high-throughput approach, named 2OMe-seq, that enables fast and accurate mapping at single-base resolution, and relative quantitation, of 2΄-OMe modified residues. We compare our method to other state-of-art approaches, and show that it achieves higher sensitivity and specificity. By applying 2OMe-seq to HeLa cells, we show that it is able to recover the majority of the annotated 2΄-OMe sites on ribosomal RNA. By performing knockdown of the Fibrillarin methyltransferase in mouse embryonic stem cells (ESCs) we show the ability of 2OMe-seq to capture 2΄-O-Methylation level variations. Moreover, using 2OMe-seq data we here report the discovery of 12 previously unannotated 2΄-OMe sites across 18S and 28S rRNAs, 11 of which are conserved in both human and mouse cells, and assigned the respective snoRNAs for all sites. Our approach expands the repertoire of methods for transcriptome-wide mapping of RNA post-transcriptional modifications, and promises to provide novel insights into the role of this modification.

The long intergenic non-coding RNA CCR492 functions as a let-7 competitive endogenous RNA to regulate c-Myc expression.

In mammals the cell-cycle progression through the G1 phase is a tightly regulated process mediated by the transcriptional activation of early genes in response to mitogenic stimuli, whose dysregulation often leads to tumorigenesis. We here report the discovery by RNA-seq of cell-cycle regulated (CCR) long intergenic non-coding RNAs (lincRNAs), potentially involved in the control of the cell-cycle progression. We identified 10 novel lincRNAs expressed in response to serum treatment in mouse embryonic fibroblasts (MEFs) and in BALB/c fibroblasts, comparably to early genes. By loss-of-function experiments we found that lincRNA CCR492 is required for G1/S progression, localizes in the cell cytoplasm and contains 4 let-7 microRNA recognition elements (MREs). Mechanistically, CCR492 functions as a competing endogenous RNA (ceRNA) to antagonize the function of let-7 microRNAs, leading to the de-repression of c-Myc. Moreover, we show that ectopic expression of CCR492 along with a constitutively active H-Ras promotes cell transformation. Thus, we identified a new lincRNA expressed as an early gene in mammalian cells to regulate the cell-cycle progression by upregulating c-Myc expression.