Search papers, labs, and topics across Lattice.
This paper investigates why causal methods for gene regulatory network (GRN) inference often underperform correlation-based methods on single-cell RNA-seq data, despite theoretical advantages. They introduce a controlled diagnostic framework to isolate and study the impact of seven biologically-motivated pathologies on six representative GRN inference methods. Results demonstrate that causal methods excel in clean regimes but are selectively neutralized by pathologies like dropout and latent confounders, and that error-type decomposition reveals qualitatively different errors among methods with similar aggregate accuracy.
Causal gene regulatory network inference methods only dominate when data is pristine; common single-cell data pathologies like dropout and latent confounders selectively negate their advantages.
Despite theoretical advantages, causal methods for Gene Regulatory Network (GRN) inference from single-cell RNA-seq data consistently fail to match or outperform correlation-based baselines in many realistic benchmarks, a persistent puzzle which casts doubt on the value of causality for this task. We argue that existing benchmarks are insufficiently controlled to answer this question because they evaluate on real or semi-real data where multiple pathologies co-occur, confounding failure modes, and obscuring the specific conditions under which different inference methods excel or fail. To address this gap, we introduce a controlled diagnostic framework that isolates seven biologically motivated pathologies (dropout, latent confounders, cell-type mixing, feedback loops, network density, sample size, and pseudotime drift) and measure how six representative methods spanning three inference paradigms degrade as each pathology intensifies. Across 6,120 controlled experiments, we find that causal methods genuinely dominate in clean and structurally favorable regimes, but specific pathologies (notably dropout and latent confounders) selectively neutralize their advantages. We further introduce an error-type decomposition that reveals methods with similar aggregate accuracy commit qualitatively different errors. To probe whether single-pathology effects persist when multiple stressors co-occur, we perform an interaction sweep over the three most impactful pathologies and find that their joint effects are sub-additive, while also exposing density-conditional cross-overs invisible to single-dial analysis. Our findings offer a nuanced understanding of when and why different methods succeed or fail for GRN inference, providing actionable insights for method development and practical guidance for practitioners.