Parts-based decomposition of spatial genomics data finds distinct tissue regions

Dimension reduction is a cornerstone of exploratory data analysis; however, traditional methods fail to preserve the spatial context of spatial genomics data. In this work, we develop a nonnegative spatial factorization (NSF) model that allows interpretable, parts-based decomposition of spatial single-cell count data. NSF allows label-free annotation of regions of interest in spatial genomics data and identifies genes and cells that can be used to define those regions.