Robust de novo pathway enrichment with KeyPathwayMiner 5

Personalized Protein-Protein Interaction Networks Towards Unraveling the Molecular Mechanisms of Alzheimer's Disease

Alzheimer's disease (AD) is a highly heterogenous neurodegenerative disease, and several omic-based datasets were generated in the last decade from the patients with the disease. However, the vast majority of studies evaluate these datasets in bulk by considering all the patients as a single group, which obscures the molecular differences resulting from the heterogeneous nature of the disease. In this study, we adopted a personalized approach and analyzed the transcriptome data from 403 patients individually by mapping the data on a human protein-protein interaction network. Patient-specific subnetworks were discovered and analyzed in terms of the genes in the subnetworks, enriched functional terms, and known AD genes. We identified several affected pathways that could not be captured by the bulk comparison. We also showed that our personalized findings point to patterns of alterations consistent with the recently suggested AD subtypes.

Network-based approaches for modeling disease regulation and progression

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

Cancer driver drug interaction explorer

Cancer is a heterogeneous disease characterized by unregulated cell growth and promoted by mutations in cancer driver genes some of which encode suitable drug targets. Since the distinct set of cancer driver genes can vary between and within cancer types, evidence-based selection of drugs is crucial for targeted therapy following the precision medicine paradigm. However, many putative cancer driver genes can not be targeted directly, suggesting an indirect approach that considers alternative functionally related targets in the gene interaction network. Once potential drug targets have been identified, it is essential to consider all available drugs. Since tools that offer support for systematic discovery of drug repurposing candidates in oncology are lacking, we developed CADDIE, a web application integrating six human gene-gene and four drug-gene interaction databases, information regarding cancer driver genes, cancer-type specific mutation frequencies, gene expression information, genetically related diseases, and anticancer drugs. CADDIE offers access to various network algorithms for identifying drug targets and drug repurposing candidates. It guides users from the selection of seed genes to the identification of therapeutic targets or drug candidates, making network medicine algorithms accessible for clinical research. CADDIE is available at https://exbio.wzw.tum.de/caddie/ and programmatically via a python package at https://pypi.org/project/caddiepy/.

De novo identification of maximally deregulated subnetworks based on multi-omics data with DeRegNet

Background

With a growing amount of (multi-)omics data being available, the extraction of knowledge from these datasets is still a difficult problem. Classical enrichment-style analyses require predefined pathways or gene sets that are tested for significant deregulation to assess whether the pathway is functionally involved in the biological process under study. De novo identification of these pathways can reduce the bias inherent in predefined pathways or gene sets. At the same time, the definition and efficient identification of these pathways de novo from large biological networks is a challenging problem.

Results

We present a novel algorithm, DeRegNet, for the identification of maximally deregulated subnetworks on directed graphs based on deregulation scores derived from (multi-)omics data. DeRegNet can be interpreted as maximum likelihood estimation given a certain probabilistic model for de-novo subgraph identification. We use fractional integer programming to solve the resulting combinatorial optimization problem. We can show that the approach outperforms related algorithms on simulated data with known ground truths. On a publicly available liver cancer dataset we can show that DeRegNet can identify biologically meaningful subgraphs suitable for patient stratification. DeRegNet can also be used to find explicitly multi-omics subgraphs which we demonstrate by presenting subgraphs with consistent methylation-transcription patterns. DeRegNet is freely available as open-source software.

Conclusion

The proposed algorithmic framework and its available implementation can serve as a valuable heuristic hypothesis generation tool contextualizing omics data within biomolecular networks.

KeyPathwayMineR: De Novo Pathway Enrichment in the R Ecosystem

De novo pathway enrichment is a systems biology approach in which OMICS data are projected onto a molecular interaction network to identify subnetworks representing condition-specific functional modules and molecular pathways. Compared to classical pathway enrichment analysis methods, de novo pathway enrichment is not limited to predefined lists of pathways from (curated) databases and thus particularly suited for discovering novel disease mechanisms. While several tools have been proposed for pathway enrichment, the integration of de novo pathway enrichment in end-to-end OMICS analysis workflows in the R programming language is currently limited to a single tool. To close this gap, we have implemented an R package KeyPathwayMineR (KPM-R). The package extends the features and usability of existing versions of KeyPathwayMiner by leveraging the power, flexibility and versatility of R and by providing various novel functionalities for performing data preparation, visualization, and comparison. In addition, thanks to its interoperability with a plethora of existing R packages in e.g., Bioconductor, CRAN, and GitHub, KPM-R allows carrying out the initial preparation of the datasets and to meaningfully interpret the extracted subnetworks. To demonstrate the package's potential, KPM-R was applied to bulk RNA-Seq data of nasopharyngeal swabs from SARS-CoV-2 infected individuals, and on single cell RNA-Seq data of aging mice tissue from the Tabula Muris Senis atlas.

CTCFL regulates the PI3K-Akt pathway and it is a target for personalized ovarian cancer therapy

High-grade serous ovarian carcinoma (HGSC) is the most lethal gynecologic malignancy due to the lack of reliable biomarkers, effective treatment, and chemoresistance. Improving the diagnosis and the development of targeted therapies is still needed. The molecular pathomechanisms driving HGSC progression are not fully understood though crucial for effective diagnosis and identification of novel targeted therapy options. The oncogene CTCFL (BORIS), the paralog of CTCF, is a transcriptional factor highly expressed in ovarian cancer (but in rarely any other tissue in females) with cancer-specific characteristics and therapeutic potential. In this work, we seek to understand the regulatory functions of CTCFL to unravel new target genes with clinical relevance. We used in vitro models to evaluate the transcriptional changes due to the presence of CTCFL, followed by a selection of gene candidates using de novo network enrichment analysis. The resulting mechanistic candidates were further assessed regarding their prognostic potential and druggability. We show that CTCFL-driven genes are involved in cytoplasmic membrane functions; in particular, the PI3K-Akt initiators EGFR1 and VEGFA, as well as ITGB3 and ITGB6 are potential drug targets. Finally, we identified the CTCFL targets ACTBL2, MALT1 and PCDH7 as mechanistic biomarkers to predict survival in HGSC. Finally, we elucidated the value of CTCFL in combination with its targets as a prognostic marker profile for HGSC progression and as putative drug targets.

Alveolar Regeneration in COVID-19 Patients: A Network Perspective

A viral infection involves entry and replication of viral nucleic acid in a host organism, subsequently leading to biochemical and structural alterations in the host cell. In the case of SARS-CoV-2 viral infection, over-activation of the host immune system may lead to lung damage. Albeit the regeneration and fibrotic repair processes being the two protective host responses, prolonged injury may lead to excessive fibrosis, a pathological state that can result in lung collapse. In this review, we discuss regeneration and fibrosis processes in response to SARS-CoV-2 and provide our viewpoint on the triggering of alveolar regeneration in coronavirus disease 2019 (COVID-19) patients.

Mapping Transcriptome Data to Protein-Protein Interaction Networks of Inflammatory Bowel Diseases Reveals Disease-Specific Subnetworks

Inflammatory bowel disease (IBD) is the common name for chronic disorders associated with the inflammation of the gastrointestinal tract. IBD is triggered by environmental factors in genetically susceptible individuals and has a significant number of incidences worldwide. Crohn’s disease (CD) and ulcerative colitis (UC) are the two distinct types of IBD. While involvement in ulcerative colitis is limited to the colon, Crohn’s disease may involve the whole gastrointestinal tract. Although these two disorders differ in macroscopic inflammation patterns, they share various molecular pathogenesis, yet the diagnosis can remain unclear, and it is important to reveal their molecular signatures in the network level. Improved molecular understanding may reveal disease type-specific and even individual-specific targets. To this aim, we determine the subnetworks specific to UC and CD by mapping transcriptome data to protein–protein interaction (PPI) networks using two different approaches [KeyPathwayMiner (KPM) and stringApp] and perform the functional enrichment analysis of the resulting disease type-specific subnetworks. TP63 was identified as the hub gene in the UC-specific subnet and p63 tumor protein, being in the same family as p53 and p73, has been studied in literature for the risk associated with colorectal cancer and IBD. APP was identified as the hub gene in the CD-specific subnet, and it has an important role in the pathogenesis of Alzheimer’s disease (AD). This relation suggests that some similar genetic factors may be effective in both AD and CD. Last, in order to understand the biological meaning of these disease-specific subnets, they were functionally enriched. It is important to note that chemokines—special types of cytokines—and antibacterial response are important in UC-specific subnets, whereas cytokines and antimicrobial responses as well as cancer-related pathways are important in CD-specific subnets. Overall, these findings reveal the differences between IBD subtypes at the molecular level and can facilitate diagnosis for UC and CD as well as provide potential molecular targets that are specific to disease subtypes.

Comparative transcriptome analysis reveals key epigenetic targets in SARS-CoV-2 infection

COVID-19 is an infection caused by SARS-CoV-2 (Severe Acute Respiratory Syndrome coronavirus 2), which has caused a global outbreak. Current research efforts are focused on the understanding of the molecular mechanisms involved in SARS-CoV-2 infection in order to propose drug-based therapeutic options. Transcriptional changes due to epigenetic regulation are key host cell responses to viral infection and have been studied in SARS-CoV and MERS-CoV; however, such changes are not fully described for SARS-CoV-2. In this study, we analyzed multiple transcriptomes obtained from cell lines infected with MERS-CoV, SARS-CoV, and SARS-CoV-2, and from COVID-19 patient-derived samples. Using integrative analyses of gene co-expression networks and de-novo pathway enrichment, we characterize different gene modules and protein pathways enriched with Transcription Factors or Epifactors relevant for SARS-CoV-2 infection. We identified EP300, MOV10, RELA, and TRIM25 as top candidates, and more than 60 additional proteins involved in the epigenetic response during viral infection that has therapeutic potential. Our results show that targeting the epigenetic machinery could be a feasible alternative to treat COVID-19.

DIGGER: exploring the functional role of alternative splicing in protein interactions

Alternative splicing plays a major role in regulating the functional repertoire of the proteome. However, isoform-specific effects to protein-protein interactions (PPIs) are usually overlooked, making it impossible to judge the functional role of individual exons on a systems biology level. We overcome this barrier by integrating protein-protein interactions, domain-domain interactions and residue-level interactions information to lift exon expression analysis to a network level. Our user-friendly database DIGGER is available at https://exbio.wzw.tum.de/digger and allows users to seamlessly switch between isoform and exon-centric views of the interactome and to extract sub-networks of relevant isoforms, making it an essential resource for studying mechanistic consequences of alternative splicing.

Individuating Possibly Repurposable Drugs and Drug Targets for COVID-19 Treatment Through Hypothesis-Driven Systems Medicine Using CoVex

Coronavirus disease 2019 (COVID-19), caused by the SARS-CoV-2 virus, has developed into a pandemic causing major disruptions and hundreds of thousands of deaths in wide parts of the world. As of July 3, 2020, neither vaccines nor approved drugs for effective treatment are available. In this article, we showcase how to individuate drug targets and potentially repurposable drugs in silico using CoVex a recently presented systems medicine platform for COVID-19 drug repurposing. Starting from initial hypotheses, CoVex leverages network algorithms to individuate host proteins involved in COVID-19 disease mechanisms, as well as existing drugs targeting these potential drug targets. Our analysis reveals GLA, PLAT, and GGCX as potential drug targets, and urokinase, argatroban, dabigatran etexilate, betrixaban, ximelagatran and anisindione as potentially repurposable drugs.

Integrative biomarker detection on high-dimensional gene expression data sets: a survey on prior knowledge approaches

Gene expression data provide the expression levels of tens of thousands of genes from several hundred samples. These data are analyzed to detect biomarkers that can be of prognostic or diagnostic use. Traditionally, biomarker detection for gene expression data is the task of gene selection. The vast number of genes is reduced to a few relevant ones that achieve the best performance for the respective use case. Traditional approaches select genes based on their statistical significance in the data set. This results in issues of robustness, redundancy and true biological relevance of the selected genes. Integrative analyses typically address these shortcomings by integrating multiple data artifacts from the same objects, e.g. gene expression and methylation data. When only gene expression data are available, integrative analyses instead use curated information on biological processes from public knowledge bases. With knowledge bases providing an ever-increasing amount of curated biological knowledge, such prior knowledge approaches become more powerful. This paper provides a thorough overview on the status quo of biomarker detection on gene expression data with prior biological knowledge. We discuss current shortcomings of traditional approaches, review recent external knowledge bases, provide a classification and qualitative comparison of existing prior knowledge approaches and discuss open challenges for this kind of gene selection.

Exploring the SARS-CoV-2 virus-host-drug interactome for drug repurposing

Coronavirus Disease-2019 (COVID-19) is an infectious disease caused by the SARS-CoV-2 virus. Various studies exist about the molecular mechanisms of viral infection. However, such information is spread across many publications and it is very time-consuming to integrate, and exploit. We develop CoVex, an interactive online platform for SARS-CoV-2 host interactome exploration and drug (target) identification. CoVex integrates virus-human protein interactions, human protein-protein interactions, and drug-target interactions. It allows visual exploration of the virus-host interactome and implements systems medicine algorithms for network-based prediction of drug candidates. Thus, CoVex is a resource to understand molecular mechanisms of pathogenicity and to prioritize candidate therapeutics. We investigate recent hypotheses on a systems biology level to explore mechanistic virus life cycle drivers, and to extract drug repurposing candidates. CoVex renders COVID-19 drug research systems-medicine-ready by giving the scientific community direct access to network medicine algorithms. It is available at https://exbio.wzw.tum.de/covex/.

A Genome-Wide Integrative Association Study of DNA Methylation and Gene Expression Data and Later Life Cognitive Functioning in Monozygotic Twins

Monozygotic twins are genetically identical but rarely phenotypically identical. Epigenetic and transcriptional variation could influence this phenotypic discordance. Investigation of intra-pair differences in molecular markers and a given phenotype in monozygotic twins controls most of the genetic contribution, enabling studies of the molecular features of the phenotype. This study aimed to identify genes associated with cognition in later life using integrated enrichment analyses of the results of blood-derived intra-pair epigenome-wide and transcriptome-wide association analyses of cognition in 452 middle-aged and old-aged monozygotic twins (56–80 years). Integrated analyses were performed with an unsupervised approach using KeyPathwayMiner, and a supervised approach using the KEGG and Reactome databases. The supervised approach identified several enriched gene sets, including “neuroactive ligand receptor interaction” (p-value = 1.62^∗10-2), “Neurotrophin signaling” (p-value = 2.52^∗10-3), “Alzheimer’s disease” (p-value = 1.20^∗10-2), and “long-term depression” (p-value = 1.62^∗10-2). The unsupervised approach resulted in a 238 gene network, including the Alzheimer’s disease gene APP (Amyloid Beta Precursor Protein) as an exception node, and several novel candidate genes. The strength of the unsupervised method is that it can reveal previously uncharacterized sub-pathways and detect interplay between biological processes, which remain undetected by the current supervised methods. In conclusion, this study identified several previously reported cognition genes and pathways and, additionally, puts forward novel candidates for further verification and validation.

De Novo Pathway Enrichment with KeyPathwayMiner

Biomolecular networks such as protein-protein interaction networks provide a static picture of the interplay of genes and their products, and, consequently, they fail to capture dynamic changes taking place during the development of complex diseases. KeyPathwayMiner is a software platform designed to fill this gap by integrating previous knowledge captured in molecular interaction networks with OMICS datasets (DNA microarrays, RNA sequencing, genome-wide methylation studies, etc.) to extract connected subnetworks with a high number of deregulated genes. This protocol describes how to use KeyPathwayMiner for integrated analysis of multi-omics datasets in the network analysis tool Cytoscape and in a stand-alone web application available at https://keypathwayminer.compbio.sdu.dk .

Network-Based Approaches for Multi-omics Integration

Network-based approach is rapidly emerging as a promising strategy to integrate and interpret different -omics datasets, including metabolomics. The first section of this chapter introduces the current progresses and main concepts in multi-omics integration. The second section provides an overview of the public resources available for creation of biological networks. The third section describes three common application scenarios including subnetwork identification, network-based enrichment analysis, and systems metabolomics. The section four introduces the concept of hierarchical community network analysis. The section five discusses different tools for network visualization. The chapter ends with a future perspective on multi-omics integration.

The epigenetic factor BORIS (CTCFL) controls the androgen receptor regulatory network in ovarian cancer

The identification of prognostic biomarkers is a priority for patients suffering from high-grade serous ovarian cancer (SOC), which accounts for >70% of ovarian cancer (OC) deaths. Meanwhile, borderline ovarian cancer (BOC) is a low malignancy tumor and usually patients undergo surgery with low probabilities of recurrence. However, SOC remains the most lethal neoplasm due to the lack of biomarkers for early diagnosis and prognosis. In this regard, BORIS (CTCFL), a CTCF paralog, is a promising cancer biomarker that is overexpressed and controls transcription in several cancer types, mainly in OC. Studies suggest that BORIS has an important function in OC by altering gene expression, but the effect and extent to which BORIS influences transcription in OC from a genome-wide perspective is unclear. Here, we sought to identify BORIS target genes in an OC cell line (OVCAR3) with potential biomarker use in OC tumor samples. To achieve this, we performed in vitro knockout and knockdown experiments of BORIS in OVCAR3 cell line followed by expression microarrays and bioinformatics network enrichment analysis to identify relevant BORIS target genes. In addition, ex vivo expression data analysis of 373 ovarian cancer patients were evaluated to identify the expression patterns of BORIS target genes. In vitro, we uncovered 130 differentially expressed genes and obtained the BORIS-associated regulatory network, in which the androgen receptor (AR) acts as a major transcription factor. Also, FN1, FAM129A, and CD97 genes, which are related to chemoresistance and metastases in OC, were identified. In SOC patients, we observed that malignancy is associated with high levels of BORIS expression while BOC patients show lower levels. Our study suggests that BORIS acts as a main regulator, and has the potential to be used as a prognostic biomarker and to yield novel drug targets among the genes BORIS controls in SOC patients.

De novo pathway-based biomarker identification

Gene expression profiles have been extensively discussed as an aid to guide the therapy by predicting disease outcome for the patients suffering from complex diseases, such as cancer. However, prediction models built upon single-gene (SG) features show poor stability and performance on independent datasets. Attempts to mitigate these drawbacks have led to the development of network-based approaches that integrate pathway information to produce meta-gene (MG) features. Also, MG approaches have only dealt with the two-class problem of good versus poor outcome prediction. Stratifying patients based on their molecular subtypes can provide a detailed view of the disease and lead to more personalized therapies. We propose and discuss a novel MG approach based on de novo pathways, which for the first time have been used as features in a multi-class setting to predict cancer subtypes. Comprehensive evaluation in a large cohort of breast cancer samples from The Cancer Genome Atlas (TCGA) revealed that MGs are considerably more stable than SG models, while also providing valuable insight into the cancer hallmarks that drive them. In addition, when tested on an independent benchmark non-TCGA dataset, MG features consistently outperformed SG models. We provide an easy-to-use web service at http://pathclass.compbio.sdu.dk where users can upload their own gene expression datasets from breast cancer studies and obtain the subtype predictions from all the classifiers.

Gaining insights into cancer biology through exploration of the cancer secretome using proteomic and bioinformatic tools

INTRODUCTION

Tumor-associated proteins released by cancer cells and by tumor stroma cells, referred as 'cancer secretome', represent a valuable resource for discovery of potential cancer biomarkers. The last decade was marked by a great increase in number of studies focused on various aspects of cancer secretome including, composition and identification of components externalized by malignant cells and by the components of tumor microenvironment. Areas covered: Here, we provide an overview of achievements in the proteomic analysis of the cancer secretome, elicited through the tumor-associated interstitial fluid recovered from malignant tissues ex vivo or the protein component of conditioned media obtained from cultured cancer cells in vitro. We summarize various bioinformatic tools and approaches and critically appraise their outcomes, focusing on problems and challenges that arise when applied for the analysis of cancer secretomic databases. Expert commentary: Recent achievements in the omics- analysis of structural and metabolic aspects of altered cancer secretome contribute greatly to the various hallmarks of cancer including the identification of clinically significant biomarkers and potential targets for therapeutic intervention.

Differential Effects of Vitamins A and D on the Transcriptional Landscape of Human Monocytes during Infection

Vitamin A and vitamin D are essential nutrients with a wide range of pleiotropic effects in humans. Beyond their well-documented roles in cellular differentiation, embryogenesis, tissue maintenance and bone/calcium homeostasis, both vitamins have attracted considerable attention due to their association with-immunological traits. Nevertheless, our knowledge of their immunomodulatory potential during infection is restricted to single gene-centric studies, which do not reflect the complexity of immune processes. In the present study, we performed a comprehensive RNA-seq-based approach to define the whole immunomodulatory role of vitamins A and D during infection. Using human monocytes as host cells, we characterized the differential role of both vitamins upon infection with three different pathogens: Aspergillus fumigatus, Candida albicans and Escherichia coli. Both vitamins showed an unexpected ability to counteract the pathogen-induced transcriptional responses. Upon infection, we identified 346 and 176 immune-relevant genes that were regulated by atRA and vitD, respectively. This immunomodulatory activity was dependent on the inflammatory stimulus, allowing us to distinguish regulatory patterns which were specific for each stimulatory setting. Moreover, we explored possible direct and indirect mechanisms of vitamin-mediated regulation of the immune response. Our findings highlight the importance of vitamin-monitoring in critically ill patients. Moreover, our results underpin the potential of atRA and vitD as therapeutic options for anti-inflammatory treatment.

Robust de novo pathway enrichment with KeyPathwayMiner 5

Identifying functional modules or novel active pathways, recently termed de novo pathway enrichment, is a computational systems biology challenge that has gained much attention during the last decade. Given a large biological interaction network, KeyPathwayMiner extracts connected subnetworks that are enriched for differentially active entities from a series of molecular profiles encoded as binary indicator matrices. Since interaction networks constantly evolve, an important question is how robust the extracted results are when the network is modified. We enable users to study this effect through several network perturbation techniques and over a range of perturbation degrees. In addition, users may now provide a gold-standard set to determine how enriched extracted pathways are with relevant genes compared to randomized versions of the original network.