metaGraphite-a new layer of pathway annotation to get metabolite networks

DNA methylation changes in association with trauma-focused psychotherapy efficacy in treatment-resistant depression patients: a prospective longitudinal study

Background: Stressful events increase the risk for treatment-resistant depression (TRD), and trauma-focused psychotherapy can be useful for TRD patients exposed to early life stress (ELS). Epigenetic processes are known to be related to depression and ELS, but there is no evidence of the effects of trauma-focused psychotherapy on methylation alterations.Objective: We performed the first epigenome-wide association study to investigate methylation changes related to trauma-focused psychotherapies effects in TRD patients.Method: Thirty TRD patients assessed for ELS underwent trauma-focused psychotherapy, of those, 12 received trauma-focused cognitive behavioural therapy, and 18 Eye Movement Desensitization and Reprocessing (EMDR). DNA methylation was profiled with Illumina Infinium EPIC array at T0 (baseline), after 8 weeks (T8, end of psychotherapy) and after 12 weeks (T12 - follow-up). We examined differentially methylated CpG sites and regions, as well as pathways analysis in association with the treatment.Results: Main results obtained have shown 110 differentially methylated regions (DMRs) with a significant adjusted p-value area associated with the effects of trauma-focused psychotherapies in the entire cohort. Several annotated genes are related to inflammatory processes and psychiatric disorders, such as LTA, GFI1, ARID5B, TNFSF13, and LST1. Gene enrichment analyses revealed statistically significant processes related to tumour necrosis factor (TNF) receptor and TNF signalling pathway. Stratified analyses by type of trauma-focused psychotherapy showed statistically significant adjusted p-value area in 141 DMRs only for the group of patients receiving EMDR, with annotated genes related to inflammation and psychiatric disorders, including LTA, GFI1, and S100A8. Gene set enrichment analyses in the EMDR group indicated biological processes related to inflammatory response, particularly the TNF signalling pathway.Conclusion: We provide preliminary valuable insights into global DNA methylation changes associated with trauma-focused psychotherapies effects, in particular with EMDR treatment.

KNeXT: a NetworkX-based topologically relevant KEGG parser

Automating the recreation of gene and mixed gene-compound networks from Kyoto Encyclopedia of Genes and Genomes (KEGG) Markup Language (KGML) files is challenging because the data structure does not preserve the independent or loosely connected neighborhoods in which they were originally derived, referred to here as its topological environment. Identical accession numbers may overlap, causing neighborhoods to artificially collapse based on duplicated identifiers. This causes current parsers to create misleading or erroneous graphical representations when mixed gene networks are converted to gene-only networks. To overcome these challenges we created a python-based KEGG NetworkX Topological (KNeXT) parser that allows users to accurately recapitulate genetic networks and mixed networks from KGML map data. The software, archived as a python package index (PyPI) file to ensure broad application, is designed to ingest KGML files through built-in APIs and dynamically create high-fidelity topological representations. The utilization of NetworkX's framework to generate tab-separated files additionally ensures that KNeXT results may be imported into other graph frameworks and maintain programmatic access to the original x-y axis positions to each node in the KEGG pathway. KNeXT is a well-described Python 3 package that allows users to rapidly download and aggregate specific KGML files and recreate KEGG pathways based on a range of user-defined settings. KNeXT is platform-independent, distinctive, and it is not written on top of other Python parsers. Furthermore, KNeXT enables users to parse entire local folders or single files through command line scripts and convert the output into NCBI or UniProt IDs. KNeXT provides an ability for researchers to generate pathway visualizations while persevering the original context of a KEGG pathway. Source code is freely available at https://github.com/everest-castaneda/knext.

Protocol for generation of a time-resolved cellular interactome during tissue remodeling in adult mice

Efficient skeletal muscle regeneration necessitates fine-tuned coordination among multiple cell types through an intricate network of intercellular communication. We present a protocol for generation of a time-resolved cellular interactome during tissue remodeling. We describe steps for isolating distinct cell populations from skeletal muscle of adult mice after acute damage and extracting RNA from purified cells prior to the generation of RNA sequencing data. We then detail procedures for generating and deciphering a time- and lineage-resolved model of intercellular crosstalk. For complete details on the use and execution of this protocol, please refer to Groppa et al. (2023).1.

Temperature-driven coordination of circadian transcriptome regulation

The circadian rhythm is an evolutionarily-conserved molecular oscillator that enables species to anticipate rhythmic changes in their environment. At a molecular level, the core clock genes induce a circadian oscillation in thousands of genes in a tissue-specific manner, orchestrating myriad biological processes. While studies have investigated how the core clock circuit responds to environmental perturbations such as temperature, the downstream effects of such perturbations on circadian regulation remain poorly understood. By analyzing bulk-RNA sequencing of Drosophila fat bodies harvested from flies subjected to different environmental conditions, we demonstrate a highly condition-specific circadian transcriptome. Further employing a reference-based gene regulatory network (Reactome), we find evidence of increased gene-gene coordination at low temperatures and synchronization of rhythmic genes that are network neighbors. Our results point to the mechanisms by which the circadian clock mediates the fly's response to seasonal changes in temperature.

Transcriptional Profiling of Rat Prefrontal Cortex after Acute Inescapable Footshock Stress

Stress is a primary risk factor for psychiatric disorders such as Major Depressive Disorder (MDD) and Post Traumatic Stress Disorder (PTSD). The response to stress involves the regulation of transcriptional programs, which is supposed to play a role in coping with stress. To evaluate transcriptional processes implemented after exposure to unavoidable traumatic stress, we applied microarray expression analysis to the PFC of rats exposed to acute footshock (FS) stress that were sacrificed immediately after the 40 min session or 2 h or 24 h after. While no substantial changes were observed at the single gene level immediately after the stress session, gene set enrichment analysis showed alterations in neuronal pathways associated with glia development, glia-neuron networking, and synaptic function. Furthermore, we found alterations in the expression of gene sets regulated by specific transcription factors that could represent master regulators of the acute stress response. Of note, these pathways and transcriptional programs are activated during the early stress response (immediately after FS) and are already turned off after 2 h-while at 24 h, the transcriptional profile is largely unaffected. Overall, our analysis provided a transcriptional landscape of the early changes triggered by acute unavoidable FS stress in the PFC of rats, suggesting that the transcriptional wave is fast and mild, but probably enough to activate a cellular response to acute stress.

Spatial compartmentalization of signaling imparts source-specific functions on secreted factors

Efficient regeneration requires multiple cell types acting in coordination. To better understand the intercellular networks involved and how they change when regeneration fails, we profile the transcriptome of hematopoietic, stromal, myogenic, and endothelial cells over 14 days following acute muscle damage. We generate a time-resolved computational model of interactions and identify VEGFA-driven endothelial engagement as a key differentiating feature in models of successful and failed regeneration. In addition, the analysis highlights that the majority of secreted signals, including VEGFA, are simultaneously produced by multiple cell types. To test whether the cellular source of a factor determines its function, we delete VEGFA from two cell types residing in close proximity: stromal and myogenic progenitors. By comparing responses to different types of damage, we find that myogenic and stromal VEGFA have distinct functions in regeneration. This suggests that spatial compartmentalization of signaling plays a key role in intercellular communication networks.

Identifying Network Biomarkers for Alzheimer's Disease Using Single-Cell RNA Sequencing Data

System-level network-based approaches are an emerging field in the biomedical domain since biological networks can be used to analyze complicated biological processes and complex human disorders more efficiently. Network biomarkers are groups of interconnected molecular components causing perturbations in the entire network topology that can be used as indicators of pathogenic biological processes when studying a given disease. Although in the last years computational systems-based approaches have gained ground on the path to discovering new network biomarkers, in complex diseases like Alzheimer's disease (AD), this approach has still much to offer. Especially the adoption of single-cell RNA sequencing (scRNA-seq) has now become the dominant technology for the study of stochastic gene expression. Toward this orientation, we propose an R workflow that extracts disease-perturbed subpathways within a pathway network. We construct a gene-gene interaction network integrated with scRNA-seq expression profiles, and after network processing and pruning, the most active subnetworks are isolated from the entire network topology. The proposed methodology was applied on a real AD-based scRNA-seq data, providing already existing and new potential AD biomarkers in gene network context.

Klebsiella pneumoniae induces host metabolic stress that promotes tolerance to pulmonary infection

K. pneumoniae sequence type 258 (Kp ST258) is a major cause of healthcare-associated pneumonia. However, it remains unclear how it causes protracted courses of infection in spite of its expression of immunostimulatory lipopolysaccharide, which should activate a brisk inflammatory response and bacterial clearance. We predicted that the metabolic stress induced by the bacteria in the host cells shapes an immune response that tolerates infection. We combined in situ metabolic imaging and transcriptional analyses to demonstrate that Kp ST258 activates host glutaminolysis and fatty acid oxidation. This response creates an oxidant-rich microenvironment conducive to the accumulation of anti-inflammatory myeloid cells. In this setting, metabolically active Kp ST258 elicits a disease-tolerant immune response. The bacteria, in turn, adapt to airway oxidants by upregulating the type VI secretion system, which is highly conserved across ST258 strains worldwide. Thus, much of the global success of Kp ST258 in hospital settings can be explained by the metabolic activity provoked in the host that promotes disease tolerance.

MyoData: An expression knowledgebase at single cell/nucleus level for the discovery of coding-noncoding RNA functional interactions in skeletal muscle

•

Regulation of gene expression through non-coding RNAs at single myofiber and nucleus resolution.

•

Reinterpretation of KEGG pathways with microRNA and long non-coding RNA activities.

•

miR-149, -214, and let-7e alter mitochondrial shape.

•

The long non-coding RNA Pvt1 is a sponge for miR-27a.

•

miR-208b regulates Sox6; miR-214 regulates both Sox6 and Slc16a3.

Non-coding RNAs represent the largest part of transcribed mammalian genomes and prevalently exert regulatory functions. Long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) can modulate the activity of each other. Skeletal muscle is the most abundant tissue in mammals. It is composed of different cell types with myofibers that represent the smallest complete contractile system. Considering that lncRNAs and miRNAs are more cell type-specific than coding RNAs, to understand their function it is imperative to evaluate their expression and action within single myofibers. In this database, we collected gene expression data for coding and non-coding genes in single myofibers and used them to produce interaction networks based on expression correlations. Since biological pathways are more informative than networks based on gene expression correlation, to understand how altered genes participate in the studied phenotype, we integrated KEGG pathways with miRNAs and lncRNAs. The database also integrates single nucleus gene expression data on skeletal muscle in different patho-physiological conditions. We demonstrated that these networks can serve as a framework from which to dissect new miRNA and lncRNA functions to experimentally validate. Some interactions included in the database have been previously experimentally validated using high throughput methods. These can be the basis for further functional studies. Using database information, we demonstrate the involvement of miR-149, -214 and let-7e in mitochondria shaping; the ability of the lncRNA Pvt1 to mitigate the action of miR-27a via sponging; and the regulatory activity of miR-214 on Sox6 and Slc16a3. The MyoData is available at https://myodata.bio.unipd.it.

DysPIA: A Novel Dysregulated Pathway Identification Analysis Method

Differential co-expression-based pathway analysis is still limited and not widely used. In most current methods, the pathways were considered as gene sets, but the gene regulation relationships were not considered, and the computational speed was slow. In this article, we proposed a novel Dysregulated Pathway Identification Analysis (DysPIA) method to overcome these shortcomings. We adopted the idea of Correlation by Individual Level Product into analysis and performed a fast enrichment analysis. We constructed a combined gene-pair background which was much more sufficient than the background used in Edge Set Enrichment Analysis. In simulation study, DysPIA was able to identify the causal pathways with high AUC (0.9584 to 0.9896). In p53 mutation data, DysPIA obtained better performance than other methods. It obtained more potential dysregulated pathways that could be literature verified, and it ran much faster (∼1,700-8,000 times faster than other methods when 10,000 permutations). DysPIA was also applied to breast cancer relapse dataset and breast cancer subtype dataset. The results show that DysPIA is effective and has a great biological significance. R packages "DysPIA" and "DysPIAData" are constructed and freely available on R CRAN (https://cran.r-project.org/web/packages/DysPIA/index.html and https://cran.r-project.org/web/packages/DysPIAData/index.html), and on GitHub (https://github.com/lemonwang2020).

A computational strategy for metabolic network construction based on the overlapping ratio: Study of patients' metabolic responses to different dialysis patterns

BACKGROUND

Uremia is a worldwide epidemic disease and poses a serious threat to human health. Both maintenance hemodialysis (HD) and maintenance high flux hemodialysis (HFD) are common treatments for uremia and are generally used in clinical applications. In-depth exploration of patients' metabolic responses to different dialysis patterns can facilitate the understanding of pathological alterations associated with uremia and the effects of different dialysis methods on uremia, which may be used for future personalized therapy. However, due to variations of multiple factors (i.e., genetic, epigenetic and environment) in the process of disease treatments, identification of the similarities and differences in plasma metabolite changes in uremic patients in response to HD and HFD remains challenging.

METHODS

In this study, a computational strategy for metabolic network construction based on the overlapping ratio (MNC-OR) was proposed for disease treatment effect research. In MNC-OR, the overlapping ratio was introduced to measure metabolic reactions and to construct metabolic networks for analysis of different treatment options. Then, MNC-OR was employed to analyze HD-pattern-dependent changes in plasma metabolites to explore the pathological alterations associated with uremia and the effectiveness of different dialysis patterns (i.e., HD and HFD) on uremia. Based on the networks constructed by MNC-OR, two network analysis techniques, namely, similarity analysis and difference analysis of network topology, were used to find the similarity and differences in metabolic signals in patients under treatment with either HD or HFD, which can facilitate the understanding of pathological alterations associated with uremia and provide the guidance for personalized dialysis therapy.

RESULTS

Similarity analysis of network topology suggested that abnormal energy metabolism, gut metabolism and pyrimidine metabolism might occur in uremic patients, and maintenance of both HFD and HD therapies have beneficial effects on uremia. Then, difference analysis of network topology was employed to extract the crucial information related to HD-pattern-dependent changes in plasma metabolites. Experimental results indicated that the amino acid metabolism was closer to the normal status in HFD-treated patients; however, in HD-treated patients, the ability of antioxidation showed greater reduction, and the protein O-GlcNAcylation level was higher. Our findings demonstrate the potential of MNC-OR for explaining the metabolic similarities and differences of patients in response to different dialysis methods, thereby contributing to the guidance of personalized dialysis therapy.

Differential metabolic network construction for personalized medicine: Study of type 2 diabetes mellitus patients' response to gliclazide-modified-release-treated

Individual variation in genetic and environmental factors can cause the differences in metabolic phenotypes, which may have an effect on drug responses of patients. Deep exploration of patients' responses to therapeutic agents is a crucial and urgent event in the personalized treatment study. Using machine learning methods for the discovery of suitability evaluation biomarkers can provide deep insight into the mechanism of disease therapy and facilitate the development of personalized medicine. To find important metabolic network signals for the prediction of patients' drug responses, a novel method referred to as differential metabolic network construction (DMNC) was proposed. In DMNC, concentration changes in metabolite ratios between different pathological states are measured to construct differential metabolic networks, which can be used to advance clinical decision-making. In this study, DMNC was applied to characterize type 2 diabetes mellitus (T2DM) patients' responses against gliclazide modified-release (MR) therapy. Two T2DM metabolomics datasets from different batches of subjects treated by gliclazide MR were analyzed in depth. A network biomarker was defined to assess the patients' suitability for gliclazide MR. It can be effective in the prediction of significant responders from nonsignificant responders, achieving area under the curve values of 0.893 and 1.000 for the discovery and validation sets, respectively. Compared with the metabolites selected by the other methods, the network biomarker selected by DMNC was more stable and precise to reflect the metabolic responses in patients to gliclazide MR therapy, thereby contributing for the personalized medicine of T2DM patients. The better performance of DMNC validated its potential for the identification of network biomarkers to characterize the responses against therapeutic treatments and provide valuable information for personalized medicine.

Computational analysis to repurpose drugs for COVID-19 based on transcriptional response of host cells to SARS-CoV-2

BACKGROUND

The Coronavirus Disease 2019 (COVID-19) pandemic has infected over 10 million people globally with a relatively high mortality rate. There are many therapeutics undergoing clinical trials, but there is no effective vaccine or therapy for treatment thus far. After affected by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), molecular signaling pathways of host cells play critical roles during the life cycle of SARS-CoV-2. Thus, it is significant to identify the involved molecular signaling pathways within the host cells. Drugs targeting these molecular signaling pathways could be potentially effective for COVID-19 treatment.

METHODS

In this study, we developed a novel integrative analysis approach to identify the related molecular signaling pathways within host cells, and repurposed drugs as potentially effective treatments for COVID-19, based on the transcriptional response of host cells.

RESULTS

We identified activated signaling pathways associated with the infection caused SARS-CoV-2 in human lung epithelial cells through integrative analysis. Then, the activated gene ontologies (GOs) and super GOs were identified. Signaling pathways and GOs such as MAPK, JNK, STAT, ERK, JAK-STAT, IRF7-NFkB signaling, and MYD88/CXCR6 immune signaling were particularly activated. Based on the identified signaling pathways and GOs, a set of potentially effective drugs were repurposed by integrating the drug-target and reverse gene expression data resources. In addition to many drugs being evaluated in clinical trials, the dexamethasone was top-ranked in the prediction, which was the first reported drug to be able to significantly reduce the death rate of COVID-19 patients receiving respiratory support.

CONCLUSIONS

The integrative genomics data analysis and results can be helpful to understand the associated molecular signaling pathways within host cells, and facilitate the discovery of effective drugs for COVID-19 treatment.

A strategy to incorporate prior knowledge into correlation network cutoff selection

Correlation networks are frequently used to statistically extract biological interactions between omics markers. Network edge selection is typically based on the statistical significance of the correlation coefficients. This procedure, however, is not guaranteed to capture biological mechanisms. We here propose an alternative approach for network reconstruction: a cutoff selection algorithm that maximizes the overlap of the inferred network with available prior knowledge. We first evaluate the approach on IgG glycomics data, for which the biochemical pathway is known and well-characterized. Importantly, even in the case of incomplete or incorrect prior knowledge, the optimal network is close to the true optimum. We then demonstrate the generalizability of the approach with applications to untargeted metabolomics and transcriptomics data. For the transcriptomics case, we demonstrate that the optimized network is superior to statistical networks in systematically retrieving interactions that were not included in the biological reference used for optimization.

Identification of differentially expressed subpathways via a bilevel consensus scoring of network topology and gene expression

Identifying differentially expressed subpathways connected to the emergence of a disease that can be considered as candidates for pharmacological intervention, with minimal off-target effects, is a daunting task. In this direction, we present a bilevel subpathway analysis method to identify differentially expressed subpathways that are connected with an experimental condition, while taking into account potential crosstalks between subpathways which arise due to their connectivity in a combined multi-pathway network. The efficacy of the method is demonstrated on a hematopoietic stem cell aging dataset, with findings corroborated using recent literature.

MOSClip: multi-omic and survival pathway analysis for the identification of survival associated gene and modules

Survival analyses of gene expression data has been a useful and widely used approach in clinical applications. But, in complex diseases, such as cancer, the identification of survival-associated cell processes - rather than single genes - provides more informative results because the efficacy of survival prediction increases when multiple prognostic features are combined to enlarge the possibility of having druggable targets. Moreover, genome-wide screening in molecular medicine has rapidly grown, providing not only gene expression but also multi-omic measurements such as DNA mutations, methylation, expression, and copy number data. In cancer, virtually all these aberrations can contribute in synergy to pathological processes, and their measurements can improve a patient’s outcome and help in diagnosis and treatment decisions. Here, we present MOSClip, an R package implementing a new topological pathway analysis tool able to integrate multi-omic data and look for survival-associated gene modules. MOSClip tests the survival association of dimensionality-reduced multi-omic data using multivariate models, providing graphical devices for management, browsing and interpretation of results. Using simulated data we evaluated MOSClip performance in terms of false positives and false negatives in different settings, while the TCGA ovarian cancer dataset is used as a case study to highlight MOSClip’s potential.

Transcriptional Characterization of Stage I Epithelial Ovarian Cancer: A Multicentric Study

Stage I epithelial ovarian cancer (EOC) represents about 10% of all EOCs. It is characterized by a complex histopathological and molecular heterogeneity, and it is composed of five main histological subtypes (mucinous, endometrioid, clear cell and high, and low grade serous), which have peculiar genetic, molecular, and clinical characteristics. As it occurs less frequently than advanced-stage EOC, its molecular features have not been thoroughly investigated. In this study, using in silico approaches and gene expression data, on a multicentric cohort composed of 208 snap-frozen tumor biopsies, we explored the subtype-specific molecular alterations that regulate tumor aggressiveness in stage I EOC. We found that single genes rather than pathways are responsible for histotype specificities and that a cAMP-PKA-CREB1 signaling axis seems to play a central role in histotype differentiation. Moreover, our results indicate that immune response seems to be, at least in part, involved in histotype differences, as a higher immune-reactive behavior of serous and mucinous samples was observed with respect to other histotypes.

The Impact of Pathway Database Choice on Statistical Enrichment Analysis and Predictive Modeling

Pathway-centric approaches are widely used to interpret and contextualize -omics data. However, databases contain different representations of the same biological pathway, which may lead to different results of statistical enrichment analysis and predictive models in the context of precision medicine. We have performed an in-depth benchmarking of the impact of pathway database choice on statistical enrichment analysis and predictive modeling. We analyzed five cancer datasets using three major pathway databases and developed an approach to merge several databases into a single integrative one: MPath. Our results show that equivalent pathways from different databases yield disparate results in statistical enrichment analysis. Moreover, we observed a significant dataset-dependent impact on the performance of machine learning models on different prediction tasks. In some cases, MPath significantly improved prediction performance and also reduced the variance of prediction performances. Furthermore, MPath yielded more consistent and biologically plausible results in statistical enrichment analyses. In summary, this benchmarking study demonstrates that pathway database choice can influence the results of statistical enrichment analysis and predictive modeling. Therefore, we recommend the use of multiple pathway databases or integrative ones.

SourceSet: A graphical model approach to identify primary genes in perturbed biological pathways

Topological gene-set analysis has emerged as a powerful means for omic data interpretation. Although numerous methods for identifying dysregulated genes have been proposed, few of them aim to distinguish genes that are the real source of perturbation from those that merely respond to the signal dysregulation. Here, we propose a new method, called SourceSet, able to distinguish between the primary and the secondary dysregulation within a Gaussian graphical model context. The proposed method compares gene expression profiles in the control and in the perturbed condition and detects the differences in both the mean and the covariance parameters with a series of likelihood ratio tests. The resulting evidence is used to infer the primary and the secondary set, i.e. the genes responsible for the primary dysregulation, and the genes affected by the perturbation through network propagation. The proposed method demonstrates high specificity and sensitivity in different simulated scenarios and on several real biological case studies. In order to fit into the more traditional pathway analysis framework, SourceSet R package also extends the analysis from a single to multiple pathways and provides several graphical outputs, including Cytoscape visualization to browse the results.

The rapid increase in omic studies has created a need to understand the biological implications of their results. Gene-set analysis has emerged as a powerful means for gaining such understanding, evolving in the last decade from the classical enrichment analysis to the more powerful topological approaches. Although numerous methods for identifying dysregulated genes have been proposed, few of them aim to distinguish genes that are the real source of perturbation from those that merely respond to the signal dysregulation. This distinction is crucial for network medicine, where the prioritization of the effect of biological perturbations may help in the molecular understanding of drug treatments and diseases. Here we propose a new method, called SourceSet, able to distinguish between primary and secondary dysregulation within a graphical model context, demonstrating a high specificity and sensitivity in different simulated scenarios and on real biological case studies.

The metaRbolomics Toolbox in Bioconductor and beyond

Metabolomics aims to measure and characterise the complex composition of metabolites in a biological system. Metabolomics studies involve sophisticated analytical techniques such as mass spectrometry and nuclear magnetic resonance spectroscopy, and generate large amounts of high-dimensional and complex experimental data. Open source processing and analysis tools are of major interest in light of innovative, open and reproducible science. The scientific community has developed a wide range of open source software, providing freely available advanced processing and analysis approaches. The programming and statistics environment R has emerged as one of the most popular environments to process and analyse Metabolomics datasets. A major benefit of such an environment is the possibility of connecting different tools into more complex workflows. Combining reusable data processing R scripts with the experimental data thus allows for open, reproducible research. This review provides an extensive overview of existing packages in R for different steps in a typical computational metabolomics workflow, including data processing, biostatistics, metabolite annotation and identification, and biochemical network and pathway analysis. Multifunctional workflows, possible user interfaces and integration into workflow management systems are also reviewed. In total, this review summarises more than two hundred metabolomics specific packages primarily available on CRAN, Bioconductor and GitHub.

PathMe: merging and exploring mechanistic pathway knowledge

Background

The complexity of representing biological systems is compounded by an ever-expanding body of knowledge emerging from multi-omics experiments. A number of pathway databases have facilitated pathway-centric approaches that assist in the interpretation of molecular signatures yielded by these experiments. However, the lack of interoperability between pathway databases has hindered the ability to harmonize these resources and to exploit their consolidated knowledge. Such a unification of pathway knowledge is imperative in enhancing the comprehension and modeling of biological abstractions.

Results

Here, we present PathMe, a Python package that transforms pathway knowledge from three major pathway databases into a unified abstraction using Biological Expression Language as the pivotal, integrative schema. PathMe is complemented by a novel web application (freely available at https://pathme.scai.fraunhofer.de/) which allows users to comprehensively explore pathway crosstalk and compare areas of consensus and discrepancies.

Conclusions

This work has harmonized three major pathway databases and transformed them into a unified schema in order to gain a holistic picture of pathway knowledge. We demonstrate the utility of the PathMe framework in: i) integrating pathway landscapes at the database level, ii) comparing the degree of consensus at the pathway level, and iii) exploring pathway crosstalk and investigating consensus at the molecular level.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2863-9) contains supplementary material, which is available to authorized users.

The Impact of Pathway Database Choice on Statistical Enrichment Analysis and Predictive Modeling

Pathway-centric approaches are widely used to interpret and contextualize -omics data. However, databases contain different representations of the same biological pathway, which may lead to different results of statistical enrichment analysis and predictive models in the context of precision medicine. We have performed an in-depth benchmarking of the impact of pathway database choice on statistical enrichment analysis and predictive modeling. We analyzed five cancer datasets using three major pathway databases and developed an approach to merge several databases into a single integrative one: MPath. Our results show that equivalent pathways from different databases yield disparate results in statistical enrichment analysis. Moreover, we observed a significant dataset-dependent impact on the performance of machine learning models on different prediction tasks. In some cases, MPath significantly improved prediction performance and also reduced the variance of prediction performances. Furthermore, MPath yielded more consistent and biologically plausible results in statistical enrichment analyses. In summary, this benchmarking study demonstrates that pathway database choice can influence the results of statistical enrichment analysis and predictive modeling. Therefore, we recommend the use of multiple pathway databases or integrative ones.