MiRNATIP: a SOM-based miRNA-target interactions predictor

Machine Learning Strategies in MicroRNA Research: Bridging Genome to Phenome

MicroRNAs (miRNAs) have emerged as a prominent layer of regulation of gene expression. This article offers the salient and current aspects of machine learning (ML) tools and approaches from genome to phenome in miRNA research. First, we underline that the complexity in the analysis of miRNA function ranges from their modes of biogenesis to the target diversity in diverse biological conditions. Therefore, it is imperative to first ascertain the miRNA coding potential of genomes and understand the regulatory mechanisms of their expression. This knowledge enables the efficient classification of miRNA precursors and the identification of their mature forms and respective target genes. Second, and because one miRNA can target multiple mRNAs and vice versa, another challenge is the assessment of the miRNA-mRNA target interaction network. Furthermore, long-noncoding RNA (lncRNA)and circular RNAs (circRNAs) also contribute to this complexity. ML has been used to tackle these challenges at the high-dimensional data level. The present expert review covers more than 100 tools adopting various ML approaches pertaining to, for example, (1) miRNA promoter prediction, (2) precursor classification, (3) mature miRNA prediction, (4) miRNA target prediction, (5) miRNA- lncRNA and miRNA-circRNA interactions, (6) miRNA-mRNA expression profiling, (7) miRNA regulatory module detection, (8) miRNA-disease association, and (9) miRNA essentiality prediction. Taken together, we unpack, critically examine, and highlight the cutting-edge synergy of ML approaches and miRNA research so as to develop a dynamic and microlevel understanding of human health and diseases.

Transcriptomic and Bioinformatic Analyses Identifying a Central Mif-Cop9-Nf-kB Signaling Network in Innate Immunity Response of Ciona robusta

The Ascidian C. robusta is a powerful model for studying innate immunity. LPS induction activates inflammatory-like reactions in the pharynx and the expression of several innate immune genes in granulocyte hemocytes such as cytokines, for instance, macrophage migration inhibitory factors (CrMifs). This leads to intracellular signaling involving the Nf-kB signaling cascade that triggers downstream pro-inflammatory gene expression. In mammals, the COP9 (Constitutive photomorphogenesis 9) signalosome (CSN) complex also results in the activation of the NF-kB pathway. It is a highly conserved complex in vertebrates, mainly engaged in proteasome degradation which is essential for maintaining processes such as cell cycle, DNA repair, and differentiation. In the present study, we used bioinformatics and in-silico analyses combined with an in-vivo LPS exposure strategy, next-generation sequencing (NGS), and qRT-PCR to elucidate molecules and the temporal dynamics of Mif cytokines, Csn signaling components, and the Nf-κB signaling pathway in C. robusta. A qRT-PCR analysis of immune genes selected from transcriptome data revealed a biphasic activation of the inflammatory response. A phylogenetic and STRING analysis indicated an evolutionarily conserved functional link between the Mif-Csn-Nf-kB axis in ascidian C. robusta during LPS-mediated inflammation response, finely regulated by non-coding molecules such as microRNAs (miRNAs).

MirDIP 5.2: tissue context annotation and novel microRNA curation

MirDIP is a well-established database that aggregates microRNA-gene human interactions from multiple databases to increase coverage, reduce bias, and improve usability by providing an integrated score proportional to the probability of the interaction occurring. In version 5.2, we removed eight outdated resources, added a new resource (miRNATIP), and ran five prediction algorithms for miRBase and mirGeneDB. In total, mirDIP 5.2 includes 46 364 047 predictions for 27 936 genes and 2734 microRNAs, making it the first database to provide interactions using data from mirGeneDB. Moreover, we curated and integrated 32 497 novel microRNAs from 14 publications to accelerate the use of these novel data. In this release, we also extend the content and functionality of mirDIP by associating contexts with microRNAs, genes, and microRNA-gene interactions. We collected and processed microRNA and gene expression data from 20 resources and acquired information on 330 tissue and disease contexts for 2657 microRNAs, 27 576 genes and 123 651 910 gene-microRNA-tissue interactions. Finally, we improved the usability of mirDIP by enabling the user to search the database using precursor IDs, and we integrated miRAnno, a network-based tool for identifying pathways linked to specific microRNAs. We also provide a mirDIP API to facilitate access to its integrated predictions. Updated mirDIP is available at https://ophid.utoronto.ca/mirDIP.

mintRULS: Prediction of miRNA–mRNA Target Site Interactions Using Regularized Least Square Method

Identification of miRNA–mRNA interactions is critical to understand the new paradigms in gene regulation. Existing methods show suboptimal performance owing to inappropriate feature selection and limited integration of intuitive biological features of both miRNAs and mRNAs. The present regularized least square-based method, mintRULS, employs features of miRNAs and their target sites using pairwise similarity metrics based on free energy, sequence and repeat identities, and target site accessibility to predict miRNA-target site interactions. We hypothesized that miRNAs sharing similar structural and functional features are more likely to target the same mRNA, and conversely, mRNAs with similar features can be targeted by the same miRNA. Our prediction model achieved an impressive AUC of 0.93 and 0.92 in LOOCV and LmiTOCV settings, respectively. In comparison, other popular tools such as miRDB, TargetScan, MBSTAR, RPmirDIP, and STarMir scored AUCs at 0.73, 0.77, 0.55, 0.84, and 0.67, respectively, in LOOCV setting. Similarly, mintRULS outperformed other methods using metrics such as accuracy, sensitivity, specificity, and MCC. Our method also demonstrated high accuracy when validated against experimentally derived data from condition- and cell-specific studies and expression studies of miRNAs and target genes, both in human and mouse.

Transcriptomic Analyses Reveal 2 and 4 Family Members of Cytochromes P450 (CYP) Involved in LPS Inflammatory Response in Pharynx of Ciona robusta

Cytochromes P450 (CYP) are enzymes responsible for the biotransformation of most endogenous and exogenous agents. The expression of each CYP is influenced by a unique combination of mechanisms and factors including genetic polymorphisms, induction by xenobiotics, and regulation by cytokines and hormones. In recent years, Ciona robusta, one of the closest living relatives of vertebrates, has become a model in various fields of biology, in particular for studying inflammatory response. Using an in vivo LPS exposure strategy, next-generation sequencing (NGS) and qRT-PCR combined with bioinformatics and in silico analyses, compared whole pharynx transcripts from naïve and LPS-exposed C. robusta, and we provide the first view of cytochrome genes expression and miRNA regulation in the inflammatory response induced by LPS in a hematopoietic organ. In C. robusta, cytochromes belonging to 2B,2C, 2J, 2U, 4B and 4F subfamilies were deregulated and miRNA network interactions suggest that different conserved and species-specific miRNAs are involved in post-transcriptional regulation of cytochrome genes and that there could be an interplay between specific miRNAs regulating both inflammation and cytochrome molecules in the inflammatory response in C. robusta.

Bioinformatics analysis of mRNA and miRNA microarray to identify the key miRNA-mRNA pairs in cisplatin-resistant ovarian cancer

Background

Ovarian cancer (OC) is a gynecological malignancy with the highest mortality rate. Cisplatin (DDP) based chemotherapy is a standard strategy for ovarian cancer. Despite good response rates for initial chemotherapy, almost 80% of the patients treated with DDP based chemotherapy will experience recurrence due to drug-resistant, which will ultimately result in fatality. The aim of the present study was to examine the pathogenesis and potential molecular markers of cisplatin-resistant OC by studying the differential expression of mRNAs and miRNAs between cisplatin resistant OC cell lines and normal cell lines.

Methods

Two mRNA datasets (GSE58470 and GSE45553) and two miRNA sequence datasets (GSE58469 and GSE148251) were downloaded from the Gene expression omnibus (GEO) database. Differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs) were screened by the NetworkAnalyst. Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted to analyze the biological functions of DEGs. The protein-protein interaction network was constructed using STRING and Cytoscape software to identify the molecular mechanisms of key signaling pathways and cellular activities. FunRich and MiRNATip databases were used to identify the target genes of the DEMs.

Results

A total of 380 DEGs, and 5 DEMs were identified. Protein–protein interaction (PPI) network of DEGs containing 379 nodes and 1049 edges was constructed, and 4 key modules and 24 hub genes related to cisplatin-resistant OC were screened. Two hundred ninety-nine target genes of the 5 DEMs were found out. Subsequently, one of these 299 target genes (UBB) belonging to the hub genes of GSE58470 and GSE45553 was identified by MCODE and CytoHubba,which was regulated by one miRNA (mir-454).

Conclusions

One miRNA–mRNA regulatory pairs (mir-454-UBB) was established. Taken together, our study provided evidence concerning the alteration genes involved in cisplatin-resistant OC, which will help to unravel the mechanisms underlying drug resistant.

ceRNA Network Regulation of TGF-β, WNT, FOXO, Hedgehog Pathways in the Pharynx of Ciona robusta

The transforming growth factor-β (TGF-β) family of cytokines performs a multifunctional signaling, which is integrated and coordinated in a signaling network that involves other pathways, such as Wintless, Forkhead box-O (FOXO) and Hedgehog and regulates pivotal functions related to cell fate in all tissues. In the hematopoietic system, TGF-β signaling controls a wide spectrum of biological processes, from immune system homeostasis to the quiescence and self-renewal of hematopoietic stem cells (HSCs). Recently an important role in post-transcription regulation has been attributed to two type of ncRNAs: microRNAs and pseudogenes. Ciona robusta, due to its philogenetic position close to vertebrates, is an excellent model to investigate mechanisms of post-transcriptional regulation evolutionarily highly conserved in immune homeostasis. The combined use of NGS and bioinformatic analyses suggests that in the pharynx, the hematopoietic organ of Ciona robusta, the Tgf-β, Wnt, Hedgehog and FoxO pathways are involved in tissue homeostasis, as they are in human. Furthermore, ceRNA network interactions and 3'UTR elements analyses of Tgf-β, Wnt, Hedgehog and FoxO pathways genes suggest that different miRNAs conserved (cin-let-7d, cin-mir-92c, cin-mir-153), species-specific (cin-mir-4187, cin-mir-4011a, cin-mir-4056, cin-mir-4150, cin-mir-4189, cin-mir-4053, cin-mir-4016, cin-mir-4075), pseudogenes (ENSCING00000011392, ENSCING00000018651, ENSCING00000007698) and mRNA 3'UTR elements are involved in post-transcriptional regulation in an integrated way in C. robusta.

BITS2019: the sixteenth annual meeting of the Italian society of bioinformatics

The 16th Annual Meeting of the Bioinformatics Italian Society was held in Palermo, Italy, on June 26-28, 2019. More than 80 scientific contributions were presented, including 4 keynote lectures, 31 oral communications and 49 posters. Also, three workshops were organised before and during the meeting. Full papers from some of the works presented in Palermo were submitted for this Supplement of BMC Bioinformatics. Here, we provide an overview of meeting aims and scope. We also shortly introduce selected papers that have been accepted for publication in this Supplement, for a complete presentation of the outcomes of the meeting.

miRNA target identification and prediction as a function of time in gene expression data

The understanding of miRNA target interactions is still limited due to conflicting data and the fact that high-quality validation of targets is a time-consuming process. Faster methods like high-throughput screens and bioinformatics predictions are employed but suffer from several problems. One of these, namely the potential occurrence of downstream (i.e. secondary) effects in high-throughput screens has been only little discussed so far. However, such effects limit usage for both the identification of interactions and for the training of bioinformatics tools. In order to analyse this problem more closely, we performed time-dependent microarray screening experiments overexpressing human miR-517a-3p, and, together with published time-dependent datasets of human miR-17-5p, miR-135b and miR-124 overexpression, we analysed the dynamics of deregulated genes. We show that the number of deregulated targets increases over time, whereas seed sequence content and performance of several miRNA target prediction algorithms actually decrease over time. Bioinformatics recognition success of validated miR-17 targets was comparable to that of data gained only 12 h post-transfection. We therefore argue that the timing of microarray experiments is of critical importance for detecting direct targets with high confidence and for the usability of these data for the training of bioinformatics prediction tools.