Predicting disease genes using protein-protein interactions

Strategic pathway analysis for dual management of epilepsy and comorbid depression: a systems biology perspective

Depression is a common psychiatric comorbidity among patients with epilepsy (PWE), affecting more than a third of PWE. Management of depression may improve quality of life of epileptic patients. Unfortunately, available antidepressants worsen epilepsy by reducing the seizure threshold. This situation demands search of new safer target for combined directorate of epilepsy and comorbid depression. A system biology approach may be useful to find novel pathways/markers for the cure of both epilepsy and associated depression via analyzing available genomic and proteomic information. Hence, the system biology approach using curated 64 seed genes involved in temporal lobe epilepsy and mental depression was applied. The interplay of 600 potential proteins was revealed by the Disease Module Detection (DIAMOnD) Algorithm for the treatment of both epilepsy and comorbid depression using these seed genes. The gene enrichment analysis of seed and diamond genes through DAVID suggested 95 pathways. Selected pathways were refined based on their syn or anti role in epilepsy and depression. In conclusion, total 8 pathways and 27 DIAMOnD genes/proteins were finally deduced as potential new targets for modulation of selected pathways to manage epilepsy and comorbid depression.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-024-00208-1.

Graph Neural Network Model for Prediction of Non-Small Cell Lung Cancer Lymph Node Metastasis Using Protein-Protein Interaction Network and 18F-FDG PET/CT Radiomics

The image texture features obtained from 18F-fluorodeoxyglucose positron emission tomography/computed tomography (18F-FDG PET/CT) images of non-small cell lung cancer (NSCLC) have revealed tumor heterogeneity. A combination of genomic data and radiomics may improve the prediction of tumor prognosis. This study aimed to predict NSCLC metastasis using a graph neural network (GNN) obtained by combining a protein-protein interaction (PPI) network based on gene expression data and image texture features. 18F-FDG PET/CT images and RNA sequencing data of 93 patients with NSCLC were acquired from The Cancer Imaging Archive. Image texture features were extracted from 18F-FDG PET/CT images and area under the curve receiver operating characteristic curve (AUC) of each image feature was calculated. Weighted gene co-expression network analysis (WGCNA) was used to construct gene modules, followed by functional enrichment analysis and identification of differentially expressed genes. The PPI of each gene module and genes belonging to metastasis-related processes were converted via a graph attention network. Images and genomic features were concatenated. The GNN model using PPI modules from WGCNA and metastasis-related functions combined with image texture features was evaluated quantitatively. Fifty-five image texture features were extracted from 18F-FDG PET/CT, and radiomic features were selected based on AUC (n = 10). Eighty-six gene modules were clustered by WGCNA. Genes (n = 19) enriched in the metastasis-related pathways were filtered using DEG analysis. The accuracy of the PPI network, derived from WGCNA modules and metastasis-related genes, improved from 0.4795 to 0.5830 (p < 2.75 × 10-12). Integrating PPI of four metastasis-related genes with 18F-FDG PET/CT image features in a GNN model elevated its accuracy over a without image feature model to 0.8545 (95% CI = 0.8401-0.8689, p-value < 0.02). This model demonstrated significant enhancement compared to the model using PPI and 18F-FDG PET/CT derived from WGCNA (p-value < 0.02), underscoring the critical role of metastasis-related genes in prediction model. The enhanced predictive capability of the lymph node metastasis prediction GNN model for NSCLC, achieved through the integration of comprehensive image features with genomic data, demonstrates promise for clinical implementation.

In vivo mapping of protein-protein interactions of schizophrenia risk factors generates an interconnected disease network

Genetic analyses of Schizophrenia (SCZ) patients have identified thousands of risk factors. In silico protein-protein interaction (PPI) network analysis has provided strong evidence that disrupted PPI networks underlie SCZ pathogenesis. In this study, we performed in vivo PPI analysis of several SCZ risk factors in the rodent brain. Using endogenous antibody immunoprecipitations coupled to mass spectrometry (MS) analysis, we constructed a SCZ network comprising 1612 unique PPI with a 5% FDR. Over 90% of the PPI were novel, reflecting the lack of previous PPI MS studies in brain tissue. Our SCZ PPI network was enriched with known SCZ risk factors, which supports the hypothesis that an accumulation of disturbances in selected PPI networks underlies SCZ. We used Stable Isotope Labeling in Mammals (SILAM) to quantitate phencyclidine (PCP) perturbations in the SCZ network and found that PCP weakened most PPI but also led to some enhanced or new PPI. These findings demonstrate that quantitating PPI in perturbed biological states can reveal alterations to network biology.

Evaluation of input data modality choices on functional gene embeddings

Functional gene embeddings, numerical vectors capturing gene function, provide a promising way to integrate functional gene information into machine learning models. These embeddings are learnt by applying self-supervised machine-learning algorithms on various data types including quantitative omics measurements, protein-protein interaction networks and literature. However, downstream evaluations comparing alternative data modalities used to construct functional gene embeddings have been lacking. Here we benchmarked functional gene embeddings obtained from various data modalities for predicting disease-gene lists, cancer drivers, phenotype-gene associations and scores from genome-wide association studies. Off-the-shelf predictors trained on precomputed embeddings matched or outperformed dedicated state-of-the-art predictors, demonstrating their high utility. Embeddings based on literature and protein-protein interactions inferred from low-throughput experiments outperformed embeddings derived from genome-wide experimental data (transcriptomics, deletion screens and protein sequence) when predicting curated gene lists. In contrast, they did not perform better when predicting genome-wide association signals and were biased towards highly-studied genes. These results indicate that embeddings derived from literature and low-throughput experiments appear favourable in many existing benchmarks because they are biased towards well-studied genes and should therefore be considered with caution. Altogether, our study and precomputed embeddings will facilitate the development of machine-learning models in genetics and related fields.

Noncoding RNAs improve the predictive power of network medicine

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

Exploring pathway interactions to detect molecular mechanisms of disease: 22q11.2 deletion syndrome

BACKGROUND

22q11.2 Deletion Syndrome (22q11DS) is a genetic disorder characterized by the deletion of adjacent genes at a location specified as q11.2 of chromosome 22, resulting in an array of clinical phenotypes including autistic spectrum disorder, schizophrenia, congenital heart defects, and immune deficiency. Many characteristics of the disorder are known, such as the phenotypic variability of the disease and the biological processes associated with it; however, the exact and systemic molecular mechanisms between the deleted area and its resulting clinical phenotypic expression, for example that of neuropsychiatric diseases, are not yet fully understood.

RESULTS

Using previously published transcriptomics data (GEO:GSE59216), we constructed two datasets: one set compares 22q11DS patients experiencing neuropsychiatric diseases versus healthy controls, and the other set 22q11DS patients without neuropsychiatric diseases versus healthy controls. We modified and applied the pathway interaction method, originally proposed by Kelder et al. (2011), on a network created using the WikiPathways pathway repository and the STRING protein-protein interaction database. We identified genes and biological processes that were exclusively associated with the development of neuropsychiatric diseases among the 22q11DS patients. Compared with the 22q11DS patients without neuropsychiatric diseases, patients experiencing neuropsychiatric diseases showed significant overrepresentation of regulated genes involving the natural killer cell function and the PI3K/Akt signalling pathway, with affected genes being closely associated with downregulation of CRK like proto-oncogene adaptor protein. Both the pathway interaction and the pathway overrepresentation analysis observed the disruption of the same biological processes, even though the exact lists of genes collected by the two methods were different.

CONCLUSIONS

Using the pathway interaction method, we were able to detect a molecular network that could possibly explain the development of neuropsychiatric diseases among the 22q11DS patients. This way, our method was able to complement the pathway overrepresentation analysis, by filling the knowledge gaps on how the affected pathways are linked to the original deletion on chromosome 22. We expect our pathway interaction method could be used for problems with similar contexts, where complex genetic mechanisms need to be identified to explain the resulting phenotypic plasticity.

Identification of disease modules using higher-order network structure

Motivation

Higher-order interaction patterns among proteins have the potential to reveal mechanisms behind molecular processes and diseases. While clustering methods are used to identify functional groups within molecular interaction networks, these methods largely focus on edge density and do not explicitly take into consideration higher-order interactions. Disease genes in these networks have been shown to exhibit rich higher-order structure in their vicinity, and considering these higher-order interaction patterns in network clustering have the potential to reveal new disease-associated modules.

Results

We propose a higher-order community detection method which identifies community structure in networks with respect to specific higher-order connectivity patterns beyond edges. Higher-order community detection on four different protein-protein interaction networks identifies biologically significant modules and disease modules that conventional edge-based clustering methods fail to discover. Higher-order clusters also identify disease modules from genome-wide association study data, including new modules that were not discovered by top-performing approaches in a Disease Module DREAM Challenge. Our approach provides a more comprehensive view of community structure that enables us to predict new disease-gene associations.

Availability and implementation

https://github.com/Reed-CompBio/graphlet-clustering.

Combination of antidepressants and antipsychotics as a novel treatment option for psychosis in Alzheimer's disease

Psychotic symptoms are reported as one of the most common complications of Alzheimer's disease (AD), in whom they are associated with more rapid deterioration and increased mortality. Empiric treatments, namely first and second-generation antipsychotics, confer modest efficacy in patients with AD and with psychosis (AD+P) and themselves increase mortality. Recent studies suggested the use and beneficial effects of antidepressants among patients with AD+P. This motivates our rationale for exploring their potential as a novel combination therapy option among these patients. We included electronic medical records of 10,260 patients with AD in our study. Survival analysis was performed to assess the effects of the combination of antipsychotics and antidepressants on the mortality of these patients. A protein-protein interaction network representing AD+P was built, and network analysis methods were used to quantify the efficacy of these drugs on AD+P. A combined score was developed to measure the potential synergetic effect against AD+P. Our survival analyses showed that the co-administration of antidepressants with antipsychotics have a significant beneficial effect in reducing mortality. Our network analysis showed that the targets of antipsychotics and antidepressants are well-separated, and antipsychotics and antidepressants have similar Signed Jaccard Index (SJI) scores to AD+P. Eight drug pairs, including some popular recommendations like aripiprazole/sertraline, showed higher than average scores which suggest their potential in treating AD+P via strong synergetic effects. Our proposed combinations of antipsychotic and antidepressant therapy showed a strong superiority over current antipsychotics treatment for AD+P. The observed beneficial effects can be further strengthened by optimizing drug-pair selection based on our systems pharmacology analysis.

The Role of Apoptotic Genes and Protein-Protein Interactions in Triple-negative Breast Cancer

BACKGROUND/AIM

Compared to other breast cancer types, triple-negative breast cancer (TNBC) has historically had few treatment alternatives. Therefore, exploring and pinpointing potentially implicated genes could be used for treating and managing TNBC. By doing this, we will provide essential data to comprehend how the genes are involved in the apoptotic pathways of the cancer cells to identify potential therapeutic targets. Analysis of a single genetic alteration may not reveal the pathogenicity driving TNBC due to the high genomic complexity and heterogeneity of TNBC. Therefore, searching through a large variety of gene interactions enabled the identification of molecular therapeutic genes.

MATERIALS AND METHODS

This study used integrated bioinformatics methods such as UALCAN, TNM plotter, PANTHER, GO-KEEG and PPIs to assess the gene expression, protein-protein interaction (PPI), and transcription factor interaction of apoptosis-regulated genes.

RESULTS

Compared to normal breast tissue, gene expressions of BNIP3, TNFRSF10B, MCL1, and CASP4 were downregulated in UALCAN. At the same time, BIK, AKT1, BAD, FADD, DIABLO, and CASP9 was down-regulated in bc-GeneExMiner v4.5 mRNA expression (BCGM) databases. Based on GO term enrichment analysis, the cellular process (GO:0009987), which has about 21 apoptosis-regulated genes, is the top category in the biological processes (BP), followed by biological regulation (GO:0065007). We identified 29 differentially regulated pathways, including the p53 pathway, angiogenesis, apoptosis signaling pathway, and the Alzheimer's disease presenilin pathway. We examined the PPIs between the genes that regulate apoptosis; CASP3 and CASP9 interact with FADD, MCL1, TNF, TNFRSRF10A, and TNFRSF10; additionally, CASP3 significantly forms PPIs with CASP9, DFFA, and TP53, and CASP9 with DIABLO. In the top 10 transcription factors, the androgen receptor (AR) interacts with five apoptosis-regulated genes (p<0.0001; q<0.01), followed by retinoic acid receptor alpha (RARA) (p<0.0001; q<0.01) and ring finger protein (RNF2) (p<0.0001; q<0.01). Overall, the gene expression profile, PPIs, and the apoptosis-TF interaction findings suggest that the 27 apoptosis-regulated genes might be used as promising targets in treating and managing TNBC. Furthermore, from a total of 27 key genes, CASP2, CASP3, DAPK1, TNF, TRAF2, and TRAF3 were significantly correlated with poor overall survival in TNBC (p-value <0.05); they could play important roles in the progression of TNBC and provide attractive therapeutic targets that may offer new candidate molecules for targeted therapy.

CONCLUSION

Our findings demonstrate that CASP2, CASP3, DAPK1, TNF, TRAF2, and TRAF3 were substantially associated with the overall survival rate (OS) difference of TNBC patients out of a total of 27 specific genes used in this study, which may play crucial roles in the development of TNBC and offer promising therapeutic interventions.

DGH-GO: dissecting the genetic heterogeneity of complex diseases using gene ontology

BACKGROUND

Complex diseases such as neurodevelopmental disorders (NDDs) exhibit multiple etiologies. The multi-etiological nature of complex-diseases emerges from distinct but functionally similar group of genes. Different diseases sharing genes of such groups show related clinical outcomes that further restrict our understanding of disease mechanisms, thus, limiting the applications of personalized medicine approaches to complex genetic disorders.

RESULTS

Here, we present an interactive and user-friendly application, called DGH-GO. DGH-GO allows biologists to dissect the genetic heterogeneity of complex diseases by stratifying the putative disease-causing genes into clusters that may contribute to distinct disease outcome development. It can also be used to study the shared etiology of complex-diseases. DGH-GO creates a semantic similarity matrix for the input genes by using Gene Ontology (GO). The resultant matrix can be visualized in 2D plots using different dimension reduction methods (T-SNE, Principal component analysis, umap and Principal coordinate analysis). In the next step, clusters of functionally similar genes are identified from genes functional similarities assessed through GO. This is achieved by employing four different clustering methods (K-means, Hierarchical, Fuzzy and PAM). The user may change the clustering parameters and explore their effect on stratification immediately. DGH-GO was applied to genes disrupted by rare genetic variants in Autism Spectrum Disorder (ASD) patients. The analysis confirmed the multi-etiological nature of ASD by identifying four clusters of genes that were enriched for distinct biological mechanisms and clinical outcome. In the second case study, the analysis of genes shared by different NDDs showed that genes causing multiple disorders tend to aggregate in similar clusters, indicating a possible shared etiology.

CONCLUSION

DGH-GO is a user-friendly application that allows biologists to study the multi-etiological nature of complex diseases by dissecting their genetic heterogeneity. In summary, functional similarities, dimension reduction and clustering methods, coupled with interactive visualization and control over analysis allows biologists to explore and analyze their datasets without requiring expert knowledge on these methods. The source code of proposed application is available at https://github.com/Muh-Asif/DGH-GO.

Network medicine: an approach to complex kidney disease phenotypes

Scientific reductionism has been the basis of disease classification and understanding for more than a century. However, the reductionist approach of characterizing diseases from a limited set of clinical observations and laboratory evaluations has proven insufficient in the face of an exponential growth in data generated from transcriptomics, proteomics, metabolomics and deep phenotyping. A new systematic method is necessary to organize these datasets and build new definitions of what constitutes a disease that incorporates both biological and environmental factors to more precisely describe the ever-growing complexity of phenotypes and their underlying molecular determinants. Network medicine provides such a conceptual framework to bridge these vast quantities of data while providing an individualized understanding of disease. The modern application of network medicine principles is yielding new insights into the pathobiology of chronic kidney diseases and renovascular disorders by expanding the understanding of pathogenic mediators, novel biomarkers and new options for renal therapeutics. These efforts affirm network medicine as a robust paradigm for elucidating new advances in the diagnosis and treatment of kidney disorders.

A new view on functions of the lysine demalonylase activity of SIRT5

AIMS

The specificity of the lysine demalonylation substrates of the pharmaceutically attractive tumor promoter/suppressor SIRT5 is not comprehensively clarified. The present study re-analyses publicly available data and highlights potentially pharmaceutically interesting outcomes by the use of bioinformatics.

MATERIALS AND METHODS

The interaction networks of SIRT5 malonylome from the wild-type and ob/ob (obese pre-diabetic type) mice were subjected to the pathway enrichment and gene function prediction analysis using GeneMania (3.5.2) application run under Cytoscape (3.9.1) environment.

KEY FINDINGS

The analysis in the wild-type mice revealed the involvement of SIRT5 malonylome in Eukaryotic translation elongation (ETE; the nodes EF1A1, EEF2, EEF1D, and EEF1G), Amino acid and derivative metabolism (AADM), and Selenoamino acid metabolism (SAM). The tumor promoter/suppressor activity of SIRT5 is mediated through the tumor promoter substrates included in AADM (GLUD1, SHMT1, ACAT1), and the tumor suppressor substrates involved in AADM and SAM (ALDH9A1, BHMT, GNMT). Selen stimulates the expression of SIRT5 and other sirtuins. SIRT5 in turn regulates the selenocysteine synthesis, which creates a regulatory loop. The analysis of SIRT5 malonylome in pre-diabetic ob/ob mice identifies the mTORC1 pathway as a mechanism, which facilitates SIRT5 functions. The comparison of the outcomes of SIRT5 malonylome, succinylome, and glutarylome analysis disclosed several differences.

SIGNIFICANCE

The analysis showed additional aspects of SIRT5 malonylome functions besides the control of glucose metabolism. It defined several unique substrates and pathways, and it showed differences compared to other enzymatic activities of SIRT5, which could be used for pharmaceutical benefits.

Network expansion of genetic associations defines a pleiotropy map of human cell biology

Interacting proteins tend to have similar functions, influencing the same organismal traits. Interaction networks can be used to expand the list of candidate trait-associated genes from genome-wide association studies. Here, we performed network-based expansion of trait-associated genes for 1,002 human traits showing that this recovers known disease genes or drug targets. The similarity of network expansion scores identifies groups of traits likely to share an underlying genetic and biological process. We identified 73 pleiotropic gene modules linked to multiple traits, enriched in genes involved in processes such as protein ubiquitination and RNA processing. In contrast to gene deletion studies, pleiotropy as defined here captures specifically multicellular-related processes. We show examples of modules linked to human diseases enriched in genes with known pathogenic variants that can be used to map targets of approved drugs for repurposing. Finally, we illustrate the use of network expansion scores to study genes at inflammatory bowel disease genome-wide association study loci, and implicate inflammatory bowel disease-relevant genes with strong functional and genetic support.

Neurodegenerative disorders: From clinicopathology convergence to systems biology divergence

Neurodegenerative diseases are multifactorial. This means that several genetic, epigenetic, and environmental factors contribute to their emergence. Therefore, for the future management of these highly prevalent diseases, it is necessary to change perspective. If a holistic viewpoint is assumed, the phenotype (the clinicopathological convergence) emerges from the perturbation of a complex system of functional interactions among proteins (systems biology divergence). The systems biology top-down approach starts with the unbiased collection of sets of data generated through one or more -omics techniques and has the aim to identify the networks and the components that participate in the generation of a phenotype (disease), often without any available a priori knowledge. The principle behind the top-down method is that the molecular components that respond similarly to experimental perturbations are somehow functionally related. This allows the study of complex and relatively poorly characterized diseases without requiring extensive knowledge of the processes under investigation. In this chapter, the use of a global approach will be applied to the comprehension of neurodegeneration, with a particular focus on the two most prevalent ones, Alzheimer's and Parkinson's diseases. The final purpose is to distinguish disease subtypes (even with similar clinical manifestations) to launch a future of precision medicine for patients with these disorders.

PGAGP: Predicting pathogenic genes based on adaptive network embedding algorithm

The study of disease-gene associations is an important topic in the field of computational biology. The accumulation of massive amounts of biomedical data provides new possibilities for exploring potential relations between diseases and genes through computational strategy, but how to extract valuable information from the data to predict pathogenic genes accurately and rapidly is currently a challenging and meaningful task. Therefore, we present a novel computational method called PGAGP for inferring potential pathogenic genes based on an adaptive network embedding algorithm. The PGAGP algorithm is to first extract initial features of nodes from a heterogeneous network of diseases and genes efficiently and effectively by Gaussian random projection and then optimize the features of nodes by an adaptive refining process. These low-dimensional features are used to improve the disease-gene heterogenous network, and we apply network propagation to the improved heterogenous network to predict pathogenic genes more effectively. By a series of experiments, we study the effect of PGAGP's parameters and integrated strategies on predictive performance and confirm that PGAGP is better than the state-of-the-art algorithms. Case studies show that many of the predicted candidate genes for specific diseases have been implied to be related to these diseases by literature verification and enrichment analysis, which further verifies the effectiveness of PGAGP. Overall, this work provides a useful solution for mining disease-gene heterogeneous network to predict pathogenic genes more effectively.

Network analysis of genome-wide association studies for drug target prioritisation

Over the past decades, genome-wide association studies (GWAS) have led to a dramatic expansion of genetic variants implicated with human traits and diseases. These advances are expected to result in new drug targets but the identification of causal genes and the cell biology underlying human diseases from GWAS remains challenging. Here, we review protein interaction network-based methods to analyse GWAS data. These approaches can rank candidate drug targets at GWAS-associated loci or among interactors of disease genes without direct genetic support. These methods identify the cell biology affected in common across diseases, offering opportunities for drug repurposing, as well as be combined with expression data to identify focal tissues and cell types. Going forward, we expect that these methods will further improve from advances in the characterisation of context specific interaction networks and the joint analysis of rare and common genetic signals.

Efficacy difference of antipsychotics in Alzheimer's disease and schizophrenia: explained with network efficiency and pathway analysis methods

Approximately 50% of Alzheimer's disease (AD) patients will develop psychotic symptoms and these patients will experience severe rapid cognitive decline compared with those without psychosis (AD-P). Currently, no medication has been approved by the Food and Drug Administration for AD with psychosis (AD+P) specifically, although atypical antipsychotics are widely used in clinical practice. These drugs have demonstrated modest efficacy in managing psychosis in individuals with AD, with an increased frequency of adverse events, including excess mortality. We compared the differences between the genetic variations/genes associated with AD+P and schizophrenia from existing Genome-Wide Association Study and differentially expressed genes (DEGs). We also constructed disease-specific protein-protein interaction networks for AD+P and schizophrenia. Network efficiency was then calculated to characterize the topological structures of these two networks. The efficiency of antipsychotics in these two networks was calculated. A weight adjustment based on binding affinity to drug targets was later applied to refine our results, and 2013 and 2123 genes were identified as related to AD+P and schizophrenia, respectively, with only 115 genes shared. Antipsychotics showed a significantly lower efficiency in the AD+P network than in the schizophrenia network (P < 0.001) indicating that antipsychotics may have less impact in AD+P than in schizophrenia. AD+P may be caused by mechanisms distinct from those in schizophrenia which result in a decreased efficacy of antipsychotics in AD+P. In addition, the network analysis methods provided quantitative explanations of the lower efficacy of antipsychotics in AD+P.

Inferring Latent MicroRNA-Disease Associations on a Gene-Mediated Tripartite Heterogeneous Multiplexing Network

MicroRNA (miRNA) is a class of non-coding single-stranded RNA molecules encoded by endogenous genes with a length of about 22 nucleotides. MiRNAs have been successfully identified as differentially expressed in various cancers. There is evidence that disorders of miRNAs are associated with a variety of complex diseases. Therefore, inferring potential miRNA-disease associations (MDAs) is very important for understanding the aetiology and pathogenesis of many diseases and is useful to disease diagnosis, prognosis and treatment. First, We creatively fused multiple similarity subnetworks from multi-sources for miRNAs, genes and diseases by multiplexing technology, respectively. Then, three multiplexed biological subnetworks are connected through the extended binary association to form a tripartite complete heterogeneous multiplexed network (Tri-HM). Finally, because the constructed Tri-HM network can retain subnetworks' original topology and biological functions and expands the binary association and dependence between the three biological entities, rich neighbourhood information is obtained iteratively from neighbours by a non-equilibrium random walk. Through cross-validation, our tri-HM-RWR model obtained an AUC value of 0.8657, and an AUPR value of 0.2139 in the global 5-fold cross-validation, which shows that our model can more fully speculate disease-related miRNAs.

Computational systems biology in disease modeling and control, review and perspectives

Omics-based approaches have become increasingly influential in identifying disease mechanisms and drug responses. Considering that diseases and drug responses are co-expressed and regulated in the relevant omics data interactions, the traditional way of grabbing omics data from single isolated layers cannot always obtain valuable inference. Also, drugs have adverse effects that may impair patients, and launching new medicines for diseases is costly. To resolve the above difficulties, systems biology is applied to predict potential molecular interactions by integrating omics data from genomic, proteomic, transcriptional, and metabolic layers. Combined with known drug reactions, the resulting models improve medicines' therapeutical performance by re-purposing the existing drugs and combining drug molecules without off-target effects. Based on the identified computational models, drug administration control laws are designed to balance toxicity and efficacy. This review introduces biomedical applications and analyses of interactions among gene, protein and drug molecules for modeling disease mechanisms and drug responses. The therapeutical performance can be improved by combining the predictive and computational models with drug administration designed by control laws. The challenges are also discussed for its clinical uses in this work.

Identification of differentially expressed genes and their major pathways among the patient with COVID-19, cystic fibrosis, and chronic kidney disease

The SARS-CoV-2 virus causes Coronavirus disease, an infectious disease. The majority of people who are infected with this virus will have mild to moderate respiratory symptoms. Multiple studies have proved that there is a substantial pathophysiological link between COVID-19 disease and patients having comorbidities such as cystic fibrosis and chronic kidney disease. In this study, we attempted to identify differentially expressed genes as well as genes that intersected among them in order to comprehend their compatibility. Gene expression profiling indicated that 849 genes were mutually exclusive and functional analysis was done within the context of gene ontology and key pathways involvement. Three genes (PRPF31, FOXN2, and RIOK3) were commonly upregulated in the analysed datasets of three disease categories. These genes could be potential biomarkers for patients with COVID-19 and cystic fibrosis, and COVID-19 and chronic kidney disease. Further extensive analyses have been performed to describe how these genes are regulated by various transcription factors and microRNAs. Then, our analyses revealed six hub genes (PRPF31, FOXN2, RIOK3, UBC, HNF4A, and ELAVL). As they were involved in the interaction between COVID-19 and the patient with CF and CKD, they could help researchers identify potential therapeutic molecules. Some drugs have been predicted based on the upregulated genes, which may have a significant impact on reducing the burden of these diseases in the future.

Diabetes mellitus susceptibility with varied diseased phenotypes and its comparison with phenome interactome networks

An emerging area of interest in understanding disease phenotypes is systems genomics. Complex diseases such as diabetes have played an important role towards understanding the susceptible genes and mutations. A wide number of methods have been employed and strategies such as polygenic risk score and allele frequencies have been useful, but understanding the candidate genes harboring those mutations is an unmet goal. In this perspective, using systems genomic approaches, we highlight the application of phenome-interactome networks in diabetes and provide deep insights. LINC01128, which we previously described as candidate for diabetes, is shown as an example to discuss the approach.

Network-Based Methods for Approaching Human Pathologies from a Phenotypic Point of View

Network and systemic approaches to studying human pathologies are helping us to gain insight into the molecular mechanisms of and potential therapeutic interventions for human diseases, especially for complex diseases where large numbers of genes are involved. The complex human pathological landscape is traditionally partitioned into discrete “diseases”; however, that partition is sometimes problematic, as diseases are highly heterogeneous and can differ greatly from one patient to another. Moreover, for many pathological states, the set of symptoms (phenotypes) manifested by the patient is not enough to diagnose a particular disease. On the contrary, phenotypes, by definition, are directly observable and can be closer to the molecular basis of the pathology. These clinical phenotypes are also important for personalised medicine, as they can help stratify patients and design personalised interventions. For these reasons, network and systemic approaches to pathologies are gradually incorporating phenotypic information. This review covers the current landscape of phenotype-centred network approaches to study different aspects of human diseases.

Discovering miRNAs Associated With Multiple Sclerosis Based on Network Representation Learning and Deep Learning Methods

Identifying biomarkers of Multiple Sclerosis is important for the diagnosis and treatment of Multiple Sclerosis. The existing study has shown that miRNA is one of the most important biomarkers for diseases. However, few existing methods are designed for predicting Multiple Sclerosis-related miRNAs. To fill this gap, we proposed a novel computation framework for predicting Multiple Sclerosis-associated miRNAs. The proposed framework uses a network representation model to learn the feature representation of miRNA and uses a deep learning-based model to predict the miRNAs associated with Multiple Sclerosis. The evaluation result shows that the proposed model can predict the miRNAs associated with Multiple Sclerosis precisely. In addition, the proposed model can outperform several existing methods in a large margin.

Cordycepin: a bioactive metabolite of Cordyceps militaris and polyadenylation inhibitor with therapeutic potential against COVID-19

Spike protein and main proteases of SARS-CoV-2 have been identified as potential therapeutic targets and their inhibition may lead to the reticence of viral entry and replication in the host body. Despite several efforts; till now no specific drugs are available to treat SARS-CoV-2. Considering all these challenges, the main objective of the present study was to establish therapeutic potential of cordycepin against COVID-19 as a conventional therapeutic strategy. In the present study; molecular interaction study was performed to assess potential binding affinity of cordycepin with SARS-CoV-2 target proteins using computational approach. Additionally, network pharmacology was used to understand cordycepin-protein interactions and their associated pathways in human body. Cordycepin is under clinical trial (NCT00709215) and possesses structural similarity with adenosine except that, it lacks a 3' hydroxyl group in its ribose moiety and hence it served as a poly(A) polymerase inhibitor and terminate premature protein synthesis. Additionally, it is known that functional RNAs of SARS-CoV-2 genome are highly 3'-plyadenylated and leading to synthesis of all viral proteins and if cordycepin can destabilize SARS-CoV-2 RNAs by inhibiting polyadenylation process then it may step forward in terms of inhibition of viral replication and multiplication in the host. Moreover, cordycepin showed strong binding affinity with SARS-CoV-2 spike protein (-145.3) and main proteases (-180.5) that further corroborate therapeutic potential against COVID-19. Since cordycepin has both pre-clinical and clinical information about antiviral activities, therefore; it is suggested to the world community to undertake repurposing cordycepin to test efficacy and safety for the treatment of COVID-19.

Protein interaction networks provide insight into fetal origins of chronic obstructive pulmonary disease

BACKGROUND

Chronic obstructive pulmonary disease (COPD) is a leading cause of death in adults that may have origins in early lung development. It is a complex disease, influenced by multiple factors including genetic variants and environmental factors. Maternal smoking during pregnancy may influence the risk for diseases during adulthood, potentially through epigenetic modifications including methylation.

METHODS

In this work, we explore the fetal origins of COPD by utilizing lung DNA methylation marks associated with in utero smoke (IUS) exposure, and evaluate the network relationships between methylomic and transcriptomic signatures associated with adult lung tissue from former smokers with and without COPD. To identify potential pathobiological mechanisms that may link fetal lung, smoke exposure and adult lung disease, we study the interactions (physical and functional) of identified genes using protein-protein interaction networks.

RESULTS

We build IUS-exposure and COPD modules, which identify connected subnetworks linking fetal lung smoke exposure to adult COPD. Studying the relationships and connectivity among the different modules for fetal smoke exposure and adult COPD, we identify enriched pathways, including the AGE-RAGE and focal adhesion pathways.

CONCLUSIONS

The modules identified in our analysis add new and potentially important insights to understanding the early life molecular perturbations related to the pathogenesis of COPD. We identify AGE-RAGE and focal adhesion as two biologically plausible pathways that may reveal lung developmental contributions to COPD. We were not only able to identify meaningful modules but were also able to study interconnections between smoke exposure and lung disease, augmenting our knowledge about the fetal origins of COPD.

Finding New Ways How to Control BACE1

Recently, all applications of BACE1 inhibitors failed as therapeutical targets for Alzheimer´s disease (AD) due to severe side effects. Therefore, alternative ways for treatment development are a hot research topic. The present analysis investigates BACE1 protein-protein interaction networks and attempts to solve the absence of complete knowledge about pathways involving BACE1. A bioinformatics analysis matched the functions of the non-substrate interaction network with Voltage-gated potassium channels, which also appear as top priority protein nodes. Targeting BACE1 interactions with PS1 and GGA-s, blocking of BACE1 access to APP by BRI3 and RTN-s, activation of Wnt signaling and upregulation of β-catenin, and brain delivery of the extracellular domain of p75^NTR, are the main alternatives to the use of BACE 1 inhibitors highlighted by the analysis. The pathway enrichment analysis also emphasized substrates and substrate candidates with essential biological functions, which cleavage must remain controlled. They include ephrin receptors, ROBO1, ROBO2, CNTN-s, CASPR-s, CD147, CypB, TTR, APLP1/APLP2, NRXN-s, and PTPR-s. The analysis of the interaction subnetwork of BACE1 functionally related to inflammation identified a connection to three cardiomyopathies, which supports the hypothesis of the common molecular mechanisms with AD. A lot of potential shows the regulation of BACE1 activity through post-translational modifications. The interaction network of BACE1 and its phosphorylation enzyme CSNK1D functionally match the Circadian clock, p53, and Hedgehog signaling pathways. The regulation of BACE1 glycosylation could be achieved through N-acetylglucosamine transferases, α-(1→6)-fucosyltransferase, β-galactoside α-(2→6)-sialyltransferases, galactosyltransferases, and mannosidases suggested by the interaction network analysis of BACE1-MGAT3. The present analysis proposes possibilities for the alternative control of AD pathology.

Scalable Non-Linear Graph Fusion for Prioritizing Cancer-Causing Genes

In the past few decades, both gene expression data and protein-protein interaction (PPI)networks have been extensively studied, due to their ability to depict important characteristics of disease-associated genes. In this regard, the paper presents a new gene prioritization algorithm to identify and prioritize cancer-causing genes, integrating judiciously the complementary information obtained from two data sources. The proposed algorithm selects disease-causing genes by maximizing the importance of selected genes and functional similarity among them. A new quantitative index is introduced to evaluate the importance of a gene. It considers whether a gene exhibits a differential expression pattern across sick and healthy individuals, and has a strong connectivity in the PPI network, which are the important characteristics of a potential biomarker. As disease-associated genes are expected to have similar expression profiles and topological structures, a scalable non-linear graph fusion technique, termed as ScaNGraF, is proposed to learn a disease-dependent functional similarity network from the co-expression and common neighbor based similarity networks. The proposed ScaNGraF, which is based on message passing algorithm, efficiently combines the shared and complementary information provided by different data sources with significantly lower computational cost. A new measure, termed as DiCoIN, is introduced to evaluate the quality of a learned affinity network. The performance of the proposed graph fusion technique and gene selection algorithm is extensively compared with that of some existing methods, using several cancer data sets.

Network-Based Approach to Repurpose Approved Drugs for COVID-19 by Integrating GWAS and Text Mining Data

The coronavirus disease 19 (COVID-19) is a global pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which has a rapidly increasing prevalence and has caused significant morbidity/mortality. Despite the availability of many vaccines that can offer widespread immunization, it is also important to reach effective treatment for COVID-19 patients. However, the development of novel drug therapeutics is usually a time-consuming and costly process, and therefore, repositioning drugs that were previously approved for other purposes could have a major impact on the fight against COVID-19. Here, we first identified lung-specific gene regulatory/interaction subnetworks (COVID-19-related genes modules) enriched for COVID-19-associated genes obtained from GWAS and text mining. We then screened the targets of 220 approved drugs from DrugBank, obtained their drug-induced gene expression profiles in the LINCS database, and constructed lung-specific drug-related gene modules. By applying an integrated network-based approach to quantify the interactions of the COVID-19-related gene modules and drug-related gene modules, we prioritized 13 approved drugs (e.g., alitretinoin, clocortolone, terazosin, doconexent, and pergolide) that could potentially be repurposed for the treatment of COVID-19. These findings provide important and timely insights into alternative therapeutic options that should be further explored as COVID-19 continues to spread.

HPODNets: deep graph convolutional networks for predicting human protein-phenotype associations

MOTIVATION

Deciphering the relationship between human genes/proteins and abnormal phenotypes is of great importance in the prevention, diagnosis and treatment against diseases. The Human Phenotype Ontology (HPO) is a standardized vocabulary that describes the phenotype abnormalities encountered in human disorders. However, the current HPO annotations are still incomplete. Thus, it is necessary to computationally predict human protein-phenotype associations. In terms of current, cutting-edge computational methods for annotating proteins (such as functional annotation), three important features are (i) multiple network input, (ii) semi-supervised learning and (iii) deep graph convolutional network (GCN), whereas there are no methods with all these features for predicting HPO annotations of human protein.

RESULTS

We develop HPODNets with all above three features for predicting human protein-phenotype associations. HPODNets adopts a deep GCN with eight layers which allows to capture high-order topological information from multiple interaction networks. Empirical results with both cross-validation and temporal validation demonstrate that HPODNets outperforms seven competing state-of-the-art methods for protein function prediction. HPODNets with the architecture of deep GCNs is confirmed to be effective for predicting HPO annotations of human protein and, more generally, node label ranking problem with multiple biomolecular networks input in bioinformatics.

AVAILABILITY AND IMPLEMENTATION

https://github.com/liulizhi1996/HPODNets.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Focus on Molecular Functions of Anti-Aging Deacetylase SIRT3

SIRT3 is a protein lysine deacetylase with a prominent role in the maintenance of mitochondrial integrity, which is a vulnerable target in many diseases. Intriguingly, cellular aging is reversible just by SIRT3 overexpression, which raises many questions about the role of SIRT3 in the molecular anti-aging mechanisms. Therefore, functions of SIRT3 were analyzed through the interaction network of 407 substrates collected by data mining. Results of the pathway enrichment and gene function prediction confirmed functions in the primary metabolism and mitochondrial ATP production. However, it also suggested involvement in thermogenesis, brain-related neurodegenerative diseases Alzheimer's (AD), Parkinson's, Huntington's disease (HD), and non-alcoholic fatty liver disease. The protein node prioritization analysis identified subunits of the complex I of the mitochondrial respiratory chain (MRC) as the nodes with the main regulatory effect within the entire interaction network. Additional high-ranked nodes were succinate dehydrogenase subunit B (SDHB), complex II, and ATP5F1, complex V of MRC. The analysis supports existence of the NADH/NAD+ driven regulatory feedback loop between SIRT3, complex I (MRC), and acetyl-CoA synthetases, and existence of the nuclear substrates of SIRT3. Unexplored functions of SIRT3 substrates such as LMNA and LMNB; HIF-1a, p53, DNA-PK, and PARK7 are highlighted for further scientific advances. SIRT3 acts as a repressor of BACE1 through the SIRT3-LKB1-AMPK-CREB-PGC1A-PPARG-BACE1 (SIRT3-BACE1), which functions are fitted the best by the Circadian Clock pathway. It forms a new working hypothesis as the therapeutical target for AD treatment. Other important pathways linked to SIRT3 activity are highlighted for therapeutical interventions.

GPCRs Are Optimal Regulators of Complex Biological Systems and Orchestrate the Interface between Health and Disease

GPCRs arguably represent the most effective current therapeutic targets for a plethora of diseases. GPCRs also possess a pivotal role in the regulation of the physiological balance between healthy and pathological conditions; thus, their importance in systems biology cannot be underestimated. The molecular diversity of GPCR signaling systems is likely to be closely associated with disease-associated changes in organismal tissue complexity and compartmentalization, thus enabling a nuanced GPCR-based capacity to interdict multiple disease pathomechanisms at a systemic level. GPCRs have been long considered as controllers of communication between tissues and cells. This communication involves the ligand-mediated control of cell surface receptors that then direct their stimuli to impact cell physiology. Given the tremendous success of GPCRs as therapeutic targets, considerable focus has been placed on the ability of these therapeutics to modulate diseases by acting at cell surface receptors. In the past decade, however, attention has focused upon how stable multiprotein GPCR superstructures, termed receptorsomes, both at the cell surface membrane and in the intracellular domain dictate and condition long-term GPCR activities associated with the regulation of protein expression patterns, cellular stress responses and DNA integrity management. The ability of these receptorsomes (often in the absence of typical cell surface ligands) to control complex cellular activities implicates them as key controllers of the functional balance between health and disease. A greater understanding of this function of GPCRs is likely to significantly augment our ability to further employ these proteins in a multitude of diseases.

Discovering Cerebral Ischemic Stroke Associated Genes Based on Network Representation Learning

Cerebral ischemic stroke (IS) is a complex disease caused by multiple factors including vascular risk factors, genetic factors, and environment factors, which accentuates the difficulty in discovering corresponding disease-related genes. Identifying the genes associated with IS is critical for understanding the biological mechanism of IS, which would be significantly beneficial to the diagnosis and clinical treatment of cerebral IS. However, existing methods to predict IS-related genes are mainly based on the hypothesis of guilt-by-association (GBA). These methods cannot capture the global structure information of the whole protein-protein interaction (PPI) network. Inspired by the success of network representation learning (NRL) in the field of network analysis, we apply NRL to the discovery of disease-related genes and launch the framework to identify the disease-related genes of cerebral IS. The utilized framework contains three main parts: capturing the topological information of the PPI network with NRL, denoising the gene feature with the participation of a stacked autoencoder (SAE), and optimizing a support vector machine (SVM) classifier to identify IS-related genes. Superior to the existing methods on IS-related gene prediction, our framework presents more accurate results. The case study also shows that the proposed method can identify IS-related genes.

Unravelling the Shared Genetic Mechanisms Underlying 18 Autoimmune Diseases Using a Systems Approach

Autoimmune diseases (AiDs) are complex heterogeneous diseases characterized by hyperactive immune responses against self. Genome-wide association studies have identified thousands of single nucleotide polymorphisms (SNPs) associated with several AiDs. While these studies have identified a handful of pleiotropic loci that confer risk to multiple AiDs, they lack the power to detect shared genetic factors residing outside of these loci. Here, we integrated chromatin contact, expression quantitative trait loci and protein-protein interaction (PPI) data to identify genes that are regulated by both pleiotropic and non-pleiotropic SNPs. The PPI analysis revealed complex interactions between the shared and disease-specific genes. Furthermore, pathway enrichment analysis demonstrated that the shared genes co-occur with disease-specific genes within the same biological pathways. In conclusion, our results are consistent with the hypothesis that genetic risk loci associated with multiple AiDs converge on a core set of biological processes that potentially contribute to the emergence of polyautoimmunity.

Computational Methods for Prediction of Human Protein-Phenotype Associations: A Review

Deciphering the relationship between human proteins (genes) and phenotypes is one of the fundamental tasks in phenomics research. The Human Phenotype Ontology (HPO) builds upon a standardized logical vocabulary to describe the abnormal phenotypes encountered in human diseases and paves the way towards the computational analysis of their genetic causes. To date, many computational methods have been proposed to predict the HPO annotations of proteins. In this paper, we conduct a comprehensive review of the existing approaches to predicting HPO annotations of novel proteins, identifying missing HPO annotations, and prioritizing candidate proteins with respect to a certain HPO term. For each topic, we first give the formalized description of the problem, and then systematically revisit the published literatures highlighting their advantages and disadvantages, followed by the discussion on the challenges and promising future directions. In addition, we point out several potential topics to be worthy of exploration including the selection of negative HPO annotations and detecting HPO misannotations. We believe that this review will provide insight to the researchers in the field of computational phenotype analyses in terms of comprehending and developing novel prediction algorithms.

Optimized Molecular Interaction Networks for the Study of Skeletal Muscle

BACKGROUND

Molecular interaction networks (MINs) aim to capture the complex relationships between interacting molecules within a biological system. MINs can be constructed from existing knowledge of molecular functional associations, such as protein-protein binding interactions (PPI) or gene co-expression, and these different sources may be combined into a single MIN. A given MIN may be more or less optimal in its representation of the important functional relationships of molecules in a tissue.

OBJECTIVE

The aim of this study was to establish whether a combined MIN derived from different types of functional association could better capture muscle-relevant biology compared to its constituent single-source MINs.

METHODS

MINs were constructed from functional association databases for both protein-binding and gene co-expression. The networks were then compared based on the capture of muscle-relevant genes and gene ontology (GO) terms, tested in two different ways using established biological network clustering algorithms. The top performing MINs were combined to test whether an optimal MIN for skeletal muscle could be constructed.

RESULTS

The STRING PPI network was the best performing single-source MIN among those tested. Combining STRING with interactions from either the MyoMiner or CoXPRESSdb gene co-expression sources resulted in a combined network with improved performance relative to its constituent networks.

CONCLUSION

MINs constructed from multiple types of functional association can better represent the functional relationships of molecules in a given tissue. Such networks may be used to improve the analysis and interpretation of functional genomics data in the study of skeletal muscle and neuromuscular diseases. Networks and clusters described by this study, including the combinations of STRING with MyoMiner or with CoXPRESSdb, are available for download from https://www.sys-myo.com/myominer/download.php.

Detection of Metabolite-Protein Interactions in Complex Biological Samples by High-Resolution Relaxometry: Toward Interactomics by NMR

Metabolomics, the systematic investigation of metabolites in biological fluids, cells, or tissues, reveals essential information about metabolism and diseases. Metabolites have functional roles in a myriad of biological processes, as substrates and products of enzymatic reactions but also as cofactors and regulators of large numbers of biochemical mechanisms. These functions involve interactions of metabolites with macromolecules. Yet, methods to systematically investigate these interactions are still scarce to date. In particular, there is a need for techniques suited to identify and characterize weak metabolite-macromolecule interactions directly in complex media such as biological fluids. Here, we introduce a method to investigate weak interactions between metabolites and macromolecules in biological fluids. Our approach is based on high-resolution NMR relaxometry and does not require any invasive procedure or separation step. We show that we can detect interactions between small and large molecules in human blood serum and quantify the size of the complex. Our work opens the way for investigations of metabolite (or other small molecules)-protein interactions in biological fluids for interactomics or pharmaceutical applications.

HPOFiller: identifying missing protein-phenotype associations by graph convolutional network

MOTIVATION

Exploring the relationship between human proteins and abnormal phenotypes is of great importance in the prevention, diagnosis and treatment of diseases. The human phenotype ontology (HPO) is a standardized vocabulary that describes the phenotype abnormalities encountered in human diseases. However, the current HPO annotations of proteins are not complete. Thus, it is important to identify missing protein-phenotype associations.

RESULTS

We propose HPOFiller, a graph convolutional network (GCN)-based approach, for predicting missing HPO annotations. HPOFiller has two key GCN components for capturing embeddings from complex network structures: 1) S-GCN for both protein-protein interaction (PPI) network and HPO semantic similarity network to utilize network weights; 2) Bi-GCN for the protein-phenotype bipartite graph to conduct message passing between proteins and phenotypes. The core idea of HPOFiller is to repeat run these two GCN modules consecutively over the three networks, to refine the embeddings. Empirical results of extremely stringent evaluation avoiding potential information leakage including cross-validation and temporal validation demonstrates that HPOFiller significantly outperforms all other state-of-the-art methods. In particular, the ablation study shows that batch normalization contributes the most to the performance. The further examination offers literature evidence for highly ranked predictions. Finally using known disease-HPO term associations, HPOFiller could suggest promising, unknown disease-gene associations, presenting possible genetic causes of human disorders.

AVAILABILITY

https://github.com/liulizhi1996/HPOFiller.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

SIGNAL: A web-based iterative analysis platform integrating pathway and network approaches optimizes hit selection from genome-scale assays

Hit selection from high-throughput assays remains a critical bottleneck in realizing the potential of omic-scale studies in biology. Widely used methods such as setting of cutoffs, prioritizing pathway enrichments, or incorporating predicted network interactions offer divergent solutions yet are associated with critical analytical trade-offs. The specific limitations of these individual approaches and the lack of a systematic way by which to integrate their rankings have contributed to limited overlap in the reported results from comparable genome-wide studies and costly inefficiencies in secondary validation efforts. Using comparative analysis of parallel independent studies as a benchmark, we characterize the specific complementary contributions of each approach and demonstrate an optimal framework to integrate these methods. We describe selection by iterative pathway group and network analysis looping (SIGNAL), an integrated, iterative approach that uses both pathway and network methods to optimize gene prioritization. SIGNAL is accessible as a rapid user-friendly web-based application (https://signal.niaid.nih.gov). A record of this paper's transparent peer review is included in the Supplemental information.

Network-based analysis using chromosomal microdeletion syndromes as a model

Microdeletion syndromes (MSs) are a heterogeneous group of genetic diseases that can virtually affect all functions and organs in humans. Although systems biology approaches integrating multiomics and database information into biological networks have expanded our knowledge of genetic disorders, cytogenomic network-based analysis has rarely been applied to study MSs. In this study, we analyzed data of 28 MSs, using network-based approaches, to investigate the associations between the critical chromosome regions and the respective underlying biological network systems. We identified MSs-associated proteins that were organized in a network of linked modules within the human interactome. Certain MSs formed highly interlinked self-contained disease modules. Furthermore, we observed disease modules involving proteins from other disease groups in the MSs interactome. Moreover, analysis of integrated data from 564 genes located in known chromosomal critical regions, including those contributing to topological parameters, shared pathways, and gene-disease associations, indicated that complex biological systems and cellular networks may underlie many genotype to phenotype associations in MSs. In conclusion, we used a network-based analysis to provide resources that may contribute to better understanding of the molecular pathways involved in MSs.

PRYNT: a tool for prioritization of disease candidates from proteomics data using a combination of shortest-path and random walk algorithms

The urinary proteome is a promising pool of biomarkers of kidney disease. However, the protein changes observed in urine only partially reflect the deregulated mechanisms within kidney tissue. In order to improve on the mechanistic insight based on the urinary protein changes, we developed a new prioritization strategy called PRYNT (PRioritization bY protein NeTwork) that employs a combination of two closeness-based algorithms, shortest-path and random walk, and a contextualized protein-protein interaction (PPI) network, mainly based on clique consolidation of STRING network. To assess the performance of our approach, we evaluated both precision and specificity of PRYNT in prioritizing kidney disease candidates. Using four urinary proteome datasets, PRYNT prioritization performed better than other prioritization methods and tools available in the literature. Moreover, PRYNT performed to a similar, but complementary, extent compared to the upstream regulator analysis from the commercial Ingenuity Pathway Analysis software. In conclusion, PRYNT appears to be a valuable freely accessible tool to predict key proteins indirectly from urinary proteome data. In the future, PRYNT approach could be applied to other biofluids, molecular traits and diseases. The source code is freely available on GitHub at: https://github.com/Boizard/PRYNT and has been integrated as an interactive web apps to improved accessibility ( https://github.com/Boizard/PRYNT/tree/master/AppPRYNT ).

A Dual Systems Genetics Approach Identifies Common Genes, Networks, and Pathways for Type 1 and 2 Diabetes in Human Islets

Type 1 and 2 diabetes (T1/2D) are complex metabolic diseases caused by absolute or relative loss of functional β-cell mass, respectively. Both diseases are influenced by multiple genetic loci that alter disease risk. For many of the disease-associated loci, the causal candidate genes remain to be identified. Remarkably, despite the partially shared phenotype of the two diabetes forms, the associated loci for T1D and T2D are almost completely separated. We hypothesized that some of the genes located in risk loci for T1D and T2D interact in common pancreatic islet networks to mutually regulate important islet functions which are disturbed by disease-associated variants leading to β-cell dysfunction. To address this, we took a dual systems genetics approach. All genes located in 57 T1D and 243 T2D established genome-wide association studies (GWAS) loci were extracted and filtered for genes expressed in human islets using RNA sequencing data, and then integrated with; (1) human islet expression quantitative trait locus (eQTL) signals in linkage disequilibrium (LD) with T1D- and T2D-associated variants; or (2) with genes transcriptionally regulated in human islets by pro-inflammatory cytokines or palmitate as in vitro models of T1D and T2D, respectively. Our in silico systems genetics approaches created two interaction networks consisting of densely-connected T1D and T2D loci genes. The "T1D-T2D islet eQTL interaction network" identified 9 genes (GSDMB, CARD9, DNLZ, ERAP1, PPIP5K2, TMEM69, SDCCAG3, PLEKHA1, and HEMK1) in common T1D and T2D loci that harbor islet eQTLs in LD with disease-associated variants. The "cytokine and palmitate islet interaction network" identified 4 genes (ASCC2, HIBADH, RASGRP1, and SRGAP2) in common T1D and T2D loci whose expression is mutually regulated by cytokines and palmitate. Functional annotation analyses of the islet networks revealed a number of significantly enriched pathways and molecular functions including cell cycle regulation, inositol phosphate metabolism, lipid metabolism, and cell death and survival. In summary, our study has identified a number of new plausible common candidate genes and pathways for T1D and T2D.

Identification of KIF4A and its effect on the progression of lung adenocarcinoma based on the bioinformatics analysis

BACKGROUND

Lung adenocarcinoma (LUAD) is the most frequent histological type of lung cancer, and its incidence has displayed an upward trend in recent years. Nevertheless, little is known regarding effective biomarkers for LUAD.

METHODS

The robust rank aggregation method was used to mine differentially expressed genes (DEGs) from the gene expression omnibus (GEO) datasets. The Search Tool for the Retrieval of Interacting Genes (STRING) database was used to extract hub genes from the protein-protein interaction (PPI) network. The expression of the hub genes was validated using expression profiles from TCGA and Oncomine databases and was verified by real-time quantitative PCR (qRT-PCR). The module and survival analyses of the hub genes were determined using Cytoscape and Kaplan-Meier curves. The function of KIF4A as a hub gene was investigated in LUAD cell lines.

RESULTS

The PPI analysis identified seven DEGs including BIRC5, DLGAP5, CENPF, KIF4A, TOP2A, AURKA, and CCNA2, which were significantly upregulated in Oncomine and TCGA LUAD datasets, and were verified by qRT-PCR in our clinical samples. We determined the overall and disease-free survival analysis of the seven hub genes using GEPIA. We further found that CENPF, DLGAP5, and KIF4A expressions were positively correlated with clinical stage. In LUAD cell lines, proliferation and migration were inhibited and apoptosis was promoted by knocking down KIF4A expression.

CONCLUSION

We have identified new DEGs and functional pathways involved in LUAD. KIF4A, as a hub gene, promoted the progression of LUAD and might represent a potential therapeutic target for molecular cancer therapy.

Repurposing potential of FDA-approved and investigational drugs for COVID-19 targeting SARS-CoV-2 spike and main protease and validation by machine learning algorithm

The present study aimed to assess the repurposing potential of existing antiviral drug candidates (FDA-approved and investigational) against SARS-CoV-2 target proteins that facilitates viral entry and replication into the host body. To evaluate molecular affinities between antiviral drug candidates and SARS-CoV-2 associated target proteins such as spike protein (S) and main protease (M^pro ), a molecular interaction simulation was performed by docking software (MVD) and subsequently the applicability score was calculated by machine learning algorithm. Furthermore, the STITCH algorithm was used to predict the pharmacology network involving multiple pathways of active drug candidate(s). Pharmacophore features of active drug(s) molecule was also determined to predict structure-activity relationship (SAR). The molecular interaction analysis showed that cordycepin has strong binding affinities with S protein (-180) and M^pro proteins (-205) which were relatively highest among other drug candidates used. Interestingly, compounds with low IC₅₀ showed high binding energy. Furthermore, machine learning algorithm also revealed high applicability scores (0.42-0.47) of cordycepin. It is worth mentioning that the pharmacology network depicted the involvement of cordycepin in different pathways associated with bacterial and viral diseases including tuberculosis, hepatitis B, influenza A, viral myocarditis, and herpes simplex infection. The embedded pharmacophore features with cordycepin also suggested strong SAR. Cordycepin's anti-SARS-CoV-2 activity indicated 65% (E-gene) and 42% (N-gene) viral replication inhibition after 48h of treatment. Since, cordycepin has both preclinical and clinical evidences on antiviral activity, in addition the present findings further validate and suggest repurposing potential of cordycepin against COVID-19.

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).

Genome-wide pathogenesis interpretation using a heat diffusion-based systems genetics method and implications for gene function annotation

Background

Genetics is best dedicated to interpreting pathogenesis and revealing gene functions. The past decade has witnessed unprecedented progress in genetics, particularly in genome‐wide identification of disorder variants through Genome‐Wide Association Studies (GWAS) and Phenome‐Wide Association Studies (PheWAS). However, it is still a great challenge to use GWAS/PheWAS‐derived data to elucidate pathogenesis.

Methods

In this study, we used HotNet2, a heat diffusion‐based systems genetics algorithm, to calculate the networks for disease genes obtained from GWAS and PheWAS, with an attempt to get deeper insights into disease pathogenesis at a molecular level.

Results

Through HotNet2 calculation, significant networks for 202 (for GWAS) and 167 (for PheWAS) types of diseases were identified and evaluated, respectively. The GWAS‐derived disease networks exhibit a stronger biomedical relevance than PheWAS counterparts. Therefore, the GWAS‐derived networks were used for pathogenesis interpretation by integrating the accumulated biomedical information. As a result, the pathogenesis for 64 diseases was elucidated in terms of mutation‐caused abnormal transcriptional regulation, and 47 diseases were preliminarily interpreted in terms of mutation‐caused varied protein‐protein interactions. In addition, 3,802 genes (including 46 function‐unknown genes) were assigned with new functions by disease network information, some of which were validated through mice gene knockout experiments.

Conclusions

Systems genetics algorithm HotNet2 can efficiently establish genotype‐phenotype links at the level of biological networks. Compared with original GWAS/PheWAS results, HotNet2‐calculated disease‐gene associations have stronger biomedical significance, hence provide better interpretations for the pathogenesis of genome‐wide variants, and offer new insights into gene functions as well. These results are also helpful in drug development.

Proteinarium: Multi-sample protein-protein interaction analysis and visualization tool

We posit the likely architecture of complex diseases is that subgroups of patients share variants in genes in specific networks which are sufficient to give rise to a shared phenotype. We developed Proteinarium, a multi-sample protein-protein interaction (PPI) tool, to identify clusters of patients with shared gene networks. Proteinarium converts user defined seed genes to protein symbols and maps them onto the STRING interactome. A PPI network is built for each sample using Dijkstra's algorithm. Pairwise similarity scores are calculated to compare the networks and cluster the samples. A layered graph of PPI networks for the samples in any cluster can be visualized. To test this newly developed analysis pipeline, we reanalyzed publicly available data sets, from which modest outcomes had previously been achieved. We found significant clusters of patients with unique genes which enhanced the findings in the original study.

Protein-Protein interactions uncover candidate 'core genes' within omnigenic disease networks

Genome wide association studies (GWAS) of human diseases have generally identified many loci associated with risk with relatively small effect sizes. The omnigenic model attempts to explain this observation by suggesting that diseases can be thought of as networks, where genes with direct involvement in disease-relevant biological pathways are named ‘core genes’, while peripheral genes influence disease risk via their interactions or regulatory effects on core genes. Here, we demonstrate a method for identifying candidate core genes solely from genes in or near disease-associated SNPs (GWAS hits) in conjunction with protein-protein interaction network data. Applied to 1,381 GWAS studies from 5 ancestries, we identify a total of 1,865 candidate core genes in 343 GWAS studies. Our analysis identifies several well-known disease-related genes that are not identified by GWAS, including BRCA1 in Breast Cancer, Amyloid Precursor Protein (APP) in Alzheimer’s Disease, INS in A1C measurement and Type 2 Diabetes, and PCSK9 in LDL cholesterol, amongst others. Notably candidate core genes are preferentially enriched for disease relevance over GWAS hits and are enriched for both Clinvar pathogenic variants and known drug targets—consistent with the predictions of the omnigenic model. We subsequently use parent term annotations provided by the GWAS catalog, to merge related GWAS studies and identify candidate core genes in over-arching disease processes such as cancer–where we identify 109 candidate core genes.

A recent theory suggests that only a small number of genes underpin the biology of a disease, these genes are called ‘core genes’, and for most diseases, these core genes remain unknown. The suggested methods for finding them requires complex and expensive experiments. We reasoned that if we merge currently available datasets in smart ways, we may be able to uncover these ‘core genes’. Our method finds “hub” proteins by merging lists of genes previously linked with disease to information on how proteins interact with each other. We found that many of these hub proteins have central roles in disease, such as insulin for both A1C measurement and Type 2 Diabetes, BRCA1 in Breast cancer, and Amyloid Precursor Protein in Alzheimer’s Disease. We think these ‘hub’ proteins are candidate ‘core genes’, and offer our method as a way to find ‘core genes’ by utilizing publicly available reference datasets.

A Cancer Gene Module Mining Method Based on Bio-Network of Multi-Omics Gene Groups

The initiation, promotion and progression of cancer are highly associated to the environment a human lives in as well as individual genetic factors. In view of the dangers to life and health caused by this abnormally complex systemic disease, many top scientific research institutions around the world have been actively carrying out research in order to discover the pathogenic mechanisms driving cancer occurrence and development. The emergence of high-throughput sequencing technology has greatly advanced oncology research and given rise to the revelation of important oncogenes and the interrelationship among them. Here, we have studied heterogeneous multi-level data within a context of integrated data, and scientifically introduced lncRNA omics data to construct multi-omics bio-network models, allowing the screening of key cancer-related gene groups. We propose a compactness clustering algorithm based on corrected cumulative rank scores, which uses the functional similarity between groups of genes as a distance measure to excavate key gene modules for abnormal regulation contained in gene groups through clustering. We also conducted a survival analysis using our results and found that our model could divide groups of different levels very well. The results also demonstrate that the integration of multi-omics biological data, key gene modules and their dysregulated gene groups can be discovered, which is crucial for cancer research.