APADB: a database for alternative polyadenylation and microRNA regulation events

Lost in translation: The 3'-UTR of IGF1R as an ancient long noncoding RNA

Background and objectives

The insulin-like growth factor (IGF) signaling system is a major arena of intragenomic conflict over embryonic growth between imprinted genes of maternal and paternal origin and the IGF type 1 receptor (IGF1R) promotes proliferation of many human cancers. The 3'-untranslated region (3'-UTR) of the mouse Igf1r mRNA is targeted by miR-675-3p derived from the imprinted H19 long noncoding RNA. We undertook a comparative sequence analysis of vertebrate IGF1R 3'-UTRs to determine the evolutionary history of miR-675 target sequences and to identify conserved features that are likely to be involved in post-transcriptional regulation of IGF1R translation.

Methodology

Sequences of IGF1R 3'-UTRs were obtained from public databases and analyzed using publicly available algorithms.

Results

A very long 3'-UTR is a conserved feature of vertebrate IGF1R mRNAs. We found that some ancient microRNAs, such as let-7 and mir-182, have predicted binding sites that are conserved between cartilaginous fish and mammals. One very conserved region is targeted by multiple, maternally expressed imprinted microRNAs that appear to have evolved more recently than the targeted sequences.

Conclusions and implications

The conserved structures we identify in the IGF1R 3'-UTR are strong candidates for regulating cell proliferation during development and carcinogenesis. These conserved structures are now targeted by multiple imprinted microRNAs. These observations emphasize the central importance of IGF signaling pathways in the mediation of intragenomic conflicts over embryonic growth and identify possible targets for therapeutic interventions in cancer.

MiRIAD update: using alternative polyadenylation, protein interaction network analysis and additional species to enhance exploration of the role of intragenic miRNAs and their host genes

MicroRNAs have established their role as potent regulators of the epigenome. Interestingly, most miRNAs are located within protein-coding genes with functional consequences that have yet to be fully investigated. MiRIAD is a database with an interactive and user-friendly online interface that has been facilitating research on intragenic miRNAs. In this article, we present a major update. First, data for five additional species (chimpanzee, rat, dog, cow and frog) were added to support the exploration of evolutionary aspects of the relationship between host genes and intragenic miRNAs. Moreover, we integrated data from two different sources to generate a comprehensive alternative polyadenylation dataset. The miRIAD interface was therefore redesigned and provides a completely new gene model representation, including an interactive visualization of the 3′ untranslated region (UTR) with alternative polyadenylation sites, corresponding signals and potential miRNA binding sites. Furthermore, we expanded on functional host gene network analysis. Although the previous version solely reported protein interactions, the update features a separate network analysis view that can either be accessed through the submission of a list of genes of interest or directly from a gene’s list of protein interactions. In addition to statistical properties of the submitted gene set, the interaction network graph is presented and miRNAs with seed site over- and underrepresentation are identified. In summary, the update of miRIAD provides novel datasets and bioinformatics resources with a significant increase in functionality to facilitate intragenic miRNA research in a user-friendly and interactive way.

Database URL:http://www.miriad-database.org

Untranslated Parts of Genes Interpreted: Making Heads or Tails of High-Throughput Transcriptomic Data via Computational Methods: Computational methods to discover and quantify isoforms with alternative untranslated regions

In this review we highlight the importance of defining the untranslated parts of transcripts, and present a number of computational approaches for the discovery and quantification of alternative transcription start and poly-adenylation events in high-throughput transcriptomic data. The fate of eukaryotic transcripts is closely linked to their untranslated regions, which are determined by the position at which transcription starts and ends at a genomic locus. Although the extent of alternative transcription starts and alternative poly-adenylation sites has been revealed by sequencing methods focused on the ends of transcripts, the application of these methods is not yet widely adopted by the community. We suggest that computational methods applied to standard high-throughput technologies are a useful, albeit less accurate, alternative to the expertise-demanding 5' and 3' sequencing and they are the only option for analysing legacy transcriptomic data. We review these methods here, focusing on technical challenges and arguing for the need to include better normalization of the data and more appropriate statistical models of the expected variation in the signal.

Dysregulation of Alternative Poly-adenylation as a Potential Player in Autism Spectrum Disorder

We present here the hypothesis that alternative poly-adenylation (APA) is dysregulated in the brains of individuals affected by Autism Spectrum Disorder (ASD), due to disruptions in the calcium signaling networks. APA, the process of selecting different poly-adenylation sites on the same gene, yielding transcripts with different-length 3′ untranslated regions (UTRs), has been documented in different tissues, stages of development and pathologic conditions. Differential use of poly-adenylation sites has been shown to regulate the function, stability, localization and translation efficiency of target RNAs. However, the role of APA remains rather unexplored in neurodevelopmental conditions. In the human brain, where transcripts have the longest 3′ UTRs and are thus likely to be under more complex post-transcriptional regulation, erratic APA could be particularly detrimental. In the context of ASD, a condition that affects individuals in markedly different ways and whose symptoms exhibit a spectrum of severity, APA dysregulation could be amplified or dampened depending on the individual and the extent of the effect on specific genes would likely vary with genetic and environmental factors. If this hypothesis is correct, dysregulated APA events might be responsible for certain aspects of the phenotypes associated with ASD. Evidence supporting our hypothesis is derived from standard RNA-seq transcriptomic data but we suggest that future experiments should focus on techniques that probe the actual poly-adenylation site (3′ sequencing). To address issues arising from the use of post-mortem tissue and low numbers of heterogeneous samples affected by confounding factors (such as the age, gender and health of the individuals), carefully controlled in vitro systems will be required to model the effect of calcium signaling dysregulation in the ASD brain.

Transcriptomic analyses of the anti-adipogenic effects of oleuropein in human mesenchymal stem cells

Extra virgin olive oil has positive effects on health. Oleuropein is a polyphenolic compound present in olive-tree leaves, fruits (olives) and olive oil. It is responsible for the relevant organoleptic and biological properties of olive oil, including antiadipogenic properties. Thus, the effects of oleuropein on the adipogenesis of human bone-marrow mesenchymal stem cells were studied by transcriptomics and differential gene-expression analyses. Oleuropein could upregulate expression of 60% of adipogenesis-repressed genes. Besides, it could activate signaling pathways such as Rho and β-catenin, maintaining cells at an undifferentiated stage. Our data suggest that mitochondrial activity is reduced by oleuropein, mostly during adipogenic differentiation. These results shed light on oleuropein activity on cells, with potential application as a "nutraceutical" for the prevention and treatment of diseases such as obesity and osteoporosis.

Transcriptomic analysis of molecular responses in Malus domestica 'M26' roots affected by apple replant disease

Gene expression studies in roots of apple replant disease affected plants suggested defense reactions towards biotic stress to occur which did not lead to adequate responses to the biotic stressors. Apple replant disease (ARD) leads to growth inhibition and fruit yield reduction in replanted populations and results in economic losses for tree nurseries and fruit producers. The etiology is not well understood on a molecular level and causal agents show a great diversity indicating that no definitive cause, which applies to the majority of cases, has been found out yet. Hence, it is pivotal to gain a better understanding of the molecular and physiological reactions of the plant when affected by ARD and later to overcome the disease, for example by developing tolerant rootstocks. For the first time, gene expression was investigated in roots of ARD affected plants employing massive analysis of cDNA ends (MACE) and RT-qPCR. In reaction to ARD, genes in secondary metabolite production as well as plant defense, regulatory and signaling genes were upregulated whereas for several genes involved in primary metabolism lower expression was detected. For internal verification of MACE data, candidate genes were tested via RT-qPCR and a strong positive correlation between both datasets was observed. Comparison of apple 'M26' roots cultivated in ARD soil or γ-irradiated ARD soil suggests that typical defense reactions towards biotic stress take place in ARD affected plants but they did not allow responding to the biotic stressors attack adequately, leading to the observed growth depressions in ARD variants.

Prediction of miRNA-mRNA Interactions Using miRGate

miRGate ( http://mirgate.bioinfo.cnio.es /) is a freely available database that contains predicted and experimentally validated microRNA-messenger RNA (miRNA-mRNA) target pairs. This resource includes novel predictions from five well-established algorithms, but recalculated from a common and comprehensive sequence dataset. It includes all 3'-UTR sequences of all known genes of the three more widely employed genomes (human, mouse, and rat), and all annotated miRNA sequences from those genomes. Besides, it also contains predictions for all genes in human targeted by miRNA viruses such as Epstein-Barr and Kaposi sarcoma-associated herpes virus.The approach intends to circumvent one of the main drawbacks in this area, as diverse sequences and gene database versions cause poor overlap among different target prediction methods even with experimentally confirmed targets. As a result, miRGate predictions have been successfully validated using functional assays in several laboratories.This chapter describes how a user can access target information via miRGate's web interface. It also shows how automatically access the database through the programmatic interface based on representational state transfer services (REST), using the application programming interface (API) available at http://mirgate.bioinfo.cnio.es/API .

Microscopic and Molecular Characterization of the Prehaustorial Resistance against Wheat Leaf Rust (Puccinia triticina) in Einkorn (Triticum monococcum)

Puccinia triticina f. sp. tritici (Eriks.), the causal agent of leaf rust, causes substantial yield losses in wheat production. In wheat many major leaf rust resistance genes have been overcome by virulent races. In contrast, the prehaustorial resistance (phr) against wheat leaf rust detected in the diploid wheat Einkorn (Triticum monoccocum var. monococcum) accession PI272560 confers race-independent resistance against isolates virulent on accessions harboring resistance genes located on the A-genome of Triticum aestivum. Phr in PI272560 leads to abortion of fungal development during the formation of haustorial mother cells and to increased hydrogen peroxide concentration in comparison to the susceptible accession 36554 (Triticum boeoticum ssp. thaoudar var. reuteri). Increased peroxidase and endochitinase activity was detected in PI272560 within 6 h after inoculation (hai). Comparative transcriptome profiling using Massive Analysis of cDNA Ends (MACE) in infected and non-infected leaves detected 14220 differentially expressed tags in PI272560 and 15472 in accession 36554. Of these 2908 and 3004, respectively, could be assigned to Gene Ontology (GO) categories of which 463 were detected in both accessions and 311 were differentially expressed between the accessions. In accordance with the concept of non-host resistance in PI272560, genes with similarity to peroxidases, chitinases, β-1,3-glucanases and other pathogenesis-related genes were up-regulated within the first 8 hai, whereas up-regulation of such genes was delayed in 36554. Moreover, a Phosphoribulokinase gene contributing to non-host resistance in rice against stripe rust was exclusively expressed in the resistant accession PI272560. Gene expression underpinned physiological and phenotypic observations at the site of infection and are in accordance with the concept of non-host resistance.

Transcriptomic Analyses of Adipocyte Differentiation From Human Mesenchymal Stromal-Cells (MSC)

Adipogenesis is a physiological process required for fat-tissue development, mainly involved in regulating the organism energetic-state. Abnormal distribution-changes and dysfunctions in such tissue are associated to different pathologies. Adipocytes are generated from progenitor cells, via a complex differentiating process not yet well understood. Therefore, we investigated differential mRNA and miRNA expression patterns of human mesenchymal stromal-cells (MSC) induced and not induced to differentiate into adipocytes by next (second)-generation sequencing. A total of 2,866 differentially expressed genes (101 encoding miRNA) were identified, with 705 (46 encoding miRNA) being upregulated in adipogenesis. They were related to different pathways, including PPARG, lipid, carbohydrate and energy metabolism, redox, membrane-organelle biosynthesis, and endocrine system. Downregulated genes were related to extracellular matrix and cell migration, proliferation, and differentiation. Analyses of mRNA-miRNA interaction showed that repressed miRNA-encoding genes can act downregulating PPARG-related genes; mostly the PPARG activator (PPARGC1A). Induced miRNA-encoding genes regulate downregulated genes related to TGFB1. These results shed new light to understand adipose-tissue differentiation and physiology, increasing our knowledge about pathologies like obesity, type-2 diabetes and osteoporosis. J. Cell. Physiol. 232: 771-784, 2017. © 2016 Wiley Periodicals, Inc.

m(6)A-LAIC-seq reveals the census and complexity of the m(6)A epitranscriptome

N(6)-Methyladenosine (m(6)A) is a widespread, reversible chemical modification of RNA molecules, implicated in many aspects of RNA metabolism. Little quantitative information exists as to either how many transcript copies of particular genes are m(6)A modified ('m(6)A levels') or the relationship of m(6)A modification(s) to alternative RNA isoforms. To deconvolute the m(6)A epitranscriptome, we developed m(6)A-level and isoform-characterization sequencing (m(6)A-LAIC-seq). We found that cells exhibit a broad range of nonstoichiometric m(6)A levels with cell-type specificity. At the level of isoform characterization, we discovered widespread differences in the use of tandem alternative polyadenylation (APA) sites by methylated and nonmethylated transcript isoforms of individual genes. Strikingly, there is a strong bias for methylated transcripts to be coupled with proximal APA sites, resulting in shortened 3' untranslated regions, while nonmethylated transcript isoforms tend to use distal APA sites. m(6)A-LAIC-seq yields a new perspective on transcriptome complexity and links APA usage to m(6)A modifications.

PlantAPA: A Portal for Visualization and Analysis of Alternative Polyadenylation in Plants

Alternative polyadenylation (APA) is an important layer of gene regulation that produces mRNAs that have different 3' ends and/or encode diverse protein isoforms. Up to 70% of annotated genes in plants undergo APA. Increasing numbers of poly(A) sites collected in various plant species demand new methods and tools to access and mine these data. We have created an open-access web service called PlantAPA (http://bmi.xmu.edu.cn/plantapa) to visualize and analyze genome-wide poly(A) sites in plants. PlantAPA provides various interactive and dynamic graphics and seamlessly integrates a genome browser that can profile heterogeneous cleavage sites and quantify expression patterns of poly(A) sites across different conditions. Particularly, through PlantAPA, users can analyze poly(A) sites in extended 3' UTR regions, intergenic regions, and ambiguous regions owing to alternative transcription or RNA processing. In addition, it also provides tools for analyzing poly(A) site selections, 3' UTR lengthening or shortening, non-canonical APA site switching, and differential gene expression between conditions, making it more powerful for the study of APA-mediated gene expression regulation. More importantly, PlantAPA offers a bioinformatics pipeline that allows users to upload their own short reads or ESTs for poly(A) site extraction, enabling users to further explore poly(A) site selection using stored PlantAPA poly(A) sites together with their own poly(A) site datasets. To date, PlantAPA hosts the largest database of APA sites in plants, including Oryza sativa, Arabidopsis thaliana, Medicago truncatula, and Chlamydomonas reinhardtii. As a user-friendly web service, PlantAPA will be a valuable addition to the community of biologists studying APA mechanisms and gene expression regulation in plants.

Introduction to Bioinformatics Resources for Post-transcriptional Regulation of Gene Expression

Untranslated regions (UTRs) and, to a lesser extent, coding sequences of mRNAs are involved in defining the fate of the mature transcripts through the modulation of three primary control processes, mRNA localization, degradation and translation; the action of trans-factors such as RNA-binding proteins (RBPs) and noncoding RNAs (ncRNAs) combined with the presence of defined sequence and structural cis-elements ultimately determines translation levels. Identifying functional regions in UTRs and uncovering post-transcriptional regulators acting upon these regions is thus of paramount importance to understand the spectrum of regulatory possibilities for any given mRNA. This tasks can now be approached computationally, to reduce the space of testable hypotheses and to drive experimental validation.This chapter focuses on presenting databases and tools allowing to study the various aspects of post-transcriptional regulation, including motif search (sequence and secondary structure), prediction of regulatory networks (e.g., RBP and ncRNA binding sites), profiling of the mRNAs translational state, and other aspects of this level of gene expression regulation. Two analysis pipelines are also presented as practical examples of how the described tools could be integrated and effectively employed.

Principles of miRNA-mRNA interactions: beyond sequence complementarity

MicroRNAs (miRNAs) are small non-coding RNAs that post-transcriptionally regulate gene expression by altering the translation efficiency and/or stability of targeted mRNAs. In vertebrates, more than 50% of all protein-coding RNAs are assumed to be subject to miRNA-mediated control, but current high-throughput methods that reliably measure miRNA-mRNA interactions either require prior knowledge of target mRNAs or elaborate preparation procedures. Consequently, experimentally validated interactions are relatively rare. Furthermore, in silico prediction based on sequence complementarity of miRNAs and their corresponding target sites suffers from extremely high false positive rates. Apparently, sequence complementarity alone is often insufficient to reflect the complex post-transcriptional regulation of mRNAs by miRNAs, which is especially true for animals. Therefore, combined analysis of small non-coding and protein-coding RNAs is indispensable to better understand and predict the complex dynamics of miRNA-regulated gene expression. Single-nucleotide polymorphisms (SNPs) and alternative polyadenylation (APA) can affect miRNA binding of a given transcript from different individuals and tissues, and especially APA is currently emerging as a major factor that contributes to variations in miRNA-mRNA interplay in animals. In this review, we focus on the influence of APA and SNPs on miRNA-mediated gene regulation and discuss the computational approaches that take these mechanisms into account.

Next-generation sequencing reveals novel differentially regulated mRNAs, lncRNAs, miRNAs, sdRNAs and a piRNA in pancreatic cancer

Background

Previous studies identified microRNAs (miRNAs) and messenger RNAs with significantly different expression between normal pancreas and pancreatic cancer (PDAC) tissues. Due to technological limitations of microarrays and real-time PCR systems these studies focused on a fixed set of targets. Expression of other RNA classes such as long intergenic non-coding RNAs or sno-derived RNAs has rarely been examined in pancreatic cancer. Here, we analysed the coding and non-coding transcriptome of six PDAC and five control tissues using next-generation sequencing.

Results

Besides the confirmation of several deregulated mRNAs and miRNAs, miRNAs without previous implication in PDAC were detected: miR-802, miR-2114 or miR-561. SnoRNA-derived RNAs (e.g. sno-HBII-296B) and piR-017061, a piwi-interacting RNA, were found to be differentially expressed between PDAC and control tissues. In silico target analysis of miR-802 revealed potential binding sites in the 3′ UTR of TCF4, encoding a transcription factor that controls Wnt signalling genes. Overexpression of miR-802 in MiaPaCa pancreatic cancer cells reduced TCF4 protein levels. Using Massive Analysis of cDNA Ends (MACE) we identified differential expression of 43 lincRNAs, long intergenic non-coding RNAs, e.g. LINC00261 and LINC00152 as well as several natural antisense transcripts like HNF1A-AS1 and AFAP1-AS1. Differential expression was confirmed by qPCR on the mRNA/miRNA/lincRNA level and by immunohistochemistry on the protein level.

Conclusions

Here, we report a novel lncRNA, sncRNA and mRNA signature of PDAC. In silico prediction of ncRNA targets allowed for assigning potential functions to differentially regulated RNAs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12943-015-0358-5) contains supplementary material, which is available to authorized users.

Dual 3'Seq using deepSuperSAGE uncovers transcriptomes of interacting Salmonella enterica Typhimurium and human host cells

Background

The interaction of eukaryotic host and prokaryotic pathogen cells is linked to specific changes in the cellular proteome, and consequently to infection-related gene expression patterns of the involved cells. To simultaneously assess the transcriptomes of both organisms during their interaction we developed dual 3’Seq, a tag-based sequencing protocol that allows for exact quantification of differentially expressed transcripts in interacting pro- and eukaryotic cells without prior fixation or physical disruption of the interaction.

Results

Human epithelial cells were infected with Salmonella enterica Typhimurium as a model system for invasion of the intestinal epithelium, and the transcriptional response of the infected host cells together with the differential expression of invading and intracellular pathogen cells was determined by dual 3’Seq coupled with the next-generation sequencing-based transcriptome profiling technique deepSuperSAGE (deep Serial Analysis of Gene Expression). Annotation to reference transcriptomes comprising the operon structure of the employed S. enterica Typhimurium strain allowed for in silico separation of the interacting cells including quantification of polycistronic RNAs. Eighty-nine percent of the known loci are found to be transcribed in prokaryotic cells prior or subsequent to infection of the host, while 75% of all protein-coding loci are represented in the polyadenylated transcriptomes of human host cells.

Conclusions

Dual 3’Seq was alternatively coupled to MACE (Massive Analysis of cDNA ends) to assess the advantages and drawbacks of a library preparation procedure that allows for sequencing of longer fragments. Additionally, the identified expression patterns of both organisms were validated by qRT-PCR using three independent biological replicates, which confirmed that RELB along with NFKB1 and NFKB2 are involved in the initial immune response of epithelial cells after infection with S. enterica Typhimurium.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1489-1) contains supplementary material, which is available to authorized users.

miRGate: a curated database of human, mouse and rat miRNA-mRNA targets

MicroRNAs (miRNAs) are small non-coding elements involved in the post-transcriptional down-regulation of gene expression through base pairing with messenger RNAs (mRNAs). Through this mechanism, several miRNA-mRNA pairs have been described as critical in the regulation of multiple cellular processes, including early embryonic development and pathological conditions. Many of these pairs (such as miR-15 b/BCL2 in apoptosis or BART-6/BCL6 in diffuse large B-cell lymphomas) were experimentally discovered and/or computationally predicted. Available tools for target prediction are usually based on sequence matching, thermodynamics and conservation, among other approaches. Nevertheless, the main issue on miRNA-mRNA pair prediction is the little overlapping results among different prediction methods, or even with experimentally validated pairs lists, despite the fact that all rely on similar principles. To circumvent this problem, we have developed miRGate, a database containing novel computational predicted miRNA-mRNA pairs that are calculated using well-established algorithms. In addition, it includes an updated and complete dataset of sequences for both miRNA and mRNAs 3'-Untranslated region from human (including human viruses), mouse and rat, as well as experimentally validated data from four well-known databases. The underlying methodology of miRGate has been successfully applied to independent datasets providing predictions that were convincingly validated by functional assays. miRGate is an open resource available at http://mirgate.bioinfo.cnio.es. For programmatic access, we have provided a representational state transfer web service application programming interface that allows accessing the database at http://mirgate.bioinfo.cnio.es/API/ Database URL: http://mirgate.bioinfo.cnio.es

APASdb: a database describing alternative poly(A) sites and selection of heterogeneous cleavage sites downstream of poly(A) signals

Increasing amounts of genes have been shown to utilize alternative polyadenylation (APA) 3'-processing sites depending on the cell and tissue type and/or physiological and pathological conditions at the time of processing, and the construction of genome-wide database regarding APA is urgently needed for better understanding poly(A) site selection and APA-directed gene expression regulation for a given biology. Here we present a web-accessible database, named APASdb (http://mosas.sysu.edu.cn/utr), which can visualize the precise map and usage quantification of different APA isoforms for all genes. The datasets are deeply profiled by the sequencing alternative polyadenylation sites (SAPAS) method capable of high-throughput sequencing 3'-ends of polyadenylated transcripts. Thus, APASdb details all the heterogeneous cleavage sites downstream of poly(A) signals, and maintains near complete coverage for APA sites, much better than the previous databases using conventional methods. Furthermore, APASdb provides the quantification of a given APA variant among transcripts with different APA sites by computing their corresponding normalized-reads, making our database more useful. In addition, APASdb supports URL-based retrieval, browsing and display of exon-intron structure, poly(A) signals, poly(A) sites location and usage reads, and 3'-untranslated regions (3'-UTRs). Currently, APASdb involves APA in various biological processes and diseases in human, mouse and zebrafish.

Kleat: cleavage site analysis of transcriptomes

In eukaryotic cells, alternative cleavage of 3' untranslated regions (UTRs) can affect transcript stability, transport and translation. For polyadenylated (poly(A)) transcripts, cleavage sites can be characterized with short-read sequencing using specialized library construction methods. However, for large-scale cohort studies as well as for clinical sequencing applications, it is desirable to characterize such events using RNA-seq data, as the latter are already widely applied to identify other relevant information, such as mutations, alternative splicing and chimeric transcripts. Here we describe KLEAT, an analysis tool that uses de novo assembly of RNA-seq data to characterize cleavage sites on 3' UTRs. We demonstrate the performance of KLEAT on three cell line RNA-seq libraries constructed and sequenced by the ENCODE project, and assembled using Trans-ABySS. Validating the KLEAT predictions with matched ENCODE RNA-seq and RNA-PET libraries, we show that the tool has over 90% positive predictive value when there are at least three RNA-seq reads supporting a poly(A) tail and requiring at least three RNA-PET reads mapping within 100 nucleotides as validation. We also compare the performance of KLEAT with other popular RNA-seq analysis pipelines that reconstruct 3' UTR ends, and show that it performs favourably, based on an ROC-like curve.