RNA molecular recording with an engineered RNA deaminase.

RNA deaminases are powerful tools for base editing and RNA molecular recording. However, the enzymes used in currently available RNA molecular recorders such as TRIBE, DART or STAMP have limitations due to RNA structure and sequence dependence. We designed a platform for directed evolution of RNA molecular recorders. We engineered an RNA A-to-I deaminase (an RNA adenosine base editor, rABE) that has high activity, low bias and low background. Using rABE, we present REMORA (RNA-encoded molecular recording in adenosines), wherein deamination by rABE writes a molecular record of RNA-protein interactions. By combining rABE with the C-to-U deaminase APOBEC1 and long-read RNA sequencing, we measured binding by two RNA-binding proteins on single messenger RNAs. Orthogonal RNA molecular recording of mammalian Pumilio proteins PUM1 and PUM2 shows that PUM1 competes with PUM2 for a subset of sites in cells. Furthermore, we identify transcript isoform-specific RNA-protein interactions driven by isoform changes distal to the binding site. The genetically encodable RNA deaminase rABE enables single-molecule identification of RNA-protein interactions with cell type specificity.

Impacts of uORF codon identity and position on translation regulation.

Translation regulation plays an important role in eukaryotic gene expression. Upstream open reading frames (uORFs) are potent regulatory elements located in 5' mRNA transcript leaders. Translation of uORFs usually inhibit the translation of downstream main open reading frames, but some enhance expression. While a minority of uORFs encode conserved functional peptides, the coding regions of most uORFs are not conserved. Thus, the importance of uORF coding sequences on their regulatory functions remains largely unknown. We investigated the impact of an uORF coding region on gene regulation by assaying the functions of thousands of variants in the yeast YAP1 uORF. Varying uORF codons resulted in a wide range of functions, including repressing and enhancing expression of the downstream ORF. The presence of rare codons resulted in the most inhibitory YAP1 uORF variants. Inhibitory functions of such uORFs were abrogated by overexpression of complementary tRNA. Finally, regression analysis of our results indicated that both codon identity and position impact uORF function. Our results support a model in which a uORF coding sequence impacts its regulatory functions by altering the speed of uORF translation.

Structural analyses of NEAT1 lncRNAs suggest long-range RNA interactions that may contribute to paraspeckle architecture.

Paraspeckles are nuclear bodies that regulate multiple aspects of gene expression. The long non-coding RNA (lncRNA) NEAT1 is essential for paraspeckle formation. NEAT1 has a highly ordered spatial organization within the paraspeckle, such that its 5' and 3' ends localize on the periphery of paraspeckle, while central sequences of NEAT1 are found within the paraspeckle core. As such, the structure of NEAT1 RNA may be important as a scaffold for the paraspeckle. In this study, we used SHAPE probing and computational analyses to investigate the secondary structure of human and mouse NEAT1. We propose a secondary structural model of the shorter (3,735 nt) isoform hNEAT1_S, in which the RNA folds into four separate domains. The secondary structures of mouse and human NEAT1 are largely different, with the exception of several short regions that have high structural similarity. Long-range base-pairing interactions between the 5' and 3' ends of the long isoform NEAT1 (NEAT1_L) were predicted computationally and verified using an in vitro RNA-RNA interaction assay. These results suggest that the conserved role of NEAT1 as a paraspeckle scaffold does not require extensively conserved RNA secondary structure and that long-range interactions among NEAT1 transcripts may have an important architectural function in paraspeckle formation.

Mod-seq: high-throughput sequencing for chemical probing of RNA structure.

The functions of RNA molecules are intimately linked to their ability to fold into complex secondary and tertiary structures. Thus, understanding how these molecules fold is essential to determining how they function. Current methods for investigating RNA structure often use small molecules, enzymes, or ions that cleave or modify the RNA in a solvent-accessible manner. While these methods have been invaluable to understanding RNA structure, they can be fairly labor intensive and often focus on short regions of single RNAs. Here we present a new method (Mod-seq) and data analysis pipeline (Mod-seeker) for assaying the structure of RNAs by high-throughput sequencing. This technique can be utilized both in vivo and in vitro, with any small molecule that modifies RNA and consequently impedes reverse transcriptase. As proof-of-principle, we used dimethyl sulfate (DMS) to probe the in vivo structure of total cellular RNAs in Saccharomyces cerevisiae. Mod-seq analysis simultaneously revealed secondary structural information for all four ribosomal RNAs and 32 additional noncoding RNAs. We further show that Mod-seq can be used to detect structural changes in 5.8S and 25S rRNAs in the absence of ribosomal protein L26, correctly identifying its binding site on the ribosome. While this method is applicable to RNAs of any length, its high-throughput nature makes Mod-seq ideal for studying long RNAs and complex RNA mixtures.

ROS-mediated SRMS activation confers platinum resistance in ovarian cancer.

Ovarian cancer is the leading cause of death among gynecological malignancies. Checkpoint blockade immunotherapy has so far only shown modest efficacy in ovarian cancer and platinum-based chemotherapy remains the front-line treatment. Development of platinum resistance is one of the most important factors contributing to ovarian cancer recurrence and mortality. Through kinome-wide synthetic lethal RNAi screening combined with unbiased datamining of cell line platinum response in CCLE and GDSC databases, here we report that Src-Related Kinase Lacking C-Terminal Regulatory Tyrosine And N-Terminal Myristylation Sites (SRMS), a non-receptor tyrosine kinase, is a novel negative regulator of MKK4-JNK signaling under platinum treatment and plays an important role in dictating platinum efficacy in ovarian cancer. Suppressing SRMS specifically sensitizes p53-deficient ovarian cancer cells to platinum in vitro and in vivo. Mechanistically, SRMS serves as a "sensor" for platinum-induced ROS. Platinum treatment-induced ROS activates SRMS, which inhibits MKK4 kinase activity by directly phosphorylating MKK4 at Y269 and Y307, and consequently attenuates MKK4-JNK activation. Suppressing SRMS leads to enhanced MKK4-JNK-mediated apoptosis by inhibiting MCL1 transcription, thereby boosting platinum efficacy. Importantly, through a "drug repurposing" strategy, we uncovered that PLX4720, a small molecular selective inhibitor of B-RafV600E, is a novel SRMS inhibitor that can potently boost platinum efficacy in ovarian cancer in vitro and in vivo. Therefore, targeting SRMS with PLX4720 holds the promise to improve the efficacy of platinum-based chemotherapy and overcome chemoresistance in ovarian cancer.

Genome-wide identification, characterization and functional analysis of the chitianse and chitinase-like gene family in Diaphorina citri.

BACKGROUND: Insect chitinases play vital roles in postembryonic development, especially during the molting process, and are potential targets for the RNA interference (RNAi)-based insecticidal strategy. Systematic functional analyses of chitinase genes have already been conducted on numerous insect pests, but similar analyses have not been carried out on Diaphorina citri. RESULTS: Eleven chitinase/chitinase-like genes and one endo-ß-N-acetylglucosaminidase (ENGase) gene were identified in the Diaphorina citri genome using various bioinformatic tools. Transcriptomes of the integument and midgut from fifth-instar nymphs and freshly-emerged adults of Diaphorina citri were generated and sequenced. Potential functions of 12 chitinase/chitinase-like genes were examined during nymph-adult transitions. Four chitinase genes, including DcCht5, DcCht7, DcCht10-1 and DcCht10-2, were mainly expressed in the integument of fifth-instar nymphs. These four genes were also up-regulated significantly under 20-hydroxyecdysone (20E) treatments. RNAi-mediated knockdown of these four genes suggests that they are essential for nymph-adult transition. CONCLUSION: Our results demonstrated essential roles of the chitinase/chitinase-like genes during the nymph-adult transition in Diaphorina citri, which are potentially useful targets for controlling the Diaphorina citri pest. © 2022 Society of Chemical Industry.

Generalized Score Functions for Causal Discovery.

Discovery of causal relationships from observational data is a fundamental problem. Roughly speaking, there are two types of methods for causal discovery, constraint-based ones and score-based ones. Score-based methods avoid the multiple testing problem and enjoy certain advantages compared to constraint-based ones. However, most of them need strong assumptions on the functional forms of causal mechanisms, as well as on data distributions, which limit their applicability. In practice the precise information of the underlying model class is usually unknown. If the above assumptions are violated, both spurious and missing edges may result. In this paper, we introduce generalized score functions for causal discovery based on the characterization of general (conditional) independence relationships between random variables, without assuming particular model classes. In particular, we exploit regression in RKHS to capture the dependence in a non-parametric way. The resulting causal discovery approach produces asymptotically correct results in rather general cases, which may have nonlinear causal mechanisms, a wide class of data distributions, mixed continuous and discrete data, and multidimensional variables. Experimental results on both synthetic and real-world data demonstrate the efficacy of our proposed approach.

Mod-seq: A High-Throughput Method for Probing RNA Secondary Structure.

It has become increasingly clear that large RNA molecules, especially long noncoding RNAs, function in almost all gene regulatory processes (Cech & Steitz, 2014). Many large RNAs appear to be structural scaffolds for assembly of important RNA/protein complexes. However, the structures of most large cellular RNA molecules are currently unknown (Hennelly & Sanbonmatsu, 2012). While chemical probing can reveal single-stranded regions of RNA, traditional approaches to identify sites of chemical modification are time consuming. Mod-seq is a high-throughput method used to map chemical modification sites on RNAs of any size, including complex mixtures of RNA. In this protocol, we describe preparation of Mod-seq high-throughput sequencing libraries from chemically modified RNA. We also describe a software package "Mod-seeker," which is a compilation of scripts written in Python, for the analysis of Mod-seq data. Mod-seeker returns statistically significant modification sites, which can then be used to aid in secondary structure prediction.