Gene co-expression in the interactome: moving from correlation toward causation via an integrated approach to disease module discovery

Network-based drug repurposing for psychiatric disorders using single-cell genomics

Neuropsychiatric disorders lack effective treatments due to a limited understanding of underlying cellular and molecular mechanisms. To address this, we integrated population-scale single-cell genomics data and analyzed cell-type-level gene regulatory networks across schizophrenia, bipolar disorder, and autism (23 cell classes/subclasses). Our analysis revealed potential druggable transcription factors co-regulating known risk genes that converge into cell-type-specific co-regulated modules. We applied graph neural networks on those modules to prioritize novel risk genes and leveraged them in a network-based drug repurposing framework to identify 220 drug molecules with the potential for targeting specific cell types. We found evidence for 37 of these drugs in reversing disorder-associated transcriptional phenotypes. Additionally, we discovered 335 drug-associated cell-type eQTLs, revealing genetic variation's influence on drug target expression at the cell-type level. Our results provide a single-cell network medicine resource that provides mechanistic insights for advancing treatment options for neuropsychiatric disorders.

Transforming drug discovery: the impact of AI and molecular simulation on R&D efficiency

The process of developing new drugs in the pharmaceutical industry is both time-consuming and costly, making efficiency crucial. Recent advances in hardware and computational methods have led to the widespread application of computational science approaches in drug discovery. These approaches, including artificial intelligence and molecular simulations, span from target identification to pharmacokinetics research, aiming to reduce the likelihood of failure and present lower costs. Machine learning-based methods predict new applications for developing new drugs based on accumulated knowledge, while molecular simulations estimate interactions between drugs and target proteins at the atomic level based on physical laws. Each approach has its advantages and disadvantages, and they complement each other. As a result, the future of computational science approaches in drug discovery is expected to focus on developing new methodologies that integrate these two techniques to enhance the efficiency of drug discovery.

Network medicine based approach for identifying the type 2 diabetes, osteoarthritis and triple negative breast cancer interactome: Finding the hub of hub genes

The increasing prevalence of multi-morbidities, particularly the incidence of breast cancer in diabetic/osteoarthritic patients emphasize on the need for exploring the underlying molecular mechanisms resulting in carcinogenesis. To address this, present study employed a systems biology approach to identify switch genes pivotal to the crosstalk between diseased states resulting in multi-morbid conditions. Hub genes previously reported for type 2 diabetes mellitus (T2DM), osteoarthritis (OA), and triple negative breast cancer (TNBC), were extracted from published literature and fed into an integrated bioinformatics analyses pipeline. Thirty-one hub genes common to all three diseases were identified. Functional enrichment analyses showed these were mainly enriched for immune and metabolism associated terms including advanced glycation end products (AGE) pathways, cancer pathways, particularly breast neoplasm, immune system signalling and adipose tissue. The T2DM-OA-TNBC interactome was subjected to protein-protein interaction network analyses to identify meta hub/clustered genes. These were prioritized and wired into a three disease signalling map presenting the enriched molecular crosstalk on T2DM-OA-TNBC axes to gain insight into the molecular mechanisms underlying disease-disease interactions. Deciphering the molecular bases for the intertwined metabolic and immune states may potentiate the discovery of biomarkers critical for identifying and targeting the immuno-metabolic origin of disease.

ECM Microenvironment in Vascular Homeostasis: New Targets for Atherosclerosis

Alterations in vascular extracellular matrix (ECM) components, interactions, and mechanical properties influence both the formation and stability of atherosclerotic plaques. This review discusses the contribution of the ECM microenvironment in vascular homeostasis and remodeling in atherosclerosis, highlighting Cartilage oligomeric matrix protein (COMP) and its degrading enzyme ADAMTS7 as examples, and proposes potential avenues for future research aimed at identifying novel therapeutic targets for atherosclerosis based on the ECM microenvironment.

Immune cell landscapes are associated with high-grade serous ovarian cancer survival

High-grade serous ovarian cancer (HGSOC) is an aggressive disease known to develop resistance to chemotherapy. We investigated the prognostic significance of tumor cell states and potential mechanisms underlying chemotherapy resistance in HGSOC. Transcriptome deconvolution was performed to address cellular heterogeneity. Kaplan-Meier survival curves were plotted to illustrate the outcomes of patients with varying cellular abundances. The association between gene expression and chemotherapy response was tested. After adjusting for surgery status and grading, several cell states exhibited a significant correlation with patient survival. Cell states can organize into carcinoma ecotypes (CE). CE9 and CE10 were proinflammatory, characterized by higher immunoreactivity, and were associated with favorable survival outcomes. Ratios of cell states and ecotypes had better prognostic abilities than a single cell state or ecotype. A total of 1265 differentially expressed genes were identified between samples with high and low levels of C9 or CE10. These genes were partitioned into three co-expressed modules, which were associated with tumor cells and immune cells. Pogz was identified to be linked with immune cell genes and the chemotherapy response of paclitaxel. Collectively, the survival of HGSOC patients is correlated with specific cell states and ecotypes.

scRank infers drug-responsive cell types from untreated scRNA-seq data using a target-perturbed gene regulatory network

Cells respond divergently to drugs due to the heterogeneity among cell populations. Thus, it is crucial to identify drug-responsive cell populations in order to accurately elucidate the mechanism of drug action, which is still a great challenge. Here, we address this problem with scRank, which employs a target-perturbed gene regulatory network to rank drug-responsive cell populations via in silico drug perturbations using untreated single-cell transcriptomic data. We benchmark scRank on simulated and real datasets, which shows the superior performance of scRank over existing methods. When applied to medulloblastoma and major depressive disorder datasets, scRank identifies drug-responsive cell types that are consistent with the literature. Moreover, scRank accurately uncovers the macrophage subpopulation responsive to tanshinone IIA and its potential targets in myocardial infarction, with experimental validation. In conclusion, scRank enables the inference of drug-responsive cell types using untreated single-cell data, thus providing insights into the cellular-level impacts of therapeutic interventions.

Exploration of functional relations among differentially co-expressed genes identifies regulators in glioblastoma

The conventional computational approaches to investigating a disease confront inherent constraints as they often need to improve in delving beyond protein functional associations and grasping their deeper contextual significance within the disease framework. Such context-specificity can be explored using clinical data by evaluating the change in interaction between the biological entities in different conditions by investigating the differential co-expression relationships. We believe that the integration and analysis of differential co-expression and the functional relationships, primarily focusing on the source nodes, will open novel insights about disease progression as the source proteins could trigger signaling cascades, mostly because they are transcription factors, cell surface receptors, or enzymes that respond instantly to a particular stimulus. A thorough contextual investigation of these nodes could lead to a helpful beginning point for identifying potential causal linkages and guiding subsequent scientific investigations to uncover mechanisms underlying observed associations. Our methodology includes functional protein-protein Interaction (PPI) data and co-expression information and filters functional linkages through a series of critical steps, culminating in the identification of a robust set of regulators. Our analysis identified eleven key regulators-AKT1, BRCA1, CAMK2G, CUL1, FGFR3, KIF3A, NUP210, PRKACB, RAB8A, RPS6KA2 and TGFB3-in glioblastoma. These regulators play a pivotal role in disease classification, cell growth control, and patient survivability and exhibit associations with immune infiltrations and disease hallmarks. This underscores the importance of assessing correlation towards causality in unraveling complex biological insights.

In-depth analysis of immune cell landscapes reveals differences between lung adenocarcinoma and lung squamous cell carcinoma

Background

Lung cancer is the leading cause of cancer deaths globally, with lung adenocarcinoma (LUAD) and squamous cell carcinoma (LUSC) being major subtypes. Immunotherapy has emerged as a promising approach for the treatment of lung cancer, but understanding the underlying mechanisms of immune dysregulation is crucial for the development of effective therapies. This study aimed to investigate the distinctive cellular features of LUAD and LUSC and identify potential biomarkers associated with the pathogenesis and clinical outcomes of each subtype.

Methods

We used digital cytometry techniques to analyze the RNA-Seq data of 1128 lung cancer patients from The Cancer Genome Atlas (TCGA) database. The abundance of cell subtypes and ecotypes in LUAD and LUSC patients was quantified. Univariate survival analysis was used to investigate their associations with patient overall survival (OS). Differential gene expression analysis and gene co-expression network construction were carried out to explore the gene expression patterns of LUSC patients with distinct survival outcomes. Scratch wound-healing assay, colony formation assay, and transwell assay were used to validate the candidate drugs for LUSC treatment.

Results

We found differential expression of cell subtypes between LUAD and LUSC, with certain cell subtypes being prognostic for survival in both subtypes. We also identified differential gene expression and gene co-expression modules associated with macrophages.3/PCs.2 ratio in LUSC patients with distinct survival outcomes. Furthermore, ecotype ratios were found to be prognostic in both subtypes and machine learning models showed that certain cell subtypes, such as epithelial.cells.1, epithelial.cells.5, and endothelial.cells.2 are important for predicting LUSC. Ginkgolide B and triamterene can inhibit the proliferation, invasion, and migration of LUSC cell lines.

Conclusion

We provide insight into the distinctive cellular features of LUAD and LUSC, and identify potential biomarkers associated with the pathogenesis and clinical outcomes of each subtype. Ginkgolide B and triamterene could be promising drugs for LUSC treatment.

TTD: Therapeutic Target Database describing target druggability information

Target discovery is one of the essential steps in modern drug development, and the identification of promising targets is fundamental for developing first-in-class drug. A variety of methods have emerged for target assessment based on druggability analysis, which refers to the likelihood of a target being effectively modulated by drug-like agents. In the therapeutic target database (TTD), nine categories of established druggability characteristics were thus collected for 426 successful, 1014 clinical trial, 212 preclinical/patented, and 1479 literature-reported targets via systematic review. These characteristic categories were classified into three distinct perspectives: molecular interaction/regulation, human system profile and cell-based expression variation. With the rapid progression of technology and concerted effort in drug discovery, TTD and other databases were highly expected to facilitate the explorations of druggability characteristics for the discovery and validation of innovative drug target. TTD is now freely accessible at: https://idrblab.org/ttd/.

Cross-linking breast tumor transcriptomic states and tissue histology

Identification of the gene expression state of a cancer patient from routine pathology imaging and characterization of its phenotypic effects have significant clinical and therapeutic implications. However, prediction of expression of individual genes from whole slide images (WSIs) is challenging due to co-dependent or correlated expression of multiple genes. Here, we use a purely data-driven approach to first identify groups of genes with co-dependent expression and then predict their status from WSIs using a bespoke graph neural network. These gene groups allow us to capture the gene expression state of a patient with a small number of binary variables that are biologically meaningful and carry histopathological insights for clinical and therapeutic use cases. Prediction of gene expression state based on these gene groups allows associating histological phenotypes (cellular composition, mitotic counts, grading, etc.) with underlying gene expression patterns and opens avenues for gaining biological insights from routine pathology imaging directly.

Regulation of expression quantitative trait loci by SVA retrotransposons within the major histocompatibility complex

Genomic and transcriptomic studies of expression quantitative trait loci (eQTL) revealed that SINE-VNTR-Alu (SVA) retrotransposon insertion polymorphisms (RIPs) within human genomes markedly affect the co-expression of many coding and noncoding genes by coordinated regulatory processes. This study examined the polymorphic SVA modulation of gene co-expression within the major histocompatibility complex (MHC) genomic region where more than 160 coding genes are involved in innate and adaptive immunity. We characterized the modulation of SVA RIPs utilizing the genomic and transcriptomic sequencing data obtained from whole blood of 1266 individuals in the Parkinson's Progression Markers Initiative (PPMI) cohort that included an analysis of human leukocyte antigen (HLA)-A regulation in a subpopulation of the cohort. The regulatory properties of eight SVAs located within the class I and class II MHC regions were associated with differential co-expression of 71 different genes within and 75 genes outside the MHC region. Some of the same genes were affected by two or more different SVA. Five SVA are annotated in the human genomic reference sequence GRCh38.p14/hg38, whereas the other three were novel insertions within individuals. We also examined and found distinct structural effects (long and short variants and the CT internal variants) for one of the SVA (R_SVA_24) insertions on the differential expression of the HLA-A gene within a subpopulation (550 individuals) of the PPMI cohort. This is the first time that many HLA and non-HLA genes (multilocus expression units) and splicing mechanisms have been shown to be regulated by eight structurally polymorphic SVA within the MHC genomic region by applying precise statistical analysis of RNA data derived from the blood samples of a human cohort population. This study shows that SVA within the MHC region are important regulators or rheostats of gene co-expression that might have potential roles in diversity, health, and disease.

Navigating the Gene Co-Expression Network and Drug Repurposing Opportunities for Brain Disorders Associated with Neurocognitive Impairment

Neurocognitive impairment refers to a spectrum of disorders characterized by a decline in cognitive functions such as memory, attention, and problem-solving, which are often linked to structural or functional abnormalities in the brain. While its exact etiology remains elusive, genetic factors play a pivotal role in disease onset and progression. This study aimed to identify highly correlated gene clusters (modules) and key hub genes shared across neurocognition-impairing diseases, including Alzheimer's disease (AD), Parkinson's disease with dementia (PDD), HIV-associated neurocognitive disorders (HAND), and glioma. Herein, the microarray datasets AD (GSE5281), HAND (GSE35864), glioma (GSE15824), and PD (GSE7621) were used to perform Weighted Gene Co-expression Network Analysis (WGCNA) to identify highly preserved modules across the studied brain diseases. Through gene set enrichment analysis, the shared modules were found to point towards processes including neuronal transcriptional dysregulation, neuroinflammation, protein aggregation, and mitochondrial dysfunction, hallmarks of many neurocognitive disorders. These modules were used in constructing protein-protein interaction networks to identify hub genes shared across the diseases of interest. These hub genes were found to play pivotal roles in processes including protein homeostasis, cell cycle regulation, energy metabolism, and signaling, all associated with brain and CNS diseases, and were explored for their drug repurposing experiments. Drug repurposing based on gene signatures highlighted drugs including Dorzolamide and Oxybuprocaine, which were found to modulate the expression of the hub genes in play and may have therapeutic implications in neurocognitive disorders. While both drugs have traditionally been used for other medical purposes, our study underscores the potential of a combined WGCNA and drug repurposing strategy for searching for new avenues in the simultaneous treatment of different diseases that have similarities in gene co-expression networks.

Ecological network analysis reveals cancer-dependent chaperone-client interaction structure and robustness

Cancer cells alter the expression levels of metabolic enzymes to fuel proliferation. The mitochondrion is a central hub of metabolic reprogramming, where chaperones service hundreds of clients, forming chaperone-client interaction networks. How network structure affects its robustness to chaperone targeting is key to developing cancer-specific drug therapy. However, few studies have assessed how structure and robustness vary across different cancer tissues. Here, using ecological network analysis, we reveal a non-random, hierarchical pattern whereby the cancer type modulates the chaperones' ability to realize their potential client interactions. Despite the low similarity between the chaperone-client interaction networks, we highly accurately predict links in one cancer type based on another. Moreover, we identify groups of chaperones that interact with similar clients. Simulations of network robustness show that this group structure affects cancer-specific response to chaperone removal. Our results open the door for new hypotheses regarding the ecology and evolution of chaperone-client interaction networks and can inform cancer-specific drug development strategies.

Circadian Regulation of the Neuroimmune Environment Across the Lifespan: From Brain Development to Aging

Circadian clocks confer 24-h periodicity to biological systems, to ultimately maximize energy efficiency and promote survival in a world with regular environmental light cycles. In mammals, circadian rhythms regulate myriad physiological functions, including the immune, endocrine, and central nervous systems. Within the central nervous system, specialized glial cells such as astrocytes and microglia survey and maintain the neuroimmune environment. The contributions of these neuroimmune cells to both homeostatic and pathogenic demands vary greatly across the day. Moreover, the function of these cells changes across the lifespan. In this review, we discuss circadian regulation of the neuroimmune environment across the lifespan, with a focus on microglia and astrocytes. Circadian rhythms emerge in early life concurrent with neuroimmune sculpting of brain circuits and wane late in life alongside increasing immunosenescence and neurodegeneration. Importantly, circadian dysregulation can alter immune function, which may contribute to susceptibility to neurodevelopmental and neurodegenerative diseases. In this review, we highlight circadian neuroimmune interactions across the lifespan and share evidence that circadian dysregulation within the neuroimmune system may be a critical component in human neurodevelopmental and neurodegenerative diseases.

Persistent homology analysis of type 2 diabetes genome-wide association studies in protein-protein interaction networks

Genome-wide association studies (GWAS) involving increasing sample sizes have identified hundreds of genetic variants associated with complex diseases, such as type 2 diabetes (T2D); however, it is unclear how GWAS hits form unique topological structures in protein-protein interaction (PPI) networks. Using persistent homology, this study explores the evolution and persistence of the topological features of T2D GWAS hits in the PPI network with increasing p-value thresholds. We define an n-dimensional persistent disease module as a higher-order generalization of the largest connected component (LCC). The 0-dimensional persistent T2D disease module is the LCC of the T2D GWAS hits, which is significantly detected in the PPI network (196 nodes and 235 edges, P< 0.05). In the 1-dimensional homology group analysis, all 18 1-dimensional holes (loops) of the T2D GWAS hits persist over all p-value thresholds. The 1-dimensional persistent T2D disease module comprising these 18 persistent 1-dimensional holes is significantly larger than that expected by chance (59 nodes and 83 edges, P< 0.001), indicating a significant topological structure in the PPI network. Our computational topology framework potentially possesses broad applicability to other complex phenotypes in identifying topological features that play an important role in disease pathobiology.

R-Mixup: Riemannian Mixup for Biological Networks

Biological networks are commonly used in biomedical and healthcare domains to effectively model the structure of complex biological systems with interactions linking biological entities. However, due to their characteristics of high dimensionality and low sample size, directly applying deep learning models on biological networks usually faces severe overfitting. In this work, we propose R-Mixup, a Mixup-based data augmentation technique that suits the symmetric positive definite (SPD) property of adjacency matrices from biological networks with optimized training efficiency. The interpolation process in R-Mixup leverages the log-Euclidean distance metrics from the Riemannian manifold, effectively addressing the swelling effect and arbitrarily incorrect label issues of vanilla Mixup. We demonstrate the effectiveness of R-Mixup with five real-world biological network datasets on both regression and classification tasks. Besides, we derive a commonly ignored necessary condition for identifying the SPD matrices of biological networks and empirically study its influence on the model performance. The code implementation can be found in Appendix E.

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.

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.

Globally invariant behavior of oncogenes and random genes at population but not at single cell level

Cancer is widely considered a genetic disease. Notably, recent works have highlighted that every human gene may possibly be associated with cancer. Thus, the distinction between genes that drive oncogenesis and those that are associated to the disease, but do not play a role, requires attention. Here we investigated single cells and bulk (cell-population) datasets of several cancer transcriptomes and proteomes in relation to their healthy counterparts. When analyzed by machine learning and statistical approaches in bulk datasets, both general and cancer-specific oncogenes, as defined by the Cancer Genes Census, show invariant behavior to randomly selected gene sets of the same size for all cancers. However, when protein-protein interaction analyses were performed, the oncogenes-derived networks show higher connectivity than those relative to random genes. Moreover, at single-cell scale, we observe variant behavior in a subset of oncogenes for each considered cancer type. Moving forward, we concur that the role of oncogenes needs to be further scrutinized by adopting protein causality and higher-resolution single-cell analyses.

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.

Landscape of co-expressed genes between the myocardium and blood in sepsis and ceRNA network construction: a bioinformatic approach

Septic cardiomyopathy is a serious complication of sepsis. The mechanism of disease pathogenesis, which is caused by infection, is well researched. Despite ongoing efforts, there are no viable biological markers in the peripheral blood for early detection and diagnosis of septic cardiomyopathy. We aimed to uncover potential biomarkers of septic cardiomyopathy by comparing the covaried genes and pathways in the blood and myocardium of sepsis patients. Gene expression profiling of GSE79962, GSE65682, GSE54514, and GSE134364 was retrieved from the GEO database. Student's t-test was used for differential expression analysis. K-means clustering analysis was applied for subgroup identification. Least absolute shrinkage and selection operator (LASSO) and logistic regression were utilized for screening characteristic genes and model construction. Receiver operating characteristic (ROC) curves were generated for estimating the diagnostic efficacy. For ceRNA information prediction, miWalk and lncBase were applied. Cytoscape was used for ceRNA network construction. Inflammation-associated genes were upregulated, while genes related to mitochondria and aerobic metabolism were downregulated in both blood and the myocardium. Three groups with a significantly different mortality were identified by these covaried genes, using clustering analysis. Five characteristic genes-BCL2A1, CD44, ADGRG1, TGIF1, and ING3-were identified, which enabled the prediction of mortality of sepsis. The pathophysiological changes in the myocardium of patients with sepsis were also reflected in peripheral blood to some extent. The co-occurring pathological processes can affect the prognosis of sepsis. Thus, the genes we identified have the potential to become biomarkers for septic cardiomyopathy.

Exploring novel disease-disease associations based on multi-view fusion network

Established taxonomy system based on disease symptom and tissue characteristics have provided an important basis for physicians to correctly identify diseases and treat them successfully. However, these classifications tend to be based on phenotypic observations, lacking a molecular biological foundation. Therefore, there is an urgent to integrate multi-dimensional molecular biological information or multi-omics data to redefine disease classification in order to provide a powerful perspective for understanding the molecular structure of diseases. Therefore, we offer a flexible disease classification that integrates the biological process, gene expression, and symptom phenotype of diseases, and propose a disease-disease association network based on multi-view fusion. We applied the fusion approach to 223 diseases and divided them into 24 disease clusters. The contribution of internal and external edges of disease clusters were analyzed. The results of the fusion model were compared with Medical Subject Headings, a traditional and commonly used disease taxonomy. Then, experimental results of model performance comparison show that our approach performs better than other integration methods. As it was observed, the obtained clusters provided more interesting and novel disease-disease associations. This multi-view human disease association network describes relationships between diseases based on multiple molecular levels, thus breaking through the limitation of the disease classification system based on tissues and organs. This approach which motivates clinicians and researchers to reposition the understanding of diseases and explore diagnosis and therapy strategies, extends the existing disease taxonomy.

AVAILABILITY OF DATA AND MATERIALS

The preprocessed dataset and source code supporting the conclusions of this article are available at GitHub repository https://github.com/yangxiaoxi89/mvHDN.

Networks as Biomarkers: Uses and Purposes

Networks-based approaches are often used to analyze gene expression data or protein-protein interactions but are not usually applied to study the relationships between different biomarkers. Given the clinical need for more comprehensive and integrative biomarkers that can help to identify personalized therapies, the integration of biomarkers of different natures is an emerging trend in the literature. Network analysis can be used to analyze the relationships between different features of a disease; nodes can be disease-related phenotypes, gene expression, mutational events, protein quantification, imaging-derived features and more. Since different biomarkers can exert causal effects between them, describing such interrelationships can be used to better understand the underlying mechanisms of complex diseases. Networks as biomarkers are not yet commonly used, despite being proven to lead to interesting results. Here, we discuss in which ways they have been used to provide novel insights into disease susceptibility, disease development and severity.

A Network of MicroRNAs and mRNAs Involved in Melanosome Maturation and Trafficking Defines the Lower Response of Pigmentable Melanoma Cells to Targeted Therapy

BACKGROUND

The ability to increase their degree of pigmentation is an adaptive response that confers pigmentable melanoma cells higher resistance to BRAF inhibitors (BRAFi) compared to non-pigmentable melanoma cells.

METHODS

Here, we compared the miRNome and the transcriptome profile of pigmentable 501Mel and SK-Mel-5 melanoma cells vs. non-pigmentable A375 melanoma cells, following treatment with the BRAFi vemurafenib (vem). In depth bioinformatic analyses (clusterProfiler, WGCNA and SWIMmeR) allowed us to identify the miRNAs, mRNAs and biological processes (BPs) that specifically characterize the response of pigmentable melanoma cells to the drug. Such BPs were studied using appropriate assays in vitro and in vivo (xenograft in zebrafish embryos).

RESULTS

Upon vem treatment, miR-192-5p, miR-211-5p, miR-374a-5p, miR-486-5p, miR-582-5p, miR-1260a and miR-7977, as well as GPR143, OCA2, RAB27A, RAB32 and TYRP1 mRNAs, are differentially expressed only in pigmentable cells. These miRNAs and mRNAs belong to BPs related to pigmentation, specifically melanosome maturation and trafficking. In fact, an increase in the number of intracellular melanosomes-due to increased maturation and/or trafficking-confers resistance to vem.

CONCLUSION

We demonstrated that the ability of pigmentable cells to increase the number of intracellular melanosomes fully accounts for their higher resistance to vem compared to non-pigmentable cells. In addition, we identified a network of miRNAs and mRNAs that are involved in melanosome maturation and/or trafficking. Finally, we provide the rationale for testing BRAFi in combination with inhibitors of these biological processes, so that pigmentable melanoma cells can be turned into more sensitive non-pigmentable cells.

Molecular interaction networks and drug development: Novel approach to drug target identification and drug repositioning

Conventional drug discovery requires identifying a protein target believed to be important for disease mechanism and screening compounds for those that beneficially alter the target's function. While this approach has been an effective one for decades, recent data suggest that its continued success is limited largely owing to the highly prevalent irreducibility of biologically complex systems that govern disease phenotype to a single primary disease driver. Network medicine, a new discipline that applies network science and systems biology to the analysis of complex biological systems and disease, offers a novel approach to overcoming these limitations of conventional drug discovery. Using the comprehensive protein-protein interaction network (interactome) as the template through which subnetworks that govern specific diseases are identified, potential disease drivers are unveiled and the effect of novel or repurposed drugs, used alone or in combination, is studied. This approach to drug discovery offers new and exciting unbiased possibilities for advancing our knowledge of disease mechanisms and precision therapeutics.

Gene Differential Co-Expression Networks Based on RNA-Seq: Construction and Its Applications

Gene co-expression network (GCN) becomes an increasingly important tool in omics data analysis. A great challenge for GCN construction is that the sample size is far lower than the number of genes. Traditional methods rely on considerable samples. Moreover, association signals are likely weak, nonlinear and stochastic, which are difficult to be identified among thousands of candidates. In this paper, the gray correlation coefficient (GCC) is introduced, and a novel method to construct gene differential co-expression networks (GDCNs) is proposed. Based on the GDCNs, three measures are proposed to explore informative genes. The proposed method can make full use of the information provided by a handful of samples and overcome the shortages of GCNs, which can evaluate the changes of co-expression relationships that are possibly triggered by treatments. Based on RNA-seq data of Brassica napus, GDCNs under multiple experimental conditions are constructed and investigated. It is found that the GCC-based method is very robust to data processing. The GDCNs facilitate the inference of gene functions and the identification of informative genes that are responsible for stress responsiveness. The GDCN-based approaches integrate the 'guilt by association' and the 'guilt by rewiring' rules, which provide alternative tools for omics data analysis.

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.

Repurposing Histaminergic Drugs in Multiple Sclerosis

Multiple sclerosis is an autoimmune disease with a strong neuroinflammatory component that contributes to severe demyelination, neurodegeneration and lesions formation in white and grey matter of the spinal cord and brain. Increasing attention is being paid to the signaling of the biogenic amine histamine in the context of several pathological conditions. In multiple sclerosis, histamine regulates the differentiation of oligodendrocyte precursors, reduces demyelination, and improves the remyelination process. However, the concomitant activation of histamine H1–H4 receptors can sustain either damaging or favorable effects, depending on the specifically activated receptor subtype/s, the timing of receptor engagement, and the central versus peripheral target district. Conventional drug development has failed so far to identify curative drugs for multiple sclerosis, thus causing a severe delay in therapeutic options available to patients. In this perspective, drug repurposing offers an exciting and complementary alternative for rapidly approving some medicines already approved for other indications. In the present work, we have adopted a new network-medicine-based algorithm for drug repurposing called SAveRUNNER, for quantifying the interplay between multiple sclerosis-associated genes and drug targets in the human interactome. We have identified new histamine drug-disease associations and predicted off-label novel use of the histaminergic drugs amodiaquine, rupatadine, and diphenhydramine among others, for multiple sclerosis. Our work suggests that selected histamine-related molecules might get to the root causes of multiple sclerosis and emerge as new potential therapeutic strategies for the disease.

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.

Protein-protein interaction networks as miners of biological discovery

Protein-protein interactions (PPIs) form the basis of a myriad of biological pathways and mechanism, such as the formation of protein-complexes or the components of signaling cascades. Here, we reviewed experimental methods for identifying PPI pairs, including yeast two-hybrid, mass spectrometry, co-localization, and co-immunoprecipitation. Furthermore, a range of computational methods leveraging biochemical properties, evolution history, protein structures and more have enabled identification of additional PPIs. Given the wealth of known PPIs, we reviewed important network methods to construct and analyze networks of PPIs. These methods aid biological discovery through identifying hub genes and dynamic changes in the network, and have been thoroughly applied in various fields of biological research. Lastly, we discussed the challenges and future direction of research utilizing the power of PPI networks. This article is protected by copyright. All rights reserved.

Artificial intelligence in cancer target identification and drug discovery

Artificial intelligence is an advanced method to identify novel anticancer targets and discover novel drugs from biology networks because the networks can effectively preserve and quantify the interaction between components of cell systems underlying human diseases such as cancer. Here, we review and discuss how to employ artificial intelligence approaches to identify novel anticancer targets and discover drugs. First, we describe the scope of artificial intelligence biology analysis for novel anticancer target investigations. Second, we review and discuss the basic principles and theory of commonly used network-based and machine learning-based artificial intelligence algorithms. Finally, we showcase the applications of artificial intelligence approaches in cancer target identification and drug discovery. Taken together, the artificial intelligence models have provided us with a quantitative framework to study the relationship between network characteristics and cancer, thereby leading to the identification of potential anticancer targets and the discovery of novel drug candidates.

Comprehensive network medicine-based drug repositioning via integration of therapeutic efficacy and side effects

Despite advances in modern medicine that led to improvements in cardiovascular outcomes, cardiovascular disease (CVD) remains the leading cause of mortality and morbidity globally. Thus, there is an urgent need for new approaches to improve CVD drug treatments. As the development time and cost of drug discovery to clinical application are excessive, alternate strategies for drug development are warranted. Among these are included computational approaches based on omics data for drug repositioning, which have attracted increasing attention. In this work, we developed an adjusted similarity measure implemented by the algorithm SAveRUNNER to reposition drugs for cardiovascular diseases while, at the same time, considering the side effects of drug candidates. We analyzed nine cardiovascular disorders and two side effects. We formulated both disease disorders and side effects as network modules in the human interactome, and considered those drug candidates that are proximal to disease modules but far from side-effects modules as ideal. Our method provides a list of drug candidates for cardiovascular diseases that are unlikely to produce common, adverse side-effects. This approach incorporating side effects is applicable to other diseases, as well.

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.

GCEN: An Easy-to-Use Toolkit for Gene Co-Expression Network Analysis and lncRNAs Annotation

Gene co-expression network analysis has been widely used in gene function annotation, especially for long noncoding RNAs (lncRNAs). However, there is a lack of effective cross-platform analysis tools. For biologists to easily build a gene co-expression network and to predict gene function, we developed GCEN, a cross-platform command-line toolkit developed with C++. It is an efficient and easy-to-use solution that will allow everyone to perform gene co-expression network analysis without the requirement of sophisticated programming skills, especially in cases of RNA-Seq research and lncRNAs function annotation. Because of its modular design, GCEN can be easily integrated into other pipelines.

scTenifoldKnk: An efficient virtual knockout tool for gene function predictions via single-cell gene regulatory network perturbation

Gene knockout (KO) experiments are a proven, powerful approach for studying gene function. However, systematic KO experiments targeting a large number of genes are usually prohibitive due to the limit of experimental and animal resources. Here, we present scTenifoldKnk, an efficient virtual KO tool that enables systematic KO investigation of gene function using data from single-cell RNA sequencing (scRNA-seq). In scTenifoldKnk analysis, a gene regulatory network (GRN) is first constructed from scRNA-seq data of wild-type samples, and a target gene is then virtually deleted from the constructed GRN. Manifold alignment is used to align the resulting reduced GRN to the original GRN to identify differentially regulated genes, which are used to infer target gene functions in analyzed cells. We demonstrate that the scTenifoldKnk-based virtual KO analysis recapitulates the main findings of real-animal KO experiments and recovers the expected functions of genes in relevant cell types.

Network analysis of atherosclerotic genes elucidates druggable targets

Background

Atherosclerosis is one of the major causes of cardiovascular disease. It is characterized by the accumulation of atherosclerotic plaque in arteries under the influence of inflammatory responses, proliferation of smooth muscle cell, accumulation of modified low density lipoprotein. The pathophysiology of atherosclerosis involves the interplay of a number of genes and metabolic pathways. In traditional translation method, only a limited number of genes and pathways can be studied at once. However, the new paradigm of network medicine can be explored to study the interaction of a large array of genes and their functional partners and their connections with the concerned disease pathogenesis. Thus, in our study we employed a branch of network medicine, gene network analysis as a tool to identify the most crucial genes and the miRNAs that regulate these genes at the post transcriptional level responsible for pathogenesis of atherosclerosis.

Result

From NCBI database 988 atherosclerotic genes were retrieved. The protein–protein interaction using STRING database resulted in 22,693 PPI interactions among 872 nodes (genes) at different confidence score. The cluster analysis of the 872 genes using MCODE, a plug-in of Cytoscape software revealed a total of 18 clusters, the topological parameter and gene ontology analysis facilitated in the selection of four influential genes viz., AGT, LPL, ITGB2, IRS1 from cluster 3. Further, the miRNAs (miR-26, miR-27, and miR-29 families) targeting these genes were obtained by employing MIENTURNET webtool.

Conclusion

Gene network analysis assisted in filtering out the 4 probable influential genes and 3 miRNA families in the pathogenesis of atherosclerosis. These genes, miRNAs can be targeted to restrict the occurrence of atherosclerosis. Given the importance of atherosclerosis, any approach in the understanding the genes involved in its pathogenesis can substantially enhance the health care system.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01195-y.

In silico recognition of a prognostic signature in basal-like breast cancer patients

Background

Triple-negative breast cancers (TNBCs) display poor prognosis, have a high risk of tumour recurrence, and exhibit high resistance to drug treatments. Based on their gene expression profiles, the majority of TNBCs are classified as basal-like breast cancers. Currently, there are not available widely-accepted prognostic markers to predict outcomes in basal-like subtype, so the selection of new prognostic indicators for this BC phenotype represents an unmet clinical challenge.

Results

Here, we attempted to address this challenging issue by exploiting a bioinformatics pipeline able to integrate transcriptomic, genomic, epigenomic, and clinical data freely accessible from public repositories. This pipeline starts from the application of the well-established network-based SWIM methodology on the transcriptomic data to unveil important (switch) genes in relation with a complex disease of interest. Then, survival and linear regression analyses are performed to associate the gene expression profiles of the switch genes with both the patients’ clinical outcome and the disease aggressiveness. This allows us to identify a prognostic gene signature that in turn is fed to the last step of the pipeline consisting of an analysis at DNA level, to investigate whether variations in the expression of identified prognostic switch genes could be related to genetic (copy number variations) or epigenetic (DNA methylation differences) alterations in their gene loci, or to the activities of transcription factors binding to their promoter regions. Finally, changes in the protein expression levels corresponding to the so far identified prognostic switch genes are evaluated by immunohistochemical staining results taking advantage of the Human Protein Atlas.

Conclusion

The application of the proposed pipeline on the dataset of The Cancer Genome Atlas (TCGA)-Breast Invasive Carcinoma (BRCA) patients affected by basal-like subtype led to an in silico recognition of a basal-like specific gene signature composed of 11 potential prognostic biomarkers to be further investigated.

Artificial intelligence framework identifies candidate targets for drug repurposing in Alzheimer's disease

BACKGROUND

Genome-wide association studies (GWAS) have identified numerous susceptibility loci for Alzheimer's disease (AD). However, utilizing GWAS and multi-omics data to identify high-confidence AD risk genes (ARGs) and druggable targets that can guide development of new therapeutics for patients suffering from AD has heretofore not been successful.

METHODS

To address this critical problem in the field, we have developed a network-based artificial intelligence framework that is capable of integrating multi-omics data along with human protein-protein interactome networks to accurately infer accurate drug targets impacted by GWAS-identified variants to identify new therapeutics. When applied to AD, this approach integrates GWAS findings, multi-omics data from brain samples of AD patients and AD transgenic animal models, drug-target networks, and the human protein-protein interactome, along with large-scale patient database validation and in vitro mechanistic observations in human microglia cells.

RESULTS

Through this approach, we identified 103 ARGs validated by various levels of pathobiological evidence in AD. Via network-based prediction and population-based validation, we then showed that three drugs (pioglitazone, febuxostat, and atenolol) are significantly associated with decreased risk of AD compared with matched control populations. Pioglitazone usage is significantly associated with decreased risk of AD (hazard ratio (HR) = 0.916, 95% confidence interval [CI] 0.861-0.974, P = 0.005) in a retrospective case-control validation. Pioglitazone is a peroxisome proliferator-activated receptor (PPAR) agonist used to treat type 2 diabetes, and propensity score matching cohort studies confirmed its association with reduced risk of AD in comparison to glipizide (HR = 0.921, 95% CI 0.862-0.984, P = 0.0159), an insulin secretagogue that is also used to treat type 2 diabetes. In vitro experiments showed that pioglitazone downregulated glycogen synthase kinase 3 beta (GSK3β) and cyclin-dependent kinase (CDK5) in human microglia cells, supporting a possible mechanism-of-action for its beneficial effect in AD.

CONCLUSIONS

In summary, we present an integrated, network-based artificial intelligence methodology to rapidly translate GWAS findings and multi-omics data to genotype-informed therapeutic discovery in AD.

SWIMmeR: an R-based software to unveiling crucial nodes in complex biological networks

SUMMARY

We present SWIMmeR, an open-source version of its predecessor SWIM (SWitchMiner) that is a network-based tool for mining key (switch) genes that are associated with intriguing patterns of molecular co-abundance and may play a crucial role in phenotypic transitions in various biological settings. SWIM was originally written in MATLAB®, a proprietary programming language that requires the purchase of a license to install, manipulate, operate and run the software. Over the last years, SWIM has sparked a widespread interest within the scientific community thanks to the promising results obtained through its application in a broad range of phenotype-specific scenarios, spanning from complex diseases to grapevine berry maturation. This success has created the call for it to be distributed in a freely accessible, open-source, runtime environment, such as R, aimed at a general audience of non-expert users that cannot afford the leading proprietary solution. SWIMmeR is provided as a comprehensive collection of R functions and it also includes several additional features that make it less intensive in terms of computer time and more efficient in terms of usability and further implementation and extension.

AVAILABILITY AND IMPLEMENTATION

The SWIMmeR source code is freely available at https://github.com/sportingCode/SWIMmeR.git, along with a practical user guide, including a usage example of its application on breast cancer dataset.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Changes in gene-gene interactions associated with cancer onset and progression are largely independent of changes in gene expression

Recent findings indicate that changes underlying cancer onset and progression are not only attributable to changes in DNA structure and expression of individual genes but to changes in interactions among these genes as well. We examined co-expression changes in gene-network structure occurring during the onset and progression of nine different cancer types. Network complexity is generally reduced in the transition from normal precursor tissues to corresponding primary tumors. Cross-tissue cancer network similarity generally increases in early-stage cancers followed by a subsequent loss in cross-tissue cancer similarity as tumors reacquire cancer-specific network complexity. Gene-gene connections remaining stable through cancer development are enriched for “housekeeping” gene functions, whereas newly acquired interactions are associated with established cancer-promoting functions. Surprisingly, >90% of changes in gene-gene network interactions in cancers are not associated with changes in the expression of network genes relative to normal precursor tissues.

•

Gene-gene network complexity is reduced in the transition from normal to cancer

•

Network similarity across cancer types is higher in early-stage versus late-stage cancers

•

Network interactions among housekeeping genes are stable through cancer development

•

<10% of changes in network interactions in cancer involve changes in gene expression

Gene Network Analysis of Hepatocellular Carcinoma Identifies Modules Associated with Disease Progression, Survival, and Chemo Drug Resistance

Background

Hepatocellular carcinoma (HCC) is the second leading cause of cancer-related mortality worldwide. HCC transcriptome has been extensively studied; however, the progress in disease mechanisms, prognosis, and treatment is still slow.

Methods

A rank-based module-centric workflow was introduced to analyze important modules associated with HCC development, prognosis, and drug resistance. The currently largest HCC cell line RNA-Seq dataset from the LIMORE database was used to construct the reference modules by weighted gene co-expression network analysis.

Results

Thirteen reference modules were identified with validated reproducibility. These modules were all associated with specific biological functions. Differentially expressed module analysis revealed the crucial modules during HCC development. Modules and hub genes are indicative of patient survival. Modules can differentiate patients in different HCC stages. Furthermore, drug resistance was revealed by drug-module association analysis. Based on differentially expressed modules and hub genes, six candidate drugs were screened. The hub genes of those modules merit further investigation.

Conclusion

We proposed a reference module-based analysis of the HCC transcriptome. The identified modules are associated with HCC development, survival, and drug resistance. M3 and M6 may play important roles during HCV to HCC development. M1, M3, M5, and M7 are associated with HCC survival. High M4, high M9, low M1, and low M3 may be associated with dasatinib, doxorubicin, CD532, and simvastatin resistance. Our analysis provides useful information for HCC diagnosis and treatment.

Novel cancer subtyping method based on patient-specific gene regulatory network

The identification of cancer subtypes is important for the understanding of tumor heterogeneity. In recent years, numerous computational methods have been proposed for this problem based on the multi-omics data of patients. It is widely accepted that different cancer subtypes are induced by different molecular regulatory networks. However, only a few incorporate the differences between their molecular systems into the identification processes. In this study, we present a novel method to identify cancer subtypes based on patient-specific molecular systems. Our method realizes this by quantifying patient-specific gene networks, which are estimated from their transcriptome data, and by clustering their quantified networks. Comprehensive analyses of The Cancer Genome Atlas (TCGA) datasets applied to our method confirmed that they were able to identify more clinically meaningful cancer subtypes than the existing subtypes and found that the identified subtypes comprised different molecular features. Our findings also show that the proposed method can identify the novel cancer subtypes even with single omics data, which cannot otherwise be captured by existing methods using multi-omics data.

Modular network inference between miRNA-mRNA expression profiles using weighted co-expression network analysis

Connecting transcriptional and post-transcriptional regulatory networks solves an important puzzle in the elucidation of gene regulatory mechanisms. To decipher the complexity of these connections, we build co-expression network modules for mRNA as well as miRNA expression profiles of breast cancer data. We construct gene and miRNA co-expression modules using the weighted gene co-expression network analysis (WGCNA) method and establish the significance of these modules (Genes/miRNAs) for cancer phenotype. This work also infers an interaction network between the genes of the turquoise module from mRNA expression data and hubs of the turquoise module from miRNA expression data. A pathway enrichment analysis using a miRsystem web tool for miRNA hubs and some of their targets, reveal their enrichment in several important pathways associated with the progression of cancer.

Computational approaches leveraging integrated connections of multi-omic data toward clinical applications

In line with the advances in high-throughput technologies, multiple omic datasets have accumulated to study biological systems and diseases coherently. No single omics data type is capable of fully representing cellular activity. The complexity of the biological processes arises from the interactions between omic entities such as genes, proteins, and metabolites. Therefore, multi-omic data integration is crucial but challenging. The impact of the molecular alterations in multi-omic data is not local in the neighborhood of the altered gene or protein; rather, the impact diffuses in the network and changes the functionality of multiple signaling pathways and regulation of the gene expression. Additionally, multi-omic data is high-dimensional and has background noise. Several integrative approaches have been developed to accurately interpret the multi-omic datasets, including machine learning, network-based methods, and their combination. In this review, we overview the most recent integrative approaches and tools with a focus on network-based methods. We then discuss these approaches according to their specific applications, from disease-network and biomarker identification to patient stratification, drug discovery, and repurposing.

MOSES: A New Approach to Integrate Interactome Topology and Functional Features for Disease Gene Prediction

Disease gene prediction is to date one of the main computational challenges of precision medicine. It is still uncertain if disease genes have unique functional properties that distinguish them from other non-disease genes or, from a network perspective, if they are located randomly in the interactome or show specific patterns in the network topology. In this study, we propose a new method for disease gene prediction based on the use of biological knowledge-bases (gene-disease associations, genes functional annotations, etc.) and interactome network topology. The proposed algorithm called MOSES is based on the definition of two somewhat opposing sets of genes both disease-specific from different perspectives: warm seeds (i.e., disease genes obtained from databases) and cold seeds (genes far from the disease genes on the interactome and not involved in their biological functions). The application of MOSES to a set of 40 diseases showed that the suggested putative disease genes are significantly enriched in their reference disease. Reassuringly, known and predicted disease genes together, tend to form a connected network module on the human interactome, mitigating the scattered distribution of disease genes which is probably due to both the paucity of disease-gene associations and the incompleteness of the interactome.

Identification of fibronectin 1 (FN1) and complement component 3 (C3) as immune infiltration-related biomarkers for diabetic nephropathy using integrated bioinformatic analysis

Immune cell infiltration (ICI) plays a pivotal role in the development of diabetic nephropathy (DN). Evidence suggests that immune-related genes play an important role in the initiation of inflammation and the recruitment of immune cells. However, the underlying mechanisms and immune-related biomarkers in DN have not been elucidated. Therefore, this study aimed to explore immune-related biomarkers in DN and the underlying mechanisms using bioinformatic approaches. In this study, four DN glomerular datasets were downloaded, merged, and divided into training and test cohorts. First, we identified 55 differentially expressed immune-related genes; their biological functions were mainly enriched in leukocyte chemotaxis and neutrophil migration. The CIBERSORT algorithm was then used to evaluate the infiltrated immune cells; macrophages M1/M2, T cells CD8, and resting mast cells were strongly associated with DN. The ICI-related gene modules as well as 25 candidate hub genes were identified to construct a protein-protein interactive network and conduct molecular complex detection using the GOSemSim algorithm. Consequently, FN1, C3, and VEGFC were identified as immune-related biomarkers in DN, and a related transcription factor-miRNA-target network was constructed. Receiver operating characteristic curve analysis was estimated in the test cohort; FN1 and C3 had large area under the curve values (0.837 and 0.824, respectively). Clinical validation showed that FN1 and C3 were negatively related to the glomerular filtration rate in patients with DN. Six potential therapeutic small molecule compounds, such as calyculin, phenamil, and clofazimine, were discovered in the connectivity map. In conclusion, FN1 and C3 are immune-related biomarkers of DN.

In silico drug repurposing in COVID-19: A network-based analysis

The SARS-CoV-2 pandemic is a worldwide public health emergency. Despite the beginning of a vaccination campaign, the search for new drugs to appropriately treat COVID-19 patients remains a priority. Drug repurposing represents a faster and cheaper method than de novo drug discovery. In this study, we examined three different network-based approaches to identify potentially repurposable drugs to treat COVID-19. We analyzed transcriptomic data from whole blood cells of patients with COVID-19 and 21 other related conditions, as compared with those of healthy subjects. In addition to conventionally used drugs (e.g., anticoagulants, antihistaminics, anti-TNFα antibodies, corticosteroids), unconventional candidate compounds, such as SCN5A inhibitors and drugs active in the central nervous system, were identified. Clinical judgment and validation through clinical trials are always mandatory before use of the identified drugs in a clinical setting.

Expression levels of chemokine (C-X-C motif) ligands CXCL1 and CXCL3 as prognostic biomarkers in rectal adenocarcinoma: evidence from Gene Expression Omnibus (GEO) analyses

Rectal cancer is a life‑threatening disease worldwide. Chemotherapy resistance is common in rectal adenocarcinoma patients and has unfavorable survival outcomes; however, its related molecular mechanisms remain unknown. To identify genes related to the initiation and progression of rectal adenocarcinoma, three datasets were obtained from the Gene Expression Omnibus database. In total, differentially expressed genes were analyzed from 294 tumor and 277 para-carcinoma samples from patients with rectal cancer. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes functions were investigated. Cytoscape software and MicroRNA Enrichment Turned Network were applied to construct a protein-protein interaction network of the dependent hub genes and related microRNAs. The Oncomine database was used to identify hub genes. Additionally, Gene Expression Profiling Interactive Analysis was applied to determine the RNA expression level. Tumor immune infiltration was assessed using the Tumor Immune Estimation Resource database. The expression profiles of hub genes between stages, and their prognostic value, were also evaluated. During this study, data from The Cancer Genome Atlas were utilized. In rectal adenocarcinoma, four hub genes including CXCL1, CXCL2, CXCL3, and GNG4 were highly expressed at the gene and RNA levels. The expression of CXCL1, CXCL2, and CXCL3 was regulated by has-miR-1-3p and had a strong positive correlation with macrophage and neutrophil. CXCL2 and CXCL3 were differentially expressed at different tumor stages. High expression levels of CXCL1 and CXCL3 predicted poor survival. In conclusion, the CXCL1 and CXCL3 genes may have potential for prognosis and molecular targeted therapy of rectal adenocarcinoma.