signatureSearch: environment for gene expression signature searching and functional interpretation

Integrated edge information and pathway topology for drug-disease associations

Drug repurposing is a promising approach to find new therapeutic indications for approved drugs. Many computational approaches have been proposed to prioritize candidate anticancer drugs by gene or pathway level. However, these methods neglect the changes in gene interactions at the edge level. To address the limitation, we develop a computational drug repurposing method (iEdgePathDDA) based on edge information and pathway topology. First, we identify drug-induced and disease-related edges (the changes in gene interactions) within pathways by using the Pearson correlation coefficient. Next, we calculate the inhibition score between drug-induced edges and disease-related edges. Finally, we prioritize drug candidates according to the inhibition score on all disease-related edges. Case studies show that our approach successfully identifies new drug-disease pairs based on CTD database. Compared to the state-of-the-art approaches, the results demonstrate our method has the superior performance in terms of five metrics across colorectal, breast, and lung cancer datasets.

Signature reversion of three disease-associated gene signatures prioritizes cancer drug repurposing candidates

Drug repurposing is promising because approving a drug for a new indication requires fewer resources than approving a new drug. Signature reversion detects drug perturbations most inversely related to the disease-associated gene signature to identify drugs that may reverse that signature. We assessed the performance and biological relevance of three approaches for constructing disease-associated gene signatures (i.e., limma, DESeq2, and MultiPLIER) and prioritized the resulting drug repurposing candidates for four low-survival human cancers. Our results were enriched for candidates that had been used in clinical trials or performed well in the PRISM drug screen. Additionally, we found that pamidronate and nimodipine, drugs predicted to be efficacious against the brain tumor glioblastoma (GBM), inhibited the growth of a GBM cell line and cells isolated from a patient-derived xenograft (PDX). Our results demonstrate that by applying multiple disease-associated gene signature methods, we prioritized several drug repurposing candidates for low-survival cancers.

Unraveling athletic performance: Transcriptomics and external load monitoring in handball competition

OBJECTIVE

This study aims to comprehend the impact of handball practice on sub-elite athletes by investigating transcriptomic changes that occur during a match. The primary focus encompasses a dual objective: firstly, to identify and characterize these transcriptomic alterations, and secondly, to establish correlations between internal factors (gene expression), and external loads measured through Electronic Performance and Tracking Systems (EPTS variables). Ultimately, this comprehensive analysis seeks to evaluate both acute and chronic responses to exercise within the context of handball training.

METHODS

The study included sixteen elite male athletes from the FC Barcelona handball second team. Blood samples were extracted at three different time points: before the match at baseline levels (T1), immediately upon completion (T2), and 24 hours after completion (T3). Differential gene expression, Gene Ontology Term and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were conducted in two comparisons: Comparison 1 (T1 vs T2) and Comparison 2 (T1 vs T3). Further, the correlation between gene expression levels and training variables (external load) was conducted.

RESULTS

In T1 vs T2, 3717 of the 14632 genes detected were differentially expressed (adjusted p-value < 0.05), and enrichment of terms related to the immune system, mitochondria, and metabolic processes was found. Further, significant linear correlations were obtained between High-Speed running (HSR) and high-intensity variables such as acceleration ACC and deceleration DEC values with amino acids, and inflammatory and oxidative environment-related pathways, both in chronic and acute response.

CONCLUSIONS

This research highlights the effects of external workload on elite athletes during a handball match and throughout the season. The study identifies deregulation in the immune system, mitochondrial functions, and various metabolic pathways during the match. Additionally, it establishes correlations between the external load and pathways associated with amino acids, inflammation, oxidative environment, and regulation. These findings offer insights into the immediate and chronic responses of athletes to physical effort.

A network-based drug prioritization and combination analysis for the MEK5/ERK5 pathway in breast cancer

BACKGROUND

Prioritizing candidate drugs based on genome-wide expression data is an emerging approach in systems pharmacology due to its holistic perspective for preclinical drug evaluation. In the current study, a network-based approach was proposed and applied to prioritize plant polyphenols and identify potential drug combinations in breast cancer. We focused on MEK5/ERK5 signalling pathway genes, a recently identified potential drug target in cancer with roles spanning major carcinogenesis processes.

RESULTS

By constructing and identifying perturbed protein-protein interaction networks for luminal A breast cancer, plant polyphenols and drugs from transcriptome data, we first demonstrated their systemic effects on the MEK5/ERK5 signalling pathway. Subsequently, we applied a pathway-specific network pharmacology pipeline to prioritize plant polyphenols and potential drug combinations for use in breast cancer. Our analysis prioritized genistein among plant polyphenols. Drug combination simulations predicted several FDA-approved drugs in breast cancer with well-established pharmacology as candidates for target network synergistic combination with genistein. This study also highlights the concept of target network enhancer drugs, with drugs previously not well characterised in breast cancer being prioritized for use in the MEK5/ERK5 pathway in breast cancer.

CONCLUSION

This study proposes a computational framework for drug prioritization and combination with the MEK5/ERK5 signaling pathway in breast cancer. The method is flexible and provides the scientific community with a robust method that can be applied to other complex diseases.

Prioritized polycystic kidney disease drug targets and repurposing candidates from pre-cystic and cystic mouse Pkd2 model gene expression reversion

BACKGROUND

Autosomal dominant polycystic kidney disease (ADPKD) is one of the most prevalent monogenic human diseases. It is mostly caused by pathogenic variants in PKD1 or PKD2 genes that encode interacting transmembrane proteins polycystin-1 (PC1) and polycystin-2 (PC2). Among many pathogenic processes described in ADPKD, those associated with cAMP signaling, inflammation, and metabolic reprogramming appear to regulate the disease manifestations. Tolvaptan, a vasopressin receptor-2 antagonist that regulates cAMP pathway, is the only FDA-approved ADPKD therapeutic. Tolvaptan reduces renal cyst growth and kidney function loss, but it is not tolerated by many patients and is associated with idiosyncratic liver toxicity. Therefore, additional therapeutic options for ADPKD treatment are needed.

METHODS

As drug repurposing of FDA-approved drug candidates can significantly decrease the time and cost associated with traditional drug discovery, we used the computational approach signature reversion to detect inversely related drug response gene expression signatures from the Library of Integrated Network-Based Cellular Signatures (LINCS) database and identified compounds predicted to reverse disease-associated transcriptomic signatures in three publicly available Pkd2 kidney transcriptomic data sets of mouse ADPKD models. We focused on a pre-cystic model for signature reversion, as it was less impacted by confounding secondary disease mechanisms in ADPKD, and then compared the resulting candidates' target differential expression in the two cystic mouse models. We further prioritized these drug candidates based on their known mechanism of action, FDA status, targets, and by functional enrichment analysis.

RESULTS

With this in-silico approach, we prioritized 29 unique drug targets differentially expressed in Pkd2 ADPKD cystic models and 16 prioritized drug repurposing candidates that target them, including bromocriptine and mirtazapine, which can be further tested in-vitro and in-vivo.

CONCLUSION

Collectively, these results indicate drug targets and repurposing candidates that may effectively treat pre-cystic as well as cystic ADPKD.

Repurposing Approved Drugs for Sarcopenia Based on Transcriptomics Data in Humans

Sarcopenia, characterized by age-related loss of muscle mass, strength, and decreased physical performance, is a growing public health challenge amid the rapidly ageing population. As there are no approved drugs that target sarcopenia, it has become increasingly urgent to identify promising pharmacological interventions. In this study, we conducted an integrative drug repurposing analysis utilizing three distinct approaches. Firstly, we analyzed skeletal muscle transcriptomic sequencing data in humans and mice using gene differential expression analysis, weighted gene co-expression analysis, and gene set enrichment analysis. Subsequently, we employed gene expression profile similarity assessment, hub gene expression reversal, and disease-related pathway enrichment to identify and repurpose candidate drugs, followed by the integration of findings with rank aggregation algorithms. Vorinostat, the top-ranking drug, was also validated in an in vitro study, which demonstrated its efficacy in promoting muscle fiber formation. Although still requiring further validation in animal models and human clinical trials, these results suggest a promising drug repurposing prospect in the treatment and prevention of sarcopenia.

Optimizing Translational Research for Exceptional Health and Life Span: A Systematic Narrative of Studies to Identify Translatable Therapeutic Target(s) for Exceptional Health Span in Humans

BACKGROUND

Exceptional longevity as manifested by the lower incidence and delayed onset of age-related disabilities/diseases that include cardiovascular disease, Alzheimer's disease, cancer is believed to be influenced by inherent protective molecular factors in exceptionally long-lived individuals. Unraveling these protective factors could lead to the discovery of therapeutic target(s) and interventions to promote healthy aging.

METHODS

In this context, the National Institute on Aging has established a collection of translational longevity research projects (ie, the Long-Life Family Study, the Longevity Consortium, Longevity Genomics, and the Integrative Longevity Omics) which are generating large omics data sets spanning the human genome to phenome and have embarked on cross-species multiomic data analyses integrating human and nonhuman species that display wide variation in their life spans.

RESULTS

It is expected that these studies will discover key signaling pathways that influence exceptional health span and identify therapeutic targets for translation to enhance health and life span. Other efforts related to translational longevity research include the "Comprehensive Evaluation of Aging-Related Clinical Outcomes and Geroproteins study," which focuses on potential effects in humans of polypeptides/proteins whose circulating levels change with age, and for which experimental evidence indicates reversal or acceleration of aging changes. The "Predictive Human Mechanistic Markers Network" is devoted to the development of predictive markers of aging, for target engagement when testing novel interventions for healthy aging.

CONCLUSION

We describe here the significance, the unique study design, categories of data sets, analytical strategies, and a data portal to facilitate open science and sharing of resources from these longevity studies to identify and validate potential therapeutic targets for healthy aging.

Single-cell multiomics analysis reveals regulatory programs in clear cell renal cell carcinoma

The clear cell renal cell carcinoma (ccRCC) microenvironment consists of many different cell types and structural components that play critical roles in cancer progression and drug resistance, but the cellular architecture and underlying gene regulatory features of ccRCC have not been fully characterized. Here, we applied single-cell RNA sequencing (scRNA-seq) and single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq) to generate transcriptional and epigenomic landscapes of ccRCC. We identified tumor cell-specific regulatory programs mediated by four key transcription factors (TFs) (HOXC5, VENTX, ISL1, and OTP), and these TFs have prognostic significance in The Cancer Genome Atlas (TCGA) database. Targeting these TFs via short hairpin RNAs (shRNAs) or small molecule inhibitors decreased tumor cell proliferation. We next performed an integrative analysis of chromatin accessibility and gene expression for CD8⁺ T cells and macrophages to reveal the different regulatory elements in their subgroups. Furthermore, we delineated the intercellular communications mediated by ligand–receptor interactions within the tumor microenvironment. Taken together, our multiomics approach further clarifies the cellular heterogeneity of ccRCC and identifies potential therapeutic targets.

Considerations and challenges for sex-aware drug repurposing

Sex differences are essential factors in disease etiology and manifestation in many diseases such as cardiovascular disease, cancer, and neurodegeneration [33]. The biological influence of sex differences (including genomic, epigenetic, hormonal, immunological, and metabolic differences between males and females) and the lack of biomedical studies considering sex differences in their study design has led to several policies. For example, the National Institute of Health's (NIH) sex as a biological variable (SABV) and Sex and Gender Equity in Research (SAGER) policies to motivate researchers to consider sex differences [204]. However, drug repurposing, a promising alternative to traditional drug discovery by identifying novel uses for FDA-approved drugs, lacks sex-aware methods that can improve the identification of drugs that have sex-specific responses [7, 11, 14, 33]. Sex-aware drug repurposing methods either select drug candidates that are more efficacious in one sex or deprioritize drug candidates based on if they are predicted to cause a sex-bias adverse event (SBAE), unintended therapeutic effects that are more likely to occur in one sex. Computational drug repurposing methods are encouraging approaches to develop for sex-aware drug repurposing because they can prioritize sex-specific drug candidates or SBAEs at lower cost and time than traditional drug discovery. Sex-aware methods currently exist for clinical, genomic, and transcriptomic information [1, 7, 155]. They have not expanded to other data types, such as DNA variation, which has been beneficial in other drug repurposing methods that do not consider sex [114]. Additionally, some sex-aware methods suffer from poorer performance because a disproportionate number of male and female samples are available to train computational methods [7]. However, there is development potential for several different categories (i.e., data mining, ligand binding predictions, molecular associations, and networks). Low-dimensional representations of molecular association and network approaches are also especially promising candidates for future sex-aware drug repurposing methodologies because they reduce the multiple hypothesis testing burden and capture sex-specific variation better than the other methods [151, 159]. Here we review how sex influences drug response, the current state of drug repurposing including with respect to sex-bias drug response, and how model organism study design choices influence drug repurposing validation.

Metabolism-Associated Gene Signatures for FDG Avidity on PET/CT and Prognostic Validation in Hepatocellular Carcinoma

Introduction

The prognostic value of F-18 fluorodeoxyglucose positron emission tomography/computed tomography (FDG PET/CT) in hepatocellular carcinoma (HCC) was established in previous reports. However, there is no evidence suggesting the prognostic value of transcriptomes associated with tumor FDG uptake in HCC. It was aimed to elucidate metabolic genes and functions associated with FDG uptake, followed by assessment of those prognostic value.

Methods

Sixty HCC patients with Edmondson–Steiner grade II were included. FDG PET/CT scans were performed before any treatment. RNA sequencing data were obtained from tumor and normal liver tissue. Associations between each metabolism-associated gene and tumor FDG uptake were investigated by Pearson correlation analyses. A novel score between glucose and lipid metabolism-associated gene expression was calculated. In The Cancer Genome Atlas Liver Hepatocellular Carcinoma dataset, the prognostic power of selected metabolism-associated genes and a novel score was evaluated for external validation.

Results

Nine genes related to glycolysis and the HIF-1 signaling pathway showed positive correlations with tumor FDG uptake; 21 genes related to fatty acid metabolism and the PPAR signaling pathway demonstrated negative correlations. Seven potential biomarker genes, PFKFB4, ALDOA, EGLN3, EHHADH, GAPDH, HMGCS2, and ENO2 were identified. A metabolic gene expression balance score according to the dominance between glucose and lipid metabolism demonstrated good prognostic value in HCC.

Conclusions

The transcriptomic evidence of this study strongly supports the prognostic power of FDG PET/CT and indicates the potential usefulness of FDG PET/CT imaging biomarkers to select appropriate patients for metabolism-targeted therapy in HCC.

Host metabolic reprogramming in response to SARS-CoV-2 infection: A systems biology approach

Understanding the pathogenesis of SARS-CoV-2 is essential for developing effective treatment strategies. Viruses hijack the host metabolism to redirect the resources for their replication and survival. The influence of SARS-CoV-2 on host metabolism is yet to be fully understood. In this study, we analyzed the transcriptomic data obtained from different human respiratory cell lines and patient samples (nasopharyngeal swab, peripheral blood mononuclear cells, lung biopsy, bronchoalveolar lavage fluid) to understand metabolic alterations in response to SARS-CoV-2 infection. We explored the expression pattern of metabolic genes in the comprehensive genome-scale network model of human metabolism, Recon3D, to extract key metabolic genes, pathways, and reporter metabolites under each SARS-CoV-2-infected condition. A SARS-CoV-2 core metabolic interactome was constructed for network-based drug repurposing. Our analysis revealed the host-dependent dysregulation of glycolysis, mitochondrial metabolism, amino acid metabolism, nucleotide metabolism, glutathione metabolism, polyamine synthesis, and lipid metabolism. We observed different pro- and antiviral metabolic changes and generated hypotheses on how the host metabolism can be targeted for reducing viral titers and immunomodulation. These findings warrant further exploration with more samples and in vitro studies to test predictions.

Reconciling multiple connectivity scores for drug repurposing

The basis of several recent methods for drug repurposing is the key principle that an efficacious drug will reverse the disease molecular ‘signature’ with minimal side effects. This principle was defined and popularized by the influential ‘connectivity map’ study in 2006 regarding reversal relationships between disease- and drug-induced gene expression profiles, quantified by a disease-drug ‘connectivity score.’ Over the past 15 years, several studies have proposed variations in calculating connectivity scores toward improving accuracy and robustness in light of massive growth in reference drug profiles. However, these variations have been formulated inconsistently using various notations and terminologies even though they are based on a common set of conceptual and statistical ideas. Therefore, we present a systematic reconciliation of multiple disease-drug similarity metrics (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$ES$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$css$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$Sum$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$Cosine$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$XSum$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$XCor$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$XSpe$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$XCos$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$EWCos$\end{document}) and connectivity scores (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$CS$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$RGES$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$NCS$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$WCS$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$Tau$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$CSS$\end{document}, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$EMUDRA$\end{document}) by defining them using consistent notation and terminology. In addition to providing clarity and deeper insights, this coherent definition of connectivity scores and their relationships provides a unified scheme that newer methods can adopt, enabling the computational drug-development community to compare and investigate different approaches easily. To facilitate the continuous and transparent integration of newer methods, this article will be available as a live document (https://jravilab.github.io/connectivity_scores) coupled with a GitHub repository (https://github.com/jravilab/connectivity_scores) that any researcher can build on and push changes to.

Identification of 37 Heterogeneous Drug Candidates for Treatment of COVID-19 via a Rational Transcriptomics-Based Drug Repurposing Approach

A year after the initial outbreak, the COVID-19 pandemic caused by SARS-CoV-2 virus remains a serious threat to global health, while current treatment options are insufficient to bring major improvements. The aim of this study is to identify repurposable drug candidates with a potential to reverse transcriptomic alterations in the host cells infected by SARS-CoV-2. We have developed a rational computational pipeline to filter publicly available transcriptomic datasets of SARS-CoV-2-infected biosamples based on their responsiveness to the virus, to generate a list of relevant differentially expressed genes, and to identify drug candidates for repurposing using LINCS connectivity map. Pathway enrichment analysis was performed to place the results into biological context. We identified 37 structurally heterogeneous drug candidates and revealed several biological processes as druggable pathways. These pathways include metabolic and biosynthetic processes, cellular developmental processes, immune response and signaling pathways, with steroid metabolic process being targeted by half of the drug candidates. The pipeline developed in this study integrates biological knowledge with rational study design and can be adapted for future more comprehensive studies. Our findings support further investigations of some drugs currently in clinical trials, such as itraconazole and imatinib, and suggest 31 previously unexplored drugs as treatment options for COVID-19.