Zhipeng Lu

RNA Duplex Map in Living Cells Reveals Higher-Order Transcriptome Structure

RNA has the intrinsic property to base pair, forming complex structures fundamental to its diverse functions. Here, we develop PARIS, a method based on reversible psoralen crosslinking for global mapping of RNA duplexes with near base-pair resolution in living cells. PARIS analysis in three human and mouse cell types reveals frequent long-range structures, higher-order architectures, and RNA-RNA interactions in trans across the transcriptome. PARIS determines base-pairing interactions on an individual-molecule level, revealing pervasive alternative conformations. We used PARIS-determined helices to guide phylogenetic analysis of RNA structures and discovered conserved long-range and alternative structures. XIST, a long noncoding RNA (lncRNA) essential for X chromosome inactivation, folds into evolutionarily conserved RNA structural domains that span many kilobases. XIST A-repeat forms complex inter-repeat duplexes that nucleate higher-order assembly of the key epigenetic silencing protein SPEN. PARIS is a generally applicable and versatile method that provides novel insights into the RNA structurome and interactome. VIDEO ABSTRACT.

REV3L confers chemoresistance to cisplatin in human gliomas: The potential of its RNAi for synergistic therapy

The REV3L gene, encoding the catalytic subunit of human polymerase zeta, plays a significant role in the cytotoxicity, mutagenicity, and chemoresistance of certain tumors. However, the role of REV3L in regulating the sensitivity of glioma cells to chemotherapy remains unknown. In this study, we investigated the expression of the REV3L gene in 10 normal brain specimens and 30 human glioma specimens and examined the value of REV3L as a potential modulator of cellular response to various DNA-damaging agents. Reverse transcriptase PCR/real-time PCR analysis revealed that REV3L was overexpressed in human gliomas compared with normal brain tissues. A glioma cell model with stable overexpression of REV3L was used to probe the role of REV3L in cisplatin treatment; upregulation of REV3L markedly attenuated cisplatin-induced apoptosis of the mitochondrial apoptotic pathway. We therefore assessed the REV3L-targeted treatment modality that combines suppression of REV3L expression using RNA interference (RNAi) with the cytotoxic effects of DNA-damaging agents. Downregulation of REV3L expression significantly enhanced the sensitivity of glioma cells to cisplatin, as evidenced by the increased apoptosis rate and marked alterations in the anti-apoptotic proteins B-cell lymphoma 2 (Bcl-2) and B-cell lymphoma-extra large (Bcl-xl) and proapoptotic Bcl-2-associated x protein (Bax) expression levels, and reduced mutation frequencies in surviving glioma cells. These results suggest that REV3L may potentially contribute to gliomagenesis and play a crucial role in regulating cellular response to the DNA cross-linking agent cisplatin. Our findings indicate that RNAi targeting REV3L combined with chemotherapy has synergistic therapeutic effects on glioma cells, which warrants further investigation as an effective novel therapeutic regimen for patients with this malignancy.

Developmental arrest of Drosophila survival motor neuron (Smn) mutants accounts for differences in expression of minor intron-containing genes

Reduced levels of survival motor neuron (SMN) protein lead to a neuromuscular disease called spinal muscular atrophy (SMA). Animal models of SMA recapitulate many aspects of the human disease, including locomotion and viability defects, but have thus far failed to uncover the causative link between a lack of SMN protein and neuromuscular dysfunction. While SMN is known to assemble small nuclear ribonucleoproteins (snRNPs) that catalyze pre-mRNA splicing, it remains unclear whether disruptions in splicing are etiologic for SMA. To investigate this issue, we carried out RNA deep-sequencing (RNA-seq) on age-matched Drosophila Smn-null and wild-type larvae. Comparison of genome-wide mRNA expression profiles with publicly available data sets revealed the timing of a developmental arrest in the Smn mutants. Furthermore, genome-wide differences in splicing between wild-type and Smn animals did not correlate with changes in mRNA levels. Specifically, we found that mRNA levels of genes that contain minor introns vary more over developmental time than they do between wild-type and Smn mutants. An analysis of reads mapping to minor-class intron-exon junctions revealed only small changes in the splicing of minor introns in Smn larvae, within the normal fluctuations that occur throughout development. In contrast, Smn mutants displayed a prominent increase in levels of stress-responsive transcripts, indicating a systemic response to the developmental arrest induced by loss of SMN protein. These findings not only provide important mechanistic insight into the developmental arrest displayed by Smn mutants, but also argue against a minor-intron-dependent etiology for SMA.

Decoding the RNA structurome.

Structures of RNA molecules are essential for their architectural, regulatory, and catalytic functions. Recent advances in high throughput sequencing enabled the development of methods for probing RNA structures on a transcriptome-wide scale-termed the RNA structurome. Here we review the state-of-the-art technologies for probing the RNA structurome, and highlight insights gained from these studies. We also point out the limits of current methods and discuss potential directions for future improvements.

RIP-seq analysis of eukaryotic Sm proteins identifies three major categories of Sm-containing ribonucleoproteins

BackgroundSm proteins are multimeric RNA-binding factors, found in all three domains of life. Eukaryotic Sm proteins, together with their associated RNAs, form small ribonucleoprotein (RNP) complexes important in multiple aspects of gene regulation. Comprehensive knowledge of the RNA components of Sm RNPs is critical for understanding their functions.ResultsWe developed a multi-targeting RNA-immunoprecipitation sequencing (RIP-seq) strategy to reliably identify Sm-associated RNAs from Drosophila ovaries and cultured human cells. Using this method, we discovered three major categories of Sm-associated transcripts: small nuclear (sn)RNAs, small Cajal body (sca)RNAs and mRNAs. Additional RIP-PCR analysis showed both ubiquitous and tissue-specific interactions. We provide evidence that the mRNA-Sm interactions are mediated by snRNPs, and that one of the mechanisms of interaction is via base pairing. Moreover, the Sm-associated mRNAs are mature, indicating a splicing-independent function for Sm RNPs.ConclusionsThis study represents the first comprehensive analysis of eukaryotic Sm-containing RNPs, and provides a basis for additional functional analyses of Sm proteins and their associated snRNPs outside of the context of pre-mRNA splicing. Our findings expand the repertoire of eukaryotic Sm-containing RNPs and suggest new functions for snRNPs in mRNA metabolism.

Developmental Analysis of Spliceosomal snRNA Isoform Expression

Pre-mRNA splicing is a critical step in eukaryotic gene expression that contributes to proteomic, cellular, and developmental complexity. Small nuclear (sn)RNAs are core spliceosomal components; however, the extent to which differential expression of snRNA isoforms regulates splicing is completely unknown. This is partly due to difficulties in the accurate analysis of the spatial and temporal expression patterns of snRNAs. Here, we use high-throughput RNA-sequencing (RNA-seq) data to profile expression of four major snRNAs throughout Drosophila development. This analysis shows that individual isoforms of each snRNA have distinct expression patterns in the embryo, larva, and pharate adult stages. Expression of these isoforms is more heterogeneous during embryogenesis; as development progresses, a single isoform from each snRNA subtype gradually dominates expression. Despite the lack of stable snRNA orthologous groups during evolution, this developmental switching of snRNA isoforms also occurs in distantly related vertebrate species, such as Xenopus, mouse, and human. Our results indicate that expression of snRNA isoforms is regulated and lays the foundation for functional studies of individual snRNA isoforms.

Vicinal: a method for the determination of ncRNA ends using chimeric reads from RNA-seq experiments

Non-coding (nc)RNAs are important structural and regulatory molecules. Accurate determination of the primary sequence and secondary structure of ncRNAs is important for understanding their functions. During cDNA synthesis, RNA 3′ end stem-loops can self-prime reverse transcription, creating RNA–cDNA chimeras. We found that chimeric RNA–cDNA fragments can also be detected at 5′ end stem-loops, although at much lower frequency. Using the Gubler–Hoffman method, both types of chimeric fragments can be converted to cDNA during library construction, and they are readily detectable in high-throughput RNA sequencing (RNA-seq) experiments. Here, we show that these chimeric reads contain valuable information about the boundaries of ncRNAs. We developed a bioinformatic method, called Vicinal, to precisely map the ends of numerous fruitfly, mouse and human ncRNAs. Using this method, we analyzed chimeric reads from over 100 RNA-seq datasets, the results of which we make available for users to find RNAs of interest. In summary, we show that Vicinal is a useful tool for determination of the precise boundaries of uncharacterized ncRNAs, facilitating further structure/function studies.

Mechanistic insights in X-chromosome inactivation

X-chromosome inactivation (XCI) is a critical epigenetic mechanism for balancing gene dosage between XY males and XX females in eutherian mammals. A long non-coding RNA (lncRNA), XIST, and its associated proteins orchestrate this multi-step process, resulting in the inheritable silencing of one of the two X-chromosomes in females. The XIST RNA is large and complex, exemplifying the unique challenges associated with the structural and functional analysis of lncRNAs. Recent technological advances in the analysis of macromolecular structure and interactions have enabled us to systematically dissect the XIST ribonucleoprotein complex, which is larger than the ribosome, and its place of action, the inactive X-chromosome. These studies shed light on key mechanisms of XCI, such as XIST coating of the X-chromosome, recruitment of DNA, RNA and histone modification enzymes, and compaction and compartmentalization of the inactive X. Here, we summarize recent studies on XCI, highlight the critical contributions of new technologies and propose a unifying model for XIST function in XCI where modular domains serve as the structural and functional units in both lncRNA–protein complexes and DNA–protein complexes in chromatin. This article is part of the themed issue ‘X-chromosome inactivation: a tribute to Mary Lyon’.

RISE: a database of RNA interactome from sequencing experiments

Abstract We present RISE (http://rise.zhanglab.net), a database of RNA Interactome from Sequencing Experiments. RNA-RNA interactions (RRIs) are essential for RNA regulation and function. RISE provides a comprehensive collection of RRIs that mainly come from recent transcriptome-wide sequencing-based experiments like PARIS, SPLASH, LIGR-seq, and MARIO, as well as targeted studies like RIA-seq, RAP-RNA and CLASH. It also includes interactions aggregated from other primary databases and publications. The RISE database currently contains 328,811 RNA-RNA interactions mainly in human, mouse and yeast. While most existing RNA databases mainly contain interactions of miRNA targeting, notably, more than half of the RRIs in RISE are among mRNA and long non-coding RNAs. We compared different RRI datasets in RISE and found limited overlaps in interactions resolved by different techniques and in different cell lines. It may suggest technology preference and also dynamic natures of RRIs. We also analyzed the basic features of the human and mouse RRI networks and found that they tend to be scale-free, small-world, hierarchical and modular. The analysis may nominate important RNAs or RRIs for further investigation. Finally, RISE provides a Circos plot and several table views for integrative visualization, with extensive molecular and functional annotations to facilitate exploration of biological functions for any RRI of interest.

PARIS: Psoralen Analysis of RNA Interactions and Structures with High Throughput and Resolution

RNA has the intrinsic propensity to form base pairs, leading to complex intramolecular and intermolecular helices. Direct measurement of base pairing interactions in living cells is critical to solving transcriptome structure and interactions, and investigating their functions (Lu and Chang, Curr Opin Struct Biol 36:142-148, 2016). Toward this goal, we developed an experimental method, PARIS (Psoralen Analysis of RNA Interactions and Structures), to directly determine transcriptome-wide base pairing interactions (Lu et al., Cell 165(5):1267-1279, 2016). PARIS combines four critical steps, in vivo cross-linking, 2D gel purification, proximity ligation, and high-throughput sequencing to achieve high-throughput and near-base pair resolution determination of the RNA structurome and interactome in living cells. In this chapter, we aim to provide a comprehensive discussion on the principles behind the experimental and computational strategies, and a step-by-step description of the experiment and analysis.