Integrating spatial transcriptomics and single-cell RNA-seq dissects immune microenvironment in fatty liver regeneration

Abstract

Dear Editor,

Liver zonal regeneration in healthy states involves a dynamic interplay between parenchymal and nonparenchymal cells, whereas fatty liver chronicity disrupts immune-niche coordination, altering intercellular crosstalk.^{1, 2} While single-cell technologies resolve cellular heterogeneity, they often overlook spatial regulation of cellular functions. We performed PHx (partial hepatectomy) on healthy mice and mice on a high-fat diet, and sampling was performed on postoperative days 0, 2, 4, and 6. Our study integrates spatial transcriptomics with single-cell profiling, bulk RNA-seq, and smFISH to construct a spatiotemporal atlas of liver regeneration post-PHx.

We constructed single-cell maps (105 442 cells after quality control) and spatial maps (25 995 points) depicting the dynamics of liver regeneration after hepatectomy in normal and fatty livers (Figure 1A,B). Batch-effect-free integration (Figure S1A–C) confirmed impaired ecological niche coordination in the fatty microenvironment. Clustering identified 14 major cell types (annotated by literature calibration markers; Figure 1C; Table S1), in which Kupffer/endothelial cells—key regenerative mediators—were reduced in the fatty liver state and parenchymal cells (hepatocytes/BECs) were diminished (Figure 1D; Figure S1B).^{3, 4} Spatial validation confirmed lipid-laden in the fatty model (Figure 1E). The temporal analysis highlighted a different immune response: fatty livers exhibited an early neutrophil/monocyte surge at day 2 (Figure S1D), in contrast to the 48 h regeneration peak in normal livers.⁵ These data establish a spatiotemporal map that identifies fatty liver-specific defects in parenchymal-immune crosstalk and delayed regenerative activation.

Meanwhile, spatial mapping validated the classical 1, 2, and 3 zonation markers (Sds/Igfbp2/Oat; Figure 1E–G) and identified 7 molecular niches after PHx (Figure S1E–H).^6-8 To further characterize the features of each molecular niche, we performed differential analysis to identify the representative markers (Table S2). Niche 2 and niche 1/6 expressed Cyp2f2/Sds and Cyp2e1/Oat, respectively, whereas midlobular niche 3/4 was enriched in Igfbp2.⁶ The molecular niches might be representative of the spatial structure (Figure 2A,C). Spatial projection confirmed niche-structure alignment: the zonal structure remains clear and significant in different disease states and at different points in time (Figure 2B). Pathway analysis showed attenuated Wnt activation (Figure 2E,F) and similarly attenuated TGF-β inhibition (Figure S2A,B) in fatty liver compared with normal liver.^{9, 10} Single-cell profiling showed that regeneration of hepatocytes in the fatty liver was delayed—Wnt signalling and proliferative activity peaked on day 4, whereas in the normal group, it peaked on day 2 (Figure 2D,G; Figure S2C,D). This temporal variation was also confirmed by bulk RNA-seq and western blot results of cell cycle markers such as Pcna (Figure 2H,I). These data establish that fatty liver retains regenerative capacity but with impaired spatiotemporal coordination, evidenced by disorganized niche architecture, desynchronized Wnt/TGF-β signalling, and delayed proliferation kinetics.

We integrated single-cell data from hepatocytes and cholangiocytes and noted a specific subpopulation, Vim⁺Cd44⁺ interface hepatocytes (Figure 3A). In the fatty state, interface hepatocytes had an enhanced proliferative capacity compared with other hepatocytes. However, Hippo/Wnt signalling remained strong in the normal state (Figure 3B,C). Spatial mapping with smFISH showed significant enrichment of Vim⁺Cd44⁺ interface hepatocytes in fatty liver at the peak of regeneration (day 4, Figures 3D, S3A–D). These interface cells showed regional redistribution: at the peak of normal regeneration (day 2) they were enriched in zone 3, whereas at the peak of fatty state regeneration (day 4) they were enriched in zone 1 (Figure 3E). Fatty liver is less regenerative, and Vim⁺Cd44⁺ interface hepatocytes may be a compensatory response to stress, acting as a transient intermediate to compensate for impaired regenerative signalling.

Single-cell analysis identified 12 T cell subsets and characterized the cytotoxic, immunosuppressive, and inflammatory capacities of the T subpopulations (Figure 4A–C), with Foxp3⁺Ctla4⁺CD4⁺ Tregs and Ccl5⁺CD8⁺ T cells exhibiting dynamic regulation during regeneration. In normal liver, both subsets decreased at day 2 (pro-regeneration phase) before rebounding post-repair, while fatty liver displayed dysregulation with early CD8⁺ T cell accumulation and insufficient Treg expansion (Figure 4D). Compared with the other timepoint, the Cd4⁺ Tregs specifically re-located itself to Zone 3 to conduct the immunosuppressive function, while Ccl5⁺ Cd8⁺ T seemed to be relatively sparse over Zone 2 and Zone 3 to harness the regeneration at both regenerations’ turning point (Day 2 in normal status, Day 4 in fatty status, Figure 4E, Figure S3E–G).

Moreover, the single-cell analysis identified six endothelial subtypes (Figure S4A,B) with zonal specialization. The endothelial cells were well located with Kit⁺ LSCEs and Bmp4⁺ LVECs located in Zone 3, Mup11⁺ LSECs and Ccl21a⁺ LVECs located in Zone 1, while Ccr1⁺ LSECs and Mki67⁺ LSEC located between Zone 1 and Zone 3 (Figure S4C,D). Zone 3 endothelial cells were enhanced in regenerative pathways such as PI3K-Akt and Hippo signalling pathways (Figure S4E). We noticed the immune suppressive function of Mup11⁺ LSCEs showed a similar trend to liver regeneration ability, indicating this LSCE might be another essential endothelial cell in liver regeneration (Figure S4F).

Myeloid profiling identified 13 subsets (Figure S5A,B). We focused on a subgroup termed Ace⁺ monocytes (Figure S5C). This subgroup is more abundant on day 2 in normal status and on day 4 in fatty status, suggesting that it may have some contribution to liver regeneration (Figure S5D). smFISH validated the results (Figure S5F–I). We explore the function of Ace⁺ monocytes showing this subtype is mainly involved in injury repair functions including focal adhesion, platelet activation, and phagocytosis to remove dead cells and avoid immune response (Figure S5E).

In summary, the fatty liver shows a delay in regeneration compared with normal liver after PHx. The intermediate fetal state (Vim⁺Cd44⁺ interface hepatocytes), however, may be an important source of fatty liver regeneration. The microenvironment of fatty liver was reshaped to hamper the liver regeneration with excessive immune function by insufficient Cd4⁺ Tregs and enhanced injury response by early accumulated Ccl5⁺ Cd8⁺ T at day 2. We also identified several subgroups that potentially promote liver regeneration such as Mup11⁺ LSEC located in Zone 1 and Ace⁺ monocyte more enriched at liver regeneration turning point (Figure 4F). Spatial positioning dictates niche-specific cellular functions across parenchymal and nonparenchymal compartments. Our spatiotemporal mapping of fatty liver microenvironments reveals therapeutic targets for cell-type or location-based regenerative strategies.

Chenhao Xu and Renyi Su: Conceptualization, methodology, data curation, formal analysis, software, resources, project administration, visualization, and writing—original draft. Yisu Song: Writing—review & editing. Wenzhi Shu: Data Curation and writing—review & editing. Mengfan Yang: Data curation. Zhe Yang: Conceptualization and supervision. Xiao Xu: Conceptualization, writing—review & editing, funding acquisition, and supervision. Xuyong Wei: Conceptualization and funding acquisition. All authors have approved the submission of the manuscript.

The authors declare no conflict of interest.