Inferring microRNA-disease associations by random walk on a heterogeneous network with multiple data sources

M6APred-EL: A Sequence-Based Predictor for Identifying N6-methyladenosine Sites Using Ensemble Learning

N6-methyladenosine (m⁶A) modification is the most abundant RNA methylation modification and involves various biological processes, such as RNA splicing and degradation. Recent studies have demonstrated the feasibility of identifying m⁶A peaks using high-throughput sequencing techniques. However, such techniques cannot accurately identify specific methylated sites, which is important for a better understanding of m⁶A functions. In this study, we develop a novel machine learning-based predictor called M6APred-EL for the identification of m⁶A sites. To predict m⁶A sites accurately within genomic sequences, we trained an ensemble of three support vector machine classifiers that explore the position-specific information and physical chemical information from position-specific k-mer nucleotide propensity, physical-chemical properties, and ring-function-hydrogen-chemical properties. We examined and compared the performance of our predictor with other state-of-the-art methods of benchmarking datasets. Comparative results showed that the proposed M6APred-EL performed more accurately for m⁶A site identification. Moreover, a user-friendly web server that implements the proposed M6APred-EL is well established and is currently available at http://server.malab.cn/M6APred-EL/. It is expected to be a practical and effective tool for the investigation of m⁶A functional mechanisms.

IRWNRLPI: Integrating Random Walk and Neighborhood Regularized Logistic Matrix Factorization for lncRNA-Protein Interaction Prediction

Long non-coding RNA (lncRNA) plays an important role in many important biological processes and has attracted widespread attention. Although the precise functions and mechanisms for most lncRNAs are still unknown, we are certain that lncRNAs usually perform their functions by interacting with the corresponding RNA- binding proteins. For example, lncRNA-protein interactions play an important role in post transcriptional gene regulation, such as splicing, translation, signaling, and advances in complex diseases. However, experimental verification of lncRNA-protein interactions prediction is time-consuming and laborious. In this work, we propose a computational method, named IRWNRLPI, to find the potential associations between lncRNAs and proteins. IRWNRLPI integrates two algorithms, random walk and neighborhood regularized logistic matrix factorization, which can optimize a lot more than using an algorithm alone. Moreover, the method is semi-supervised and does not require negative samples. Based on the leave-one-out cross validation, we obtain the AUC of 0.9150 and the AUPR of 0.7138, demonstrating its reliable performance. In addition, by means of case study in the “Mus musculus,” many lncRNA-protein interactions which are predicted by our method can be successfully confirmed by experiments. This suggests that IRWNRLPI will be a useful bioinformatics resource in biomedical research.

TLHNMDA: Triple Layer Heterogeneous Network Based Inference for MiRNA-Disease Association Prediction

In recent years, microRNAs (miRNAs) have been confirmed to be involved in many important biological processes and associated with various kinds of human complex diseases. Therefore, predicting potential associations between miRNAs and diseases with the huge number of verified heterogeneous biological datasets will provide a new perspective for disease therapy. In this article, we developed a novel computational model of Triple Layer Heterogeneous Network based inference for MiRNA-Disease Association prediction (TLHNMDA) by using the experimentally verified miRNA-disease associations, miRNA-long noncoding RNA (lncRNA) interactions, miRNA function similarity information, disease semantic similarity information and Gaussian interaction profile kernel similarity for lncRNAs into an triple layer heterogeneous network to predict new miRNA-disease associations. As a result, the AUCs of TLHNMDA are 0.8795 and 0.8795 ± 0.0010 based on leave-one-out cross validation (LOOCV) and 5-fold cross validation, respectively. Furthermore, TLHNMDA was implemented on three complex human diseases to evaluate predictive ability. As a result, 84% (kidney neoplasms), 78% (lymphoma) and 76% (prostate neoplasms) of top 50 predicted miRNAs for the three complex diseases can be verified by biological experiments. In addition, based on the HMDD v1.0 database, 98% of top 50 potential esophageal neoplasms-associated miRNAs were confirmed by experimental reports. It is expected that TLHNMDA could be a useful model to predict potential miRNA-disease associations with high prediction accuracy and stability.

A Computational Systems Biology Approach for Identifying Candidate Drugs for Repositioning for Cardiovascular Disease

We report an in silico method to screen for receptors or pathways that could be targeted to elicit beneficial transcriptional changes in a cellular model of a disease of interest. In our method, we integrate: (1) a dataset of transcriptome responses of a cell line to a panel of drugs; (2) two sets of genes for the disease; and (3) mappings between drugs and the receptors or pathways that they target. We carried out a gene set enrichment analysis (GSEA) test for each of the two gene sets against a list of genes ordered by fold-change in response to a drug in a relevant cell line (HL60), with the overall score for a drug being the difference of the two enrichment scores. Next, we applied GSEA for drug targets based on drugs that have been ranked by their differential enrichment scores. The method ranks drugs by the degree of anti-correlation of their gene-level transcriptional effects on the cell line with the genes in the disease gene sets. We applied the method to data from (1) CMap 2.0; (2) gene sets from two transcriptome profiling studies of atherosclerosis; and (3) a combined dataset of drug/target information. Our analysis recapitulated known targets related to CVD (e.g., PPARγ; HMG-CoA reductase, HDACs) and novel targets (e.g., amine oxidase A, δ-opioid receptor). We conclude that combining disease-associated gene sets, drug-transcriptome-responses datasets and drug-target annotations can potentially be useful as a screening tool for diseases that lack an accepted cellular model for in vitro screening.

Integrated multifactor analysis explores core dysfunctional modules in autism spectrum disorder

Autism spectrum disorder (ASD) is a complex neurodevelopmental disease in early childhood, and growing up to be a major cause of disability in children. However, the underlying molecular mechanism of ASD remains elusive. Hence, we represented integrated multifactor analysis exploring dysfunctional modules based on RNA-Seq data from corpus callosum in 6 patients with ASD and 6 normal individuals. According to protein-protein interactions (PPIs) and WGCNA, we performed co-expression modules analysis for ASD-associated genes, and identified 25 modules with differentially expressed genes (DEGs), observing that genes in these modules were significantly involved in various biological processes in nervous system, sensory system, phylogenetic system and variety of signaling pathways. Then, based on transcriptional and post-transcriptional regulations, integrating transcription factor (TF)-target and RNA-associated interactions, significant regulators of co-expression modules were identified as pivot regulators, including 67 pivot TFs, 13 pivot miRNAs and 6 pivot lncRNAs. GO and KEGG pathway enrichment analysis demonstrated that the pivot miRNAs significantly enriched in neural or mental-associated biological progresses. The pivot TFs were mainly involved in various regulation of transcription, immune system and organs development. Finally, our work deciphered a multifactor dysfunctional co-expression subnetwork involved in ASD, helps uncover core dysfunctional modules for this disease and improves our understanding of its underlying molecular mechanism.

Predicting Potential Drugs for Breast Cancer based on miRNA and Tissue Specificity

Network-based computational method, with the emphasis on biomolecular interactions and biological data integration, has succeeded in drug development and created new directions, such as drug repositioning and drug combination. Drug repositioning, that is finding new uses for existing drugs to treat more patients, offers time, cost and efficiency benefits in drug development, especially when in silico techniques are used. MicroRNAs (miRNAs) play important roles in multiple biological processes and have attracted much scientific attention recently. Moreover, cumulative studies demonstrate that the mature miRNAs as well as their precursors can be targeted by small molecular drugs. At the same time, human diseases result from the disordered interplay of tissue- and cell lineage-specific processes. However, few computational researches predict drug-disease potential relationships based on miRNA data and tissue specificity. Therefore, based on miRNA data and the tissue specificity of diseases, we propose a new method named as miTS to predict the potential treatments for diseases. Firstly, based on miRNAs data, target genes and information of FDA (Food and Drug Administration) approved drugs, we evaluate the relationships between miRNAs and drugs in the tissue-specific PPI (protein-protein) network. Then, we construct a tripartite network: drug-miRNA-disease Finally, we obtain the potential drug-disease associations based on the tripartite network. In this paper, we take breast cancer as case study and focus on the top-30 predicted drugs. 25 of them (83.3%) are found having known connections with breast cancer in CTD (Comparative Toxicogenomics Database) benchmark and the other 5 drugs are potential drugs for breast cancer. We further evaluate the 5 newly predicted drugs from clinical records, literature mining, KEGG pathways enrichment analysis and overlapping genes between enriched pathways. For each of the 5 new drugs, strongly supported evidences can be found in three or more aspects. In particular, Regorafenib (DB08896) has 15 overlapping KEGG pathways with breast cancer and their p-values are all very small. In addition, whether in the literature curation or clinical validation, Regorafenib has a strong correlation with breast cancer. All the facts show that Regorafenib is likely to be a truly effective drug, worthy of our further study. It further follows that our method miTS is effective and practical for predicting new drug indications, which will provide potential values for treatments of complex diseases.

Regularized Multi-View Subspace Clustering for Common Modules Across Cancer Stages

Discovering the common modules that are co-expressed across various stages can lead to an improved understanding of the underlying molecular mechanisms of cancers. There is a shortage of efficient tools for integrative analysis of gene expression and protein interaction networks for discovering common modules associated with cancer progression. To address this issue, we propose a novel regularized multi-view subspace clustering (rMV-spc) algorithm to obtain a representation matrix for each stage and a joint representation matrix that balances the agreement across various stages. To avoid the heterogeneity of data, the protein interaction network is incorporated into the objective of rMV-spc via regularization. Based on the interior point algorithm, we solve the optimization problem to obtain the common modules. By using artificial networks, we demonstrate that the proposed algorithm outperforms state-of-the-art methods in terms of accuracy. Furthermore, the rMV-spc discovers common modules in breast cancer networks based on the breast data, and these modules serve as biomarkers to predict stages of breast cancer. The proposed model and algorithm effectively integrate heterogeneous data for dynamic modules.

Global Similarity Method Based on a Two-tier Random Walk for the Prediction of microRNA-Disease Association

microRNAs (miRNAs) mutation and maladjustment are related to the occurrence and development of human diseases. Studies on disease-associated miRNA have contributed to disease diagnosis and treatment. To address the problems, such as low prediction accuracy and failure to predict the relationship between new miRNAs and diseases and so on, we design a Laplacian score of graphs to calculate the global similarity of networks and propose a Global Similarity method based on a Two-tier Random Walk for the prediction of miRNA-disease association (GSTRW) to reveal the correlation between miRNAs and diseases. This method is a global approach that can simultaneously predict the correlation between all diseases and miRNAs in the absence of negative samples. Experimental results reveal that this method is better than existing approaches in terms of overall prediction accuracy and ability to predict orphan diseases and novel miRNAs. A case study on GSTRW for breast cancer and conlon cancer is also conducted, and the majority of miRNA-disease association can be verified by our experiment. This study indicates that this method is feasible and effective.

Comprehensive analysis of aberrantly expressed lncRNAs and construction of ceRNA network in gastric cancer

Gastric cancer remains fifth most common cancer often diagnosed at an advanced stage and is the second leading cause of cancer-related death worldwide. Long non-coding RNAs (lncRNAs) involved in various cellular pathways are essential for tumor occurrence and progression and they have high potential to promote or suppress the expression of many genes. In this study, we profiled 19 selected cancer-associated lncRNAs in thirty gastric adenocarcinomas and matching normal tissues by qRT-PCR. Our results showed that most of the lncRNAs were significantly upregulated (12/19). Further, we performed bioinformatic screening of miRNAs that share common miRNA response elements (MREs) with lncRNAs and their downstream mRNA targets. The prediction identified three microRNAs (miR-21, miR-145 and miR-148a) and five gastric cancer-specific target genes (EGFR, KLF4, DNMT1 and AGO4) which also showed strong correlation with lncRNAs in regression analysis. Finally, we constructed an integrated lncRNA-miRNA-mRNA interaction network of the candidate genes to understand the post-transcriptional gene regulation. The ceRNA network analysis revealed that the differentially regulated miR-21 and miR-148a were playing as central candidates coordinating sponging activity of the lncRNAs analyzed (H19, TUG1 and MALAT1) in this study and the overexpression of H19 and miR-21 could be a signature event of gastric tumorigenesis that could serve as prognostic indicators and therapeutic targets.

SRMDAP: SimRank and Density-Based Clustering Recommender Model for miRNA-Disease Association Prediction

Aberrant expression of microRNAs (miRNAs) can be applied for the diagnosis, prognosis, and treatment of human diseases. Identifying the relationship between miRNA and human disease is important to further investigate the pathogenesis of human diseases. However, experimental identification of the associations between diseases and miRNAs is time-consuming and expensive. Computational methods are efficient approaches to determine the potential associations between diseases and miRNAs. This paper presents a new computational method based on the SimRank and density-based clustering recommender model for miRNA-disease associations prediction (SRMDAP). The AUC of 0.8838 based on leave-one-out cross-validation and case studies suggested the excellent performance of the SRMDAP in predicting miRNA-disease associations. SRMDAP could also predict diseases without any related miRNAs and miRNAs without any related diseases.

MiroRNA-188 Acts as Tumor Suppressor in Non-Small-Cell Lung Cancer by Targeting MAP3K3

Non-small cell lung cancer (NSCLC) is the most prevalent form of lung cancer. MicroRNAs have been increasingly implicated in NSCLC and may serve as novel therapeutic targets to combat cancer. Here we investigated the functional implication of miR-188 in NSCLC. We first analyzed miR-188 expression in both NSCLC clinical samples and cancer cell lines. Next we investigated its role in A549 and H2126 cells with cell proliferation, migration, and apoptosis assays. To extend the in vitro study, we employed both xenograft model and LSL- K-ras G12D lung cancer model to examine the role of miR-188 in tumorigenesis. Last we tested MAP3K3 as miR-188 target in NSCLC model. MiR-188 expression was significantly downregulated at the NSCLC tumor sites and lung cancer cells. In vitro transfection of miR-188 reduced cell proliferation and migration potential and promoted cell apoptosis. In xenograft model, miR-188 inhibited tumor growth derived from cancer cells. Intranasal miR-188 administration reduced tumor formation in NSCLC animal model. MAP3K3 was validated as direct target of miR-188. Knocking down MAP3K3 in mice also inhibited tumorigenesis in LSL- K-ras G12D model. Our results demonstrate that miR-188 and its downstream target MAP3K3 could be a potential therapeutic target for NSCLC.

Inferring MicroRNA Targets Based on Restricted Boltzmann Machines

Predicting the miRNA-target interactions (MTIs) is a critical task for elucidating mechanistic roles of miRNAs in pathophysiology. However, most existing techniques have a higher false positive because the precise miRNA target mechanisms are poorly known. Considering that ensemble methods can take advantage of the complementary knowledge in different methods, we propose an alternative optimization framework, Inferring MiRNA Targets based on Restricted Boltzmann Machines (IMTRBM), to enhance the accuracy of previous prediction results. First, the proposed method directly constructs a weighted MTI network though the results predicted by individual methods and each miRNA target pair is weighted based on the frequency appearing in these results. Second, we transform the miRNA-target prediction problem into a complete bipartite graph model, named restricted Boltzmann machine, and utilize a practical learning procedure to train our model and make predictions. Our results show that the algorithm outperforms individual miRNA-target prediction approach in the number of validated miRNA targets at cutoffs of top list. Moreover, our framework can tolerate the decrease and increase of predicted MTIs and even discover new miRNA targets, which have been a challenge to predict for any individual methods. Finally, for the miRNAs that are not appearing in IMTRBM, we design a new method to supplement IMTRBM based on the intuition that similar miRNAs have similar functions, which also achieves a comparable result. The source code of IMTRBM is available at https://github.com/liuying201705/IMTRBM.

Prognostic and clinicopathological role of long non-coding RNA UCA1 in various carcinomas

Urothelial cancer associated 1 (UCA1) as an oncogenic long non-coding RNA (LncRNA) was aberrantly upregulated in various solid tumors. Numerous studies have demonstrated overexpression of UCA1 is an unfavorable prognostic indicator in cancer patients. This study aimed to further explore the prognosis role and clinical significance of UCA1 in cancer. Eligible studies were recruited by a systematic search in PubMed, Embase, Cochrane Library and Web of Science databases. A total of 19/16 studies with 1587/1291 cancer patients were included to evaluate the association between UCA1 expression and overall survival (OS) and clinicopathological factors of malignancies by computing hazard ratio (HR), odds ratios (OR) and confidence interval (CI). The meta-analysis indicated overexpression of UCA1 was significantly correlated with unexpected OS in patients with cancer (pooled HR = 1.85, 95% CI 1.62-2.10, p < 0.001). There was also a significantly negative association between high level of UCA1 and poor grade cancer (pooled OR = 2.74, 95% CI 2.04-3.70, p < 0.001) and positive lymphatic metastasis (pooled OR = 2.43, 95% CI 1.72-3.41, p < 0.001). In conclusion, our study suggested that UCA1 was correlated with more advanced clinicopathological features and poor prognosis as a novel predictive biomarker of patients with various tumors.

Which statistical significance test best detects oncomiRNAs in cancer tissues? An exploratory analysis

MicroRNAs(miRNAs) often exert their oncogenic and tumor suppressor functions by suppressing protein-coding genes expressions in cancers and thus have a strong association with cancers' generation, development and metastasis. Through comprehensively understanding differentially expressed miRNAs (oncomiRNA) in tumor tissues, we can elucidate the underlying molecular mechanisms in tumorigenesis and develop novel strategies for cancer diagnosis and treatment. The differential expression of miRNAs can now be analyzed through numerous statistical significance tests based on different principles, which are also available in various R packages. However, the results can be notably different. In this study, we compared miRNAs obtained from 6 common significance tests/R packages (t-test, Limma, DESeq, edgeR, LRT and MARS) with the miRNAs archived in two databases; HMDD 2.0 database, which collects experimentally validated differentially expressed miRNAs, and Infer microRNA-disease association database, which contains the potential disease-associated miRNAs by network forecasting. Finally, we sought the MARS method in DEGseq package more effectively searched out differentially expressed miRNAs than other common methods.

MiR-93-5p Promotes Cell Proliferation through Down-Regulating PPARGC1A in Hepatocellular Carcinoma Cells by Bioinformatics Analysis and Experimental Verification

Peroxisome proliferator-activated receptor gamma coactivator-1 alpha (PPARGC1A, formerly known as PGC-1a) is a transcriptional coactivator and metabolic regulator. Previous studies are mainly focused on the association between PPARGC1A and hepatoma. However, the regulatory mechanism remains unknown. A microRNA associated with cancer (oncomiR), miR-93-5p, has recently been found to play an essential role in tumorigenesis and progression of various carcinomas, including liver cancer. Therefore, this paper aims to explore the regulatory mechanism underlying these two proteins in hepatoma cells. Firstly, an integrative analysis was performed with miRNA–mRNA modules on microarray and The Cancer Genome Atlas (TCGA) data and obtained the core regulatory network and miR-93-5p/PPARGC1A pair. Then, a series of experiments were conducted in hepatoma cells with the results including miR-93-5p upregulated and promoted cell proliferation. Thirdly, the inverse correlation between miR-93-5p and PPARGC1A expression was validated. Finally, we inferred that miR-93-5p plays an essential role in inhibiting PPARGC1A expression by directly targeting the 3′-untranslated region (UTR) of its mRNA. In conclusion, these results suggested that miR-93-5p overexpression contributes to hepatoma development by inhibiting PPARGC1A. It is anticipated to be a promising therapeutic strategy for patients with liver cancer in the future.

Identification and Analysis of Rice Yield-Related Candidate Genes by Walking on the Functional Network

Rice (Oryza sativa L.) is one of the most important staple foods in the world. It is possible to identify candidate genes associated with rice yield using the model of random walk with restart on a functional similarity network. We demonstrated the high performance of this approach by a five-fold cross-validation experiment, as well as the robustness of the parameter r. We also assessed the strength of associations between known seeds and candidate genes in the light of the results scores. The candidates ranking at the top of the results list were considered to be the most relevant rice yield-related genes. This study provides a valuable alternative for rice breeding and biology research. The relevant dataset and script can be downloaded at the website: http://lab.malab.cn/jj/rice.htm.

Prioritizing chronic obstructive pulmonary disease (COPD) candidate genes in COPD-related networks

Chronic obstructive pulmonary disease (COPD) is a multi-factor disease, which could be caused by many factors, including disturbances of metabolism and protein-protein interactions (PPIs). In this paper, a weighted COPD-related metabolic network and a weighted COPD-related PPI network were constructed base on COPD disease genes and functional information. Candidate genes in these weighted COPD-related networks were prioritized by making use of a gene prioritization method, respectively. Literature review and functional enrichment analysis of the top 100 genes in these two networks suggested the correlation of COPD and these genes. The performance of our gene prioritization method was superior to that of ToppGene and ToppNet for genes from the COPD-related metabolic network or the COPD-related PPI network after assessing using leave-one-out cross-validation, literature validation and functional enrichment analysis. The top-ranked genes prioritized from COPD-related metabolic and PPI networks could promote the better understanding about the molecular mechanism of this disease from different perspectives. The top 100 genes in COPD-related metabolic network or COPD-related PPI network might be potential markers for the diagnosis and treatment of COPD.

miRDDCR: a miRNA-based method to comprehensively infer drug-disease causal relationships

Revealing the cause-and-effect mechanism behind drug-disease relationships remains a challenging task. Recent studies suggested that drugs can target microRNAs (miRNAs) and alter their expression levels. In the meanwhile, the inappropriate expression of miRNAs will lead to various diseases. Therefore, targeting specific miRNAs by small-molecule drugs to modulate their activities provides a promising approach to human disease treatment. However, few studies attempt to discover drug-disease causal relationships through the molecular level of miRNAs. Here, we developed a miRNA-based inference method miRDDCR to comprehensively predict drug-disease causal relationships. We first constructed a three-layer drug-miRNA-disease heterogeneous network by combining similarity measurements, existing drug-miRNA associations and miRNA-disease associations. Then, we extended the algorithm of Random Walk to the three-layer heterogeneous network and ranked the potential indications for drugs. Leave-one-out cross-validations and case studies demonstrated that our method miRDDCR can achieve excellent prediction power. Compared with related methods, our causality discovery-based algorithm showed superior prediction ability and highlighted the molecular basis miRNAs, which can be used to assist in the experimental design for drug development and disease treatment. Finally, comprehensively inferred drug-disease causal relationships were released for further studies.

RNA-sequencing reveals genome-wide long non-coding RNAs profiling associated with early development of diabetic nephropathy

Background

Diabetic nephropathy (DN) seriously threatens the lives of patients, and the mechanism of DN remains largely unknown because of the complex regulation between long non-coding RNA (lncRNA) and protein-coding genes. In early development of diabetic nephropathy (DN), pathogenesis remains largely unknown.

Results

We used RNA-sequencing to profile protein-coding and lncRNA gene transcriptome of mouse kidney proximal tubular cells during early stage of DN at various time points. Over 7000 protein-coding and lncRNA genes were differentially expressed, and most of them were time-specific. Nearly 40% of lncRNA genes overlapped with functional element signals using CHIP-Seq data from ENCODE database. Disease progression was characterized by lncRNA expression patterns, rather than protein-coding genes, indicating that the lncRNA genes are potential biomarkers for DN. For gene ontologies related to kidney, enrichment was observed in protein-coding genes co-expressed with neighboring lncRNA genes. Based on protein-coding and lncRNA gene profiles, clustering analysis reveals dynamic expression patterns for kidney, suggesting that they are highly correlated during disease progression. To evaluate translation of mouse model to human conditions, we experimentally validated orthologous genes in human cells in vitro diabetic model. In mouse model, most gene expression patterns were repeated in human cell lines.

Conclusions

These results define dynamic transcriptome and novel functional roles for lncRNAs in diabetic kidney cells; these roles may result in lncRNA-based diagnosis and therapies for DN.

A novel method for multifactorial bio-chemical experiments design based on combinational design theory

Experimental design focuses on describing or explaining the multifactorial interactions that are hypothesized to reflect the variation. The design introduces conditions that may directly affect the variation, where particular conditions are purposely selected for observation. Combinatorial design theory deals with the existence, construction and properties of systems of finite sets whose arrangements satisfy generalized concepts of balance and/or symmetry. In this work, borrowing the concept of "balance" in combinatorial design theory, a novel method for multifactorial bio-chemical experiments design is proposed, where balanced templates in combinational design are used to select the conditions for observation. Balanced experimental data that covers all the influencing factors of experiments can be obtianed for further processing, such as training set for machine learning models. Finally, a software based on the proposed method is developed for designing experiments with covering influencing factors a certain number of times.

A deep ensemble model to predict miRNA-disease association

Cumulative evidence from biological experiments has confirmed that microRNAs (miRNAs) are related to many types of human diseases through different biological processes. It is anticipated that precise miRNA-disease association prediction could not only help infer potential disease-related miRNA but also boost human diagnosis and disease prevention. Considering the limitations of previous computational models, a more effective computational model needs to be implemented to predict miRNA-disease associations. In this work, we first constructed a human miRNA-miRNA similarity network utilizing miRNA-miRNA functional similarity data and heterogeneous miRNA Gaussian interaction profile kernel similarities based on the assumption that similar miRNAs with similar functions tend to be associated with similar diseases, and vice versa. Then, we constructed disease-disease similarity using disease semantic information and heterogeneous disease-related interaction data. We proposed a deep ensemble model called DeepMDA that extracts high-level features from similarity information using stacked autoencoders and then predicts miRNA-disease associations by adopting a 3-layer neural network. In addition to five-fold cross-validation, we also proposed another cross-validation method to evaluate the performance of the model. The results show that the proposed model is superior to previous methods with high robustness.

Identify Huntington's disease associated genes based on restricted Boltzmann machine with RNA-seq data

BACKGROUND

Predicting disease-associated genes is helpful for understanding the molecular mechanisms during the disease progression. Since the pathological mechanisms of neurodegenerative diseases are very complex, traditional statistic-based methods are not suitable for identifying key genes related to the disease development. Recent studies have shown that the computational models with deep structure can learn automatically the features of biological data, which is useful for exploring the characteristics of gene expression during the disease progression.

RESULTS

In this paper, we propose a deep learning approach based on the restricted Boltzmann machine to analyze the RNA-seq data of Huntington's disease, namely stacked restricted Boltzmann machine (SRBM). According to the SRBM, we also design a novel framework to screen the key genes during the Huntington's disease development. In this work, we assume that the effects of regulatory factors can be captured by the hierarchical structure and narrow hidden layers of the SRBM. First, we select disease-associated factors with different time period datasets according to the differentially activated neurons in hidden layers. Then, we select disease-associated genes according to the changes of the gene energy in SRBM at different time periods.

CONCLUSIONS

The experimental results demonstrate that SRBM can detect the important information for differential analysis of time series gene expression datasets. The identification accuracy of the disease-associated genes is improved to some extent using the novel framework. Moreover, the prediction precision of disease-associated genes for top ranking genes using SRBM is effectively improved compared with that of the state of the art methods.

Using computer simulation models to investigate the most promising microRNAs to improve muscle regeneration during ageing

MicroRNAs (miRNAs) regulate gene expression through interactions with target sites within mRNAs, leading to enhanced degradation of the mRNA or inhibition of translation. Skeletal muscle expresses many different miRNAs with important roles in adulthood myogenesis (regeneration) and myofibre hypertrophy and atrophy, processes associated with muscle ageing. However, the large number of miRNAs and their targets mean that a complex network of pathways exists, making it difficult to predict the effect of selected miRNAs on age-related muscle wasting. Computational modelling has the potential to aid this process as it is possible to combine models of individual miRNA:target interactions to form an integrated network. As yet, no models of these interactions in muscle exist. We created the first model of miRNA:target interactions in myogenesis based on experimental evidence of individual miRNAs which were next validated and used to make testable predictions. Our model confirms that miRNAs regulate key interactions during myogenesis and can act by promoting the switch between quiescent/proliferating/differentiating myoblasts and by maintaining the differentiation process. We propose that a threshold level of miR-1 acts in the initial switch to differentiation, with miR-181 keeping the switch on and miR-378 maintaining the differentiation and miR-143 inhibiting myogenesis.

Chronic obstructive pulmonary disease candidate gene prioritization based on metabolic networks and functional information

Chronic obstructive pulmonary disease (COPD) is a multi-factor disease, in which metabolic disturbances played important roles. In this paper, functional information was integrated into a COPD-related metabolic network to assess similarity between genes. Then a gene prioritization method was applied to the COPD-related metabolic network to prioritize COPD candidate genes. The gene prioritization method was superior to ToppGene and ToppNet in both literature validation and functional enrichment analysis. Top-ranked genes prioritized from the metabolic perspective with functional information could promote the better understanding about the molecular mechanism of this disease. Top 100 genes might be potential markers for diagnostic and effective therapies.

Prediction of miRNA-disease Associations using an Evolutionary Tuned Latent Semantic Analysis

MicroRNAs, small non-coding elements implied in gene regulation, are very interesting biomarkers for various diseases such as cancers. They represent potential prodigious biotechnologies for early diagnosis and gene therapies. However, experimental verification of microRNA-disease associations are time-consuming and costly, so that computational modeling is a proper solution. Previously, we designed MiRAI, a predictive method based on distributional semantics, to identify new associations between microRNA molecules and human diseases. Our preliminary results showed very good prediction scores compared to other available methods. However, MiRAI performances depend on numerous parameters that cannot be tuned manually. In this study, a parallel evolutionary algorithm is proposed for finding an optimal configuration of our predictive method. The automatically parametrized version of MiRAI achieved excellent performance. It highlighted new miRNA-disease associations, especially the potential implication of mir-188 and mir-795 in various diseases. In addition, our method allowed to detect several putative false associations contained in the reference database.

A non-negative matrix factorization based method for predicting disease-associated miRNAs in miRNA-disease bilayer network

MOTIVATION

Identification of disease-associated miRNAs (disease miRNAs) is critical for understanding disease etiology and pathogenesis. Since miRNAs exert their functions by regulating the expression of their target mRNAs, several methods based on the target genes were proposed to predict disease miRNA candidates. They achieved only limited success as they all suffered from the high false-positive rate of target prediction results. Alternatively, other prediction methods were based on the observation that miRNAs with similar functions tend to be associated with similar diseases and vice versa. The methods exploited the information about miRNAs and diseases, including the functional similarities between miRNAs, the similarities between diseases, and the associations between miRNAs and diseases. However, how to integrate the multiple kinds of information completely and consider the biological characteristic of disease miRNAs is a challenging problem.

RESULTS

We constructed a bilayer network to represent the complex relationships among miRNAs, among diseases and between miRNAs and diseases. We proposed a non-negative matrix factorization based method to rank, so as to predict, the disease miRNA candidates. The method integrated the miRNA functional similarity, the disease similarity and the miRNA-disease associations seamlessly, which exploited the complex relationships within the bilayer network and the consensus relationship between multiple kinds of information. Considering the correlation between the candidates related to various diseases, it predicted their respective candidates for all the diseases simultaneously. In addition, the sparseness characteristic of disease miRNAs was introduced to generate more reliable prediction model that excludes those noisy candidates. The results on 15 common diseases showed a superior performance of the new method for not only well-characterized diseases but also new ones. A detailed case study on breast neoplasms, colorectal neoplasms, lung neoplasms and 32 other diseases demonstrated the ability of the method for discovering potential disease miRNAs.

AVAILABILITY AND IMPLEMENTATION

The web service for the new method and the list of predicted candidates for all the diseases are available at http://www.bioinfolab.top.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A network-based method using a random walk with restart algorithm and screening tests to identify novel genes associated with Menière's disease

As a chronic illness derived from hair cells of the inner ear, Menière’s disease (MD) negatively influences the quality of life of individuals and leads to a number of symptoms, such as dizziness, temporary hearing loss, and tinnitus. The complete identification of novel genes related to MD would help elucidate its underlying pathological mechanisms and improve its diagnosis and treatment. In this study, a network-based method was developed to identify novel MD-related genes based on known MD-related genes. A human protein-protein interaction (PPI) network was constructed using the PPI information reported in the STRING database. A classic ranking algorithm, the random walk with restart (RWR) algorithm, was employed to search for novel genes using known genes as seed nodes. To make the identified genes more reliable, a series of screening tests, including a permutation test, an interaction test and an enrichment test, were designed to select essential genes from those obtained by the RWR algorithm. As a result, several inferred genes, such as CD4, NOTCH2 and IL6, were discovered. Finally, a detailed biological analysis was performed on fifteen of the important inferred genes, which indicated their strong associations with MD.

Integrated analysis of dosage effect lncRNAs in lung adenocarcinoma based on comprehensive network

Accumulating evidences indicate that cancer-related lncRNAs occur frequent somatic copy number alternation (SCNA). Although individual SCNA lncRNAs have been implicated in tumor biology, their regulatory mechanism has not been assessed in a systematic way. In order to explore the expression characteristics and biological functions of SCNA lncRNAs in cancer, we built a computational framework based on lncRNA expression profiles, lncRNA copy numbers and dosage sensitivity score (DSS). First, we found that the lncRNAs with different DSS were involved in distinct biological processes, while those with the same DSS had similar functions. Second, some of the lncRNAs participated in the progression and metastasis of lung adenocarcinoma (LUAD) through cis-acting regulation. In lncRNA-TF-mRNA network, lncRNAs interacted with 4 TFs and affected the immune system, and further influenced LUAD progression. Third, competing endogenous RNA network analysis inferred that lncRNA ENSG00000240990 competed with HOXA10 to absorb hsa-let-7a/b/f/g-5p and affected patient prognosis in LUAD. Last but not least, by integrating target information of miRNA we also provided a new perspective for the discovery of potential small molecule drugs. In summary, we systematically analyzed the regulatory role of SCNA lncRNAs. This work may facilitate cancer research and serve as the basis for future efforts to understand the role of SCNA lncRNAs, develop novel biomarkers and improve knowledge of tumor biology.

A 3-year observation of testosterone deficiency in Chinese patients with chronic heart failure

Testosterone deficiency is present in a certain proportion men with chronic heart failure (CHF). Low testosterone levels in American and European patients with CHF lead to the high mortality and readmission rates. Interestingly, this relationship has not been studied in Chinese patients. To this end, 167 Chinese men with CHF underwent clinical and laboratory evaluations associated with determinations of testosterone levels. Total testosterone (TT) levels and sex hormone-binding globulin were measured by chemiluminescence or immunoassays assays and free testosterone (FT) levels were calculated, Based upon results from these assays, patients were divided into either a low testosterone (LT; n = 93) or normal testosterone (NT; n = 74) group. Subsequently, records from each patient were reviewed over a follow-up duration of at least 3 years. Patients in the LT group experienced worse cardiac function and a higher prevalence of etiology (ischemic vs. no ischemic) and comorbidity (both P < 0.05). In addition, readmission rates of patients in the LT group were higher than that of patients in the NT group (3.32 ± 1.66 VS 1.57 ± 0.89). Overall, deficiencies in FT levels were accompanied with increased mortalities (HR = 6.301, 95% CI 3.187–12.459, P < .0001).

The landscape of DNA methylation-mediated regulation of long non-coding RNAs in breast cancer

Although systematic studies have identified a host of long non-coding RNAs (lncRNAs) which are involved in breast cancer, the knowledge about the methyla-tion-mediated dysregulation of those lncRNAs remains limited. Here, we integrated multi-omics data to analyze the methylated alteration of lncRNAs in breast invasive carcinoma (BRCA). We found that lncRNAs showed diverse methylation patterns on promoter regions in BRCA. LncRNAs were divided into two categories and four subcategories based on their promoter methylation patterns and expression levels be-tween tumor and normal samples. Through cis-regulatory analysis and gene ontology network, abnormally methylated lncRNAs were identified to be associated with can-cer regulation, proliferation or expression of transcription factors. Competing endog-enous RNA network and functional enrichment analysis of abnormally methylated lncRNAs showed that lncRNAs with different methylation patterns were involved in several hallmarks and KEGG pathways of cancers significantly. Finally, survival analysis based on mRNA modules in networks revealed that lncRNAs silenced by high methylation were associated with prognosis significantly in BRCA. This study enhances the understanding of aberrantly methylated patterns of lncRNAs and pro-vides a novel insight for identifying cancer biomarkers and potential therapeutic tar-gets in breast cancer.

BRWLDA: bi-random walks for predicting lncRNA-disease associations

Increasing efforts have been done to figure out the association between lncRNAs and complex diseases. Many computational models construct various lncRNA similarity networks, disease similarity networks, along with known lncRNA-disease associations to infer novel associations. However, most of them neglect the structural difference between lncRNAs network and diseases network, hierarchical relationships between diseases and pattern of newly discovered associations. In this study, we developed a model that performs Bi-Random Walks to predict novel LncRNA-Disease Associations (BRWLDA in short). This model utilizes multiple heterogeneous data to construct the lncRNA functional similarity network, and Disease Ontology to construct a disease network. It then constructs a directed bi-relational network based on these two networks and available lncRNAs-disease associations. Next, it applies bi-random walks on the network to predict potential associations. BRWLDA achieves reliable and better performance than other comparing methods not only on experiment verified associations, but also on the simulated experiments with masked associations. Case studies further demonstrate the feasibility of BRWLDA in identifying new lncRNA-disease associations.

Improved low-rank matrix recovery method for predicting miRNA-disease association

MicroRNAs (miRNAs) performs crucial roles in various human diseases, but miRNA-related pathogenic mechanisms remain incompletely understood. Revealing the potential relationship between miRNAs and diseases is a critical problem in biomedical research. Considering limitation of existing computational approaches, we develop improved low-rank matrix recovery (ILRMR) for miRNA-disease association prediction. ILRMR is a global method that can simultaneously prioritize potential association for all diseases and does not require negative samples. ILRMR can also identify promising miRNAs for investigating diseases without any known related miRNA. By integrating miRNA-miRNA similarity information, disease-disease similarity information, and miRNA family information to matrix recovery, ILRMR performs better than other methods in cross validation and case studies.

Disease-related gene module detection based on a multi-label propagation clustering algorithm

Detecting disease-related gene modules by analyzing gene expression data is of great significance. It is helpful for exploratory analysis of the interaction mechanisms of genes under complex disease phenotypes. The multi-label propagation algorithm (MLPA) has been widely used in module detection for its fast and easy implementation. The accuracy of MLPA greatly depends on the connections between nodes, and most existing research focuses on measuring the similarity between nodes. However, MLPA does not perform well with loose connections between disease-related genes. Moreover, the biological significance of modules obtained by MLPA has not been demonstrated. To solve these problems, we designed a double label propagation clustering algorithm (DLPCA) based on MLPA to study Huntington's disease. In DLPCA, in addition to category labels, we introduced pathogenic labels to supervise the process of multi-label propagation clustering. The pathogenic labels contain pathogenic information about disease genes and the hierarchical structure of gene expression data. Experimental results demonstrated the superior performance of DLPCA compared with other conventional gene-clustering algorithms.

Identifying novel fruit-related genes in Arabidopsis thaliana based on the random walk with restart algorithm

Fruit is essential for plant reproduction and is responsible for protection and dispersal of seeds. The development and maturation of fruit is tightly regulated by numerous genetic factors that respond to environmental and internal stimulation. In this study, we attempted to identify novel fruit-related genes in a model organism, Arabidopsis thaliana, using a computational method. Based on validated fruit-related genes, the random walk with restart (RWR) algorithm was applied on a protein-protein interaction (PPI) network using these genes as seeds. The identified genes with high probabilities were filtered by the permutation test and linkage tests. In the permutation test, the genes that were selected due to the structure of the PPI network were discarded. In the linkage tests, the importance of each candidate gene was measured from two aspects: (1) its functional associations with validated genes and (2) its similarity with validated genes on gene ontology (GO) terms and KEGG pathways. Finally, 255 inferred genes were obtained, subsequent extensive analysis of important genes revealed that they mainly contribute to ubiquitination (UBQ9, UBQ8, UBQ11, UBQ10), serine hydroxymethyl transfer (SHM7, SHM5, SHM6) or glycol-metabolism (HXKL2_ARATH, CSY5, GAPCP1), suggesting essential roles during the development and maturation of fruit in Arabidopsis thaliana.

A large-scale benchmark of gene prioritization methods

In order to maximize the use of results from high-throughput experimental studies, e.g. GWAS, for identification and diagnostics of new disease-associated genes, it is important to have properly analyzed and benchmarked gene prioritization tools. While prospective benchmarks are underpowered to provide statistically significant results in their attempt to differentiate the performance of gene prioritization tools, a strategy for retrospective benchmarking has been missing, and new tools usually only provide internal validations. The Gene Ontology(GO) contains genes clustered around annotation terms. This intrinsic property of GO can be utilized in construction of robust benchmarks, objective to the problem domain. We demonstrate how this can be achieved for network-based gene prioritization tools, utilizing the FunCoup network. We use cross-validation and a set of appropriate performance measures to compare state-of-the-art gene prioritization algorithms: three based on network diffusion, NetRank and two implementations of Random Walk with Restart, and MaxLink that utilizes network neighborhood. Our benchmark suite provides a systematic and objective way to compare the multitude of available and future gene prioritization tools, enabling researchers to select the best gene prioritization tool for the task at hand, and helping to guide the development of more accurate methods.

Computational prediction of human disease-related microRNAs by path-based random walk

MicroRNAs (miRNAs) are a class of small, endogenous RNAs that are 21–25 nucleotides in length. In animals and plants, miRNAs target specific genes for degradation or translation repression. Discovering disease-related miRNA is fundamental for understanding the pathogenesis of diseases. The association between miRNA and a disease is mainly determined via biological investigation, which is complicated by increased biological information due to big data from different databases. Researchers have utilized different computational methods to harmonize experimental approaches to discover miRNA that articulates restrictively in specific environmental situations. In this work, we present a prediction model that is based on the theory of path features and random walk to obtain a relevancy score of miRNA-related disease. In this model, highly ranked scores are potential miRNA-disease associations. Features were extracted from positive and negative samples of miRNA-disease association. Then, we compared our method with other presented models using the five-fold cross-validation method, which obtained an area under the receiver operating characteristic curve of 88.6%. This indicated that our method has a better performance compared to previous methods and will help future biological investigations.

Identification of exosomes and its signature miRNAs of male and female Cynoglossus semilaevis

Exosomes are small membrane particles which are widely found in various cell lines and physiological fluids in mammalian. MicroRNAs (miRNAs) enclosed in exosomes have been identified as proper signatures for many diseases and response to therapies. However, the composition of exosomes and enclosed miRNAs in fishes has not been investigated. Cynoglossus semilaevis is an important commercial flatfish with ambiguous distinction between males and females before sex maturation, which leads to screening difficulty in reproduction and cultivation. An effective detection method was required for sex differentiation of C. semilaevis. In this work, we successfully identified exosomes in C. semilaevis serum. The analysis of nucleotide composition showed that miRNA dominated in exosomes. Thereafter the miRNA profiles in exosomes from males and females were sequenced and compared to identify the signature miRNAs corresponding to sex differentiation. The functions of signature miRNAs were analyzed by target matching and annotation. Furthermore, 7 miRNAs with high expression in males were selected from signature miRNAs as the markers for sex identification with their expression profiles verified by real time quantitative PCR. Exosomes were first found in fish serum in this work. Investigation of marker miRNAs supplies an effective index for the filtration of male and female C. semilaevis in cultivation.

Cross disease analysis of co-functional microRNA pairs on a reconstructed network of disease-gene-microRNA tripartite

BACKGROUND

MicroRNAs always function cooperatively in their regulation of gene expression. Dysfunctions of these co-functional microRNAs can play significant roles in disease development. We are interested in those multi-disease associated co-functional microRNAs that regulate their common dysfunctional target genes cooperatively in the development of multiple diseases. The research is potentially useful for human disease studies at the transcriptional level and for the study of multi-purpose microRNA therapeutics.

METHODS AND RESULTS

We designed a computational method to detect multi-disease associated co-functional microRNA pairs and conducted cross disease analysis on a reconstructed disease-gene-microRNA (DGR) tripartite network. The construction of the DGR tripartite network is by the integration of newly predicted disease-microRNA associations with those relationships of diseases, microRNAs and genes maintained by existing databases. The prediction method uses a set of reliable negative samples of disease-microRNA association and a pre-computed kernel matrix instead of kernel functions. From this reconstructed DGR tripartite network, multi-disease associated co-functional microRNA pairs are detected together with their common dysfunctional target genes and ranked by a novel scoring method. We also conducted proof-of-concept case studies on cancer-related co-functional microRNA pairs as well as on non-cancer disease-related microRNA pairs.

CONCLUSIONS

With the prioritization of the co-functional microRNAs that relate to a series of diseases, we found that the co-function phenomenon is not unusual. We also confirmed that the regulation of the microRNAs for the development of cancers is more complex and have more unique properties than those of non-cancer diseases.

Large-scale prediction of microRNA-disease associations by combinatorial prioritization algorithm

Identification of the associations between microRNA molecules and human diseases from large-scale heterogeneous biological data is an important step for understanding the pathogenesis of diseases in microRNA level. However, experimental verification of microRNA-disease associations is expensive and time-consuming. To overcome the drawbacks of conventional experimental methods, we presented a combinatorial prioritization algorithm to predict the microRNA-disease associations. Importantly, our method can be used to predict microRNAs (diseases) associated with the diseases (microRNAs) without the known associated microRNAs (diseases). The predictive performance of our proposed approach was evaluated and verified by the internal cross-validations and external independent validations based on standard association datasets. The results demonstrate that our proposed method achieves the impressive performance for predicting the microRNA-disease association with the Area Under receiver operation characteristic Curve (AUC), 86.93%, which is indeed outperform the previous prediction methods. Particularly, we observed that the ensemble-based method by integrating the predictions of multiple algorithms can give more reliable and robust prediction than the single algorithm, with the AUC score improved to 92.26%. We applied our combinatorial prioritization algorithm to lung neoplasms and breast neoplasms, and revealed their top 30 microRNA candidates, which are in consistent with the published literatures and databases.

Extracting microRNA-gene relations from biomedical literature using distant supervision

Many biomedical relation extraction approaches are based on supervised machine learning, requiring an annotated corpus. Distant supervision aims at training a classifier by combining a knowledge base with a corpus, reducing the amount of manual effort necessary. This is particularly useful for biomedicine because many databases and ontologies have been made available for many biological processes, while the availability of annotated corpora is still limited. We studied the extraction of microRNA-gene relations from text. MicroRNA regulation is an important biological process due to its close association with human diseases. The proposed method, IBRel, is based on distantly supervised multi-instance learning. We evaluated IBRel on three datasets, and the results were compared with a co-occurrence approach as well as a supervised machine learning algorithm. While supervised learning outperformed on two of those datasets, IBRel obtained an F-score 28.3 percentage points higher on the dataset for which there was no training set developed specifically. To demonstrate the applicability of IBRel, we used it to extract 27 miRNA-gene relations from recently published papers about cystic fibrosis. Our results demonstrate that our method can be successfully used to extract relations from literature about a biological process without an annotated corpus. The source code and data used in this study are available at https://github.com/AndreLamurias/IBRel.

Integrative analyses of transcriptome sequencing identify novel functional lncRNAs in esophageal squamous cell carcinoma

Long non-coding RNAs (lncRNAs) have a critical role in cancer initiation and progression, and thus may mediate oncogenic or tumor suppressing effects, as well as be a new class of cancer therapeutic targets. We performed high-throughput sequencing of RNA (RNA-seq) to investigate the expression level of lncRNAs and protein-coding genes in 30 esophageal samples, comprised of 15 esophageal squamous cell carcinoma (ESCC) samples and their 15 paired non-tumor tissues. We further developed an integrative bioinformatics method, denoted URW-LPE, to identify key functional lncRNAs that regulate expression of downstream protein-coding genes in ESCC. A number of known onco-lncRNA and many putative novel ones were effectively identified by URW-LPE. Importantly, we identified lncRNA625 as a novel regulator of ESCC cell proliferation, invasion and migration. ESCC patients with high lncRNA625 expression had significantly shorter survival time than those with low expression. LncRNA625 also showed specific prognostic value for patients with metastatic ESCC. Finally, we identified E1A-binding protein p300 (EP300) as a downstream executor of lncRNA625-induced transcriptional responses. These findings establish a catalog of novel cancer-associated functional lncRNAs, which will promote our understanding of lncRNA-mediated regulation in this malignancy.

Network-based collaborative filtering recommendation model for inferring novel disease-related miRNAs

According to the miRNA and disease similarity network, the unknown associations are predicted by combining the known miRNA-disease association network based on collaborative filtering recommendation algorithm.

A Reconfiguration Strategy of Distribution Networks Considering Node Importance

Node importance degree is a vital index in distribution network reconfiguration because it reflects the robustness of the network structure by evaluating node importance. Since the traditional reconfiguration ignores this index, the reconstructed network structure may be vulnerable which would reduce the security and stability of the distribution systems. This paper presents a novel reconfiguration strategy considering the node importance. The optimization objectives are the improvement of the node importance degree and the reduction of power loss. To balance the objectives, the reconfiguration mathematical model is formulated as a compound objective function with weight coefficients. Then the quantum particle swarm algorithm is employed to address this compound objective optimization problem. The strategy can model different scenarios network reconfiguration by adjusting the weight vector based on the tendencies of the utility decision maker. Illustrative examples verify the effectiveness of the proposed strategy.

The Human Kinome Targeted by FDA Approved Multi-Target Drugs and Combination Products: A Comparative Study from the Drug-Target Interaction Network Perspective

The human kinome is one of the most productive classes of drug target, and there is emerging necessity for treating complex diseases by means of polypharmacology (multi-target drugs and combination products). However, the advantages of the multi-target drugs and the combination products are still under debate. A comparative analysis between FDA approved multi-target drugs and combination products, targeting the human kinome, was conducted by mapping targets onto the phylogenetic tree of the human kinome. The approach of network medicine illustrating the drug-target interactions was applied to identify popular targets of multi-target drugs and combination products. As identified, the multi-target drugs tended to inhibit target pairs in the human kinome, especially the receptor tyrosine kinase family, while the combination products were able to against targets of distant homology relationship. This finding asked for choosing the combination products as a better solution for designing drugs aiming at targets of distant homology relationship. Moreover, sub-networks of drug-target interactions in specific disease were generated, and mechanisms shared by multi-target drugs and combination products were identified. In conclusion, this study performed an analysis between approved multi-target drugs and combination products against the human kinome, which could assist the discovery of next generation polypharmacology.

A Survey of Methods for Constructing Rooted Phylogenetic Networks

Rooted phylogenetic networks are primarily used to represent conflicting evolutionary information and describe the reticulate evolutionary events in phylogeny. So far a lot of methods have been presented for constructing rooted phylogenetic networks, of which the methods based on the decomposition property of networks and by means of the incompatible graph (such as the CASS, the LNETWORK and the BIMLR) are more efficient than other available methods. The paper will discuss and compare these methods by both the practical and artificial datasets, in the aspect of the running time of the methods and the effective of constructed phylogenetic networks. The results show that the LNETWORK can construct much simper networks than the others.

Network Consistency Projection for Human miRNA-Disease Associations Inference

Prediction and confirmation of the presence of disease-related miRNAs is beneficial to understand disease mechanisms at the miRNA level. However, the use of experimental verification to identify disease-related miRNAs is expensive and time-consuming. Effective computational approaches used to predict miRNA-disease associations are highly specific. In this study, we develop the Network Consistency Projection for miRNA-Disease Associations (NCPMDA) method to reveal the potential associations between miRNAs and diseases. NCPMDA is a non-parametric universal network-based method that can simultaneously predict miRNA-disease associations in all diseases but does not require negative samples. NCPMDA can also confirm the presence of miRNAs in isolated diseases (diseases without any known miRNA association). Leave-one-out cross validation and case studies have shown that the predictive performance of NCPMDA is superior over that of previous method.

lncScore: alignment-free identification of long noncoding RNA from assembled novel transcripts

RNA-Seq based transcriptome assembly has been widely used to identify novel lncRNAs. However, the best-performing transcript reconstruction methods merely identified 21% of full-length protein-coding transcripts from H. sapiens. Those partial-length protein-coding transcripts are more likely to be classified as lncRNAs due to their incomplete CDS, leading to higher false positive rate for lncRNA identification. Furthermore, potential sequencing or assembly error that gain or abolish stop codons also complicates ORF-based prediction of lncRNAs. Therefore, it remains a challenge to identify lncRNAs from the assembled transcripts, particularly the partial-length ones. Here, we present a novel alignment-free tool, lncScore, which uses a logistic regression model with 11 carefully selected features. Compared to other state-of-the-art alignment-free tools (e.g. CPAT, CNCI, and PLEK), lncScore outperforms them on accurately distinguishing lncRNAs from mRNAs, especially partial-length mRNAs in the human and mouse datasets. In addition, lncScore also performed well on transcripts from five other species (Zebrafish, Fly, C. elegans, Rat, and Sheep). To speed up the prediction, multithreading is implemented within lncScore, and it only took 2 minute to classify 64,756 transcripts and 54 seconds to train a new model with 21,000 transcripts with 12 threads, which is much faster than other tools. lncScore is available at https://github.com/WGLab/lncScore.

OAHG: an integrated resource for annotating human genes with multi-level ontologies

OAHG, an integrated resource, aims to establish a comprehensive functional annotation resource for human protein-coding genes (PCGs), miRNAs, and lncRNAs by multi-level ontologies involving Gene Ontology (GO), Disease Ontology (DO), and Human Phenotype Ontology (HPO). Many previous studies have focused on inferring putative properties and biological functions of PCGs and non-coding RNA genes from different perspectives. During the past several decades, a few of databases have been designed to annotate the functions of PCGs, miRNAs, and lncRNAs, respectively. A part of functional descriptions in these databases were mapped to standardize terminologies, such as GO, which could be helpful to do further analysis. Despite these developments, there is no comprehensive resource recording the function of these three important types of genes. The current version of OAHG, release 1.0 (Jun 2016), integrates three ontologies involving GO, DO, and HPO, six gene functional databases and two interaction databases. Currently, OAHG contains 1,434,694 entries involving 16,929 PCGs, 637 miRNAs, 193 lncRNAs, and 24,894 terms of ontologies. During the performance evaluation, OAHG shows the consistencies with existing gene interactions and the structure of ontology. For example, terms with more similar structure could be associated with more associated genes (Pearson correlation γ² = 0.2428, p < 2.2e-16).

Differential network analysis from cross-platform gene expression data

Understanding how the structure of gene dependency network changes between two patient-specific groups is an important task for genomic research. Although many computational approaches have been proposed to undertake this task, most of them estimate correlation networks from group-specific gene expression data independently without considering the common structure shared between different groups. In addition, with the development of high-throughput technologies, we can collect gene expression profiles of same patients from multiple platforms. Therefore, inferring differential networks by considering cross-platform gene expression profiles will improve the reliability of network inference. We introduce a two dimensional joint graphical lasso (TDJGL) model to simultaneously estimate group-specific gene dependency networks from gene expression profiles collected from different platforms and infer differential networks. TDJGL can borrow strength across different patient groups and data platforms to improve the accuracy of estimated networks. Simulation studies demonstrate that TDJGL provides more accurate estimates of gene networks and differential networks than previous competing approaches. We apply TDJGL to the PI3K/AKT/mTOR pathway in ovarian tumors to build differential networks associated with platinum resistance. The hub genes of our inferred differential networks are significantly enriched with known platinum resistance-related genes and include potential platinum resistance-related genes.