SyntDB: defining orthologues of human long noncoding RNAs across primates.

SyntDB (http://syntdb.amu.edu.pl/) is a collection of data on long noncoding RNAs (lncRNAs) and their evolutionary relationships in twelve primate species, including humans. This is the first database dedicated to primate lncRNAs, thousands of which are uniquely stored in SyntDB. The lncRNAs were predicted with our computational pipeline using publicly available RNA-Seq data spanning diverse tissues and organs. Most of the species included in SyntDB still lack lncRNA annotations in public resources. In addition to providing users with unique sets of lncRNAs and their characteristics, SyntDB provides data on orthology relationships between the lncRNAs of humans and other primates, which are not available on this scale elsewhere. Keeping in mind that only a small fraction of currently known human lncRNAs have been functionally characterized and that lncRNA conservation is frequently used to identify the most relevant lncRNAs for functional studies, we believe that SyntDB will contribute to ongoing research aimed at deciphering the biological roles of lncRNAs.

lncEvo: automated identification and conservation study of long noncoding RNAs.

BACKGROUND: Long noncoding RNAs represent a large class of transcripts with two common features: they exceed an arbitrary length threshold of 200 nt and are assumed to not encode proteins. Although a growing body of evidence indicates that the vast majority of lncRNAs are potentially nonfunctional, hundreds of them have already been revealed to perform essential gene regulatory functions or to be linked to a number of cellular processes, including those associated with the etiology of human diseases. To better understand the biology of lncRNAs, it is essential to perform a more in-depth study of their evolution. In contrast to protein-encoding transcripts, however, they do not show the strong sequence conservation that usually results from purifying selection; therefore, software that is typically used to resolve the evolutionary relationships of protein-encoding genes and transcripts is not applicable to the study of lncRNAs. RESULTS: To tackle this issue, we developed lncEvo, a computational pipeline that consists of three modules: (1) transcriptome assembly from RNA-Seq data, (2) prediction of lncRNAs, and (3) conservation study-a genome-wide comparison of lncRNA transcriptomes between two species of interest, including search for orthologs. Importantly, one can choose to apply lncEvo solely for transcriptome assembly or lncRNA prediction, without calling the conservation-related part. CONCLUSIONS: lncEvo is an all-in-one tool built with the Nextflow framework, utilizing state-of-the-art software and algorithms with customizable trade-offs between speed and sensitivity, ease of use and built-in reporting functionalities. The source code of the pipeline is freely available for academic and nonacademic use under the MIT license at https://gitlab.com/spirit678/lncrna_conservation_nf .

A chromatin-associated splicing isoform of OIP5-AS1 acts in cis to regulate the OIP5 oncogene.

A large portion of the human genome is transcribed into long noncoding RNAs that can range from 200 nucleotides to several kilobases in length. The number of identified lncRNAs is still growing, but only a handful of them have been functionally characterized. However, it is known that the functions of lncRNAs are closely related to their subcellular localization. Cytoplasmic lncRNAs can regulate mRNA stability, affect translation and act as miRNA sponges, while nuclear-retained long noncoding RNAs have been reported to be involved in transcriptional control, chromosome scaffolding, modulation of alternative splicing and chromatin remodelling. Through these processes, lncRNAs have diverse regulatory roles in cell biology and diseases. OIP5-AS1 (also known as Cyrano), a poorly characterized lncRNA expressed antisense to the OIP5 oncogene, is deregulated in multiple cancers. We showed that one of the OIP5-AS1 splicing forms (ENST00000501665.2) is retained in the cell nucleus where it associates with chromatin, thus narrowing down the spectrum of its possible mechanisms of action. Its knockdown with antisense LNA gapmeRs led to inhibited expression of a sense partner, OIP5, strongly suggesting a functional coupling between OIP5 and ENST00000501665.2. A subsequent bioinformatics analysis followed by RAP-MS and RNA Immunoprecipitation experiments suggested its possible mode of action; in particular, we found that ENST00000501665.2 directly binds to a number of nuclear proteins, including SMARCA4, a component of the SWI/SNF chromatin remodelling complex, whose binding motif is located in the promoter of the OIP5 oncogene.

Identification of Novel miRNAs and Their Target Genes in the Response to Abscisic Acid in Arabidopsis.

miRNAs are involved in various biological processes, including adaptive responses to abiotic stress. To understand the role of miRNAs in the response to ABA, ABA-responsive miRNAs were identified by small RNA sequencing in wild-type Arabidopsis, as well as in abi1td, mkkk17, and mkkk18 mutants. We identified 10 novel miRNAs in WT after ABA treatment, while in abi1td, mkkk17, and mkkk18 mutants, three, seven, and nine known miRNAs, respectively, were differentially expressed after ABA treatment. One novel miRNA (miRn-8) was differentially expressed in the mkkk17 mutant. Potential target genes of the miRNA panel were identified using psRNATarget. Sequencing results were validated by quantitative RT-PCR of several known and novel miRNAs in all genotypes. Of the predicted targets of novel miRNAs, seven target genes of six novel miRNAs were further validated by 5' RLM-RACE. Gene ontology analyses showed the potential target genes of ABA-responsive known and novel miRNAs to be involved in diverse cellular processes in plants, including development and stomatal movement. These outcomes suggest that a number of the identified miRNAs have crucial roles in plant responses to environmental stress, as well as in plant development, and might have common regulatory roles in the core ABA signaling pathway.

Waterlogging-Stress-Responsive LncRNAs, Their Regulatory Relationships with miRNAs and Target Genes in Cucumber (Cucumis sativus L.).

Low oxygen level is a phenomenon often occurring during the cucumber cultivation period. Genes involved in adaptations to stress can be regulated by non-coding RNA. The aim was the identification of long non-coding RNAs (lncRNAs) involved in the response to long-term waterlogging stress in two cucumber haploid lines, i.e., DH2 (waterlogging tolerant-WL-T) and DH4 (waterlogging sensitive-WL-S). Plants, at the juvenile stage, were waterlogged for 7 days (non-primed, 1xH), and after a 14-day recovery period, plants were stressed again for another 7 days (primed, 2xH). Roots were collected for high-throughput RNA sequencing. Implementation of the bioinformatic pipeline made it possible to determine specific lncRNAs for non-primed and primed plants of both accessions, highlighting differential responses to hypoxia stress. In total, 3738 lncRNA molecules were identified. The highest number (1476) of unique lncRNAs was determined for non-primed WL-S plants. Seventy-one lncRNAs were depicted as potentially being involved in acquiring tolerance to hypoxia in cucumber. Understanding the mechanism of gene regulation under long-term waterlogging by lncRNAs and their interactions with miRNAs provides sufficient information in terms of adaptation to the oxygen deprivation in cucumber. To the best of our knowledge, this is the first report concerning the role of lncRNAs in the regulation of long-term waterlogging tolerance by priming application in cucumber.

miRNEST database: an integrative approach in microRNA search and annotation.

Despite accumulating data on animal and plant microRNAs and their functions, existing public miRNA resources usually collect miRNAs from a very limited number of species. A lot of microRNAs, including those from model organisms, remain undiscovered. As a result there is a continuous need to search for new microRNAs. We present miRNEST (http://mirnest.amu.edu.pl), a comprehensive database of animal, plant and virus microRNAs. The core part of the database is built from our miRNA predictions conducted on Expressed Sequence Tags of 225 animal and 202 plant species. The miRNA search was performed based on sequence similarity and as many as 10,004 miRNA candidates in 221 animal and 199 plant species were discovered. Out of them only 299 have already been deposited in miRBase. Additionally, miRNEST has been integrated with external miRNA data from literature and 13 databases, which includes miRNA sequences, small RNA sequencing data, expression, polymorphisms and targets data as well as links to external miRNA resources, whenever applicable. All this makes miRNEST a considerable miRNA resource in a sense of number of species (544) that integrates a scattered miRNA data into a uniform format with a user-friendly web interface.

HuntMi: an efficient and taxon-specific approach in pre-miRNA identification.

BACKGROUND: Machine learning techniques are known to be a powerful way of distinguishing microRNA hairpins from pseudo hairpins and have been applied in a number of recognised miRNA search tools. However, many current methods based on machine learning suffer from some drawbacks, including not addressing the class imbalance problem properly. It may lead to overlearning the majority class and/or incorrect assessment of classification performance. Moreover, those tools are effective for a narrow range of species, usually the model ones. This study aims at improving performance of miRNA classification procedure, extending its usability and reducing computational time. RESULTS: We present HuntMi, a stand-alone machine learning miRNA classification tool. We developed a novel method of dealing with the class imbalance problem called ROC-select, which is based on thresholding score function produced by traditional classifiers. We also introduced new features to the data representation. Several classification algorithms in combination with ROC-select were tested and random forest was selected for the best balance between sensitivity and specificity. Reliable assessment of classification performance is guaranteed by using large, strongly imbalanced, and taxon-specific datasets in 10-fold cross-validation procedure. As a result, HuntMi achieves a considerably better performance than any other miRNA classification tool and can be applied in miRNA search experiments in a wide range of species. CONCLUSIONS: Our results indicate that HuntMi represents an effective and flexible tool for identification of new microRNAs in animals, plants and viruses. ROC-select strategy proves to be superior to other methods of dealing with class imbalance problem and can possibly be used in other machine learning classification tasks. The HuntMi software as well as datasets used in the research are freely available at http://lemur.amu.edu.pl/share/HuntMi/.

ERISdb: a database of plant splice sites and splicing signals.

Splicing is one of the major contributors to observed spatiotemporal diversification of transcripts and proteins in metazoans. There are numerous factors that affect the process, but splice sites themselves along with the adjacent splicing signals are critical here. Unfortunately, there is still little known about splicing in plants and, consequently, further research in some fields of plant molecular biology will encounter difficulties. Keeping this in mind, we performed a large-scale analysis of splice sites in eight plant species, using novel algorithms and tools developed by us. The analyses included identification of orthologous splice sites, polypyrimidine tracts and branch sites. Additionally we identified putative intronic and exonic cis-regulatory motifs, U12 introns as well as splice sites in 45 microRNA genes in five plant species. We also provide experimental evidence for plant splice sites in the form of expressed sequence tag and RNA-Seq data. All the data are stored in a novel database called ERISdb and are freely available at http://lemur.amu.edu.pl/share/ERISdb/.

The Bro1-like domain-containing protein, AtBro1, modulates growth and abiotic stress responses in Arabidopsis.

Abscisic acid (ABA) affects plant physiology by altering gene expression, enabling plants to adapt to a wide range of environments. Plants have evolved protective mechanisms to allow seed germination in harsh conditions. Here, we explore a subset of these mechanisms involving the AtBro1 gene, which encodes one of a small family of poorly characterised Bro1-like domain-containing proteins, in Arabidopsis thaliana plants subjected to multiple abiotic stresses. AtBro1 transcripts were upregulated by salt, ABA and mannitol stress, while AtBro1-overexpression lines demonstrated robust tolerance to drought and salt stress. Furthermore, we found that ABA elicits stress-resistance responses in loss-of-function bro1-1 mutant plants and AtBro1 regulates drought resistance in Arabidopsis. When the AtBro1 promoter was fused to the ß-glucuronidase (GUS) gene and introduced into plants, GUS was expressed mainly in rosette leaves and floral clusters, especially in anthers. Using a construct expressing an AtBro1-GFP fusion protein, AtBro1 was found to be localized in the plasma membrane in Arabidopsis protoplasts. A broad RNA-sequencing analysis revealed specific quantitative differences in the early transcriptional responses to ABA treatment between wild-type and loss-of-function bro1-1 mutant plants, suggesting that ABA stimulates stress-resistance responses via AtBro1. Additionally, transcripts levels of MOP9.5, MRD1, HEI10, and MIOX4 were altered in bro1-1 plants exposed to different stress conditions. Collectively, our results show that AtBro1 plays a significant role in the regulation of the plant transcriptional response to ABA and the induction of resistance responses to abiotic stress.

[MicroRNA databases]. / Bazy danych mikroRNA.

microRNAs (miRNAs) are small RNAs that play key roles in regulation of cellular processes and therefore could largely contribute to solving many problems in medicine, biotechnology, and other biological sciences. As a result, the numbers of research projects and publications on miRNAs are constantly growing, which is accompanied by increasing amounts of new data and databases need to be created for data storage. There are 51 dedicated miRNA databases at the moment, what make it quite difficult for the users to find relevant data. Moreover, such problems as insufficient documentation, low quality of data or flaws in the graphical interface make the things even worse. However, there are positive signs, including standardization of database interfaces, a tendency to create integrated systems that collect data from a number of databases and present it in a uniform format, and emergence of systems for automated data search and download.

Emerging Roles of Long Noncoding RNAs in Breast Cancer Epigenetics and Epitranscriptomics.

Breast carcinogenesis is a multistep process that involves both genetic and epigenetic changes. Epigenetics refers to reversible changes in gene expression that are not accompanied by changes in gene sequence. In breast cancer (BC), dysregulated epigenetic changes, such as DNA methylation and histone modifications, are accompanied by epitranscriptomic changes, in particular adenine to inosine modifications within RNA molecules. Factors that trigger these phenomena are largely unknown, but there is evidence for widespread participation of long noncoding RNAs (lncRNAs) that already have been linked to virtually any aspect of BC biology, making them promising biomarkers and therapeutic targets in BC patients. Here, we provide a systematic review of known and possible roles of lncRNAs in epigenetic and epitranscriptomic processes, along with methods and tools to study them, followed by a brief overview of current challenges regarding the use of lncRNAs in medical applications.

Comparative genomics in the search for conserved long noncoding RNAs.

Long noncoding RNAs (lncRNAs) have emerged as prominent regulators of gene expression in eukaryotes. The identification of lncRNA orthologs is essential in efforts to decipher their roles across model organisms, as homologous genes tend to have similar molecular and biological functions. The relatively high sequence plasticity of lncRNA genes compared with protein-coding genes, makes the identification of their orthologs a challenging task. This is why comparative genomics of lncRNAs requires the development of specific and, sometimes, complex approaches. Here, we briefly review current advancements and challenges associated with four levels of lncRNA conservation: genomic sequences, splicing signals, secondary structures and syntenic transcription.

Long-Term Waterlogging as Factor Contributing to Hypoxia Stress Tolerance Enhancement in Cucumber: Comparative Transcriptome Analysis of Waterlogging Sensitive and Tolerant Accessions.

Waterlogging (WL), excess water in the soil, is a phenomenon often occurring during plant cultivation causing low oxygen levels (hypoxia) in the soil. The aim of this study was to identify candidate genes involved in long-term waterlogging tolerance in cucumber using RNA sequencing. Here, we also determined how waterlogging pre-treatment (priming) influenced long-term memory in WL tolerant (WL-T) and WL sensitive (WL-S) i.e., DH2 and DH4 accessions, respectively. This work uncovered various differentially expressed genes (DEGs) activated in the long-term recovery in both accessions. De novo assembly generated 36,712 transcripts with an average length of 2236 bp. The results revealed that long-term waterlogging had divergent impacts on gene expression in WL-T DH2 and WL-S DH4 cucumber accessions: after 7 days of waterlogging, more DEGs in comparison to control conditions were identified in WL-S DH4 (8927) than in WL-T DH2 (5957). Additionally, 11,619 and 5007 DEGs were identified after a second waterlogging treatment in the WL-S and WL-T accessions, respectively. We identified genes associated with WL in cucumber that were especially related to enhanced glycolysis, adventitious roots development, and amino acid metabolism. qRT-PCR assay for hypoxia marker genes i.e., alcohol dehydrogenase (adh), 1-aminocyclopropane-1-carboxylate oxidase (aco) and long chain acyl-CoA synthetase 6 (lacs6) confirmed differences in response to waterlogging stress between sensitive and tolerant cucumbers and effectiveness of priming to enhance stress tolerance.

MEF2C shapes the microtranscriptome during differentiation of skeletal muscles.

Myocyte enhancer factor 2C (MEF2C) is a transcription factor that regulates heart and skeletal muscle differentiation and growth. Several protein-encoding genes were identified as targets of this factor; however, little is known about its contribution to the microtranscriptome composition and dynamics in myogenic programs. In this report, we aimed to address this question. Deep sequencing of small RNAs of human muscle cells revealed a set of microRNAs (miRNAs), including several muscle-specific miRNAs, that are sensitive to MEF2C depletion. As expected, in cells with knockdown of MEF2C, we found mostly downregulated miRNAs; nevertheless, as much as one-third of altered miRNAs were upregulated. The majority of these changes are driven by transcription efficiency. Moreover, we found that MEF2C affects nontemplated 3'-end nucleotide addition of miRNAs, mainly oligouridylation. The rate of these modifications is associated with the level of TUT4 which mediates RNA 3'-uridylation. Finally, we found that a quarter of miRNAs which significantly changed upon differentiation of human skeletal myoblasts is inversely altered in MEF2C deficient cells. We concluded that MEF2C is an essential factor regulating both the quantity and quality of the microtranscriptome, leaving an imprint on the stability and perhaps specificity of many miRNAs during the differentiation of muscle cells.

Complex Analysis of Retroposed Genes' Contribution to Human Genome, Proteome and Transcriptome.

Gene duplication is a major driver of organismal evolution. One of the main mechanisms of gene duplications is retroposition, a process in which mRNA is first transcribed into DNA and then reintegrated into the genome. Most gene retrocopies are depleted of the regulatory regions. Nevertheless, examples of functional retrogenes are rapidly increasing. These functions come from the gain of new spatio-temporal expression patterns, imposed by the content of the genomic sequence surrounding inserted cDNA and/or by selectively advantageous mutations, which may lead to the switch from protein coding to regulatory RNA. As recent studies have shown, these genes may lead to new protein domain formation through fusion with other genes, new regulatory RNAs or other regulatory elements. We utilized existing data from high-throughput technologies to create a complex description of retrogenes functionality. Our analysis led to the identification of human retroposed genes that substantially contributed to transcriptome and proteome. These retrocopies demonstrated the potential to encode proteins or short peptides, act as cis- and trans- Natural Antisense Transcripts (NATs), regulate their progenitors' expression by competing for the same microRNAs, and provide a sequence to lncRNA and novel exons to existing protein-coding genes. Our study also revealed that retrocopies, similarly to retrotransposons, may act as recombination hot spots. To our best knowledge this is the first complex analysis of these functions of retrocopies.

Towards a deeper annotation of human lncRNAs.

A substantial fraction of the human transcriptome is composed of the so-called long noncoding RNAs (lncRNAs), yet the available catalogs of known lncRNAs are far from complete. Moreover, functional studies of these RNAs are challenged by several factors, such as their tissue-specific expression and functional heterogeneity, resulting in only ca. 1% of them being well characterized. Here, we describe a set of 41,400 novel lncRNAs discovered with RNA-Seq data from 1463 samples encompassing diverse tissues and cell lines. We utilized publicly available transcriptomic and genomic data to provide their characteristics, such as tissue specificity, cellular abundance, polyA status, cellular localization, evolutionary conservation and transcript stability, which allowed us to speculate on their possible biological roles. We also pinpointed 24 novel lncRNAs as candidates for breast cancer biomarkers. The results bring us closer to a comprehensive annotation of human lncRNAs, though vast amounts of further work are needed to validate the predictions and fully decipher their biology. This article is part of a Special Issue entitled: ncRNA in control of gene expression edited by Kotb Abdelmohsen.

CANTATAdb 2.0: Expanding the Collection of Plant Long Noncoding RNAs.

Long non-coding RNAs (lncRNAs) are a class of potent regulators of gene expression that are found in a wide array of eukaryotes; however, our knowledge about these molecules in plants is very limited. In particular, a number of plant species with important roles in biotechnology, agriculture and basic research still lack comprehensively identified and annotated sets of lncRNAs. To address these shortcomings, we previously created a database of lncRNAs in 10 model species, called CANTATAdb, and now we are expanding this online resource to encompass 39 species, including three algae. The lncRNAs were identified computationally using publicly available RNA sequencing (RNA-Seq) data. Expression values, coding potential calculations and other types of information were used to provide annotations for the identified lncRNAs. The data are freely available for searching, browsing and downloading from an online database called CANTATAdb 2.0 ( http://cantata.amu.edu.pl , http://yeti.amu.edu.pl/CANTATA/ ).

Natural antisense transcripts in diseases: From modes of action to targeted therapies.

Antisense transcription is a widespread phenomenon in mammalian genomes, leading to production of RNAs molecules referred to as natural antisense transcripts (NATs). NATs apply diverse transcriptional and post-transcriptional regulatory mechanisms to carry out a wide variety of biological roles that are important for the normal functioning of living cells, but their dysfunctions can be associated with human diseases. In this review, we attempt to provide a molecular basis for the involvement of NATs in the etiology of human disorders such as cancers and neurodegenerative and cardiovascular diseases. We also discuss the pros and cons of oligonucleotide-based therapies targeted against NATs, and we comment on state-of-the-art progress in this promising area of clinical research. WIREs RNA 2018, 9:e1461. doi: 10.1002/wrna.1461 This article is categorized under: RNA in Disease and Development > RNA in Disease Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Interactions with Proteins and Other Molecules > Small Molecule-RNA Interactions.

lncRNA-RNA Interactions across the Human Transcriptome.

Long non-coding RNAs (lncRNAs) represent a numerous class of non-protein coding transcripts longer than 200 nucleotides. There is possibility that a fraction of lncRNAs are not functional and represent mere transcriptional noise but a growing body of evidence shows they are engaged in a plethora of molecular functions and contribute considerably to the observed diversification of eukaryotic transcriptomes and proteomes. Still, however, only ca. 1% of lncRNAs have well established functions and much remains to be done towards decipherment of their biological roles. One of the least studied aspects of lncRNAs biology is their engagement in gene expression regulation through RNA-RNA interactions. By hybridizing with mate RNA molecules, lncRNAs could potentially participate in modulation of pre-mRNA splicing, RNA editing, mRNA stability control, translation activation, or abrogation of miRNA-induced repression. Here, we implemented a similarity-search based method for transcriptome-wide identification of RNA-RNA interactions, which enabled us to find 18,871,097 lncRNA-RNA base-pairings in human. Further analyses showed that the interactions could be involved in processing, stability control and functions of 57,303 transcripts. An extensive use of RNA-Seq data provided support for approximately one third of the interactions, at least in terms of the two RNA components being co-expressed. The results suggest that lncRNA-RNA interactions are broadly used to regulate and diversify the human transcriptome.

Retroposition as a source of antisense long non-coding RNAs with possible regulatory functions.

Long non-coding RNAs (lncRNAs) are a class of intensely studied, yet enigmatic molecules that make up a substantial portion of the human transcriptome. In this work, we link the origins and functions of some lncRNAs to retroposition, a process resulting in the creation of intronless copies (retrocopies) of the so-called parental genes. We found 35 human retrocopies transcribed in antisense and giving rise to 58 lncRNA transcripts. These lncRNAs share sequence similarity with the corresponding parental genes but in the sense/antisense orientation, meaning they have the potential to interact with each other and to form RNA:RNA duplexes. We took a closer look at these duplexes and found that 10 of the lncRNAs might regulate parental gene expression and processing at the pre-mRNA and mRNA levels. Further analysis of the co-expression and expression correlation provided support for the existence of functional coupling between lncRNAs and their mate parental gene transcripts.