
==== Front
Int J Biol Sci
Int J Biol Sci
ijbs
International Journal of Biological Sciences
1449-2288
Ivyspring International Publisher Sydney

10.7150/ijbs.98362
ijbsv20p4438
Research Paper
Multiomics analyses decipher intricate changes in the cellular and metabolic landscape of steatotic livers upon dietary restriction and sleeve gastrectomy
Chen Shuai 1†
Zeng Qinghe 23†
Cai Xiurong 4
Xue Jiaming 1
Yin Guo 5
Song Peng 1
Tang Liming 1
Klein Christophe 3
Tacke Frank 5
Guillot Adrien 5✉
Liu Hanyang 15✉
1 Department of General Surgery, The Third Affiliated Hospital of Nanjing Medical University, Changzhou Medical Center, Nanjing Medical University, Changzhou 213000, China.
2 Laboratoire d'Informatique Paris Descartes (LIPADE), Université Paris Cité, Paris 75014, France.
3 Centre d'Histologie, d'Imagerie et de Cytométrie (CHIC), Centre de Recherche des Cordeliers, INSERM, Sorbonne Université, Université de Paris, Paris 75014, France.
4 The Jackson Laboratory, Bar Harbor, ME, USA.
5 Department of Hepatology and Gastroenterology, Charité Universitätsmedizin Berlin, Campus Virchow-Klinikum and Campus Charité Mitte, Berlin 13353, Germany.
✉ Corresponding authors: Hanyang Liu, Department of General Surgery, The Third Affiliated Hospital of Nanjing Medical University, Changzhou Medical Center, Nanjing Medical University, Changzhou 213000, China; Department of Hepatology and Gastroenterology, Charité Universitätsmedizin Berlin, Campus Virchow-Klinikum and Campus Charité Mitte, 13353 Berlin, Germany. Email: brandenliu@live.com. Adrien Guillot, Department of Hepatology and Gastroenterology, Charité Universitätsmedizin Berlin, Campus Virchow-Klinikum and Campus Charité Mitte, 13353 Berlin, Germany. Email: adrien.guillot@charite.de.
† These authors contributed equally.

Competing Interests: FT's lab has received research funding from Gilead, AstraZeneca and MSD (funding to the institution). FT has received honoraria for consulting or lectures from AstraZeneca, Gilead, AbbVie, BMS, Boehringer, Madrigal, Intercept, Falk, Inventiva, MSD, GSK, Orphalan, Merz, Pfizer, Alnylam, CSL Behring, Novo Nordisk, Sanofi and Novartis.

2024
19 8 2024
20 11 44384457
12 5 2024
4 8 2024
© The author(s)
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/). See http://ivyspring.com/terms for full terms and conditions.
Metabolic dysfunction-associated steatotic liver disease (MASLD) is a chronic, progressive liver disease that encompasses a spectrum of steatosis, steatohepatitis (or MASH), and fibrosis. Evidence suggests that dietary restriction (DR) and sleeve gastrectomy (SG) can lead to remission of hepatic steatosis and inflammation through weight loss, but it is unclear whether these procedures induce distinct metabolic or immunological changes in MASLD livers. This study aims to elucidate the intricate hepatic changes following DR, SG or sham surgery in rats fed a high-fat diet as a model of obesity-related MASLD, in comparison to a clinical cohort of patients undergoing SG. Single-cell and single-nuclei transcriptome analysis, spatial metabolomics, and immunohistochemistry revealed the liver landscape, while circulating biomarkers were measured in serum samples. Artificial intelligence (AI)-assisted image analysis characterized the spatial distribution of hepatocytes, myeloid cells and lymphocytes. In patients and experimental MASLD rats, SG improved body mass index, circulating liver injury biomarkers and triglyceride levels. Both DR and SG attenuated liver steatosis and fibrosis in rats. Metabolism-related genes (Ppara, Cyp2e1 and Cyp7a1) were upregulated in hepatocytes upon DR and SG, while SG broadly upregulated lipid metabolism on cholangiocytes, monocytes, macrophages, and neutrophils. Furthermore, SG promoted restorative myeloid cell accumulation in the liver not only ameliorating inflammation but activating liver repair processes. Regions with potent myeloid infiltration were marked with enhanced metabolic capacities upon SG. Additionally, a disruption of periportal hepatocyte functions was observed upon DR. In conclusion, this study indicates a dynamic cellular crosstalk in steatotic livers of patients undergoing SG. Notably, PPARα- and gut-liver axis-related processes, and metabolically active myeloid cell infiltration indicate intervention-related mechanisms supporting the indication of SG for the treatment of MASLD.

MASLD
Bariatric surgery
Immunometabolism
Liver zonation
Myeloid cells
==== Body
pmcIntroduction

Obesity and overweight represent significant global health challenges, affecting approximately 2 billion people worldwide 1. Obesity-related complications encompass a spectrum of metabolic syndromes, including metabolic dysfunction-associated liver diseases (MASLD), metabolic dysfunction-associated steatohepatitis (MASH) and multi-organ disorders, contributing to elevated mortality rates among affected individuals 2, 3. MASLD encompasses a range of chronic liver pathologies affecting more than a quarter of the global population, characterized predominantly by steatosis and fibrosis, with potential progression to cirrhosis and hepatocellular carcinoma and potent immune cell involvement at all disease stages 4. Despite considerable research efforts, pharmacotherapeutic targets for MASLD remain largely elusive, with only a limited number of candidates demonstrating efficacy. Only recently, resmetirom (a selective thyroid hormone receptor beta agonist) was approved by the U.S. Food and Drug Administration (FDA), as the first treatment for patients with MASH with moderate to advanced liver fibrosis, along with diet management and exercise 5.

Lifestyle modifications, such as dietary restriction (DR) and exercise, remain the mainstay of recommended therapeutic interventions in overweight and obese individuals with MASLD, with the aim of reducing body weight by at least 10% 6, 7. However, the sustainability of long-term (hypocaloric) dietary management poses several challenges, particularly for individuals with severe obesity 8. Consequently, bariatric surgery (BS) has emerged as a viable option for obese patients, particularly when conventional nutritional and behavioral therapy fails to achieve the desired therapeutic outcomes. In a randomized controlled clinical trial, BS-based interventions were superior to lifestyle intervention alone in weight loss as well as in improving histological features of MASLD such as resolution of MASH and improvement of liver fibrosis by at least 1 stage 9. Furthermore, observational studies demonstrated long-term benefits regarding cardiovascular and liver-related morbidity as well as mortality 10.

Among various surgical modalities, sleeve gastrectomy (SG) is the most frequently performed procedure, supported by robust evidence of long-term efficacy and safety 11. SG restricts nutritional intake by reducing the stomach volume and capacity, without reconstructing the gastrointestinal tract 12. Emerging studies have suggested that BS potentially mitigates lipid accumulation and halts MASLD progression 13-15. Although the primary consequence of SG is the restriction of food intake (akin to DR), SG generates remarkable influences on systemic inflammation and gut microbiota composition 16. Recent advancements in omics technologies allow for sophisticated inquiries into liver diseases from diverse angles by profiling gene expression, protein function, metabolism, and microbiome activity 17. The latest study demonstrated that SG promotes liver regenerative capacities according to results from mouse models 18. Simultaneously, our previous study demonstrated that BS exerts superior protective effects in the livers of patients with MASLD compared to DR, a finding associated with increased macrophage infiltration in BS 14. However, little has been known about metabolic and inflammatory alterations of steatotic livers after SG and DR interventions.

Here, we sought to compare the molecular and cellular consequences of DR and SG related to the amelioration of liver steatosis, using a rodent obesity-MASLD model. We implemented single-cell (sc)/single-nuclei (sn) transcriptome analysis to describe intracellular alterations and extracellular interactions, and spatial metabolomics with the assistance of artificial intelligence (AI)-based image analysis assisted in deciphering metabolite production collating immune cell zonation. This study provides novel insights into multiomics alterations in steatotic livers from patients undergoing DR and SG, which may hopefully increase our understanding of mechanistic consequences of surgery-based therapies in MASLD.

Results

SG induces weight loss and ameliorates non-invasive biomarkers in human patients with MASLD

BMI and serological indexes were analyzed retrospectively in 18 patients who underwent SG (Fig. 1A). In comparison to baseline at one day pre-SG, significant declines of body mass index (BMI), glucose, postprandial blood glucose (PPG), aspartate aminotransferase (AST), alanine aminotransferase (ALT), gamma-glutamyltransferase (γ-GT), triglycerides (TG), high-density lipoprotein-cholesterol (HDL-C) and total bilirubin (TBIL) occurred at 6 months post-SG, while total bile acids (TBA) (p = 0.051) and fibroblast growth factor (FGF)-19 (p = 0.058) tended to increase (Fig. 1B-D, Table S2). This data hence revealed profound changes in varying processes, including inflammatory, metabolic and gut-derived endocrine signals. Intriguingly, male patients benefitted more from SG particularly on serological TG and low-density lipoprotein-cholesterol (LDL-C) levels than female patients (Table S2). Thus, to further investigate the cellular mechanisms involved in the liver response to SG and compare with DR, healthy wild-type male rats were used to recapitulate interventions on an obesity-related MASLD model (Fig. 2A, Fig. S1). All rats gained weight after 12 weeks of high-fat diet (HFD) feeding, whereas Sham + DR and SG rats had a reduced weight gain in the 8 weeks following the intervention compared to HFD alone. Moreover, DR and SG strikingly suppressed weight gain of rats from 13 to 17 weeks, while ameliorative effects compromised afterward in contrast to Sham (Fig. 2B). Serum analysis revealed that in comparison to Sham animals on HFD, Sham + DR suppressed the level of LDL but increased the level of HDL; SG significantly suppressed the levels of TG, LDL and AST, but increased the levels of HDL, TBA and FGF-19. Additionally, SG drastically suppressed cholesterol levels (p = 0.06) but led to increased serum FGF-21 (p = 0.06) (Fig. 2C). HFD exaggerated hepatic lipid accumulation in comparison to rats fed a normal diet.

Yet, both Sham + DR and SG intervention dramatically diminished hepatic lipid accumulation compared to the HFD-fed sham group (Fig. 2D and E). Simultaneously, fibrosis was increased after HFD, while it was significantly attenuated by both DR and SG (Fig. 2F and G). Results imply that both DR and SG ameliorated obesity progression, liver steatosis and fibrosis, while SG was the most efficient intervention for improving systemic levels of lipid metabolism and bile acid (BA) production.

Sc/sn-transcriptome analysis evidences an altered hepatic cellular landscape after DR and SG

Rat liver samples [n=3 each per Sham (on HFD), Sham + DR and SG group] were analyzed using sc/snRNA-sequencing analysis, with an explicit classification of main hepatic and immune cell populations [including two populations of hepatocytes: HepaC(a) and HepaC(b), cholangiocytes, hepatic stellate cells (HSC), endothelial cells, macrophages, monocytes, neutrophils, dendritic cells (DC, two populations), B lymphocytes, T lymphocytes and natural killer (NK) cells] (Fig. 3A, Fig. S2A). Reconstituted cell composition suggested higher numbers of NK cells, B cells, neutrophils and cholangiocytes leading to a lower proportion of hepatocytes and hepatic stellate cells upon Sham + DR (vs. Sham). In contrast, SG intervention led to an increased proportion of macrophages, while DR led to an increased proportion of NK cells but decreased proportions of hepatocytes and HSCs (Fig. 3A, Fig. S2B). The gene expression of metabolism-related markers (Cyp2e1, Cyp7a1, Ppara and Nr1h4) and key MASLD markers (Pnpla3 and Gcg) was analyzed in hepatocytes and immune cells in each group (Fig. 3B). Essentially, both Sham + DR and SG suppressed the gene expression of Pnpla3 (mainly in hepatocytes) and Gcg (in total level, hepatocytes and non-hepatocytes). Nonetheless, SG increased the gene expression of Cyp2e1, Ppara and Nr1h4 (in the non-hepatocytes) in contrast to Sham + DR (Fig. 3C). Immunohistochemistry (IHC) evidenced that patatin like phospholipase domain containing 3 (PNPLA3) protein was significantly downregulated in Sham + DR and SG groups (vs. Sham), while peroxisome proliferator activated receptor (PPAR)α-expressing liver cells increased in the SG group (vs. Sham) (Fig. 3D and E), thereby confirming the profound regulation towards beneficial metabolic pathways in hepatocytes on a protein level after SG.

Furthermore, differentially expressed genes (DEGs) were screened in parenchymal (Fig. 4A) and immune cells (Fig. 4B) by comparing Sham + DR and SG to Sham. Transcriptomic influences of diverse cell types in the liver upon Sham + DR and SG were assessed based on cell proportions and the amount of significantly up-/down-regulated DEGs. Hence, the distance of dots to the origin points (value '0') are considered as the magnitude of influences (Fig. 4C). As shown in the rainbow plot, hepatocytes stand distinctly out of other cell types and as the most influential ones by SG for both gene up- and down-regulation, implying that SG impacts steatotic livers mostly through altering hepatocytic functions (Fig. 4C). Contrastingly, Sham + DR appears to result in relatively uniform effects on liver cells. In accordance with the KEGG database, enriched signaling pathways were identified based on significantly upregulated DEGs. Up- and down-regulated DEGs at the bulk levels were compared between Sham, Sham + DR and SG groups, indicating enhanced metabolic processes (SG vs. Sham and SG vs. Sham + DR) and declined metabolic processes (Sham + DR vs. Sham) (Fig. S3), which were in line with findings on liver samples from patients from our earlier study 14.

Next, regulated signaling pathways were compared between the groups (Sham + DR vs. Sham, SG vs. Sham, and SG vs. Sham + DR) for each cell population (Fig. 5A). Fatty acid degradation-related gene expression was increased not only in hepatocytes (Sham + DR vs. Sham, SG vs. Sham, SG vs. Sham + DR) but also in cholangiocytes and endothelial cells (SG vs. Sham). Furthermore, PPAR signaling pathways, fatty acid degradation and cholesterol metabolism were more pronounced (not only in hepatocytes) by both Sham + DR and SG interventions (vs. Sham). Interestingly, in comparison to Sham + DR, the PPAR signaling pathway was promoted in cholangiocytes, endothelial cells and myeloid cell populations (monocytes, macrophages, DCs and neutrophils). Of note, butanoate metabolism was enhanced in hepatocytes upon both Sham + DR and SG (Fig. 5A). In addition, cellular interactions (total and ligand-receptor) were elucidated between each gradual pair of cell types among all three groups using the CellChat algorithm (Fig. S4A and B). Key signaling pathways that enriched according to total cellular interactions were described (Fig. S5A). In contrast to Sham, both Sham + DR and SG elevated ligand reception of neutrophils from immune cells, while SG elevated ligand release from cholangiocytes to other cells. Of note, in contrast to Sham + DR groups, SG promoted interaction levels of non-hepatocyte parenchymal cells (cholangiocytes, HSCs and endothelial cells) as well as immune cells (mainly for macrophages and neutrophils) (Fig. 5B). Further dissecting cholangiocyte-associated cellular interactions, significant signaling pathways from cholangiocytes to monocytes, macrophages, neutrophils, B lymphocytes, T lymphocytes and NK cells, as well as from monocytes, macrophages, DCs, hepatocytes, T lymphocytes and NK cells to cholangiocytes were illustrated (Fig. S5B and C).

Taken together, both Sham + DR and SG improve liver lipid metabolism (remarkably PPAR-associated processes) in steatotic livers. In particular, SG triggers metabolic enhancement broadly in liver cells, including hepatocytes, cholangiocytes, endothelial cells and myeloid cells. Noteworthy, cholangiocytes appear to actively interact with multiple liver cell types upon SG.

SG promotes liver metabolic capacities through myeloid cell accumulation

To explore the metabolic landscape of the liver, three representative liver samples were freshly acquired from the Sham, Sham + DR and SG groups and then analyzed using an ambient mass spectrometry imaging (MSI) metabolomics approach (Fig. 6A). The spectrum and distribution of total metabolites (with positive and negative ions) were illustrated (Fig. S6A - E). Histological features of liver samples together with spatially distributed metabolite abundance were displayed, indicating strongly reduced production of total metabolites in steatotic livers undergoing DR or SG, indicating the alleviation of global liver metabolism activities (Fig. 6B). Contrastingly, function enrichment analysis based on the total abundance of differentially upregulated metabolites elucidated that both Sham + DR and SG enhanced lipid metabolism in steatotic livers in contrast to Sham groups, potentially attenuating steatotic stress in liver. It is worth noting that SG significantly enhanced the capacities of linoleic acid metabolism, fat digestion and absorption, and cholesterol metabolism compared to Sham + DR (Fig. 6C).

AI-based pathological exploration has evolved as a promising and precise tool for MASLD investigations 19. Hence, in this study, a deep-learning model was established based on cell annotations (hepatocytes, lymphocytes and myeloid cells) that were supervised by liver pathology experts (Fig. 7A). Accordingly, the numbers of hepatocytes, lymphocytes and myeloid cells were assessed on a total of 18 rat liver samples. Importantly, hepatocyte proportions (70 - 80%) of total liver cells were not significantly altered by different interventions (Fig. 7B). Lymphocyte proportions were significantly increased in Sham with HFD (1.5 - 6.1%) vs. healthy control (0.05 - 2.1%) but declined upon Sham + DR (0.7 - 2.0%) and SG (0.5 - 2.2%) (vs. Sham) (Fig. 7B). Simultaneously, SG exacerbated the infiltration of myeloid cells (vs. healthy control, vs. Sham + DR) (Fig. 7B). In comparison, the numbers of lymphocytes and myeloid cells were assessed based on the sc/sn RNA-seq dataset. Myeloid cell proportions tended to increase by SG (Fig. S7A and B). Besides, the cell distribution of lymphocytes and myeloid cells was depicted to correspondingly map the metabolites in diverse cell regions (Fig. 7C). The abundance of total metabolites between regions with high or low myeloid cell infiltration was compared (high-/low-infiltration region selection displayed in Fig. 7C). SG remarkably increased metabolite abundance in regions with high myeloid cell infiltration (Fig. 7D). In addition, upregulated metabolites in regions rich in myeloid cells significantly refer to processes of autophagy, together with the metabolism of choline, linoleic acids, cholesterol and glycerophospholipids (Fig. 7E). Given the diverse phenotypic functions of macrophages in liver, a list of markers (Trem2, Spp1, Gpnmb, Mmp9, Mmp13, Plin2, Ppara, Pparg, Cd207, Cx3cr1, Clec4f, Nos2, Cd80, Nlrp3, Il1b, Tnfa, Cd163, Mrc1, Il10 and Tgfb1) was selected from the sc/sn RNA-seq dataset to assess phenotypes of scar-associated macrophage (SAM), lipid-associated macrophage (LAM), liver capsular macrophage (LCM), Kupffer cell (KC), pro- and anti-inflammatory macrophages 20-22. Comparing macrophage clusters (SG vs. Sham), markers of SAM, LCM and KC were upregulated (SG vs. Sham), whereas pro-inflammatory macrophage markers were downregulated (SG vs. Sham and Sham + DR vs. Sham) (Fig. 8A, Fig. S7C). Simultaneously, SG promoted classical phenotypes of monocytes but suppressed neutrophil maturation. In contrast, DR tended to promote conventional DC (cDC) 1 phenotype and neutrophil maturation (Fig. 8A, Fig. S7C). In addition, myeloid cells were further divided into subtypes (SAM, LAM, non-/classical monocytes, KC, cDC1/2, mature/naïve neutrophils and other macrophages) by phenotypic characteristics and annotated (Fig. 8B). In line, proportions of myeloid cell subtypes were generated for Sham, Sham + DR and SG groups. Intriguingly, SG led to a lower frequency of LAM and mature neutrophils but a higher frequency of SAM and KC in comparison to Sham and Sham + DR (Fig. 8C). Taken together, SG and DR enhanced lipid metabolism in the liver, while SG promoted myeloid cell accumulation in the liver, facilitating metabolic improvement and liver repair processes. Nonetheless, SG may fuel fibrosis via enhanced SAM population in the liver.

SG facilitates liver repair processes in steatotic livers

The strong accumulation of macrophages in SG-treated livers, in conjunction with the weight development and metabolic improvement, suggested the induction of myeloid cell-driven hepatic repair processes 23. To further assess the cellular phenotypic alterations in the liver, a list of markers was explored (for hepatocytes: Pcna, Ppara, Cyp2e1, Cyp7a1 and Afp; for cholangiocytes: Pcna, Sox9, Slc4a2, Hnf4a and Hnf1b; for HSCs: Pdgfrb, Tgfb1, Col1a1, Acta2 and Vim) (Fig. 9A). The expression of Pcna, Ppara, Cyp2e1 and Cyp7a1 in hepatocytes, the expression of Pcna and Hnf1b in cholangiocytes significantly upregulated by SG (vs. Sham). In contrast, the expression of Pcna in hepatocytes and cholangiocytes in Sham + DR, and the expression of Pdgfrb in HSCs in Sham + DR and SG were downregulated (vs. Sham) (Fig. 8A). Besides, cell differentiation was assessed and depicted using a cell trajectory analysis model (Fig. 8B). Cell trajectory of hepatocytes implied higher expression of Pcna, Cyp2e1 and Ppara in Sham + DR and SG groups than the Sham group, while the expression distinctly varied with the cell trajectory (Fig. 9B, Fig. S8). Cholangiocytes were classified into three groups. Cholangiocytes from Sham and SG groups were more associated with the expression of Pcna and Slc4a2 but not with the expression of Sox9 (Fig. 9B, Fig. S8). Furthermore, HSCs of Sham + DR groups were shown in relatively inactivated clusters associated with low Pdgfrb expression (Fig. 9B), but not associated with the expression of Tgfb1 and Col1a1 (Fig. S8). In addition, IHC analyses revealed the occurrence of ductular reaction by immunostaining of CK7, and interleukin (IL)-17A+ cells, typically associated with liver fibrosis 24, 25. The presence of CK7+ cells was significantly increased by steatosis induction, while it was alleviated by Sham + DR (vs. Sham) (Fig. 9C and D). Simultaneously, the presence of IL-17A+ cells was significantly increased by HFD but not in the Sham + DR and SG (vs. Sham) (Fig. 9C and D). Hence, our data indicated that Sham + DR and SG can attenuate the induction of fibrogenic immune response, while Sham + DR repressed ductular reaction. Taken together, SG and especially DR tend to reduce the involvement of portal area injury (particularly ductular reaction and fibrogenesis). Notably, DR appears to compromise regenerative capacities in hepatocytes, which requires further investigation.

DR affects zonation characteristics of periportal hepatocytes

Given the remarkable cell trajectory with variable expression of phenotypic markers, differential characteristic alterations in hepatocyte subtypes were hypothesized. Hence, hepatocytes were further dissected into subpopulations based on the sc/sn RNA-seq dataset (Fig. 10A). Metabolism-related function enrichment analysis revealed two distinct groups: metabolism one group (Met1, including cluster 6, 7 and 10) enriched metabolism processes of glycolysis, lipids and drugs, while metabolism two group (Met2, including clusters 1, 2, 3, 4, 5, 8 and 9) enriched fatty acid biosynthesis (Fig. 10B).

Interestingly, illustrations of the gene expression (pericentral hepatocyte markers: Ctnnb1, Ccne1, Lrp5, Lrp6, Myc and Glul; periportal hepatocyte markers: Cyp2f2 and Alb) further identified that the Met1 group harbors more periportal hepatocytes, while the Met2 group harbors more pericentral hepatocytes (Fig. 10C and D). Furthermore, cell trajectory analysis indicated the explicit differentiation status of these two metabolism groups, which supports the zonal classification (Fig. 10E). In comparison to periportal (Met1) hepatocytes, pericentral (Met2) hepatocytes illustrated higher numbers and strength of interactions with immune cells, while myeloid cells, T lymphocytes and NK cells exerted particularly high involvement in cellular interactions (Fig. S9A and B). Particularly, periportal phenotypes accounted for 15 - 20% of total hepatocytes in both Sham and SG groups but were almost absent (0.4%) in the Sham + DR group, indicating a phenotypic reshaping of hepatocytes induced by DR (Fig. 10F). In addition, metabolism markers [cytochrome (CYP)2E1 and CYP7A1] and cell proliferation (Ki67) were investigated in situ (Fig. 10G). Strikingly, proliferative cells (Ki67+) were more present in the liver tissues of the SG group than in other groups, indicating a strong potential for liver regeneration led by SG (Fig. 10G). CYP2E1 and CYP7A1 were determined to be more enriched surrounding central veins in Ctrl, Sham and SG groups, but to be homogeneously distributed in Sham + DR groups (Fig. 10G and H). In addition, the spatial distribution of fatty acid, BA and sphingolipid metabolites appeared to be enriched around periportal areas in Sham and SG groups but rather widely spread in Sham + DR groups (Fig. 10I). Taken together, DR deprives hepatocyte zonation features potentially by reshaping periportal hepatocyte characteristics, which may compromise the capacities of metabolism and regeneration in steatotic livers in comparison to SG.

Discussion

In this study, we employed multiomics methods to unravel relevant alterations of obesity-related parameters and liver steatosis following SG and DR in both patient cohorts and rat models. Specifically, sc/sn RNA-seq data unveiled an augmentation of lipid metabolism in both hepatocytes and immune cells, with a notable increase in myeloid cells following SG. This data was supported by spatially-resolved metabolomics analysis that revealed the strengthened metabolic capacities of liver cells induced by SG and DR. Notably, SG elicited elevated lipid metabolic capacities in regions exhibiting intense myeloid cell infiltration, underscoring the intricate interplay between metabolic and immune responses in the liver microenvironment. PPAR-α emerged as a prominent candidate mediating the observed metabolic enhancements induced by both SG and DR. Additionally, the increased serological levels of bile acids and gut-derived FGF-19, together with hepatic butanoate metabolism underscored the putative pivotal roles of the gut-liver axis in modulating liver metabolism post-SG.

To date, only a few studies have employed the bulk transcriptomic analysis of patients' liver samples post-BS to identify alterations of inflammation and lipid metabolism, yet those studies suggested potential therapeutic approaches for MASLD 26. Essentially, PPARα was elucidated as a key nuclear receptor responsible for weight-loss-induced MASLD amelioration 27, 28. Consequently, PPAR agonism has been regarded as a promising approach to tackling MASLD 29, 30. Robust preclinical data have determined the potential of targeting PPARα to improve lipid metabolism, inflammation, and fibrosis in MASLD 31, 32. In line, our previous study has proved that lanifibranor (a pan-PPAR agonist) significantly improves liver steatosis and fibrosis in mice and in an in-vitro liver-on-a-chip model of MASLD 33. However, numerous PPAR dual/pan agonists have failed to reach clinical practices, because they could not satisfy clinical endpoints (amelioration of e.g., liver injury, fibrosis, steatosis) 34. In addition, SG effectively enhanced the PPARα level and metabolic pathways broadly on hepatic parenchymal and immune cells, to a higher degree than DR. It has been demonstrated that PPARα activation in the liver can ameliorate the accumulation of circulating monocytes upon DR intervention 35. In addition, cell types were predicted from hematoxylin-eosin-stained rat liver sections using an AI-assisted digital method, regarded as an advanced tool for MASLD diagnosis 19, 36. Intriguingly, in rat MASLD models, SG can strengthen the macrophage pool despite a broad PPARα activation in the liver. Further characterization of the macrophage population indicates that SG endows the anti-inflammatory and repair-associated potential of macrophages rather than their detrimental pro-inflammatory roles in the liver.

It has been reported that lipid metabolism is critical in immune response during steatotic liver diseases 37-39. Accordingly, our study outlined plausible direct links between metabolism restoration and myeloid cell infiltration in post-SG steatotic rat livers. The diverse and contradictory roles of liver macrophages make them able to promote liver homeostasis, disease progression and remission depending on a plethora of possible phenotypic activation profiles 30, 40-42. Liver macrophages can mediate hepatic lipid metabolism during steatotic stress 43, 44, while lipid metabolism is closely associated with phenotypic alterations in macrophages 45, 46. Thus, the liver macrophage compartment remains an intriguing target for MASLD therapies 47-49. Moreover, the functions of liver macrophages vary distinctly by spatial distribution 50. In this study, we revealed that upon SG, the liver regions with high myeloid cell infiltration exert strong lipid metabolism capacity (e.g., autophagy and choline metabolism), while the causes are poorly understood. As reported, autophagy and choline metabolism play crucial roles in mediating myeloid cell (especially macrophage) phenotypes. For instance, lipid autophagy (term 'lipophagy') can potentiate their activation towards macrophage foam cells and a lipid-associated macrophage phenotype 51. Contrastingly, choline uptake and metabolism modulate macrophage inflammasome processes and promote a pro-inflammatory phenotype 52. Therefore, we hypothesize that SG-stimulated liver myeloid cells could better sense and assist metabolic stress in the liver, thereby leading to their heterogeneous distribution in the liver.

The hepatic metabolic and regenerative functions vary by hepatocyte zonation (periportal and pericentral), leading to diverse responses to injuries and therapies 53, 54. Emerging studies demonstrated that absent metabolic function in the periportal zone provokes cholestasis and ductular reaction 55, 56. On the other hand, compromised hepatic periportal characteristics may be attributed to dysregulated WNT signaling and cell regeneration 55, 57. Nowadays, the advances in spatially resolved multiomics approaches (e.g., sc/sn RNA-seq, spatial metabolomics, and spatial proteomics) provide both broader views and deeper understandings of heterogenic pathomechanisms of liver diseases 53, 58, 59. Our study elucidated that hepatic periportal characteristics were dramatically subsided by DR but preserved by SG. Since DR is rarely associated with adverse events in the liver, we hypothesized that periportal hepatocytes might acquire 'pericentral' features to aid hepatic metabolism. Previous studies have shown that fasting can influence hepatocyte zonation features, potentially by altering endoplasmic reticulum (ER)-mitochondria function in pericentral and mild-lobular hepatocytes 55, 60. However, a causative link between liver zonation and dietary intake or surgical interventions remains unknown. Likely, enhanced immunometabolism levels may interpret the preservation of zonal features upon the SG intervention. Furthermore, both DR and SG appeared to ameliorate stress in the liver portal area, illustrating less presence of ductular cells and fibrosis-promoters. In addition, our data indicated that the regenerative potentials of hepatocytes were enhanced by SG but not DR, which strengthened the findings from the latest study 18.

As observed in multiple clinical studies, FGF19 levels were increased in patients' circulation post-BS, thereby implying a potential regulatory mechanism of BS 61. Intriguingly, we also spotted that female patients tended to gain a higher FGF-19 increase than male patients post-SG. However, it has to be verified based on future studies with larger cohorts. Noteworthy, our study determined improved capacities of BA secretion and butanoate metabolism, as well as elevated FGF19 levels in serum upon the SG, suggesting that SG influences liver metabolism by regulating the gut-liver axis. The current understanding of MASLD pathogenesis is deficient in regard to insights into gut-liver crosstalk. This encompasses a number of key areas, including hepatic cell metabolism, inflammation, fibrogenesis, and complex cellular interactions, which play pivotal roles in MASLD progression 23, 30, 62. Increasing evidence indicates critical roles that BS may play in regulating the gut-liver axis, including GM, and BA metabolism 63. As suggested by recent studies, post-BS changes in the GM and BA circulation, as well as a decrease in the portal influx of free fatty acids, are beneficial to steatotic livers and obesity 64-66. Furthermore, BS can drastically repopulate GM and reverse the primary/secondary BA ratio in the intraluminal ileal, inducing improvement in metabolic syndrome in patients with MASLD 67. On the other hand, the butanoate production by GM was reported to be associated with bile acids and to mediate liver metabolism 68. Notably, butanoate/butyrate production can be enhanced after bariatric surgery 69, while favorable effects of butyrate have been proved on MASLD, including attenuating hepatic steatosis and inflammation while enhancing insulin sensitivity 70. Given the limited amount of literature focusing on the connection between PPAR-α, FGF-19, and butanoate metabolism, it is suggested that butyrate can upregulate PPAR-α in hepatocytes, thereby alleviating HFD-induced MASLD in rats through the activation of β-oxidation and inflammation suppression 71. Our study not only supports the fundamental roles of the gut-liver axis and gut-derived signals in protecting against liver steatosis, but also reveals disease ameliorations led by BS in contrast to simple DR.

Nonetheless, the limitations inherent in the current study should be acknowledged. The exclusion of patients with MASLD and obesity who solely underwent DR stemmed from challenges associated with dietary normalization to SG and inherent biases in dietary quantification. Furthermore, the sample size utilized for spatial metabolomics analysis was relatively modest, potentially constraining the generalizability of the findings. In addition, the in vitro modeling of BS remains challenging, while it is promising to employ multicellular systems (e.g., liver-on-a-chip 72) to break down mechanisms of immunometabolism in MASLD. Despite these limitations, this study highlights a prominent enhancement in liver immunometabolism upon SG. Notably, multi-omics approaches offer a comprehensive exploration of the cellular and metabolic landscape within steatotic livers, marking a pioneering endeavor in this domain. Beyond traditional analyses, the incorporation of sc/sn RNA transcriptome and spatially resolved metabolite distribution analyses provides unprecedented insights into the immunometabolism dynamics of MASLD post-DR and -SG.

Materials and Methods

Human cohort

The clinical retrospective study was conducted on 18 patients who received SG in The Third Affiliated Hospital of Nanjing Medical University/Changzhou Second People's Hospital (from January 2023 to July 2023). This cohort comprised n = 8 male (44.4%) and n = 10 female (55.6%) individuals at a median age of 34 years (range: 29 - 40). We included patients that fulfilled the following criteria: 1) BMI > 32.5 kg/m2; 2) Age: 16 - 65; 3) clinically diagnosed with MASLD (based on cardiometabolic risk factors and presence of steatosis on liver imaging; 4) diagnosed with type 2 diabetes or sleep apnea syndrome; 5) without a history of other hepatic or genetic illness; 6) received sleeve gastrectomy, without postoperative complications. Clinical indexes (e.g., BMI, serological indexes) were evaluated on the day before SG and 6 months post-surgery. This clinical retrospective study was approved by the ethics committee of The Third Affiliated Hospital of Nanjing Medical University/Changzhou Second People's Hospital (ID: [2023]KY124-01).

Animal models and surgical interventions

Male wild-type SPRAGUE DAWLEY® rats aged 8 weeks were obtained from Laboratory Animal Co., Ltd. (Changzhou, China). Rats were bred under specific-pathogen-free conditions and fed with HFD (D12492, Research Diets, Inc., USA) for 12 weeks. Five rats received sham surgery followed by uncontrolled HFD, 9 rats received sleeve gastrectomy followed by uncontrolled HFD, and 8 rats received sham surgery followed by restricted HFD feeding (comparable to dietary intake of post-SG rats) for 8 weeks. Procedures of sleeve gastrectomy were illustrated in supplementary data (Supplementary methods and Fig. S1) and were based on our previous research 73. Five healthy control rats were set with normal diet feeding but without any specific interventions. All these procedures were performed and supported by Kerbio Co., Ltd., Changzhou, China. Animal experiments were approved by the Institutional Animal Care and Use Committee (ID: IACUC23-0068).

Serological tests

Patients' blood samples were acquired from the median cubital vein. Rat blood samples were harvested through rat hearts immediately after rats were sacrificed. Serum was harvested using centrifugation (450×g, 10 min), and stored at -80℃. Circulating markers including AST, ALT, PPG, γ-GT, TG, cholesterol (CHO), HDL-C, LDL-C, alkaline phosphatase (ALP) and TBA were measured using Indiko analyzer (ThermoFisher, USA). Cytokines including IL-6, IL-10, FGF-19, FGF-21 and transferrin (TRF) were measured using ELISA kits and a microplate reader (Rayto, China). All measurements were conducted following the manufacturer's instructions. Antibodies that were used in this study are listed in Table S1.

Immunohistochemistry and histological exploration

Liver tissues were harvested immediately after rats were sacrificed. Tissues were fixed in 10%-formalin for 24 hours, followed by paraffin embedding. Sections of 5 µm thickness were obtained using a HistoCore BIOCUT microtome (Leica, Germany). Hematoxylin and eosin (HE) staining were performed using standard procedures. Extracellular collagen area was indicated using the Picro Sirius Red Stain Kit (Abcam, UK). The IHC method was employed using the kits purchased from Beyotime Biotechnology, China. Antigen retrieval was conducted using a heat-induced epitope retrieval method within a citrate buffer. Primary antibodies against the target proteins were applied, followed by appropriate secondary antibodies. Diaminobenzidine (DAB) was used for chromogenic detection, and counterstaining was performed with hematoxylin. Slides were examined under a microscope to evaluate protein localization and expression levels. In situ staining was observed under a light microscope (Olympus, Japan). Antibodies used in this study are listed in Table S1. For the quantitation of cell nucleus and protein levels, FIJI software was utilized 74. In particular, relative staining densities of proteins per µm2 were assessed and compared between peri-central and -portal regions in each tissue slide.

Single-cell and single-nuclei sequencing and analysis

Single cells and cell nuclei were isolated and extracted from freshly collected rat liver samples using the dissection kit purchased from OE Biotech Co., Ltd., Shanghai, China. The libraries were acquired following the manufacturer's protocol and then sequenced using the MobiNova platform (Shanghai, China). The integration of sc/sn-RNA sequencing data was operated following the common protocol with eliminating batch effects 75. Cell Ranger was used to filter the data, while the downstream analyses were conducted using Seurat (V5.0.1) 76, CellChat (1.6.0) 77, and Monocle2 78 software packages for R. Cell clusters were annotated based on the expression of representative cell makers and verified using the Enrichr database 79. Function enrichment was performed based on Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. All the technical procedures including sampling, quality control, library build-up and sequencing analysis were supported by the genome platform of OE Biotech Co., Ltd., Shanghai, China.

Spatial metabolomics analysis

Liver samples were freshly acquired from rats and then frozen at -80℃ immediately after optimal cutting temperature (OCT) embedding. Embedded tissue samples were sliced into 5 µm using a cryo-microtome CM 1950 (Leica, Germany). In-situ profiles of total metabolites were detected with positive and negative ion binding approaches using the MSI metabolomics method 80. Tissue regions were annotated after H&E staining. Metabolite enrichment was assessed based on the KEGG database. All the technical procedures including sample preparation, quality control, detection and analysis were supported by the metabolomics platform of Luming Biotechnology Co., Ltd. Shanghai, China.

AI-guided imaging analysis

The AI-guided imaging analysis performed herein is based on a previously published approach 81-83. H&E-stained rat liver sections were scanned using a light microscope (Olympus, Japan). Processing of the whole slide images (WSIs) was operated using QuPath 0.4.4 following a pre-defined protocol 84. Sparse annotations, comprising approximately 15 cells for each of the cell populations, were manually labeled across three WSIs representing different conditions (Sham, Sham + DR and SG). A simple artificial neural network classifier with two hidden layers, consisting of 20 and 10 neurons, respectively, was trained to distinguish between the three classes versus the rest. Example regions without these classes were used as negative examples during training. Subsequently, cell classification and measurements were exported, and for each cell population, 2D hexagonal binning plots of the number of cells were generated in Python (using the matplotlib function 'hexbin') 85. The model was trained for representative cells (at least 15 cells/sample' of hepatocytes, lymphocyte-like cells and myeloid-like cells) based on three liver tissue slides, which were defined with the assistance of pathological expertise. The trained model was applied to the sample cohort of rat liver sections (including 5 Ctrl, 5 Sham, 8 Sham + DR and 8 SG samples). Detailed methodology is described in the Supplementary methods.

Statistical analysis

The GraphPad Prism 9.0 software (GraphPad Software, USA) and the ggplot2 package (version: 3.4.4) for R were used to generate plots 86. The statistical tests used in this study are indicated in the figure legend for each panel. Data are presented as the mean ± S.D. 'p < 0.05' was considered a significant difference.

Supplementary Material

Supplementary methods, figures and tables.

Xiurong Cai and Hanyang Liu were supported by China Scholarship Council Foundation. The authors thank Shanghai OE Biotech CO., Ltd. and Shanghai Luming Biological Technology Co., Ltd. for the access to the MobiNova single-cell system and AFADESI-MSI platforms. Schematic figures were generated via BioRender.

Funding

National Nature Science Foundation of China (82300730); Changzhou Health Commission of China (CZQM2022007); Changzhou Foreign Expert Project from Changzhou Science and Technology Bureau of China (CQ20214053); Changzhou Medical Center of Nanjing Medical University (CZKYCMCB202221); German Research Foundation (DFG Ta434/8-1, SFB/TRR 296 and CRC1382, Project-ID 403224013).

Author contributions

H.L. conceptualized the research. H.L. and A.G. supervised the research. S.C., Q.Z. and H.L. developed methodology. S.C. and P.S. performed animal experiments. Q.Z. and C.K. performed computational analysis. S.C. and J.X. performed experiments and clinical data collection. H.L., X.C., G.Y. and S.C. performed data analysis. The manuscript was written by H.L. and revised by A.G. and F.T. H.L., H.L., A.G. and F.T. acquired funding.

Ethics committee approval and patient consent

This clinical retrospective study was approved by the ethics committee of The Third Affiliated Hospital of Nanjing Medical University/Changzhou Second People's Hospital (ID: [2023]KY124-01). Animal experiments in this study were approved by the Institutional Animal Care and Use Committee (ID: IACUC23-0068).

Data and materials availability

Raw and processed data of sc/sn-RNA sequencing analysis can be obtained from the GEO database (https://www.ncbi.nlm.nih.gov/geo/) with the reference number GSE254841. Raw data of spatial metabolomics can be obtained from the OMIX, China National Center for Bioinformation (https://ngdc.cncb.ac.cn/omix/) with the reference number OMIX005793. Technical materials for the deep-learning-based image analysis model may be provided by the corresponding authors upon reasonable requests. All data are available in the main text or the supplementary materials.

Figure 1 Sleeve gastrectomy induces weight loss and MASLD amelioration in patients. (A) Description of the patient cohort receiving SG. Comparisons of (B) BMI and (C) PPG, together with (D) ALT, AST, TG, γ-GT, ALP, CHO, HDL-C, LDL-C, TBA, TBIL, TRF, FGF-19, FGF-21, IL-6 and IL-10 in patients between one day pre-SG and 6 months post-SG. Abbreviations: SG: sleeve gastrectomy; BMI: body mass index; PPG: postprandial blood glucose; AST: aspartate aminotransferase; ALT: alanine aminotransferase; γ-GT: gamma-glutamyltransferase; TG, triglycerides; CHO, cholesterol; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol; TBIL, total bilirubin; ALP, alkaline phosphatase; TBA, total bile acid. TRF: transferrin; IL: interleukin; FGF, fibroblast growth factor. The paired t-test was performed. '*' represents 'p < 0.05' and statistical significance.

Figure 2 Sleeve gastrectomy and dietary restriction induce weight loss and MASLD amelioration in rat models. (A) The technical scheme of the rat MASLD model and surgical interventions. All animals except “Ctrl” were on HFD during this experiment. (B) Body weight was recorded for 20 weeks of modeling. (C) Serological indexes (glucose, triglycerides, cholesterol, HDL-C, LDL-C, TBA, ALT, AST, ALP, IL-6, IL-10, FGF-19 and FGF-21) were measured. (D) H&E-stained rat liver tissue and (E) quantification of steatosis (area of hepatocyte lipid vesicles) from Ctrl, Sham, Sham + DR and SG groups. (F) Sirius red-stained rat liver tissue and (G) quantification of fibrosis from Ctrl, Sham, Sham + DR and SG groups. Abbreviations: Sham/S: sham surgery; DR: dietary restriction; SG: sleeve gastrectomy; Ctrl: control fed a normal diet; HDL-C: high-density lipoprotein-cholesterol; LDL-C: low-density lipoprotein-cholesterol; TG: triglycerides; TBA: total bile acids; ALT: alanine transaminase; AST: aspartate aminotransferase; ALP: alkaline phosphatase; IL: interleukin; FGF: fibroblast growth factor. The one-way ANOVA test followed by Tukey's multiple comparisons was performed. '*' represents 'p < 0.05' and statistical significance.

Figure 3 Sleeve gastrectomy improves hepatic metabolism according to sc/sn-transcriptome analysis and in-situ exploration. (A) UMAPs describe the distribution of B cells, cholangiocytes, DCs, endothelial cells, hepatic stellate cells, hepatocytes, macrophages, monocytes, neutrophils, T cells and NK cells. Histograms describing cell proportions. (B) Gene expressions of Cyp2e1, Cyp7a1, Ppara, Nr1h4, Pnpla3 and Gcg were illustrated on rat liver cells upon Sham, Sham + DR and SG. (C) Gene expression changes of Cyp2e1, Cyp7a1, Ppara, Nr1h4, Pnpla3 and Gcg were illustrated on rat liver cells upon Sham, Sham + DR and SG. Protein levels of PNPLA3, PPARα and FXR were illustrated (D) and quantified (E) on rat liver tissue of Ctrl, Sham, Sham + DR and SG groups. Abbreviations: Ctrl: control fed a normal diet; Sham: sham surgery; DR: dietary restriction; SG: sleeve gastrectomy; PPAR: peroxisome proliferator-activated receptor; PNPLA: patatin-like phospholipase domain-containing protein; FXR: farnesoid X receptor; CK: cytokeratin; IL: interleukin; DC: dendritic cells; NK: natural killer cells; B: B lymphocytes; T: T lymphocytes; Neu: neutrophils; Mac: macrophages; Mono: monocytes; Endo: endothelial cells; HSC: hepatic stellate cells; Chol: cholangiocytes; HepaC: hepatocytes. The one-way ANOVA test followed by Tukey's multiple comparisons was performed. '*' represents 'p < 0.05' and statistical significance.

Figure 4 Sc/sn-transcriptome analysis illustrates differentially expressed genes in diverse liver cells. (A & B) Differentially expressed genes of hepatocytes, cholangiocytes, endothelial cells, hepatic stellate cells, B, T, NK, DCs, macrophages, monocytes and neutrophils (Sham + DR vs. Sham and SG vs. Sham) were illustrated in volcano plots. (C) Transcriptomic influences of hepatocytes, cholangiocytes, endothelial cells, hepatic stellate cells, B, T, NK, DCs, macrophages, monocytes and neutrophils in the liver were assessed with DEG numbers and cell proportion of total cells Abbreviations: Sham: sham surgery; DR: dietary restriction; SG: sleeve gastrectomy; DC: dendritic cells; NK: natural killers; B: B lymphocytes; T: T lymphocytes; Neu: neutrophils; Mac: macrophages; Mono: monocytes; Endo: endothelial cells; HSC: hepatic stellate cells; Chol: cholangiocytes; HepaC: hepatocytes; DEG: differentially expressed gene. The one-way ANOVA test followed by Tukey's multiple comparisons was performed. 'Adjust p < 0.05' represents statistical significance.

Figure 5 Sc/sn-transcriptome analysis illustrates up-regulated enriched pathways and cellular interactions in diverse liver cells. Significantly up-regulated (A) KEGG signaling pathways and (B) cellular interactions (numbers and strength) in hepatocytes, cholangiocytes, endothelial cells, hepatic stellate cells, B, T, NK, DCs, macrophages, monocytes and neutrophils were illustrated (Sham + DR vs. Sham and SG vs. Sham). Abbreviations: Sham: sham surgery; DR: dietary restriction; SG: sleeve gastrectomy; HepaC: hepatocytes; Chol: cholangiocytes; HSC: hepatic stellate cells; Endo: endothelial cells; Mono: monocytes; Mac: macrophages; DC: dendritic cells; Neu: neutrophils; NK: natural killer cells; B: B lymphocytes; T: T lymphocytes; KEGG: Kyoto encyclopedia of genes and genomes. The Fisher's exact test and the one-way ANOVA test followed by Tukey's multiple comparisons were performed. 'Adjust p<0.05' represents statistical significance.

Figure 6 Sleeve gastrectomy and dietary restriction enhance capacities of lipid metabolism in steatotic livers. (A) Technical scheme of the ambient mass spectrometry imaging metabolomics method. (B) H&E-stained rat liver sample section was mapped with the abundance of total metabolites in Sham, Sham + DR and SG groups. (C) Function enrichment of upregulated differentially generated metabolites (Sham + DR vs. Sham, SG vs. Sham and SG vs. Sham + DR) was analyzed based on the KEGG database. Abbreviations: Sham: sham surgery; DR: dietary restriction; SG: sleeve gastrectomy; KEGG: Kyoto encyclopedia of genes and genomes. The Fisher's exact test was performed. 'Adjusted p < 0.05' represents statistical significance.

Figure 7 Sleeve gastrectomy facilitates lipid metabolism associated with myeloid cell infiltration in the liver. (A) The technical scheme of the deep-learning-based image analysis. The deep-learning approach assessed (B) proportions of hepatocytes, lymphocytes and myeloid cells in total liver cells per tissue section and (C) spatial distribution in representative liver samples. Regions with higher (red circles) and lower (green circles) myeloid cell densities were selected (D) Differentially generated metabolites were illustrated (higher myeloid cell infiltrated regions vs. lower myeloid cell infiltrated regions) in Sham, Sham + DR and SG groups. (E) Function enrichment of upregulated differentially generated metabolites (higher myeloid cell infiltrated regions vs. lower myeloid cell infiltrated regions) in the liver sample of SG group was analyzed based on the KEGG database. Abbreviations: C: control; Sham/S: sham surgery; DR: dietary restriction; SG: sleeve gastrectomy; KEGG: Kyoto encyclopedia of genes and genomes. The unpaired t-test, the Fisher's exact test and the one-way ANOVA test followed by Tukey's multiple comparison were performed. '*' represents 'p < 0.05' and statistical significance.

Figure 8 Sleeve gastrectomy induces phenotypic alterations in myeloid cells. (A) The expression of phenotype-associated markers in myeloid cell clusters was illustrated based on sc/sn RNA-seq data (SG vs. Sham). (B) Myeloid cells were characterized into phenotypic subtypes (SAM, LAM, non-/classical monocytes, KC, cDC1/2, mature/naïve neutrophils and other macrophages) with annotations. (C) The proportion of myeloid cell subtypes in Sham, Sham + DR and SG groups were illustrated. Abbreviations: Sham: sham surgery; DR: dietary restriction; SG: sleeve gastrectomy; sc/sn: single-cell/single-nuclei; Mono: monocytes; Mac: macrophages; (c)DC: (conventional) dendritic cells; Neu: neutrophils; SAM: scar-associated macrophage; LAM: lipid-associated macrophage; LCM: liver capsular macrophage; KC: Kupffer cell. The unpaired t-test and Fisher's exact test were performed. '*' represents 'p < 0.05' and statistical significance.

Figure 9 Dietary restriction and sleeve gastrectomy alter the cellular differentiation of hepatocytes, cholangiocytes and hepatic stellate cells. (A) Cellular characteristic markers of HepaC (Pcna, Ppara, Cyp2e1, Cyp7a1 and Afp), Chol (Pcna, Sox9, Slc4a2, Hnf4a and Hnf1b) and HSC (Pdgfrb, Tgfb1, Col1a1, Acta2 and Vim) were displayed. (B) Cell trajectories of HepaC, Chol and HSC together with Cyp2e1 expression in HepaC, Pcna expression in Chol, and Pdgfrb expression in HSC. CK7- and IL-17A- expressing cells were (C) illustrated in situ and (D) quantified from total rat liver cells. Abbreviations: C: control; Sham/S: sham surgery; DR: dietary restriction; SG: sleeve gastrectomy; HepaC: hepatocyte; Chol: cholangiocyte; HSC: hepatic stellate cell; KEGG: Kyoto encyclopedia of genes and genomes. The one-way ANOVA test was performed. '*' represents 'p < 0.05' and statistical significance.

Figure 10 Dietary restriction compromises zonation characteristics of periportal hepatocytes. (A) Subclusters of rat hepatocytes from Sham, Sham + DR and SG are illustrated in UMAPs. (B) Metabolism enrichment on each hepatocyte subcluster is displayed in a matrix. (C) Gene expression of Ctnnb1, Cnne1, Lrp5, Lrp6, Myc, Glul, Cyp2f2 and Cyp7a1 illustrated in rat hepatocytes. (D) Metabolism group-1 and -2 in rat hepatocytes. (E) Metabolism group-1 and -2 in the cell trajectory of rat hepatocytes. (F) Proportions of metabolism group-1 and -2 in rat hepatocytes of Sham, Sham + DR and SG groups. (G) Protein levels of Ki67, CYP2E1 and CYP7E1 illustrated in rat liver tissue (each view contains one or more portal and central areas) of Sham, Sham + DR and SG groups. (H) Relative differences of CYP2E1 and CYP7E1 levels between peri-central and -portal regions in each rat liver tissue of Sham, Sham + DR and SG groups. (I) Spatial distribution of fatty acid, bile acid and sphingolipid metabolites. In H&E- and IHC-characterized rat liver sample sections, periportal regions were outlined in red circles and pericentral regions were outlined in blue circles. Abbreviations: Sham: sham surgery; DR: dietary restriction; SG: sleeve gastrectomy; HepaC: hepatocyte; Met: metabolism; CYP: cytochrome. H&E: hematoxylin and eosin. '*' represents 'p < 0.05' and statistical significance.
==== Refs
1 Collaboration NCDRF Worldwide trends in underweight and obesity from 1990 to 2022: a pooled analysis of 3663 population-representative studies with 222 million children, adolescents, and adults Lancet 2024 403 1027 50 38432237
2 Pillon NJ Loos RJF Marshall SM Zierath JR Metabolic consequences of obesity and type 2 diabetes: Balancing genes and environment for personalized care Cell 2021 184 1530 44 33675692
3 Quek J Chan KE Wong ZY Tan C Tan B Lim WH Global prevalence of non-alcoholic fatty liver disease and non-alcoholic steatohepatitis in the overweight and obese population: a systematic review and meta-analysis Lancet Gastroenterol Hepatol 2023 8 20 30 36400097
4 Devarbhavi H Asrani SK Arab JP Nartey YA Pose E Kamath PS Global burden of liver disease: 2023 update J Hepatol 2023 79 516 37 36990226
5 Harrison SA Bedossa P Guy CD Schattenberg JM Loomba R Taub R A Phase 3, Randomized, Controlled Trial of Resmetirom in NASH with Liver Fibrosis N Engl J Med 2024 390 497 509 38324483
6 Younossi ZM Zelber-Sagi S Henry L Gerber LH Lifestyle interventions in nonalcoholic fatty liver disease Nat Rev Gastroenterol Hepatol 2023 20 708 22 37402873
7 European Association for the Study of the Liver, et al EASL-EASD-EASO Clinical Practice Guidelines on the management of metabolic dysfunction-associated steatotic liver disease (MASLD) J Hepatol. 2024: S0168-8278(24)00329-5
8 Greenway FL Physiological adaptations to weight loss and factors favouring weight regain Int J Obes (Lond) 2015 39 1188 96 25896063
9 Verrastro O Panunzi S Castagneto-Gissey L De Gaetano A Lembo E Capristo E Bariatric-metabolic surgery versus lifestyle intervention plus best medical care in non-alcoholic steatohepatitis (BRAVES): a multicentre, open-label, randomised trial Lancet 2023 401 1786 97 37088093
10 Aminian A Al-Kurd A Wilson R Bena J Fayazzadeh H Singh T Association of Bariatric Surgery With Major Adverse Liver and Cardiovascular Outcomes in Patients With Biopsy-Proven Nonalcoholic Steatohepatitis JAMA 2021 326 2031 42 34762106
11 Zhang L Tan WH Chang R Eagon JC Perioperative risk and complications of revisional bariatric surgery compared to primary Roux-en-Y gastric bypass Surg Endosc 2015 29 1316 20 25294534
12 Brajcich BC Hungness ES Sleeve Gastrectomy JAMA 2020 324 908 32870299
13 Lopez Lopez AP Tuli S Lauze M Becetti I Pedreira CC Huber FA Changes in Hepatic Fat Content by CT 1 Year After Sleeve Gastrectomy in Adolescents and Young Adults With Obesity J Clin Endocrinol Metab 2023 108 e1489 e95 37403207
14 Chen S Tang L Guillot A Liu H Bariatric Surgery Associates with Nonalcoholic Steatohepatitis/Hepatocellular Carcinoma Amelioration via SPP1 Suppression Metabolites 2022 13 11 36676937
15 Chen S Jiang Y Qi X Song P Tang L Liu H Bioinformatics analysis to obtain critical genes regulated in subcutaneous adipose tissue after bariatric surgery Adipocyte 2022 11 550 61 36036283
16 Cerreto M Santopaolo F Gasbarrini A Pompili M Ponziani FR Bariatric Surgery and Liver Disease: General Considerations and Role of the Gut-Liver Axis Nutrients 2021 13 2649 34444807
17 Thiele M Villesen IF Niu L Johansen S Sulek K Nishijima S Opportunities and barriers for omics-based biomarker discovery in steatotic liver diseases Journal of Hepatology 2024 81 345 359 38552880
18 Yu T Ma X Cheng Y Wang Z Zhang G Ding H Amelioration of NAFLD by sleeve gastrectomy-triggered hepatocyte regeneration in mice - experimental research Int J Surg 2024 110 3307 3325 38573134
19 Ratziu V Hompesch M Petitjean M Serdjebi C Iyer JS Parwani AV Artificial intelligence-assisted digital pathology for non-alcoholic steatohepatitis: current status and future directions J Hepatol 2024 80 335 51 37879461
20 Guilliams M Bonnardel J Haest B Vanderborght B Wagner C Remmerie A Spatial proteogenomics reveals distinct and evolutionarily conserved hepatic macrophage niches Cell 2022 185 379 96 e38 35021063
21 Peiseler M Araujo David B Zindel J Surewaard BGJ Lee WY Heymann F Kupffer cell-like syncytia replenish resident macrophage function in the fibrotic liver Science 2023 381 eabq5202 37676943
22 Ramachandran P Dobie R Wilson-Kanamori JR Dora EF Henderson BEP Luu NT Resolving the fibrotic niche of human liver cirrhosis at single-cell level Nature 2019 575 512 8 31597160
23 Peiseler M Schwabe R Hampe J Kubes P Heikenwalder M Tacke F Immune mechanisms linking metabolic injury to inflammation and fibrosis in fatty liver disease - novel insights into cellular communication circuits J Hepatol 2022 77 1136 60 35750137
24 Li N Yamamoto G Fuji H Kisseleva T Interleukin-17 in Liver Disease Pathogenesis Semin Liver Dis 2021 41 507 15 34130335
25 Guillot A Hamdaoui N Bizy A Zoltani K Souktani R Zafrani ES Cannabinoid receptor 2 counteracts interleukin-17-induced immune and fibrogenic responses in mouse liver Hepatology 2014 59 296 306 23813495
26 Haas JT Vonghia L Mogilenko DA Verrijken A Molendi-Coste O Fleury S Transcriptional Network Analysis Implicates Altered Hepatic Immune Function in NASH development and resolution Nat Metab 2019 1 604 14 31701087
27 Lefebvre P Lalloyer F Bauge E Pawlak M Gheeraert C Dehondt H Interspecies NASH disease activity whole-genome profiling identifies a fibrogenic role of PPARalpha-regulated dermatopontin JCI Insight 2017 2 e92264 28679947
28 Puengel T Liu H Guillot A Heymann F Tacke F Peiseler M Nuclear Receptors Linking Metabolism, Inflammation, and Fibrosis in Nonalcoholic Fatty Liver Disease Int J Mol Sci 2022 23 2668 35269812
29 Kamata S Honda A Ishikawa R Akahane M Fujita A Kaneko C Functional and Structural Insights into the Human PPARalpha/delta/gamma Targeting Preferences of Anti-NASH Investigational Drugs, Lanifibranor, Seladelpar, and Elafibranor Antioxidants (Basel) 2023 12 1523 37627519
30 Tacke F Puengel T Loomba R Friedman SL An integrated view of anti-inflammatory and antifibrotic targets for the treatment of NASH J Hepatol 2023 79 552 66 37061196
31 Hu P Li K Peng X Kan Y Li H Zhu Y Nuclear Receptor PPARalpha as a Therapeutic Target in Diseases Associated with Lipid Metabolism Disorders Nutrients 2023 15
32 Gong L Wei F Gonzalez FJ Li G Hepatic fibrosis: Targeting peroxisome proliferator-activated receptor alpha from mechanism to medicines Hepatology 2023 78 1625 53 36626642
33 Lefere S Puengel T Hundertmark J Penners C Frank AK Guillot A Differential effects of selective- and pan-PPAR agonists on experimental steatohepatitis and hepatic macrophages☆ J Hepatol 2020 73 757 70 32360434
34 Kamata S Honda A Ishii I Current Clinical Trial Status and Future Prospects of PPAR-Targeted Drugs for Treating Nonalcoholic Fatty Liver Disease Biomolecules 2023 13 1264 37627329
35 Jordan S Tung N Casanova-Acebes M Chang C Cantoni C Zhang D Dietary Intake Regulates the Circulating Inflammatory Monocyte Pool Cell 2019 178 1102 14 e17 31442403
36 Uche-Anya E Anyane-Yeboa A Berzin TM Ghassemi M May FP Artificial intelligence in gastroenterology and hepatology: how to advance clinical practice while ensuring health equity Gut 2022 71 1909 15 35688612
37 Abimannan T Parthibane V Le SH Vijaykrishna N Fox SD Karim B Sphingolipid biosynthesis is essential for metabolic rewiring during T(H)17 cell differentiation Sci Adv 2024 10 eadk1045 38657065
38 Xu JC Chen ZY Huang XJ Wu J Huang H Niu LF Multi-omics analysis reveals that linoleic acid metabolism is associated with variations of trained immunity induced by distinct BCG strains Sci Adv 2024 10 eadk8093 38578989
39 Gao B Tsukamoto H Inflammation in Alcoholic and Nonalcoholic Fatty Liver Disease: Friend or Foe? Gastroenterology 2016 150 1704 9 26826669
40 Vonderlin J Chavakis T Sieweke M Tacke F The Multifaceted Roles of Macrophages in NAFLD Pathogenesis Cell Mol Gastroenterol Hepatol 2023 15 1311 24 36907380
41 Kohlhepp MS Liu H Tacke F Guillot A The contradictory roles of macrophages in non-alcoholic fatty liver disease and primary liver cancer-Challenges and opportunities Front Mol Biosci 2023 10 1129831 36845555
42 Guillot A Tacke F Liver macrophages revisited: The expanding universe of versatile responses in a spatiotemporal context Hepatol Commun 2024 8 e0491 38967563
43 Lefere S Tacke F Macrophages in obesity and non-alcoholic fatty liver disease: Crosstalk with metabolism JHEP Rep 2019 1 30 43 32149275
44 Ji PX Chen YX Ni XX Miao Q Hua J Effect of triggering receptor expressed on myeloid cells 2-associated alterations on lipid metabolism in macrophages in the development of non-alcoholic fatty liver disease J Gastroenterol Hepatol 2024 39 369 80 38012119
45 Yuan Y Zhang Y Lu X Li J Wang M Zhang W Novel insights into macrophage immunometabolism in nonalcoholic steatohepatitis Int Immunopharmacol 2024 131 111833 38503012
46 Sangineto M Ciarnelli M Colangelo T Moola A Bukke VN Duda L Monocyte bioenergetics: An immunometabolic perspective in metabolic dysfunction-associated steatohepatitis Cell Rep Med 2024 5 101564 38733988
47 Wen Y Lambrecht J Ju C Tacke F Hepatic macrophages in liver homeostasis and diseases-diversity, plasticity and therapeutic opportunities Cell Mol Immunol 2021 18 45 56 33041338
48 Tacke F Targeting hepatic macrophages to treat liver diseases J Hepatol 2017 66 1300 12 28267621
49 Ju C Tacke F Hepatic macrophages in homeostasis and liver diseases: from pathogenesis to novel therapeutic strategies Cell Mol Immunol 2016 13 316 27 26908374
50 Adrien G Frank T Spatial dimension of macrophage heterogeneity in liver diseases eGastroenterology 2023 1 e000003
51 Robichaud S Fairman G Vijithakumar V Mak E Cook DP Pelletier AR Identification of novel lipid droplet factors that regulate lipophagy and cholesterol efflux in macrophage foam cells Autophagy 2021 17 3671 89 33590792
52 Sanchez-Lopez E Zhong Z Stubelius A Sweeney SR Booshehri LM Antonucci L Choline Uptake and Metabolism Modulate Macrophage IL-1beta and IL-18 Production Cell Metab 2019 29 1350 62 e7 30982734
53 Paris J Henderson NC Liver zonation, revisited Hepatology 2022 76 1219 30 35175659
54 Porat-Shliom N Compartmentalization, cooperation, and communication: The 3Cs of Hepatocyte zonation Curr Opin Cell Biol 2024 86 102292 38064779
55 Martini T Naef F Tchorz JS Spatiotemporal Metabolic Liver Zonation and Consequences on Pathophysiology Annu Rev Pathol 2023 18 439 66 36693201
56 Kang SWS Cunningham RP Miller CB Brown LA Cultraro CM Harned A A spatial map of hepatic mitochondria uncovers functional heterogeneity shaped by nutrient-sensing signaling Nat Commun 2024 15 1799 38418824
57 Hoshida Y Nijman SM Kobayashi M Chan JA Brunet JP Chiang DY Integrative transcriptome analysis reveals common molecular subclasses of human hepatocellular carcinoma Cancer Res 2009 69 7385 92 19723656
58 Santos AA Delgado TC Marques V Ramirez-Moncayo C Alonso C Vidal-Puig A Spatial metabolomics and its application in the liver Hepatology 2024 79 1158 79 36811413
59 Guillot A Winkler M Silva Afonso M Aggarwal A Lopez D Berger H Mapping the hepatic immune landscape identifies monocytic macrophages as key drivers of steatohepatitis and cholangiopathy progression Hepatology 2023 78 150 66 36630995
60 Parlakgul G Pang S Artico LL Min N Cagampan E Villa R Spatial mapping of hepatic ER and mitochondria architecture reveals zonated remodeling in fasting and obesity Nat Commun 2024 15 3982 38729945
61 Ryan PM Hayward NE Sless RT Garwood P Rahmani J Effect of bariatric surgery on circulating FGF-19: A systematic review and meta-analysis Obes Rev 2020 21 e13038 32329176
62 Hammerich L Tacke F Hepatic inflammatory responses in liver fibrosis Nat Rev Gastroenterol Hepatol 2023 20 633 46 37400694
63 Aron-Wisnewsky J Dore J Clement K The importance of the gut microbiota after bariatric surgery Nat Rev Gastroenterol Hepatol 2012 9 590 8 22926153
64 Musso G Gambino R Cassader M Interactions between gut microbiota and host metabolism predisposing to obesity and diabetes Annu Rev Med 2011 62 361 80 21226616
65 le Roux CW Aylwin SJ Batterham RL Borg CM Coyle F Prasad V Gut hormone profiles following bariatric surgery favor an anorectic state, facilitate weight loss, and improve metabolic parameters Ann Surg 2006 243 108 14 16371744
66 Koulas SG Stefanou CK Stefanou SK Tepelenis K Zikos N Tepetes K Gut Microbiota in Patients with Morbid Obesity Before and After Bariatric Surgery: a Ten-Year Review Study (2009-2019) Obes Surg 2021 31 317 26 33130944
67 Talavera-Urquijo E Beisani M Balibrea JM Alverdy JC Is bariatric surgery resolving NAFLD via microbiota-mediated bile acid ratio reversal?. A comprehensive review Surg Obes Relat Dis 2020 16 1361 9 32336663
68 Zhou F Mai T Wang Z Zeng Z Shi J Zhang F The improvement of intestinal dysbiosis and hepatic metabolic dysfunction in dextran sulfate sodium-induced colitis mice: effects of curcumin J Gastroenterol Hepatol 2023 38 1333 45 37210613
69 Park YS Ahn K Yun K Jeong J Baek KW Lee J Alterations in gastric and gut microbiota following sleeve gastrectomy in high-fat diet-induced obese rats Sci Rep 2023 13 21294 38042896
70 Coppola S Avagliano C Calignano A Berni Canani R The Protective Role of Butyrate against Obesity and Obesity-Related Diseases Molecules 2021 26 682 33525625
71 Sun B Jia Y Hong J Sun Q Gao S Hu Y Sodium Butyrate Ameliorates High-Fat-Diet-Induced Non-alcoholic Fatty Liver Disease through Peroxisome Proliferator-Activated Receptor alpha-Mediated Activation of beta Oxidation and Suppression of Inflammation J Agric Food Chem 2018 66 7633 42 29961332
72 Liu H Yin G Kohlhepp MS Schumacher F Hundertmark J Hassan MIA Dissecting Acute Drug-Induced Hepatotoxicity and Therapeutic Responses of Steatotic Liver Disease Using Primary Mouse Liver and Blood Cells in a Liver-On-A-Chip Model Adv Sci (Weinh) 2024 11 e2403516 38868948
73 Qi X Li X Jiao Y Chen S Song P Qian Z The Effects of Antral Preservation and Antral Resection on Body Composition, Glycemic Control and Bone Mineral Density Following Vertical Sleeve Gastrectomy in C57BL/6J Mice with Obesity and Type 2 Diabetes Diabetes Metab Syndr Obes 2022 15 297 305 35140488
74 Albersheim M Fedje-Johnston W Carlson C Arnoczky SP Toth F Shea K Cell Count and Cell Density Decrease as Age Increases in Cadaveric Pediatric Medial Menisci Arthrosc Sports Med Rehabil 2023 5 100795 37868658
75 Slyper M Porter CBM Ashenberg O Waldman J Drokhlyansky E Wakiro I A single-cell and single-nucleus RNA-Seq toolbox for fresh and frozen human tumors Nat Med 2020 26 792 802 32405060
76 Slovin S Carissimo A Panariello F Grimaldi A Bouche V Gambardella G Single-Cell RNA Sequencing Analysis: A Step-by-Step Overview Methods Mol Biol 2021 2284 343 65 33835452
77 Jin S Guerrero-Juarez CF Zhang L Chang I Ramos R Kuan CH Inference and analysis of cell-cell communication using CellChat Nat Commun 2021 12 1088 33597522
78 Qiu X Mao Q Tang Y Wang L Chawla R Pliner HA Reversed graph embedding resolves complex single-cell trajectories Nat Methods 2017 14 979 82 28825705
79 Kuleshov MV Jones MR Rouillard AD Fernandez NF Duan Q Wang Z Enrichr: a comprehensive gene set enrichment analysis web server 2016 update Nucleic Acids Res 2016 44 W90 7 27141961
80 Unsihuay D Mesa Sanchez D Laskin J Quantitative Mass Spectrometry Imaging of Biological Systems Annu Rev Phys Chem 2021 72 307 29 33441032
81 Calderaro J Ghaffari Laleh N Zeng Q Maille P Favre L Pujals A Deep learning-based phenotyping reclassifies combined hepatocellular-cholangiocarcinoma Nat Commun 2023 14 8290 38092727
82 Zeng Q Klein C Caruso S Maille P Allende DS Minguez B Artificial intelligence-based pathology as a biomarker of sensitivity to atezolizumab-bevacizumab in patients with hepatocellular carcinoma: a multicentre retrospective study Lancet Oncol 2023 24 1411 22 37951222
83 Zeng Q Klein C Caruso S Maille P Laleh NG Sommacale D Artificial intelligence predicts immune and inflammatory gene signatures directly from hepatocellular carcinoma histology J Hepatol 2022 77 116 27 35143898
84 Muti HS Loeffler CML Echle A Heij LR Bülow RD Krause J The Aachen Protocol for Deep Learning Histopathology: A hands-on guide for data preprocessing 2020
85 Han S Kwak IY Mastering data visualization with Python: practical tips for researchers J Minim Invasive Surg 2023 26 167 75 38098348
86 Ito K Murphy D Application of ggplot2 to Pharmacometric Graphics CPT Pharmacometrics Syst Pharmacol 2013 2 e79 24132163
