Splicing signature database development to delineate cancer pathways using literature mining and transcriptome machine learning

Alternative splicing (AS) events modulate certain pathways and phenotypic plasticity in cancer. Although previous studies have computationally analyzed splicing events, it is still a challenge to uncover biological functions induced by reliable AS events from tremendous candidates. To provide essential splicing event signatures to assess pathway regulation, we developed a database by collecting two datasets: (i) reported literature and (ii) cancer transcriptome profile. The former includes knowledge-based splicing signatures collected from 63,229 PubMed abstracts using natural language processing, extracted for 202 pathways. The latter is the machine learning-based splicing signatures identified from pan-cancer transcriptome for 16 cancer types and 42 pathways. We established six different learning models to classify pathway activities from splicing profiles as a learning dataset. Top-ranked AS events by learning model feature importance became the signature for each pathway. To validate our learning results, we performed evaluations by (i) performance metrics, (ii) differential AS sets acquired from external datasets, and (iii) our knowledge-based signatures. The area under the receiver operating characteristic values of the learning models did not exhibit any drastic difference. However, random-forest distinctly presented the best performance to compare with the AS sets identified from external datasets and our knowledge-based signatures. Therefore, we used the signatures obtained from the random-forest model. Our database provided the clinical characteristics of the AS signatures, including survival test, molecular subtype, and tumor microenvironment. The regulation by splicing factors was additionally investigated. Our database for developed signatures supported retrieval and visualization system.

Analysis of aberrant miRNA-mRNA interaction networks in prostate cancer to conjecture its molecular mechanisms

BACKGROUND

MicroRNAs (miRNAs) capable of post-transcriptionally regulating mRNA expression are essential to tumor occurrence and progression.

OBJECTIVE

This study aims to find negatively regulatory miRNA-mRNA pairs in prostate adenocarcinoma (PRAD).

METHODS

Combining The Cancer Genome Atlas (TCGA) RNA-Seq data with Gene Expression Omnibus (GEO) mRNA/miRNA expression profiles, differently expressed miRNA/mRNA (DE-miRNAs/DE-mRNAs) were identified. MiRNA-mRNA pairs were screened by miRTarBase and TarBase, databases collecting experimentally confirmed miRNA-mRNA pairs, and verified in 30 paired prostate specimens by real-time reverse transcription polymerase chain reaction (RT-qPCR). The diagnostic values of miRNA-mRNA pairs were measured by receiver operation characteristic (ROC) curve and Decision Curve Analysis (DCA). DAVID-mirPath database and Connectivity Map were employed in GO/KEGG analysis and compounds research. Interactions between miRNA-mRNA pairs and phenotypic features were analyzed with correlation heatmap in hiplot.

RESULTS

Based on TCGA RNA-Seq data, 22 miRNA and 14 mRNA GEO datasets, 67 (20 down and 47 up) miRNAs and 351 (139 up and 212 down) mRNAs were selected. After screening from 2 databases, 8 miRNA (up)-mRNA (down) and 7 miRNA (down)-mRNA (up) pairs were identified with Pearson's correlation in TCGA. By external validation, miR-221-3p (down)/GALNT3 (up) and miR-20a-5p (up)/FRMD6 (down) were chosen. The model combing 4 signatures possessed better diagnostic value. These two miRNA-mRNA pairs were significantly connected with immune cells fraction and tumor immune microenvironment.

CONCLUSIONS

The diagnostic model containing 2 negatively regulatory miRNA-mRNA pairs was established to distinguish PRADs from normal controls.

In Silico and In Vivo Evaluation of microRNA-181c-5p's Role in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is a fatal disease, accounting for 75-85% of primary liver cancers. The conclusive research on miR-181c-5p's role in hepatocarcinogenesis, whether it has oncogenic effects or acts as a tumor repressor, is limited and fluctuating. Therefore, the current study aimed to elucidate the role of miR-181c-5p in HCC in silico and in vivo. The bioinformatics analysis of miR-181c-5p expression data in HCC using several databases strongly shed light on its involvement in HCC development, but also confirmed the fluctuating data around its role. miR-181c-5p was proven here to have an oncogenic role by increasing HepG2 cells' viability as confirmed by MTT analysis. In addition, miR-181c-5p was upregulated in the HCC positive control group and progressed the HCC development and malignant features by its forced expression in an HCC mouse model by targeted delivery using a LA-PAMAM polyplex. This is indicated by the cancerous gross and histological features, and the significant increase in liver function biomarkers. The functional enrichment bioinformatics analyses of miR-181c-5p-downregulated targets in HCC indicated that miR-181c-5p targets were significantly enriched in multiple pathways and biological processes involved in HCC development. Fbxl3, an example for miR-181c-5p potential targets, downregulation and its correlation with miR-181c-5p were validated by qPCR. In conclusion, miR-181c-5p is upregulated in HCC and has an oncogenic role enhancing HCC progression.

SGAEMDA: Predicting miRNA-Disease Associations Based on Stacked Graph Autoencoder

MicroRNA (miRNA)-disease association (MDA) prediction is critical for disease prevention, diagnosis, and treatment. Traditional MDA wet experiments, on the other hand, are inefficient and costly.Therefore, we proposed a multi-layer collaborative unsupervised training base model called SGAEMDA (Stacked Graph Autoencoder-Based Prediction of Potential miRNA-Disease Associations). First, from the original miRNA and disease data, we defined two types of initial features: similarity features and association features. Second, stacked graph autoencoder is then used to learn unsupervised low-dimensional representations of meaningful higher-order similarity features, and we concatenate the association features with the learned low-dimensional representations to obtain the final miRNA-disease pair features. Finally, we used a multilayer perceptron (MLP) to predict scores for unknown miRNA-disease associations. SGAEMDA achieved a mean area under the ROC curve of 0.9585 and 0.9516 in 5-fold and 10-fold cross-validation, which is significantly higher than the other baseline methods. Furthermore, case studies have shown that SGAEMDA can accurately predict candidate miRNAs for brain, breast, colon, and kidney neoplasms.

Cis-regulatory mutations associate with transcriptional and post-transcriptional deregulation of gene regulatory programs in cancers

Most cancer alterations occur in the noncoding portion of the human genome, where regulatory regions control gene expression. The discovery of noncoding mutations altering the cells' regulatory programs has been limited to few examples with high recurrence or high functional impact. Here, we show that transcription factor binding sites (TFBSs) have similar mutation loads to those in protein-coding exons. By combining cancer somatic mutations in TFBSs and expression data for protein-coding and miRNA genes, we evaluate the combined effects of transcriptional and post-transcriptional alterations on the regulatory programs in cancers. The analysis of seven TCGA cohorts culminates with the identification of protein-coding and miRNA genes linked to mutations at TFBSs that are associated with a cascading trans-effect deregulation on the cells' regulatory programs. Our analyses of cis-regulatory mutations associated with miRNAs recurrently predict 12 mature miRNAs (derived from 7 precursors) associated with the deregulation of their target gene networks. The predictions are enriched for cancer-associated protein-coding and miRNA genes and highlight cis-regulatory mutations associated with the dysregulation of key pathways associated with carcinogenesis. By combining transcriptional and post-transcriptional regulation of gene expression, our method predicts cis-regulatory mutations related to the dysregulation of key gene regulatory networks in cancer patients.

Updated review of advances in microRNAs and complex diseases: experimental results, databases, webservers and data fusion

MicroRNAs (miRNAs) are gene regulators involved in the pathogenesis of complex diseases such as cancers, and thus serve as potential diagnostic markers and therapeutic targets. The prerequisite for designing effective miRNA therapies is accurate discovery of miRNA-disease associations (MDAs), which has attracted substantial research interests during the last 15 years, as reflected by more than 55 000 related entries available on PubMed. Abundant experimental data gathered from the wealth of literature could effectively support the development of computational models for predicting novel associations. In 2017, Chen et al. published the first-ever comprehensive review on MDA prediction, presenting various relevant databases, 20 representative computational models, and suggestions for building more powerful ones. In the current review, as the continuation of the previous study, we revisit miRNA biogenesis, detection techniques and functions; summarize recent experimental findings related to common miRNA-associated diseases; introduce recent updates of miRNA-relevant databases and novel database releases since 2017, present mainstream webservers and new webserver releases since 2017 and finally elaborate on how fusion of diverse data sources has contributed to accurate MDA prediction.

Updated review of advances in microRNAs and complex diseases: towards systematic evaluation of computational models

Currently, there exist no generally accepted strategies of evaluating computational models for microRNA-disease associations (MDAs). Though K-fold cross validations and case studies seem to be must-have procedures, the value of K, the evaluation metrics, and the choice of query diseases as well as the inclusion of other procedures (such as parameter sensitivity tests, ablation studies and computational cost reports) are all determined on a case-by-case basis and depending on the researchers' choices. In the current review, we include a comprehensive analysis on how 29 state-of-the-art models for predicting MDAs were evaluated. Based on the analytical results, we recommend a feasible evaluation workflow that would suit any future model to facilitate fair and systematic assessment of predictive performance.

PMMS: Predicting essential miRNAs based on multi-head self-attention mechanism and sequences

Increasing evidence has proved that miRNA plays a significant role in biological progress. In order to understand the etiology and mechanisms of various diseases, it is necessary to identify the essential miRNAs. However, it is time-consuming and expensive to identify essential miRNAs by using traditional biological experiments. It is critical to develop computational methods to predict potential essential miRNAs. In this study, we provided a new computational method (called PMMS) to identify essential miRNAs by using multi-head self-attention and sequences. First, PMMS computes the statistic and structure features and extracts the static feature by concatenating them. Second, PMMS extracts the deep learning original feature (BiLSTM-based feature) by using bi-directional long short-term memory (BiLSTM) and pre-miRNA sequences. In addition, we further obtained the multi-head self-attention feature (MS-based feature) based on BiLSTM-based feature and multi-head self-attention mechanism. By considering the importance of the subsequence of pre-miRNA to the static feature of miRNA, we obtained the deep learning final feature (WA-based feature) based on the weighted attention mechanism. Finally, we concatenated WA-based feature and static feature as an input to the multilayer perceptron) model to predict essential miRNAs. We conducted five-fold cross-validation to evaluate the prediction performance of PMMS. The areas under the ROC curves (AUC), the F1-score, and accuracy (ACC) are used as performance metrics. From the experimental results, PMMS obtained best prediction performances (AUC: 0.9556, F1-score: 0.9030, and ACC: 0.9097). It also outperformed other compared methods. The experimental results also illustrated that PMMS is an effective method to identify essential miRNA.

MicroRNAs in Cancer and Cardiovascular Disease

Although cardiac tumor formation is rare, accumulating evidence suggests that the two leading causes of deaths, cancers, and cardiovascular diseases are similar in terms of pathogenesis, including angiogenesis, immune responses, and fibrosis. These similarities have led to the creation of new exciting field of study called cardio-oncology. Here, we review the similarities between cancer and cardiovascular disease from the perspective of microRNAs (miRNAs). As miRNAs are well-known regulators of translation by binding to the 3'-untranslated regions (UTRs) of messenger RNAs (mRNAs), we carefully dissect how a specific set of miRNAs are both oncomiRs (miRNAs in cancer) and myomiRs (muscle-related miRNAs). Furthermore, from the standpoint of similar pathogenesis, miRNAs categories related to the similar pathogenesis are discussed; namely, angiomiRs, Immune-miRs, and fibromiRs.

CmirC: an integrated database of clustered miRNAs co-localized with copy number variations in cancer

Genomic rearrangements and copy number variations (CNVs) are the major regulators of clustered microRNAs (miRNAs) expression. Several clustered miRNAs are harbored in and around chromosome fragile sites (CFSs) and cancer-associated genomic hotspots. Aberrant expression of such clusters can lead to oncogenic or tumor suppressor activities. Here, we developed CmirC (Clustered miRNAs co-localized with CNVs), a comprehensive database of clustered miRNAs co-localized with CNV regions. The database consists of 481 clustered miRNAs co-localized with CNVs and their expression patterns in 35 cancer types of the TCGA. The portal also provides information on CFSs, miRNA cluster candidates, genomic coordinates, target gene networks, and gene functionality. The web portal is integrated with advanced tools such as JBrowse, NCBI-BLAST, GeneSCF, visNetwork, and NetworkD3 to help the researchers in data analysis, visualization, and browsing. This portal provides a promising avenue for integrated data analytics and offers additional evidence for the complex regulation of clustered miRNAs in cancer. The web portal is freely accessible at http://slsdb.manipal.edu/cmirclust to explore clinically significant miRNAs.

An Exosomal miRNA Biomarker for the Detection of Pancreatic Ductal Adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) remains a difficult tumor to diagnose and treat. To date, PDAC lacks routine screening with no markers available for early detection. Exosomes are 40-150 nm-sized extracellular vesicles that contain DNA, RNA, and proteins. These exosomes are released by all cell types into circulation and thus can be harvested from patient body fluids, thereby facilitating a non-invasive method for PDAC detection. A bioinformatics analysis was conducted utilizing publicly available miRNA pancreatic cancer expression and genome databases. Through this analysis, we identified 18 miRNA with strong potential for PDAC detection. From this analysis, 10 (MIR31, MIR93, MIR133A1, MIR210, MIR330, MIR339, MIR425, MIR429, MIR1208, and MIR3620) were chosen due to high copy number variation as well as their potential to differentiate patients with chronic pancreatitis, neoplasms, and PDAC. These 10 were examined for their mature miRNA expression patterns, giving rise to 18 mature miRs for further analysis. Exosomal RNA from cell culture media was analyzed via RTqPCR and seven mature miRs exhibited statistical significance (miR-31-5p, miR-31-3p, miR-210-3p, miR-339-5p, miR-425-5p, miR-425-3p, and miR-429). These identified biomarkers can potentially be used for early detection of PDAC.

Identification of neoplasm-specific signatures of miRNA interactions by employing a systems biology approach

MicroRNAs represent major regulatory components of the disease epigenome and they constitute powerful biomarkers for the accurate diagnosis and prognosis of various diseases, including cancers. The advent of high-throughput technologies facilitated the generation of a vast amount of miRNA-cancer association data. Computational approaches have been utilized widely to effectively analyze and interpret these data towards the identification of miRNA signatures for diverse types of cancers. Herein, a novel computational workflow was applied to discover core sets of miRNA interactions for the major groups of neoplastic diseases by employing network-based methods. To this end, miRNA-cancer association data from four comprehensive publicly available resources were utilized for constructing miRNA-centered networks for each major group of neoplasms. The corresponding miRNA-miRNA interactions were inferred based on shared functionally related target genes. The topological attributes of the generated networks were investigated in order to detect clusters of highly interconnected miRNAs that form core modules in each network. Those modules that exhibited the highest degree of mutual exclusivity were selected from each graph. In this way, neoplasm-specific miRNA modules were identified that could represent potential signatures for the corresponding diseases.

A Survey on Computational Methods for Investigation on ncRNA-Disease Association through the Mode of Action Perspective

Molecular and sequencing technologies have been successfully used in decoding biological mechanisms of various diseases. As revealed by many novel discoveries, the role of non-coding RNAs (ncRNAs) in understanding disease mechanisms is becoming increasingly important. Since ncRNAs primarily act as regulators of transcription, associating ncRNAs with diseases involves multiple inference steps. Leveraging the fast-accumulating high-throughput screening results, a number of computational models predicting ncRNA-disease associations have been developed. These tools suggest novel disease-related biomarkers or therapeutic targetable ncRNAs, contributing to the realization of precision medicine. In this survey, we first introduce the biological roles of different ncRNAs and summarize the databases containing ncRNA-disease associations. Then, we suggest a new trend in recent computational prediction of ncRNA-disease association, which is the mode of action (MoA) network perspective. This perspective includes integrating ncRNAs with mRNA, pathway and phenotype information. In the next section, we describe computational methodologies widely used in this research domain. Existing computational studies are then summarized in terms of their coverage of the MoA network. Lastly, we discuss the potential applications and future roles of the MoA network in terms of integrating biological mechanisms for ncRNA-disease associations.

A critical review of datasets and computational suites for improving cancer theranostics and biomarker discovery

Cancer has been constantly evolving and so is the research pertaining to cancer diagnosis and therapeutic regimens. Early detection and specific therapeutics are the key features of modern cancer therapy. These requirements can only be fulfilled with the integration of diverse high-throughput technologies. Integration of advanced omics methodology involving genomics, epigenomics, proteomics, and transcriptomics provide a clear understanding of multi-faceted cancer. In the past few years, tremendous high-throughput data have been generated from cancer genomics and epigenomic analyses, which on further methodological analyses can yield better biological insights. The major epigenetic alterations reported in cancer are DNA methylation levels, histone post-translational modifications, and epi-miRNA regulating the oncogenes and tumor suppressor genes. While the genomic analyses like gene expression profiling, cancer gene prediction, and genome annotation divulge the genetic alterations in oncogenes or tumor suppressor genes. Also, systems biology approach using biological networks is being extensively used to identify novel cancer biomarkers. Therefore, integration of these multi-dimensional approaches will help to identify potential diagnostic and therapeutic biomarkers. Here, we reviewed the critical databases and tools dedicated to various epigenomic and genomic alterations in cancer. The review further focuses on the multi-omics resources available for further validating the identified cancer biomarkers. We also highlighted the tools for cancer biomarker discovery using a systems biology approach utilizing genomic and epigenomic data. Biomarkers predicted using such integrative approaches are shown to be more clinically relevant.

RNADisease v4.0: an updated resource of RNA-associated diseases, providing RNA-disease analysis, enrichment and prediction

Numerous studies have shown that RNA plays an important role in the occurrence and development of diseases, and RNA-disease associations are not limited to noncoding RNAs in mammals but also exist for protein-coding RNAs. Furthermore, RNA-associated diseases are found across species including plants and nonmammals. To better analyze diseases at the RNA level and facilitate researchers in exploring the pathogenic mechanism of diseases, we decided to update and change MNDR v3.0 to RNADisease v4.0, a repository for RNA-disease association (http://www.rnadisease.org/ or http://www.rna-society.org/mndr/). Compared to the previous version, new features include: (i) expanded data sources and categories of species, RNA types, and diseases; (ii) the addition of a comprehensive analysis of RNAs from thousands of high-throughput sequencing data of cancer samples and normal samples; (iii) the addition of an RNA-disease enrichment tool and (iv) the addition of four RNA-disease prediction tools. In summary, RNADisease v4.0 provides a comprehensive and concise data resource of RNA-disease associations which contains a total of 3 428 058 RNA-disease entries covering 18 RNA types, 117 species and 4090 diseases to meet the needs of biological research and lay the foundation for future therapeutic applications of diseases.

AntiVIRmiR: A repository of host antiviral miRNAs and their expression along with experimentally validated viral miRNAs and their targets

miRNAs play an essential role in promoting viral infections as well as modulating the antiviral defense. Several miRNA repositories have been developed for different species, e.g., human, mouse, and plant. However, ‘VIRmiRNA’ is the only existing resource for experimentally validated viral miRNAs and their targets. We have developed a ‘AntiVIRmiR’ resource encompassing data on host/virus miRNA expression during viral infection. This resource with 22,741 entries is divided into four sub-databases viz., ‘DEmiRVIR’, ‘AntiVmiR’, ‘VIRmiRNA2’ and ‘VIRmiRTar2’. ‘DEmiRVIR’ has 10,033 differentially expressed host-viral miRNAs for 21 viruses. ‘AntiVmiR’ incorporates 1,642 entries for host miRNAs showing antiviral activity for 34 viruses. Additionally, ‘VIRmiRNA2’ includes 3,340 entries for experimentally validated viral miRNAs from 50 viruses along with 650 viral isomeric sequences for 14 viruses. Further, ‘VIRmiRTar2’ has 7,726 experimentally validated targets for viral miRNAs against 21 viruses. Furthermore, we have also performed network analysis for three sub-databases. Interactions between up/down-regulated human miRNAs and viruses are displayed for ‘AntiVmiR’ as well as ‘DEmiRVIR’. Moreover, ‘VIRmiRTar2’ interactions are shown among different viruses, miRNAs, and their targets. We have provided browse, search, external hyperlinks, data statistics, and useful analysis tools. The database available at https://bioinfo.imtech.res.in/manojk/antivirmir would be beneficial for understanding the host-virus interactions as well as viral pathogenesis.

Inferring human miRNA–disease associations via multiple kernel fusion on GCNII

Increasing evidence shows that the occurrence of human complex diseases is closely related to the mutation and abnormal expression of microRNAs(miRNAs). MiRNAs have complex and fine regulatory mechanisms, which makes it a promising target for drug discovery and disease diagnosis. Therefore, predicting the potential miRNA-disease associations has practical significance. In this paper, we proposed an miRNA–disease association predicting method based on multiple kernel fusion on Graph Convolutional Network via Initial residual and Identity mapping (GCNII), called MKFGCNII. Firstly, we built a heterogeneous network of miRNAs and diseases to extract multi-layer features via GCNII. Secondly, multiple kernel fusion method was applied to weight fusion of embeddings at each layer. Finally, Dual Laplacian Regularized Least Squares was used to predict new miRNA–disease associations by the combined kernel in miRNA and disease spaces. Compared with the other methods, MKFGCNII obtained the highest AUC value of 0.9631. Code is available at https://github.com/cuntjx/bioInfo.

Exploration of shared TF-miRNA‒mRNA and mRNA-RBP-pseudogene networks in type 2 diabetes mellitus and breast cancer

Type 2 diabetes mellitus (T2DM) has been confirmed to be closely associated with breast cancer (BC). However, the shared mechanisms between these diseases remain unclear. By comparing different datasets, we identified shared differentially expressed (DE) RNAs in T2DM and BC, including 427 mRNAs and 6 miRNAs from the GEO(Gene Expression Omnibus) database. We used databases to predict interactions to construct two critical networks. The transcription factor (TF)-miRNA‒mRNA network contained 236 TFs, while the RNA binding protein (RBP)-pseudogene-mRNA network showed that the pseudogene S-phase kinase associated protein 1 pseudogene 1 (SKP1P1) might play a key role in regulating gene expression. The shared mRNAs between T2DM and BC were enriched in cytochrome (CYP) pathways, and further analysis of CPEB1 and COLEC12 expression in cell lines, single cells and other cancers showed that they were strongly correlated with the survival and prognosis of patients with BC. This result suggested that patients with T2DM presenting the downregulation of CPEB1 and COLEC12 might have a higher risk of developing BC. Overall, our work revealed that high expression of CYPs in patients with T2DM might be a susceptibility factor for BC and identified novel gene candidates and immune features that are promising targets for immunotherapy in patients with BC.

Updated review of advances in microRNAs and complex diseases: taxonomy, trends and challenges of computational models

Since the problem proposed in late 2000s, microRNA-disease association (MDA) predictions have been implemented based on the data fusion paradigm. Integrating diverse data sources gains a more comprehensive research perspective, and brings a challenge to algorithm design for generating accurate, concise and consistent representations of the fused data. After more than a decade of research progress, a relatively simple algorithm like the score function or a single computation layer may no longer be sufficient for further improving predictive performance. Advanced model design has become more frequent in recent years, particularly in the form of reasonably combing multiple algorithms, a process known as model fusion. In the current review, we present 29 state-of-the-art models and introduce the taxonomy of computational models for MDA prediction based on model fusion and non-fusion. The new taxonomy exhibits notable changes in the algorithmic architecture of models, compared with that of earlier ones in the 2017 review by Chen et al. Moreover, we discuss the progresses that have been made towards overcoming the obstacles to effective MDA prediction since 2017 and elaborated on how future models can be designed according to a set of new schemas. Lastly, we analysed the strengths and weaknesses of each model category in the proposed taxonomy and proposed future research directions from diverse perspectives for enhancing model performance.

Comparison of plasma- and saliva-derived exosomal miRNA profiles reveals diagnostic potential in head and neck cancer

Background: Head and neck squamous cell carcinomas (HNSCC) lack tumor-specific biomarkers. Exosomes from HNSCC patients carry immunomodulatory molecules, and correlate with clinical parameters. We compared miRNA profiles of plasma- and saliva-derived exosomes to reveal liquid biomarker candidates for HNSCC.

Methods: Exosomes were isolated by differential ultracentrifugation from corresponding plasma and saliva samples from 11 HNSCC patients and five healthy donors (HD). Exosomal miRNA profiles, as determined by nCounter^® SPRINT technology, were analyzed regarding their diagnostic and prognostic potential, correlated to clinical data and integrated into network analysis.

Results: 119 miRNAs overlapped between plasma- and saliva-derived exosomes of HNSCC patients, from which 29 tumor-exclusive miRNAs, associated with TP53, TGFB1, PRDM1, FOX O 1 and CDH1 signaling, were selected. By intra-correlation of tumor-exclusive miRNAs from plasma and saliva, top 10 miRNA candidates with the strongest correlation emerged as diagnostic panels to discriminate cancer and healthy as well as potentially prognostic panels for disease-free survival (DFS). Further, exosomal miRNAs were differentially represented in human papillomavirus (HPV) positive and negative as well as low and high stage disease.

Conclusion: A plasma- and a saliva-derived panel of tumor-exclusive exosomal miRNAs hold great potential as liquid biopsy for discrimination between cancer and healthy as well as HPV status and disease stage. Exosomal miRNAs from both biofluids represent a promising tool for future biomarker studies, emphasizing the possibility to substitute plasma by less-invasive saliva collection.

Prediction of biomarker-disease associations based on graph attention network and text representation

MOTIVATION

The associations between biomarkers and human diseases play a key role in understanding complex pathology and developing targeted therapies. Wet lab experiments for biomarker discovery are costly, laborious and time-consuming. Computational prediction methods can be used to greatly expedite the identification of candidate biomarkers.

RESULTS

Here, we present a novel computational model named GTGenie for predicting the biomarker-disease associations based on graph and text features. In GTGenie, a graph attention network is utilized to characterize diverse similarities of biomarkers and diseases from heterogeneous information resources. Meanwhile, a pretrained BERT-based model is applied to learn the text-based representation of biomarker-disease relation from biomedical literature. The captured graph and text features are then integrated in a bimodal fusion network to model the hybrid entity representation. Finally, inductive matrix completion is adopted to infer the missing entries for reconstructing relation matrix, with which the unknown biomarker-disease associations are predicted. Experimental results on HMDD, HMDAD and LncRNADisease data sets showed that GTGenie can obtain competitive prediction performance with other state-of-the-art methods.

AVAILABILITY

The source code of GTGenie and the test data are available at: https://github.com/Wolverinerine/GTGenie.

High-dimension to high-dimension screening for detecting genome-wide epigenetic and noncoding RNA regulators of gene expression

MOTIVATION

The advancement of high-throughput technology characterizes a wide variety of epigenetic modifications and noncoding RNAs across the genome involved in disease pathogenesis via regulating gene expression. The high dimensionality of both epigenetic/noncoding RNA and gene expression data make it challenging to identify the important regulators of genes. Conducting univariate test for each possible regulator-gene pair is subject to serious multiple comparison burden, and direct application of regularization methods to select regulator-gene pairs is computationally infeasible. Applying fast screening to reduce dimension first before regularization is more efficient and stable than applying regularization methods alone.

RESULTS

We propose a novel screening method based on robust partial correlation to detect epigenetic and noncoding RNA regulators of gene expression over the whole genome, a problem that includes both high-dimensional predictors and high-dimensional responses. Compared to existing screening methods, our method is conceptually innovative that it reduces the dimension of both predictor and response, and screens at both node (regulators or genes) and edge (regulator-gene pairs) levels. We develop data-driven procedures to determine the conditional sets and the optimal screening threshold, and implement a fast iterative algorithm. Simulations and applications to long noncoding RNA and microRNA regulation in Kidney cancer and DNA methylation regulation in Glioblastoma Multiforme illustrate the validity and advantage of our method.

AVAILABILITY AND IMPLEMENTATION

The R package, related source codes and real datasets used in this article are provided at https://github.com/kehongjie/rPCor.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

idenMD-NRF: a ranking framework for miRNA-disease association identification

Identifying miRNA-disease associations is an important task for revealing pathogenic mechanism of complicated diseases. Different computational methods have been proposed. Although these methods obtained encouraging performance for detecting missing associations between known miRNAs and diseases, how to accurately predict associated diseases for new miRNAs is still a difficult task. In this regard, a ranking framework named idenMD-NRF is proposed for miRNA-disease association identification. idenMD-NRF treats the miRNA-disease association identification as an information retrieval task. Given a novel query miRNA, idenMD-NRF employs Learning to Rank algorithm to rank associated diseases based on high-level association features and various predictors. The experimental results on two independent test datasets indicate that idenMD-NRF is superior to other compared predictors. A user-friendly web server of idenMD-NRF predictor is freely available at http://bliulab.net/idenMD-NRF/.

MicroRNA Methylome Signature and Their Functional Roles in Colorectal Cancer Diagnosis, Prognosis, and Chemoresistance

Colorectal cancer (CRC) is one of the leading causes of cancer-related deaths worldwide. Despite significant advances in the diagnostic services and patient care, several gaps remain to be addressed, from early detection, to identifying prognostic variables, effective treatment for the metastatic disease, and the implementation of tailored treatment strategies. MicroRNAs, the short non-coding RNA species, are deregulated in CRC and play a significant role in the occurrence and progression. Nevertheless, microRNA research has historically been based on expression levels to determine its biological significance. The exact mechanism underpinning microRNA deregulation in cancer has yet to be elucidated, but several studies have demonstrated that epigenetic mechanisms play important roles in the regulation of microRNA expression, particularly DNA methylation. However, the methylation profiles of microRNAs remain unknown in CRC patients. Methylation is the next major paradigm shift in cancer detection since large-scale epigenetic alterations are potentially better in identifying and classifying cancers at an earlier stage than somatic mutations. This review aims to provide insight into the current state of understanding of microRNA methylation in CRC. The new knowledge from this study can be utilized for personalized health diagnostics, disease prediction, and monitoring of treatment.

DNRLCNN: A CNN Framework for Identifying MiRNA-Disease Associations Using Latent Feature Matrix Extraction with Positive Samples

Emerging evidence indicates that miRNAs have strong relationships with many human diseases. Investigating the associations will contribute to elucidating the activities of miRNAs and pathogenesis mechanisms, and providing new opportunities for disease diagnosis and drug discovery. Therefore, it is of significance to identify potential associations between miRNAs and diseases. The existing databases about the miRNA-disease associations (MDAs) only provide the known MDAs, which can be regarded as positive samples. However, the unknown MDAs are not sufficient to regard as reliable negative samples. To deal with this uncertainty, we proposed a convolutional neural network (CNN) framework, named DNRLCNN, based on a latent feature matrix extracted by only positive samples to predict MDAs. First, by only considering the positive samples into the calculation process, we captured the latent feature matrix for complex interactions between miRNAs and diseases in low-dimensional space. Then, we constructed a feature vector for each miRNA and disease pair based on the feature representation. Finally, we adopted a modified CNN for the feature vector to predict MDAs. As a result, our model achieves better performance than other state-of-the-art methods which based CNN in fivefold cross-validation on both miRNA-disease association prediction task (average AUC of 0.9030) and miRNA-phenotype association prediction task (average AUC of 0. 9442). In addition, we carried out case studies on two human diseases, and all the top-50 predicted miRNAs for lung neoplasms are confirmed by HMDD v3.2 and dbDEMC 2.0 databases, 98% of the top-50 predicted miRNAs for heart failure are confirmed. The experiment results show that our model has the capability of inferring potential disease-related miRNAs.

dbDEMC 3.0: Functional Exploration of Differentially Expressed miRNAs in Cancers of Human and Model Organisms

MicroRNAs (miRNAs) are important regulators in gene expression. The dysregulation of miRNA expression is widely reported in the transformation from physiological to pathological states of cells. A large number of differentially expressed miRNAs (DEMs) have been identified in various human cancers by using high-throughput technologies, such as microarray and miRNA-seq. Through mining of published studies with high-throughput experiment information, the database of DEMs in human cancers (dbDEMC) was constructed with the aim of providing a systematic resource for the storage and query of the DEMs. Here we report an update of the dbDEMC to version 3.0, which contains two-fold more data entries than the second version and now includes also data from mice and rats. The dbDEMC 3.0 contains 3268 unique DEMs in 40 different cancer types. The current datasets for differential expression analysis have expanded to 9 generalized categories. Moreover, the current release integrates functional annotations of DEMs obtained by using experimentally validated targets. The annotations can be of great benefit to the intensive analysis of the roles of DEMs in cancer. In summary, dbDEMC 3.0 provides a valuable resource for characterizing molecular functions and regulatory mechanisms of DEMs in human cancers. The dbDEMC 3.0 is freely accessible at https://www.biosino.org/dbDEMC.

In silico Methods for Identification of Potential Therapeutic Targets

At the initial stage of drug discovery, identifying novel targets with maximal efficacy and minimal side effects can improve the success rate and portfolio value of drug discovery projects while simultaneously reducing cycle time and cost. However, harnessing the full potential of big data to narrow the range of plausible targets through existing computational methods remains a key issue in this field. This paper reviews two categories of in silico methods—comparative genomics and network-based methods—for finding potential therapeutic targets among cellular functions based on understanding their related biological processes. In addition to describing the principles, databases, software, and applications, we discuss some recent studies and prospects of the methods. While comparative genomics is mostly applied to infectious diseases, network-based methods can be applied to infectious and non-infectious diseases. Nonetheless, the methods often complement each other in their advantages and disadvantages. The information reported here guides toward improving the application of big data-driven computational methods for therapeutic target discovery.

Graphical abstract

Construction of the miRNA-mRNA Regulatory Networks and Explore Their Role in the Development of Lung Squamous Cell Carcinoma

Purpose: MicroRNA (miRNA) binds to target mRNA and inhibit post-transcriptional gene expression. It plays an essential role in regulating gene expression, cell cycle, and biological development. This study aims to identify potential miRNA-mRNA regulatory networks that contribute to the pathogenesis of lung squamous cell carcinoma (LUSC).

Patients and Methods: MiRNA microarray and RNA-Seq datasets were obtained from the gene expression omnibus (GEO) databases, the cancer genome atlas (TCGA), miRcancer, and dbDEMC. The GEO2R tool, “limma” and “DEseq” R packages were used to perform differential expression analysis. Gene enrichment analysis was conducted using the DAVID, DIANA, and Hiplot tools. The miRNA-mRNA regulatory networks were screened from the experimentally validated miRNA-target interactions databases (miRTarBase and TarBase). External validation was carried out in 30 pairs of LUSC tissues by Real-Time Quantitative Reverse Transcription PCR (qRT-PCR). Receiver operating characteristic curve (ROC) and decision curve analysis (DCA) were conducted to evaluate the diagnostic value. Clinical, survival and phenotypic analysis of miRNA-mRNA regulatory networks were further explored.

Results: We screened 5 miRNA and 10 mRNA expression datasets from GEO and identified 7 DE-miRNAs and 270 DE-mRNAs. After databases screening and correlation analysis, four pairs of miRNA-mRNA regulatory networks were screened out. The miRNA-mRNA network of miR-205-5p (up) and PTPRM (down) was validated in 30 pairs of LUSC tissues. MiR-205-5p and PTPRM have good diagnostic efficacy and are expressed differently in different clinical features and are related to tumor immunity.

Conclusion: The research identified a potential miRNA-mRNA regulatory network, providing a new way to explore the genesis and development of LUSC.

iUMRG: multi-layered network-guided propagation modeling for the inference of susceptibility genes and potential drugs against uveal melanoma

Uveal melanoma (UM) is the most common primary malignant intraocular tumor. The use of precision medicine for UM to enable personalized diagnosis, prognosis, and treatment require the development of computer-aided strategies and predictive tools that can identify novel high-confidence susceptibility genes (HSGs) and potential therapeutic drugs. In the present study, a computational framework via propagation modeling on integrated multi-layered molecular networks (abbreviated as iUMRG) was proposed for the systematic inference of HSGs in UM. Under the leave-one-out cross-validation experiments, the iUMRG achieved superior predictive performance and yielded a higher area under the receiver operating characteristic curve value (0.8825) for experimentally verified SGs. In addition, using the experimentally verified SGs as seeds, genome-wide screening was performed to detect candidate HSGs using the iUMRG. Multi-perspective validation analysis indicated that most of the top 50 candidate HSGs were indeed markedly associated with UM carcinogenesis, progression, and outcome. Finally, drug repositioning experiments performed on the HSGs revealed 17 potential targets and 10 potential drugs, of which six have been approved for UM treatment. In conclusion, the proposed iUMRG is an effective supplementary tool in UM precision medicine, which may assist the development of new medical therapies and discover new SGs.

miRCom: Tensor Completion Integrating Multi-View Information to Deduce the Potential Disease-Related miRNA-miRNA Pairs

MicroRNAs (miRNAs) are consistently capable of regulating gene expression synergistically in a combination mode and play a key role in various biological processes associated with the initiation and development of human diseases, which indicate that comprehending the synergistic molecular mechanism of miRNAs may facilitate understanding the pathogenesis of diseases or even overcome it. However, most existing computational methods had an incomprehensive acknowledge of the miRNA synergistic effect on the pathogenesis of complex diseases, or were hard to be extended to a large-scale prediction task of miRNA synergistic combinations for different diseases. In this article, we propose a novel tensor completion framework integrating multi-view miRNAs and diseases information, called miRCom, for the discovery of potential disease-associated miRNA-miRNA pairs. We first construct an incomplete three-order association tensor and several types of similarity matrices based on existing biological knowledge. Then, we formulate an objective function via performing the factorizations of coupled tensor and matrices simultaneously. Finally, we build an optimization schema by adopting the ADMM algorithm. After that, we obtain the prediction of miRNA-miRNA pairs for different diseases from the full tensor. The contrastive experimental results with other approaches verified that miRCom effectively identify the potential disease-related miRNA-miRNA pairs. Moreover, case study results further illustrated that miRNA-miRNA pairs have more biologically significance and prognostic value than single miRNAs.

BioTAGME: A Comprehensive Platform for Biological Knowledge Network Analysis

The inference of novel knowledge and new hypotheses from the current literature analysis is crucial in making new scientific discoveries. In bio-medicine, given the enormous amount of literature and knowledge bases available, the automatic gain of knowledge concerning relationships among biological elements, in the form of semantically related terms (or entities), is rising novel research challenges and corresponding applications. In this regard, we propose BioTAGME, a system that combines an entity-annotation framework based on Wikipedia corpus (i.e., TAGME tool) with a network-based inference methodology (i.e., DT-Hybrid). This integration aims to create an extensive Knowledge Graph modeling relations among biological terms and phrases extracted from titles and abstracts of papers available in PubMed. The framework consists of a back-end and a front-end. The back-end is entirely implemented in Scala and runs on top of a Spark cluster that distributes the computing effort among several machines. The front-end is released through the Laravel framework, connected with the Neo4j graph database to store the knowledge graph.

Investigating the interference of single nucleotide polymorphisms with miRNA mediated gene regulation in pancreatic ductal adenocarcinoma: An in silico approach

BACKGROUND

Pancreatic ductal adenocarcinoma (PDAC) has a strong genetic component and single nucleotide polymorphisms (SNPs) in key genes have been found to modulate the susceptibility of the individuals to the disease. SNPs in 3'-UTR of the target genes or in miRNA seed region has gained much importance as this may lead to impairment of miRNA-mRNA interaction. Not much information about this phenomenon is available with respect to PDAC and we wanted to predict such SNPs which could affect miRNA function in the disease using bioinformatics tools.

METHODS

After identifying the deregulated miRNAs and genes in PDAC, we determined how many of those altered genes are among experimentally validated targets of those miRNAs. Subsequently, SNPs which could alter these miRNA-mRNA interactions were detected using multiple webtools following high stringent conditions. Disease relevance of the SNPs were also evaluated.

RESULTS

We identified a total of 2492 experimentally validated target genes for 303 miRNAs deregulated in PDAC. Our meta-analysis from 363 PDAC patients and 162 control individuals resulted in a set of differentially expressed genes in pancreatic cancer, which was further compared with the miRNA target genes to get targets differentially expressed in pancreatic cancer. We further detected SNPs either in 'seed' region of miRNAs or 'seed-match' sequence of mRNAs either having disruption or creation of miRNA binding site, correlated the expression for each miRNA-SNP-mRNA interaction. Selected SNPs were found to be in LD with important GWAS identified SNPs.

CONCLUSION

Our study, hereby, explores the probable effects of SNPs on miRNA-target mRNA interactions. Through stringent analytical methods, we have identified 3 common variants and 13other rare variants possibly interfering with miRNA mediated gene regulation in PDAC.

miRModuleNet: Detecting miRNA-mRNA Regulatory Modules

Increasing evidence that microRNAs (miRNAs) play a key role in carcinogenesis has revealed the need for elucidating the mechanisms of miRNA regulation and the roles of miRNAs in gene-regulatory networks. A better understanding of the interactions between miRNAs and their mRNA targets will provide a better understanding of the complex biological processes that occur during carcinogenesis. Increased efforts to reveal these interactions have led to the development of a variety of tools to detect and understand these interactions. We have recently described a machine learning approach miRcorrNet, based on grouping and scoring (ranking) groups of genes, where each group is associated with a miRNA and the group members are genes with expression patterns that are correlated with this specific miRNA. The miRcorrNet tool requires two types of -omics data, miRNA and mRNA expression profiles, as an input file. In this study we describe miRModuleNet, which groups mRNA (genes) that are correlated with each miRNA to form a star shape, which we identify as a miRNA-mRNA regulatory module. A scoring procedure is then applied to each module to further assess their contribution in terms of classification. An important output of miRModuleNet is that it provides a hierarchical list of significant miRNA-mRNA regulatory modules. miRModuleNet was further validated on external datasets for their disease associations, and functional enrichment analysis was also performed. The application of miRModuleNet aids the identification of functional relationships between significant biomarkers and reveals essential pathways involved in cancer pathogenesis. The miRModuleNet tool and all other supplementary files are available at https://github.com/malikyousef/miRModuleNet/

PDMDA: predicting deep-level miRNA-disease associations with graph neural networks and sequence features

MOTIVATION

Many studies have shown that microRNAs (miRNAs) play a key role in human diseases. Meanwhile, traditional experimental methods for miRNA-disease association identification are extremely costly, time-consuming and challenging. Therefore, many computational methods have been developed to predict potential associations between miRNAs and diseases. However, those methods mainly predict the existence of miRNA-disease associations, and they cannot predict the deep-level miRNA-disease association types.

RESULTS

In this study, we propose a new end-to-end deep learning method (called PDMDA) to predict deep-level miRNA-disease associations with graph neural networks (GNNs) and miRNA sequence features. Based on the sequence and structural features of miRNAs, PDMDA extracts the miRNA feature representations by a fully connected network (FCN). The disease feature representations are extracted from the disease-gene network and gene-gene interaction network by GNN model. Finally, a multilayer with three fully connected layers and a softmax layer is designed to predict the final miRNA-disease association scores based on the concatenated feature representations of miRNAs and diseases. Note that PDMDA does not take the miRNA-disease association matrix as input to compute the Gaussian interaction profile similarity. We conduct three experiments based on six association type samples (including circulations, epigenetics, target, genetics, known association of which their types are unknown and unknown association samples). We conduct fivefold cross-validation validation to assess the prediction performance of PDMDA. The area under the receiver operating characteristic curve scores is used as metric. The experiment results show that PDMDA can accurately predict the deep-level miRNA-disease associations.

AVAILABILITY AND IMPLEMENTATION

Data and source codes are available at https://github.com/27167199/PDMDA.

MLRDFM: a multi-view Laplacian regularized DeepFM model for predicting miRNA-disease associations

MOTIVATION

MicroRNAs (miRNAs), as critical regulators, are involved in various fundamental and vital biological processes, and their abnormalities are closely related to human diseases. Predicting disease-related miRNAs is beneficial to uncovering new biomarkers for the prevention, detection, prognosis, diagnosis and treatment of complex diseases.

RESULTS

In this study, we propose a multi-view Laplacian regularized deep factorization machine (DeepFM) model, MLRDFM, to predict novel miRNA-disease associations while improving the standard DeepFM. Specifically, MLRDFM improves DeepFM from two aspects: first, MLRDFM takes the relationships among items into consideration by regularizing their embedding features via their similarity-based Laplacians. In this study, miRNA Laplacian regularization integrates four types of miRNA similarity, while disease Laplacian regularization integrates two types of disease similarity. Second, to judiciously train our model, Laplacian eigenmaps are utilized to initialize the weights in the dense embedding layer. The experimental results on the latest HMDD v3.2 dataset show that MLRDFM improves the performance and reduces the overfitting phenomenon of DeepFM. Besides, MLRDFM is greatly superior to the state-of-the-art models in miRNA-disease association prediction in terms of different evaluation metrics with the 5-fold cross-validation. Furthermore, case studies further demonstrate the effectiveness of MLRDFM.

Production and Stabilization of Specific Upregulated Long Noncoding RNA HOXD-AS2 in Glioblastomas Are Mediated by TFE3 and miR-661, Respectively

Differential expression of long noncoding RNAs (lncRNA) plays a key role in the development of gliomas. Because gliomas are the most common primary central nervous system tumor and glioblastomas have poor prognosis, it is urgent to develop new diagnostic methods. We have previously reported that lncRNA HOXD-AS2, which is specifically up-regulated in gliomas, can activate cell cycle and promote the development of gliomas. It is expected to be a new marker for molecular diagnosis of gliomas, but little is known about HOXD-AS2. Here, we demonstrate that TFE3 and miR-661 maintain the high expression level of HOXD-AS2 by regulating its production and degradation. We found that TFE3 acted as a transcription factor binding to the HOXD-AS2 promoter region and raised H3K27ac to activate HOXD-AS2. As the cytoplasmic-located lncRNA, HOXD-AS2 could be degraded by miR-661. This process was inhibited in gliomas due to the low expression of miR-661. Our study explains why HOXD-AS2 was specifically up-regulated in gliomas, helps to understand the molecular characteristics of gliomas, and provids insights for the search for specific markers in gliomas.

Integrated multi-omics analysis reveals miR-20a as a regulator for metabolic colorectal cancer

Single-driver molecular events specific to the metabolic colorectal cancer (CRC) have not been clearly elucidated. Herein, we identified 12 functional miRNAs linked to activated metabolism by integrating multi-omics features in metabolic CRC. These miRNAs exhibited significantly enriched CRC driver miRNAs, significant impacts on CRC cell growth and significantly correlated metabolites. Importantly, miR-20a is minimally expressed in normal colorectal tissues but highly expressed in metabolic CRC, suggesting the potential therapeutic target. Bioinformatics analyses further revealed miR-20a as the most powerful determinant that regulates a cascade of dysregulated events, including Wnt signaling pathway, core enzymes involved in FA metabolism program and triacylglycerol abundances. In vitro assays demonstrated that elevated miR-20a up-regulated FA synthesis enzymes via Wnt/β-catenin signaling, and finally promoted proliferative and migration of metabolic CRC cells. Overall, our study revealed that miR-20a promoted progression of metabolic CRC by regulating FA metabolism and served as a potential target for preventing tumor metastasis.

A miRNA-Disease Association Identification Method Based on Reliable Negative Sample Selection and Improved Single-Hidden Layer Feedforward Neural Network

miRNAs are a category of important endogenous non-coding small RNAs and are ubiquitous in eukaryotes. They are widely involved in the regulatory process of post-transcriptional gene expression and play a critical part in the development of human diseases. By utilizing recent advancements in big data technology, using bioinformatics methods to identify causative miRNA becomes a hot spot. In this paper, a method called RNSSLFN is proposed to identify the miRNA-disease associations by reliable negative sample selection and an improved single-hidden layer feedforward neural network (SLFN). It involves, firstly, obtaining integrated similarity for miRNAs and diseases; next, selecting reliable negative samples from unknown miRNA-disease associations via distinguishing up-regulated or down-regulated miRNAs; then, introducing an improved SLFN to solve the prediction task. The experimental results on the latest data sets HMDD v3.2 and the framework of 5-fold cross-validation (CV) show that the average AUC and AUPR of RNSSLFN achieve 0.9316 and 0.9065 m, respectively, which are superior to the other three state-of-the-art methods. Furthermore, in the case studies of 10 common cancers, more than 70% of the top 30 predicted miRNA-disease association pairs are verified in the databases, which further confirms the reliability and effectiveness of the RNSSLFN model. Generally, RNSSLFN in predicting miRNA-disease associations has prodigious potential and extensive foreground.

Investigation of miRNA expression profiles using cohort samples reveals potential early detectability of colorectal cancers by serum miR-26a-5p before clinical diagnosis

Previous studies have investigated the usefulness of microRNA (miRNA/miR) expression data for the early detection of colorectal cancer (CRC). However, limited data are available regarding miRNAs that detect CRC before clinical diagnoses. Accordingly, the present study investigated the early detectability of CRC by miRNAs using the preserved serum samples of the cohort participants affected with CRC within 2 years of study enrollment. First, the significant miRNAs were revealed using clinical CRC samples for a (seven early CRCs and seven controls) microarray analysis based on significance analysis of microarrays. Next, replicability was verified by reverse transcription-quantitative (RT-q)PCR (eight early CRCs and eight controls, together with 12 CRCs and 12 controls). Finally, early detectability was tested using the cohort samples of Japan Multi-Institutional Collaborative Cohort Study (17 CRCs and 17 controls) to reveal how a certain number of patients developed CRC within 2 years after participation. In the discovery phase, miRNA expression measurements were conducted using a 3D-Gene Human miRNA Oligo Chip for 2,555 miRNAs, and RT-qPCR analyses were performed to validate the replicability. In the first validation set with eight CRCs with early clinical stage and eight age- and gender-matched controls, miR-26a-5p and miR-223-3p demonstrated the highest diagnostic accuracy of area under the curve (AUC)=1.000 (sensitivity and specificity 100%). In an examination of the predictability of CRC incidence using pre-clinical cohort samples, miR-26a-5p demonstrated good predictability of advanced CRC incidence with an AUC of 0.840. Overall, the present study revealed serum miR-26a-5p as a potential early detection marker for CRC.

Biomedical data, computational methods and tools for evaluating disease-disease associations

In recent decades, exploring potential relationships between diseases has been an active research field. With the rapid accumulation of disease-related biomedical data, a lot of computational methods and tools/platforms have been developed to reveal intrinsic relationship between diseases, which can provide useful insights to the study of complex diseases, e.g. understanding molecular mechanisms of diseases and discovering new treatment of diseases. Human complex diseases involve both external phenotypic abnormalities and complex internal molecular mechanisms in organisms. Computational methods with different types of biomedical data from phenotype to genotype can evaluate disease-disease associations at different levels, providing a comprehensive perspective for understanding diseases. In this review, available biomedical data and databases for evaluating disease-disease associations are first summarized. Then, existing computational methods for disease-disease associations are reviewed and classified into five groups in terms of the usages of biomedical data, including disease semantic-based, phenotype-based, function-based, representation learning-based and text mining-based methods. Further, we summarize software tools/platforms for computation and analysis of disease-disease associations. Finally, we give a discussion and summary on the research of disease-disease associations. This review provides a systematic overview for current disease association research, which could promote the development and applications of computational methods and tools/platforms for disease-disease associations.

Screening Differential CircRNAs Expression Profiles Reveals the Regulatory Role of the has_circTPT1_003-has-miR-218-5p-CCNE2/SMC4 Signaling Axis in Bladder Carcinoma Progression

Circular RNAs (circRNAs) are a class of noncoding RNAs closely related to the development and progression of various human cancers. However, it is unclear whether circRNAs play an important role in the development of bladder cancer. We utilized human circRNA array V2 microarrays to screen circRNA expression profiles in bladder cancer tissues. Bioinformatic tools including circBank, dbDEMC 2.0, miRCancer, TarBase v7.0, miRtarbase, TCGA-BLCA, Cytoscape-MCODE, String, ENCORI, and Venny 2.1 were then employed to construct the circRNA-miRNA-mRNA regulatory networks. In total, 105 upregulated circRNAs and 167 downregulated circRNAs (fold change >2 and p < 0.001) were filtered out. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of filtered dysregulated circRNAs disclosed that the circRNAs regulatory network was closely related with mRNA processing and cell cycle, etc. Further excavation analysis showed that seven differentially overexpressed circRNAs including hsa_circ_0000133, hsa_circ_0023610, hsa_circ_0005615, hsa_circ_0030162, hsa_circ_0077007, hsa_circ_0001140, and hsa_circ_0107031 were associated with bladder cancer invasiveness, and the cell cycle signal axis. has_circTPT1_003-has-miR-218-5p-CCNE2/SMC4 was finally clarified as a possible mechanism for bladder cancer progression. Based on results derived from multiple approaches, we identified that has_circTPT1_003-has-miR-218-5p-CCNE2/SMC4 signal axis may be involved in the invasion process of bladder cancer.

Comprehensive bioinformatics analysis of functional molecules in colorectal cancer

BACKGROUND

Colorectal cancer (CRC) is the 3rd most common cancer and the 2nd leading cause of cancer-related death. Numerous studies have found that aberrations in cellular molecules play an important role in the development of tumors. Studying and determining the interactions between these molecules can contribute to the diagnosis, treatment, and prognosis of tumors.

METHODS

The GSE151021, GSE156720, and GSE156719 data sets were analyzed to screen the differentially expressed messenger RNAs (DEmRNAs), long non-coding RNAs (DElncRNAs), and microRNAs (DEmiRNAs) in CRC. Database for Annotation, Visualization and Integrated Discovery (DAVID) and the Search Tool for the Retrieval of Interacting Genes/Proteins software were used to examine gene enrichment and the hub genes. Gene Expression Profiling Interactive Analysis 2 (GEPIA2) and UALCAN was used to verify the expression of the hub genes. To analyze the overall survival (OS) of the hub genes, Kaplan-Meier plotter (KM plotter) was performed. Finally, the miRCancer database, TargetScan, and GSE156719 were used to identify the targets of the identified miRNAs. To predict the lncRNA-miRNA interactions, we used DIANA-LncBase v2 and GSE156720. Finally, the visualization protein‑protein interaction (PPI), competitive endogenous RNA (ceRNA) network was constructed using Cytoscape v3.1.

RESULTS

By analyzing GSE151021 and GSE156720, 23 upregulated mRNAs and 10 downregulated mRNAs were identified as sharing the differentially expressed genes (DEGs) between CRC and adjacent tissues. Furthermore, nucleolar protein 14 (NOP14), the sonic hedgehog (SHH) signaling molecule, phorbol-12-myristate-13-acetate-induced protein 1 (PMAIP1), the BCL2 apoptosis regulator (BCL2), and zinc finger E-box binding homeobox 2 (ZEB2) were considered hub genes. The constructed lncRNA-miRNA-mRNA network revealed 7 intersecting miRNAs (4 upregulated and 3 downregulated), 79 lncRNAs (40 upregulated and 39 downregulated), and 5 mRNAs (3 upregulated and 2 downregulated). Finally, we determined that the dysregulation of lncRNAs, such as HCG16, CASC9, SNHG16, HAND2-AS1, and NR2F1-AS1, secluded altered the expression of several miRNAs, such as hsa-miR-193a-5p, hsa-miR-485-5p, hsa-miR-17-5p, and hsa-miR-92a-3p, and affected the occurrence and development of CRC.

CONCLUSIONS

We identified a series of DElncRNAs, DEmRNAs, and DEmiRNAs in CRC that might be considered potential biomarkers in understanding the complex molecular pathways leading to CRC development.

MIMRDA: A Method Incorporating the miRNA and mRNA Expression Profiles for Predicting miRNA-Disease Associations to Identify Key miRNAs (microRNAs)

Identifying cancer-related miRNAs (or microRNAs) that precisely target mRNAs is important for diagnosis and treatment of cancer. Creating novel methods to identify candidate miRNAs becomes an imminent Frontier of researches in the field. One major obstacle lies in the integration of the state-of-the-art databases. Here, we introduce a novel method, MIMRDA, which incorporates the miRNA and mRNA expression profiles for predicting miRNA-disease associations to identify key miRNAs. As a proof-of-principle study, we use the MIMRDA method to analyze TCGA datasets of 20 types (BLCA, BRCA, CESE, CHOL, COAD, ESCA, HNSC, KICH, KIRC, KIRP, LIHC, LUAD, LUSC, PAAD, PRAD, READ, SKCM, STAD, THCA and UCEC) of cancer, which identified hundreds of top-ranked miRNAs. Some (as Category 1) of them are endorsed by public databases including TCGA, miRTarBase, miR2Disease, HMDD, MISIM, ncDR and mTD; others (as Category 2) are supported by literature evidences. miR-21 (representing Category 1) and miR-1258 (representing Category 2) display the excellent characteristics of biomarkers in multi-dimensional assessments focusing on the function similarity analysis, overall survival analysis, and anti-cancer drugs’ sensitivity or resistance analysis. We compare the performance of the MIMRDA method over the Limma and SPIA packages, and estimate the accuracy of the MIMRDA method in classifying top-ranked miRNAs via the Random Forest simulation test. Our results indicate the superiority and effectiveness of the MIMRDA method, and recommend some top-ranked key miRNAs be potential biomarkers that warrant experimental validations.

MiRNA-disease association prediction based on meta-paths

Since miRNAs can participate in the posttranscriptional regulation of gene expression, they may provide ideas for the development of new drugs or become new biomarkers for drug targets or disease diagnosis. In this work, we propose an miRNA-disease association prediction method based on meta-paths (MDPBMP). First, an miRNA-disease-gene heterogeneous information network was constructed, and seven symmetrical meta-paths were defined according to different semantics. After constructing the initial feature vector for the node, the vector information carried by all nodes on the meta-path instance is extracted and aggregated to update the feature vector of the starting node. Then, the vector information obtained by the nodes on different meta-paths is aggregated. Finally, miRNA and disease embedding feature vectors are used to calculate their associated scores. Compared with the other methods, MDPBMP obtained the highest AUC value of 0.9214. Among the top 50 predicted miRNAs for lung neoplasms, esophageal neoplasms, colon neoplasms and breast neoplasms, 49, 48, 49 and 50 have been verified. Furthermore, for breast neoplasms, we deleted all the known associations between breast neoplasms and miRNAs from the training set. These results also show that for new diseases without known related miRNA information, our model can predict their potential miRNAs. Code and data are available at https://github.com/LiangYu-Xidian/MDPBMP.

Role of BCAR4 in prostate cancer cell autophagy

Background

Increased autophagy of prostate cancer (PC) cells contributes to their resistance to chemotherapy. Recently, we reported that a long non-coding RNA (lncRNA)-breast-cancer anti-estrogen resistance 4 (BCAR4)-is highly expressed in PC and contributes to castration resistance through activation of GLI2 signaling. However, the role of BCAR4 in the regulation of PC cell autophagy is unknown and is the subject of the current study.

Methods

BCAR4 and Beclin-1 levels and the alteration in autophagy pathway genes were assessed in PC using a public database and in our own clinical specimens. The correlation between BCAR4 and Beclin-1 levels in PC and PC cell lines was determined and their regulatory relationship was assessed by overexpression and knockout assay. The final effect on autophagy was measured by microtubule-associated protein 1A/1B-light chain 3 (LC3) levels. The mechanism that underlies the control of Beclin-1 by BCAR4 was analyzed by cancer database and gain-of-function and loss-of-function approaches.

Results

BCAR4 and Beclin-1 were both upregulated in PC and were positively correlated. BCAR4 directly activated Beclin-1 at transcriptional level, which subsequently increased the ratio of LC3 II to LC3I to augment PC cell autophagy. Beclin-1 did not control levels of BCAR4. Mechanically, BCAR4 and Beclin-1 shared several targeting microRNAs, among which miR-15 and miR-146 appeared to be the mediators of the effects of BACR4 on Beclin-1.

Conclusions

BCAR4 may enhance PC cell autophagy through altering miRNA-regulated Beclin-1 expression in PC.