SEMgsa: topology-based pathway enrichment analysis with structural equation models

Background

Pathway enrichment analysis is extensively used in high-throughput experimental studies to gain insight into the functional roles of pre-defined subsets of genes, proteins and metabolites. Methods that leverages information on the topology of the underlying pathways outperform simpler methods that only consider pathway membership, leading to improved performance. Among all the proposed software tools, there’s the need to combine high statistical power together with a user-friendly framework, making it difficult to choose the best method for a particular experimental environment.

Results

We propose SEMgsa, a topology-based algorithm developed into the framework of structural equation models. SEMgsa combine the SEM p values regarding node-specific group effect estimates in terms of activation or inhibition, after statistically controlling biological relations among genes within pathways. We used SEMgsa to identify biologically relevant results in a Coronavirus disease (COVID-19) RNA-seq dataset (GEO accession: GSE172114) together with a frontotemporal dementia (FTD) DNA methylation dataset (GEO accession: GSE53740) and compared its performance with some existing methods. SEMgsa is highly sensitive to the pathways designed for the specific disease, showing low p values (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$< 0.001$$\end{document}<0.001) and ranking in high positions, outperforming existing software tools. Three pathway dysregulation mechanisms were used to generate simulated expression data and evaluate the performance of methods in terms of type I error followed by their statistical power. Simulation results confirm best overall performance of SEMgsa.

Conclusions

SEMgsa is a novel yet powerful method for identifying enrichment with regard to gene expression data. It takes into account topological information and exploits pathway perturbation statistics to reveal biological information. SEMgsa is implemented in the R package SEMgraph, easily available at https://CRAN.R-project.org/package=SEMgraph.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04884-8.

BZD9L1 sirtuin inhibitor: Identification of key molecular targets and their biological functions in HCT 116 colorectal cancer cells

BZD9L1 was previously described as a SIRT1/2 inhibitor with anti-cancer activities in colorectal cancer (CRC), either as a standalone chemotherapy or in combination with 5-fluorouracil. BZD9L1 was reported to induce apoptosis in CRC cells; however, the network of intracellular pathways and crosstalk between molecular players mediated by BZD9L1 is not fully understood. This study aimed to uncover the mechanisms involved in BZD9L1-mediated cytotoxicity based on previous and new findings for the prediction and identification of related pathways and key molecular players. BZD9L1-regulated candidate targets (RCTs) were identified using a range of molecular, cell-based and biochemical techniques on the HCT 116 cell line. BZD9L1 regulated major cancer pathways including Notch, p53, cell cycle, NFκB, Myc/MAX, and MAPK/ERK signalling pathways. BZD9L1 also induced reactive oxygen species (ROS), regulated apoptosis-related proteins, and altered cell polarity and adhesion profiles. In silico analyses revealed that most RCTs were interconnected, and were involved in the modulation of catalytic activity, metabolism and transcription regulation, response to cytokines, and apoptosis signalling pathways. These RCTs were implicated in p53-dependent apoptosis pathway. This study provides the first assessment of possible associations of molecular players underlying the cytotoxic activity of BZD9L1, and establishes the links between RCTs and apoptosis through the p53 pathway.

A comparative study of topology-based pathway enrichment analysis methods

BACKGROUND

Pathway enrichment extensively used in the analysis of Omics data for gaining biological insights into the functional roles of pre-defined subsets of genes, proteins and metabolites. A large number of methods have been proposed in the literature for this task. The vast majority of these methods use as input expression levels of the biomolecules under study together with their membership in pathways of interest. The latest generation of pathway enrichment methods also leverages information on the topology of the underlying pathways, which as evidence from their evaluation reveals, lead to improved sensitivity and specificity. Nevertheless, a systematic empirical comparison of such methods is still lacking, making selection of the most suitable method for a specific experimental setting challenging. This comparative study of nine network-based methods for pathway enrichment analysis aims to provide a systematic evaluation of their performance based on three real data sets with different number of features (genes/metabolites) and number of samples.

RESULTS

The findings highlight both methodological and empirical differences across the nine methods. In particular, certain methods assess pathway enrichment due to differences both across expression levels and in the strength of the interconnectedness of the members of the pathway, while others only leverage differential expression levels. In the more challenging setting involving a metabolomics data set, the results show that methods that utilize both pieces of information (with NetGSA being a prototypical one) exhibit superior statistical power in detecting pathway enrichment.

CONCLUSION

The analysis reveals that a number of methods perform equally well when testing large size pathways, which is the case with genomic data. On the other hand, NetGSA that takes into consideration both differential expression of the biomolecules in the pathway, as well as changes in the topology exhibits a superior performance when testing small size pathways, which is usually the case for metabolomics data.