Molecular networks in Network Medicine: Development and applications

Mechanistic prediction and validation of Brevilin A Therapeutic effects in Lung Cancer

BACKGROUND

Traditional Chinese medicine (TCM) has been found widespread application in neoplasm treatment, yielding promising therapeutic candidates. Previous studies have revealed the anti-cancer properties of Brevilin A, a naturally occurring sesquiterpene lactone derived from Centipeda minima (L.) A.Br. (C. minima), a TCM herb, specifically against lung cancer. However, the underlying mechanisms of its effects remain elusive. This study employs network pharmacology and experimental analyses to unravel the molecular mechanisms of Brevilin A in lung cancer.

METHODS

The Batman-TCM, Swiss Target Prediction, Pharmmapper, SuperPred, and BindingDB databases were screened to identify Brevilin A targets. Lung cancer-related targets were sourced from GEO, Genecards, OMIM, TTD, and Drugbank databases. Utilizing Cytoscape software, a protein-protein interaction (PPI) network was established. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene set enrichment analysis (GSEA), and gene-pathway correlation analysis were conducted using R software. To validate network pharmacology results, molecular docking, molecular dynamics simulations, and in vitro experiments were performed.

RESULTS

We identified 599 Brevilin A-associated targets and 3864 lung cancer-related targets, with 155 overlapping genes considered as candidate targets for Brevilin A against lung cancer. The PPI network highlighted STAT3, TNF, HIF1A, PTEN, ESR1, and MTOR as potential therapeutic targets. GO and KEGG analyses revealed 2893 enriched GO terms and 157 enriched KEGG pathways, including the PI3K-Akt signaling pathway, FoxO signaling pathway, and HIF-1 signaling pathway. GSEA demonstrated a close association between hub genes and lung cancer. Gene-pathway correlation analysis indicated significant associations between hub genes and the cellular response to hypoxia pathway. Molecular docking and dynamics simulations confirmed Brevilin A's interaction with PTEN and HIF1A, respectively. In vitro experiments demonstrated Brevilin A-induced dose- and time-dependent cell death in A549 cells. Notably, Brevilin A treatment significantly reduced HIF-1α mRNA expression while increasing PTEN mRNA levels.

CONCLUSIONS

This study demonstrates that Brevilin A exerts anti-cancer effects in treating lung cancer through a multi-target and multi-pathway manner, with the HIF pathway potentially being involved. These results lay a theoretical foundation for the prospective clinical application of Brevilin A.

Transcriptomics analysis of long non-coding RNAs in smooth muscle cells from patients with peripheral artery disease and diabetes mellitus

Diabetes mellitus (DM) is a significant risk factor for peripheral arterial disease (PAD), and PAD is an independent predictor of cardiovascular disorders (CVDs). Growing evidence suggests that long non-coding RNAs (lncRNAs) significantly contribute to disease development and underlying complications, particularly affecting smooth muscle cells (SMCs). So far, no study has focused on transcriptome analysis of lncRNAs in PAD patients with and without DM. Tissue samples were obtained from our Vascular Biobank. Due to the sample's heterogeneity, expression analysis of lncRNAs in whole tissue detected only ACTA2-AS1 with a 4.9-fold increase in PAD patients with DM. In contrast, transcriptomics of SMCs revealed 28 lncRNAs significantly differentially expressed between PAD with and without DM (FDR < 0.1). Sixteen lncRNAs were of unknown function, six were described in cancer, one connected with macrophages polarisation, and four were associated with CVDs, mainly with SMC function and phenotypic switch (NEAT1, MIR100HG, HIF1A-AS3, and MRI29B2CHG). The enrichment analysis detected additional lncRNAs H19, CARMN, FTX, and MEG3 linked with DM. Our study revealed several lncRNAs in diabetic PAD patients associated with the physiological function of SMCs. These lncRNAs might serve as potential therapeutic targets to improve the function of SMCs within the diseased tissue and, thus, the clinical outcome.

Finding new analgesics: Computational pharmacology faces drug discovery challenges

Despite the worldwide prevalence and huge burden of pain, pain is an undertreated phenomenon. Currently used analgesics have several limitations regarding their efficacy and safety. The discovery of analgesics possessing a novel mechanism of action has faced multiple challenges, including a limited understanding of biological processes underpinning pain and analgesia and poor animal-to-human translation. Computational pharmacology is currently employed to face these challenges. In this review, we discuss the theory, methods, and applications of computational pharmacology in pain research. Computational pharmacology encompasses a wide variety of theoretical concepts and practical methodological approaches, with the overall aim of gaining biological insight through data acquisition and analysis. Data are acquired from patients or animal models with pain or analgesic treatment, at different levels of biological organization (molecular, cellular, physiological, and behavioral). Distinct methodological algorithms can then be used to analyze and integrate data. This helps to facilitate the identification of biological molecules and processes associated with pain phenotype, build quantitative models of pain signaling, and extract translatable features between humans and animals. However, computational pharmacology has several limitations, and its predictions can provide false positive and negative findings. Therefore, computational predictions are required to be validated experimentally before drawing solid conclusions. In this review, we discuss several case study examples of combining and integrating computational tools with experimental pain research tools to meet drug discovery challenges.

Prioritization of genes for translation: a computational approach

INTRODUCTION

Gene identification for genetic diseases is critical for the development of new diagnostic approaches and personalized treatment options. Prioritization of gene translation is an important consideration in the molecular biology field, allowing researchers to focus on the most promising candidates for further investigation.

AREAS COVERED

In this paper, we discussed different approaches to prioritize genes for translation, including the use of computational tools and machine learning algorithms, as well as experimental techniques such as knockdown and overexpression studies. We also explored the potential biases and limitations of these approaches and proposed strategies to improve the accuracy and reliability of gene prioritization methods. Although numerous computational methods have been developed for this purpose, there is a need for computational methods that incorporate tissue-specific information to enable more accurate prioritization of candidate genes. Such methods should provide tissue-specific predictions, insights into underlying disease mechanisms, and more accurate prioritization of genes.

EXPERT OPINION

Using advanced computational tools and machine learning algorithms to prioritize genes, we can identify potential targets for therapeutic intervention of complex diseases. This represents an up-and-coming method for drug development and personalized medicine.

HHIP protein interactions in lung cells provide insight into COPD pathogenesis

Chronic obstructive pulmonary disease (COPD) is the third leading cause of death worldwide. The primary causes of COPD are environmental, including cigarette smoking; however, genetic susceptibility also contributes to COPD risk. Genome-Wide Association Studies (GWASes) have revealed more than 80 genetic loci associated with COPD, leading to the identification of multiple COPD GWAS genes. However, the biological relationships between the identified COPD susceptibility genes are largely unknown. Genes associated with a complex disease are often in close network proximity, i.e. their protein products often interact directly with each other and/or similar proteins. In this study, we use affinity purification mass spectrometry (AP-MS) to identify protein interactions with HHIP , a well-established COPD GWAS gene which is part of the sonic hedgehog pathway, in two disease-relevant lung cell lines (IMR90 and 16HBE). To better understand the network neighborhood of HHIP , its proximity to the protein products of other COPD GWAS genes, and its functional role in COPD pathogenesis, we create HUBRIS, a protein-protein interaction network compiled from 8 publicly available databases. We identified both common and cell type-specific protein-protein interactors of HHIP. We find that our newly identified interactions shorten the network distance between HHIP and the protein products of several COPD GWAS genes, including DSP, MFAP2, TET2 , and FBLN5 . These new shorter paths include proteins that are encoded by genes involved in extracellular matrix and tissue organization. We found and validated interactions to proteins that provide new insights into COPD pathobiology, including CAVIN1 (IMR90) and TP53 (16HBE). The newly discovered HHIP interactions with CAVIN1 and TP53 implicate HHIP in response to oxidative stress.

Transgenerational Epigenetic Inheritance of Cardiovascular Diseases: A Network Medicine Perspective

INTRODUCTION

The ability to identify early epigenetic signatures underlying the inheritance of cardiovascular risk, including trans- and intergenerational effects, may help to stratify people before cardiac symptoms occur.

METHODS

Prospective and retrospective cohorts and case-control studies focusing on DNA methylation and maternal/paternal effects were searched in Pubmed from 1997 to 2023 by using the following keywords: DNA methylation, genomic imprinting, and network analysis in combination with transgenerational/intergenerational effects.

RESULTS

Maternal and paternal exposures to traditional cardiovascular risk factors during critical temporal windows, including the preconceptional period or early pregnancy, may perturb the plasticity of the epigenome (mainly DNA methylation) of the developing fetus especially at imprinted loci, such as the insulin-like growth factor type 2 (IGF2) gene. Thus, the epigenome is akin to a "molecular archive" able to memorize parental environmental insults and predispose an individual to cardiovascular diseases onset in later life. Direct evidence for human transgenerational epigenetic inheritance (at least three generations) of cardiovascular risk is lacking but it is supported by epidemiological studies. Several blood-based association studies showed potential intergenerational epigenetic effects (single-generation studies) which may mediate the transmittance of cardiovascular risk from parents to offspring.

DISCUSSION

In this narrative review, we discuss some relevant examples of trans- and intergenerational epigenetic associations with cardiovascular risk. In our perspective, we propose three network-oriented approaches which may help to clarify the unsolved issues regarding transgenerational epigenetic inheritance of cardiovascular risk and provide potential early biomarkers for primary prevention.

Construction and evaluation of endometriosis diagnostic prediction model and immune infiltration based on efferocytosis-related genes

Background: Endometriosis (EM) is a long-lasting inflammatory disease that is difficult to treat and prevent. Existing research indicates the significance of immune infiltration in the progression of EM. Efferocytosis has an important immunomodulatory function. However, research on the identification and clinical significance of efferocytosis-related genes (EFRGs) in EM is sparse. Methods: The EFRDEGs (differentially expressed efferocytosis-related genes) linked to datasets associated with endometriosis were thoroughly examined utilizing the Gene Expression Omnibus (GEO) and GeneCards databases. The construction of the protein-protein interaction (PPI) and transcription factor (TF) regulatory network of EFRDEGs ensued. Subsequently, machine learning techniques including Univariate logistic regression, LASSO, and SVM classification were applied to filter and pinpoint diagnostic biomarkers. To establish and assess the diagnostic model, ROC analysis, multivariate regression analysis, nomogram, and calibration curve were employed. The CIBERSORT algorithm and single-cell RNA sequencing (scRNA-seq) were employed to explore immune cell infiltration, while the Comparative Toxicogenomics Database (CTD) was utilized for the identification of potential therapeutic drugs for endometriosis. Finally, immunohistochemistry (IHC) and reverse transcription quantitative polymerase chain reaction (RT-qPCR) were utilized to quantify the expression levels of biomarkers in clinical samples of endometriosis. Results: Our findings revealed 13 EFRDEGs associated with EM, and the LASSO and SVM regression model identified six hub genes (ARG2, GAS6, C3, PROS1, CLU, and FGL2). Among these, ARG2, GAS6, and C3 were confirmed as diagnostic biomarkers through multivariate logistic regression analysis. The ROC curve analysis of GSE37837 (AUC = 0.627) and GSE6374 (AUC = 0.635), along with calibration and DCA curve assessments, demonstrated that the nomogram built on these three biomarkers exhibited a commendable predictive capacity for the disease. Notably, the ratio of nine immune cell types exhibited significant differences between eutopic and ectopic endometrial samples, with scRNA-seq highlighting M0 Macrophages, Fibroblasts, and CD8 Tex cells as the cell populations undergoing the most substantial changes in the three biomarkers. Additionally, our study predicted seven potential medications for EM. Finally, the expression levels of the three biomarkers in clinical samples were validated through RT-qPCR and IHC, consistently aligning with the results obtained from the public database. Conclusion: we identified three biomarkers and constructed a diagnostic model for EM in this study, these findings provide valuable insights for subsequent mechanistic research and clinical applications in the field of endometriosis.

A framework for considering prior information in network-based approaches to omics data analysis

For decades, molecular biologists have been uncovering the mechanics of biological systems. Efforts to bring their findings together have led to the development of multiple databases and information systems that capture and present pathway information in a computable network format. Concurrently, the advent of modern omics technologies has empowered researchers to systematically profile cellular processes across different modalities. Numerous algorithms, methodologies, and tools have been developed to use prior knowledge networks (PKNs) in the analysis of omics datasets. Interestingly, it has been repeatedly demonstrated that the source of prior knowledge can greatly impact the results of a given analysis. For these methods to be successful it is paramount that their selection of PKNs is amenable to the data type and the computational task they aim to accomplish. Here we present a five-level framework that broadly describes network models in terms of their scope, level of detail, and ability to inform causal predictions. To contextualize this framework, we review a handful of network-based omics analysis methods at each level, while also describing the computational tasks they aim to accomplish.

The potential DNA methylation markers of cardiovascular disease in patients with type 2 diabetes

BACKGROUND

DNA methylation is associated with cardiovascular (CV) disease. However, in type 2 diabetes (T2D) patients, the role of gene methylation in the development of CV disease is under-studied. We aimed to identify the CV disease-related DNA methylation loci in patients with T2D and to explore the potential pathways underlying the development of CV disease using a two-stage design.

METHODS

The participants were from the Jinan Diabetes Cohort Study (JNDCS), an ongoing longitudinal study designed to evaluate the development of CV risk in patients with T2D. In the discovery cohort, 10 diabetic patients with CV events at baseline were randomly selected as the case group, and another 10 diabetic patients without CV events were matched for sex, age, smoking status, and body mass index as the control group. In 1438 T2D patients without CV disease at baseline, 210 patients with CV events were identified after a mean 6.5-year follow-up. Of whom, 100 patients who experienced CV events during the follow-up were randomly selected as cases, and 100 patients who did not have CV events were randomly selected as the control group in the validation cohort. Reduced representation bisulfite sequencing and Targeted Bisulfite Sequencing were used to measure the methylation profiles in the discovery and validation cohort, respectively.

RESULTS

In the discover cohort, 127 DMRs related to CV disease were identified in T2D patients. Further, we validated 23 DMRs mapped to 25 genes, of them, 4 genes (ARSG, PNPLA6, NEFL, and CRYGEP) for the first time were reported. There was evidence that the addition of DNA methylation data improved the prediction performance of CV disease in T2D patients. Pathway analysis identified some significant signaling pathways involved in CV comorbidities, T2D, and inflammation.

CONCLUSIONS

In this study, we identified 23 DMRs mapped to 25 genes associated with CV disease in T2D patients, of them, 4 DMRs for the first time were reported. DNA methylation testing may help identify a high CV-risk population in T2D patients.

Adjustment of spurious correlations in co-expression measurements from RNA-Sequencing data

MOTIVATION

Gene co-expression measurements are widely used in computational biology to identify coordinated expression patterns across a group of samples. Coordinated expression of genes may indicate that they are controlled by the same transcriptional regulatory program, or involved in common biological processes. Gene co-expression is generally estimated from RNA-Sequencing data, which are commonly normalized to remove technical variability. Here, we demonstrate that certain normalization methods, in particular quantile-based methods, can introduce false-positive associations between genes. These false-positive associations can consequently hamper downstream co-expression network analysis. Quantile-based normalization can, however, be extremely powerful. In particular, when preprocessing large-scale heterogeneous data, quantile-based normalization methods such as smooth quantile normalization can be applied to remove technical variability while maintaining global differences in expression for samples with different biological attributes.

RESULTS

We developed SNAIL (Smooth-quantile Normalization Adaptation for the Inference of co-expression Links), a normalization method based on smooth quantile normalization specifically designed for modeling of co-expression measurements. We show that SNAIL avoids formation of false-positive associations in co-expression as well as in downstream network analyses. Using SNAIL, one can avoid arbitrary gene filtering and retain associations to genes that only express in small subgroups of samples. This highlights the method's potential future impact on network modeling and other association-based approaches in large-scale heterogeneous data.

AVAILABILITY AND IMPLEMENTATION

The implementation of the SNAIL algorithm and code to reproduce the analyses described in this work can be found in the GitHub repository https://github.com/kuijjerlab/PySNAIL.

Epigenetic hallmarks in pulmonary fibrosis: New advances and perspectives

Epigenetics indicates that certain phenotypes of an organism can undergo heritable changes in the absence of changes in the genetic DNA sequence. Many studies have shown that epigenetic patterns play an important role in the lung and lung diseases. Pulmonary fibrosis (PF) is also a type of lung disease. PF is an end-stage change of a large group of lung diseases, characterized by fibroblast proliferation and massive accumulation of extracellular matrix, accompanied by inflammatory injury and histological destruction, that is, structural abnormalities caused by abnormal repair of normal alveolar tissue. It causes loss of lung function in patients with multiple complex diseases, leading to respiratory failure and subsequent death. However, current treatment options for IPF are very limited and no drugs have been shown to significantly prolong the survival of patients. Therefore, based on a systematic understanding of the disease mechanisms of PF, this review integrates the role of epigenetics in the development and course of PF, describes preventive and potential therapeutic targets for PF, and provides a theoretical basis for further exploration of the mechanisms of PF.

DapBCH: a disease association prediction model Based on Cross-species and Heterogeneous graph embedding

The study of comorbidity can provide new insights into the pathogenesis of the disease and has important economic significance in the clinical evaluation of treatment difficulty, medical expenses, length of stay, and prognosis of the disease. In this paper, we propose a disease association prediction model DapBCH, which constructs a cross-species biological network and applies heterogeneous graph embedding to predict disease association. First, we combine the human disease-gene network, mouse gene-phenotype network, human-mouse homologous gene network, and human protein-protein interaction network to reconstruct a heterogeneous biological network. Second, we apply heterogeneous graph embedding based on meta-path aggregation to generate the feature vector of disease nodes. Finally, we employ link prediction to obtain the similarity of disease pairs. The experimental results indicate that our model is highly competitive in predicting the disease association and is promising for finding potential disease associations.

Mediator complex in neurological disease

Neurological diseases, including traumatic brain injuries, stroke (haemorrhagic and ischemic), and inherent neurodegenerative diseases cause acquired disability in humans, representing a leading cause of death worldwide. The Mediator complex (MED) is a large, evolutionarily conserved multiprotein that facilities the interaction between transcription factors and RNA Polymerase II in eukaryotes. Some MED subunits have been found altered in the brain, although their specific functions in neurodegenerative diseases are not fully understood. Mutations in MED subunits were associated with a wide range of genetic diseases for MED12, MED13, MED13L, MED20, MED23, MED25, and CDK8 genes. In addition, MED12 and MED23 were deregulated in the Alzheimer's Disease. Interestingly, most of the genomic mutations have been found in the subunits of the kinase module. To date, there is only one evidence on MED1 involvement in post-stroke cognitive deficits. Although the underlying neurodegenerative disorders may be different, we are confident that the signal cascades of the biological-cognitive mechanisms of brain adaptation, which begin after brain deterioration, may also differ. Here, we analysed relevant studies in English published up to June 2023. They were identified through a search of electronic databases including PubMed, Medline, EMBASE and Scopus, including search terms such as "Mediator complex", "neurological disease", "brains". Thematic content analysis was conducted to collect and summarize all studies demonstrating MED alteration to understand the role of this central transcriptional regulatory complex in the brain. Improved and deeper knowledge of the regulatory mechanisms in neurological diseases can increase the ability of physicians to predict onset and progression, thereby improving diagnostic care and providing appropriate treatment decisions.

Multidimensional Landscape of SA-AKI Revealed by Integrated Proteomics and Metabolomics Analysis

Sepsis-associated acute kidney injury (SA-AKI) is a severe and life-threatening condition with high morbidity and mortality among emergency patients, and it poses a significant risk of chronic renal failure. Clinical treatments for SA-AKI remain reactive and non-specific, lacking effective diagnostic biomarkers or treatment targets. In this study, we established an SA-AKI mouse model using lipopolysaccharide (LPS) and performed proteomics and metabolomics analyses. A variety of bioinformatic analyses, including gene set enrichment analysis (GSEA), weighted gene co-expression network analysis (WGCNA), protein and protein interactions (PPI), and MetaboAnalyst analysis, were conducted to investigate the key molecules of SA-AKI. Integrated proteomics and metabolomics analysis revealed that sepsis led to impaired renal mitochondrial function and metabolic disorders. Immune-related pathways were found to be activated in kidneys upon septic infection. The catabolic products of polyamines accumulated in septic kidneys. Overall, our integrated analysis provides a multidimensional understanding of SA-AKI and identifies potential pathways for this condition.

From time-series transcriptomics to gene regulatory networks: A review on inference methods

Inference of gene regulatory networks has been an active area of research for around 20 years, leading to the development of sophisticated inference algorithms based on a variety of assumptions and approaches. With the ever increasing demand for more accurate and powerful models, the inference problem remains of broad scientific interest. The abstract representation of biological systems through gene regulatory networks represents a powerful method to study such systems, encoding different amounts and types of information. In this review, we summarize the different types of inference algorithms specifically based on time-series transcriptomics, giving an overview of the main applications of gene regulatory networks in computational biology. This review is intended to give an updated reference of regulatory networks inference tools to biologists and researchers new to the topic and guide them in selecting the appropriate inference method that best fits their questions, aims, and experimental data.

A review of multi-omics data integration through deep learning approaches for disease diagnosis, prognosis, and treatment

Accurate diagnosis is the key to providing prompt and explicit treatment and disease management. The recognized biological method for the molecular diagnosis of infectious pathogens is polymerase chain reaction (PCR). Recently, deep learning approaches are playing a vital role in accurately identifying disease-related genes for diagnosis, prognosis, and treatment. The models reduce the time and cost used by wet-lab experimental procedures. Consequently, sophisticated computational approaches have been developed to facilitate the detection of cancer, a leading cause of death globally, and other complex diseases. In this review, we systematically evaluate the recent trends in multi-omics data analysis based on deep learning techniques and their application in disease prediction. We highlight the current challenges in the field and discuss how advances in deep learning methods and their optimization for application is vital in overcoming them. Ultimately, this review promotes the development of novel deep-learning methodologies for data integration, which is essential for disease detection and treatment.

A network approach to define the predictive role of immune profile on tumor response and toxicity of anti PD-1 single agent immunotherapy in patients with solid tumors

Background

The immune profile of each patient could be considered as a portrait of the fitness of his/her own immune system. The predictive role of the immune profile in immune-related toxicities (irAEs) development and tumour response to treatment was investigated.

Methods

A prospective, multicenter study evaluating, through a multiplex assay, the soluble immune profile at the baseline of 53 patients with advanced cancer, treated with immunotherapy as single agent was performed. Four connectivity heat maps and networks were obtained by calculating the Spearman correlation coefficients for each group: responder patients who developed cumulative toxicity (R-T), responders who did not develop cumulative toxicity (R-NT), non-responders who developed cumulative toxicity (NR-T), non-responders who did not develop cumulative toxicity (NR-NT).

Results

A statistically significant up-regulation of IL-17A, sCTLA4, sCD80, I-CAM-1, sP-Selectin and sEselectin in NR-T was detected. A clear loss of connectivity of most of the soluble immune checkpoints and cytokines characterized the immune profile of patients with toxicity, while an inversion of the correlation for ICAM-1 and sP-selectin was observed in NR-T. Four connectivity networks were built for each group. The highest number of connections characterized the NR-T.

Conclusions

A connectivity network of immune dysregulation was defined for each subgroup of patients, regardless of tumor type. In patients with the worst prognosis (NR-T) the peculiar connectivity model could facilitate their early and timely identification, as well as the design of a personalized treatment approach to improve outcomes or prevent irAEs.

Gauging innovation and health impact from biomedical research: survey results and interviews with recipients of EU-funding in the fields of Alzheimer's disease, breast cancer and prostate cancer

Biomedical research on Alzheimer's disease (AD), breast cancer (BC) and prostate cancer (PC) has globally improved our understanding of the etiopathological mechanisms underlying the onset of these diseases, often with the goal to identify associated genetic and environmental risk factors and develop new medicines. However, the prevalence of these diseases and failure rate in drug development remain high. Being able to retrospectively monitor the major scientific breakthroughs and impact of such investment endeavors is important to re-address funding strategies if and when needed. The EU has supported research into those diseases via its successive framework programmes for research, technological development and innovation. The European Commission (EC) has already undertaken several activities to monitor research impact. As an additional contribution, the EC Joint Research Centre (JRC) launched in 2020 a survey addressed to former and current participants of EU-funded research projects in the fields of AD, BC and PC, with the aim to understand how EU-funded research has contributed to scientific innovation and societal impact, and how the selection of the experimental models may have underpinned the advances made. Further feedback was also gathered through in-depth interviews with some selected survey participants representative of the diverse pre-clinical models used in the EU-funded projects. A comprehensive analysis of survey replies, complemented with the information derived from the interviews, has recently been published in a Synopsis report. Here we discuss the main findings of this analysis and propose a set of priority actions that could be considered to help improving the translation of scientific innovation of biomedical research into societal impact.

Using knowledge graphs to infer gene expression in plants

Introduction

Climate change is already affecting ecosystems around the world and forcing us to adapt to meet societal needs. The speed with which climate change is progressing necessitates a massive scaling up of the number of species with understood genotype-environment-phenotype (G×E×P) dynamics in order to increase ecosystem and agriculture resilience. An important part of predicting phenotype is understanding the complex gene regulatory networks present in organisms. Previous work has demonstrated that knowledge about one species can be applied to another using ontologically-supported knowledge bases that exploit homologous structures and homologous genes. These types of structures that can apply knowledge about one species to another have the potential to enable the massive scaling up that is needed through in silico experimentation.

Methods

We developed one such structure, a knowledge graph (KG) using information from Planteome and the EMBL-EBI Expression Atlas that connects gene expression, molecular interactions, functions, and pathways to homology-based gene annotations. Our preliminary analysis uses data from gene expression studies in Arabidopsis thaliana and Populus trichocarpa plants exposed to drought conditions.

Results

A graph query identified 16 pairs of homologous genes in these two taxa, some of which show opposite patterns of gene expression in response to drought. As expected, analysis of the upstream cis-regulatory region of these genes revealed that homologs with similar expression behavior had conserved cis-regulatory regions and potential interaction with similar trans-elements, unlike homologs that changed their expression in opposite ways.

Discussion

This suggests that even though the homologous pairs share common ancestry and functional roles, predicting expression and phenotype through homology inference needs careful consideration of integrating cis and trans-regulatory components in the curated and inferred knowledge graph.

Mitochondrial quality control in cardiac fibrosis: Epigenetic mechanisms and therapeutic strategies

Cardiac fibrosis (CF) is considered an ultimate common pathway of a wide variety of heart diseases in response to diverse pathological and pathophysiological stimuli. Mitochondria are characterized as isolated organelles with a double-membrane structure, and they primarily contribute to and maintain highly dynamic energy and metabolic networks whose distribution and structure exert potent support for cellular properties and performance. Because the myocardium is a highly oxidative tissue with high energy demands to continuously pump blood, mitochondria are the most abundant organelles within mature cardiomyocytes, accounting for up to one-third of the total cell volume, and play an essential role in maintaining optimal performance of the heart. Mitochondrial quality control (MQC), including mitochondrial fusion, fission, mitophagy, mitochondrial biogenesis, and mitochondrial metabolism and biosynthesis, is crucial machinery that modulates cardiac cells and heart function by maintaining and regulating the morphological structure, function and lifespan of mitochondria. Certain investigations have focused on mitochondrial dynamics, including manipulating and maintaining the dynamic balance of energy demand and nutrient supply, and the resultant findings suggest that changes in mitochondrial morphology and function may contribute to bioenergetic adaptation during cardiac fibrosis and pathological remodeling. In this review, we discuss the function of epigenetic regulation and molecular mechanisms of MQC in the pathogenesis of CF and provide evidence for targeting MQC for CF. Finally, we discuss how these findings can be applied to improve the treatment and prevention of CF.

Recent advances in mass spectrometry-based computational metabolomics

The computational metabolomics field brings together computer scientists, bioinformaticians, chemists, clinicians, and biologists to maximize the impact of metabolomics across a wide array of scientific and medical disciplines. The field continues to expand as modern instrumentation produces datasets with increasing complexity, resolution, and sensitivity. These datasets must be processed, annotated, modeled, and interpreted to enable biological insight. Techniques for visualization, integration (within or between omics), and interpretation of metabolomics data have evolved along with innovation in the databases and knowledge resources required to aid understanding. In this review, we highlight recent advances in the field and reflect on opportunities and innovations in response to the most pressing challenges. This review was compiled from discussions from the 2022 Dagstuhl seminar entitled "Computational Metabolomics: From Spectra to Knowledge".

Three genes expressed in relation to lipid metabolism considered as potential biomarkers for the diagnosis and treatment of diabetic peripheral neuropathy

Diabetic neuropathy is one of the most common chronic complications and is present in approximately 50% of diabetic patients. A bioinformatic approach was used to analyze candidate genes involved in diabetic distal symmetric polyneuropathy and their potential mechanisms. GSE95849 was downloaded from the Gene Expression Omnibus database for differential analysis, together with the identified diabetic peripheral neuropathy-associated genes and the three major metabolism-associated genes in the CTD database to obtain overlapping Differentially Expressed Genes (DEGs). Gene Set Enrichment Analysis and Functional Enrichment Analysis were performed. Protein-Protein Interaction and hub gene networks were constructed using the STRING database and Cytoscape software. The expression levels of target genes were evaluated using GSE24290 samples, followed by Receiver operating characteristic, curve analysis. And Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were performed on the target genes. Finally, mRNA-miRNA networks were constructed. A total of 442 co-expressed DEGs were obtained through differential analysis, of which 353 expressed up-regulated genes and 89 expressed down-regulated genes. The up-regulated DEGs were involved in 742 GOs and 10 KEGG enrichment results, mainly associated with lipid metabolism-related pathways, TGF-β receptor signaling pathway, lipid transport, and PPAR signaling pathway. A total of 4 target genes (CREBBP, EP300, ME1, CD36) were identified. Analysis of subject operating characteristic curves indicated that CREBBP (AUC = 1), EP300 (AUC = 0.917), ME1 (AUC = 0.944) and CD36 (AUC = 1) may be candidate serum biomarkers for DPN. Conclusion: Diabetic peripheral neuropathy pathogenesis and progression is caused by multiple pathways, which also provides clinicians with potential therapeutic tools.

Network modeling helps to tackle the complexity of drug-disease systems

From the (patho)physiological point of view, diseases can be considered as emergent properties of living systems stemming from the complexity of these systems. Complex systems display some typical features, including the presence of emergent behavior and the organization in successive hierarchic levels. Drug treatments increase this complexity scenario, and from some years the use of network models has been introduced to describe drug-disease systems and to make predictions about them with regard to several aspects related to drug discovery. Here, we review some recent examples thereof with the aim to illustrate how network science tools can be very effective in addressing both tasks. We will examine the use of bipartite networks that lead to the important concept of "disease module", as well as the introduction of more articulated models, like multi-scale and multiplex networks, able to describe disease systems at increasing levels of organization. Examples of predictive models will then be discussed, considering both those that exploit approaches purely based on graph theory and those that integrate machine learning methods. A short account of both kinds of methodological applications will be provided. Finally, the point will be made on the present situation of modeling complex drug-disease systems highlighting some open issues. This article is categorized under: Neurological Diseases > Computational Models Infectious Diseases > Computational Models Cardiovascular Diseases > Computational Models.

Pulmonary arterial hypertension associated with congenital heart disease: An omics study

Pulmonary arterial hypertension associated with congenital heart disease (PAH-CHD) is a severely progressive condition with uncertain physiological course. Hence, it has become increasingly relevant to clarify the specific mechanisms of molecular modification, which is crucial to identify more treatment strategies. With the rapid development of high-throughput sequencing, omics technology gives access to massive experimental data and advanced techniques for systems biology, permitting comprehensive assessment of disease occurrence and progression. In recent years, significant progress has been made in the study of PAH-CHD and omics. To provide a comprehensive description and promote further in-depth investigation of PAH-CHD, this review attempts to summarize the latest developments in genomics, transcriptomics, epigenomics, proteomics, metabolomics, and multi-omics integration.

Immune-related toxicity and soluble profile in patients affected by solid tumors: a network approach

BACKGROUND

Immune checkpoint inhibitors (ICIs) have particular, immune-related adverse events (irAEs), as a consequence of interfering with self-tolerance mechanisms. The incidence of irAEs varies depending on ICI class, administered dose and treatment schedule. The aim of this study was to define a baseline (T0) immune profile (IP) predictive of irAE development.

METHODS

A prospective, multicenter study evaluating the immune profile (IP) of 79 patients with advanced cancer and treated with anti-programmed cell death protein 1 (anti-PD-1) drugs as a first- or second-line setting was performed. The results were then correlated with irAEs onset. The IP was studied by means of multiplex assay, evaluating circulating concentration of 12 cytokines, 5 chemokines, 13 soluble immune checkpoints and 3 adhesion molecules. Indoleamine 2, 3-dioxygenase (IDO) activity was measured through a modified liquid chromatography-tandem mass spectrometry using the high-performance liquid chromatography-mass spectrometry (HPLC-MS/MS) method. A connectivity heatmap was obtained by calculating Spearman correlation coefficients. Two different networks of connectivity were constructed, based on the toxicity profile.

RESULTS

Toxicity was predominantly of low/moderate grade. High-grade irAEs were relatively rare, while cumulative toxicity was high (35%). Positive and statistically significant correlations between the cumulative toxicity and IP10 and IL8, sLAG3, sPD-L2, sHVEM, sCD137, sCD27 and sICAM-1 serum concentration were found. Moreover, patients who experienced irAEs had a markedly different connectivity pattern, characterized by disruption of most of the paired connections between cytokines, chemokines and connections of sCD137, sCD27 and sCD28, while sPDL-2 pair-wise connectivity values seemed to be intensified. Network connectivity analysis identified a total of 187 statistically significant interactions in patients without toxicity and a total of 126 statistically significant interactions in patients with toxicity. Ninety-eight interactions were common to both networks, while 29 were specifically observed in patients who experienced toxicity.

CONCLUSIONS

A particular, common pattern of immune dysregulation was defined in patients developing irAEs. This immune serological profile, if confirmed in a larger patient population, could lead to the design of a personalized therapeutic strategy in order to prevent, monitor and treat irAEs at an early stage.

Normalized L3-based link prediction in protein-protein interaction networks

BACKGROUND

Protein-protein interaction (PPI) data is an important type of data used in functional genomics. However, high-throughput experiments are often insufficient to complete the PPI interactome of different organisms. Computational techniques are thus used to infer missing data, with link prediction being one such approach that uses the structure of the network of PPIs known so far to identify non-edges whose addition to the network would make it more sound, according to some underlying assumptions. Recently, a new idea called the L3 principle introduced biological motivation into PPI link predictions, yielding predictors that are superior to general-purpose link predictors for complex networks. Interestingly, the L3 principle can be interpreted in another way, so that other signatures of PPI networks can also be characterized for PPI predictions. This alternative interpretation uncovers candidate PPIs that the current L3-based link predictors may not be able to fully capture, underutilizing the L3 principle.

RESULTS

In this article, we propose a formulation of link predictors that we call NormalizedL3 (L3N) which addresses certain missing elements within L3 predictors in the perspective of network modeling. Our computational validations show that the L3N predictors are able to find missing PPIs more accurately (in terms of true positives among the predicted PPIs) than the previously proposed methods on several datasets from the literature, including BioGRID, STRING, MINT, and HuRI, at the cost of using more computation time in some of the cases. In addition, we found that L3-based link predictors (including L3N) ranked a different pool of PPIs higher than the general-purpose link predictors did. This suggests that different types of PPIs can be predicted based on different topological assumptions, and that even better PPI link predictors may be obtained in the future by improved network modeling.

Monitoring of Serological, Cellular and Genomic Biomarkers in Transplantation, Computational Prediction Models and Role of Cell-Free DNA in Transplant Outcome

The process and evolution of an organ transplant procedure has evolved in terms of the prevention of immunological rejection with the improvement in the determination of immune response genes. These techniques include considering more important genes, more polymorphism detection, more refinement of the response motifs, as well as the analysis of epitopes and eplets, its capacity to fix complement, the PIRCHE algorithm and post-transplant monitoring with promising new biomarkers that surpass the classic serum markers such as creatine and other similar parameters of renal function. Among these new biomarkers, we analyze new serological, urine, cellular, genomic and transcriptomic biomarkers and computational prediction, with particular attention to the analysis of donor free circulating DNA as an optimal marker of kidney damage.

NIAPU: network-informed adaptive positive-unlabeled learning for disease gene identification

MOTIVATION

Gene-disease associations are fundamental for understanding disease etiology and developing effective interventions and treatments. Identifying genes not yet associated with a disease due to a lack of studies is a challenging task in which prioritization based on prior knowledge is an important element. The computational search for new candidate disease genes may be eased by positive-unlabeled learning, the machine learning (ML) setting in which only a subset of instances are labeled as positive while the rest of the dataset is unlabeled. In this work, we propose a set of effective network-based features to be used in a novel Markov diffusion-based multi-class labeling strategy for putative disease gene discovery.

RESULTS

The performances of the new labeling algorithm and the effectiveness of the proposed features have been tested on 10 different disease datasets using three ML algorithms. The new features have been compared against classical topological and functional/ontological features and a set of network- and biological-derived features already used in gene discovery tasks. The predictive power of the integrated methodology in searching for new disease genes has been found to be competitive against state-of-the-art algorithms.

AVAILABILITY AND IMPLEMENTATION

The source code of NIAPU can be accessed at https://github.com/AndMastro/NIAPU. The source data used in this study are available online on the respective websites.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

"Transplantomics" for predicting allograft rejection: real-life applications and new strategies from Network Medicine

Although decades of the reductionist approach achieved great milestones in optimizing the immunosuppression therapy, traditional clinical parameters still fail in predicting both acute and chronic (mainly) rejection events leading to higher rates across all solid organ transplants. To clarify the underlying immune-related cellular and molecular mechanisms, current biomedical research is increasingly focusing on "transplantomics" which relies on a huge quantity of big data deriving from genomics, transcriptomics, epigenomics, proteomics, and metabolomics platforms. The AlloMap (gene expression) and the AlloSure (donor-derived cell-free DNA) tests represent two successful examples of how omics and liquid biopsy can really improve the precision medicine of heart and kidney transplantation. One of the major challenges in translating big data in clinically useful biomarkers is the integration and interpretation of the different layers of omics datasets. Network Medicine offers advanced bioinformatic-molecular strategies which were widely used to integrate large omics datasets and clinical information in end-stage patients to prioritize potential biomarkers and drug targets. The application of network-oriented approaches to clarify the complex nature of graft rejection is still in its infancy. Here, we briefly discuss the real-life clinical applications derived from omics datasets as well as novel opportunities for establishing predictive tests in solid organ transplantation. Also, we provide an original "graft rejection interactome" and propose network-oriented strategies which can be useful to improve precision medicine of solid organ transplantation.

Multi-omics and machine learning for the prevention and management of female reproductive health

Females typically carry most of the burden of reproduction in mammals. In humans, this burden is exacerbated further, as the evolutionary advantage of a large and complex human brain came at a great cost of women's reproductive health. Pregnancy thus became a highly demanding phase in a woman's life cycle both physically and emotionally and therefore needs monitoring to assure an optimal outcome. Moreover, an increasing societal trend towards reproductive complications partly due to the increasing maternal age and global obesity pandemic demands closer monitoring of female reproductive health. This review first provides an overview of female reproductive biology and further explores utilization of large-scale data analysis and -omics techniques (genomics, transcriptomics, proteomics, and metabolomics) towards diagnosis, prognosis, and management of female reproductive disorders. In addition, we explore machine learning approaches for predictive models towards prevention and management. Furthermore, mobile apps and wearable devices provide a promise of continuous monitoring of health. These complementary technologies can be combined towards monitoring female (fertility-related) health and detection of any early complications to provide intervention solutions. In summary, technological advances (e.g., omics and wearables) have shown a promise towards diagnosis, prognosis, and management of female reproductive disorders. Systematic integration of these technologies is needed urgently in female reproductive healthcare to be further implemented in the national healthcare systems for societal benefit.

Association Between Circulating CD4+ T Cell Methylation Signatures of Network-Oriented SOCS3 Gene and Hemodynamics in Patients Suffering Pulmonary Arterial Hypertension

Pathogenic DNA methylation changes may be involved in pulmonary arterial hypertension (PAH) onset and its progression, but there is no data on potential associations with patient-derived hemodynamic parameters. The reduced representation bisulfite sequencing (RRBS) platform identified N = 631 differentially methylated CpG sites which annotated to N = 408 genes (DMGs) in circulating CD4⁺ T cells isolated from PAH patients vs. healthy controls (CTRLs). A promoter-restricted network analysis established the PAH subnetwork that included 5 hub DMGs (SOCS3, GNAS, ITGAL, NCOR2, NFIC) and 5 non-hub DMGs (NR4A2, GRM2, PGK1, STMN1, LIMS2). The functional analysis revealed that the SOCS3 gene was the most recurrent among the top ten significant pathways enriching the PAH subnetwork, including the growth hormone receptor and the interleukin-6 signaling. Correlation analysis showed that the promoter methylation levels of each network-oriented DMG were associated individually with hemodynamic parameters. In particular, SOCS3 hypomethylation was negatively associated with right atrial pressure (RAP) and positively associated with cardiac index (CI) (|r|≥ 0.6). A significant upregulation of the SOCS3, ITGAL, NFIC, NCOR2, and PGK1 mRNA levels (qRT-PCR) in peripheral blood mononuclear cells from PAH patients vs. CTRLs was found (P ≤ 0.05). By immunoblotting, a significant upregulation of the SOCS3 protein was confirmed in PAH patients vs. CTRLs (P < 0.01). This is the first network-oriented study which integrates circulating CD4⁺ T cell DNA methylation signatures, hemodynamic parameters, and validation experiments in PAH patients at first diagnosis or early follow-up. Our data suggests that SOCS3 gene might be involved in PAH pathogenesis and serve as potential prognostic biomarker.

Semi-automated socio-anthropologic analysis of the medical discourse on rheumatoid arthritis: Potential impact on public health

BACKGROUND

The debilitating effects of noncommunicable diseases (NCDs) and the accompanying chronic inflammation represent a significant obstacle for the sustainability of our development, with efforts spreading worldwide to counteract the diffusion of NCDs, as per the United Nations Sustainable Development Goals (SDG 3). In fact, despite efforts of varied intensity in numerous directions (from innovations in biotechnology to lifestyle modifications), the incidence of NCDs remains pandemic. The present work wants to contribute to addressing this major concern, with a specific focus on the fragmentation of medical approaches, via an interdisciplinary analysis of the medical discourse, i.e. the heterogenous reporting that biomedical scientific literature uses to describe the anti-inflammatory therapeutic landscape in NCDs. The aim is to better capture the roots of this compartmentalization and the power relations existing among three segregated pharmacological, experimental and unstandardized biomedical approaches to ultimately empower collaboration beyond medical specialties and possibly tap into a more ample and effective reservoir of integrated therapeutic opportunities.

METHOD

Using rheumatoid arthritis (RA) as an exemplar disease, twenty-eight articles were manually translated into a nine-dimensional categorical variable of medical socio-anthropological relevance, relating in particular (but not only) to legitimacy, temporality and spatialization. This digitalized picture (9 x 28 table) of the medical discourse was further analyzed by simple automated learning approaches to identify differences and highlight commonalities among the biomedical categories.

RESULTS

Interpretation of these results provides original insights, including suggestions to: empower scientific communication between unstandardized approaches and basic biology; promote the repurposing of non-pharmacological therapies to enhance robustness of experimental approaches; and align the spatial representation of diseases and therapies in pharmacology to effectively embrace the systemic approach promoted by modern personalized and preventive medicines. We hope this original work can expand and foster interdisciplinarity among public health stakeholders, ultimately contributing to the achievement of SDG3.

Network-based approaches for modeling disease regulation and progression

Molecular interaction networks lay the foundation for studying how biological functions are controlled by the complex interplay of genes and proteins. Investigating perturbed processes using biological networks has been instrumental in uncovering mechanisms that underlie complex disease phenotypes. Rapid advances in omics technologies have prompted the generation of high-throughput datasets, enabling large-scale, network-based analyses. Consequently, various modeling techniques, including network enrichment, differential network extraction, and network inference, have proven to be useful for gaining new mechanistic insights. We provide an overview of recent network-based methods and their core ideas to facilitate the discovery of disease modules or candidate mechanisms. Knowledge generated from these computational efforts will benefit biomedical research, especially drug development and precision medicine. We further discuss current challenges and provide perspectives in the field, highlighting the need for more integrative and dynamic network approaches to model disease development and progression.

Network-based methods for psychometric data of eating disorders: A systematic review

BACKGROUND

Network science represents a powerful and increasingly promising method for studying complex real-world problems. In the last decade, it has been applied to psychometric data in the attempt to explain psychopathologies as complex systems of causally interconnected symptoms. One category of mental disorders, relevant for their severity, incidence and multifaceted structure, is that of eating disorders (EDs), serious disturbances that negatively affect a person's eating behavior.

AIMS

We aimed to review the corpus of psychometric network analysis methods by scrutinizing a large sample of network-based studies that exploit psychometric data related to EDs. A particular focus is given to the description of the methodologies for network estimation, network description and network stability analysis providing also a review of the statistical software packages currently used to carry out each phase of the network estimation and analysis workflow. Moreover, we try to highlight aspects with potential clinical impact such as core symptoms, influences of external factors, comorbidities, and related changes in network structure and connectivity across both time and subpopulations.

METHODS

A systematic search was conducted (February 2022) on three different literature databases to identify 57 relevant research articles. The exclusion criteria comprehended studies not based on psychometric data, studies not using network analysis, studies with different aims or not focused on ED, and review articles.

RESULTS

Almost all the selected 57 papers employed the same analytical procedures implemented in a collection of R packages specifically designed for psychometric network analysis and are mostly based on cross-sectional data retrieved from structured psychometric questionnaires, with just few exemptions of panel data. Most of them used the same techniques for all phases of their analysis. In particular, a pervasive use of the Gaussian Graphical Model with LASSO regularization was registered for in network estimation step. Among the clinically relevant results, we can include the fact that all papers found strong symptom interconnections between specific and nonspecific ED symptoms, suggesting that both types should therefore be addressed by clinical treatment.

CONCLUSIONS

We here presented the largest and most comprehensive review to date about psychometric network analysis methods. Although these methods still need solid validation in the clinical setting, they have already been able to show many strengths and important results, as well as great potentials and perspectives, which have been analyzed here to provide suggestions on their use and their possible improvement.

Biological Random Walks: multi-omics integration for disease gene prioritization

MOTIVATION

Over the past decade, network-based approaches have proven useful in identifying disease modules within the human interactome, often providing insights into key mechanisms and guiding the quest for therapeutic targets. This is all the more important, since experimental investigation of potential gene candidates is an expensive task, thus not always a feasible option. On the other hand, many sources of biological information exist beyond the interactome and an important research direction is the design of effective techniques for their integration.

RESULTS

In this work, we introduce the Biological Random Walks (BRW) approach for disease gene prioritization in the human interactome. The proposed framework leverages multiple biological sources within an integrated framework. We perform an extensive, comparative study of BRW's performance against well-established baselines.

AVAILABILITY AND IMPLEMENTATION

All codes are publicly available and can be downloaded at https://github.com/LeoM93/BiologicalRandomWalks. We used publicly available datasets, details on their retrieval and preprocessing are provided in the Supplementary Material.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Network Pharmacology of Adaptogens in the Assessment of Their Pleiotropic Therapeutic Activity

The reductionist concept, based on the ligand–receptor interaction, is not a suitable model for adaptogens, and herbal preparations affect multiple physiological functions, revealing polyvalent pharmacological activities, and are traditionally used in many conditions. This review, for the first time, provides a rationale for the pleiotropic therapeutic efficacy of adaptogens based on evidence from recent gene expression studies in target cells and where the network pharmacology and systems biology approaches were applied. The specific molecular targets and adaptive stress response signaling mechanisms involved in nonspecific modes of action of adaptogens are identified.

Molecular basis of epigenetic regulation in cancer diagnosis and treatment

The global cancer cases and mortality rates are increasing and demand efficient biomarkers for accurate screening, detection, diagnosis, and prognosis. Recent studies have demonstrated that variations in epigenetic mechanisms like aberrant promoter methylation, altered histone modification and mutations in ATP-dependent chromatin remodelling complexes play an important role in the development of carcinogenic events. However, the influence of other epigenetic alterations in various cancers was confirmed with evolving research and the emergence of high throughput technologies. Therefore, alterations in epigenetic marks may have clinical utility as potential biomarkers for early cancer detection and diagnosis. In this review, an outline of the key epigenetic mechanism(s), and their deregulation in cancer etiology have been discussed to decipher the future prospects in cancer therapeutics including precision medicine. Also, this review attempts to highlight the gaps in epigenetic drug development with emphasis on integrative analysis of epigenetic biomarkers to establish minimally non-invasive biomarkers with clinical applications.

Identification of pathogenic genes associated with CKD: An integrated bioinformatics approach

Chronic kidney disease (CKD) is defined as a persistent abnormality in the structure and function of kidneys and leads to high morbidity and mortality in individuals across the world. Globally, approximately 8%–16% of the population is affected by CKD. Proper screening, staging, diagnosis, and the appropriate management of CKD by primary care clinicians are essential in preventing the adverse outcomes associated with CKD worldwide. In light of this, the identification of biomarkers for the appropriate management of CKD is urgently required. Growing evidence has suggested the role of mRNAs and microRNAs in CKD, however, the gene expression profile of CKD is presently uncertain. The present study aimed to identify diagnostic biomarkers and therapeutic targets for patients with CKD. The human microarray profile datasets, consisting of normal samples and treated samples were analyzed thoroughly to unveil the differentially expressed genes (DEGs). After selection, the interrelationship among DEGs was carried out to identify the overlapping DEGs, which were visualized using the Cytoscape program. Furthermore, the PPI network was constructed from the String database using the selected DEGs. Then, from the PPI network, significant modules and sub-networks were extracted by applying the different centralities methods (closeness, betweenness, stress, etc.) using MCODE, Cytohubba, and Centiserver. After sub-network analysis we identified six overlapped hub genes (RPS5, RPL37A, RPLP0, CXCL8, HLA-A, and ANXA1). Additionally, the enrichment analysis was undertaken on hub genes to determine their significant functions. Furthermore, these six genes were used to find their associated miRNAs and targeted drugs. Finally, two genes CXCL8 and HLA-A were common for Ribavirin drug (the gene-drug interaction), after docking studies HLA-A was selected for further investigation. To conclude our findings, we can say that the identified hub genes and their related miRNAs can serve as potential diagnostic biomarkers and therapeutic targets for CKD treatment strategies.

Network Proximity-Based Drug Repurposing Strategy for Early and Late Stages of Primary Biliary Cholangitis

Primary biliary cholangitis (PBC) is a chronic, cholestatic, immune-mediated, and progressive liver disorder. Treatment to preventing the disease from advancing into later and irreversible stages is still an unmet clinical need. Accordingly, we set up a drug repurposing framework to find potential therapeutic agents targeting relevant pathways derived from an expanded pool of genes involved in different stages of PBC. Starting with updated human protein-protein interaction data and genes specifically involved in the early and late stages of PBC, a network medicine approach was used to provide a PBC "proximity" or "involvement" gene ranking using network diffusion algorithms and machine learning models. The top genes in the proximity ranking, when combined with the original PBC-related genes, resulted in a final dataset of the genes most involved in PBC disease. Finally, a drug repurposing strategy was implemented by mining and utilizing dedicated drug-gene interaction and druggable genome information knowledge bases (e.g., the DrugBank repository). We identified several potential drug candidates interacting with PBC pathways after performing an over-representation analysis on our initial 1121-seed gene list and the resulting disease-associated (algorithm-obtained) genes. The mechanism and potential therapeutic applications of such drugs were then thoroughly discussed, with a particular emphasis on different stages of PBC disease. We found that interleukin/EGFR/TNF-alpha inhibitors, branched-chain amino acids, geldanamycin, tauroursodeoxycholic acid, genistein, antioestrogens, curcumin, antineovascularisation agents, enzyme/protease inhibitors, and antirheumatic agents are promising drugs targeting distinct stages of PBC. We developed robust and transparent selection mechanisms for prioritizing already approved medicinal products or investigational products for repurposing based on recognized unmet medical needs in PBC, as well as solid preliminary data to achieve this goal.

Astragalus–Scorpion Drug Pair Inhibits the Development of Prostate Cancer by Regulating GDPD4-2/PI3K/AKT/mTOR Pathway and Autophagy

Objective: Prostate cancer (PCa) is an epithelial malignancy of the prostate that currently lacks effective treatment. Traditional Chinese medicine (TCM) can play an anticancer role through regulating the immune system, anti-tumor angiogenesis, regulating tumor cell apoptosis, autophagy dysfunction, and other mechanisms. This study attempted to explore the active ingredients and potential mechanism of action of the Astragalus–Scorpion (A–S) drug pair in PCa, in order to provide new insights into the treatment of PCa.

Methods: Network pharmacology was used to analyze the A–S drug pair and PCa targets. Bioinformatics analysis was used to analyze the LncRNAs with significant differences in PCa. The expression of LC3 protein was detected by immunofluorescence. CCK8 was used to detect cell proliferation. The expressions of GDPD4-2, AC144450.1, LINC01513, AC004009.2, AL096869.1, AP005210.1, and BX119924.1 were detected by RT-qPCR. The expression of the PI3K/AKT/mTOR pathway and autophagy-related proteins were detected by western blot. LC-MS/MS was used to identify the active components of Astragalus and Scorpion.

Results: A–S drug pair and PCa have a total of 163 targets, which were mainly related to the prostate cancer and PI3K/AKT pathways. A–S drug pair inhibited the formation of PCa, promoted the expression of LC3Ⅱ and Beclin1 proteins, and inhibited the expression of P62 and PI3K–AKT pathway proteins in PCa mice. Astragaloside IV and polypeptide extract from scorpion venom (PESV) were identified as the main active components of the A–S drug pair. GDPD4-2 was involved in the treatment of PCa by Astragaloside IV-PESV. Silencing GDPD4-2 reversed the therapeutic effects of Astragaloside IV-PESV by regulating the PI3K/AKT/mTOR pathway.

Conclusion: Astragaloside IV-PESV is the main active components of A–S drug pair treated PCa by regulating the GDPD4-2/PI3K–AKT/mTOR pathway and autophagy.

Identification of Five Hub Genes Based on Single-Cell RNA Sequencing Data and Network Pharmacology in Patients With Acute Myocardial Infarction

Acute myocardial infarction (AMI) has a high mortality. The single-cell RNA sequencing (scRNA-seq) method was used to analyze disease heterogeneity at the single-cell level. From the Gene Expression Omnibus (GEO) database (GSE180678), AMI scRNA-seq were downloaded and preprocessed by the Seurat package. Gene expression data came from GSE182923. Cell cluster analysis was conducted. Cell types were identified. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analyses were performed on hub genes. Drugs were predicted by protein–protein interaction (PPI) and molecular docking. In total, 7 cell clusters were defined based on the scRNA-seq dataset, and the clusters were labeled as 5 cell types by marker genes. Hematopoietic stem cell types as a differential subgroups were higher in AMI than in healthy tissues. From available databases and PPI analysis, 52 common genets were identified. Based on 52 genes, 5 clusters were obtained using the MCODE algorithm, and genes in these 5 clusters involved in immune and inflammatory pathways were determined. Correlation analysis showed that hematopoietic stem cell types were negatively correlated with ATM, CARM1, and CASP8 but positively correlated with CASP3 and PPARG. This was reversed with immune cells. Molecular docking analysis showed that DB05490 had the lowest docking score with PPARG. We identified 5 hub genes (ATM, CARM1, CASP8, CASP3, and PPARG) involved in AMI progression. Compound DB05490 was a potential inhibitor of PPAG.

Mechanisms of Quercetin against atrial fibrillation explored by network pharmacology combined with molecular docking and experimental validation

Atrial fibrillation (AF) is a common atrial arrhythmia for which there is no specific therapeutic drug. Quercetin (Que) has been used to treat cardiovascular diseases such as arrhythmias. In this study, we explored the mechanism of action of Que in AF using network pharmacology and molecular docking. The chemical structure of Que was obtained from Pubchem. TCMSP, Swiss Target Prediction, Drugbank, STITCH, Pharmmapper, CTD, GeneCards, DISGENET and TTD were used to obtain drug component targets and AF-related genes, and extract AF and normal tissue by GEO database differentially expressed genes by GEO database. The top targets were IL6, VEGFA, JUN, MMP9 and EGFR, and Que for AF treatment might involve the role of AGE-RAGE signaling pathway in diabetic complications, MAPK signaling pathway and IL-17 signaling pathway. Molecular docking showed that Que binds strongly to key targets and is differentially expressed in AF. In vivo results showed that Que significantly reduced the duration of AF fibrillation and improved atrial remodeling, reduced p-MAPK protein expression, and inhibited the progression of AF. Combining network pharmacology and molecular docking approaches with in vivo studies advance our understanding of the intensive mechanisms of Quercetin, and provide the targeted basis for clinical Atrial fibrillation treatment.

A comprehensive review of Artificial Intelligence and Network based approaches to drug repurposing in Covid-19

Conventional drug discovery and development is tedious and time-taking process; because of which it has failed to keep the required pace to mitigate threats and cater demands of viral and re-occurring diseases, such as Covid-19. The main reasons of this delay in traditional drug development are: high attrition rates, extensive time requirements, and huge financial investment with significant risk. The effective solution to de novo drug discovery is drug repurposing. Previous studies have shown that the network-based approaches and analysis are versatile platform for repurposing as the network biology is used to model the interactions between variety of biological concepts. Herein, we provide a comprehensive background of machine learning and deep learning in drug repurposing while specifically focusing on the applications of network-based approach to drug repurposing in Covid-19, data sources, and tools used. Furthermore, use of network proximity, network diffusion, and AI on network-based drug repurposing for Covid-19 is well-explained. Finally, limitations of network-based approaches in general and specific to network are stated along with future recommendations for better network-based models.

Connections for Matters of the Heart: Network Medicine in Cardiovascular Diseases

Cardiovascular diseases (CVD) are diverse disorders affecting the heart and vasculature in millions of people worldwide. Like other fields, CVD research has benefitted from the deluge of multiomics biomedical data. Current CVD research focuses on disease etiologies and mechanisms, identifying disease biomarkers, developing appropriate therapies and drugs, and stratifying patients into correct disease endotypes. Systems biology offers an alternative to traditional reductionist approaches and provides impetus for a comprehensive outlook toward diseases. As a focus area, network medicine specifically aids the translational aspect of in silico research. This review discusses the approach of network medicine and its application to CVD research.

Network Approaches for Charting the Transcriptomic and Epigenetic Landscape of the Developmental Origins of Health and Disease

The early developmental phase is of critical importance for human health and disease later in life. To decipher the molecular mechanisms at play, current biomedical research is increasingly relying on large quantities of diverse omics data. The integration and interpretation of the different datasets pose a critical challenge towards the holistic understanding of the complex biological processes that are involved in early development. In this review, we outline the major transcriptomic and epigenetic processes and the respective datasets that are most relevant for studying the periconceptional period. We cover both basic data processing and analysis steps, as well as more advanced data integration methods. A particular focus is given to network-based methods. Finally, we review the medical applications of such integrative analyses.

Unveiling the effect of Withania somnifera on neuronal cytoarchitecture and synaptogenesis: A combined in vitro and network pharmacology approach

Withania somnifera (WS), is known for its remarkable contribution in herbal medicine and Ayurveda, which is therapeutically applied to improve memory and anxiety in patients. However, the pharmacological details of this plant on memory boosting yet remained undefined. This study provides mechanistic insights on the effect of ethanol solution extract of the whole plant of WS (WSEE) on neuritogenesis by combining in vitro and in silico network pharmacology approaches. WSEE promoted significant neuronal growth through early differentiation, axodendritic arborization, and synaptogenesis on primary hippocampal neurons. The network pharmacological study confirmed that the neuritogenic activity is potentially mediated by modulating the neurotrophin signaling pathway, where NRTK1 (TrkA) was revealed as the primary target of WS secondary metabolites. This neurotrophic activity of WSEE was significantly stifled by the presence of TrkA inhibitor, which further confirms the TrkA-dependent activity of WSEE. In addition, a molecular docking study suggested steroidal lactones present in the WS might act as nerve growth factor (NGF)-mimetics, activating TrkA by binding to the NGF-binding domain. As a whole, the findings of the study suggest a significant role of WSEE on neuritogenesis and its potential to function as a therapeutic agent and in drug designing for the prevention and treatment of memory-related neurological disorders.

A Comparison of Network-Based Methods for Drug Repurposing along with an Application to Human Complex Diseases

Drug repurposing strategy, proposing a therapeutic switching of already approved drugs with known medical indications to new therapeutic purposes, has been considered as an efficient approach to unveil novel drug candidates with new pharmacological activities, significantly reducing the cost and shortening the time of de novo drug discovery. Meaningful computational approaches for drug repurposing exploit the principles of the emerging field of Network Medicine, according to which human diseases can be interpreted as local perturbations of the human interactome network, where the molecular determinants of each disease (disease genes) are not randomly scattered, but co-localized in highly interconnected subnetworks (disease modules), whose perturbation is linked to the pathophenotype manifestation. By interpreting drug effects as local perturbations of the interactome, for a drug to be on-target effective against a specific disease or to cause off-target adverse effects, its targets should be in the nearby of disease-associated genes. Here, we used the network-based proximity measure to compute the distance between the drug module and the disease module in the human interactome by exploiting five different metrics (minimum, maximum, mean, median, mode), with the aim to compare different frameworks for highlighting putative repurposable drugs to treat complex human diseases, including malignant breast and prostate neoplasms, schizophrenia, and liver cirrhosis. Whilst the standard metric (that is the minimum) for the network-based proximity remained a valid tool for efficiently screening off-label drugs, we observed that the other implemented metrics specifically predicted further interesting drug candidates worthy of investigation for yielding a potentially significant clinical benefit.

Integrated Analysis of ceRNA Network to Reveal Potential Prognostic Biomarkers for Glioblastoma

Glioblastoma (GBM), originating in the brain, is a universally aggressive malignant tumor with a particularly poor prognosis. Therefore, insight into the critical role of underlying genetic mechanisms is essential to developing new therapeutic approaches. This study aims to identify potential markers with clinical and prognostic significance in GBM. To this end, increasing numbers of differentially expressed RNA have been identified used to construct competitive endogenous RNA networks for prognostic analysis via comparison and analysis of RNA expression levels of tumor and normal tissues in glioblastoma. This analysis demonstrated that the RNA expression patterns of normal and tumor samples were significantly different. Thus, the resulting differentially expressed RNAs were used to construct competitive endogenous RNA (competing endogenous RNA, ceRNA) networks. The functional enrichment indicated mRNAs in the network are critically involved in a variety of biological functions. Additionally, the prognostic analysis suggested 27 lncRNAs, including LOXL1-AS1, AL356414.1, etc., were significantly associated with patient survival. Given the prognostic significance of these 27 lncRNAs in GBM, we sought to classify the samples. Importantly, Kaplan-Meier analysis revealed that survival times varied significantly among the different categories. Overall, these results identify that the candidate lncRNAs are potential prognostic markers of GBM and its corresponding mRNAs may be a potential target for therapy.