Refining the pool of RNA-binding domains advances the classification and prediction of RNA-binding proteins.

From transcription to decay, RNA-binding proteins (RBPs) influence RNA metabolism. Using the RBP2GO database that combines proteome-wide RBP screens from 13 species, we investigated the RNA-binding features of 176 896 proteins. By compiling published lists of RNA-binding domains (RBDs) and RNA-related protein family (Rfam) IDs with lists from the InterPro database, we analyzed the distribution of the RBDs and Rfam IDs in RBPs and non-RBPs to select RBDs and Rfam IDs that were enriched in RBPs. We also explored proteins for their content in intrinsically disordered regions (IDRs) and low complexity regions (LCRs). We found a strong positive correlation between IDRs and RBDs and a co-occurrence of specific LCRs. Our bioinformatic analysis indicated that RBDs/Rfam IDs were strong indicators of the RNA-binding potential of proteins and helped predicting new RBP candidates, especially in less investigated species. By further analyzing RBPs without RBD, we predicted new RBDs that were validated by RNA-bound peptides. Finally, we created the RBP2GO composite score by combining the RBP2GO score with new quality factors linked to RBDs and Rfam IDs. Based on the RBP2GO composite score, we compiled a list of 2018 high-confidence human RBPs. The knowledge collected here was integrated into the RBP2GO database at https://RBP2GO-2-Beta.dkfz.de.

Proteome-Wide Identification of RNA-Dependent Proteins in Lung Cancer Cells.

Following the concept of RNA dependence and exploiting its application in the R-DeeP screening approach, we have identified RNA-dependent proteins in A549 lung adenocarcinoma cells. RNA-dependent proteins are defined as proteins whose interactome depends on RNA and thus entails RNA-binding proteins (RBPs) as well as proteins in ribonucleoprotein complexes (RNPs) without direct RNA interaction. With this proteome-wide technique based on sucrose density gradient ultracentrifugation and fractionation followed by quantitative mass spectrometry and bioinformatic analysis, we have identified 1189 RNA-dependent proteins including 170 proteins which had never been linked to RNA before. R-DeeP provides quantitative information on the fraction of a protein being RNA-dependent as well as it allows the reconstruction of protein complexes based on co-segregation. The RNA dependence of three newly identified RNA-dependent proteins, DOCK5, ELMO2, also known as CED12A, and ABRAXAS1, also known as CCDC98, was validated using western blot analysis, and the direct RNA interaction was verified by iCLIP2 for the migration-related protein DOCK5 and the mitosis-related protein ABRAXAS1. The R-DeeP 2.0 database provides proteome-wide and cell line-specific information from A549 and HeLa S3 cells on proteins and their RNA dependence to contribute to understanding the functional role of RNA and RNA-binding proteins in cancer cells.

RBP2GO: a comprehensive pan-species database on RNA-binding proteins, their interactions and functions.

RNA-protein complexes have emerged as central players in numerous key cellular processes with significant relevance in health and disease. To further deepen our knowledge of RNA-binding proteins (RBPs), multiple proteome-wide strategies have been developed to identify RBPs in different species leading to a large number of studies contributing experimentally identified as well as predicted RBP candidate catalogs. However, the rapid evolution of the field led to an accumulation of isolated datasets, hampering the access and comparison of their valuable content. Moreover, tools to link RBPs to cellular pathways and functions were lacking. Here, to facilitate the efficient screening of the RBP resources, we provide RBP2GO (https://RBP2GO.DKFZ.de), a comprehensive database of all currently available proteome-wide datasets for RBPs across 13 species from 53 studies including 105 datasets identifying altogether 22 552 RBP candidates. These are combined with the information on RBP interaction partners and on the related biological processes, molecular functions and cellular compartments. RBP2GO offers a user-friendly web interface with an RBP scoring system and powerful advanced search tools allowing forward and reverse searches connecting functions and RBPs to stimulate new research directions.

circ2GO: A Database Linking Circular RNAs to Gene Function.

Circular RNAs (circRNAs) play critical roles in a broad spectrum of physiological and pathological processes, including cancer. Here, we provide a comprehensive database-circ2GO-systematically linking circRNAs to the functions and processes of their linear counterparts. circ2GO contains 148,811 circular human RNAs originating from 12,251 genes, which we derived from deep transcriptomics after rRNA depletion in a panel of 60 lung cancer and non-transformed cell lines. The broad circRNA expression dataset is mapped to all isoforms of the respective gene. The data are visualized in transcript maps and in heatmaps, to intuitively display a comprehensive portrait for the abundance of circRNAs across transcripts and cell lines. By integrating gene ontology (GO) information for all genes in our dataset, circ2GO builds a connection between circRNAs and their host genes' biological functions and molecular mechanisms. Additionally, circ2GO offers target predictions for circRNA-microRNA (miRNA) pairs for 25,166 highly abundant circRNAs from 6578 genes and 897 high-confidence human miRNAs. Visualization, user-friendliness, intuitive and advanced forward and reverse search options, batch processing and download options make circ2GO a comprehensive source for circRNA information to build hypotheses on their function, processes, and miRNA targets.

A pan-cancer analysis reveals nonstop extension mutations causing SMAD4 tumour suppressor degradation.

Nonstop or stop-loss mutations convert a stop into a sense codon, resulting in translation into the 3' untranslated region as a nonstop extension mutation to the next in-frame stop codon or as a readthrough mutation into the poly-A tail. Nonstop mutations have been characterized in hereditary diseases, but not in cancer genetics. In a pan-cancer analysis, we curated and analysed 3,412 nonstop mutations from 62 tumour entities, generating a comprehensive database at http://NonStopDB.dkfz.de. Six different nonstop extension mutations affected the tumour suppressor SMAD4, extending its carboxy terminus by 40 amino acids. These caused rapid degradation of the SMAD4 mutants via the ubiquitin-proteasome system. A hydrophobic degron signal sequence of ten amino acids within the carboxy-terminal extension was required to induce complete loss of the SMAD4 protein. Thus, we discovered that nonstop mutations can be functionally important in cancer and characterize their loss-of-function impact on the tumour suppressor SMAD4.

Identification, quantification and bioinformatic analysis of RNA-dependent proteins by RNase treatment and density gradient ultracentrifugation using R-DeeP.

Analysis of RNA-protein complexes is central to understanding the molecular circuitry governing cellular processes. In recent years, several proteome-wide studies have been dedicated to the identification of RNA-binding proteins. Here, we describe in detail R-DeeP, an approach built on RNA dependence, defined as the ability of a protein to engage in protein complexes only in the presence of RNA, involving direct or indirect interaction with RNA. This approach provides-for the first time, to our knowledge-quantitative information on the fraction of a protein associated with RNA-protein complexes. R-DeeP is independent of any potentially biased purification procedures. It is based on cellular lysate fractionation by density gradient ultracentrifugation and subsequent analysis by proteome-wide mass spectrometry (MS) or individual western blotting. The comparison of lysates with and without previous RNase treatment enables the identification of differences in the apparent molecular weight and, hence, the size of the complexes. In combination with information from databases of protein-protein complexes, R-DeeP facilitates the computational reconstruction of protein complexes from proteins migrating in the same fraction. In addition, we developed a pipeline for the statistical analysis of the MS dataset to automatically identify RNA-dependent proteins (proteins whose interactome depends on RNA). With this protocol, the individual analysis of proteins of interest by western blotting can be completed within 1-2 weeks. For proteome-wide studies, additional time is needed for the integration of the proteomic and statistical analyses. In the future, R-DeeP can be extended to other fractionation techniques, such as chromatography.

A pan-cancer analysis of synonymous mutations.

Synonymous mutations have been viewed as silent mutations, since they only affect the DNA and mRNA, but not the amino acid sequence of the resulting protein. Nonetheless, recent studies suggest their significant impact on splicing, RNA stability, RNA folding, translation or co-translational protein folding. Hence, we compile 659194 synonymous mutations found in human cancer and characterize their properties. We provide the user-friendly, comprehensive resource for synonymous mutations in cancer, SynMICdb ( http://SynMICdb.dkfz.de ), which also contains orthogonal information about gene annotation, recurrence, mutation loads, cancer association, conservation, alternative events, impact on mRNA structure and a SynMICdb score. Notably, synonymous and missense mutations are depleted at the 5'-end of the coding sequence as well as at the ends of internal exons independent of mutational signatures. For patient-derived synonymous mutations in the oncogene KRAS, we indicate that single point mutations can have a relevant impact on expression as well as on mRNA secondary structure.

R-DeeP: Proteome-wide and Quantitative Identification of RNA-Dependent Proteins by Density Gradient Ultracentrifugation.

The comprehensive but specific identification of RNA-binding proteins as well as the discovery of RNA-associated protein functions remain major challenges in RNA biology. Here we adapt the concept of RNA dependence, defining a protein as RNA dependent when its interactome depends on RNA. We converted this concept into a proteome-wide, unbiased, and enrichment-free screen called R-DeeP (RNA-dependent proteins), based on density gradient ultracentrifugation. Quantitative mass spectrometry identified 1,784 RNA-dependent proteins, including 537 lacking known links to RNA. Exploiting the quantitative nature of R-DeeP, proteins were classified as not, partially, or completely RNA dependent. R-DeeP identified the transcription factor CTCF as completely RNA dependent, and we uncovered that RNA is required for the CTCF-chromatin association. Additionally, R-DeeP allows reconstruction of protein complexes based on co-segregation. The whole dataset is available at http://R-DeeP.dkfz.de, providing proteome-wide, specific, and quantitative identification of proteins with RNA-dependent interactions and aiming at future functional discovery of RNA-protein complexes.

RNA motifs and combinatorial prediction of interactions, stability and localization of noncoding RNAs.

Although the number of documented noncoding RNAs (ncRNAs) is rapidly increasing, knowledge of their molecular function is lagging behind. The identification of specific RNA motifs that mediate transcript stability, interactions and localization may aid in the prediction of these features in new transcripts and may have potential implications for ncRNA function. Here, we review RNA motifs, focusing on four recent studies identifying nuclear-retention motifs, and discuss the limited specificity of short-RNA motifs and the resulting challenge for effective functional prediction. Future approaches may succeed by integrating combinatorial and cooperative effects of additional partially sequence-based properties.

Mitochondrial mutations in human cancer: Curation of translation.

As a genetic disease, cancer is caused by the activation of oncogenes and the inhibition of tumor suppressor genes via genetic and epigenetic mechanisms. Given the important role of energy metabolism in tumors, we analyzed the cancer-derived mutations occurring in the DNA of the mitochondrion. Mutations in the mitochondrial DNA (mtDNA) compared to nuclear DNA are 62% decreased relative to the coding length per chromosome. We find that the majority of these mutations affects highly conserved nucleotides - significantly exceeding the conservation of the mtDNA - and are devoid of single nucleotide polymorphisms (SNPs). Surprisingly, the leading resources for tumor genetics information universally use the standard genetic code for translation of nucleotide into amino acid sequences in their online resources. However, the nuclear and mitochondrial genetic codes differ for four codons and the usage of incomplete STOP codons. Hence, we analyze and curate the consequences for all mutations in the mtDNA and comprehensively reclassify missense, nonsense and synonymous mutations accordingly. In total, 10% of the mutations are incorrectly translated leading to significant changes in the distribution of mutation types with tripling of nonsense and 69% loss of nonstop extension mutations. Lastly, we provide a curated dataset of coding and non-coding mitochondrial mutations in cancer merged, standardized, duplicate-free and aggregated from two databases as a resource including orthogonal data on their high conservation and SNPs. This study generally highlights the need to universally regard the important differences between the standard and mitochondrial genetic code in life science research.

A cautionary tale of sense-antisense gene pairs: independent regulation despite inverse correlation of expression.

Long non-coding RNAs (lncRNAs) have been proven to play important roles in diverse cellular processes including the DNA damage response. Nearly 40% of annotated lncRNAs are transcribed in antisense direction to other genes and have often been implicated in their regulation via transcript- or transcription-dependent mechanisms. However, it remains unclear whether inverse correlation of gene expression would generally point toward a regulatory interaction between the genes. Here, we profiled lncRNA and mRNA expression in lung and liver cancer cells after exposure to DNA damage. Our analysis revealed two pairs of mRNA-lncRNA sense-antisense transcripts being inversely expressed upon DNA damage. The lncRNA NOP14-AS1 was strongly upregulated upon DNA damage, while the mRNA for NOP14 was downregulated, both in a p53-dependent manner. For another pair, the lncRNA LIPE-AS1 was downregulated, while its antisense mRNA CEACAM1 was upregulated. To test whether as expected the antisense genes would regulate each other resulting in this highly significant inverse correlation, we employed antisense oligonucleotides and RNAi to study transcript-dependent effects as well as dCas9-based transcriptional modulation by CRISPRi/CRISPRa for transcription-dependent effects. Surprisingly, despite the strong stimulus-dependent inverse correlation, our data indicate that neither transcript- nor transcription-dependent mechanisms explain the inverse regulation of NOP14-AS1:NOP14 or LIPE-AS1:CEACAM1 expression. Hence, sense-antisense pairs whose expression is strongly-positively or negatively-correlated can be nonetheless regulated independently. This highlights the requirement of individual experimental studies for each antisense pair and prohibits drawing conclusions on regulatory mechanisms from expression correlations.

Tethering RNA to chromatin for fluorescence microscopy based analysis of nuclear organization.

Nuclear RNAs emerge as important factors to orchestrate the dynamic organization of the nucleus into functional subcompartments. By tethering RNAs to distinct genomic loci, RNA-dependent chromatin changes can be dissected by fluorescence microscopic analysis. Here we describe how this approach is implemented in mammalian cells. It involves two high-affinity protein-nucleic acid interactions that can be established with a number of different protein domains and DNA and RNA sequences. A prototypic system is described here in detail: It consists of the binding of MS2 bacteriophage coat protein to its RNA recognition sequence and the interaction between the bacterial LacI repressor protein to its target lacO operator DNA sequence. Via these interactions RNAs tagged with the MS2 recognition sequences can be recruited to a locus with integrated lacO repeats. By inducing RNA-chromatin binding a number of RNA-dependent activities can be dissected: (i) The RNA-induced compaction or decondensation of chromatin, (ii) identification of RNA-interacting chromatin modifiers that set epigenetic signals such as posttranslational histone modifications, and (iii) nuclear relocation of a genomic locus targeted by the tethered RNA. Thus, a variety of RNA-dependent activities can be evaluated with the MS2-LacI system, which are crucial for understanding how RNA shapes nuclear organization.

The lncRNA VELUCT strongly regulates viability of lung cancer cells despite its extremely low abundance.

Little is known about the function of most non-coding RNAs (ncRNAs). The majority of long ncRNAs (lncRNAs) is expressed at very low levels and it is a matter of intense debate whether these can be of functional relevance. Here, we identified lncRNAs regulating the viability of lung cancer cells in a high-throughput RNA interference screen. Based on our previous expression profiling, we designed an siRNA library targeting 638 lncRNAs upregulated in human cancer. In a functional siRNA screen analyzing the viability of lung cancer cells, the most prominent hit was a novel lncRNA which we called Viability Enhancing LUng Cancer Transcript (VELUCT). In silico analyses confirmed the non-coding properties of the transcript. Surprisingly, VELUCT was below the detection limit in total RNA from NCI-H460 cells by RT-qPCR as well as RNA-Seq, but was robustly detected in the chromatin-associated RNA fraction. It is an extremely low abundant lncRNA with an RNA copy number of less than one copy per cell. Blocking transcription with actinomycin D revealed that VELUCT RNA was highly unstable which may partially explain its low steady-state concentration. Despite its extremely low abundance, loss-of-function of VELUCT with three independent experimental approaches in three different lung cancer cell lines led to a significant reduction of cell viability: Next to four individual siRNAs, also two complex siPOOLs as well as two antisense oligonucleotides confirmed the strong and specific phenotype. In summary, the extremely low abundant lncRNA VELUCT is essential for regulation of cell viability in several lung cancer cell lines. Hence, VELUCT is the first example for a lncRNA that is expressed at a very low level, but has a strong loss-of-function phenotype. Thus, our study proves that at least individual low-abundant lncRNAs can play an important functional role.

Regulation of nucleolus assembly by non-coding RNA polymerase II transcripts.

The nucleolus is a nuclear subcompartment for tightly regulated rRNA production and ribosome subunit biogenesis. It also acts as a cellular stress sensor and can release enriched factors in response to cellular stimuli. Accordingly, the content and structure of the nucleolus change dynamically, which is particularly evident during cell cycle progression: the nucleolus completely disassembles during mitosis and reassembles in interphase. Although the mechanisms that drive nucleolar (re)organization have been the subject of a number of studies, they are only partly understood. Recently, we identified Alu element-containing RNA polymerase II transcripts (aluRNAs) as important for nucleolar structure and rRNA synthesis. Integrating these findings with studies on the liquid droplet-like nature of the nucleolus leads us to propose a model on how RNA polymerase II transcripts could regulate the assembly of the nucleolus in response to external stimuli and during cell cycle progression.

Isolation of the protein and RNA content of active sites of transcription from mammalian cells.

Mammalian cell nuclei contain three RNA polymerases (RNAP I, RNAP II and RNAP III), which transcribe different gene subsets, and whose active forms are contained in supramolecular complexes known as 'transcription factories.' These complexes are difficult to isolate because they are embedded in the 3D structure of the nucleus. Factories exchange components with the soluble nucleoplasmic pool over time as gene expression programs change during development or disease. Analysis of their content can provide information on the nascent transcriptome and its regulators. Here we describe a protocol for the isolation of large factory fragments under isotonic salt concentrations in <72 h. It relies on DNase I-mediated detachment of chromatin from the nuclear substructure of freshly isolated, unfixed cells, followed by caspase treatment to release multi-megadalton factory complexes. These complexes retain transcriptional activity, and isolation of their contents is compatible with downstream analyses by mass spectrometry (MS) or RNA-sequencing (RNA-seq) to catalog the proteins and RNA associated with sites of active transcription.

Alu element-containing RNAs maintain nucleolar structure and function.

Non-coding RNAs play a key role in organizing the nucleus into functional subcompartments. By combining fluorescence microscopy and RNA deep-sequencing-based analysis, we found that RNA polymerase II transcripts originating from intronic Alu elements (aluRNAs) were enriched in the nucleolus. Antisense-oligo-mediated depletion of aluRNAs or drug-induced inhibition of RNA polymerase II activity disrupted nucleolar structure and impaired RNA polymerase I-dependent transcription of rRNA genes. In contrast, overexpression of a prototypic aluRNA sequence increased both nucleolus size and levels of pre-rRNA, suggesting a functional link between aluRNA, nucleolus integrity and pre-rRNA synthesis. Furthermore, we show that aluRNAs interact with nucleolin and target ectopic genomic loci to the nucleolus. Our study suggests an aluRNA-based mechanism that links RNA polymerase I and II activities and modulates nucleolar structure and rRNA production.

Dissecting the nascent human transcriptome by analysing the RNA content of transcription factories.

While mapping total and poly-adenylated human transcriptomes has now become routine, characterizing nascent transcripts remains challenging, largely because nascent RNAs have such short half-lives. Here, we describe a simple, fast and cost-effective method to isolate RNA associated with transcription factories, the sites responsible for the majority of nuclear transcription. Following stimulation of human endothelial cells with the pro-inflammatory cytokine TNFα, we isolate and analyse the RNA content of factories by sequencing. Comparison with total, poly(A)(+) and chromatin RNA fractions reveals that sequencing of purified factory RNA maps the complete nascent transcriptome; it is rich in intronic unprocessed transcript, as well as long intergenic non-coding (lincRNAs) and enhancer-associated RNAs (eRNAs), micro-RNA precursors and repeat-derived RNAs. Hence, we verify that transcription factories produce most nascent RNA and confer a regulatory role via their association with a set of specifically-retained non-coding transcripts.

Loss of cooperativity of secreted CD40L and increased dose-response to IL4 on CLL cell viability correlates with enhanced activation of NF-kB and STAT6.

Chronic lymphocytic leukemia (CLL) cells fail to enter apoptosis in vivo as opposed to their non-malignant B-lymphocyte counterparts. The ability of CLL cells to escape apoptosis is highly dependent on their microenvironment. Compared to non-malignant B cells, CLL cells are more responsive to complex stimuli that can be reproduced in vitro by the addition of cytokines. To understand the molecular mechanism of the environment-dependent anti-apoptotic signaling circuitry of CLL cells, we quantified the effect of the SDF-1, BAFF, APRIL, anti-IgM, interleukin-4 (IL4) and secreted CD40L (sCD40L) on the survival of in vitro cultured CLL cells and found IL4 and sCD40L to be most efficient in rescuing CLL cells from apoptosis. In quantitative dose-response experiments using cell survival as readout, the binding affinity of IL4 to its receptor was similar between malignant and non-malignant cells. However, the downstream signaling in terms of the amount of STAT6 and its degree of phosphorylation was highly stimulated in CLL cells. In contrast, the response to sCD40L showed a loss of cooperative binding in CLL cells but displayed a largely increased ligand binding affinity. Although a high-throughput microscopy analysis did not reveal a significant difference in the spatial CD40 receptor organization, the downstream signaling showed an enhanced activation of the NF-kB pathway in the malignant cells. Thus, we propose that the anti-apoptotic phenotype of CLL involves a sensitized response for IL4 dependent STAT6 phosphorylation, and an activation of NF-kB signaling due to an increased affinity of sCD40L to its receptor.

Epigenetic upregulation of lncRNAs at 13q14.3 in leukemia is linked to the In Cis downregulation of a gene cluster that targets NF-kB.

Non-coding RNAs are much more common than previously thought. However, for the vast majority of non-coding RNAs, the cellular function remains enigmatic. The two long non-coding RNA (lncRNA) genes DLEU1 and DLEU2 map to a critical region at chromosomal band 13q14.3 that is recurrently deleted in solid tumors and hematopoietic malignancies like chronic lymphocytic leukemia (CLL). While no point mutations have been found in the protein coding candidate genes at 13q14.3, they are deregulated in malignant cells, suggesting an epigenetic tumor suppressor mechanism. We therefore characterized the epigenetic makeup of 13q14.3 in CLL cells and found histone modifications by chromatin-immunoprecipitation (ChIP) that are associated with activated transcription and significant DNA-demethylation at the transcriptional start sites of DLEU1 and DLEU2 using 5 different semi-quantitative and quantitative methods (aPRIMES, BioCOBRA, MCIp, MassARRAY, and bisulfite sequencing). These epigenetic aberrations were correlated with transcriptional deregulation of the neighboring candidate tumor suppressor genes, suggesting a coregulation in cis of this gene cluster. We found that the 13q14.3 genes in addition to their previously known functions regulate NF-kB activity, which we could show after overexpression, siRNA-mediated knockdown, and dominant-negative mutant genes by using Western blots with previously undescribed antibodies, by a customized ELISA as well as by reporter assays. In addition, we performed an unbiased screen of 810 human miRNAs and identified the miR-15/16 family of genes at 13q14.3 as the strongest inducers of NF-kB activity. In summary, the tumor suppressor mechanism at 13q14.3 is a cluster of genes controlled by two lncRNA genes that are regulated by DNA-methylation and histone modifications and whose members all regulate NF-kB. Therefore, the tumor suppressor mechanism in 13q14.3 underlines the role both of epigenetic aberrations and of lncRNA genes in human tumorigenesis and is an example of colocalization of a functionally related gene cluster.