AURA: Atlas of UTR Regulatory Activity

Connecting the dots: A practical evaluation of web-tools for describing protein dynamics as networks

Protein Structure Networks (PSNs) are a well-known mathematical model for estimation and analysis of the three-dimensional protein structure. Investigating the topological architecture of PSNs may help identify the crucial amino acid residues for protein stability and protein-protein interactions, as well as deduce any possible mutational effects. But because proteins go through conformational changes to give rise to essential biological functions, this has to be done dynamically over time. The most effective method to describe protein dynamics is molecular dynamics simulation, with the most popular software programs for manipulating simulations to infer interaction networks being RING, MD-TASK, and NAPS. Here, we compare the computational approaches used by these three tools-all of which are accessible as web servers-to understand the pathogenicity of missense mutations and talk about their potential applications as well as their advantages and disadvantages.

Comprehensive characterization of posttranscriptional impairment-related 3'-UTR mutations in 2413 whole genomes of cancer patients

The 3' untranslated region (3'-UTR) is the vital element regulating gene expression, but most studies have focused on variations in RNA-binding proteins (RBPs), miRNAs, alternative polyadenylation (APA) and RNA modifications. To explore the posttranscriptional function of 3'-UTR somatic mutations in tumorigenesis, we collected whole-genome data from 2413 patients across 18 cancer types. Our updated algorithm, PIVar, revealed 25,216 3'-UTR posttranscriptional impairment-related SNVs (3'-UTR piSNVs) spanning 2930 genes; 24 related RBPs were significantly enriched. The somatic 3'-UTR piSNV ratio was markedly increased across all 18 cancer types, which was associated with worse survival for four cancer types. Several cancer-related genes appeared to facilitate tumorigenesis at the protein and posttranscriptional regulation levels, whereas some 3'-UTR piSNV-affected genes functioned mainly via posttranscriptional mechanisms. Moreover, we assessed immune cell and checkpoint characteristics between the high/low 3'-UTR piSNV ratio groups and predicted 80 compounds associated with the 3'-UTR piSNV-affected gene expression signature. In summary, our study revealed the prevalence and clinical relevance of 3'-UTR piSNVs in cancers, and also demonstrates that in addition to affecting miRNAs, 3'-UTR piSNVs perturb RBPs binding, APA and m6A RNA modification, which emphasized the importance of considering 3'-UTR piSNVs in cancer biology.

Genome-wide dynamics of RNA synthesis, processing, and degradation without RNA metabolic labeling

The quantification of the kinetic rates of RNA synthesis, processing, and degradation are largely based on the integrative analysis of total and nascent transcription, the latter being quantified through RNA metabolic labeling. We developed INSPEcT−, a computational method based on the mathematical modeling of premature and mature RNA expression that is able to quantify kinetic rates from steady-state or time course total RNA-seq data without requiring any information on nascent transcripts. Our approach outperforms available solutions, closely recapitulates the kinetic rates obtained through RNA metabolic labeling, improves the ability to detect changes in transcript half-lives, reduces the cost and complexity of the experiments, and can be adopted to study experimental conditions in which nascent transcription cannot be readily profiled. Finally, we applied INSPEcT− to the characterization of post-transcriptional regulation landscapes in dozens of physiological and disease conditions. This approach was included in the INSPEcT Bioconductor package, which can now unveil RNA dynamics from steady-state or time course data, with or without the profiling of nascent RNA.

Prevalence and architecture of posttranscriptionally impaired synonymous mutations in 8,320 genomes across 22 cancer types

Somatic synonymous mutations are one of the most frequent genetic variants occurring in the coding region of cancer genomes, while their contributions to cancer development remain largely unknown. To assess whether synonymous mutations involved in post-transcriptional regulation contribute to the genetic etiology of cancers, we collected whole exome data from 8,320 patients across 22 cancer types. By employing our developed algorithm, PIVar, we identified a total of 22,948 posttranscriptionally impaired synonymous SNVs (pisSNVs) spanning 2,042 genes. In addition, 35 RNA binding proteins impacted by these identified pisSNVs were significantly enriched. Remarkably, we discovered markedly elevated ratio of somatic pisSNVs across all 22 cancer types, and a high pisSNV ratio was associated with worse patient survival in five cancer types. Intriguing, several well-established cancer genes, including PTEN, RB1 and PIK3CA, appeared to contribute to tumorigenesis at both protein function and posttranscriptional regulation levels, whereas some pisSNV-hosted genes, including UBR4, EP400 and INTS1, exerted their function during carcinogenesis mainly via posttranscriptional mechanisms. Moreover, we predicted three drugs associated with two pisSNVs, and numerous compounds associated with expression signature of pisSNV-hosted genes. Our study reveals the prevalence and clinical relevance of pisSNVs in cancers, and emphasizes the importance of considering posttranscriptional impaired synonymous mutations in cancer biology.

Discovering the 3' UTR-mediated regulation of alpha-synuclein

Recent evidence indicates a link between Parkinson's Disease (PD) and the expression of a-synuclein (SNCA) isoforms with different 3′ untranslated regions (3′UTRs). Yet, the post-transcriptional mechanisms regulating SNCA expression are unknown. Using a large-scale in vitro /in silico screening we identified RNA-binding proteins (RBPs) that interact with SNCA 3′ UTRs. We identified two RBPs, ELAVL1 and TIAR, that bind with high affinity to the most abundant and translationally active 3′ UTR isoform (575 nt). Knockdown and overexpression experiments indicate that both ELAVL1 and TIAR positively regulate endogenous SNCA in vivo. The mechanism of regulation implies mRNA stabilization as well as enhancement of translation in the case of TIAR. We observed significant alteration of both TIAR and ELAVL1 expression in motor cortex of post-mortem brain donors and primary cultured fibroblast from patients affected by PD and Multiple System Atrophy (MSA). Moreover, trans expression quantitative trait loci (trans-eQTLs) analysis revealed that a group of single nucleotide polymorphisms (SNPs) in TIAR genomic locus influences SNCA expression in two different brain areas, nucleus accumbens and hippocampus. Our study sheds light on the 3′ UTR-mediated regulation of SNCA and its link with PD pathogenesis, thus opening up new avenues for investigation of post-transcriptional mechanisms in neurodegeneration.

Associating transcription factors and conserved RNA structures with gene regulation in the human brain

Anatomical subdivisions of the human brain can be associated with different neuronal functions. This functional diversification is reflected by differences in gene expression. By analyzing post-mortem gene expression data from the Allen Brain Atlas, we investigated the impact of transcription factors (TF) and RNA secondary structures on the regulation of gene expression in the human brain. First, we modeled the expression of a gene as a linear combination of the expression of TFs. We devised an approach to select robust TF-gene interactions and to determine localized contributions to gene expression of TFs. Among the TFs with the most localized contributions, we identified EZH2 in the cerebellum, NR3C1 in the cerebral cortex and SRF in the basal forebrain. Our results suggest that EZH2 is involved in regulating ZIC2 and SHANK1 which have been linked to neurological diseases such as autism spectrum disorder. Second, we associated enriched regulatory elements inside differentially expressed mRNAs with RNA secondary structure motifs. We found a group of purine-uracil repeat RNA secondary structure motifs plus other motifs in neuron related genes such as ACSL4 and ERLIN2.

The 5'-untranslated region of p16INK4a melanoma tumor suppressor acts as a cellular IRES, controlling mRNA translation under hypoxia through YBX1 binding

CDKN2A/p16INK4a is an essential tumor suppressor gene that controls cell cycle progression and replicative senescence. It is also the main melanoma susceptibility gene. Here we report that p16INK4a 5'UTR mRNA acts as a cellular Internal Ribosome Entry Site (IRES). The potential for p16INK4a 5'UTR to drive cap-independent translation was evaluated by dual-luciferase assays using bicistronic and monocistronic vectors. Results of reporters' relative activities coupled to control analyses for actual bicistronic mRNA transcription, indicated that the wild type p16INK4a 5'UTR could stimulate cap-independent translation. Notably, hypoxic stress and the treatment with mTOR inhibitors enhanced the translation-stimulating property of p16INK4a 5'UTR. RNA immunoprecipitation performed in melanoma-derived SK-Mel-28 and in a patient-derived lymphoblastoid cell line indicated that YBX1 can bind the wild type p16INK4a mRNA increasing its translation efficiency, particularly during hypoxic stress. Modulation of YBX1 expression further supported its involvement in cap-independent translation of the wild type p16INK4a but not a c.-42T>A variant. RNA SHAPE assays revealed local flexibility changes for the c.-42T>A variant at the predicted YBX1 binding site region. Our results indicate that p16INK4a 5'UTR contains a cellular IRES that can enhance mRNA translation efficiency, in part through YBX1.

RBP-Var: a database of functional variants involved in regulation mediated by RNA-binding proteins

Transcription factors bind to the genome by forming specific contacts with the primary DNA sequence; however, RNA-binding proteins (RBPs) have greater scope to achieve binding specificity through the RNA secondary structure. It has been revealed that single nucleotide variants (SNVs) that alter RNA structure, also known as RiboSNitches, exhibit 3-fold greater local structure changes than replicates of the same DNA sequence, demonstrated by the fact that depletion of RiboSNitches could result in the alteration of specific RNA shapes at thousands of sites, including 3' UTRs, binding sites of microRNAs and RBPs. However, the network between SNVs and post-transcriptional regulation remains unclear. Here, we developed RBP-Var, a database freely available at http://www.rbp-var.biols.ac.cn/, which provides annotation of functional variants involved in post-transcriptional interaction and regulation. RBP-Var provides an easy-to-use web interface that allows users to rapidly find whether SNVs of interest can transform the secondary structure of RNA and identify RBPs whose binding may be subsequently disrupted. RBP-Var integrates DNA and RNA biology to understand how various genetic variants and post-transcriptional mechanisms cooperate to orchestrate gene expression. In summary, RBP-Var is useful in selecting candidate SNVs for further functional studies and exploring causal SNVs underlying human diseases.

Optimizing RNA structures by sequence extensions using RNAcop

A key aspect of RNA secondary structure prediction is the identification of novel functional elements. This is a challenging task because these elements typically are embedded in longer transcripts where the borders between the element and flanking regions have to be defined. The flanking sequences impact the folding of the functional elements both at the level of computational analyses and when the element is extracted as a transcript for experimental analysis. Here, we analyze how different flanking region lengths impact folding into a constrained structure by computing probabilities of folding for different sizes of flanking regions. Our method, RNAcop (RNA context optimization by probability), is tested on known and de novo predicted structures. In vitro experiments support the computational analysis and suggest that for a number of structures, choosing proper lengths of flanking regions is critical. RNAcop is available as web server and stand-alone software via http://rth.dk/resources/rnacop.

DoRiNA 2.0--upgrading the doRiNA database of RNA interactions in post-transcriptional regulation

The expression of almost all genes in animals is subject to post-transcriptional regulation by RNA binding proteins (RBPs) and microRNAs (miRNAs). The interactions between both RBPs and miRNAs with mRNA can be mapped on a whole-transcriptome level using experimental and computational techniques established in the past years. The combined action of RBPs and miRNAs is thought to form a post-transcriptional regulatory code. Here we present doRiNA 2.0, available at http://dorina.mdc-berlin.de. In this highly improved new version, we have completely reworked the user interface and expanded the database to improve the usability of the website. Taking into account user feedback over the past years, the input forms for both the simple and the combinatorial search function have been streamlined and combined into a single web page that will also display the search results. Especially, custom uploads is one of the key new features in doRiNA 2.0. To enable the inclusion of doRiNA into third-party analysis pipelines, all operations are accessible via a REST API. Alternatively, local installations can be queried using a Python API. Both the web application and the APIs are available under an OSI-approved Open Source license that allows research and commercial access and re-use.

Fingerprints of a message: integrating positional information on the transcriptome

The recent explosion of high-throughput sequencing methods applied to RNA molecules is allowing us to go beyond the description of sequence variants and their relative abundances, as measured by RNA-seq. We can now probe for RNA engagement in polysomes, for ribosomes, RNA binding proteins and microRNAs binding sites, for RNA secondary structure and for RNA methylation. These descriptors produce a steadily growing multidimensional array of positional information on RNA sequences, whose effective integration only would bring to decipher the regulatory interplay occurring between proteins, RNAs and their modifications on the transcriptome. This interplay ultimately dictates the degree of mRNA availability to translation, and thus the occurrence of cell phenotypes. However, several issues in data presentation are slowing down effective integration. A standardization effort for new dataset types produced should be urgently undertaken to solve these issues. Providing uniformed experimental details along with datasets processed to be directly usable and employing shared formats would greatly simplify integration efforts, strengthening hypotheses stemming from correlative observations and eventually bringing to mechanistic understanding.

APADB: a database for alternative polyadenylation and microRNA regulation events

Alternative polyadenylation (APA) is a widespread mechanism that contributes to the sophisticated dynamics of gene regulation. Approximately 50% of all protein-coding human genes harbor multiple polyadenylation (PA) sites; their selective and combinatorial use gives rise to transcript variants with differing length of their 3′ untranslated region (3′UTR). Shortened variants escape UTR-mediated regulation by microRNAs (miRNAs), especially in cancer, where global 3′UTR shortening accelerates disease progression, dedifferentiation and proliferation. Here we present APADB, a database of vertebrate PA sites determined by 3′ end sequencing, using massive analysis of complementary DNA ends. APADB provides (A)PA sites for coding and non-coding transcripts of human, mouse and chicken genes. For human and mouse, several tissue types, including different cancer specimens, are available. APADB records the loss of predicted miRNA binding sites and visualizes next-generation sequencing reads that support each PA site in a genome browser. The database tables can either be browsed according to organism and tissue or alternatively searched for a gene of interest. APADB is the largest database of APA in human, chicken and mouse. The stored information provides experimental evidence for thousands of PA sites and APA events. APADB combines 3′ end sequencing data with prediction algorithms of miRNA binding sites, allowing to further improve prediction algorithms. Current databases lack correct information about 3′UTR lengths, especially for chicken, and APADB provides necessary information to close this gap.

Database URL:http://tools.genxpro.net/apadb/

Scan for Motifs: a webserver for the analysis of post-transcriptional regulatory elements in the 3' untranslated regions (3' UTRs) of mRNAs

Background

Gene expression in vertebrate cells may be controlled post-transcriptionally through regulatory elements in mRNAs. These are usually located in the untranslated regions (UTRs) of mRNA sequences, particularly the 3′UTRs.

Results

Scan for Motifs (SFM) simplifies the process of identifying a wide range of regulatory elements on alignments of vertebrate 3′UTRs. SFM includes identification of both RNA Binding Protein (RBP) sites and targets of miRNAs. In addition to searching pre-computed alignments, the tool provides users the flexibility to search their own sequences or alignments. The regulatory elements may be filtered by expected value cutoffs and are cross-referenced back to their respective sources and literature. The output is an interactive graphical representation, highlighting potential regulatory elements and overlaps between them. The output also provides simple statistics and links to related resources for complementary analyses. The overall process is intuitive and fast. As SFM is a free web-application, the user does not need to install any software or databases.

Conclusions

Visualisation of the binding sites of different classes of effectors that bind to 3′UTRs will facilitate the study of regulatory elements in 3′ UTRs.

Protein-specific prediction of mRNA binding using RNA sequences, binding motifs and predicted secondary structures

Background

RNA-binding proteins interact with specific RNA molecules to regulate important cellular processes. It is therefore necessary to identify the RNA interaction partners in order to understand the precise functions of such proteins. Protein-RNA interactions are typically characterized using in vivo and in vitro experiments but these may not detect all binding partners. Therefore, computational methods that capture the protein-dependent nature of such binding interactions could help to predict potential binding partners in silico.

Results

We have developed three methods to predict whether an RNA can interact with a particular RNA-binding protein using support vector machines and different features based on the sequence (the Oli method), the motif score (the OliMo method) and the secondary structure (the OliMoSS method). We applied these approaches to different experimentally-derived datasets and compared the predictions with RNAcontext and RPISeq. Oli outperformed OliMoSS and RPISeq, confirming our protein-specific predictions and suggesting that tetranucleotide frequencies are appropriate discriminative features. Oli and RNAcontext were the most competitive methods in terms of the area under curve. A precision-recall curve analysis achieved higher precision values for Oli. On a second experimental dataset including real negative binding information, Oli outperformed RNAcontext with a precision of 0.73 vs. 0.59.

Conclusions

Our experiments showed that features based on primary sequence information are sufficiently discriminating to predict specific RNA-protein interactions. Sequence motifs and secondary structure information were not necessary to improve these predictions. Finally we confirmed that protein-specific experimental data concerning RNA-protein interactions are valuable sources of information that can be used for the efficient training of models for in silico predictions. The scripts are available upon request to the corresponding author.

Constitutive patterns of gene expression regulated by RNA-binding proteins

Background

RNA-binding proteins regulate a number of cellular processes, including synthesis, folding, translocation, assembly and clearance of RNAs. Recent studies have reported that an unexpectedly large number of proteins are able to interact with RNA, but the partners of many RNA-binding proteins are still uncharacterized.

Results

We combined prediction of ribonucleoprotein interactions, based on catRAPID calculations, with analysis of protein and RNA expression profiles from human tissues. We found strong interaction propensities for both positively and negatively correlated expression patterns. Our integration of in silico and ex vivo data unraveled two major types of protein–RNA interactions, with positively correlated patterns related to cell cycle control and negatively correlated patterns related to survival, growth and differentiation. To facilitate the investigation of protein–RNA interactions and expression networks, we developed the catRAPID express web server.

Conclusions

Our analysis sheds light on the role of RNA-binding proteins in regulating proliferation and differentiation processes, and we provide a data exploration tool to aid future experimental studies.

Searching for repetitions in biological networks: methods, resources and tools

We present here a compact overview of the data, models and methods proposed for the analysis of biological networks based on the search for significant repetitions. In particular, we concentrate on three problems widely studied in the literature: 'network alignment', 'network querying' and 'network motif extraction'. We provide (i) details of the experimental techniques used to obtain the main types of interaction data, (ii) descriptions of the models and approaches introduced to solve such problems and (iii) pointers to both the available databases and software tools. The intent is to lay out a useful roadmap for identifying suitable strategies to analyse cellular data, possibly based on the joint use of different interaction data types or analysis techniques.

tRanslatome: an R/Bioconductor package to portray translational control

Summary: High-throughput technologies have led to an explosion of genomic data available for automated analysis. The consequent possibility to simultaneously sample multiple layers of variation along the gene expression flow requires computational methods integrating raw information from different ‘-omics’. It has been recently demonstrated that translational control is a widespread phenomenon, with profound and still underestimated regulation capabilities. Although detecting changes in the levels of total messenger RNAs (mRNAs; the transcriptome), of polysomally loaded mRNAs (the translatome) and of proteins (the proteome) is experimentally feasible in a high-throughput way, the integration of these levels is still far from being robustly approached. Here we introduce tRanslatome, a new R/Bioconductor package, which is a complete platform for the simultaneous pairwise analysis of transcriptome, translatome and proteome data. The package includes most of the available statistical methods developed for the analysis of high-throughput data, allowing the parallel comparison of differentially expressed genes and the corresponding differentially enriched biological themes. Notably, it also enables the prediction of translational regulatory elements on mRNA sequences. The utility of this tool is demonstrated with two case studies.

Availability and implementation: tRanslatome is available in Bioconductor.

Contact: t.tebaldi@unitn.it

Supplementary information:Supplementary data are available at Bioinformatics online.

Postnatal dynamics of Zeb2 expression in rat brain: analysis of novel 3' UTR sequence reveals a miR-9 interacting site

ZEB2 is a transcription factor with established roles in neurogenesis but no defined function in postnatal brain despite extensive neuronal expression in telencephalic structures. Multiple, incompletely annotated transcripts derive from the Zeb2 locus; the purpose of the present study was to structurally characterize rat brain Zeb2 transcripts with respect to 3' untranslated (UTR) sequence in order to understand Zeb2 transcript regulation including possible interactions with regulatory molecules such as neuronal miRNAs. We cloned a 5054-nucleotide Zeb2 3' UTR that is included in the most abundant Zeb2 transcript in neonatal rat brain. Unique features of the distal 3' UTR region included a number of brain-specific miRNA target sites; a highly conserved miR-9 target site at 3' UTR position 4097 was selected for functional verification in transfection experiments. Parallel analysis of Zeb2 transcript, ZEB2 protein and miR-9 levels across postnatal brain cortical development revealed a significant accumulation of ZEB2 protein levels only between postnatal days P2 and P5 that was associated with an acute loss of postnatal miR-9 expression in cortex. These studies demonstrate novel features of Zeb2 gene expression in postnatal rat brain and highlight the importance of full transcript annotation for identifying the complement of potential transcript-interacting regulators.

Affinity analysis of differentially expressed genes in hepatocytes expressing HCV core genotype 1b or 3a

Chronic hepatitis C patients display many genotype-specific clinical features of HCV infection. The core proteins encoded by different genotypes dysregulate numerous sets of distinct host genes. In this study we tested the hypothesis that HCV core proteins 1b and 3a would actually act on a limited number of independent cellular players, as well as on several functionally linked gene products. Structural and functional tests identified a core set of host genes dysregulated by HCV core genotypes 1b and 3a. The core proteins of HCV genotypes 1b and 3a target specifically limited sets of functionally related gene products, which may be responsible for the variations in the clinical spectra associated with HCV infection.

Hyper conserved elements in vertebrate mRNA 3'-UTRs reveal a translational network of RNA-binding proteins controlled by HuR

Little is known regarding the post-transcriptional networks that control gene expression in eukaryotes. Additionally, we still need to understand how these networks evolve, and the relative role played in them by their sequence-dependent regulatory factors, non-coding RNAs (ncRNAs) and RNA-binding proteins (RBPs). Here, we used an approach that relied on both phylogenetic sequence sharing and conservation in the whole mapped 3'-untranslated regions (3'-UTRs) of vertebrate species to gain knowledge on core post-transcriptional networks. The identified human hyper conserved elements (HCEs) were predicted to be preferred binding sites for RBPs and not for ncRNAs, namely microRNAs and long ncRNAs. We found that the HCE map identified a well-known network that post-transcriptionally regulates histone mRNAs. We were then able to discover and experimentally confirm a translational network composed of RNA Recognition Motif (RRM)-type RBP mRNAs that are positively controlled by HuR, another RRM-type RBP. HuR shows a preference for these RBP mRNAs bound in stem-loop motifs, confirming its role as a 'regulator of regulators'. Analysis of the transcriptome-wide HCE distribution revealed a profile of prevalently small clusters separated by unconserved intercluster RNA stretches, which predicts the formation of discrete small ribonucleoprotein complexes in the 3'-UTRs.

Tuning the engine: an introduction to resources on post-transcriptional regulation of gene expression

In the last years post-transcriptional regulation (PTR) of gene expression has been increasingly recognized to be a powerful and general determinant of the quantitative changes in proteomes, and therefore a driving force for cell phenotypes. By means of networks of trans-factors on one hand, and cis-elements found primarily in untranslated regions (UTRs) of mRNA on the other hand, mRNA availability to translation and translation rates are tightly controlled and can be rapidly tuned according to the changing state of the cell. A number of dedicated resources and tools, including databases and predictive algorithms, have been proposed as bioinformatics aids for the study of this fundamental layer of gene expression regulation. Their use, however, is rendered difficult by heterogeneity and fragmentation.   This review aims to locate these resources in their proper space, classifying them according to their goals, limitations and integration capabilities and, in the end, to provide the user with an initial toolbox for the bioinformatic analysis of post-transcriptional regulation of gene expression. The accompanying website, available at www.ptrguide.org, lists all resources, provides summary and features for each one and will be regularly updated in the future.