By performing snATAC-seq in a rat model of NAFLD and employing novel statistical methods, we elucidated a global picture of in vivo TF regulation in each cell type as a set of modules and discovered core genes for NAFLD-relevant biological processes. This evidence concerns the gene TF and metabolic dysfunction-associated steatotic liver disease.