Accurate estimation of pathway activity in single cells for clustering and differential analysis

Inferring which and how biological pathways and gene sets change is a key question in many studies that utilize single-cell RNA sequencing. Typically, these questions are addressed by quantifying the enrichment of known gene sets in lists of genes derived from global analysis. Here we offer SiPSiC, a new method to infer pathway activity in every single cell. This allows more sensitive differential analysis and utilization of pathway scores to cluster cells and compute UMAP or other similar projections. We apply our method to COVID-19, lung adenocarcinoma and glioma data sets, and demonstrate its utility. SiPSiC analysis results are consistent with findings reported in previous studies in many cases, but SiPSiC also reveals the differential activity of novel pathways, enabling us to suggest new mechanisms underlying the pathophysiology of these diseases and demonstrating SiPSiC's high accuracy and sensitivity in detecting biological function and traits. In addition, we demonstrate how it can be used to better classify cells based on activity of biological pathways instead of single genes and its ability to overcome patient-specific artifacts.