Abstract
Single nucleus RNA sequencing (snRNA-seq), an alternative to single cell RNA sequencing (scRNA-seq), encounters technical challenges in obtaining high-quality nuclei and RNA, persistently hindering its applications. Here, we present a robust technique for isolating nuclei across various tissue types, remarkably enhancing snRNA-seq data quality. Employing this approach, we comprehensively characterize the depot-dependent cellular dynamics of various cell types underlying adipose tissue remodeling during obesity. By integrating bulk nuclear RNA-seq from adipocyte nuclei of different sizes, we identify distinct adipocyte subpopulations categorized by size and functionality. These subpopulations follow two divergent trajectories, adaptive and pathological, with their prevalence varying by depot. Specifically, we identify a key molecular feature of dysfunctional hypertrophic adipocytes, a global shutdown in gene expression, along with elevated stress and inflammatory responses. Furthermore, our differential gene expression analysis reveals distinct contributions of adipocyte subpopulations to the overall pathophysiology of adipose tissue. Our study establishes a robust snRNA-seq method, providing novel insights into the biological processes involved in adipose tissue remodeling during obesity, with broader applicability across diverse biological systems.