Cited by 102Liver physiology and pathologySingle-cell and spatial transcriptomicsTissue Engineering and Regenerative Medicine

The DOI matches its Crossref record; 148 of the 148 references it cites were found in a registry or an index.

Cellular plasticity balances the metabolic and proliferation dynamics of a regenerating liver

Open access
Chembazhi, Ullas V., Bangru, Sushant, Hernaez, Mikel, Kalsotra, Auinash·2021·Genome Research

Published version of a preprint

Abstract

The adult liver has an exceptional ability to regenerate, but how it maintains its specialized functions during regeneration is unclear. Here, we used partial hepatectomy (PHx) in tandem with single-cell transcriptomics to track cellular transitions and heterogeneities of ∼22,000 liver cells through the initiation, progression, and termination phases of mouse liver regeneration. Our results uncovered that, following PHx, a subset of hepatocytes transiently reactivates an early-postnatal-like gene expression program to proliferate, while a distinct population of metabolically hyperactive cells appears to compensate for any temporary deficits in liver function. Cumulative EdU labeling and immunostaining of metabolic, portal, and central vein–specific markers revealed that hepatocyte proliferation after PHx initiates in the midlobular region before proceeding toward the periportal and pericentral areas. We further demonstrate that portal and central vein proximal hepatocytes retain their metabolically active state to preserve essential liver functions while midlobular cells proliferate nearby. Through combined analysis of gene regulatory networks and cell–cell interaction maps, we found that regenerating hepatocytes redeploy key developmental regulons, which are guided by extensive ligand-receptor-mediated signaling events between hepatocytes and nonparenchymal cells. Altogether, our study offers a detailed blueprint of the intercellular crosstalk and cellular reprogramming that balances the metabolic and proliferative requirements of a regenerating liver.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1101/gr.267013.120

Cites

148 of 148 references verified0 citation errors · 0 not checkable

Statistics

Figures

Figure 1. Single-cell analysis of resident hepatic cell populations from immature, adult, and regenerating mouse livers. (A) Time course plot showing the restoration of liver-to-body weight ratio after partial hepatectomy (PHx). The liver recovers its original mass within 7 d after PHx (n = 5). (B) Fluorescent imaging of hepatocyte proliferation measured by in vivo EdU incorporation in post-PHx and Sham livers. White arrows indicate proliferating hepatocytes (colabeled for HNF4A in green, incorporated EdU in red, and nuclei in blue). Images taken under 20× resolution are shown. (C) Overview schematic demonstrating workflow for isolation of mouse liver cells for single-cell RNA sequencing (scRNA-seq). Portal vein perfusion of collagenase containing buffer was used to isolate single liver cells from uninjured P14 pups and adults as well as mice at 24, 48, and 96 h after 2/3rd PHx (n = 1/time point). Single-cell library preparation was performed with whole-cell suspensions individually for each mouse using a 10x Genomics Chromium Single Cell 3′ Reagent Kit (V3 chemistry) after magnetic-activated cell sorting to remove dead cells. The inset details our PHx procedure, showing the position of two knots before excision of the respective liver lobes. (D) Combined UMAP projection of all 22,068 cells identified after QC cutoffs and batch correction. Cells are colored by the batch of origin, and the total number of cells identified from each batch is given in parentheses. (E) Identification of hepatocyte and nonparenchymal cell (NPC) subpopulations. Graph-based clustering in Seurat v3.1 followed by marker gene analysis revealed broad epithelial and nonepithelial cellular identities. Feature plots shown as insets show higher expression of expression of Hnf4a (a hepatocyte marker) and Vim (a nonepithelial marker) especially in populations identified as hepatocytes and NPCs, respectively. (F) Combined UMAP projection of all cells, colored by the annotated cell type. (G) PHATE projection of the ∼22,000 cells from different stages of liver development and regeneration. Cells are colored by annotated cell types from F.
Figure 1. Single-cell analysis of resident hepatic cell populations from immature, adult, and regenerating mouse livers. (A) Time course plot showing the restoration of liver-to-body weight ratio after partial hepatectomy (PHx). The liver recovers its original mass within 7 d after PHx (n = 5). (B) Fluorescent imaging of hepatocyte proliferation measured by in vivo EdU incorporation in post-PHx and Sham livers. White arrows indicate proliferating hepatocytes (colabeled for HNF4A in green, incorporated EdU in red, and nuclei in blue). Images taken under 20× resolution are shown. (C) Overview schematic demonstrating workflow for isolation of mouse liver cells for single-cell RNA sequencing (scRNA-seq). Portal vein perfusion of collagenase containing buffer was used to isolate single liver cells from uninjured P14 pups and adults as well as mice at 24, 48, and 96 h after 2/3rd PHx (n = 1/time point). Single-cell library preparation was performed with whole-cell suspensions individually for each mouse using a 10x Genomics Chromium Single Cell 3′ Reagent Kit (V3 chemistry) after magnetic-activated cell sorting to remove dead cells. The inset details our PHx procedure, showing the position of two knots before excision of the respective liver lobes. (D) Combined UMAP projection of all 22,068 cells identified after QC cutoffs and batch correction. Cells are colored by the batch of origin, and the total number of cells identified from each batch is given in parentheses. (E) Identification of hepatocyte and nonparenchymal cell (NPC) subpopulations. Graph-based clustering in Seurat v3.1 followed by marker gene analysis revealed broad epithelial and nonepithelial cellular identities. Feature plots shown as insets show higher expression of expression of Hnf4a (a hepatocyte marker) and Vim (a nonepithelial marker) especially in populations identified as hepatocytes and NPCs, respectively. (F) Combined UMAP projection of all cells, colored by the annotated cell type. (G) PHATE projection of the ∼22,000 cells from different stages of liver development and regeneration. Cells are colored by annotated cell types from F.CC BY-NC · PubMed Central
Figure 2. Specific hepatocyte population reversibly reprograms to an immature postnatal-like state during regeneration. (A) Ridge plots showing relative scoring on hepatocyte subpopulation using Seurat3.1 demonstrate extensive rewiring of metabolic genes during regeneration. Relative scores were computed based on the lists of genes for each pathway obtained from the Rat Genome Database (RGD). (B) Heat map showing relative scores of the top differentially regulated metabolic pathways. (C) Pseudotime plots demonstrating cellular trajectories during postnatal maturation (including P14 and adult hepatocytes), initiation, and progression (including adult, PHx24 and PHx48 hepatocytes), and termination and rematuration (including PHx48 and PHx96 hepatocytes). Single-cell trajectories were constructed and pseudotime values calculated using Monocle 2. Trajectories are colored by pseudotime (left) and sample identity (right). (D) Heat maps representing modules of genes that co-vary along the pseudotime during postnatal maturation, initiation–progression, and termination–rematuration phases. DAVID-based Gene Ontology (GO) analysis revealed that reversible reprogramming of developmentally regulated gene expression programs essentially reverts postnatal maturation, and this is followed by transitions that reinstate mature hepatic program. The top up-regulated and down-regulated GO terms are described below the respective heat maps. (E) Pseudotime plot indicating cellular trajectories of hepatocytes from all samples. Single-cell trajectories were constructed and pseudotime values were calculated using Monocle 2. Trajectories are colored by pseudotime. (F) Pseudotime plots showing distribution of each sample along combined cellular trajectories shown in C. The adult and P14 hepatocytes present distinct distributions along the trajectory; however, the distribution shifts toward P14 at PHx24–48 and back toward the adult at PHx96.
Figure 2. Specific hepatocyte population reversibly reprograms to an immature postnatal-like state during regeneration. (A) Ridge plots showing relative scoring on hepatocyte subpopulation using Seurat3.1 demonstrate extensive rewiring of metabolic genes during regeneration. Relative scores were computed based on the lists of genes for each pathway obtained from the Rat Genome Database (RGD). (B) Heat map showing relative scores of the top differentially regulated metabolic pathways. (C) Pseudotime plots demonstrating cellular trajectories during postnatal maturation (including P14 and adult hepatocytes), initiation, and progression (including adult, PHx24 and PHx48 hepatocytes), and termination and rematuration (including PHx48 and PHx96 hepatocytes). Single-cell trajectories were constructed and pseudotime values calculated using Monocle 2. Trajectories are colored by pseudotime (left) and sample identity (right). (D) Heat maps representing modules of genes that co-vary along the pseudotime during postnatal maturation, initiation–progression, and termination–rematuration phases. DAVID-based Gene Ontology (GO) analysis revealed that reversible reprogramming of developmentally regulated gene expression programs essentially reverts postnatal maturation, and this is followed by transitions that reinstate mature hepatic program. The top up-regulated and down-regulated GO terms are described below the respective heat maps. (E) Pseudotime plot indicating cellular trajectories of hepatocytes from all samples. Single-cell trajectories were constructed and pseudotime values were calculated using Monocle 2. Trajectories are colored by pseudotime. (F) Pseudotime plots showing distribution of each sample along combined cellular trajectories shown in C. The adult and P14 hepatocytes present distinct distributions along the trajectory; however, the distribution shifts toward P14 at PHx24–48 and back toward the adult at PHx96.CC BY-NC · PubMed Central
Figure 3. Bifurcation of hepatocyte trajectory during regeneration produces hepatocytes enriched with complementary functions in proliferation and metabolism. (A) Heat map showing bifurcating of gene expression programs executed along the pseudotime after branching. The top GO terms enriched in each class of genes are listed with their corresponding adjusted P-values. (B) Trajectory demonstrating the three distinct states of hepatocytes. The branch point under evaluation is shown in red. (C) Box plots demonstrating cell cycle phase scores calculated from Seurat v3.1 for hepatocytes belonging to different cell states. Q.S., T.S., P.S., and M.H.A. denote quiescent, transition, proliferating, and metabolically hyperactive states, respectively. P-values were derived from a parametric t-test (unpaired). (*) P ≤ 0.05, (****) P ≤ 0.0001. (D) “Proliferating” and “metabolically hyperactive” states uniquely up-regulate proliferation- or metabolism-related functions, respectively. Box plots showing relative scoring of indicated pathways in hepatocytes belonging to different states. P-values were derived from a parametric t-test (unpaired). (*) P ≤ 0.05, (**) P ≤ 0.01, (****) P ≤ 0.0001, (ns) P > 0.05. (E) Metabolically hyperactive state transiently up-regulates metabolism-related functions during regeneration. Box plots showing time point-based scoring of hepatocytes from the metabolically hyperactive state for the indicated pathways. P-values were derived from a parametric t-test (unpaired). (*) P ≤ 0.05, (**) P ≤ 0.01, (****) P ≤ 0.0001, (ns) P > 0.05.
Figure 3. Bifurcation of hepatocyte trajectory during regeneration produces hepatocytes enriched with complementary functions in proliferation and metabolism. (A) Heat map showing bifurcating of gene expression programs executed along the pseudotime after branching. The top GO terms enriched in each class of genes are listed with their corresponding adjusted P-values. (B) Trajectory demonstrating the three distinct states of hepatocytes. The branch point under evaluation is shown in red. (C) Box plots demonstrating cell cycle phase scores calculated from Seurat v3.1 for hepatocytes belonging to different cell states. Q.S., T.S., P.S., and M.H.A. denote quiescent, transition, proliferating, and metabolically hyperactive states, respectively. P-values were derived from a parametric t-test (unpaired). (*) P ≤ 0.05, (****) P ≤ 0.0001. (D) “Proliferating” and “metabolically hyperactive” states uniquely up-regulate proliferation- or metabolism-related functions, respectively. Box plots showing relative scoring of indicated pathways in hepatocytes belonging to different states. P-values were derived from a parametric t-test (unpaired). (*) P ≤ 0.05, (**) P ≤ 0.01, (****) P ≤ 0.0001, (ns) P > 0.05. (E) Metabolically hyperactive state transiently up-regulates metabolism-related functions during regeneration. Box plots showing time point-based scoring of hepatocytes from the metabolically hyperactive state for the indicated pathways. P-values were derived from a parametric t-test (unpaired). (*) P ≤ 0.05, (**) P ≤ 0.01, (****) P ≤ 0.0001, (ns) P > 0.05.CC BY-NC · PubMed Central
Figure 4. Gene regulatory networks are rewired to a postnatal-like state during regeneration. (A) UMAP projection of all hepatocytes based on the AUC scores for each regulon calculated with SCENIC. Cells are colored according to the sample of origin. (B) AUC score-based UMAP projection, grouped according to the sample of origin. Cells are colored according to sample of origin. Adult and PHx96 hepatocytes cluster together, whereas PHx24 and PHx48 hepatocytes cluster together with P14 hepatocytes. (C) Heat map depicting the activities of different regulons that show time point-dependent variations. (D) Violin plot showing distribution of AUC scores for RELA, E2F1, GABP1, and ETS2 regulons across hepatocytes from each time point demonstrating their high activity in PHx24, PHx48, and P14 hepatocytes. (E) Violin plots showing distribution of AUC scores for HNF4A, DBP, CBPA, and HES6 regulons across hepatocytes from each time point demonstrating their high activity in adult and PHx96 hepatocytes. (F) Violin plots showing distribution of AUC scores of representative regulons across hepatocytes showing their up-regulation in the proliferative state. (G) Violin plots showing distribution of AUC scores of representative regulons across hepatocytes showing their up-regulation in quiescent and metabolically active states. (H) Pseudotime plots of hepatocyte cellular trajectories colored by the AUC scores of representative regulons showing high activity in the proliferative state. (I) Pseudotime plots of hepatocyte cellular trajectories colored by the AUC scores of representative regulons showing high activity in quiescent and metabolically active states.
Figure 4. Gene regulatory networks are rewired to a postnatal-like state during regeneration. (A) UMAP projection of all hepatocytes based on the AUC scores for each regulon calculated with SCENIC. Cells are colored according to the sample of origin. (B) AUC score-based UMAP projection, grouped according to the sample of origin. Cells are colored according to sample of origin. Adult and PHx96 hepatocytes cluster together, whereas PHx24 and PHx48 hepatocytes cluster together with P14 hepatocytes. (C) Heat map depicting the activities of different regulons that show time point-dependent variations. (D) Violin plot showing distribution of AUC scores for RELA, E2F1, GABP1, and ETS2 regulons across hepatocytes from each time point demonstrating their high activity in PHx24, PHx48, and P14 hepatocytes. (E) Violin plots showing distribution of AUC scores for HNF4A, DBP, CBPA, and HES6 regulons across hepatocytes from each time point demonstrating their high activity in adult and PHx96 hepatocytes. (F) Violin plots showing distribution of AUC scores of representative regulons across hepatocytes showing their up-regulation in the proliferative state. (G) Violin plots showing distribution of AUC scores of representative regulons across hepatocytes showing their up-regulation in quiescent and metabolically active states. (H) Pseudotime plots of hepatocyte cellular trajectories colored by the AUC scores of representative regulons showing high activity in the proliferative state. (I) Pseudotime plots of hepatocyte cellular trajectories colored by the AUC scores of representative regulons showing high activity in quiescent and metabolically active states.CC BY-NC · PubMed Central
Figure 5. Proliferating and metabolically active hepatocytes are discretely localized within regenerating livers. (A) Schematic showing the experimental strategy for cumulative labeling of proliferating cells by continued feeding of EdU during liver regeneration. (B) Representative immunofluorescence images of proliferating hepatocytes measured by cumulative EdU incorporation after PHx or Sham surgeries (n = 4 mice/time point). Proliferating hepatocytes were colabeled for HNF4A (green) and EdU (red). White arrowheads point to the hepatocytes without any EdU incorporation at each time point, representing cells that had not proliferated. Nuclei were stained with To-Pro-3 dye (blue). (C) Percentages of proliferated (EdU+) and nonproliferated (EdU−) hepatocytes in regenerating livers at PHx36, 48, and 96 h. Data are mean ± SD; n = 4 mice/time point. (D) Representative immunofluorescence image demonstrating down-regulation of HNF4A protein specifically in proliferating nuclei (EdU+) of a regenerating liver. Liver sections were colabeled for HNF4A in green, EdU in red, and nuclei in blue. (E) Representative images showing an overlay of hepatic glycogen content by PAS staining (purple) and fluorescently detected EdU incorporation (yellow) in PHx48 livers. Nonproliferating hepatocytes exhibited minimal glycogen depletion (n = 4 mice/time point). (F) Representative immunofluorescence images showing spatial segregation of metabolic and proliferating hepatocytes through different phases of liver regeneration. Periportal (CHD1+) or pericentral (GLUL+) hepatocytes were colabeled (green) along with EdU (red). Nuclei were stained with To-Pro-3 dye (blue) (n = 4 mice/time point).
Figure 5. Proliferating and metabolically active hepatocytes are discretely localized within regenerating livers. (A) Schematic showing the experimental strategy for cumulative labeling of proliferating cells by continued feeding of EdU during liver regeneration. (B) Representative immunofluorescence images of proliferating hepatocytes measured by cumulative EdU incorporation after PHx or Sham surgeries (n = 4 mice/time point). Proliferating hepatocytes were colabeled for HNF4A (green) and EdU (red). White arrowheads point to the hepatocytes without any EdU incorporation at each time point, representing cells that had not proliferated. Nuclei were stained with To-Pro-3 dye (blue). (C) Percentages of proliferated (EdU+) and nonproliferated (EdU−) hepatocytes in regenerating livers at PHx36, 48, and 96 h. Data are mean ± SD; n = 4 mice/time point. (D) Representative immunofluorescence image demonstrating down-regulation of HNF4A protein specifically in proliferating nuclei (EdU+) of a regenerating liver. Liver sections were colabeled for HNF4A in green, EdU in red, and nuclei in blue. (E) Representative images showing an overlay of hepatic glycogen content by PAS staining (purple) and fluorescently detected EdU incorporation (yellow) in PHx48 livers. Nonproliferating hepatocytes exhibited minimal glycogen depletion (n = 4 mice/time point). (F) Representative immunofluorescence images showing spatial segregation of metabolic and proliferating hepatocytes through different phases of liver regeneration. Periportal (CHD1+) or pericentral (GLUL+) hepatocytes were colabeled (green) along with EdU (red). Nuclei were stained with To-Pro-3 dye (blue) (n = 4 mice/time point).CC BY-NC · PubMed Central
Figure 6. Dynamics of cell–cell interactions during liver regeneration. (A) Heat map showing expression of various ligand molecules and cellular receptors from different liver cell types (left) and from hepatocytes belonging to different cell states (right). (B) Network diagrams showing cell–cell interactions indicated by arrows (edges) pointing in the source-to-target direction. Thickness indicates the sum of weighted paths between populations, and the color of arrows corresponds to the source. Network diagrams for Adult, PHx24, PHx48, and PHx96 are shown. (C) Dot plot of representative inbound signals to hepatocytes at PHx48. Size of each dot indicates the weight of the corresponding ligand–receptor interaction and the color indicates negative log10 P-value of the source-to-target interaction. Empirical P-values were calculated and Benjamini–Hochberg correction was performed. (D) Dot plot of representative outbound signals from hepatocytes to various liver cells at PHx48.
Figure 6. Dynamics of cell–cell interactions during liver regeneration. (A) Heat map showing expression of various ligand molecules and cellular receptors from different liver cell types (left) and from hepatocytes belonging to different cell states (right). (B) Network diagrams showing cell–cell interactions indicated by arrows (edges) pointing in the source-to-target direction. Thickness indicates the sum of weighted paths between populations, and the color of arrows corresponds to the source. Network diagrams for Adult, PHx24, PHx48, and PHx96 are shown. (C) Dot plot of representative inbound signals to hepatocytes at PHx48. Size of each dot indicates the weight of the corresponding ligand–receptor interaction and the color indicates negative log10 P-value of the source-to-target interaction. Empirical P-values were calculated and Benjamini–Hochberg correction was performed. (D) Dot plot of representative outbound signals from hepatocytes to various liver cells at PHx48.CC BY-NC · PubMed Central
Figure 7. Division-of-labor model for liver regeneration. Following surgical resection (PHx), the remnant liver tissue regenerates quickly and restores its original size and function. Hepatocyte proliferation initiates in the midlobular region before proceeding toward the periportal and pericentral areas. We propose that a subset of residual hepatocytes in the midlobular area reversibly activate an early-postnatal-like gene program to enter a proliferative state. Simultaneously, a distinct population near the portal and central vein proximal regions up-regulates their metabolic gene program to offset any regeneration-induced deficits in liver function. These reversible cell state transitions are guided by distinct ligand-receptor-mediated signaling events between hepatocytes and nonparenchymal cells. Thus, the division of labor maximizes the benefit-cost ratio of regeneration for an organism, ensuring quick and robust replenishment of the hepatic parenchyma while sustaining adequate metabolic and detoxification activities. (NPCs) PA: portal artery; PV: portal vein; CV: central vein.
Figure 7. Division-of-labor model for liver regeneration. Following surgical resection (PHx), the remnant liver tissue regenerates quickly and restores its original size and function. Hepatocyte proliferation initiates in the midlobular region before proceeding toward the periportal and pericentral areas. We propose that a subset of residual hepatocytes in the midlobular area reversibly activate an early-postnatal-like gene program to enter a proliferative state. Simultaneously, a distinct population near the portal and central vein proximal regions up-regulates their metabolic gene program to offset any regeneration-induced deficits in liver function. These reversible cell state transitions are guided by distinct ligand-receptor-mediated signaling events between hepatocytes and nonparenchymal cells. Thus, the division of labor maximizes the benefit-cost ratio of regeneration for an organism, ensuring quick and robust replenishment of the hepatic parenchyma while sustaining adequate metabolic and detoxification activities. (NPCs) PA: portal artery; PV: portal vein; CV: central vein.CC BY-NC · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Research works printed but not linked: 1

  1. 10.1016/0270-9139(93)90234-E · DOI

Cited by

Sorted by citations

See all 102

Metadata

Authors
Chembazhi, U. V., Bangru, S., Hernaez, M., & Kalsotra, A.
Affiliations
University of Illinois Urbana-Champaign (US)
Corresponding author
Kalsotra, Auinash
Venue
Genome Research
Published
2021-03-01
Volume
31
Issue
4
Pages
576-591
ISSN
  • 1088-9051 (ISSN-L)
  • 1549-5469
OpenAlex
W3133820991
Cited by
102
References
147
Citations by year
2026: 13 · 2025: 19 · 2024: 22 · 2023: 26 · 2022: 9 · 2021: 10
Type
posted content
Language
en
Open access
Locations
  • Genome Research (publishedVersion)
  • bioRxiv (Cold Spring Harbor Laboratory) (acceptedVersion, cc-by-nc-nd)
  • Europe PMC (PubMed Central) (publishedVersion, cc-by-nc)
  • PubMed Central (publishedVersion, cc-by-nc)
  • PubMed (publishedVersion)
Funders
  • National Institutes of Health
  • NIH (R01HL126845, R01AA010154, T32-EB019944)
  • Muscular Dystrophy Association (MDA514335)
Topics
  • Health Sciences › Medicine › Hepatology › Liver physiology and pathology
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Single-cell and spatial transcriptomics
  • Health Sciences › Medicine › Surgery › Tissue Engineering and Regenerative Medicine
Keywords
liver regeneration, partial hepatectomy, hepatocyte proliferation, cellular plasticity, cell-to-cell communication, non-parenchymal cells, developmental gene regulatory networks, single-cell transcriptomics
MeSH terms
Cell Plasticity, Animals, Hepatectomy, Liver, Liver Regeneration, Hepatocytes, Cell Proliferation, Mice, Single-Cell Analysis, Transcriptome

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review