PathExt: a general framework for path-based mining of omics-integrated biological networks

Protocol for identifying key genes using network-based approach as an alternative to differential expression analysis

In a variety of biological contexts, characterizing genes associated with disease etiology and mediating global transcriptomic change is a key initial step. Here, we present a protocol to identify such key genes using our tool "PathExt," a tool that implements a network-based approach. We describe steps for installing libraries, preparing input data and detailed procedures for running PathExt, and characterizing differential pathways and key genes based on ripple centrality scores. For complete details on the use and execution of this protocol, please refer to Agrawal et al.1,2.

Network-based approach elucidates critical genes in BRCA subtypes and chemotherapy response in triple negative breast cancer

Breast cancers (BRCA) exhibit substantial transcriptional heterogeneity, posing a significant clinical challenge. The global transcriptional changes in a disease context, however, are likely mediated by few key genes which reflect disease etiology better than the differentially expressed genes (DEGs). We apply our network-based tool PathExt to 1,059 BRCA tumors across 4 subtypes to identify key mediator genes in each subtype. Compared to conventional differential expression analysis, PathExt-identified genes exhibit greater concordance across tumors, revealing shared and subtype-specific biological processes; better recapitulate BRCA-associated genes in multiple benchmarks, and are more essential in BRCA subtype-specific cell lines. Single-cell transcriptomic analysis reveals a subtype-specific distribution of PathExt-identified genes in multiple cell types from the tumor microenvironment. Application of PathExt to a TNBC chemotherapy response dataset identified subtype-specific key genes and biological processes associated with resistance. We described putative drugs that target key genes potentially mediating drug resistance.

PathTracer Comprehensively Identifies Hypoxia-Induced Dormancy Adaptations in Mycobacterium tuberculosis

Mining large-scale data to discover biologically relevant information remains a challenge despite the rapid development of bioinformatics tools. Here, we have developed a new tool, PathTracer, to identify biologically relevant information flows by mining genome-wide protein-protein interaction networks following integration of gene expression data. PathTracer successfully mines interactions between genes and traces the most perturbed paths of perceived activities under the conditions of the study. We further demonstrated the utility of this tool by identifying adaptation mechanisms of hypoxia-induced dormancy in Mycobacterium tuberculosis (Mtb).

Systems-level profiling of early peripheral host-response landscape variations across COVID-19 severity states in an Indian cohort

Host immune response to COVID-19 plays a significant role in regulating disease severity. Although big data analysis has provided significant insights into the host biology of COVID-19 across the world, very few such studies have been performed in the Indian population. This study utilizes a transcriptome-integrated network analysis approach to compare the immune responses between asymptomatic or mild and moderate-severe COVID-19 patients in an Indian cohort. An immune suppression phenotype is observed in the early stages of moderate-severe COVID-19 manifestation. A number of pathways are identified that play crucial roles in the host control of the disease such as the type I interferon response and classical complement pathway which show different activity levels across the severity spectrum. This study also identifies two transcription factors, IRF7 and ESR1, to be important in regulating the severity of COVID-19. Overall this study provides a deep understanding of the peripheral immune landscape in the COVID-19 severity spectrum in the Indian genetic background and opens up future research avenues to compare immune responses across global populations.

Network-based approach elucidates critical genes in BRCA subtypes and chemotherapy response in Triple Negative Breast Cancer

Breast cancers exhibit substantial transcriptional heterogeneity, posing a significant challenge to the prediction of treatment response and prognostication of outcomes. Especially, translation of TNBC subtypes to the clinic remains a work in progress, in part because of a lack of clear transcriptional signatures distinguishing the subtypes. Our recent network-based approach, PathExt, demonstrates that global transcriptional changes in a disease context are likely mediated by a small number of key genes, and these mediators may better reflect functional or translationally relevant heterogeneity. We apply PathExt to 1059 BRCA tumors and 112 healthy control samples across 4 subtypes to identify frequent, key-mediator genes in each BRCA subtype. Compared to conventional differential expression analysis, PathExt-identified genes (1) exhibit greater concordance across tumors, revealing shared as well as BRCA subtype-specific biological processes, (2) better recapitulate BRCA-associated genes in multiple benchmarks, and (3) exhibit greater dependency scores in BRCA subtype-specific cancer cell lines. Single cell transcriptomes of BRCA subtype tumors reveal a subtype-specific distribution of PathExt-identified genes in multiple cell types from the tumor microenvironment. Application of PathExt to a TNBC chemotherapy response dataset identified TNBC subtype-specific key genes and biological processes associated with resistance. We described putative drugs that target top novel genes potentially mediating drug resistance. Overall, PathExt applied to breast cancer refines previous views of gene expression heterogeneity and identifies potential mediators of TNBC subtypes, including potential therapeutic targets.

Identification of Probucol as a candidate for combination therapy with Metformin for Type 2 diabetes

Type 2 Diabetes (T2D) is often managed with metformin as the drug of choice. While it is effective overall, many patients progress to exhibit complications. Strategic drug combinations to tackle this problem would be useful. We constructed a genome-wide protein-protein interaction network capturing a global perspective of perturbations in diabetes by integrating T2D subjects' transcriptomic data. We computed a 'frequently perturbed subnetwork' in T2D that captures common perturbations across tissue types and mapped the possible effects of Metformin onto it. We then identified a set of remaining T2D perturbations and potential drug targets among them, related to oxidative stress and hypercholesterolemia. We then identified Probucol as the potential co-drug for adjunct therapy with Metformin and evaluated the efficacy of the combination in a rat model of diabetes. We find Metformin-Probucol at 5:0.5 mg/kg effective in restoring near-normal serum glucose, lipid, and cholesterol levels.

Profiling system-wide variations and similarities between Rheumatic Heart Disease and Acute Rheumatic Fever-A pilot analysis

Rheumatic heart disease (RHD) continues to affect developing countries with low income due to the lack of resources and effective diagnostic techniques. Understanding the genetic basis common to both the diseases and that of progression from its prequel disease state, Acute Rheumatic Fever (ARF), would aid in developing predictive biomarkers and improving patient care. To gain system-wide molecular insights into possible causes for progression, in this pilot study, we collected blood transcriptomes from ARF (5) and RHD (5) patients. Using an integrated transcriptome and network analysis approach, we identified a subnetwork comprising the most significantly differentially expressed genes and most perturbed pathways in RHD compared to ARF. For example, the chemokine signaling pathway was seen to be upregulated, while tryptophan metabolism was found to be downregulated in RHD. The subnetworks of variation between the two conditions provide unbiased molecular-level insights into the host processes that may be linked with the progression of ARF to RHD, which has the potential to inform future diagnostics and therapeutic strategies. We also found a significantly raised neutrophil/lymphocyte ratio in both ARF and RHD cohorts. Activated neutrophils and inhibited Natural Killer cell gene signatures reflected the drivers of the inflammatory process typical to both disease conditions.

DrDimont: explainable drug response prediction from differential analysis of multi-omics networks

MOTIVATION

While it has been well established that drugs affect and help patients differently, personalized drug response predictions remain challenging. Solutions based on single omics measurements have been proposed, and networks provide means to incorporate molecular interactions into reasoning. However, how to integrate the wealth of information contained in multiple omics layers still poses a complex problem.

RESULTS

We present DrDimont, Drug response prediction from Differential analysis of multi-omics networks. It allows for comparative conclusions between two conditions and translates them into differential drug response predictions. DrDimont focuses on molecular interactions. It establishes condition-specific networks from correlation within an omics layer that are then reduced and combined into heterogeneous, multi-omics molecular networks. A novel semi-local, path-based integration step ensures integrative conclusions. Differential predictions are derived from comparing the condition-specific integrated networks. DrDimont's predictions are explainable, i.e. molecular differences that are the source of high differential drug scores can be retrieved. We predict differential drug response in breast cancer using transcriptomics, proteomics, phosphosite and metabolomics measurements and contrast estrogen receptor positive and receptor negative patients. DrDimont performs better than drug prediction based on differential protein expression or PageRank when evaluating it on ground truth data from cancer cell lines. We find proteomic and phosphosite layers to carry most information for distinguishing drug response.

AVAILABILITY AND IMPLEMENTATION

DrDimont is available on CRAN: https://cran.r-project.org/package=DrDimont.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A Path-Based Analysis of Infected Cell Line and COVID-19 Patient Transcriptome Reveals Novel Potential Targets and Drugs Against SARS-CoV-2

Most transcriptomic studies of SARS-CoV-2 infection have focused on differentially expressed genes, which do not necessarily reveal the genes mediating the transcriptomic changes. In contrast, exploiting curated biological network, our PathExt tool identifies central genes from the differentially active paths mediating global transcriptomic response. Here we apply PathExt to multiple cell line infection models of SARS-CoV-2 and other viruses, as well as to COVID-19 patient-derived PBMCs. The central genes mediating SARS-CoV-2 response in cell lines were uniquely enriched for ATP metabolic process, G1/S transition, leukocyte activation and migration. In contrast, PBMC response reveals dysregulated cell-cycle processes. In PBMC, the most frequently central genes are associated with COVID-19 severity. Importantly, relative to differential genes, PathExt-identified genes show greater concordance with several benchmark anti-COVID-19 target gene sets. We propose six novel anti-SARS-CoV-2 targets ADCY2, ADSL, OCRL, TIAM1, PBK, and BUB1, and potential drugs targeting these genes, such as Bemcentinib, Phthalocyanine, and Conivaptan.

A new blood-based RNA signature (R9), for monitoring effectiveness of tuberculosis treatment in a South Indian longitudinal cohort

Tuberculosis (TB) treatment involves a multidrug regimen for six months, and until two months, it is unclear if treatment is effective. This delay can lead to the evolution of drug resistance, lung damage, disease spread, and transmission. We identify a blood-based 9-gene signature using a computational pipeline that constructs and interrogates a genome-wide transcriptome-integrated protein-interaction network. The identified signature is able to determine treatment response at week 1-2 in three independent public datasets. Signature-based R9-score correctly detected treatment response at individual timepoints (204 samples) from a newly developed South Indian longitudinal cohort involving 32 patients with pulmonary TB. These results are consistent with conventional clinical metrics and can discriminate good from poor treatment responders at week 2 (AUC 0.93(0.81-1.00)). In this work, we provide proof of concept that the R9-score can determine treatment effectiveness, making a case for designing a larger clinical study.