Walking the interactome for prioritization of candidate disease genes

Identification of Carcinogenic Chemicals with Network Embedding and Deep Learning Methods

Background:

Cancer is the second leading cause of human death in the world. To date, many factors have been confirmed to be the cause of cancer. Among them, carcinogenic chemicals have been widely accepted as the important ones. Traditional methods for detecting carcinogenic chemicals are of low efficiency and high cost.

Objective:

The aim of this study was to design an efficient computational method for the identification of carcinogenic chemicals.

Methods:

A new computational model was proposed for detecting carcinogenic chemicals. As a data-driven model, carcinogenic and non-carcinogenic chemicals were obtained from Carcinogenic Potency Database (CPDB). These chemicals were represented by features extracted from five chemical networks, representing five types of chemical associations, via a network embedding method, Mashup. Obtained features were fed into a powerful deep learning method, recurrent neural network, to build the model.

Results:

The jackknife test on such model provided the F-measure of 0.971 and AUROC of 0.971.

Conclusion:

The proposed model was quite effective and was superior to the models with traditional machine learning algorithms, classic chemical encoding schemes or direct usage of chemical associations.

PRSS1 Upregulation Predicts Platinum Resistance in Ovarian Cancer Patients

Ovarian cancer is the most frequent cause of death among gynecologic malignancies. A total of 80% of patients who have completed platinum-based chemotherapy suffer from relapse and develop resistance within 2 years. In the present study, we obtained patients' complete platinum (cisplatin and carboplatin) medication information from The Cancer Genome Atlas database and then divided them into two categories: resistance and sensitivity. Difference analysis was performed to screen differentially expressed genes (DEgenes) related to platinum response. Subsequently, we annotated DEgenes into the protein–protein interaction network as seed nodes and analyzed them by random walk. Finally, second-ranking protease serine 1 gene (PRSS1) was selected as a candidate gene for verification analysis. PRSS1's expression pattern was continuously studied in Oncomine and cBio Cancer Genomic Portal databases, revealing the key roles of PRSS1 in ovarian cancer formation. Hereafter, we conducted in-depth explorations on PRSS1's platinum response to ovarian cancer through tissue and cytological experiments. Quantitative real-time polymerase chain reaction and Western blot assay results indicated that PRSS1 expression levels in platinum-resistant samples (tissue/cell) were significantly higher than in samples sensitive to platinum. By cell transfection assay, we observed that knockdown of PRSS1 reduced the resistance of ovarian cancer cells to cisplatin. Meanwhile, overexpression of PRSS1 increased the resistance to cisplatin. In conclusion, we identified a novel risk gene PRSS1 related to ovarian cancer platinum response and confirmed its key roles using multiple levels of low-throughput experiments, revealing a new treatment strategy based on a novel target factor for overcoming cisplatin resistance in ovarian cancer.

Identification of Functional lncRNAs Associated With Ovarian Endometriosis Based on a ceRNA Network

Background

Endometriosis is a common gynecological disease affecting women of reproductive age; however, the mechanisms underlying this condition are not fully clear. The aim of this study was to identify functional long non-coding RNAs (lncRNAs) associated with ovarian endometriosis for potential use as biomarkers and therapeutic targets.

Methods

RNA-seq profiles of paired ectopic (EC) and eutopic (EU) endometrial samples from patients with ovarian endometriosis were downloaded from the publicly available Gene Expression Omnibus (GEO) database. Bioinformatics algorithms were used to construct a network of ovarian endometriosis-related competing endogenous RNAs (ceRNAs) and to detect functional lncRNAs.

Results

A total of 4,213 mRNAs, 1,474 lncRNAs, and 221 miRNAs were identified as being differentially expressed between EC and EU samples, and an ovarian endometriosis-related ceRNA network was constructed through analysis of these differentially expressed RNAs. H19 and GS1-358P8.4 were identified as key ovarian endometriosis-related lncRNAs through topological feature analysis, and RP11-96D1.10 was identified using a random walk with restart algorithm.

Conclusion

Based on bioinformatics analysis of a ceRNA network, we identified the lncRNAs H19, GS1-358P8.4, and RP11-96D1.10 as being strongly associated with ovarian endometriosis. These three lncRNAs hold potential as targets for medical therapy and as diagnostic biomarkers. Further studies are needed to elucidate the detailed biological function of these lncRNAs in the pathogenesis of endometriosis.

CNA2Subpathway: identification of dysregulated subpathway driven by copy number alterations in cancer

Biological pathways reflect the key cellular mechanisms that dictate disease states, drug response and altered cellular function. The local areas of pathways are defined as subpathways (SPs), whose dysfunction has been reported to be associated with the occurrence and development of cancer. With the development of high-throughput sequencing technology, identifying dysfunctional SPs by using multi-omics data has become possible. Moreover, the SPs are not isolated in the biological system but interact with each other. Here, we propose a network-based calculated method, CNA2Subpathway, to identify dysfunctional SPs is driven by somatic copy number alterations (CNAs) in cancer through integrating pathway topology information, multi-omics data and SP crosstalk. This provides a novel way of SP analysis by using the SP interactions in the system biological level. Using data sets from breast cancer and head and neck cancer, we validate the effectiveness of CNA2Subpathway in identifying cancer-relevant SPs driven by the somatic CNAs, which are also shown to be associated with cancer immune and prognosis of patients. We further compare our results with five pathway or SP analysis methods based on CNA and gene expression data without considering SP crosstalk. With these analyses, we show that CNA2Subpathway could help to uncover dysfunctional SPs underlying cancer via the use of SP crosstalk. CNA2Subpathway is developed as an R-based tool, which is freely available on GitHub (https://github.com/hanjunwei-lab/CNA2Subpathway).

A network-based machine-learning framework to identify both functional modules and disease genes

Disease gene identification is a critical step towards uncovering the molecular mechanisms of diseases and systematically investigating complex disease phenotypes. Despite considerable efforts to develop powerful computing methods, candidate gene identification remains a severe challenge owing to the connectivity of an incomplete interactome network, which hampers the discovery of true novel candidate genes. We developed a network-based machine-learning framework to identify both functional modules and disease candidate genes. In this framework, we designed a semi-supervised non-negative matrix factorization model to obtain the functional modules related to the diseases and genes. Of note, we proposed a disease gene-prioritizing method called MapGene that integrates the correlations from both functional modules and network closeness. Our framework identified a set of functional modules with highly functional homogeneity and close gene interactions. Experiments on a large-scale benchmark dataset showed that MapGene performs significantly better than the state-of-the-art algorithms. Further analysis demonstrates MapGene can effectively relieve the impact of the incompleteness of interactome networks and obtain highly reliable rankings of candidate genes. In addition, disease cases on Parkinson's disease and diabetes mellitus confirmed the generalization of MapGene for novel candidate gene identification. This work proposed, for the first time, an integrated computing framework to predict both functional modules and disease candidate genes. The methodology and results support that our framework has the potential to help discover underlying functional modules and reliable candidate genes in human disease.

iMPTCE-Hnetwork: A Multilabel Classifier for Identifying Metabolic Pathway Types of Chemicals and Enzymes with a Heterogeneous Network

Metabolic pathway is an important type of biological pathways. It produces essential molecules and energies to maintain the life of living organisms. Each metabolic pathway consists of a chain of chemical reactions, which always need enzymes to participate in. Thus, chemicals and enzymes are two major components for each metabolic pathway. Although several metabolic pathways have been uncovered, the metabolic pathway system is still far from complete. Some hidden chemicals or enzymes are not discovered in a certain metabolic pathway. Besides the traditional experiments to detect hidden chemicals or enzymes, an alternative pipeline is to design efficient computational methods. In this study, we proposed a powerful multilabel classifier, called iMPTCE-Hnetwork, to uniformly assign chemicals and enzymes to metabolic pathway types reported in KEGG. Such classifier adopted the embedding features derived from a heterogeneous network, which defined chemicals and enzymes as nodes and the interactions between chemicals and enzymes as edges, through a powerful network embedding algorithm, Mashup. The popular RAndom k-labELsets (RAKEL) algorithm was employed to construct the classifier, which incorporated the support vector machine (polynomial kernel) as the basic classifier. The ten-fold cross-validation results indicated that such a classifier had good performance with accuracy higher than 0.800 and exact match higher than 0.750. Several comparisons were done to indicate the superiority of the iMPTCE-Hnetwork.

Complex Networks in Health and Disease

From protein interactions to signal transduction, from metabolism to the nervous system: Virtually all processes in health and disease rely on the careful orchestration of a large number of diverse individual components ranging from molecules to cells and entire organs. Networks provide a powerful framework for describing and understanding these complex systems in a wholistic fashion. They offer a unique combination of a highly intuitive, qualitative description, and a plethora of analytical, quantitative tools. Here we provide a brief introduction to the emerging field of network medicine. After an overview of the core concepts for connecting network characteristics to biological functions, we review commonly used networks, ranging from the molecular interaction networks that form the basis of all biological processes in the cell to the global transportation networks that govern the spread of global epidemics. Lastly, we highlight current conceptual and practical challenges.

Prioritizing Type 2 Diabetes Genes by Weighted PageRank on Bilayer Heterogeneous Networks

The prevalence of diabetes mellitus has been increasing rapidly in recent years. Type 2 diabetes makes up about 90 percent cases of diabetes. The interacting mixed effects of genetics and environments build possible interpretable pathogenesis. Thus, finding the causal disease genes is crucial in its clinical diagnosis and medical treatment. Currently, network-based computational method becomes a powerful tool of systematically analyzing complex diseases, such as the identification of candidate disease genes from networks. In this paper, we propose a bioinformatics framework of prioritizing type 2 diabetes genes by leveraging the modified PageRank algorithm on bilayer biomolecular networks consisting an ensemble gene-gene regulatory network and an integrative protein-protein interaction network. We specifically weigh the networks by differential mutual information for measuring the context specificities between genes and between proteins by transcriptomic and proteomic datasets, respectively. After formulating the network into two components of known disease genes and the other normal healthy genes, we rank the diabetes genes and others by bringing the orders in the bilayer network via an improved PageRank algorithm. We conclude that these known disease genes achieve significantly higher ranks compared to these randomly-selected normal genes, and the ranks are robust and consistent in multiple validation scenarios. In functional analysis, these high-ranked genes are identified to perform relevant risks and dysfunctions of type 2 diabetes.

Frequent Subgraph Mining of Functional Interaction Patterns Across Multiple Cancers

Molecular mechanisms characterizing cancer development and progression are complex and process through thousands of interacting elements in the cell. Understanding the underlying structure of interactions requires the integration of cellular networks with extensive combinations of dysregulation patterns. Recent pan-cancer studies focused on identifying common dysregulation patterns in a confined set of pathways or targeting a manually curated set of genes. However, the complex nature of the disease presents a challenge for finding pathways that would constitute a basis for tumor progression and requires evaluation of subnetworks with functional interactions. Uncovering these relationships is critical for translational medicine and the identification of future therapeutics. We present a frequent subgraph mining algorithm to find functional dysregulation patterns across the cancer spectrum. We mined frequent subgraphs coupled with biased random walks utilizing genomic alterations, gene expression profiles, and protein-protein interaction networks. In this unsupervised approach, we have recovered expert-curated pathways previously reported for explaining the underlying biology of cancer progression in multiple cancer types. Furthermore, we have clustered the genes identified in the frequent subgraphs into highly connected networks using a greedy approach and evaluated biological significance through pathway enrichment analysis. Gene clusters further elaborated on the inherent heterogeneity of cancer samples by both suggesting specific mechanisms for cancer type and common dysregulation patterns across different cancer types. Survival analysis of sample level clusters also revealed significant differences among cancer types (p < 0.001). These results could extend the current understanding of disease etiology by identifying biologically relevant interactions.Supplementary Information: Supplementary methods, figures, tables and code are available at https://github.com/bebeklab/FSM_Pancancer.

Pathogenic gene prediction based on network embedding

In disease research, the study of gene-disease correlation has always been an important topic. With the emergence of large-scale connected data sets in biology, we use known correlations between the entities, which may be from different sets, to build a biological heterogeneous network and propose a new network embedded representation algorithm to calculate the correlation between disease and genes, using the correlation score to predict pathogenic genes. Then, we conduct several experiments to compare our method to other state-of-the-art methods. The results reveal that our method achieves better performance than the traditional methods.

Cancer subtype classification and modeling by pathway attention and propagation

MOTIVATION

Biological pathway is an important curated knowledge of biological processes. Thus, cancer subtype classification based on pathways will be very useful to understand differences in biological mechanisms among cancer subtypes. However, pathways include only a fraction of the entire gene set, only one-third of human genes in KEGG, and pathways are fragmented. For this reason, there are few computational methods to use pathways for cancer subtype classification.

RESULTS

We present an explainable deep-learning model with attention mechanism and network propagation for cancer subtype classification. Each pathway is modeled by a graph convolutional network. Then, a multi-attention-based ensemble model combines several hundreds of pathways in an explainable manner. Lastly, network propagation on pathway-gene network explains why gene expression profiles in subtypes are different. In experiments with five TCGA cancer datasets, our method achieved very good classification accuracies and, additionally, identified subtype-specific pathways and biological functions.

AVAILABILITY AND IMPLEMENTATION

The source code is available at http://biohealth.snu.ac.kr/software/GCN_MAE.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Predicting human microbe-drug associations via graph convolutional network with conditional random field

Motivation

Human microbes play critical roles in drug development and precision medicine. How to systematically understand the complex interaction mechanism between human microbes and drugs remains a challenge nowadays. Identifying microbe-drug associations can not only provide great insights into understanding the mechanism, but also boost the development of drug discovery and repurposing. Considering the high cost and risk of biological experiments, the computational approach is an alternative choice. However, at present, few computational approaches have been developed to tackle this task.

Results

In this work, we leveraged rich biological information to construct a heterogeneous network for drugs and microbes, including a microbe similarity network, a drug similarity network, and a microbe-drug interaction network. We then proposed a novel Graph Convolutional Network (GCN) based framework for predicting human Microbe-Drug Associations, named GCNMDA. In the hidden layer of GCN, we further exploited the Conditional Random Field (CRF), which can ensure that similar nodes (i.e., microbes or drugs) have similar representations. To more accurately aggregate representations of neighborhoods, an attention mechanism was designed in the CRF layer. Moreover, we performed a random walk with restart (RWR) based scheme on both drug and microbe similarity networks to learn valuable features for drugs and microbes respectively. Experimental results on three different datasets showed that our GCNMDA model consistently achieved better performance than seven state-of-the-art methods. Case studies for three microbes including SARS-CoV-2 and two antimicrobial drugs (i.e., Ciprofloxacin and Moxifloxacin) further confirmed the effectiveness of GCNMDA in identifying potential microbe-drug associations.

Availability

Python codes and dataset are available at: https://github.com/longyahui/GCNMDA.

Supplementary information

Supplementary data are available at Bioinformatics online.

Recent advances in network-based methods for disease gene prediction

Disease-gene association through genome-wide association study (GWAS) is an arduous task for researchers. Investigating single nucleotide polymorphisms that correlate with specific diseases needs statistical analysis of associations. Considering the huge number of possible mutations, in addition to its high cost, another important drawback of GWAS analysis is the large number of false positives. Thus, researchers search for more evidence to cross-check their results through different sources. To provide the researchers with alternative and complementary low-cost disease-gene association evidence, computational approaches come into play. Since molecular networks are able to capture complex interplay among molecules in diseases, they become one of the most extensively used data for disease-gene association prediction. In this survey, we aim to provide a comprehensive and up-to-date review of network-based methods for disease gene prediction. We also conduct an empirical analysis on 14 state-of-the-art methods. To summarize, we first elucidate the task definition for disease gene prediction. Secondly, we categorize existing network-based efforts into network diffusion methods, traditional machine learning methods with handcrafted graph features and graph representation learning methods. Thirdly, an empirical analysis is conducted to evaluate the performance of the selected methods across seven diseases. We also provide distinguishing findings about the discussed methods based on our empirical analysis. Finally, we highlight potential research directions for future studies on disease gene prediction.

NECo: A node embedding algorithm for multiplex heterogeneous networks

Complex diseases such as hypertension, cancer, and diabetes cause nearly 70% of the deaths in the U.S. and involve multiple genes and their interactions with environmental factors. Therefore, identification of genetic factors to understand and decrease the morbidity and mortality from complex diseases is an important and challenging task. With the generation of an unprecedented amount of multi-omics datasets, network-based methods have become popular to represent the multilayered complex molecular interactions. Particularly node embeddings, the low-dimensional representations of nodes in a network are utilized for gene function prediction. Integrated network analysis of multi-omics data alleviates the issues related to missing data and lack of context-specific datasets. Most of the node embedding methods, however, are unable to integrate multiple types of datasets from genes and phenotypes. To address this limitation, we developed a node embedding algorithm called Node Embeddings of Complex networks (NECo) that can utilize multilayered heterogeneous networks of genes and phenotypes. We evaluated the performance of NECo using genotypic and phenotypic datasets from rat (Rattus norvegicus) disease models to classify hypertension disease-related genes. Our method significantly outperformed the state-of-the-art node embedding methods, with AUC of 94.97% compared 85.98% in the second-best performer, and predicted genes not previously implicated in hypertension.

Understanding oncogenicity of cancer driver genes and mutations in the cancer genomics era

One of the key challenges of cancer biology is to catalogue and understand the somatic genomic alterations leading to cancer. Although alternative definitions and search methods have been developed to identify cancer driver genes and mutations, analyses of thousands of cancer genomes return a remarkably similar catalogue of around 300 genes that are mutated in at least one cancer type. Yet, many features of these genes and their role in cancer remain unclear, first and foremost when a somatic mutation is truly oncogenic. In this review, we first summarize some of the recent efforts in completing the catalogue of cancer driver genes. Then, we give an overview of different aspects that influence the oncogenicity of somatic mutations in the core cancer driver genes, including their interactions with the germline genome, other cancer driver mutations, the immune system, or their potential role in healthy tissues. In the coming years, this research holds promise to illuminate how, when, and why cancer driver genes and mutations are really drivers, and thereby move personalized cancer medicine and targeted therapies forward.

Transcriptional insights into pathogenesis of cutaneous systemic sclerosis using pathway driven meta-analysis assisted by machine learning methods

Pathophysiology of systemic sclerosis (SSc, Scleroderma), an autoimmune rheumatic disease, comprises of mechanisms that drive vasculopathy, inflammation and fibrosis. Understanding of the disease and associated clinical heterogeneity has advanced considerably in the past decade, highlighting the necessity of more specific targeted therapy. While many of the recent trials in SSc failed to meet the primary end points that predominantly relied on changes in modified Rodnan skin scores (MRSS), sub-group analysis, especially those focused on the basal skin transcriptomic data have provided insights into patient subsets that respond to therapies. These findings suggest that deeper understanding of the molecular changes in pathways is very important to define disease drivers in various patient subgroups. In view of these challenges, we performed meta-analysis on 9 public available SSc microarray studies using a novel pathway pivoted approach combining consensus clustering and machine learning assisted feature selection. Selected pathway modules were further explored through cluster specific topological network analysis in search of novel therapeutic concepts. In addition, we went beyond previously described SSc class divisions of 3 clusters (e.g. inflammation, fibro-proliferative, normal-like) and expanded into a much finer stratification in order to profile SSc patients more accurately. Our analysis unveiled an important 80 pathway signatures that differentiated SSc patients into 8 unique subtypes. The 5 pathway modules derived from such signature successfully defined the 8 SSc subsets and were validated by in-silico cellular deconvolution analysis. Myeloid cells and fibroblasts involvement in different clusters were confirmed and linked to corresponding pathway activities. Collectively, our findings revealed more complex disease subtypes in SSc; Key gene mediators such as IL6, FGFR1, TLR7, PLCG2, IRK2 identified by network analysis underscored the scientific rationale for exploring additional targets in treatment of SSc.

A Deep Learning-Based Approach for Identifying the Medicinal Uses of Plant-Derived Natural Compounds

Medicinal plants and their extracts have been used as important sources for drug discovery. In particular, plant-derived natural compounds, including phytochemicals, antioxidants, vitamins, and minerals, are gaining attention as they promote health and prevent disease. Although several in vitro methods have been developed to confirm the biological activities of natural compounds, there is still considerable room to reduce time and cost. To overcome these limitations, several in silico methods have been proposed for conducting large-scale analysis, but they are still limited in terms of dealing with incomplete and heterogeneous natural compound data. Here, we propose a deep learning-based approach to identify the medicinal uses of natural compounds by exploiting massive and heterogeneous drug and natural compound data. The rationale behind this approach is that deep learning can effectively utilize heterogeneous features to alleviate incomplete information. Based on latent knowledge, molecular interactions, and chemical property features, we generated 686 dimensional features for 4,507 natural compounds and 2,882 approved and investigational drugs. The deep learning model was trained using the generated features and verified drug indication information. When the features of natural compounds were applied as input to the trained model, potential efficacies were successfully predicted with high accuracy, sensitivity, and specificity.

NetCore: a network propagation approach using node coreness

We present NetCore, a novel network propagation approach based on node coreness, for phenotype–genotype associations and module identification. NetCore addresses the node degree bias in PPI networks by using node coreness in the random walk with restart procedure, and achieves improved re-ranking of genes after propagation. Furthermore, NetCore implements a semi-supervised approach to identify phenotype-associated network modules, which anchors the identification of novel candidate genes at known genes associated with the phenotype. We evaluated NetCore on gene sets from 11 different GWAS traits and showed improved performance compared to the standard degree-based network propagation using cross-validation. Furthermore, we applied NetCore to identify disease genes and modules for Schizophrenia GWAS data and pan-cancer mutation data. We compared the novel approach to existing network propagation approaches and showed the benefits of using NetCore in comparison to those. We provide an easy-to-use implementation, together with a high confidence PPI network extracted from ConsensusPathDB, which can be applied to various types of genomics data in order to obtain a re-ranking of genes and functionally relevant network modules.

Connectivity Significance for Disease Gene Prioritization in an Expanding Universe

A fundamental topic in network medicine is disease genes prioritization. The underlying hypothesis is that disease genes are organized as modules confined within the interactome. Here, we propose a novel algorithm called DiaBLE (DIAMOnD Background Local Expansion) which is a modified version of DIAMOnD, a successful algorithm based on the concept of connectivity significance. Instead of taking the whole interactome as the background model, DiaBLE considers as gene universe the smallest local expansion of the current seeds set at each iteration step. We show that DiaBLE significantly increases the overall DIAMOnD ranking quality of genes prioritization both in terms of cross-validation and biological consistency. Here, we focus on the two algorithms only since a comparative analysis among gene prioritization methods is beyond the scope of this study. Finally, we briefly discuss the improvement of biological insight provided by DiaBLE for two cancers (head and neck squamous cell carcinoma and kidney renal clear cell carcinoma).

Supervised learning is an accurate method for network-based gene classification

Background

Assigning every human gene to specific functions, diseases and traits is a grand challenge in modern genetics. Key to addressing this challenge are computational methods, such as supervised learning and label propagation, that can leverage molecular interaction networks to predict gene attributes. In spite of being a popular machine-learning technique across fields, supervised learning has been applied only in a few network-based studies for predicting pathway-, phenotype- or disease-associated genes. It is unknown how supervised learning broadly performs across different networks and diverse gene classification tasks, and how it compares to label propagation, the widely benchmarked canonical approach for this problem.

Results

In this study, we present a comprehensive benchmarking of supervised learning for network-based gene classification, evaluating this approach and a classic label propagation technique on hundreds of diverse prediction tasks and multiple networks using stringent evaluation schemes. We demonstrate that supervised learning on a gene’s full network connectivity outperforms label propagaton and achieves high prediction accuracy by efficiently capturing local network properties, rivaling label propagation’s appeal for naturally using network topology. We further show that supervised learning on the full network is also superior to learning on node embeddings (derived using node2vec), an increasingly popular approach for concisely representing network connectivity. These results show that supervised learning is an accurate approach for prioritizing genes associated with diverse functions, diseases and traits and should be considered a staple of network-based gene classification workflows.

Availability and implementation

The datasets and the code used to reproduce the results and add new gene classification methods have been made freely available.

Contact

arjun@msu.edu

Supplementary information

Supplementary data are available at Bioinformatics online.

DiSNEP: a Disease-Specific gene Network Enhancement to improve Prioritizing candidate disease genes

Biological network-based strategies are useful in prioritizing genes associated with diseases. Several comprehensive human gene networks such as STRING, GIANT and HumanNet were developed and used in network-assisted algorithms to identify disease-associated genes. However, none of these networks are disease-specific and may not accurately reflect gene interactions for a specific disease. Aiming to improve disease gene prioritization using networks, we propose a Disease-Specific Network Enhancement Prioritization (DiSNEP) framework. DiSNEP first enhances a comprehensive gene network specifically for a disease through a diffusion process on a gene-gene similarity matrix derived from disease omics data. The enhanced disease-specific gene network thus better reflects true gene interactions for the disease and may improve prioritizing disease-associated genes subsequently. In simulations, DiSNEP that uses an enhanced disease-specific network prioritizes more true signal genes than comparison methods using a general gene network or without prioritization. Applications to prioritize cancer-associated gene expression and DNA methylation signal genes for five cancer types from The Cancer Genome Atlas (TCGA) project suggest that more prioritized candidate genes by DiSNEP are cancer-related according to the DisGeNET database than those prioritized by the comparison methods, consistently across all five cancer types considered, and for both gene expression and DNA methylation signal genes.

Using whole-exome sequencing and protein interaction networks to prioritize candidate genes for germline cutaneous melanoma susceptibility

Although next-generation sequencing has demonstrated great potential for novel gene discovery, confirming disease-causing genes after initial discovery remains challenging. Here, we applied a network analysis approach to prioritize candidate genes identified from whole-exome sequencing analysis of 98 cutaneous melanoma patients from 27 families. Using a network propagation method, we ranked candidate genes by their similarity to known disease genes in protein-protein interaction networks and identified gene clusters with functional connectivity. Using this approach, we identified several new candidate susceptibility genes that warrant future investigations such as NGLY1, IL1RN, FABP2, PRKDC, and PROSER2. The propagated network analysis also allowed us to link families that did not have common underlying genes but that carried variants in genes that interact on protein-protein interaction networks. In conclusion, our study provided an analysis perspective for gene prioritization in the context of genetic heterogeneity across families and prioritized top potential candidate susceptibility genes in our dataset.

Identification of Latent Oncogenes with a Network Embedding Method and Random Forest

Oncogene is a special type of genes, which can promote the tumor initiation. Good study on oncogenes is helpful for understanding the cause of cancers. Experimental techniques in early time are quite popular in detecting oncogenes. However, their defects become more and more evident in recent years, such as high cost and long time. The newly proposed computational methods provide an alternative way to study oncogenes, which can provide useful clues for further investigations on candidate genes. Considering the limitations of some previous computational methods, such as lack of learning procedures and terming genes as individual subjects, a novel computational method was proposed in this study. The method adopted the features derived from multiple protein networks, viewing proteins in a system level. A classic machine learning algorithm, random forest, was applied on these features to capture the essential characteristic of oncogenes, thereby building the prediction model. All genes except validated oncogenes were ranked with a measurement yielded by the prediction model. Top genes were quite different from potential oncogenes discovered by previous methods, and they can be confirmed to become novel oncogenes. It was indicated that the newly identified genes can be essential supplements for previous results.

RWR-algorithm-based dissection of microRNA-506-3p and microRNA-140-5p as radiosensitive biomarkers in colorectal cancer

Radiotherapy resistance is one of the main causes for treatment failure in colorectal cancer (CRC), and it is associated with the deregulation of certain microRNAs. In this study, we constructed the microRNA-mRNA network consisting of 2275 microRNAs and 7045 target genes, collected the known microRNAs related to CRC-radiosensitivity (CRCR) (n=18) as the seed nodes, and applied the algorithm of random walk with restart (RWR) to the network to identify novel CRCR-related microRNAs (n=263). In functional analysis, 263 novel microRNAs shared a high proportion of the same biological processes and pathways with the known microRNAs. In topological analysis of the sub-network of the 263 microRNAs and their targets, hsa-mir-506-3p and hsa-mir-140-5p were identified as network hub nodes. In plasma, radiosensitive patients had a higher expression level of hsa-mir-506-3p and hsa-mir-140-5p than radioresistant patients. In experimental validation, both hsa-mir-506-3p and hsa-mir-140-5p over-expression could obviously decrease the cell proliferation, survival rate and colonality in CRC cells after radiation. In conclusion, this study combined the novel network-based method with experimental validation, and identified two novel radiosensitive biomarkers of hsa-mir-506-3p and hsa-mir-140-5p in CRC.

Integrated Genomic Analysis Revealed Associated Genes for Alzheimer's Disease in APOE4 Non-Carriers

BACKGROUND

APOE4 is the strongest genetic risk factor for late-onset Alzheimer's disease (LOAD). LOAD patients carrying or not carrying APOE4 manifest distinct clinico-pathological characteristics. APOE4 has been shown to play a critical role in the pathogenesis of AD by affecting various aspects of pathological processes. However, the pathogenesis involved in LOAD not-carrying APOE4 remains elusive.

OBJECTIVE

We aimed to identify the associated genes involved in LOAD not-carrying APOE4.

METHODS

An integrated genomic analysis of datasets of genome-wide association study, genome-wide expression profiling and genome-wide linkage scan and protein-protein interaction network construction were applied to identify associated gene clusters in APOE4 non-carriers. The role of one of hub gene of an APOE4 non-carrier-associated gene cluster in tau phosphorylation was studied by knockdown and western blot.

RESULTS

We identified 12 gene clusters associated with AD APOE4 non-carriers. The hub genes associated with AD in these clusters were MAPK8, POU2F1, XRCC1, PRKCG, EXOC6, VAMP4, SIRT1, MME, NOS1, ABCA1 and LDLR. The associated genes for APOE4 non-carriers were enriched in hereditary disorder, neurological disease and psychological disorders. Moreover, knockdown of PRKCG to reduce the expression of protein kinase Cγ isoform enhanced tau phosphorylation at Thr181 and Thr231 and the expression of glycogen synthase kinase 3β and cyclin-dependent kinase 5 in the presence of APOE3 but not APOE4.

CONCLUSION

The study provides new insight into the mechanism of distinct pathogenesis of LOAD not carrying APOE4 and prompts the functional exploration of identified genes based on APOE genotypes.

The support of genetic evidence for cardiovascular risk induced by antineoplastic drugs

Cardiovascular dysfunction is one of the most common complications of long-term cancer treatment. Growing evidence has shown that antineoplastic drugs can increase cardiovascular risk during cancer therapy, seriously affecting patient survival. However, little is known about the genetic factors associated with the cardiovascular risk of antineoplastic drugs. We established a compendium of genetic evidence that supports cardiovascular risk induced by antineoplastic drugs. Most of this genetic evidence is attributed to causal alleles altering the expression of cardiovascular disease genes. We found that antineoplastic drugs predicted to induce cardiovascular risk are significantly enriched in drugs associated with cardiovascular adverse reactions, including many first-line cancer treatments. Functional experiments validated that retinoid X receptor agonists can reduce triglyceride lipolysis, thus modulating cardiovascular risk. Our results establish a link between the causal allele of cardiovascular disease genes and the direction of pharmacological modulation, which could facilitate cancer drug discovery and clinical trial design.

Computational/in silico methods in drug target and lead prediction

Drug-like compounds are most of the time denied approval and use owing to the unexpected clinical side effects and cross-reactivity observed during clinical trials. These unexpected outcomes resulting in significant increase in attrition rate centralizes on the selected drug targets. These targets may be disease candidate proteins or genes, biological pathways, disease-associated microRNAs, disease-related biomarkers, abnormal molecular phenotypes, crucial nodes of biological network or molecular functions. This is generally linked to several factors, including incomplete knowledge on the drug targets and unpredicted pharmacokinetic expressions upon target interaction or off-target effects. A method used to identify targets, especially for polygenic diseases, is essential and constitutes a major bottleneck in drug development with the fundamental stage being the identification and validation of drug targets of interest for further downstream processes. Thus, various computational methods have been developed to complement experimental approaches in drug discovery. Here, we present an overview of various computational methods and tools applied in predicting or validating drug targets and drug-like molecules. We provide an overview on their advantages and compare these methods to identify effective methods which likely lead to optimal results. We also explore major sources of drug failure considering the challenges and opportunities involved. This review might guide researchers on selecting the most efficient approach or technique during the computational drug discovery process.

Applications of Genome-Wide Screening and Systems Biology Approaches in Drug Repositioning

Modern drug discovery through de novo drug discovery entails high financial costs, low success rates, and lengthy trial periods. Drug repositioning presents a suitable approach for overcoming these issues by re-evaluating biological targets and modes of action of approved drugs. Coupling high-throughput technologies with genome-wide essentiality screens, network analysis, genome-scale metabolic modeling, and machine learning techniques enables the proposal of new drug-target signatures and uncovers unanticipated modes of action for available drugs. Here, we discuss the current issues associated with drug repositioning in light of curated high-throughput multi-omic databases, genome-wide screening technologies, and their application in systems biology/medicine approaches.

iATC-NRAKEL: an efficient multi-label classifier for recognizing anatomical therapeutic chemical classes of drugs

MOTIVATION

The anatomical therapeutic chemical (ATC) classification system plays an increasingly important role in drug repositioning and discovery. The correct identification of classes in each level of such system that a given drug may belong to is an essential problem. Several multi-label classifiers have been proposed in this regard. Although they provided satisfactory performance, the feature extraction procedures were still rough. More refined features may further improve the predicted quality.

RESULTS

In this article, we provide a novel multi-label classifier, called iATC-NRAKEL, to predict drug ATC classes in the first level. To obtain more informative drug features, we employed the drug association information in STITCH and KEGG, which was organized by seven drug networks. The powerful network embedding algorithm, Mashup, was adopted to extract informative drug features. The obtained features were fed into the RAndom k-labELsets (RAKEL) algorithm with support vector machine as the basic classification algorithm to construct the classifier. The 10-fold cross-validation of the benchmark dataset with 3883 drugs showed that the accuracy and absolute true were 76.56 and 74.51%, respectively. The comparison results indicated that iATC-NRAKEL was much superior to all previous reported classifiers. Finally, the contribution of each network was analyzed.

AVAILABILITY AND IMPLEMENTATION

The codes of iATC-NRAKEL are available at https://github.com/zhou256/iATC-NRAKEL.

Deep Collaborative Filtering for Prediction of Disease Genes

Accurate prioritization of potential disease genes is a fundamental challenge in biomedical research. Various algorithms have been developed to solve such problems. Inductive Matrix Completion (IMC) is one of the most reliable models for its well-established framework and its superior performance in predicting gene-disease associations. However, the IMC method does not hierarchically extract deep features, which might limit the quality of recovery. In this case, the architecture of deep learning, which obtains high-level representations and handles noises and outliers presented in large-scale biological datasets, is introduced into the side information of genes in our Deep Collaborative Filtering (DCF) model. Further, for lack of negative examples, we also exploit Positive-Unlabeled (PU) learning formulation to low-rank matrix completion. Our approach achieves substantially improved performance over other state-of-the-art methods on diseases from the Online Mendelian Inheritance in Man (OMIM) database. Our approach is 10 percent more efficient than standard IMC in detecting a true association, and significantly outperforms other alternatives in terms of the precision-recall metric at the top-k predictions. Moreover, we also validate the disease with no previously known gene associations and newly reported OMIM associations. The experimental results show that DCF is still satisfactory for ranking novel disease phenotypes as well as mining unexplored relationships. The source code and the data are available at https://github.com/xzenglab/DCF.

Prediction of membrane protein types by fusing protein-protein interaction and protein sequence information

Membrane proteins are gatekeepers to the cell and essential for determination of the function of cells. Identification of the types of membrane proteins is an essential problem in cell biology. It is time-consuming and expensive to identify the type of membrane proteins with traditional experimental methods. The alternative way is to design effective computational methods, which can provide quick and reliable predictions. To date, several computational methods have been proposed in this regard. Several of them used the features extracted from the sequence information of individual proteins. Recently, networks are more and more popular to tackle different protein-related problems, which can organize proteins in a system level and give an overview of all proteins. However, such form weakens the essential properties of proteins, such as their sequence information. In this study, a novel feature fusion scheme was proposed, which integrated the information of protein sequences and protein-protein interaction network. The fused features of a protein were defined as the linear combination of sequence features of all proteins in the network, where the combination coefficients were the probabilities yielded by the random walk with restart algorithm with the protein as the seed node. Several models with such fused features and different classification algorithms were built and evaluated. Their performance for predicting the type of membrane proteins was improved compared with the models only with the sequence features or network information.

Integrated Genomic Characterization of the Human Immunome in Cancer

Alterations in immune-related pathways are common hallmarks of cancer. A comprehensive understanding of how cancer mutations rewire immune signaling networks and functional output across cancer types is instrumental to realize the full potential of immunotherapy. Here, we systematically interrogated somatic mutations involved in immune signaling that alter immune responses in patients with cancer. To do so, we developed a Network-based Integrative model to Prioritize Potential immune respondER genes (NIPPER). Identified mutations were enriched in essential protein domains and genes identified by NIPPER were associated with responsiveness to multiple immunotherapy modalities. These genes were used to devise an interactome network propagation framework integrated with drug-associated gene signatures to identify potential immunomodulatory drug candidates. Together, our systems-level analysis results help interpret the heterogeneous immune responses among patients and serve as a resource for future functional studies and targeted therapeutics. SIGNIFICANCE: This study demonstrates that integration of multi-omics data can help identify critical molecular determinants for effective targeted therapeutics.

Prediction of drug response in multilayer networks based on fusion of multiomics data

Predicting the response of each individual patient to a drug is a key issue assailing personalized medicine. Our study predicted drug response based on the fusion of multiomics data with low-dimensional feature vector representation on a multilayer network model. We named this new method DREMO (Drug Response prEdiction based on MultiOmics data fusion). DREMO fuses similarities between cell lines and similarities between drugs, thereby improving the ability to predict the response of cancer cell lines to therapeutic agents. First, a multilayer similarity network related to cell lines and drugs was constructed based on gene expression profiles, somatic mutation, copy number variation (CNV), drug chemical structures, and drug targets. Next, low-dimensional feature vector representation was used to fuse the biological information in the multilayer network. Then, a machine learning model was applied to predict new drug-cell line associations. Finally, our results were validated using the well-established GDSC/CCLE databases, literature, and the functional pathway database. Furthermore, a comparison was made between DREMO and other methods. Results of the comparison showed that DREMO improves predictive capabilities significantly.

The identification of six risk genes for ovarian cancer platinum response based on global network algorithm and verification analysis

Ovarian cancer is the most lethal gynaecological cancer, and resistance of platinum‐based chemotherapy is the main reason for treatment failure. The aim of the present study was to identify candidate genes involved in ovarian cancer platinum response by analysing genes from homologous recombination and Fanconi anaemia pathways. Associations between these two functional genes were explored in the study, and we performed a random walk algorithm based on reconstructed gene‐gene network, including protein‐protein interaction and co‐expression relations. Following the random walk, all genes were ranked and GSEA analysis showed that the biological functions focused primarily on autophagy, histone modification and gluconeogenesis. Based on three types of seed nodes, the top two genes were utilized as examples. We selected a total of six candidate genes (FANCA, FANCG, POLD1, KDM1A, BLM and BRCA1) for subsequent verification. The validation results of the six candidate genes have significance in three independent ovarian cancer data sets with platinum‐resistant and platinum‐sensitive information. To explore the correlation between biomarkers and clinical prognostic factors, we performed differential analysis and multivariate clinical subgroup analysis for six candidate genes at both mRNA and protein levels. And each of the six candidate genes and their neighbouring genes with a mutation rate greater than 10% were also analysed by network construction and functional enrichment analysis. In the meanwhile, the survival analysis for platinum‐treated patients was performed in the current study. Finally, the RT‐qPCR assay was used to determine the performance of candidate genes in ovarian cancer platinum response. Taken together, this research demonstrated that comprehensive bioinformatics methods could help to understand the molecular mechanism of platinum response and provide new strategies for overcoming platinum resistance in ovarian cancer treatment.

Interactome overlap between schizophrenia and cognition

Cognitive impairments constitute a core feature of schizophrenia, and a genetic overlap between schizophrenia and cognitive functioning in healthy individuals has been identified. However, due to the high polygenicity and complex genetic architecture of both traits, overlapping biological pathways have not yet been identified between schizophrenia and normal cognitive ability. Network medicine offers a framework to study underlying biological pathways through protein-protein interactions among risk genes. Here, established network-based methods were used to characterize the biological relatedness of schizophrenia and cognition by examining the genetic link between schizophrenia risk genes and genes associated with cognitive performance in healthy individuals, through the protein interactome. First, network separation showed a profound interactome overlap between schizophrenia risk genes and genes associated with cognitive performance (S_AB = -0.22, z-score = -6.80, p = 5.38e-12). To characterize this overlap, network propagation was thereafter used to identify schizophrenia risk genes that are close to cognition-associated genes in the interactome network space (n = 140, of which 54 were part of the direct genetic overlap). Schizophrenia risk genes close to cognition were enriched for pathways including long-term potentiation and Alzheimer's disease, and included genes with a role in neurotransmitter systems important for cognitive functioning, such as glutamate and dopamine. These results pinpoint a subset of schizophrenia risk genes that are of particular interest for further examination in schizophrenia patient groups, of which some are druggable genes with potential as candidate targets for cognitive enhancing drugs.

Identification of Differentially Expressed Genes and miRNAs Associated with Esophageal Squamous Cell Carcinoma by Integrated Analysis of Microarray Data

To identify candidate key genes and miRNAs associated with esophageal squamous cell carcinoma (ESCC) development and prognosis, the gene expression profiles and miRNA microarray data including GSE20347, GSE38129, GSE23400, and GSE55856 were downloaded from the Gene Expression Omnibus (GEO) database. Clinical and survival data were retrieved from The Cancer Genome Atlas (TCGA). Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of differentially expressed genes (DEGs) was analyzed via DAVID, while the DEG-associated protein-protein interaction network (PPI) was constructed using the STRING database. Additionally, the miRNA target gene regulatory network and miRNA coregulatory network were constructed, using the Cytoscape software. Survival analysis and prognostic model construction were performed via the survival (version 2.42-6) and rbsurv R packages, respectively. The results showed a total of 2575, 2111, and 1205 DEGs, and 226 differentially expressed miRNAs (DEMs) were identified. Pathway enrichment analyses revealed that DEGs were mainly enriched in 36 pathways, such as the proteasome, p53, and beta-alanine metabolism pathways. Furthermore, 448 nodes and 1144 interactions were identified in the PPI network, with MYC having the highest random walk score. In addition, 7 DEMs in the microarray data, including miR-196a, miR-21, miR-205, miR-194, miR-103, miR-223, and miR-375, were found in the regulatory network. Moreover, several reported disease-related miRNAs, including miR-198a, miR-103, miR-223, miR-21, miR-194, and miR-375, were found to have common target genes with other DEMs. Survival analysis revealed that 85 DEMs were related to prognosis, among which hsa-miR-1248, hsa-miR-1291, hsa-miR-421, and hsa-miR-7-5p were used for a prognostic survival model. Taken together, this study revealed the important roles of DEGs and DEMs in ESCC development, as well as DEMs in the prognosis of ESCC. This will provide potential therapeutic targets and prognostic predictors for ESCC.

Multi-scale top-down approach for modelling epileptic protein-protein interaction network analysis to identify driver nodes and pathways

Protein - Protein Interaction Network (PPIN) analysis unveils molecular level mechanisms involved in disease condition. To explore the complex regulatory mechanisms behind epilepsy and to address the clinical and biological issues of epilepsy, in silico techniques are feasible in a cost- effective manner. In this work, a hierarchical procedure to identify influential genes and regulatory pathways in epilepsy prognosis is proposed. To obtain key genes and pathways causing epilepsy, integration of two benchmarked datasets which are exclusively devoted for complex disorders is done as an initial step. Using STRING database, PPIN is constructed for modelling protein-protein interactions. Further, key interactions are obtained from the established PPIN using network centrality measures followed by network propagation algorithm -Random Walk with Restart (RWR). The outcome of the method reveals some influential genes behind epilepsy prognosis, along with their associated pathways like PI3 kinase, VEGF signaling, Ras, Wnt signaling etc. In comparison with similar works, our results have shown improvement in identifying unique molecular functions, biological processes, gene co-occurrences etc. Also, CORUM provides an annotation for approximately 60% of similarity in human protein complexes with the obtained result. We believe that the formulated strategy can put-up the vast consideration of indigenous drugs towards meticulous identification of genes encoded by protein against several combinatorial disorders.

Identification of COVID-19 Infection-Related Human Genes Based on a Random Walk Model in a Virus-Human Protein Interaction Network

Coronaviruses are specific crown-shaped viruses that were first identified in the 1960s, and three typical examples of the most recent coronavirus disease outbreaks include severe acute respiratory syndrome (SARS), Middle East respiratory syndrome (MERS), and COVID-19. Particularly, COVID-19 is currently causing a worldwide pandemic, threatening the health of human beings globally. The identification of viral pathogenic mechanisms is important for further developing effective drugs and targeted clinical treatment methods. The delayed revelation of viral infectious mechanisms is currently one of the technical obstacles in the prevention and treatment of infectious diseases. In this study, we proposed a random walk model to identify the potential pathological mechanisms of COVID-19 on a virus–human protein interaction network, and we effectively identified a group of proteins that have already been determined to be potentially important for COVID-19 infection and for similar SARS infections, which help further developing drugs and targeted therapeutic methods against COVID-19. Moreover, we constructed a standard computational workflow for predicting the pathological biomarkers and related pharmacological targets of infectious diseases.

deepDR: a network-based deep learning approach to in silico drug repositioning

MOTIVATION

Traditional drug discovery and development are often time-consuming and high risk. Repurposing/repositioning of approved drugs offers a relatively low-cost and high-efficiency approach toward rapid development of efficacious treatments. The emergence of large-scale, heterogeneous biological networks has offered unprecedented opportunities for developing in silico drug repositioning approaches. However, capturing highly non-linear, heterogeneous network structures by most existing approaches for drug repositioning has been challenging.

RESULTS

In this study, we developed a network-based deep-learning approach, termed deepDR, for in silico drug repurposing by integrating 10 networks: one drug-disease, one drug-side-effect, one drug-target and seven drug-drug networks. Specifically, deepDR learns high-level features of drugs from the heterogeneous networks by a multi-modal deep autoencoder. Then the learned low-dimensional representation of drugs together with clinically reported drug-disease pairs are encoded and decoded collectively via a variational autoencoder to infer candidates for approved drugs for which they were not originally approved. We found that deepDR revealed high performance [the area under receiver operating characteristic curve (AUROC) = 0.908], outperforming conventional network-based or machine learning-based approaches. Importantly, deepDR-predicted drug-disease associations were validated by the ClinicalTrials.gov database (AUROC = 0.826) and we showcased several novel deepDR-predicted approved drugs for Alzheimer's disease (e.g. risperidone and aripiprazole) and Parkinson's disease (e.g. methylphenidate and pergolide).

AVAILABILITY AND IMPLEMENTATION

Source code and data can be downloaded from https://github.com/ChengF-Lab/deepDR.

SUPPLEMENTARY INFORMATION

Supplementary data are available online at Bioinformatics.

Stability Analysis of Biological Networks' Diffusion State

Computational knowledge acquired from noisy networks is not reliable and the network topology determines the reliability. Protein-protein interaction networks have uncertain topologies and noise that contain false positive and false negative edges at high rates. In this study, we analyze effects of the existing mutations in a network topology to the diffusion state of that network. To evaluate the sensitivity of the diffusion state, we derive the fitness measures based on the mathematically defined stability of a network. Searching for an influential set of edges in a network is a difficult problem. We handle the computational challenge by developing a novel metaheuristic optimization method and we find influential mutations time-efficiently. Our experiments, conducted on both synthetic and real networks from public databases, demonstrated that our method obtained better results than competitors for all types of network topologies. This is the first-time that the diffusion has been evaluated under topological mutations. Our analysis identifies significant biological results about the stability of biological - synthetic networks and diffusion state. In this manner, mutations in protein-protein interaction network topologies have a significant influence on the diffusion state of the network. Network stability is more affected by the network model than the network size.

PrGeFNE: Predicting disease-related genes by fast network embedding

Identifying disease-related genes is of importance for understanding of molecule mechanisms of diseases, as well as diagnosis and treatment of diseases. Many computational methods have been proposed to predict disease-related genes, but how to make full use of multi-source biological data to enhance the ability of disease-gene prediction is still challenging. In this paper, we proposed a novel method for predicting disease-related genes by using fast network embedding (PrGeFNE), which can integrate multiple types of associations related to diseases and genes. Specifically, we first constructed a heterogeneous network by using phenotype-disease, disease-gene, protein-protein and gene-GO associations; and low-dimensional representation of nodes is extracted from the network by using a fast network embedding algorithm. Then, a dual-layer heterogeneous network was reconstructed by using the low-dimensional representation, and a network propagation was applied to the dual-layer heterogeneous network to predict disease-related genes. Through cross-validation and newly added-association validation, we displayed the important roles of different types of association data in enhancing the ability of disease-gene prediction, and confirmed the excellent performance of PrGeFNE by comparing to state-of-the-art algorithms. Furthermore, we developed a web tool that can facilitate researchers to search for candidate genes of different diseases predicted by PrGeFNE, along with the enrichment analysis of GO and pathway on candidate gene set. This may be useful for investigation of diseases' molecular mechanisms as well as their experimental validations. The web tool is available at http://bioinformatics.csu.edu.cn/prgefne/.

uKIN Combines New and Prior Information with Guided Network Propagation to Accurately Identify Disease Genes

Protein interaction networks provide a powerful framework for identifying genes causal for complex genetic diseases. Here, we introduce a general framework, uKIN, that uses prior knowledge of disease-associated genes to guide, within known protein-protein interaction networks, random walks that are initiated from newly identified candidate genes. In large-scale testing across 24 cancer types, we demonstrate that our network propagation approach for integrating both prior and new information not only better identifies cancer driver genes than using either source of information alone but also readily outperforms other state-of-the-art network-based approaches. We also apply our approach to genome-wide association data to identify genes functionally relevant for several complex diseases. Overall, our work suggests that guided network propagation approaches that utilize both prior and new data are a powerful means to identify disease genes. uKIN is freely available for download at: https://github.com/Singh-Lab/uKIN.

A Knowledge-Based Machine Learning Approach to Gene Prioritisation in Amyotrophic Lateral Sclerosis

Amyotrophic lateral sclerosis is a neurodegenerative disease of the upper and lower motor neurons resulting in death from neuromuscular respiratory failure, typically within two to five years of first symptoms. Several rare disruptive gene variants have been associated with ALS and are responsible for about 15% of all cases. Although our knowledge of the genetic landscape of this disease is improving, it remains limited. Machine learning models trained on the available protein-protein interaction and phenotype-genotype association data can use our current knowledge of the disease genetics for the prediction of novel candidate genes. Here, we describe a knowledge-based machine learning method for this purpose. We trained our model on protein-protein interaction data from IntAct, gene function annotation from Gene Ontology, and known disease-gene associations from DisGeNet. Using several sets of known ALS genes from public databases and a manual review as input, we generated a list of new candidate genes for each input set. We investigated the relevance of the predicted genes in ALS by using the available summary statistics from the largest ALS genome-wide association study and by performing functional and phenotype enrichment analysis. The predicted sets were enriched for genes associated with other neurodegenerative diseases known to overlap with ALS genetically and phenotypically, as well as for biological processes associated with the disease. Moreover, using ALS genes from ClinVar and our manual review as input, the predicted sets were enriched for ALS-associated genes (ClinVar p = 0.038 and manual review p = 0.060) when used for gene prioritisation in a genome-wide association study.

Exploring the mechanism of TCM formulae in the treatment of different types of coronary heart disease by network pharmacology and machining learning

Traditional Chinese medicine (TCM) has long been used in the clinical treatment of coronary heart disease (CHD). TCM is characterized by syndrome-based medication, which is, using different TCM formulae for different syndromes. However, the underlying mode of action remains unclear. In this work, we utilized network pharmacology and machine learning to explore the mechanism of eight classic TCM formulae in the treatment of different types of CHD. First, by integrating multiple databases, a total of 669 potential bioactive compounds and 581 targets of the eight formulae were screened. Then, the effectiveness of these formulae on CHD was evaluated using two network-based indicators. The results showed that each formula's targets were significantly correlated with CHD associated genes and overlapped with the targets of 9 classes of drugs for cardio vascular diseases (CVD) to some degree. Next, from 5 different levels, i.e., herb, symptom, compound, target, and pathway level, we systematically compared the eight formulae using network clustering and hierarchical clustering. We found that all the formulae could be grouped into five clusters and the clustering results were approximately consistent at different levels. All the formulae were involved in 7 pathways closely related to CHD and may exhibit the common effect of relieving angina. Formulae in the same group collectively regulated some unique pathways and suggest further specific indications. For example, the three formulae used for Qi stagnation and blood stasis, Qi deficiency and blood stasis, and Qi-Yin deficiency syndromes acted on two special pathways (TNF signaling pathway, NF-kappa B signaling pathway) and may exert anti-inflammatory and immune-enhancing effects; the two formulae for Yin deficiency of heart and kidney, and Yang deficiency of heart and kidney syndromes regulated two special pathways (PPAR signaling pathway, thyroid hormone signaling pathway) in endocrine system and could improve renal function. Subsequently, we designed a rank algorithm, which integrated network topology with biological function, to identify important targets of these formulae. The results were consistent with the multi-level clustering results. At last, our literature mining validated about 20 % putative targets, as well as clustering results and effects of the formulae by experimental evidences. This study explained the medication patterns and scientific significance of TCM formulae on different types of CHD from perspective of systems biology. It may facilitate the understanding of different types of CHD described by traditional Chinese medicine from the perspectives of modern biology.

Dynamic rewiring of the human interactome by interferon signaling

BACKGROUND

The type I interferon (IFN) response is an ancient pathway that protects cells against viral pathogens by inducing the transcription of hundreds of IFN-stimulated genes. Comprehensive catalogs of IFN-stimulated genes have been established across species and cell types by transcriptomic and biochemical approaches, but their antiviral mechanisms remain incompletely characterized. Here, we apply a combination of quantitative proteomic approaches to describe the effects of IFN signaling on the human proteome, and apply protein correlation profiling to map IFN-induced rearrangements in the human protein-protein interaction network.

RESULTS

We identify > 26,000 protein interactions in IFN-stimulated and unstimulated cells, many of which involve proteins associated with human disease and are observed exclusively within the IFN-stimulated network. Differential network analysis reveals interaction rewiring across a surprisingly broad spectrum of cellular pathways in the antiviral response. We identify IFN-dependent protein-protein interactions mediating novel regulatory mechanisms at the transcriptional and translational levels, with one such interaction modulating the transcriptional activity of STAT1. Moreover, we reveal IFN-dependent changes in ribosomal composition that act to buffer IFN-stimulated gene protein synthesis.

CONCLUSIONS

Our map of the IFN interactome provides a global view of the complex cellular networks activated during the antiviral response, placing IFN-stimulated genes in a functional context, and serves as a framework to understand how these networks are dysregulated in autoimmune or inflammatory disease.

Identification of Conserved Gene-Regulatory Networks that Integrate Environmental Sensing and Growth in the Root Cambium

Cambium drives the lateral growth of stems and roots, contributing to diverse plant growth forms. The root crop is one of the outstanding examples of the cambium-driven growth. To understand its molecular basis, we used radish to generate a compendium of root-tissue- and stage-specific transcriptomes from two contrasting inbred lines during root growth. Expression patterns of key cambium regulators and hormone signaling components were validated. Clustering and gene ontology (GO) enrichment analyses of radish datasets followed by a comparative analysis against the newly established Arabidopsis early cambium data revealed evolutionary conserved stress-response transcription factors that may intimately control the cambium. Indeed, an in vivo network consisting of selected stress-response and cambium regulators indicated ERF-1 as a potential key checkpoint of cambial activities, explaining how cambium-driven growth is altered in response to environmental changes. The findings here provide valuable information about dynamic gene expression changes during cambium-driven root growth and have implications with regard to future engineering schemes, leading to better crop yields.