1. Medicine

Mid-zone hepatocytes trade proliferation for survival via Atf4-Chop axis in early acute liver injury

  1. Yaying Zhu
  2. Chengxiang Deng
  3. Bo Chen
  4. Jia He
  5. Yanan Liu
  6. Shan Lei
  7. Weiju Lu
  8. Cheng Peng  Is a corresponding author
  9. Zhao Shan  Is a corresponding author
  1. Yunnan Key Laboratory of Cell Metabolism and Diseases, Center for Life Sciences, School of Life Sciences, Yunnan University, China
https://doi.org/10.7554/eLife.108646.4
Share this article
Cite this article
  1. Yaying Zhu
  2. Chengxiang Deng
  3. Bo Chen
  4. Jia He
  5. Yanan Liu
  6. Shan Lei
  7. Weiju Lu
  8. Cheng Peng
  9. Zhao Shan
(2026)
Mid-zone hepatocytes trade proliferation for survival via Atf4-Chop axis in early acute liver injury
eLife 14:RP108646.
https://doi.org/10.7554/eLife.108646.4
  2. Download BibTeX
  3. Download .RIS

eLife Assessment

The manuscript is important in that it has the potential to advance our understanding of how liver regeneration is orchestrated after injury of pericentral hepatocytes, which are a target of many drugs and toxins. The findings are solid and suggest that midlobular hepatocytes are protected from toxin-induced injury by a genetic program delaying proliferation to prevent cell death. The effects of individual actors in the proposed program are convincingly shown but whether they act in concert and whether cell cycle arrest is indeed the protective mechanism remains to be determined.

https://doi.org/10.7554/eLife.108646.4.sa0

Abstract

Hepatocytes undergo extensive proliferation to facilitate liver repair after injury, yet early adaptive changes prior to proliferation remain unclear. Using a mouse model of acetaminophen (APAP)-induced liver injury, we show that hepatocytes undergo transient proliferation suppression that is most pronounced in mid-zone hepatocytes, consistent with zonal APAP metabolism. Although spatial transcriptomics (ST) alone could not independently confirm zonal proliferative restraint, ST analysis identified a unique mid-zone stress-response program centered on the Atf4-Chop axis. By integrating Cut&Run with immunohistochemistry and functional studies, we found that the Atf4-Chop axis suppresses proliferation via the cell-cycle inhibitor Btg2. Together, our findings support a model in which mid-zone hepatocytes transiently prioritize stress adaptation over proliferation, thereby preserving regenerative capacity for subsequent liver repair.

Introduction

Acetaminophen (APAP) overdose is the leading cause of drug-induced acute liver failure worldwide, resulting in over 10,000 hospitalizations and approximately 500 deaths annually in the United States (Bernal and Wendon, 2013; Larson et al., 2005; Ostapowicz et al., 2002; Lee, 2007). The hepatotoxicity of APAP is initiated by its metabolic activation through cytochrome P450 enzymes, primarily in pericentral hepatocytes, generating the reactive metabolite N-acetyl-p-benzoquinone imine (NAPQI) (Dahlin et al., 1984; Hinson et al., 2004; Jack et al., 2010; Tujios and Fontana, 2011). NAPQI depletes glutathione and forms protein adducts, leading to oxidative stress and mitochondrial dysfunction that culminate in centrilobular necrosis (Dahlin et al., 1984; Hinson et al., 2004; Jack et al., 2010; Tujios and Fontana, 2011). Although the mechanisms of APAP bioactivation and toxicity are well characterized, the early adaptive responses of hepatocytes—particularly their proliferative dynamics across lobular zones—remain poorly understood.

The liver is organized into functional hexagonal units termed lobules, where hepatocytes display spatial heterogeneity across three distinct zones: periportal (PP, zone 1), mid-zonal (Mid, zone 2), and pericentral (PC, zone 3) (Ben-Moshe et al., 2022; Wang et al., 2024; Wu et al., 2024). These zones exhibit unique gene expression profiles and metabolic specializations (Halpern et al., 2017). PC hepatocytes (zone 3), expressing high levels of cytochrome P450 enzymes, are the primary site of APAP-induced necrosis (Wang et al., 2024). However, recent spatial profiling, functional, and lineage tracing studies have identified Mid hepatocytes (zone 2) as the predominant source of new hepatocytes during both homeostasis and repair, demonstrating significant proliferative capacity and metabolic plasticity (Wei et al., 2021; He et al., 2021; Chembazhi et al., 2021). Despite these advances, how mid-zone hepatocytes dynamically adapt during the initial phase of acute stress remains unknown.

A key regulator of cellular stress adaptation is the integrated stress response (ISR), which phosphorylates eukaryotic initiation factor 2α (eIF2α) to suppress global translation while selectively upregulating stress-responsive genes, including activating transcription factor 4 (Atf4) and its downstream target pro-apoptotic factor C/EBP homologous protein (Chop, encoded by DNA damage-inducible transcript 3 [Ddit3]) (Costa-Mattioli and Walter, 2020; Vasudevan et al., 2020). Although Chop is traditionally associated with apoptosis, its transcriptional targets exhibit significant overlap with those of Atf4, encompassing genes that paradoxically support cell survival and growth (Marciniak et al., 2004; Han et al., 2013; Uzi et al., 2013). These findings underscore the context-dependent duality of the Atf4-Chop axis across different stress conditions, highlighting the need to elucidate its precise regulatory roles in specific physiological settings.

Through spatial transcriptomics (ST) and functional analyses, we demonstrate that zonal heterogeneity in APAP metabolism drives a stress-adaptive response in mid-zone hepatocytes during early APAP-induced liver injury. This response is mediated by the Atf4-Chop-Btg2 axis, which orchestrates a transient proliferation arrest to prioritize cellular survival over regenerative capacity. Our findings elucidate the spatial and molecular mechanisms governing hepatocyte adaptation to acute injury, uncovering a critical survival-regeneration trade-off in the early phase of acute liver injury.

Results

Mid-zone hepatocytes exhibit the most pronounced proliferative decline during early acute liver injury

To assess hepatocyte proliferation during early APAP-induced liver injury, wild-type (WT) C57BL/6J mice received an intraperitoneal injection of 300 mg/kg APAP, and liver sections were analyzed at 0, 3, 6, 12, and 24 hr post-injection. Immunohistochemical staining for the proliferation marker Ki67 revealed that under homeostatic conditions (0 hr), mid-zone hepatocytes exhibited approximately threefold higher basal proliferative activity than PP and PC hepatocytes, consistent with previous reports (Wei et al., 2021). During the early injury phase (3–12 hr), the number of Ki67+ mid-zone hepatocytes declined sharply from roughly 70–20 per section, reaching levels comparable to those observed in the PP and PC zones (Figure 1A and B). By 24 hr post-APAP, mid-zone proliferation gradually recovered toward baseline levels (Figure 1A and B). Although all zones exhibited reduced proliferation during the injury-initiation phase, the mid-zone showed the most substantial decline (Figure 1C), as further corroborated by value-change analysis relative to the 0 hr time point (Figure 1D). Collectively, these data identify the mid-zone as the region most profoundly affected by transient proliferative suppression during early APAP-induced acute liver injury.

Figure 1
Download asset Open asset
Mid-zone hepatocytes exhibit the most pronounced proliferative decline in early acute liver injury.

Wild-type C57BL/6J mice were intraperitoneally injected with 300 mg/kg acetaminophen (APAP) for 0, 3, 6, 12, and 24 hr. (A) Immunohistochemical staining of Ki67 to detect proliferating hepatocytes in liver sections at 0, 3, 6, 12, and 24 hr post-APAP. Representative images of pericentral (PC), middle (Mid), and periportal (PP) are shown. Red arrows indicate Ki67+ hepatocytes. (B) Quantitative analysis of Ki67+ hepatocytes across liver zones (PC, Mid, PP) at 0, 3, 6, 12, and 24 hr post-APAP. (C) Quantitative analysis of Ki67+ hepatocytes across 0, 3, 6, 12, and 24 hr post-APAP in each liver zone (PC, Mid, PP). (D) Changes of Ki67+ hepatocytes across liver zones (PC, Mid, PP) at 0, 3, 6, 12, and 24 hr post-APAP. Value change = number of Ki67+ hepatocytes at x hr post-APAP—the number of Ki67+ hepatocytes at 0 hr post-APAP. (a) denotes significance between PC and Mid regions, (b) denotes significance between PC and PP regions, and (c) denotes significance between Mid and PP regions. Data points represent mean hepatocyte counts from six high-power fields per zonal layer (n=3 mice/group) in B and C. Data are represented as means ± SD; one-way ANOVA (B, C, D). *p<0.05; **p<0.01; ***p<0.001; ****p<0.0001; ns, not significant.

Figure 1—source data 1

Numerical data of Figure 1B–D.

https://cdn.elifesciences.org/articles/108646/elife-108646-fig1-data1-v1.xlsx

Mid-zone hepatocytes show distinct gene expression profiles

To investigate the mechanism underlying proliferation pause during early APAP-induced liver injury, liver sections from WT mice treated with APAP for various time points were subjected to ST using the Visium (10× Genomics) platform (Figure 2—figure supplement 1A). Given inter-mouse variability, one representative sample per time point was selected based on H&E histology and serum ALT levels (Figure 2—figure supplement 1B). Under normal conditions (APAP 0 hr), unbiased transcriptional clustering partitioned each lobule into three distinct zones: PP (Alb, Mup20, Cyp2f2), PC (Glul, Cyp2e1, Cyp1a2), and mid-zone (intermediate between PP and PC, with elevated Igfbp2 and Hamp) (Figure 2A, Figure 2—figure supplement 1C). Following APAP-induced injury (3 hr, 6 hr), the PC zone remained identifiable by residual PC signature genes despite necrosis and reduced transcription, while PP genes remained largely unchanged. The mid-zone was distinguishable as the intervening cluster exhibiting marked transcriptional reprogramming (e.g. Sqstm1, Igfbp1) (Figure 2A, Figure 2—figure supplement 1C). Uniform Manifold Approximation and Projection (UMAP) and spatial heatmaps of Glul and Cyp2f2 confirmed persistent, though dynamic, zonal clustering across all time points (Figure 2—figure supplement 1D–F). These data indicate that the fundamental lobular architecture remains spatially resolved during the early injury phase, enabling zone-specific transcriptional analysis.

Figure 2 with 1 supplement see all
Download asset Open asset
Mid-zone hepatocytes show distinct gene expression profiles.

(A) Heatmap of the top 10 differentially expressed genes in PP, Mid, and PC hepatocyte zones at 0, 3, 6, and 24 hr post-acetaminophen (APAP). To highlight relative expression shifts across zones, this heatmap is visualized using row-wise Z-score standardization, which scales each gene to a uniform range and emphasizes spatial switching patterns rather than absolute transcript abundance. Functional categories at the bottom were assigned using gene ontology (GO) terms (biological process, cellular component, molecular function); for genes annotated with multiple GO terms, the most representative or significantly enriched term was selected. (B) Immunofluorescence staining for glutamine synthetase (GS, a PC marker) with DAPI (blue) in liver sections at 0, 3, 6, and 24 hr after APAP administration. Quantification of the zonal distribution of GS-positive cells at the indicated time points. (C) Immunohistochemical detection of p62 (Sqstm1 protein, a midzone marker during early injury) in liver sections at 0, 3, 6, and 24 hr after APAP administration. Red arrows indicate representative p62-positive hepatocytes. Quantification of the zonal distribution of p62-positive cells at the indicated time points. (D) Immunohistochemical staining of Cyp2f2 (a PP marker) in liver sections at 0, 3, 6, and 24 hr after APAP administration. Quantification of the zonal distribution of Cyp2f2-positive cells at the indicated time points. Data are presented as mean ± SD. One-way ANOVA (B–D). *p<0.05, **p<0.01, ***p<0.001, ****p<0.0001; ns, not significant.

Figure 2—source data 1

Numerical data of Figure 2B–D.

https://cdn.elifesciences.org/articles/108646/elife-108646-fig2-data1-v1.xlsx

Differential expression analysis across zones and time points revealed distinct zonal trajectories (Figure 2A). In the PP zone, gene expression remained relatively stable across 0–6 hr, with inflammatory genes (Orm1, Serpina1c) rising only at 24 hr. The mid-zone exhibited a clear temporal progression. At 0 hr, it expressed proliferative (Igfbp2), metabolic (Scd1, Cdo1), and homeostatic (Hamp, Cps1) genes. By 3–6 hr, these were replaced by stress-response and autophagy genes (Hmox1, Hspa8, Atf3, Ubc, Sqstm1). By 24 hr, wound-healing genes (Fgg, Fgl1) became prominent. The PC zone showed early downregulation of xenobiotic metabolism (Cyp2e1, Cyp1a2), followed by upregulation of cytoskeletal (Actb, Krt8) and redox (Prdx1, Srxn1) genes at 24 hr. This temporal trajectory identifies the mid-zone as the most transcriptionally dynamic region, undergoing a functional switch from homeostasis to active stress adaptation during the initial injury phase.

Focusing on the mid-zone at 6 hr versus 0 hr, volcano plot analysis identified significant upregulation of stress-related genes (Atf3, Hmox1, Ddit3, Atf4) alongside cell-cycle and apoptosis regulators (Btg2, Egr1, Sqstm1) (Figure 2—figure supplement 1G). Among these, we selected Sqstm1 (encoding p62) for validation as a representative stress-response gene with previously characterized roles in autophagy and oxidative stress. To validate our zonation assignment, we compared ST clusters with a nine-layered radial classification from the central vein (CV) to the portal vein (PV). Immunostaining for GS (PC), p62 (mid-zone), and Cyp2f2 (PP) confirmed that our PC, Mid, and PP zones corresponded to layers 1–2, 3–6, and 7–9, respectively (Figure 2B–D). These spatial correlations confirm that the mid-zone is the primary site of early stress-response gene induction, with p62 serving as a robust regional marker.

ST-based analysis of Ki67-associated spots showed a temporal pattern broadly consistent with the transient reduction in proliferation observed by immunohistochemistry (Figure 2—figure supplement 1I). To further dissect cell-cycle dynamics, we assessed S-phase and G2/M-phase gene signatures. Mid-zone hepatocytes displayed a transient increase in S-phase scores at 3 hr, followed by a decline at 6 and 24 hr; conversely, G2/M scores were reduced at 3–6 hr and returned to baseline by 24 hr (Figure 2—figure supplement 1J). These ST findings provide spatially resolved, exploratory evidence of changes in cell-cycle-associated transcription during early APAP-induced injury, consistent with the proliferative suppression observed by Ki67 immunohistochemistry.

Zonal metabolism of APAP leads to stress response in mid-zone hepatocytes

To determine what triggers the mid-zone stress response, we first examined the spatial distribution of APAP-metabolizing cytochrome P450 (Cyp) enzymes. Cyp enzymes metabolize APAP to NAPQI, which causes oxidative stress and cellular damage unless neutralized by GSH (Figure 3A). ST analysis revealed a zonal gradient in the expression of most Cyp transcripts, with highest levels in the PC zone progressively decreasing toward the PP zone (Figure 3B, Figure 3—figure supplement 1A). Immunohistochemical staining for Cyp1a2 confirmed this gradient at the protein level under homeostatic conditions (0 hr), with strong PC signal diminishing toward the portal tracts. However, following APAP treatment (3–24 hr), Cyp1a2 expression was largely lost within the necrotic CV region, while the adjacent mid-zone remained the primary area retaining detectable Cyp1a2 signal (Figure 3C and D). This spatial shift suggests that selective necrosis of Cyp-rich PC hepatocytes effectively relocates the residual metabolic capacity to the surviving mid-zone parenchyma (Figure 3E).

Figure 3 with 1 supplement see all
Download asset Open asset
Zonal metabolism of acetaminophen (APAP) leads to stress response in mid-zone hepatocytes.

(A) Schematic figure illustrating the metabolic process of liver injury caused by APAP. (B) Dot plot displays the log-normalized absolute transcript abundance of the Cyp family in hepatocyte zones (zones periportal [PP], mid-zonal [Mid], and pericentral [PC]) at 0, 3, 6, 24 hr post-APAP. (C) Immunohistochemical staining of Cyp1a2 to detect its dynamic changes in hepatocyte zones (from central vein [CV] to portal vein [PV]) at 0, 3, 6, 24 hr post-APAP. Injured area is outlined by black dashed lines. (D) Quantification of zonal distribution of Cyp1a2-positive cells in liver sections at 0, 3, 6, and 24 hr post-APAP, respectively. (E) Schematic figure illustrating expression changes of Cyp1a2 before and after APAP. (F) Immunofluorescence staining of Chop protein (green) to observe the dynamic changes in its expression across hepatocyte zones (from CV to PV) at 6 hr post-APAP with doses of 100, 300, or 500 mg/kg. Cyp1a2 protein (red) staining highlights the area around the CV. Cell nuclei are stained with DAPI (blue). (G) Quantification of zonal distribution of Chop-positive cells in liver sections from mice treated with varying doses of APAP—100, 300, or 500 mg/kg. (H) Schematic explaining why mid-lobular hepatocytes display a distinct gene expression program early after APAP injury. Cyp enzymes are expressed in a decreasing gradient from the CV to the PV. APAP is therefore metabolized primarily by pericentral Cyp-expressing hepatocytes; their necrosis shifts the main site of residual Cyp activity to the mid-lobular zone. High local APAP metabolism then drives a specific stress-response program, accounting for the unique early transcriptional signature of the mid-zone. Data are represented as means ± SD; one-way ANOVA (D, G). *p<0.05; **p<0.01; ***p<0.001; ****p<0.0001; ns, not significant.

Figure 3—source data 1

Numerical data of Figure 3D and E.

https://cdn.elifesciences.org/articles/108646/elife-108646-fig3-data1-v1.xlsx

To test whether the extent of necrosis dictates the magnitude and localization of the mid-zone stress response, we administered escalating APAP doses (100, 300, and 500 mg/kg). H&E staining confirmed a dose-dependent expansion of necrotic areas (Figure 3—figure supplement 1B and C). We then assessed Chop (the protein product of Ddit3), a key mid-zone stress marker whose transcript was upregulated at 6 hr post-APAP (Figure 2B). Chop immunoreactivity was absent at the lowest dose, where necrosis was minimal. As the dose increased and the necrotic area expanded, Chop expression emerged immediately adjacent to the necrotic front and progressively shifted toward the PV, always surrounding the injured territory (Figure 3F and G). Thus, CV necrosis establishes the adjacent peri-necrotic region as the functional epicenter of residual APAP metabolism, directly linking zonal Cyp distribution to localized stress response (Figure 3H).

The Atf4-Chop axis emerges as pivotal in mid-zone during early acute injury

To identify upstream drivers of the Mid stress program, we constructed gene regulatory networks (GRNs) from the ST data and ranked the top 10 regulons per zone at 3 and 6 hr post-APAP (Figure 4A, Figure 4—figure supplement 1A). GRN analysis revealed simultaneous activation of several transcription factors (TFs) in the mid-zone across both time points, including Atf4, Fos, Nfil3, Egr1, Maff, Atf3, and Ddit3. These TFs include known stress-response regulators (Atf4, Atf3, Ddit3), immediate-early genes (Fos, Egr1), and oxidative/metabolic regulators (Maff, Nfil3), suggesting a multi-faceted adaptive program (Woo et al., 2009; Ohoka et al., 2005; Bunch et al., 2023; von Scheidt et al., 2021; Lin et al., 2024). Among the top 10 mid-zone regulons, Ddit3 exhibited the highest inferred activity, while Atf4 ranked seventh (Figure 4B, Figure 4—figure supplement 1B). Since the Atf4-Ddit3 axis is a well-established regulator of stress adaptation (Zhang et al., 2022), the coordinated activation of this transcriptional network prompted us to investigate this pathway further. Violin plots confirmed prominent mid-zone enrichment of both Atf4 and Ddit3 at 3 and 6 hr (Figure 4C, Figure 4—figure supplement 1C). GRNs illustrated connectivity between Atf4 and its targets (e.g. Gdf15, Hspa8, Ddit3) as well as between Ddit3 and its targets (e.g. Hspa8, Hmox1, Osgin1, Gadd45a, Egr1) (Figure 4—figure supplement 1D). Heatmaps further demonstrated that the target genes of both TFs were selectively expressed in mid-zone hepatocytes at 6 hr post-APAP (Figure 4D), and many of these targets overlapped with the mid-zone differentially expressed genes (DEGs) identified earlier. This coordinated enrichment and target overlap strongly implicate the Atf4-Ddit3 axis as a principal driver of the mid-zone transcriptional reprogramming.

Figure 4 with 3 supplements see all
Download asset Open asset
The Atf4-Ddit3 axis emerges as pivotal in mid-zone during early acute injury.

(A) Heatmap displays the area under the curve (AUC) scores of transcription factor (TF) motifs, estimated per dot by single-cell regulatory network inference and clustering (SCENIC), highlighting gene regulatory networks and differentially activated TF motifs in hepatocyte zones (periportal [PP], mid-zonal [Mid], pericentral [PC]) at 6 hr post-acetaminophen (APAP). Columns represent TF motifs, rows represent dots, and color intensity indicates AUC scores. (B) Heatmap shows inferred TFs' activity across different zonal regions at 6 hr post-APAP. Activity was quantified as regulon enrichment scores using AUCell. The color scale represents regulon activity scores. (C) Violin plots show regulon activity of Atf4 and Ddit3 in each zonation at 6 hr post-APAP. Activity was quantified as regulon enrichment scores using AUCell. (D) Gene ontology (GO) pathway analysis reveals the top 15 enriched pathways for genes co-regulated by Atf4 and Ddit3 in Mid hepatocyte zones at 6 hr post-APAP. (E) Heatmap showing the target genes’ expression of Atf4 (left) and Ddit3 (right) in hepatocyte zones (zones PP, Mid, and PC) at 6 hr post-APAP. (F) Immunohistochemistry staining of p-eIF2A in liver sections from mice at 0 hr and 3 hr post-APAP. p-eIF2A-positive hepatocytes are indicated by red arrows. The number of p-eIF2A-positive cells per field of view (FOV) is quantified. n=3 mice/group. (G) Heatmap displays stress response and ER stress-related genes from differentially expressed genes (DEGs) identified in Mid hepatocyte zones at 0 and 6 hr post-APAP. Data are represented as means ± SD; unpaired Student’s t-test (F). *p<0.05; **p<0.01; ***p<0.001; ****p<0.0001; ns, not significant.

Figure 4—source data 1

Numerical data of Figure 4F.

https://cdn.elifesciences.org/articles/108646/elife-108646-fig4-data1-v1.xlsx

GO analysis of Atf4-Ddit3 co-regulated targets revealed enrichment for pathways related to unfolded protein response, cellular stress, autophagy, and apoptotic signaling (Figure 4E, Figure 4—figure supplement 1E and F). Immunohistochemical staining confirmed increased phosphorylation of eIF2α (p-eIF2α) in the mid-zone at 3 hr post-APAP, indicating ISR activation (Figure 4F). However, when we examined canonical ER-stress markers, we found that genes such as Eif2ak3, Atf6, and Xbp1 either remained unchanged or were downregulated at 3 and 6 hr post-APAP, despite strong Atf4 and Ddit3 upregulation (Figure 4G, Figure 4—figure supplement 1G). These data indicate that the mid-zone stress response engages the ISR through eIF2α phosphorylation but does not constitute a classical, full-scale ER stress/unfolded protein response. Instead, the response appears selectively routed through the Atf4-Ddit3 branch.

Reanalysis of a public APAP-injury scRNA-seq dataset (Matchett et al., 2024) confirmed zonal hepatocyte clustering (PP, Mid, PC) at 3 and 6 hr post-APAP, with mid-zone enrichment of Hmox1, Atf3, and Ubc (Figure 4—figure supplement 2A–C). GO analysis of mid-zone upregulated genes highlighted stress-response and cell-death processes (Figure 4—figure supplement 2D and E). Single-cell regulatory network inference and clustering (SCENIC) analysis independently identified activation of Ddit3 and Atf3 motifs in the mid-zone (Figure 4—figure supplement 2F and G). This external dataset validates our ST-based observations and reinforces the specificity of the mid-zone stress-regulatory program.

Analysis of human single-nucleus RNA-seq (snRNA-seq) data from healthy and APAP-induced acute liver failure (ALF) explants (GSE223561) revealed distinct PP (ALB, CYP2F2), Mid (IGFBP2, HAMP), and PC (GLUL, CYP2E1) hepatocyte populations (Figure 4—figure supplement 3A–D). In healthy human livers, Mid hepatocytes showed detectable expression of stress regulators (DDIT3, ATF4, HMOX1) alongside low cell-cycle gene expression. In APAP patients, Mid stress regulators were further elevated; however, in contrast to our early murine model, proliferation-associated genes (MKI67, CCNB1, CDK1) were also upregulated in the same Mid population (Figure 4—figure supplement 3E). Human ST data from two APAP patients revealed inter-patient heterogeneity: one specimen showed stress-gene enrichment with reduced proliferation markers in the peri-necrotic mid-zone, while the other displayed a predominantly proliferative transcriptional profile with weaker stress signaling (Figure 4—figure supplement 3F). These human data suggest that while the Mid stress program is conserved, the balance between stress adaptation and proliferation is dynamic and likely stage-dependent, consistent with the temporal progression observed in our mouse model.

The Atf4-Ddit3 axis protects hepatocytes from liver injury

To assess the functional relevance of the Atf4–Ddit3 axis, we first examined endogenous Atf4 protein dynamics. Immunohistochemical staining revealed progressive nuclear Atf4 accumulation in mid-zone hepatocytes between 0 and 6 hr post-APAP (Figure 5A). This temporal and spatial pattern aligns with the transcriptional upregulation observed in our ST data, confirming mid-zone-specific Atf4 activation.

Figure 5
Download asset Open asset
The Atf4-Ddit3 axis protects hepatocytes from liver injury.

(A) Immunohistochemical staining of Atf4 in liver sections at 0 and 6 hr post-acetaminophen (APAP). Red arrows indicate Atf4-positive hepatocytes. Zonal distribution of Atf4-positive cells in liver sections at 6 hr post-APAP is quantified. The statistic is the percentage of Atf4-positive hepatocytes in each layer over the total number of Atf4-positive hepatocytes. n=3 mice. (B) Schematic figure illustrating the overexpression of Atf4 via AAV in hepatocytes of wild-type C57BL/6J mice. Overexpression of EGFP is used as a control. (C) H&E staining showing liver morphology from mice overexpressing AAV-TBG-EGFP or AAV-TBG-Atf4 at 6 hr post-APAP. The injured area is outlined by black dashed lines. The percentage of injured area is quantified. n=3 mice/group. Necrotic areas outlined by loss of cellular architecture on H&E. (D) Serum levels of ALT are measured in mice overexpressing AAV-TBG-EGFP or AAV-TBG-Atf4 at 6, 24, 48, and 72 hr post-APAP. n=6 mice/group. (E) Immunohistochemistry staining of Atf4 in liver sections from mice overexpressing AAV-TBG-EGFP or AAV-TBG-Atf4 at 24 hr post-APAP. Atf4-positive hepatocytes are indicated by red arrows. The number of Atf4-positive cells per field of view (FOV) is quantified. n=9 mice/group. (F) Immunofluorescence staining of Ki67 (red) was performed to assess proliferating hepatocytes in mice overexpressing AAV-TBG-EGFP or AAV-TBG-Atf4 at 72 hr post-APAP. Cell nuclei were stained with DAPI (blue). The number of Ki67-positive cells per FOV is quantified. n=6 mice/group. (G) Immunohistochemistry staining of p-eIF2A in liver sections from mice overexpressing AAV-TBG-EGFP or AAV-TBG-Atf4 at 24 hr post-APAP. p-eIF2A-positive hepatocytes are indicated by red arrows. The number of p-eIF2A-positive cells per FOV is quantified. n=6 mice/group. Data are represented as means ± SD; unpaired Student’s t-test (C–G). One-way ANOVA (A). *p<0.05; **p<0.01; ***p<0.001; ****p<0.0001; ns, not significant.

Figure 5—source data 1

Numerical data of Figure 5A and C–G.

https://cdn.elifesciences.org/articles/108646/elife-108646-fig5-data1-v1.xlsx

To test causality, we overexpressed Atf4 specifically in hepatocytes using adeno-associated virus 8 (AAV8)—thyroxine-binding globulin (TBG) (Figure 5B). Atf4 overexpression markedly attenuated APAP-induced liver injury, as evidenced by reduced necrotic areas on H&E staining and significantly lower serum ALT levels (Figure 5C and D). Immunohistochemical staining confirmed nuclear Atf4 localization in hepatocytes adjacent to injured regions (Figure 5E). Furthermore, Atf4 overexpression resulted in decreased Ki67+ proliferating cells and increased p-eIF2α immunoreactivity (Figure 5F and G). Together, these data demonstrate that elevated Atf4 actively promotes hepatocyte survival by enhancing ISR signaling, suppressing proliferation, and mitigating early liver damage.

The Atf4-Ddit3 axis protects the liver by pausing proliferation via Btg2

Atf4 and Ddit3 are key transcriptional regulators of ISR, and their downstream target genes converge on overlapping stress-responsive pathways (Ohoka et al., 2005; Woo et al., 2009; Kaspar et al., 2021). To dissect the downstream effectors of the Atf4-Ddit3 axis, we mapped Ddit3 genomic occupancy by Cut&Run at 6 hr post-APAP—the peak of mid-zone Ddit3 activity during the injury-initiation phase. Compared with input control, we detected robust Ddit3-binding signals across chromatin domains (Figure 6A). Overlapping 6157 Ddit3-bound genes with the 47 mid-zone-specific DEGs identified 39 overlapping targets, termed Ddit3-binding DEGs (Figure 6B). Integrative Genomics Viewer (IGV) tracks confirmed Ddit3 occupancy at regulatory elements of Btg2, Atf3, Dnajb1, and Egr1 (Figure 6C).

Figure 6 with 2 supplements see all
Download asset Open asset
The Atf4-Ddit3 axis protects the liver by pausing proliferation via Btg2.

(A) Heatmap of Ddit3-binding signals in chromatin domains. (B) Overlap of Ddit3-binding genes by Cut&Run with differentially expressed genes (DEGs) in mid-zone was shown. These genes were referred to as Ddit3-binding DEGs. (C) Integrative Genomics Viewer (IGV) plot shows the Cut&Run peaks of Ddit3 and input on Btg2-site, Ifrd1-site, Fos-site, and Egr1-site. (D) Immunohistochemical staining of Btg2 in liver sections at 0 and 6 hr post-acetaminophen (APAP). Red arrows indicate Btg2-positive hepatocytes. (E) Schematic figure illustrating overexpression of Ddit3 via AAV in hepatocytes of wild-type C57BL/6J mice. Overexpression of EGFP is used as a control. (F) Immunohistochemical staining of Btg2 in liver sections from mice overexpressing either AAV-TBG-EGFP or AAV-TBG-Ddit3 at 12 hr post-APAP. Btg2-positive hepatocytes are indicated by red arrows. Btg2-positive hepatocytes in liver sections are quantified. n=4 mice/group. (G) Schematic figure illustrating overexpression of Btg2 via AAV in hepatocytes of wild-type C57BL/6J mice. Overexpression of EGFP is used as a control. (H) H&E staining shows morphology of livers from AAV-TBG-EGFP or AAV-TBG-Btg2 overexpressing mice at 6 hr post-APAP. Injured area is outlined by black dashed lines. The percentage of injury area is quantified. n=3 mice/group. Necrotic areas outlined by loss of cellular architecture on H&E. (I) Serum levels of ALT are measured at 6 hr post-APAP in AAV-TBG-EGFP or AAV-TBG-Btg2 overexpressing mice. n=3 mice/group. (J) Immunohistochemical staining of Ki67 in liver sections from mice overexpressing either AAV-TBG-EGFP or AAV-TBG-Btg2 at 72 hr post-APAP. Ki67-positive hepatocytes are indicated by red arrows. Ki67-positive cells in liver sections are quantified. n=4 mice/group. (K) Schematic figure illustrating knockdown of Btg2 via AAV in hepatocytes of wild-type C57BL/6J mice. IP479-EFS-sgEGFP is used as a control. (L) Serum levels of ALT are measured at 6 hr post-APAP in IP479-EFS-sgEGFP or IP479-EFS-sgBtg2 knockdown mice. n=3 mice/group. (M) H&E staining shows morphology of livers from IP479-EFS-sgEGFP or IP479-EFS-sgBtg2 knockdown mice at 6 hr post-APAP. Injured area is outlined by black dashed lines. The percentage of injury area is quantified. n=3 mice/group. Necrotic areas outlined by loss of cellular architecture on H&E. (N) Immunohistochemical staining of Ki67 in liver sections from IP479-EFS-sgEGFP or IP479-EFS-sgBtg2 knockdown mice at 72 hr post-APAP. Ki67-positive cells in liver sections are quantified. n=3 mice/group. Data are represented as means ± SD; unpaired Student’s t-test (F, H–N). **p<0.01; ***p<0.001; ****p<0.0001; ns, not significant.

Figure 6—source data 1

Numerical data of Figure 6F, H–J, and L–N.

https://cdn.elifesciences.org/articles/108646/elife-108646-fig6-data1-v1.xlsx

Among the Ddit3-binding DEGs, B-cell translocation gene 2 (Btg2) is a cell-cycle inhibitor (Stupfler et al., 2016; Park et al., 2008; Yuniati et al., 2019). Moreover, Btg2 is also co-regulated by Ddit3 and Atf4. Immunohistochemical staining further revealed increased Btg2 protein expression around the injured area at 6 hr post-APAP (Figure 6D). To determine whether Ddit3 drives Btg2 expression, we overexpressed Ddit3 in hepatocytes via AAV8-TBG (Figure 6E). Ddit3 overexpression significantly elevated Btg2 protein levels at 12 hr post-APAP compared to EGFP controls (Figure 6F). We next overexpressed Btg2 directly (Figure 6G). Btg2 overexpression significantly reduced APAP-induced injury, as shown by decreased necrotic areas and serum ALT levels (Figure 6H and I) and markedly reduced Ki67+ hepatocyte numbers (Figure 6J).

To rule out potential confounds associated with AAV transduction, we measured Cyp2e1 protein levels, which showed only a modest decrease in both Atf4- and Btg2-overexpressing groups—well below the threshold known to confer significant protection (Figure 6—figure supplement 1A; Ganetsky et al., 2019; Ye et al., 2022). AAV transduction efficiencies averaged 32% (EGFP), 18% (Atf4), and 23% (Btg2), with individual efficiency negatively correlating with ALT levels (e.g. Pearson r = –0.7681, p=0.0260 for Atf4; Figure 6—figure supplement 1B and C). These controls confirm that the protective effects are attributable to transgene expression rather than nonspecific Cyp2e1 inhibition. Forced Btg2 expression phenocopies the protective, proliferation-suppressive state observed during early injury, positioning Btg2 as a critical downstream mediator of the stress response.

To determine whether Btg2 is necessary for the protective arrest, we knocked down Btg2 using the AAV8-CasRx-sgRNA system (Figure 6K; Zhou et al., 2020). Despite a moderate (~30%) reduction in Btg2 expression (Figure 6—figure supplement 1D), knockdown resulted in a significant ~2-fold increase in serum ALT, ~1.5-fold expansion of necrotic areas, and ~1.8-fold increase in Ki67+ hepatocytes compared to sgEGFP controls (Figure 6L–N). These data suggest that loss of Btg2 abrogates the protective brake on proliferation, sensitizing the liver to injury and permitting unscheduled cell-cycle activity under stress conditions.

To assess the broader applicability of this mechanism, we tested a carbon tetrachloride (CCl4)-induced acute liver injury model, which causes centrilobular necrosis through free-radical pathways rather than cytochrome P450-dependent metabolite formation (Wang et al., 2024). At 18 hr post-CCl4 injection, we observed elevated p-eIF2α, Atf4, Chop, and Btg2 levels alongside reduced Ki67 expression near injury sites (Figure 6—figure supplement 2A–G). This conservation across mechanistically distinct toxicants suggests that ISR-mediated Btg2 induction and subsequent proliferation arrest represent a general adaptive response to centrilobular injury. Collectively, these data establish a conserved functional cascade in which the Atf4-Ddit3 axis induces Btg2 to enforce transient proliferative restraint, prioritizing survival during early acute liver injury and preserving the capacity for subsequent regeneration.

Discussion

Our study supports a model that during the early phase of APAP-induced liver injury, mid-zone hepatocytes undergo a transient proliferative restraint orchestrated by a distinctive transcriptional program centered on the Atf4-Chop-Btg2 axis. ST and functional assays demonstrate that this adaptive response promotes survival by upregulating the cell-cycle inhibitor Btg2, highlighting a trade-off between cytoprotection and proliferative capacity. Notably, zonal metabolic heterogeneity in APAP processing drives preferential activation of these stress pathways in mid-zone hepatocytes. These findings uncover a critical mechanism by which hepatocytes prioritize stress adaptation before initiating regeneration, offering new insights into the spatiotemporal regulation of liver repair (Figure 7).

Figure 7
Download asset Open asset
The Atf4-Ddit3 axis mediates integrated stress protection in mid-zone hepatocytes at the expense of proliferative capacity during early AILI.

Schematic model depicting the adaptive response of mid-zone hepatocytes to acute liver injury, wherein activation of the Atf4-Ddit3 pathway upregulates Btg2 to induce transient cell-cycle arrest, balancing survival and regeneration.

While the ISR shares upstream initiators with ER stress (e.g. eIF2α phosphorylation), its outcomes are context-dependent. ER stress specifically activates the UPR through Ire1α, Atf6, and Perk (Wang et al., 2006), with genetic studies demonstrating that knockout of Ire1α, Xbp1, or Ddit3 consistently attenuates liver damage, confirming the UPR’s contribution to hepatotoxicity (Nagy et al., 2007; Uzi et al., 2013; Hur et al., 2012; Ye et al., 2022). Unlike classical ER stress, we find that mid-zone hepatocytes leverage the Atf4-Chop axis—defined here as a functional pathway of co-regulated transcriptional targets—to transiently halt proliferation via Btg2, favoring survival. This aligns with emerging evidence that Ddit3 can paradoxically support adaptation in mitochondrial stress (Kaspar et al., 2021; Kreß et al., 2023), suggesting that Ddit3 plays a context-dependent role in fine-tuning the ISR in mammals.

Our work resolves an early adaptive phase (3–12 hr) of proliferation arrest that is temporally and functionally distinct from previously described zonal regeneration programs, including later perinecrotic proliferation (>24 hr) (Bajt et al., 2003), CXCR2+ pro-proliferative responses (Nguyen et al., 2022), and p21-mediated senescence (Umbaugh et al., 2024). Rather than acting as transitional cells, mid-zone hepatocytes deploy the Atf4-Chop-Btg2 axis to transiently redirect resources from cell division toward cytoprotective pathways, such as redox defense. This spatially restricted mechanism operates before overt necrosis expansion and, unlike earlier PP/PC stress profiles (Umbaugh et al., 2021), is causally supported by integrated ST and functional validation (Cut&Run, AAV modulation). Our findings therefore bridge early stress adaptation with subsequent liver repair, revealing that mid-zone hepatocytes may preserve regenerative capacity by prioritizing survival over proliferation during the critical injury-initiation window.

The antiproliferative activity of Btg2 has been studied through an integrated network of transcriptional and post-transcriptional regulatory mechanisms. At the post-transcriptional level, Btg2 bridges poly(A)-binding protein cytoplasmic 1 (Pabpc1) with the CCR4-NOT deadenylase complex, accelerating degradation of cell cycle-related transcripts (Stupfler et al., 2016). Transcriptionally, Btg2 disrupts the cyclin B1-Foxm1-positive feedback loop, suppressing cyclin B1 expression and inducing G2/M arrest (Park et al., 2008). Btg2 also inhibits G1/S transition through cyclin-dependent kinase interactions and modulates p53 activity (Yuniati et al., 2019). In the context of APAP-induced liver injury, the stress-responsive Atf4-Chop pathway induces Btg2 expression in Mid hepatocytes, where it likely coordinates a comprehensive proliferation arrest through these complementary mechanisms. This multilayered regulatory strategy enables hepatocytes to temporarily pause cell proliferation and redirect resources toward cellular repair and stress adaptation during the acute phase of toxic injury.

A key limitation of this study is the scarcity of early-phase APAP-induced liver injury human samples, which precludes definitive validation of the clinical relevance of these zonal adaptations. The apparent differences between the human and murine datasets are likely attributable to the stage of liver injury analyzed. Human samples were obtained from end-stage liver explants, in which regenerative programs are already active, whereas our mouse studies focused on the early (3–12 hr) injury phase preceding overt regeneration. Consistent with this interpretation, one ST specimen displayed a stress-dominant Mid transcriptional program resembling our murine observations, whereas the second exhibited a predominantly proliferative profile. These findings suggest that the balance between stress adaptation and hepatocyte proliferation is dynamic and strongly influenced by the temporal stage of liver injury. Although the available human datasets do not capture the early injury window examined in our mouse model, they support conservation of the Mid stress-response program across species.

We also acknowledge several technical limitations of our ST approach. While ST analysis and Ki67 staining consistently identified hepatocytes as the primary proliferating compartment during early APAP-induced liver injury, non-parenchymal cells (e.g. immune cells, stellate cells) may contribute to division signals. Single-cell resolution techniques would be required to exclude these contributions fully. Additionally, our ST data alone do not independently confirm zonal proliferative restraint: ST lacks overt separation in S/G2/M scores, uses non-canonical markers (e.g. Nasp, Cks1b), and has low sensitivity for rare Ki67+ events (~1% of spots). Thus, our primary evidence for mid-zone proliferation dynamics comes from Ki67 immunohistochemistry (Figure 1), with ST findings (Figure 2—figure supplement 1I–J) shown as exploratory and directionally consistent with transient impairment, but not as definitive proof of proliferative restraint.

In summary, our study uncovers a spatially coordinated stress-adaptation mechanism in mid-zone hepatocytes during early APAP-induced liver injury, where the Atf4-Chop axis transiently suppresses proliferation via Btg2 upregulation. By elucidating how hepatocytes prioritize survival before initiating repair, our work provides a framework for understanding the dynamic interplay between stress adaptation and regeneration in tissues.

Materials and methods

Animal experiments and procedures

Animals

C57BL/6J mice (strain no. N000013) were used as WT mice. All mouse colonies were maintained at the Animal Core Facility of Yunnan University. The animal studies were approved by the Yunnan University Institutional Animal Care and Use Committee (IACUC, Approval No. YNU20220262). For APAP treatment, mice (8–12 weeks of age) were fasted overnight (5:00 pm to 9:00 am) before intraperitoneally (i.p.) injected with APAP (Sigma, A7085) at a dose of 300 mg/kg for male mice, as female mice are less susceptible to APAP-induced liver injury (Munoz and Fearon, 1984). Male mice have been the choice in the vast majority of the studies of APAP-induced liver injury reported in the literature (Shan et al., 2021; Campion et al., 2008; Ben-Moshe et al., 2022; Matchett et al., 2024). Therefore, we used male mice in the majority of the experiments presented. In one experiment, mice were i.p. injected with various doses of APAP (100, 300, 500 mg/kg). For carbon tetrachloride (CCl4) treatment, mice were i.p. injected with corn oil or CCl4 (Macklin, Cat #C805329) at a dose of 0.8 µl/g (corn oil diluted to 20%).

In some experiments, APAP-treated mice were pre-injected i.p. with either AAV-overexpression plasmid (AAV-TBG-Atf4/Ddit3/Btg2) or AAV-sgRNA knockdown plasmid (IP479-EFS-sgBtg2) for 3 weeks. AAV-TBG-EGFP and IP479-EFS-sgEGFP were used as controls, respectively. ALT measurement was performed using a diagnostic assay kit (Teco Diagnostics, Anaheim, CA, USA) to assess liver injury.

Histology, immunohistochemistry, and immunofluorescence

H&E staining

Request a detailed protocol

Tissues were fixed overnight at 4°C using buffered 10% paraformaldehyde (Sangon Biotech, A500684-0500) and subsequently embedded in paraffin. For H&E staining, paraffin-embedded slides were deparaffinized and rehydrated following standard protocols. The sections were briefly immersed in hematoxylin for 30 s, rinsed in water, and then stained with Eosin solution for 3 min. After staining, the sections were dehydrated in ethanol, cleared in xylene, and mounted. Slides were examined under a microscope (Scan System, SQS-1000) and analyzed using ImageJ. Necrotic areas outlined by loss of cellular architecture on H&E.

Immunofluorescence on frozen section

Request a detailed protocol

Immunofluorescence staining was conducted on frozen sections of fresh liver tissues. Initially, the tissues were fixed using 2% paraformaldehyde for 1 hr and subsequently dehydrated overnight in a 30% sucrose solution. The following day, tissue embedding was performed using OCT (SAKURA, 4583), after which ultra-thin sections of 5 μm were sliced. Permeabilization was achieved using 0.2% Triton X-100 for 10 min at room temperature, followed by blocking with 5% normal goat serum (VivaCell, C2530-0100). Primary antibodies against mouse GS (Proteintech, 11037-2-AP, 1:200), Chop (Proteintech, 15204-1-AP, 1:1000), Cyp1a2 (Proteintech, 15204-1-AP, 1:50), and Ki67 (Abcam, ab180569, 1:1000) were then applied. Secondary antibodies (Affinipure Goat Anti-Rabbit 488-conjugated, Jackson ImmunoResearch, 111-545-003, 1:400; Goat anti-Rabbit IgG (H+L) Cross-Adsorbed Secondary Antibody, Alexa Fluor 568, Invitrogen, A11011, 1:1000, Alexa Fluor 488-conjugated Affinipure Goat Anti-Mouse IgG+IgM(H+L), Jackson ImmunoResearch, 115-545-044, 1:400; and Alexa Fluor 594 AffiniPure Goat Anti-Mouse IgG (H+L), Jackson ImmunoResearch, 115-585-003, 1:200) were used accordingly. After washing, the slides were mounted using antifade medium, and nuclei were stained with DAPI (Beyotime, C1006, prediluted). Images were captured using an Olympus BP80 microscope and analyzed using ImageJ.

Immunohistochemistry

Request a detailed protocol

Livers were fixed in 10% buffered formalin for 24 hr before paraffin embedding. For immunohistochemical staining of paraffin-embedded liver tissue, 5 μm sections were prepared, deparaffinized, and subjected to antigen retrieval by microwaving for 20 min. Antigen retrieval was performed in pH 6.0 sodium citrate buffer (Ki-67, P62, Chop, Cyp1a2, Cyp2f2), pH 6.4 sodium citrate buffer (p-eIF2A), or pH 9.0 Tris-EDTA buffer (Atf4, Btg2). Sections were blocked with 3% hydrogen peroxide and permeabilized with 0.2% Triton X-100 for 10 min at room temperature, followed by blocking with 5% normal goat serum (VivaCell, C2530-0100) for 1 hr at room temperature. Primary antibodies were incubated overnight at 4°C at the following dilutions in blocking buffer: Ki-67 (Abcam, ab15580, 1:400), P62 (Abclonal, A19700, 1:1000), Cyp2f2 (Santa Cruz, sc-374540, 1:200), Chop (Proteintech, 15204-1-AP, 1:200), Cyp1a2 (Santa Cruz, sc-53241, 1:50), p-eIF2A (CST, 3398S, 1:50), Atf4 (CST, 11815S, 1:200), and Btg2 (Proteintech, 22339-1-AP, 1:100). Secondary antibodies used were Goat anti-Rabbit IgG (H+L) Secondary Antibody, Biotin (Invitrogen, 65-6140, 1:1000) and peroxidase-AffiniPure Goat Anti-mouse IgG(H+L) (Jackson ImmunoResearch, 115-035-003, 1:1000). Avidin, NeutrAvidin, Horseradish Peroxidase conjugate (Invitrogen, A2664, 1:2000) was used as the third antibody. After three washes with PBS, color development was achieved using the DAB Peroxidase Substrate Kit (ZSGB-BIO, ZLI-9018), followed by quenching in distilled water. Slides were counterstained with hematoxylin (BBI, A600701-0050), dehydrated to xylene, and images were captured using an Olympus BP80 microscope or Scan System SQS-1000 and analyzed using ImageJ.

Western blot analysis

Request a detailed protocol

Mouse liver tissue (30 mg) was homogenized in 1 ml RIPA buffer (10 mM Tris-HCl pH 8.0, 1 mM EDTA, 0.5 mM EGTA, 1% Triton X-100, 0.1% sodium deoxycholate, 0.1% SDS, 140 mM NaCl, 1 mM PMSF, and protease inhibitor cocktail). The lysate was mixed with 5× SDS loading buffer (final 1×), boiled at 98°C for 10 min, and centrifuged at 12,000 rpm for 3 min at room temperature. Supernatants were separated by 10% SDS-PAGE, transferred to PVDF membranes, and blocked with 5% nonfat milk in TBST (Tris-buffered saline with 0.1% Tween-20). Membranes were incubated with primary antibodies against Cyp2e1 (Proteintech, 19937-1-AP, 1:5000) and GAPDH (Proteintech, 60004-1-Ig, 1:5000), followed by HRP-conjugated secondary antibodies: goat anti-mouse (Jackson ImmunoResearch, 115-035-003, 1:2000) and goat anti-rabbit (Jackson ImmunoResearch, 111-035-003, 1:2000). Protein bands were detected using enhanced chemiluminescence and a Minichemi Chemiluminescence Imaging System.

RNA extraction and quantitative real-time PCR

Request a detailed protocol

Ground mouse liver tissue was lysed in TRIzol reagent (Invitrogen, 15596018) and centrifuged at 12,000 rpm for 10 min at 4°C. The supernatant was transferred to a new tube, mixed with 200 μl chloroform, incubated at room temperature for 10 min, and centrifuged. The aqueous phase was collected, mixed with an equal volume of isopropanol, incubated at 4°C for 10 min, and centrifuged. The supernatant was discarded, and the RNA pellet was washed twice with 75% ethanol, air-dried for 10 min, and dissolved in RNase-free water. RNA concentration and purity were measured using a NanoDrop spectrophotometer. Reverse transcription was performed using a cDNA first-strand synthesis kit (Takara, 6210B). Quantitative PCR was carried out in triplicates using PowerUp SYBR Green Master Mix (Thermo Fisher, A25742) on a QuantStudio 1 Real-Time PCR System (Life Technologies) according to the manufacturer’s instructions. Primer sequences are provided in Supplementary file 1.

AAV production and purification

Request a detailed protocol

AAV8 was produced using 293T cells cultured in one or more 175 cm2 cell culture flasks. Absence of mycoplasma contamination was confirmed using PCR. Cells were plated 1 day prior to transfection, reaching 50% confluence, which allowed them to reach 80–90% confluence the next day. For transfection of one 175 cm2 flask, a mixture of 11.56 μg pHelper, 7.24 μg pAAV2/8 (Addgene, 112864), and 5.9 μg transgene plasmids was prepared in 3 ml Opti-MEM medium in a tube. PEI solution (1 μg/μl in water, pH 7.0, powder from Polysciences, 23966-2) was added to the tube at a ratio of 3 μg PEI per 1 μg DNA. The solution was mixed and incubated for 15 min before adding to the cell culture. At 24 hr post-transfection, the medium was replaced with antibiotic-containing medium (P+/S+) without FBS. At 72 hr after the medium change, cells were scraped off the flask, and 3 ml of chloroform was added, followed by vortexing for 5 min. Then, 7.6 ml of 5 M NaCl was added and vortexed for 10 s. The solution was centrifuged at 3000×g for 5 min at 4°C. The supernatant was transferred to a new 50 ml tube, and 9.4 ml of 50% (wt/vol) PEG 8000 was added, vortexed for 10 s, and left for 1 hr at 4°C. Subsequently, the mixture was centrifuged at 3000×g for 30 min at 4°C. The supernatant was decanted, and the tube was inverted for 10 min to dry. Next, 1.4 ml of 50 mM HEPES buffer (pH 8.0) was added to resuspend the pellet, and the solution was vortexed for 5 min. Then, 3.5 μl of 1 M MgCl2, 2.8 μl of 5 μg/μl DNase I, and 2.8 μl of 5 μg/μl RNase A were added, and the mixture was incubated for 20 min at 37°C. The virus was then extracted with chloroform, and the solvent was replaced with PBS. AAV was concentrated by centrifugation at 14,000×g for 5 min at 4°C using a 100 kDa ultrafiltration tube (Beyotime, FUF058). Finally, the samples were aliquoted and stored at –80°C after titer measurement.

Library preparation for 10x Visium ST

Request a detailed protocol

After administering APAP, mice were euthanized at 0 hr, 3 hr, 6 hr, and 24 hr time points using carbon dioxide. Freshly harvested liver tissues were frozen using isopentane and liquid nitrogen, embedded in OCT, sectioned into 10 μm slices, and mounted on Visium tissue optimization or spatial gene expression slides. Tissues were fixed with methanol, stained with H&E, and imaged at ×40 magnification. Tissue permeabilization was optimized for 6 min. Libraries were prepared according to the Visium spatial gene expression user guide, with cDNA amplification cycles adjusted per time point, and index PCR performed for 6 cycles. Libraries were pooled, loaded at 400 pM, and sequenced on a NovaSeq 6000 SP100 sequencer.

Processing of raw sequencing data

Request a detailed protocol

The raw data underwent processing using space ranger (version 1.1.0) to perform several analytical tasks. These included detecting tissue boundaries, aligning reads to the reference genome, generating feature-spot matrices, conducting clustering to identify spatially distinct groups of cells, and analyzing gene expression patterns. Additionally, Spaceranger facilitated the spatial placement of spots within the context of the slide image, providing a comprehensive spatial transcriptomic analysis of the tissue sample.

Basic analysis of ST data

Request a detailed protocol

For the gene-spot matrices generated by Spaceranger, we utilized Scanpy (version 1.9.3) for analysis (Patro et al., 2017). Initially, routine statistical analyses were conducted, including calculating the number of detected UMIs (nUMI) and genes (nGene) in each spot. Basic quality control (QC) measures were then applied to the data. Specifically, spots with extremely low UMI or gene counts (min_counts = 2000, min_cell = 10) were excluded using scanpy.pp.filter_cells function.

Mitochondrial and hemoglobin genes were filtered out to focus on relevant gene expression. Spot matrix was filtered out to keep only spots overlaying tissue sections. Following QC, data normalization was performed using the scanpy.pp.normalize_total function, and log-transformation was applied using the scanpy.pp.log1p function. Subsequently, we identified 2000 highly variable genes using the scanpy.pp.highly_variable_genes function based on their expression means and variances.

Principal component analysis was then applied to reduce the dimensionality of the data, projecting the spots into a lower-dimensional space using the identified principal components. Using the corrected PC matrices, we conducted unsupervised clustering based on shared nearest neighbors and visualized the data using UMAP for further exploration and analysis (Lim and Qiu, 2023).

Each sample was partitioned into three clusters. To facilitate comparison, each cluster was annotated with a region label (PC, PP, or Mid), determined by integrating information from cluster-specific marker genes and H&E staining images. Furthermore, to compare gene expression profiles across clusters, we identified DEGs using fold change analysis and the Wilcoxon rank sum test. This approach allowed us to highlight genes showing significant expression differences among all or selected clusters, providing insights into biological variations across spatial regions in the tissue samples.

Analysis of DEGs and pathway enrichment

Request a detailed protocol

For differential gene analysis across different samples, we initially integrated all samples using the combat function to mitigate batch effects (Stuart et al., 2019). Subsequently, differential expression analysis of variable genes was conducted using the FindMarkers function, employing the Wilcoxon rank sum test (Ritchie et al., 2015). Genes showing significant changes were selected based on criteria, including a log-fold change threshold of 0.1 and a minimum percentage of samples showing differential expression of 1%. The DEGs identified were enriched across various time points post-APAP injection, relative to the 0 hr time point.

Pathway analysis was performed using the Kyoto Encyclopedia of Genes and Genomes (KEGG), while gene ontology (GO) analysis was carried out using the clusterProfiler package (Yu et al., 2012) to explore functional categories associated with these DEGs. The most significantly enriched pathways were identified using the enrichGO and enrichKEGG functions within the clusterProfiler package (Yu et al., 2012). Furthermore, gene set variation analysis (GSVA) was employed to investigate the expression variations related to proliferation at different time points following APAP administration, utilizing the gsva function from the GSVA R package (Hänzelmann et al., 2013). Pathway activities were visualized using the pheatmap R package for comprehensive analysis and interpretation. To functionally classify the DEGs identified in the F2A comparison, we performed GO enrichment analysis using clusterProfiler. Each DEG was mapped to GO terms from the three ontologies: biological process, cellular component, and molecular function. When a single gene was annotated with multiple GO terms, we selected the most prominent term—defined as the one with the smallest adjusted p-value or highest enrichment score—for downstream interpretation. This approach ensured that the primary functional role of each DEG was captured without redundancy.

GRN analysis

Request a detailed protocol

The SCENIC workflow comprises three main steps: co-expression analysis, target gene motif enrichment analysis, and regulon activity assessment. The analysis employed pySCENIC (version 0.12.1) (Van de Sande et al., 2020) with default parameters, using the raw count matrix from all samples as input. Co-expression modules were initially identified, and the interaction strength between TFs and their target genes was evaluated using GRNBoost (Aibar et al., 2017). Each co-expressed module underwent motif enrichment analysis using RcisTarget with a rank threshold of 3000, to identify modules where the TF motif was significantly enriched among targets (Aibar et al., 2017). Only modules meeting these criteria were retained, establishing TFs along with their potential direct targets as regulons. The activity of each regulon within each cell was assessed using AUCell. For visualization, average regulon activity scores (AUC) were calculated for each cluster (Aibar et al., 2017). A rank plot of regulons was generated using ggplot2. Additionally, the gene expression network related to the target genes of the TFs of interest was visualized using the nx.draw function.

Analysis of snRNA-seq data (GSE223561)

Request a detailed protocol

To annotate cell clusters in the snRNA-seq data, we utilized a curated list of known marker genes corresponding to liver cell lineages, as reported in Multimodal Decoding of Human Liver Regeneration. Using the AddModuleScore function in Seurat, we calculated signature scores to assign cell lineage identities. Clusters identified as primarily composed of cycling cells were further reclustered to segregate them into their constituent lineages. This process was repeated iteratively for each identified lineage. A cleansing step was included to remove clusters enriched for nuclei identified as doublets or overexpressing marker genes from multiple lineages. Hepatocytes were subsequently extracted for further reclustering, using resolution (res)=0.2 and 80 principal components (npc). To annotate the hepatocyte clusters, we applied a curated list of known zonation-specific marker genes for hepatocytes.

Library preparation for Cut&Run samples

Request a detailed protocol

Cut&Run was conducted on mouse liver tissue samples from mice injected with APAP for 6 hr using the Cut&Run assay kit (CST, 86652) following the manufacturer’s protocol. Initially, 1 mg of fresh tissue was weighed for each antibody/MNase reaction, with an additional 1 mg for the input sample. Reactions included positive and negative controls. Tissues were fixed with 1 ml of 0.1% formaldehyde for 10 min and quenched with glycine for 5 min. After washing with PBS, tissues were resuspended in 1 ml of 1× Wash Buffer (+spermidine+protease inhibitor cocktail) and transferred to a Dounce homogenizer. Tissue was homogenized into a single-cell suspension with 10–15 strokes until no tissue chunks remained.

For the Cut&Run assay to identify Chop-binding sites, samples were bound to concanavalin A beads and incubated overnight at 4°C with the following primary antibodies: 0.5 μg CHOP (CST, 2895S, 1:100), 2 μl positive control Tri-Methyl-Histone H3 (Lys4) (CST, 9751), and 5 μl negative control Rabbit (DA1E) mAb IgG XP Isotype Control (CST, 66362). The pAG MNase enzyme was activated by adding cold calcium chloride and incubating at 4°C for 30 min. The reaction was terminated with 1× stop buffer at 37°C for 10 min without shaking to release DNA fragments into the solution.

To reverse crosslinks in fixed tissue samples, samples were brought to room temperature and treated with 10% SDS Solution (BBI, B548118-0100) to a final concentration of 0.1%, along with proteinase K (20 mg/ml). DNA Extraction Buffer (+Proteinase K+RNAse A) was added to the input sample. Cells were lysed and chromatin fragmented by sonication using a Covaris ME220 contact ultrasound apparatus on ice, with optimal conditions generating chromatin fragments ranging from 100 to 600 bp. DNA was purified using phenol/chloroform extraction followed by ethanol precipitation, and DNA quantification was performed by qPCR. For Cut&Run, purified DNA was used to prepare sequencing libraries with the NEB Next Ultra II DNA Library Prep Kit for Illumina (7645S and 7335L) and sequenced on an Illumina Nova-Seq 6000 sequencer to obtain 150 bp paired-end reads.

Cut&Run data analysis

Request a detailed protocol

For visualization purposes, BW files were generated using IGV software (Robinson et al., 2011), and bigwig files were created using Deeptools (version 3.5.1) with the ‘bamCoverage’ module (bamCoverage --binSize 100) (Ramírez et al., 2016). For each biological replicate and its corresponding IgG control, peaks were called using macs2 (version 2.2.7) with the command macs2 callpeak -f BAMPE --qvalue 0.1 --keep-dup (Zhang et al., 2008). To annotate the identified peaks, ChIPseeker (version 1.26.2) in R was utilized (Yu et al., 2015).

Quantification of zonal distribution

Request a detailed protocol

The zonation of protein (such as Ki67, Chop, Atf4)-positive nuclei was quantified using the following method: the position index (P.I.) was calculated based on distances to the nearest CV (x), PV (y), and the distance between CV and PV (z), utilizing the law of cosines. The formula employed was P.I. = (x2+z2 – y2)/(2z2). This approach aligns with the methodology described by Lin et al., 2018.

Statistical analysis

Request a detailed protocol

All experimental data are presented as the mean ± SD. Statistical analyses were performed using GraphPad Prism (version 8.0.2), with appropriate tests selected based on experimental design: unpaired two-tailed Student’s t-tests for comparisons between two groups, and one-way ANOVA for comparisons involving three or more groups. The number of animals (‘n’) used in each experiment is indicated in the Figures and corresponding legends. For quantification of liver sections, three to five random PC and PP fields of each liver sample, unless specified otherwise, were imaged and quantified using ImageJ.

Appendix 1

Appendix 1—key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
AntibodyRabbit polyclonal anti-Ki67 antibodyAbcamCat# ab15580; RRID:AB_443209IF (1:1000) IHC (1:400)
AntibodyRabbit Glutamine Synthetase Polyclonal antibodyProteintechCat# 11037–2-AP; RRID:AB_2110650IF (1:200)
AntibodyMouse monoclonal Anti-Cyp1a2 antibodySantaCruzCat# sc-53241; RRID:AB_629359IF (1:50) IHC (1:50)
AntibodyRabbi CHOP; GADD153 Polyclonal antibodyProteintechCat# 15204–1-AP; RRID:AB_2292610IF (1:1000) IHC (1:200)
AntibodyMonoclonal ATF-4 (D4B8) Rabbit mAbCell Signalling TechnologyCat# 11815 S; RRID:AB_2616025IHC (1:200)
AntibodyRabbit BTG2 Polyclonal antibodyProteintechCat# 22339–1-AP; RRID:AB_2879078IHC (1:100)
AntibodyPolyclonal Phospho-eIF2α (Ser51) (D9G8) XP Rabbit mAbCell Signalling TechnologyCat# 3398 s; RRID:AB_2096481IHC (1:50)
AntibodyGoat anti-Rabbit IgG (H+L) Secondary Antibody, BiotinInvitrogenCat# 65–6140; RRID:AB_2533969IHC (1:1000)
AntibodyPeroxidase-conjugated Affinipure Goat Anti-Mouse IgG(H+L)Jackson ImmunoResearchCat# 115-035-003; RRID:AB_10015289IHC (1:1000) WB (1:2000)
AntibodyPeroxidase-conjugated Affinipure Goat Anti-Rabbit IgG(H+L)Jackson ImmunoResearchCat# 111-035-003; RRID:AB_2313567WB (1:2000)
AntibodyAvidin, NeutrAvidin, Horseradish Peroxidase conjugateInvitrogenCat# A2664IHC (1:2000)
AntibodyAlexa Fluor 488-conjugated Affinipure Goat Anti-Mouse IgG +IgM(H+L)Jackson ImmunoResearchCat# 115-545-044; RRID:AB_2338844IHC (1:400)
AntibodyAlexa Fluor 594 AffiniPure Goat Anti-Mouse IgG (H+L)Jackson ImmunoResearchCat# 115-585-003; RRID:AB_2338871IHC (1:200)
AntibodyAlexa fluor 568-goat anti-rabbit IgG(H+L)cross-adsorbed secondary antibodyInvitrogenCat# A11011; RRID:AB_143157IHC (1:1000)
AntibodyAlexa Fluor 488 AffiniPure Goat Anti-Rabbit IgG (H+L)Jackson ImmunoResearchCat# 111-545-003; RRID:AB_2338046IHC (1:400)
AntibodyCYP2E1-Specific Polyclonal antibodyProteintechCat# 19937–1-AP; RRID:AB_10646444WB (1:5000)
AntibodyGAPDH Monoclonal antibodyProteintechCat# 60004–1-Ig; RRID:AB_2107436WB (1:5000)
Chemical compound, drugAcetaminophen(APAP)SigmaCat# A7085
Chemical compound, drugParaformaldehyde, (CH2O)nSangon BiotechCat# A500684-0500
Chemical compound, drugPhosphate Buffered SalineVivaCellCat# C3580-0500
Chemical compound, drug10% Neutral Formalin Fix SolutionSangon BiotechCat# E672001-0500
Chemical compound, drugSucrose, C12H22O11Sangon BiotechCat# A502792-0005
Chemical compound, drugHigh effect paraffinceresinShanghai Hualing Rehabilitation Equipment Manufacturing Plant.Cat# N/A
Chemical compound, drugXyleneTianjin Zhiyuan Chemical Reagents Co., Ltd.Cat# N/A
Chemical compound, drugEthanolTianjin Zhiyuan Chemical Reagents Co., Ltd.Cat# N/A
Chemical compound, drugCitrate Antigen Retrieval Solution (Powder)Sangon BiotechCat# E673002-0020
Chemical compound, drugTrisSolarbioCat# T8060
Chemical compound, drugDisodium salt dihydrate, C10H14N2O8Na2·2H2O(EDTA)Sangon BiotechCat# A500838-0500
Chemical compound, drugHydrogen peroxideCat# N/A
Chemical compound, drugTriton X-100Sangon BiotechCat# A600198-0500
Chemical compound, drugGoat serumVivaCellCat# C2530-0100
Chemical compound, drugHematoxylinSangon BiotechCat# A600701-0050
Chemical compound, drugEosin Y(water soluble)AladdinCat# E141405
Chemical compound, drugNeutral balsamSolarbioCat# G8590
Chemical compound, drugTissue-tek OCT compoundSAKURACat# REF:4583
Chemical compound, drugAcetoneChron ChemicalsCat# N/A
Chemical compound, drugTween20Sangon BiotechCat# A600560-0500
Chemical compound, drugDAPI Staining SolutionBeyotimeCat# C1006
Chemical compound, drugIsopentaneAladdinCat# M108171
Chemical compound, drugMethanolTianjin Zhiyuan Chemical Reagents Co., Ltd.Cat# N/A
Chemical compound, drug20 X TBS bufferSangon BiotechCat# B548105-0500
Chemical compound, drugUltraPure DNase/RNase-Free Distilled WaterInvitrogenCat# 10977015
Chemical compound, drugTrichloromethaneChron ChemicalsCat# N/A
Chemical compound, drugFBSVivaCellCat# C04001-500
Chemical compound, drugDMEM(High glucose)VivaCellCat# C3113-0500
Chemical compound, drugPenicillin-Streptomycin SolutionVivaCellCat# C3421-0100
Chemical compound, drugTrypsin 1:300 from Porcine pancreasSangon BiotechCat# A100260-0050
Chemical compound, drugPEIPolysciencesCat# 23966–2
Chemical compound, drugOpti-MEMGibcoCat# 11058021
Chemical compound, drugPEG-8000Sangon BiotechCat# A600433-0500
Chemical compound, drugMgCl2GhtechCat# N/A
Chemical compound, drug1 M HEPESSolarbioCat# H1095
Chemical compound, drug37% formaldehydeSigmaCat# F8775
Chemical compound, drug10% SDS SolutionSangon BiotechCat# B548118-0100
Chemical compound, drugGlycineSangon BiotechCat# A502065-0005
Chemical compound, drugPierce Protease inhibitor tablets, EDTA freeThermoCat# A32965
Chemical compound, drugDNA extract bufferSolarbioCat# P1012
Chemical compound, drug3 M Sodium acetateInvitrogenCat# AM9740
Chemical compound, drugGlycogenThermoCat# R0561
Chemical compound, drugSSC buffer 20 X concentrateSigmaCat# S6639
Chemical compound, drugPowerUp SYBR Green Master MixApplied biosystemsCat# A25742
Chemical compound, drugProteinase K Solution (20 mg/ml)Sangon BiotechCat# B600169-0002
Chemical compound, drugDNaseⅠ,RNase-freeThermoCat# EN0521
Chemical compound, drugRnase AThermoCat# R1253
Chemical compound, drugAgencourt AMPure XP magnetic beadsBeckman CoulterCat# A63880
Chemical compound, drugEthyl alcohol, Pure (200 proof, molecular biologygrade)Sigma-AldrichCat# E7023-500ML
Chemical compound, drugDEPC-treated waterBiosharpCat# 701062
Commercial assay or kitVisium Spatial Gene Expression Slide & Reagent Kit, 4 rxns10 X GenomicsCat# PN-1000187
Commercial assay or kitVisium Accessory Kit10 X GenomicsCat# PN-1000194
Commercial assay or kitDual Index TT Set APN10 X GenomicsCat# PN-1000215
Commercial assay or kitAlanine aminotransferase assay kitNanjing Jiancheng Bioengineering InstituteCat# c009-2-1
Commercial assay or kitAspartate aminotransferase Assay KitNanjing Jiancheng Bioengineering InstituteCat# C010-2-1
Commercial assay or kitCUT&RUN Assay KitCell Signalling TechnologyCat# 86652
Commercial assay or kitNEBNext Ultra II DNA Library Prep Kit for IlluminaNew England BiolabsCat# 7645 S
Commercial assay or kitNEBNext Multiplex Oligos for Illumina (Index Primers Set 1)New England BiolabsCat# 7335 L
Commercial assay or kitQubit dsDNA HS Assay KitThermo Fisher ScientificCat# Q32851
Commercial assay or kit2×Taq Master Mix (Dye Plus)VazymeCat# P112-01
Commercial assay or kit2×Rapid Taq Master MixVazymeCat# P222-01
Cell line (Homo-sapiens)HEK293TCell Bank of the Chinese Academy of SciencesATCC# CRL-3216; RRID:CVCL_0063
Recombinant DNA reagentPlasmid AAV8-TBG-EGFPThis paperpAAV2/8-based EGFP overexpression plasmid
Recombinant DNA reagentPlasmid AAV8-TBG-Atf4This paperpAAV2/8-based Atf4 overexpression plasmid
Recombinant DNA reagentPlasmid AAV8-TBG-Ddit3This paperpAAV2/8-based Ddit3 overexpression plasmid
Recombinant DNA reagentPlasmid AAV8-TBG-Btg2This paperpAAV2/8-based Btg2 overexpression plasmid
Recombinant DNA reagentPlasmid IP479-EFS-sgEGFPThis paperan EFS promoter-driven sgRNA expression plasmid targeting EGFP
Recombinant DNA reagentPlasmid IP479-EFS-sgBtg2This paperan EFS promoter-driven sgRNA expression plasmid targeting Btg2
Strain, strain background (Mus musculus, male)C57BL/6 JGemPharmate ch Co., Ltd.RRID:IMSR_GPT:N000013
Software and algorithmsGraphPad PrismGraphPad SoftwareRRID:SCR002798https://www.graphpad.com
Software and algorithmsSPSSIBM SPSS softwareRRID:SCR_002865https://www.ibm.com/
Software and algorithmsImageJNational Institutes of HealthRRID:SCR_003070https://imagej.nih.gov/ij/
Software and algorithms10 x Genomics Space Ranger (v1.1.0)10 x GenomicsRRID:SCR_025848https://10xgenomics.com
Software and algorithmsClusterProfiler R package (v4.1.4)Yu et al.RRID:SCR_016884https://bioconductor.org/packages/release/bioc/html/clusterProfiler.html
Software and algorithmsSeurat R package (v4.2.3)Stuart et al.RRID:SCR_016341https://satijalab.org/seurat/
Software and algorithmsBowtie2 (v2.2.5)Langmead et al.RRID:SCR_016368http://bowtie-bio.sourceforge.net/bowtie2/index.shtml
Software and algorithmsDeeptools (v3.5.1)Ramirez et al.RRID:SCR_016366https://deeptools.readthedocs.io/
Software and algorithmsIntegrative Genomics Viewer (IGV)Robinson et al.RRID:SCR_011793https://igv.org/
Software and algorithmsR (v4.0.5)and RstudioR ConsortiumRRID:SCR_000432https://www.rstudio.com/
Software and algorithmsPython(v3.7.1)PythonRRID:SCR_008394https://www.python.org/
Software and algorithmsGSVA R package (v1.38.2)Hanzelmann et al.RRID:SCR_021058https://www.bioconductor.org/packages/release/bioc/html/GSVA.html
Software and algorithmsLimma R package (v3.46.0)Ritchie et al.RRID:SCR_010943https://www.bioconductor.org/packages/release/bioc/html/limma.html
Software and algorithmsPheatmap R package (v.1.0.12)N/ARRID:SCR_016418https://ggplot2.tidyverse.org
Software and algorithmsGgplot2 R packageN/ARRID:SCR_014601https://ggplot2.tidyverse.org
Software and algorithmsScanpy (python package)Wolf et al.RRID:SCR_018139https://scanpy.readthedocs.io
Software and algorithmspySCENIC (version 0.12.1)Van de Sande B et al.RRID:SCR_025802https://scenic.aertslab.org
OtherSpatial transcriptomics dataThis paperGSE272564https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE272564
OtherCut&Run dataThis paperGSE272565https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE272565

Data availability

All data generated or analysed during this study are included in the manuscript and supporting files; source data files have been provided. The raw spatial transcriptomics data of mouse liver at different time points after APAP injection have been deposited in the GEO database (accession number: GSE272564). The raw Cut&Run data of mouse liver at 6 hours after APAP injection have been deposited in the GEO database (accession number: GSE272565).

The following data sets were generated
    1. Zhu Y
    2. Deng C
    3. Chen B
    4. He J
    (2024) NCBI Gene Expression Omnibus
    ID GSE272565. Mid-Lobular Hepatocytes Exhibit Transient Proliferation Pause and Stress Response through Atf4-Ddit3 axis in Early Acute Liver Injury (CUT&Run).
    https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE272565
    1. Zhu Y
    2. Deng C
    3. Chen B
    4. He J
    (2024) NCBI Gene Expression Omnibus
    ID GSE272564. Mid-Lobular Hepatocytes Exhibit Transient Proliferation Pause and Stress Response through Atf4-Ddit3 axis in Early Acute Liver Injury (10x Visium).
    https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE272564
The following previously published data sets were used

References

    1. Bernal W
    2. Wendon J
    (2013) Acute liver failure
    New England Journal of Medicine 369:2525–2534.
    https://doi.org/10.1056/NEJMra1208937
  17. Book
    1. Jack A
    2. Hinson DWR
    3. James LP
    (2010) Mechanism of acetaminophen-induced liver necrosis
    In: Uetrecht J, editors. Handbook of Experimental Pharmacology. Springer. pp. 369–405.
    https://doi.org/10.1007/978-3-642-00663-0_12

Article and author information

Author details

  1. Yaying Zhu

    Yunnan Key Laboratory of Cell Metabolism and Diseases, Center for Life Sciences, School of Life Sciences, Yunnan University, Kunming, China
    Contribution
    Conceptualization, Data curation, Software, Formal analysis, Supervision, Validation, Investigation, Visualization, Methodology, Writing – original draft, Project administration, Writing – review and editing
    Contributed equally with
    Chengxiang Deng and Bo Chen
    Competing interests
    No competing interests declared

  2. Chengxiang Deng

    Yunnan Key Laboratory of Cell Metabolism and Diseases, Center for Life Sciences, School of Life Sciences, Yunnan University, Kunming, China
    Contribution
    Conceptualization, Data curation, Software, Formal analysis, Investigation, Visualization, Methodology
    Contributed equally with
    Yaying Zhu and Bo Chen
    Competing interests
    No competing interests declared

  3. Bo Chen

    Yunnan Key Laboratory of Cell Metabolism and Diseases, Center for Life Sciences, School of Life Sciences, Yunnan University, Kunming, China
    Contribution
    Data curation, Software, Formal analysis, Visualization, Writing – review and editing
    Contributed equally with
    Yaying Zhu and Chengxiang Deng
    Competing interests
    No competing interests declared

  4. Jia He

    Yunnan Key Laboratory of Cell Metabolism and Diseases, Center for Life Sciences, School of Life Sciences, Yunnan University, Kunming, China
    Contribution
    Software, Investigation, Visualization, Methodology
    Competing interests
    No competing interests declared

  5. Yanan Liu

    Yunnan Key Laboratory of Cell Metabolism and Diseases, Center for Life Sciences, School of Life Sciences, Yunnan University, Kunming, China
    Contribution
    Data curation, Investigation
    Competing interests
    No competing interests declared

  6. Shan Lei

    Yunnan Key Laboratory of Cell Metabolism and Diseases, Center for Life Sciences, School of Life Sciences, Yunnan University, Kunming, China
    Contribution
    Investigation
    Competing interests
    No competing interests declared

  7. Weiju Lu

    Yunnan Key Laboratory of Cell Metabolism and Diseases, Center for Life Sciences, School of Life Sciences, Yunnan University, Kunming, China
    Contribution
    Software, Visualization
    Competing interests
    No competing interests declared

  8. Cheng Peng

    Yunnan Key Laboratory of Cell Metabolism and Diseases, Center for Life Sciences, School of Life Sciences, Yunnan University, Kunming, China
    Contribution
    Resources, Supervision, Funding acquisition, Methodology, Writing – review and editing
    For correspondence
    chengpeng@ynu.edu.cn
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-0448-0094

  9. Zhao Shan

    Yunnan Key Laboratory of Cell Metabolism and Diseases, Center for Life Sciences, School of Life Sciences, Yunnan University, Kunming, China
    Contribution
    Conceptualization, Resources, Supervision, Funding acquisition, Methodology, Writing – original draft, Project administration, Writing – review and editing
    For correspondence
    shanzhaolab@163.com
    Competing interests
    Reviewing editor, eLife
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-5064-1023

Funding

National Natural Science Foundation of China (82570734)

  • Zhao Shan

Yunnan Provincial Science and Technology Department (C619300A086)

  • Zhao Shan

National Natural Science Foundation of China (32170662)

  • Cheng Peng

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Bin Qi (Yunnan University) for suggestions and discussion. We thank Dr. Yin Hao (Wuhan University) for providing us the pHelper and pAAV2/8 plasmids. We thank Hui Yang (Institute of Neuroscience, SBS, CAS) for providing us the plasmids related to the CasRx system. We thank Yonglong Wei (Yunnan University) for helping in the quantification of spatial distribution. This work was supported by the National Natural Science Foundation of China (82570734 to ZS), Yunnan Provincial Science and Technology Department (C619300A086 to ZS), National Natural Science Foundation of China (32170662 to CP).

Ethics

The animal studies were approved by the Yunnan University Institutional Animal Care and Use Committee (IACUC, Approval No.YNU20220262).

Version history

  1. Sent for peer review:
  2. Preprint posted:
  3. Reviewed Preprint version 1:
  4. Reviewed Preprint version 2:
  5. Reviewed Preprint version 3:
  6. Version of Record published:

Cite all versions

You can cite all versions using the DOI https://doi.org/10.7554/eLife.108646. This DOI represents all versions, and will always resolve to the latest one.

Copyright

© 2025, Zhu, Deng, Chen et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,137
    views
  • 90
    downloads
  • 0
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Yaying Zhu
  2. Chengxiang Deng
  3. Bo Chen
  4. Jia He
  5. Yanan Liu
  6. Shan Lei
  7. Weiju Lu
  8. Cheng Peng
  9. Zhao Shan
(2026)
Mid-zone hepatocytes trade proliferation for survival via Atf4-Chop axis in early acute liver injury
eLife 14:RP108646.
https://doi.org/10.7554/eLife.108646.4

Share this article

https://doi.org/10.7554/eLife.108646

Categories and tags

Research organism