Locat: Joint enrichment and depletion testing identifies localized marker genes in single-cell transcriptomics

Abstract

Several methods identify marker genes that delineate cell populations in single-cell transcriptomic data, yet most emphasize enrichment within candidate populations without testing whether expression is significantly reduced elsewhere. We present Locat, a framework for identifying highly specific localized genes by testing whether expression is concentrated within compact regions of a cellular embedding and depleted outside them. For each gene, Locat fits weighted Gaussian mixture models to gene-specific and background densities, computes concentration and depletion statistics, and integrates them into a unified localization score. Across synthetic benchmarks with controlled ground truth, Locat detects uni-modal, multi-modal, and sparse localized patterns and loses significance when expression becomes indistinguishable from background structure. In developmental, perturbation, and differentiation datasets, Locat identifies compact marker sets that capture lineage organization, condition-specific programs, and temporal dynamics. These sets are often smaller than highly variable gene selections, while embeddings built from them preserve major cell populations and developmental programs in several cases. In murine dermis, interferon-treated PBMCs, and retinoic acid-induced embryonic stem cell differentiation, localized genes recover differentiation trajectories, stimulus-responsive programs, and reproducible stage-specific patterns. Together, these results show that jointly assessing concentration and depletion yields specific, interpretable marker genes.