RAIN: RNA-protein Association and Interaction Networks

Noncoding RNAs improve the predictive power of network medicine

Network medicine has improved the mechanistic understanding of disease, offering quantitative insights into disease mechanisms, comorbidities, and novel diagnostic tools and therapeutic treatments. Yet, most network-based approaches rely on a comprehensive map of protein-protein interactions (PPI), ignoring interactions mediated by noncoding RNAs (ncRNAs). Here, we systematically combine experimentally confirmed binding interactions mediated by ncRNA with PPI, constructing a comprehensive network of all physical interactions in the human cell. We find that the inclusion of ncRNA expands the number of genes in the interactome by 46% and the number of interactions by 107%, significantly enhancing our ability to identify disease modules. Indeed, we find that 132 diseases lacked a statistically significant disease module in the protein-based interactome but have a statistically significant disease module after inclusion of ncRNA-mediated interactions, making these diseases accessible to the tools of network medicine. We show that the inclusion of ncRNAs helps unveil disease-disease relationships that were not detectable before and expands our ability to predict comorbidity patterns between diseases. Taken together, we find that including noncoding interactions improves both the breath and the predictive accuracy of network medicine.

Exploring the landscape of tools and resources for the analysis of long non-coding RNAs

In recent years, research on long non-coding RNAs (lncRNAs) has gained considerable attention due to the increasing number of newly identified transcripts. Several characteristics make their functional evaluation challenging, which called for the urgent need to combine molecular biology with other disciplines, including bioinformatics. Indeed, the recent development of computational pipelines and resources has greatly facilitated both the discovery and the mechanisms of action of lncRNAs. In this review, we present a curated collection of the most recent computational resources, which have been categorized into distinct groups: databases and annotation, identification and classification, interaction prediction, and structure prediction. As the repertoire of lncRNAs and their analysis tools continues to expand over the years, standardizing the computational pipelines and improving the existing annotation of lncRNAs will be crucial to facilitate functional genomics studies.

The lncRNA expression profile signature of leukemia stem cells is altered upon PI3K/mTOR inhibition: an in vitro and in silico study

Genetic and/or epigenetic alterations in hematopoietic stem cells (HSCs) contribute to leukemia stem cell (LSC) formation. We aimed to identify alterations in the lncRNA expression profile signature of LSCs upon inhibition of PI3K/Akt/mTOR signaling, which provides selective advantages to LSCs. We also aimed to elucidate the potential interaction networks and functions of differentially expressed lncRNAs (DELs). We suppressed PI3K/Akt/mTOR signaling in LSC and HSC cell-lines by specific PI3K/mTOR dual-inhibitor (VS-5584) and confirmed the inhibition by antibody-array. We defined DELs by qRT-PCR. Increased SRA, ZEB2-AS1, antiPeg11, DLX6-AS, SNHG4, and decreased H19, PCGEM1, CAR-Intergenic-10, L1PA16, IGF2AS, and SNHG5 levels (|log2fold-change|>5) were strictly associated with PI3K/Akt/mTOR pathway inhibition in LSC. We performed in silico analyses for DELs. ZEB2-AS1 was found to be specifically expressed in normal bone marrow and predominantly lower in leukemic cell-lines. Three sub-clusters were identified for DELs and they were associated with "abnormality of multiple cell lineages in the bone marrow." DELs were most highly enriched for "glucuronidation" Reactome pathway and "ascorbate and aldarate metabolism" and "inositol phosphate metabolism" KEGG pathways. Transcription factors, MBD4, NANOG, PAX6, RELA, CEBPB, and CEBPA were predicted to be associated with the DEL profile. SRA was predicted to interact with CREB1, RARA, and PPARA. The possible DELs' targets were predicted to form six ontological groups, be highly enriched for phosphoprotein, and be involved in "PPAR signaling pathway" and "ChREBP regulation by carbohydrates and cAMP." These results will help to elucidate the roles of lncRNAs in the mechanisms that provide selective advantages to leukemia stem cells.

Functional Relationships between Long Non-Coding RNAs and Estrogen Receptor Alpha: A New Frontier in Hormone-Responsive Breast Cancer Management

In the complex and articulated machinery of the human genome, less than 2% of the transcriptome encodes for proteins, while at least 75% is actively transcribed into non-coding RNAs (ncRNAs). Among the non-coding transcripts, those ≥200 nucleotides long (lncRNAs) are receiving growing attention for their involvement in human diseases, particularly cancer. Genomic studies have revealed the multiplicity of processes, including neoplastic transformation and tumor progression, in which lncRNAs are involved by regulating gene expression at epigenetic, transcriptional, and post-transcriptional levels by mechanism(s) that still need to be clarified. In breast cancer, several lncRNAs were identified and demonstrated to have either oncogenic or tumor-suppressive roles. The functional understanding of the mechanisms of lncRNA action in this disease could represent a potential for translational applications, as these molecules may serve as novel biomarkers of clinical use and potential therapeutic targets. This review highlights the relationship between lncRNAs and the principal hallmark of the luminal breast cancer phenotype, estrogen receptor α (ERα), providing an overview of new potential ways to inhibit estrogenic signaling via this nuclear receptor toward escaping resistance to endocrine therapy.

MuCoMiD: A Multitask Graph Convolutional Learning Framework for miRNA-Disease Association Prediction

Growing evidence from recent studies implies that microRNAs or miRNAs could serve as biomarkers in various complex human diseases. Since wet-lab experiments for detecting miRNAs associated with a disease are expensive and time-consuming, machine learning techniques for miRNA-disease association prediction have attracted much attention in recent years. A big challenge in building reliable machine learning models is that of data scarcity. In particular, existing approaches trained on the available small datasets, even when combined with precalculated handcrafted input features, often suffer from bad generalization and data leakage problems. We overcome the limitations of existing works by proposing a novel multitask graph convolution-based approach, which we refer to as MuCoMiD. MuCoMiD allows automatic feature extraction while incorporating knowledge from five heterogeneous biological information sources (associations between miRNAs/diseases and protein-coding genes (PCGs), interactions between protein-coding genes, miRNA family information, and disease ontology) in a multitask setting which is a novel perspective and has not been studied before. To effectively test the generalization capability of our model, we conduct large-scale experiments on the standard benchmark datasets as well as on our proposed large independent testing sets and case studies. MuCoMiD obtains significantly higher Average Precision (AP) scores than all benchmarked models on three large independent testing sets, especially those with many new miRNAs, as well as in the detection of false positives. Thanks to its capability of learning directly from raw input information, MuCoMiD is easier to maintain and update than handcrafted feature-based methods, which would require recomputation of features every time there is a change in the original information sources (e.g., disease ontology, miRNA/disease-PCG associations, etc.). We share our code for reproducibility and future research at https://git.l3s.uni-hannover.de/dong/cmtt.

Detecting RNA-RNA interactome

The last decade has seen a robust increase in various types of novel RNA molecules and their complexity in gene regulation. RNA molecules play a critical role in cellular events by interacting with other biomolecules, including protein, DNA, and RNA. It has been established that RNA-RNA interactions play a critical role in several biological processes by regulating the biogenesis and function of RNA molecules. Interestingly, RNA-RNA interactions regulate the biogenesis of diverse RNA molecules, including mRNAs, microRNAs, tRNAs, and circRNAs, through splicing or backsplicing. Structured RNAs like rRNA, tRNA, and snRNAs achieve their functional conformation by intramolecular RNA-RNA interactions. In addition, functional consequences of many intermolecular RNA-RNA interactions have been extensively studied in the regulation of gene expression. Hence, it is essential to understand the mechanism and functions of RNA-RNA interactions in eukaryotes. Conventionally, RNA-RNA interactions have been identified through diverse biochemical methods for decades. The advent of high-throughput RNA-sequencing technologies has revolutionized the identification of global RNA-RNA interactome in cells and their importance in RNA structure and function in gene expression regulation. Although these technologies revealed tens of thousands of intramolecular and intermolecular RNA-RNA interactions, we further look forward to future unbiased and quantitative high-throughput technologies for detecting transcriptome-wide RNA-RNA interactions. With the ability to detect RNA-RNA interactome, we expect that future studies will reveal the higher-order structures of RNA molecules and multi-RNA hybrids impacting human health and diseases. This article is categorized under: RNA Methods > RNA Analyses In Vitro and In Silico RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

Variation in the co-expression profile highlights a loss of miRNA-mRNA regulation in multiple cancer types

Recent research provides insight into the ability of miRNA to regulate various pathways in several cancer types. Despite their involvement in the regulation of the mRNA via targeting the 3′UTR, there are relatively few studies examining the changes in these regulatory mechanisms specific to single cancer types or shared between different cancer types.

We analyzed samples where both miRNA and mRNA expression had been measured and performed a thorough correlation analysis on 7494 experimentally validated human miRNA-mRNA target-gene pairs in both healthy and tumoral samples.

We show how more than 90% of these miRNA-mRNA interactions show a loss of regulation in the tumoral samples compared with their healthy counterparts.

As expected, we found shared miRNA-mRNA dysregulated pairs among different tumors of the same tissue. However, anatomically different cancers also share multiple dysregulated interactions, suggesting that some cancer-related mechanisms are not tumor-specific. 2865 unique miRNA-mRNA pairs were identified across 13 cancer types, ≈ 40% of these pairs showed a loss of correlation in the tumoral samples in at least 2 out of the 13 analyzed cancers. Specifically, miR-200 family, miR-155 and miR-1 were identified, based on the computational analysis described below, as the miRNAs that potentially lose the highest number of interactions across different samples (only literature-based interactions were used for this analysis).

Moreover, the miR-34a/ALDH2 and miR-9/MTHFD2 pairs show a switch in their correlation between healthy and tumor kidney samples suggesting a possible change in the regulation exerted by the miRNAs. Interestingly, the expression of these mRNAs is also associated with the overall survival. The disruption of miRNA regulation on its target, therefore, suggests the possible involvement of these pairs in cell malignant functions.

The analysis reported here shows how the regulation of miRNA-mRNA interactions strongly differs between healthy and tumoral cells, based on the strong correlation variation between miRNA and its target that we obtained by analyzing the expression data of healthy and tumor tissue in highly reliable miRNA-target pairs. Finally, a go term enrichment analysis shows that the critical pairs identified are involved in cellular adhesion, proliferation, and migration.

Restoration of aged hematopoietic cells by their young counterparts through instructive microvesicles release

This study addresses the potential to reverse age-associated morbidity by establishing methods to restore the aged hematopoietic system. Parabiotic animal models indicated that young secretome could restore aged tissues, leading us to establish a heterochronic transwell system with aged mobilized peripheral blood (MPB), co-cultured with young MPB or umbilical cord blood (UCB) cells. Functional studies and omics approaches indicate that the miRNA cargo of microvesicles (MVs) restores the aged hematopoietic system. The in vitro findings were validated in immune deficient (NSG) mice carrying an aged hematopoietic system, improving aged hallmarks such as increased lymphoid:myeloid ratio, decreased inflammation and cellular senescence. Elevated MYC and E2F pathways, and decreased p53 were key to hematopoietic restoration. These processes require four restorative miRs that target the genes for transcription/differentiation, namely PAX and phosphatase PPMIF. These miRs when introduced in aged cells were sufficient to restore the aged hematopoietic system in NSG mice. The aged MPBs were the drivers of their own restoration, as evidenced by the changes from distinct baseline miR profiles in MPBs and UCB to comparable expressions after exposure to aged MPBs. Restorative natural killer cells eliminated dormant breast cancer cells in vivo, indicating the broad relevance of this cellular paradigm - preventing and reversing age-associated disorders such as clearance of early malignancies and enhanced responses to vaccine and infection.

OnTheFly2.0: a text-mining web application for automated biomedical entity recognition, document annotation, network and functional enrichment analysis

Extracting and processing information from documents is of great importance as lots of experimental results and findings are stored in local files. Therefore, extracting and analyzing biomedical terms from such files in an automated way is absolutely necessary. In this article, we present OnTheFly^2.0, a web application for extracting biomedical entities from individual files such as plain texts, office documents, PDF files or images. OnTheFly^2.0 can generate informative summaries in popup windows containing knowledge related to the identified terms along with links to various databases. It uses the EXTRACT tagging service to perform named entity recognition (NER) for genes/proteins, chemical compounds, organisms, tissues, environments, diseases, phenotypes and gene ontology terms. Multiple files can be analyzed, whereas identified terms such as proteins or genes can be explored through functional enrichment analysis or be associated with diseases and PubMed entries. Finally, protein-protein and protein-chemical networks can be generated with the use of STRING and STITCH services. To demonstrate its capacity for knowledge discovery, we interrogated published meta-analyses of clinical biomarkers of severe COVID-19 and uncovered inflammatory and senescence pathways that impact disease pathogenesis. OnTheFly^2.0 currently supports 197 species and is available at http://bib.fleming.gr:3838/OnTheFly/ and http://onthefly.pavlopouloslab.info.

A Survey of Current Resources to Study lncRNA-Protein Interactions

Phenotypes are driven by regulated gene expression, which in turn are mediated by complex interactions between diverse biological molecules. Protein-DNA interactions such as histone and transcription factor binding are well studied, along with RNA-RNA interactions in short RNA silencing of genes. In contrast, lncRNA-protein interaction (LPI) mechanisms are comparatively unknown, likely directed by the difficulties in studying LPI. However, LPI are emerging as key interactions in epigenetic mechanisms, playing a role in development and disease. Their importance is further highlighted by their conservation across kingdoms. Hence, interest in LPI research is increasing. We therefore review the current state of the art in lncRNA-protein interactions. We specifically surveyed recent computational methods and databases which researchers can exploit for LPI investigation. We discovered that algorithm development is heavily reliant on a few generic databases containing curated LPI information. Additionally, these databases house information at gene-level as opposed to transcript-level annotations. We show that early methods predict LPI using molecular docking, have limited scope and are slow, creating a data processing bottleneck. Recently, machine learning has become the strategy of choice in LPI prediction, likely due to the rapid growth in machine learning infrastructure and expertise. While many of these methods have notable limitations, machine learning is expected to be the basis of modern LPI prediction algorithms.

Effect of long noncoding RNAs on epithelial-mesenchymal transition in A549 cells and fibrotic human lungs

Long noncoding RNAs (LncRNAs) regulate epithelial-mesenchymal transition (EMT). EMT involves myofibroblast differentiation and pulmonary fibrosis (PF). We aimed to determine the expression profiles of HOTAIR, CARLo-5, and CD99P1 LncRNAs in EMT-mediated myofibroblast differentiation in A549 cells and fibrotic human lungs and to explain their roles. A group of A549s was stimulated with transforming growth factor β (TGF-β; 5 ng/ml) to induce EMT. The remaining A549s were incubated with 20 μM FH535 after 24 h of TGF-β treatment to inhibit EMT. A549s were collected at 0, 24, 36, and 48 h. Expressions of three LncRNAs and protein/genes related to EMT, myofibroblast differentiation, and PF were assayed by quantitative reverse-transcription polymerase chain reaction and Western blot analysis in A549s and fibrotic human lungs. The targets of three LncRNAs were investigated by bioinformatics methods. TGF-β stimulation resulted in increased expressions of three LncRNAs, ACTA2, COL1A1, SNAI1, CTNNB1, TCF4, LEF1, α-SMA, and active-β-catenin, and decreased E-cadherin at 24, 36, and 48 h in A549s. FH535 treatment regressed these alterations. But it increased HOTAIR expression at 36 h and did not increase E-cadherin at 48 h. Fibrotic human lungs were characterized by increased expressions of HOTAIR, CARLo-5, CD99P1, and miR-214, decreased expressions of miR-148b, miR-218-1, miR-7-1, and the presence of CARLo-5 and CD99P1 in HDAC1-LncRNAs coprecipitation products, but not HOTAIR. Bioinformatic analysis showed the interactions of three LncRNAs with both proteins and at least 13 microRNAs related to EMT and PF. In conclusion, HOTAIR, CARLo-5, and CD99P1 can regulate EMT-mediated myofibroblast differentiation through interacting with proteins and miRNAs associated with EMT and PF. These LncRNAs can be considered as potential targets to decrease EMT for treating PF.

An integrated workflow for biomarker development using microRNAs in extracellular vesicles for cancer precision medicine

EV-miRNAs are microRNA (miRNA) molecules encapsulated in extracellular vesicles (EVs), which play crucial roles in tumor pathogenesis, progression, and metastasis. Recent studies about EV-miRNAs have gained novel insights into cancer biology and have demonstrated a great potential to develop novel liquid biopsy assays for various applications. Notably, compared to conventional liquid biomarkers, EV-miRNAs are more advantageous in representing host-cell molecular architecture and exhibiting higher stability and specificity. Despite various available techniques for EV-miRNA separation, concentration, profiling, and data analysis, a standardized approach for EV-miRNA biomarker development is yet lacking. In this review, we performed a substantial literature review and distilled an integrated workflow encompassing important steps for EV-miRNA biomarker development, including sample collection and EV isolation, EV-miRNA extraction and quantification, high-throughput data preprocessing, biomarker prioritization and model construction, functional analysis, as well as validation. With the rapid growth of "big data", we highlight the importance of efficient mining of high-throughput data for the discovery of EV-miRNA biomarkers and integrating multiple independent datasets for in silico and experimental validations to increase the robustness and reproducibility. Furthermore, as an efficient strategy in systems biology, network inference provides insights into the regulatory mechanisms and can be used to select functionally important EV-miRNAs to refine the biomarker candidates. Despite the encouraging development in the field, a number of challenges still hinder the clinical translation. We finally summarize several common challenges in various biomarker studies and discuss potential opportunities emerging in the related fields.

Approaches to Identify and Characterise the Post-Transcriptional Roles of lncRNAs in Cancer

It is becoming increasingly evident that the non-coding genome and transcriptome exert great influence over their coding counterparts through complex molecular interactions. Among non-coding RNAs (ncRNA), long non-coding RNAs (lncRNAs) in particular present increased potential to participate in dysregulation of post-transcriptional processes through both RNA and protein interactions. Since such processes can play key roles in contributing to cancer progression, it is desirable to continue expanding the search for lncRNAs impacting cancer through post-transcriptional mechanisms. The sheer diversity of mechanisms requires diverse resources and methods that have been developed and refined over the past decade. We provide an overview of computational resources as well as proven low-to-high throughput techniques to enable identification and characterisation of lncRNAs in their complex interactive contexts. As more cancer research strategies evolve to explore the non-coding genome and transcriptome, we anticipate this will provide a valuable primer and perspective of how these technologies have matured and will continue to evolve to assist researchers in elucidating post-transcriptional roles of lncRNAs in cancer.

The lncRNA Toolkit: Databases and In Silico Tools for lncRNA Analysis

Long non-coding RNAs (lncRNAs) are a rapidly expanding field of research, with many new transcripts identified each year. However, only a small subset of lncRNAs has been characterized functionally thus far. To aid investigating the mechanisms of action by which new lncRNAs act, bioinformatic tools and databases are invaluable. Here, we review a selection of computational tools and databases for the in silico analysis of lncRNAs, including tissue-specific expression, protein coding potential, subcellular localization, structural conformation, and interaction partners. The assembled lncRNA toolkit is aimed primarily at experimental researchers as a useful starting point to guide wet-lab experiments, mainly containing multi-functional, user-friendly interfaces. With more and more new lncRNA analysis tools available, it will be essential to provide continuous updates and maintain the availability of key software in the future.

Versatile interactions and bioinformatics analysis of noncoding RNAs

Advances in RNA sequencing technologies and computational methodologies have provided a huge impetus to noncoding RNA (ncRNA) study. Once regarded as inconsequential results of transcriptional promiscuity, ncRNAs were later found to exert great roles in various aspects of biological functions. They are emerging as key players in gene regulatory networks by interacting with other biomolecules (DNA, RNA or protein). Here, we provide an overview of ncRNA repertoire and highlight recent discoveries of their versatile interactions. To better investigate the ncRNA-mediated regulation, it is necessary to make full use of innovative sequencing techniques and computational tools. We further describe a comprehensive workflow for in silico ncRNA analysis, providing up-to-date platforms, databases and tools dedicated to ncRNA identification and functional annotation.

RNAInter in 2020: RNA interactome repository with increased coverage and annotation

Research on RNA-associated interactions has exploded in recent years, and increasing numbers of studies are not limited to RNA-RNA and RNA-protein interactions but also include RNA-DNA/compound interactions. To facilitate the development of the interactome and promote understanding of the biological functions and molecular mechanisms of RNA, we updated RAID v2.0 to RNAInter (RNA Interactome Database), a repository for RNA-associated interactions that is freely accessible at http://www.rna-society.org/rnainter/ or http://www.rna-society.org/raid/. Compared to RAID v2.0, new features in RNAInter include (i) 8-fold more interaction data and 94 additional species; (ii) more definite annotations organized, including RNA editing/localization/modification/structure and homology interaction; (iii) advanced functions including fuzzy/batch search, interaction network and RNA dynamic expression and (iv) four embedded RNA interactome tools: RIscoper, IntaRNA, PRIdictor and DeepBind. Consequently, RNAInter contains >41 million RNA-associated interaction entries, involving more than 450 thousand unique molecules, including RNA, protein, DNA and compound. Overall, RNAInter provides a comprehensive RNA interactome resource for researchers and paves the way to investigate the regulatory landscape of cellular RNAs.

Identification of two microRNA nodes as potential cooperative modulators of liver metabolism

AIM

Hepatic insulin resistance is a hallmark of type 2 diabetes and non-alcoholic fatty liver disease. Dysregulation of microRNA (miRNA) expression in insulin-resistant livers might coordinate impaired hepatic metabolic function. Here, we aimed to discover miRNAs and their downstream targets involved in hepatic insulin resistance.

METHODS

We determined miRNA expression profiles by small RNA sequencing of two mouse models of impaired hepatic insulin action: high-fat diet-induced obesity and liver-specific insulin receptor knockout. Conversely, we assessed the hepatic miRNA expression profile after treatment with the antidiabetic hormone, fibroblast growth factor 21 (FGF21). Ontology analysis of predicted miRNA gene targets was performed to identify regulated gene pathways. Target enrichment analysis and miRNA mimic overexpression in vitro were used to identify unified protein targets of nodes of regulated miRNAs.

RESULTS

We identified an array of miRNA species regulated by impaired liver insulin action or after fibroblast growth factor 21 treatment. Ontology analysis of predicted miRNA gene targets identified pathways controlling hepatic energy metabolism and insulin sensitivity. We identified a node of two miRNAs downregulated in the livers of liver-specific insulin receptor knockout mice, miR-883b and miR-205, which positively regulate the expression of transcription factor zinc finger E-box-binding homeobox 1 (ZBED1). We found another node of two miRNAs upregulated in the livers of fibroblast growth factor 21-treated mice, miR-155-3p and miR-1968-5p, which canonically downregulates the caveola component, polymerase I and transcript release factor (PTRF), a gene previously implicated in hepatic energy metabolism.

CONCLUSIONS

This study identifies two nodes of coregulated miRNAs that might coordinately control hepatic energy metabolism in states of insulin resistance.

Non-coding RNA regulatory networks

It is well established that the vast majority of human RNA transcripts do not encode for proteins and that non-coding RNAs regulate cell physiology and shape cellular functions. A subset of them is involved in gene regulation at different levels, from epigenetic gene silencing to post-transcriptional regulation of mRNA stability. Notably, the aberrant expression of many non-coding RNAs has been associated with aggressive pathologies. Rapid advances in network biology indicates that the robustness of cellular processes is the result of specific properties of biological networks such as scale-free degree distribution and hierarchical modularity, suggesting that regulatory network analyses could provide new insights on gene regulation and dysfunction mechanisms. In this study we present an overview of public repositories where non-coding RNA-regulatory interactions are collected and annotated, we discuss unresolved questions for data integration and we recall existing resources to build and analyse networks.

RISE: a database of RNA interactome from sequencing experiments

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

Differentially Expressed MiRNAs and tRNA Genes Affect Host Homeostasis During Highly Pathogenic Porcine Reproductive and Respiratory Syndrome Virus Infections in Young Pigs

Background: Porcine respiratory and reproductive syndrome virus (PRRSV) is a single-stranded RNA virus member that infects pigs and causes losses to the commercial industry reaching upward of a billion dollars annually in combined direct and indirect costs. The virus can be separated into etiologies that contain multiple heterologous low and highly pathogenic strains. Recently, the United States has begun to see an increase in heterologous type 2 PRRSV strains of higher virulence (HP-PRRSV). The high pathogenicity of these strains can drastically alter host immune responses and the ability of the animal to maintain homeostasis. Because the loss of host homeostasis can denote underlying changes in gene and regulatory element expression profiles, the study aimed to examine the effect PRRSV infections has on miRNA and tRNA expression and the roles they play in host tolerance or susceptibility.

Results: Using transcriptomic analysis of whole blood taken from control and infected pigs at several time points (1, 3, 8 dpi), the analysis returned a total of 149 statistically significant (FDR ⫹ 0.15) miRNAs (n = 89) and tRNAs (n = 60) that were evaluated for possible pro- and anti-viral effects. The tRNA differential expression increased in both magnitude and count as dpi increased, with no statistically significant expression at 1 dpi, but increases at 3 and 8 dpi. The most abundant tRNA amino acid at 3 dpi was alanine, while glycine was the most abundant at 8 dpi. For the miRNAs, focus was put on upregulation that can inhibit gene expression. These results yielded candidates with potential anti- and pro-viral actions such as Ssc-miR-125b, which is predicted to limit PRRSV viral levels, and Ssc-miR-145-5p shown to cause alternative macrophage priming. The results also showed that both the tRNAs and miRNAs displayed expression patterns.

Conclusions: The results indicated that the HP-PRRSV infection affects host homeostasis through changes in miRNA and tRNA expression and their subsequent gene interactions that target and influence the function of host immune, metabolic, and structural pathways.

mirDIP 4.1-integrative database of human microRNA target predictions

MicroRNAs are important regulators of gene expression, achieved by binding to the gene to be regulated. Even with modern high-throughput technologies, it is laborious and expensive to detect all possible microRNA targets. For this reason, several computational microRNA-target prediction tools have been developed, each with its own strengths and limitations. Integration of different tools has been a successful approach to minimize the shortcomings of individual databases. Here, we present mirDIP v4.1, providing nearly 152 million human microRNA-target predictions, which were collected across 30 different resources. We also introduce an integrative score, which was statistically inferred from the obtained predictions, and was assigned to each unique microRNA-target interaction to provide a unified measure of confidence. We demonstrate that integrating predictions across multiple resources does not cumulate prediction bias toward biological processes or pathways. mirDIP v4.1 is freely available at http://ophid.utoronto.ca/mirDIP/.

miRandola 2017: a curated knowledge base of non-invasive biomarkers

miRandola (http://mirandola.iit.cnr.it/) is a database of extracellular non-coding RNAs (ncRNAs) that was initially published in 2012, foreseeing the relevance of ncRNAs as non-invasive biomarkers. An increasing amount of experimental evidence shows that ncRNAs are frequently dysregulated in diseases. Further, ncRNAs have been discovered in different extracellular forms, such as exosomes, which circulate in human body fluids. Thus, miRandola 2017 is an effort to update and collect the accumulating information on extracellular ncRNAs that is spread across scientific publications and different databases. Data are manually curated from 314 articles that describe miRNAs, long non-coding RNAs and circular RNAs. Fourteen organisms are now included in the database, and associations of ncRNAs with 25 drugs, 47 sample types and 197 diseases. miRandola also classifies extracellular RNAs based on their extracellular form: Argonaute2 protein, exosome, microvesicle, microparticle, membrane vesicle, high density lipoprotein and circulating. We also implemented a new web interface to improve the user experience.

Recent computational developments on CLIP-seq data analysis and microRNA targeting implications

Cross-Linking Immunoprecipitation associated to high-throughput sequencing (CLIP-seq) is a technique used to identify RNA directly bound to RNA-binding proteins across the entire transcriptome in cell or tissue samples. Recent technological and computational advances permit the analysis of many CLIP-seq samples simultaneously, allowing us to reveal the comprehensive network of RNA–protein interaction and to integrate it to other genome-wide analyses. Therefore, the design and quality management of the CLIP-seq analyses are of critical importance to extract clean and biological meaningful information from CLIP-seq experiments. The application of CLIP-seq technique to Argonaute 2 (Ago2) protein, the main component of the microRNA (miRNA)-induced silencing complex, reveals the direct binding sites of miRNAs, thus providing insightful information about the role played by miRNA(s). In this review, we summarize and discuss the most recent computational methods for CLIP-seq analysis, and discuss their impact on Ago2/miRNA-binding site identification and prediction with a regard toward human pathologies.

Systems Bioinformatics: increasing precision of computational diagnostics and therapeutics through network-based approaches

Systems Bioinformatics is a relatively new approach, which lies in the intersection of systems biology and classical bioinformatics. It focuses on integrating information across different levels using a bottom-up approach as in systems biology with a data-driven top-down approach as in bioinformatics. The advent of omics technologies has provided the stepping-stone for the emergence of Systems Bioinformatics. These technologies provide a spectrum of information ranging from genomics, transcriptomics and proteomics to epigenomics, pharmacogenomics, metagenomics and metabolomics. Systems Bioinformatics is the framework in which systems approaches are applied to such data, setting the level of resolution as well as the boundary of the system of interest and studying the emerging properties of the system as a whole rather than the sum of the properties derived from the system's individual components. A key approach in Systems Bioinformatics is the construction of multiple networks representing each level of the omics spectrum and their integration in a layered network that exchanges information within and between layers. Here, we provide evidence on how Systems Bioinformatics enhances computational therapeutics and diagnostics, hence paving the way to precision medicine. The aim of this review is to familiarize the reader with the emerging field of Systems Bioinformatics and to provide a comprehensive overview of its current state-of-the-art methods and technologies. Moreover, we provide examples of success stories and case studies that utilize such methods and tools to significantly advance research in the fields of systems biology and systems medicine.

Long noncoding RNA: a crosslink in biological regulatory network

Long noncoding RNAs (lncRNAs) had been defined as a novel class of functional RNAs longer than 200 nucleotides around a decade ago. It is widely acknowledged that lncRNAs play a significant role in regulation of gene expression, but the biological and molecular mechanisms are diverse and complex, and remain to be determined. Especially, the regulatory network of lncRNAs associated with other biological molecules is still a controversial matter, thus becoming a new frontier of the studies on transcriptome. Recent advance in high-throughput sequencing technologies and bioinformatics approaches may be an accelerator to lift the mysterious veil. In this review, we will outline well-known associations between lncRNAs and other biological molecules, demonstrate the diverse bioinformatics approaches applied in prediction and analysis of lncRNA interaction and perform a case study for lncRNA linc00460 to concretely decipher the lncRNA regulatory network.

Design, implementation, and operation of a rapid, robust named entity recognition web service

Most BioCreative tasks to date have focused on assessing the quality of text-mining annotations in terms of precision and recall. Interoperability, speed, and stability are, however, other important factors to consider for practical applications of text mining. For about a decade, we have run named entity recognition (NER) web services, which are designed to be efficient, implemented using a multi-threaded queueing system to robustly handle many simultaneous requests, and hosted at a supercomputer facility. To participate in this new task, we extended the existing NER tagging service with support for the BeCalm API. The tagger suffered no downtime during the challenge and, as in earlier tests, proved to be highly efficient, consistently processing requests of 5000 abstracts in less than half a minute. In fact, the majority of this time was spent not on the NER task but rather on retrieving the document texts from the challenge servers. The latter was found to be the main bottleneck even when hosting a copy of the tagging service on a Raspberry Pi 3, showing that local document storage or caching would be desirable features to include in future revisions of the API standard.

'One DB to rule them all'-the RING: a Regulatory INteraction Graph combining TFs, genes/proteins, SNPs, diseases and drugs

In the last decade, genomics data have been largely adopted to sketch, study and better understand the complex mechanisms that underlie biological processes. The amount of publicly available data sources has grown accordingly, and several types of regulatory interactions have been collected and documented in literature. Unfortunately, often these efforts do not follow any data naming/interoperability/formatting standards, resulting in high-quality but often uninteroperable heterogeneous data repositories. To efficiently take advantage of the large amount of available data and integrate these heterogeneous sources of information, we built the RING (Regulatory Interaction Graph), an integrative standardized multilevel database of biological interactions able to provide a comprehensive and unmatched high-level perspective on several phenomena that take place in the regulatory cascade and that researchers can use to easily build regulatory networks around entities of interest.

SFPEL-LPI: Sequence-based feature projection ensemble learning for predicting LncRNA-protein interactions

LncRNA-protein interactions play important roles in post-transcriptional gene regulation, poly-adenylation, splicing and translation. Identification of lncRNA-protein interactions helps to understand lncRNA-related activities. Existing computational methods utilize multiple lncRNA features or multiple protein features to predict lncRNA-protein interactions, but features are not available for all lncRNAs or proteins; most of existing methods are not capable of predicting interacting proteins (or lncRNAs) for new lncRNAs (or proteins), which don’t have known interactions. In this paper, we propose the sequence-based feature projection ensemble learning method, “SFPEL-LPI”, to predict lncRNA-protein interactions. First, SFPEL-LPI extracts lncRNA sequence-based features and protein sequence-based features. Second, SFPEL-LPI calculates multiple lncRNA-lncRNA similarities and protein-protein similarities by using lncRNA sequences, protein sequences and known lncRNA-protein interactions. Then, SFPEL-LPI combines multiple similarities and multiple features with a feature projection ensemble learning frame. In computational experiments, SFPEL-LPI accurately predicts lncRNA-protein associations and outperforms other state-of-the-art methods. More importantly, SFPEL-LPI can be applied to new lncRNAs (or proteins). The case studies demonstrate that our method can find out novel lncRNA-protein interactions, which are confirmed by literature. Finally, we construct a user-friendly web server, available at http://www.bioinfotech.cn/SFPEL-LPI/.

LncRNA-protein interactions play important roles in post-transcriptional gene regulation, poly-adenylation, splicing and translation. Identification of lncRNA-protein interactions helps to understand lncRNA-related activities. In this paper, we propose a novel computational method “SFPEL-LPI” to predict lncRNA-protein interactions. SFPEL-LPI makes use of lncRNA sequences, protein sequences and known lncRNA-protein associations to extract features and calculate similarities for lncRNAs and proteins, and then combines them with a feature projection ensemble learning frame. SFPEL-LPI can predict unobserved interactions between lncRNAs and proteins, and also can make predictions for new lncRNAs (or proteins), which have no interactions with any proteins (or lncRNAs). SFPEL-LPI produces high-accuracy performances on the benchmark dataset when evaluated by five-fold cross validation, and outperforms state-of-the-art methods. The case studies demonstrate that SFPEL-LPI can find out novel associations, which are confirmed by literature. To facilitate the lncRNA-protein interaction prediction, we develop a user-friendly web server, available at http://www.bioinfotech.cn/SFPEL-LPI/.

Investigation of RNA-RNA Interactions Using the RISE Database

RNA-RNA interactions (RRIs) are essential to understanding the regulatory mechanisms of RNAs. Mapping RRIs in vivo in a transcriptome-wide manner remained challenging until the recent development of several sequencing-based technologies. However, RRIs generated from large-scale studies had not been systematically collected and analyzed before. This article introduces RISE, a database of the RNA Interactome from Sequencing Experiments. RISE provides a comprehensive collection of RRIs in human, mouse, and yeast, derived from transcriptome-wide sequencing experiments, as well as targeted sequencing studies and other public databases/datasets. To facilitate better understanding of the biological roles of these RRIs, RISE also offers rich functional annotations involving RNAs, and an interactive interface to explore the analysis results. Here, we provide a brief description of the RISE website and a step-by-step protocol for using RISE to study RRIs. © 2018 by John Wiley & Sons, Inc.

Darkness in the Human Gene and Protein Function Space: Widely Modest or Absent Illumination by the Life Science Literature and the Trend for Fewer Protein Function Discoveries Since 2000

The mentioning of gene names in the body of the scientific literature 1901-2017 and their fractional counting is used as a proxy to assess the level of biological function discovery. A literature score of one has been defined as full publication equivalent (FPE), the amount of literature necessary to achieve one publication solely dedicated to a gene. It has been found that less than 5000 human genes have each at least 100 FPEs in the available literature corpus. This group of elite genes (4817 protein-coding genes, 119 non-coding RNAs) attracts the overwhelming majority of the scientific literature about genes. Yet, thousands of proteins have never been mentioned at all, ≈2000 further proteins have not even one FPE of literature and, for ≈4600 additional proteins, the FPE count is below 10. The protein function discovery rate measured as numbers of proteins first mentioned or crossing a threshold of accumulated FPEs in a given year has grown until 2000 but is in decline thereafter. This drop is partially offset by function discoveries for non-coding RNAs. The full human genome sequencing does not boost the function discovery rate. Since 2000, the fastest growing group in the literature is that with at least 500 FPEs per gene.

Genome-wide identification of clusters of predicted microRNA binding sites as microRNA sponge candidates

The number of discovered natural miRNA sponges in plants, viruses, and mammals is increasing steadily. Some sponges like ciRS-7 for miR-7 contain multiple nearby miRNA binding sites. We hypothesize that such clusters of miRNA binding sites on the genome can function together as a sponge. No systematic effort has been made in search for clusters of miRNA targets. Here, we, to our knowledge, make the first genome-wide target site predictions for clusters of mature human miRNAs. For each miRNA, we predict the target sites on a genome-wide scale, build a graph with edge weights based on the pairwise distances between sites, and apply Markov clustering to identify genomic regions with high binding site density. Significant clusters are then extracted based on cluster size difference between real and shuffled genomes preserving local properties such as the GC content. We then use conservation and binding energy to filter a final set of miRNA target site clusters or sponge candidates. Our pipeline predicts 3673 sponge candidates for 1250 miRNAs, including the experimentally verified miR-7 sponge ciRS-7. In addition, we point explicitly to 19 high-confidence candidates overlapping annotated genomic sequence. The full list of candidates is freely available at http://rth.dk/resources/mirnasponge, where detailed properties for individual candidates can be explored, such as alignment details, conservation, accessibility and target profiles, which facilitates selection of sponge candidates for further context specific analysis.

An automated identification and analysis of ontological terms in gastrointestinal diseases and nutrition-related literature provides useful insights

With an unprecedented growth in the biomedical literature, keeping up to date with the new developments presents an immense challenge. Publications are often studied in isolation of the established literature, with interpretation being subjective and often introducing human bias. With ontology-driven annotation of biomedical data gaining popularity in recent years and online databases offering metatags with rich textual information, it is now possible to automatically text-mine ontological terms and complement the laborious task of manual management, interpretation, and analysis of the accumulated literature with downstream statistical analysis. In this paper, we have formulated an automated workflow through which we have identified ontological information, including nutrition-related terms in PubMed abstracts (from 1991 to 2016) for two main types of Inflammatory Bowel Diseases: Crohn’s Disease and Ulcerative Colitis; and two other gastrointestinal (GI) diseases, namely, Coeliac Disease and Irritable Bowel Syndrome. Our analysis reveals unique clustering patterns as well as spatial and temporal trends inherent to the considered GI diseases in terms of literature that has been accumulated so far. Although automated interpretation cannot replace human judgement, the developed workflow shows promising results and can be a useful tool in systematic literature reviews. The workflow is available at https://github.com/KociOrges/pytag.

Expression and network analysis of YBX1 interactors for identification of new drug targets in lung adenocarcinoma

Y-Box Binding protein 1 (YBX-1) is known to be involved in various types of cancers. It's interactors also play major role in various cellular functions. Present work aimed to study the expression profile of the YBX-1 interactors during lung adenocarcinoma (LUAD). The differential expression analysis involved 57 genes from 95 lung adenocarcinoma samples, construction of gene network and topology analysis. A Total of 43 genes were found to be differentially expressed from which 17 genes were found to be down regulated and 26 genes were up-regulated. We observed that Polyadenylate-binding protein 1 (PABPC1), a protein involved in YBX1 translation, is highly correlated with YBX1. The interaction network analysis for a differentially expressed non-coding RNA Growth Arrest Specific 5 (GAS5) suggests that two proteins namely, Growth Arrest Specific 2 (GAS2) and Peripheral myelin protein 22 (PMP22) are potentially involved in LUAD progression. The network analysis and differential expression suggests that Collagen type 1 alpha 2 (COL1A2) can be potential biomarker and target for LUAD.

LncRNA OIP5-AS1 is overexpressed in undifferentiated oral tumors and integrated analysis identifies as a downstream effector of stemness-associated transcription factors

Long non-coding RNAs (lncRNAs) play an important role in the regulation of key cellular processes in early development and cancer. LncRNA Oip5-as1 facilitates stem cell self-renewal in mouse by sponging mmu-miR-7 and modulating NANOG level, yet its role in cancer is less understood. We analyzed OIP5-AS1 expression in oral tumors and in TCGA datasets. We observed overexpression of OIP5-AS1 in oral tumors (P < 0.001) and in tumors of epithelial origin from TCGA. OIP5-AS1 expression was strongly associated with undifferentiated tumors (P = 0.0038). In silico analysis showed miR-7 binding site is conserved in mouse and human OIP5-AS1. However, human NANOG 3'-UTR lost the binding site for hsa-miR-7a-3. Therefore, we screened for other miRNAs that can be sponged by OIP5-AS1 and identified six potential miRNAs and their downstream target genes. Expression analysis showed downregulation of miRNAs and upregulation of downstream target genes, particularly in undifferentiated tumors with high-level of OIP5-AS1 suggesting OIP5-AS1 could post-transcriptionally modulate the downstream target genes. Further, systematic epigenomic analysis of OIP5-AS1 promoter revealed binding motifs for MYC, NANOG and KLF4 suggesting that OIP5-AS1 could be transactivated by stemness-associated transcription factors in cancer. OIP5-AS1 overexpression in undifferentiated oral tumors may be suggestive of enhanced cancer stemness, and consequently, poor clinical outcome.

IFN-λ and microRNAs are important modulators of the pulmonary innate immune response against influenza A (H1N2) infection in pigs

The innate immune system is paramount in the response to and clearance of influenza A virus (IAV) infection in non-immune individuals. Known factors include type I and III interferons and antiviral pathogen recognition receptors, and the cascades of antiviral and pro- and anti-inflammatory gene expression they induce. MicroRNAs (miRNAs) are increasingly recognized to participate in post-transcriptional modulation of these responses, but the temporal dynamics of how these players of the antiviral innate immune response collaborate to combat infection remain poorly characterized. We quantified the expression of miRNAs and protein coding genes in the lungs of pigs 1, 3, and 14 days after challenge with swine IAV (H1N2). Through RT-qPCR we observed a 400-fold relative increase in IFN-λ3 gene expression on day 1 after challenge, and a strong interferon-mediated antiviral response was observed on days 1 and 3 accompanied by up-regulation of genes related to the pro-inflammatory response and apoptosis. Using small RNA sequencing and qPCR validation we found 27 miRNAs that were differentially expressed after challenge, with the highest number of regulated miRNAs observed on day 3. In contrast, the number of protein coding genes found to be regulated due to IAV infection peaked on day 1. Pulmonary miRNAs may thus be aimed at fine-tuning the initial rapid inflammatory response after IAV infection. Specifically, we found five miRNAs (ssc-miR-15a, ssc-miR-18a, ssc-miR-21, ssc-miR-29b, and hsa-miR-590-3p)-four known porcine miRNAs and one novel porcine miRNA candidate-to be potential modulators of viral pathogen recognition and apoptosis. A total of 11 miRNAs remained differentially expressed 14 days after challenge, at which point the infection had cleared. In conclusion, the results suggested a role for miRNAs both during acute infection as well as later, with the potential to influence lung homeostasis and susceptibility to secondary infections in the lungs of pigs after IAV infection.

The yeast noncoding RNA interaction network

This article describes the creation of the first expert manually curated noncoding RNA interaction networks for S. cerevisiae The RNA-RNA and RNA-protein interaction networks have been carefully extracted from the experimental literature and made available through the IntAct database (www.ebi.ac.uk/intact). We provide an initial network analysis and compare their properties to the much larger protein-protein interaction network. We find that the proteins that bind to ncRNAs in the network contain only a small proportion of classical RNA binding domains. We also see an enrichment of WD40 domains suggesting their direct involvement in ncRNA interactions. We discuss the challenges in collecting noncoding RNA interaction data and the opportunities for worldwide collaboration to fill the unmet need for this data.

Assessing the miRNA sponge potential of RUNX1T1 in t(8;21) acute myeloid leukemia

t(8;21) acute myeloid leukemia (AML) is characterized by a translocation between chromosomes 8 and 21 and formation of a distinctive RUNX1-RUNX1T1 fusion transcript. This translocation places RUNX1T1 under control of the RUNX1 promoter leading to a pronounced upregulation of RUNX1T1 transcripts in t(8;21) AML, compared to normal hematopoietic cells. We investigated the role of highly-upregulated RUNX1T1 under the hypothesis that it acts as competing endogenous RNA (ceRNA) titrating microRNAs (miRNAs) away from their target transcripts and thus contributes to AML formation. Using publicly available t(8;21) AML RNA-Seq and miRNA-Seq data available from The Cancer Genome Atlas (TCGA) project, we obtained a network consisting of 605 genes that may act as ceRNAs competing for miRNAs with the suggested RUNX1T1 miRNA sponge. Among the 605 ceRNA candidates, 121 have previously been implied in cancer development. Players in the integrin, cadherin, and Wnt signaling pathways affected by the RUNX1T1 sponge were overrepresented. Finally, among a set of 21 high interest RUNX1T1 ceRNAs we found multiple genes that have previously been linked to AML formation. In conclusion, our study offers a novel look at the role of the RUNX1-RUNX1T1 fusion transcript in t(8;21) AML beyond previously investigated genetic and epigenetic aberrations.