Geometric coherence of single-cell CRISPR perturbations reveals regulatory architecture and predicts cellular stress

TL;DR

The paper introduces a geometric stability metric, Shesha, to quantify the directional coherence of single-cell CRISPR perturbation responses, revealing insights into regulatory architecture and cellular stress.

Contribution

It presents a novel geometric stability metric, Shesha, that uncovers regulatory mechanisms and predicts cellular stress by analyzing perturbation coherence across multiple datasets.

Findings

Stability correlates strongly with effect magnitude (Spearman ρ=0.75-0.97).

Discordant cases reveal regulatory architecture, distinguishing master regulators from lineage-specific factors.

High stability associates with lower chaperone activation and cellular stress.

Abstract

Genome engineering has achieved remarkable sequence-level precision, yet predicting the transcriptomic state that a cell will occupy after perturbation remains an open problem. Single-cell CRISPR screens measure how far cells move from their unperturbed state, but this effect magnitude ignores a fundamental question: do the cells move together? Two perturbations with identical magnitude can produce qualitatively different outcomes if one drives cells coherently along a shared trajectory while the other scatters them across expression space. We introduce a geometric stability metric, Shesha, that quantifies the directional coherence of single-cell perturbation responses as the mean cosine similarity between individual cell shift vectors and the mean perturbation direction. Across five CRISPR datasets (2,200+ perturbations spanning CRISPRa, CRISPRi, and pooled screens), stability correlates strongly with effect magnitude (Spearman $\rho=0.75-0.97$), with a calibrated cross-dataset correlation of 0.97. Crucially, discordant cases where the two metrics decouple expose regulatory architecture: pleiotropic master regulators such as CEBPA and GATA1 pay a "geometric tax," producing large but incoherent shifts, while lineage-specific factors such as KLF1 produce tightly coordinated responses. After controlling for magnitude, geometric instability is independently associated with elevated chaperone activation (HSPA5/BiP; $\rho_{partial}=-0.34$ and $-0.21$ across datasets), and the high-stability/high-stress quadrant is systematically depleted. The magnitude-stability relationship persists in scGPT foundation model embeddings, confirming it is a property of biological state space rather than linear projection. Perturbation stability provides a complementary axis for hit prioritization in screens, phenotypic quality control in cell manufacturing, and evaluation of in silico perturbation predictions.