Disentangling associations between complex traits and cell types with seismic

Integrating single-cell RNA sequencing with Genome-Wide Association Studies (GWAS) can uncover cell types involved in complex traits and disease. However, current methods often lack scalability, interpretability, and robustness. We present seismic, a framework that computes a novel specificity score capturing both expression magnitude and consistency across cell types and introduces influential gene analysis, an approach to identify genes driving each cell type-trait association. Across over 1000 cell-type characterizations at different granularities and 28 polygenic traits, seismic corroborates known associations and uncovers trait-relevant cell groups not apparent through other methodologies. In Parkinson’s and Alzheimer’s, seismic unveils both cell- and brain-region-specific differences in pathology. Analyzing a pathology-based Alzheimer’s GWAS with seismic enables the identification of vulnerable neuron populations and molecular pathways implicated in their neurodegeneration. In general, seismic is a computationally efficient, powerful, and interpretable approach for mapping the relationships between polygenic traits and cell-type-specific expression, offering new insights into disease mechanisms. Here the authors present seismic, which integrates single-cell RNA-seq and GWAS to uncover trait-relevant cell types and influential genes driving their association. In Alzheimer’s, it distinguishes microglial signals in clinical GWAS versus neuronal vulnerability in tau GWAS.