Differential gene expression in periportal and perivenous mouse hepatocytes

Transcriptomic analysis of AHR wildtype and Knock-out rat livers supports TCDD's role in AHR/ARNT-mediated circadian disruption and hepatotoxicity

Single, high doses of TCDD in rats are known to cause wasting, a progressive loss of 30 to 50% body weight and death within several weeks. To identify pathway perturbations at or near doses causing wasting, we examined differentially gene expression (DGE) and pathway enrichment in centrilobular (CL) and periportal (PP) regions of female rat livers following 6 dose levels of TCDD - 0, 3, 22, 100, 300, and 1000 ng/kg/day, 5 days/week for 4 weeks. At the higher doses, rats lost weight, had increased liver/body weight ratios and nearly complete cessation of liver cell proliferation, signs consistent with wasting. DGE curves were left shifted for the CL versus the PP regions. Canonical Phase I and Phase II genes were maximally increased at lower doses and remained elevated at all doses. At lower doses, ≤ 22 ng/kg/day in the CL and ≤ 100 ng/kg/day, upregulated genes showed transcription factor (TF) enrichment for AHR and ARNT. At the mid- and high-dose doses, there was a large number of downregulated genes and pathway enrichment for DEGs which showed downregulation of many cellular metabolism processes including those for steroids, fatty acid metabolism, pyruvate metabolism and citric acid cycle. There was significant TF enrichment of the hi-dose downregulated genes for RXR, ESR1, LXR, PPARalpha. At the highest dose, there was also pathway enrichment with upregulated genes for extracellular matrix organization, collagen formation, hemostasis and innate immune system. TCDD demonstrates most of its effects through binding the aryl hydrocarbon receptor (AHR) while the downregulation of metabolism genes at higher TCDD doses is known to be independent of AHR binding to DREs. Based on our results with DEG, we provide a hypothesis for wasting in which high doses of TCDD shift circadian processes away from the resting state, leading to greatly reduced synthesis of steroids and complex lipids needed for cell growth, and producing gene expression signals consistent with an epithelial-to-mesenchymal transition in hepatocytes.

The Role of Metabolic Reprogramming in the Tumor Immune Microenvironment: Mechanisms and Opportunities for Immunotherapy in Hepatocellular Carcinoma

As one of the emerging hallmarks of tumorigenesis and tumor progression, metabolic remodeling is common in the tumor microenvironment. Hepatocellular carcinoma (HCC) is the third leading cause of global tumor-related mortality, causing a series of metabolic alterations in response to nutrient availability and consumption to fulfill the demands of biosynthesis and carcinogenesis. Despite the efficacy of immunotherapy in treating HCC, the response rate remains unsatisfactory. Recently, research has focused on metabolic reprogramming and its effects on the immune state of the tumor microenvironment, and immune response rate. In this review, we delineate the metabolic reprogramming observed in HCC and its influence on the tumor immune microenvironment. We discuss strategies aimed at enhancing response rates and overcoming immune resistance through metabolic interventions, focusing on targeting glucose, lipid, or amino acid metabolism, as well as systemic regulation.

Potential effect of the non-neuronal cardiac cholinergic system on hepatic glucose and energy metabolism

The vagus nerve belongs to the parasympathetic nervous system, which is involved in the regulation of organs throughout the body. Since the discovery of the non-neuronal cardiac cholinergic system (NNCCS), several studies have provided evidence for the positive role of acetylcholine (ACh) released from cardiomyocytes against cardiovascular diseases, such as sympathetic hyperreactivity-induced cardiac remodeling and dysfunction as well as myocardial infarction. Non-neuronal ACh released from cardiomyocytes is believed to regulate key physiological functions of the heart, such as attenuating heart rate, offsetting hypertrophic signals, maintaining action potential propagation, and modulating cardiac energy metabolism through the muscarinic ACh receptor in an auto/paracrine manner. Moreover, the NNCCS may also affect peripheral remote organs (e.g., liver) through the vagus nerve. Remote ischemic preconditioning (RIPC) and NNCCS activate the central nervous system and afferent vagus nerve. RIPC affects hepatic glucose and energy metabolism through the central nervous system and vagus nerve. In this review, we discuss the mechanisms and potential factors responsible for NNCCS in glucose and energy metabolism in the liver.

Perturbed liver gene zonation in a mouse model of non-alcoholic steatohepatitis

Liver zonation characterizes the separation of metabolic pathways along the lobules and is required for optimal hepatic function. Wnt signaling is a master regulator of spatial liver zonation. A perivenous-periportal Wnt activity gradient orchestrates metabolic zonation by activating gene expression in perivenous hepatocytes, while suppressing gene expression in their periportal counterparts. However, the understanding as to the liver gene zonation and zonation regulators in diseases is limited. Non-alcoholic steatohepatitis (NASH) is a chronic liver disease characterized by fat accumulation, inflammation, and fibrosis. Here, we investigated the perturbation of liver gene zonation in a mouse NASH model by combining spatial transcriptomics, bulk RNAseq and in situ hybridization. Wnt-target genes represented a major subset of genes showing altered spatial expression in the NASH liver. The altered Wnt-target gene expression levels and zonation spatial patterns were in line with the up regulation of Wnt regulators and the augmentation of Wnt signaling. Particularly, we found that the Wnt activator Rspo3 expression was restricted to the perivenous zone in control liver but expanded to the periportal zone in NASH liver. AAV8-mediated RSPO3 overexpression in controls resulted in zonation changes, and further amplified the disturbed zonation of Wnt-target genes in NASH, similarly Rspo3 knockdown in Rspo3+/- mice resulted in zonation changes of Wnt-target genes in both chow and HFD mouse. Interestingly, there were no impacts on steatosis, inflammation, or fibrosis NASH pathology from RSPO3 overexpression nor Rspo3 knockdown. In summary, our study demonstrated the alteration of Wnt signaling in a mouse NASH model, leading to perturbed liver zonation.

Spatial hepatocyte plasticity of gluconeogenesis during the metabolic transitions between fed, fasted and starvation states

The liver acts as a master regulator of metabolic homeostasis in part by performing gluconeogenesis. This process is dysregulated in type 2 diabetes, leading to elevated hepatic glucose output. The parenchymal cells of the liver (hepatocytes) are heterogeneous, existing on an axis between the portal triad and the central vein, and perform distinct functions depending on location in the lobule. Here, using single cell analysis of hepatocytes across the liver lobule, we demonstrate that gluconeogenic gene expression ( Pck1 and G6pc ) is relatively low in the fed state and gradually increases first in the periportal hepatocytes during the initial fasting period. As the time of fasting progresses, pericentral hepatocyte gluconeogenic gene expression increases, and following entry into the starvation state, the pericentral hepatocytes show similar gluconeogenic gene expression to the periportal hepatocytes. Similarly, pyruvate-dependent gluconeogenic activity is approximately 10-fold higher in the periportal hepatocytes during the initial fasting state but only 1.5-fold higher in the starvation state. In parallel, starvation suppresses canonical beta-catenin signaling and modulates expression of pericentral and periportal glutamine synthetase and glutaminase, resulting in an enhanced pericentral glutamine-dependent gluconeogenesis. These findings demonstrate that hepatocyte gluconeogenic gene expression and gluconeogenic activity are highly spatially and temporally plastic across the liver lobule, underscoring the critical importance of using well-defined feeding and fasting conditions to define the basis of hepatic insulin resistance and glucose production.

Hepatic nutrient and hormone signaling to mTORC1 instructs the postnatal metabolic zonation of the liver

The metabolic functions of the liver are spatially organized in a phenomenon called zonation, linked to the differential exposure of portal and central hepatocytes to nutrient-rich blood. The mTORC1 signaling pathway controls cellular metabolism in response to nutrients and insulin fluctuations. Here we show that simultaneous genetic activation of nutrient and hormone signaling to mTORC1 in hepatocytes results in impaired establishment of postnatal metabolic and zonal identity of hepatocytes. Mutant hepatocytes fail to upregulate postnatally the expression of Frizzled receptors 1 and 8, and show reduced Wnt/β-catenin activation. This defect, alongside diminished paracrine Wnt2 ligand expression by endothelial cells, underlies impaired postnatal maturation. Impaired zonation is recapitulated in a model of constant supply of nutrients by parenteral nutrition to piglets. Our work shows the role of hepatocyte sensing of fluctuations in nutrients and hormones for triggering a latent metabolic zonation program.

Region-specific cellular and molecular basis of liver regeneration after acute pericentral injury

Liver injuries often occur in a zonated manner. However, detailed regenerative responses to such zonal injuries at cellular and molecular levels remain largely elusive. By using a fate-mapping strain, Cyp2e1-DreER, to elucidate liver regeneration after acute pericentral injury, we found that pericentral regeneration is primarily compensated by the expansion of remaining pericentral hepatocytes, and secondarily by expansion of periportal hepatocytes. Employing single-cell RNA sequencing, spatial transcriptomics, immunostaining, and in vivo functional assays, we demonstrated that the upregulated expression of the mTOR/4E-BP1 axis and lactate dehydrogenase A in hepatocytes contributes to pericentral regeneration, while activation of transforming growth factor β (TGF-β1) signaling in the damaged area mediates fibrotic responses and inhibits hepatocyte proliferation. Inhibiting the pericentral accumulation of monocytes and monocyte-derived macrophages through an Arg-Gly-Asp (RGD) peptide-based strategy attenuates these cell-derived TGF-β1 signalings, thus improving pericentral regeneration. Our study provides integrated and high-resolution spatiotemporal insights into the cellular and molecular basis of pericentral regeneration.

A single-nucleus transcriptomic atlas of primate liver aging uncovers the pro-senescence role of SREBP2 in hepatocytes

Aging increases the risk of liver diseases and systemic susceptibility to aging-related diseases. However, cell type-specific changes and the underlying mechanism of liver aging in higher vertebrates remain incompletely characterized. Here, we constructed the first single-nucleus transcriptomic landscape of primate liver aging, in which we resolved cell type-specific gene expression fluctuation in hepatocytes across three liver zonations and detected aberrant cell-cell interactions between hepatocytes and niche cells. Upon in-depth dissection of this rich dataset, we identified impaired lipid metabolism and upregulation of chronic inflammation-related genes prominently associated with declined liver functions during aging. In particular, hyperactivated sterol regulatory element-binding protein (SREBP) signaling was a hallmark of the aged liver, and consequently, forced activation of SREBP2 in human primary hepatocytes recapitulated in vivo aging phenotypes, manifesting as impaired detoxification and accelerated cellular senescence. This study expands our knowledge of primate liver aging and informs the development of diagnostics and therapeutic interventions for liver aging and associated diseases.

Biomarker discovery in acetaminophen hepatotoxicity: leveraging single-cell transcriptomics and mechanistic insight

INTRODUCTION

Acetaminophen (APAP) overdose is the leading cause of drug-induced liver injury and can cause a rapid progression to acute liver failure (ALF). Therefore, the identification of prognostic biomarkers to determine which patients will require a liver transplant is critical for APAP-induced ALF.

AREAS COVERED

We begin by relating the mechanistic investigations in mouse models of APAP hepatotoxicity to the human APAP overdose pathophysiology. We draw insights from the established sequence of molecular events in mice to understand the progression of events in the APAP overdose patient. Through this mechanistic understanding, several new biomarkers, such as CXCL14, have recently been evaluated. We also explore how single-cell RNA sequencing, spatial transcriptomics, and other omics approaches have been leveraged for identifying novel biomarkers and how these approaches will continue to push the field of biomarker discovery forward.

EXPERT OPINION

Recent investigations have elucidated several new biomarkers or combination of markers such as CXCL14, a regenerative miRNA signature, a cell death miRNA signature, hepcidin, LDH, CPS1, and FABP1. While these biomarkers are promising, they all require further validation. Larger cohort studies analyzing these new biomarkers in the same patient samples, while adding these candidate biomarkers to prognostic models will further support their clinical utility.

Wnt-β-catenin in hepatobiliary homeostasis, injury, and repair

Wnt-β-catenin signaling has emerged as an important regulatory pathway in the liver, playing key roles in zonation and mediating contextual hepatobiliary repair after injuries. In this review, we will address the major advances in understanding the role of Wnt signaling in hepatic zonation, regeneration, and cholestasis-induced injury. We will also touch on some important unanswered questions and discuss the relevance of modulating the pathway to provide therapies for complex liver pathologies that remain a continued unmet clinical need.

Spatiotemporally quantitative in vivo imaging of mitochondrial fatty acid β-oxidation at cellular-level resolution in mice

Mitochondrial fatty acid β-oxidation (FAO) plays a key role in energy homeostasis. Several FAO evaluation methods are currently available, but they are not necessarily suitable for capturing the dynamics of FAO in vivo at a cellular-level spatial resolution and seconds-level time resolution. FAOBlue is a coumarin-based probe that undergoes β-oxidation to produce a fluorescent substrate, 7-hydroxycoumarin-3-(N-(2-hydroxyethyl))-carboxamide (7-HC). After confirming that 7-HC could be specifically detected using multiphoton microscopy at excitation/emission wavelength = 820/415-485 nm, wild-type C57BL/6 mice were randomly divided into control, pemafibrate, fasting (24 or 72 h), and etomoxir groups. These mice received a single intravenous injection of FAOBlue. FAO activities in the liver of these mice were visualized using multiphoton microscopy at 4.2 s/frame. These approaches could visualize the difference in FAO activities between periportal and pericentral hepatocytes in the control, pemafibrate, and fasting groups. FAO velocity, which was expressed by the maximum slope of the fluorescence intensity curve, was accelerated in the pemafibrate and 72-h fasting groups both in the periportal and the pericentral hepatocytes in comparison with the control group. Our approach revealed differences in the FAO activation mode by the two stimuli, i.e., pemafibrate and fasting, with pemafibrate accelerating the time of first detection of FAO-derived fluorescence. No increase in the fluorescence was observed in etomoxir-pretreated mice, confirming that FAOBlue specifically detected FAO in vivo. Thus, FAOBlue is useful for visualizing in vivo liver FAO dynamics at the single-cell-level spatial resolution and seconds-level time resolution.NEW & NOTEWORTHY Fatty acid β-oxidation (FAO) plays a key role in energy homeostasis. Here, the authors established a strategy for visualizing FAO activity in vivo at the cellular-level spatial resolution and seconds-level time resolution in mice. Quantitative analysis revealed spatiotemporal heterogeneity in hepatic FAO dynamics. Our method is widely applicable because it is simple and uses a multiphoton microscope to observe the FAOBlue-injected mice.

Desorption electrospray ionization mass spectrometry imaging (DESI-MSI) in disease diagnosis: an overview

Tissue analysis, which is essential to histology and is considered the benchmark for the diagnosis and prognosis of many illnesses, including cancer, is significant. During surgery, the surgical margin of the tumor is assessed using the labor-intensive, challenging, and commonly subjective technique known as frozen section histopathology. In the biopsy section, large numbers of molecules can now be visualized at once (ion images) following recent developments in [MSI] mass spectrometry imaging under atmospheric conditions. This is vastly superior to and different from the single optical tissue image processing used in traditional histopathology. This review article will focus on the advancement of desorption electrospray ionization mass spectrometry imaging [DESI-MSI] technique, which is label-free and requires little to no sample preparation. Since the proportion of molecular species in normal and abnormal tissues is different, DESI-MSI can capture ion images of the distributions of lipids and metabolites on biopsy sections, which can provide rich diagnostic information. This is not a systematic review but a summary of well-known, cutting-edge and recent DESI-MSI applications in cancer research between 2018 and 2023.

Liver in infections: a single-cell and spatial transcriptomics perspective

The liver is an immune organ that plays a vital role in the detection, capture, and clearance of pathogens and foreign antigens that invade the human body. During acute and chronic infections, the liver transforms from a tolerant to an active immune state. The defence mechanism of the liver mainly depends on a complicated network of intrahepatic and translocated immune cells and non-immune cells. Therefore, a comprehensive liver cell atlas in both healthy and diseased states is needed for new therapeutic target development and disease intervention improvement. With the development of high-throughput single-cell technology, we can now decipher heterogeneity, differentiation, and intercellular communication at the single-cell level in sophisticated organs and complicated diseases. In this concise review, we aimed to summarise the advancement of emerging high-throughput single-cell technologies and re-define our understanding of liver function towards infections, including hepatitis B virus, hepatitis C virus, Plasmodium, schistosomiasis, endotoxemia, and corona virus disease 2019 (COVID-19). We also unravel previously unknown pathogenic pathways and disease mechanisms for the development of new therapeutic targets. As high-throughput single-cell technologies mature, their integration into spatial transcriptomics, multiomics, and clinical data analysis will aid in patient stratification and in developing effective treatment plans for patients with or without liver injury due to infectious diseases.

Single-cell and spatiotemporal transcriptomic analyses reveal the effects of microorganisms on immunity and metabolism in the mouse liver

The gut-liver axis is a complex bidirectional communication pathway between the intestine and the liver in which microorganisms and their metabolites flow from the intestine through the portal vein to the liver and influence liver function. In a sterile environment, the phenotype or function of the liver is altered, but few studies have investigated the specific cellular and molecular effects of microorganisms on the liver. To this end, we constructed single-cell and spatial transcriptomic (ST) profiles of germ-free (GF) and specific-pathogen-free (SPF) mouse livers. Single-cell RNA sequencing (scRNA-seq) and single-nucleus RNA sequencing (snRNA-seq) revealed that the ratio of most immune cells was altered in the liver of GF mice; in particular, natural killer T (NKT) cells, IgA plasma cells (IgAs) and Kupffer cells (KCs) were significantly reduced in GF mice. Spatial enhanced resolution omics sequencing (Stereo-seq) confirmed that microorganisms mediated the accumulation of Kupffer cells in the periportal zone. Unexpectedly, IgA plasma cells were more numerous and concentrated in the periportal vein in liver sections from SPF mice but less numerous and scattered in GF mice. ST technology also enables the precise zonation of liver lobules into eight layers and three patterns based on the gene expression level in each layer, allowing us to further investigate the effects of microbes on gene zonation patterns and functions. Furthermore, untargeted metabolism experiments of the liver revealed that the propionic acid levels were significantly lower in GF mice, and this reduction may be related to the control of genes involved in bile acid and fatty acid metabolism. In conclusion, the combination of sc/snRNA-seq, Stereo-seq, and untargeted metabolomics revealed immune system defects as well as altered bile acid and lipid metabolic processes at the single-cell and spatial levels in the livers of GF mice. This study will be of great value for understanding host-microbiota interactions.

eSPRESSO: topological clustering of single-cell transcriptomics data to reveal informative genes for spatio-temporal architectures of cells

BACKGROUND

Bioinformatics capability to analyze spatio-temporal dynamics of gene expression is essential in understanding animal development. Animal cells are spatially organized as functional tissues where cellular gene expression data contain information that governs morphogenesis during the developmental process. Although several computational tissue reconstruction methods using transcriptomics data have been proposed, those methods have been ineffective in arranging cells in their correct positions in tissues or organs unless spatial information is explicitly provided.

RESULTS

This study demonstrates stochastic self-organizing map clustering with Markov chain Monte Carlo calculations for optimizing informative genes effectively reconstruct any spatio-temporal topology of cells from their transcriptome profiles with only a coarse topological guideline. The method, eSPRESSO (enhanced SPatial REconstruction by Stochastic Self-Organizing Map), provides a powerful in silico spatio-temporal tissue reconstruction capability, as confirmed by using human embryonic heart and mouse embryo, brain, embryonic heart, and liver lobule with generally high reproducibility (average max. accuracy = 92.0%), while revealing topologically informative genes, or spatial discriminator genes. Furthermore, eSPRESSO was used for temporal analysis of human pancreatic organoids to infer rational developmental trajectories with several candidate 'temporal' discriminator genes responsible for various cell type differentiations.

CONCLUSIONS

eSPRESSO provides a novel strategy for analyzing mechanisms underlying the spatio-temporal formation of cellular organizations.

Adenoviral Transfer of Human Aquaporin-8 Gene to Mouse Liver Improves Ammonia-Derived Ureagenesis

We previously reported that, in cultured hepatocytes, mitochondrial aquaporin-8 (AQP8) channels facilitate the conversion of ammonia to urea and that the expression of human AQP8 (hAQP8) enhances ammonia-derived ureagenesis. In this study, we evaluated whether hepatic gene transfer of hAQP8 improves detoxification of ammonia to urea in normal mice as well as in mice with impaired hepatocyte ammonia metabolism. A recombinant adenoviral (Ad) vector encoding hAQP8, AdhAQP8, or a control Ad vector was administered via retrograde infusion into the bile duct of the mice. Hepatocyte mitochondrial expression of hAQP8 was confirmed using confocal immunofluorescence and immunoblotting. The normal hAQP8-transduced mice showed decreased plasma ammonia and increased liver urea. Enhanced ureagenesis was confirmed via the NMR studies assessing the synthesis of 15N-labeled urea from 15N-labeled ammonia. In separate experiments, we made use of the model hepatotoxic agent, thioacetamide, to induce defective hepatic metabolism of ammonia in mice. The adenovirus-mediated mitochondrial expression of hAQP8 was able to restore normal ammonemia and ureagenesis in the liver of the mice. Our data suggest that hAQP8 gene transfer to mouse liver improves detoxification of ammonia to urea. This finding could help better understand and treat disorders with defective hepatic ammonia metabolism.

Characterisation of hepatic lipid signature distributed across the liver zonation using mass spectrometry imaging

Background & Aims

Lipid metabolism plays an important role in liver pathophysiology. The liver lobule asymmetrically distributes oxygen and nutrition, resulting in heterogeneous metabolic functions. Periportal and pericentral hepatocytes have different metabolic functions, which lead to generating liver zonation. We developed spatial metabolic imaging using desorption electrospray ionisation mass spectrometry to investigate lipid distribution across liver zonation with high reproducibility and accuracy.

Methods

Fresh frozen livers from healthy mice with control diet were analysed using desorption electrospray ionisation mass spectrometry imaging. Imaging was performed at 50 μm × 50 μm pixel size. Regions of interest (ROIs) were manually created by co-registering with histological data to determine the spatial hepatic lipids across liver zonation. The ROIs were confirmed by double immunofluorescence. The mass list of specific ROIs was automatically created, and univariate and multivariate statistical analysis were performed to identify statistically significant lipids across liver zonation.

Results

A wide range of lipid species was identified, including fatty acids, phospholipids, triacylglycerols, diacylglycerols, ceramides, and sphingolipids. We characterised hepatic lipid signatures in three different liver zones (periportal zone, midzone, and pericentral zone) and validated the reproducibility of our method for measuring a wide range of lipids. Fatty acids were predominantly detected in the periportal region, whereas phospholipids were distributed in both the periportal and pericentral zones. Interestingly, phosphatidylinositols, PI(36:2), PI(36:3), PI(36:4), PI(38:5), and PI(40:6) were located predominantly in the midzone (zone 2). Triacylglycerols and diacylglycerols were detected mainly in the pericentral region. De novo triacylglycerol biosynthesis appeared to be the most influenced pathway across the three zones.

Conclusions

The ability to accurately assess zone-specific hepatic lipid distribution in the liver could lead to a better understanding of lipid metabolism during the progression of liver disease.

Impact and Implications

Zone-specific hepatic lipid metabolism could play an important role in lipid homoeostasis during disease progression. Herein, we defined the zone-specific references of hepatic lipid species in the three liver zones using molecular imaging. The de novo triacylglycerol biosynthesis was highlighted as the most influenced pathway across the three zones.

Hepatic stellate cells maintain liver homeostasis through paracrine neurotrophin-3 signaling that induces hepatocyte proliferation

Organ size is maintained by the controlled proliferation of distinct cell populations. In the mouse liver, hepatocytes in the midlobular zone that are positive for cyclin D1 (CCND1) repopulate the parenchyma at a constant rate to preserve liver mass. Here, we investigated how hepatocyte proliferation is supported by hepatic stellate cells (HSCs), pericytes that are in close proximity to hepatocytes. We used T cells to ablate nearly all HSCs in the murine liver, enabling the unbiased characterization of HSC functions. In the normal liver, complete loss of HSCs persisted for up to 10 weeks and caused a gradual reduction in liver mass and in the number of CCND1+ hepatocytes. We identified neurotrophin-3 (Ntf-3) as an HSC-produced factor that induced the proliferation of midlobular hepatocytes through the activation of tropomyosin receptor kinase B (TrkB). Treating HSC-depleted mice with Ntf-3 restored CCND1+ hepatocytes in the midlobular region and increased liver mass. These findings establish that HSCs form the mitogenic niche for midlobular hepatocytes and identify Ntf-3 as a hepatocyte growth factor.

Spatiotemporal Metabolic Liver Zonation and Consequences on Pathophysiology

Hepatocytes are the main workers in the hepatic factory, managing metabolism of nutrients and xenobiotics, production and recycling of proteins, and glucose and lipid homeostasis. Division of labor between hepatocytes is critical to coordinate complex complementary or opposing multistep processes, similar to distributed tasks at an assembly line. This so-called metabolic zonation has both spatial and temporal components. Spatial distribution of metabolic function in hepatocytes of different lobular zones is necessary to perform complex sequential multistep metabolic processes and to assign metabolic tasks to the right environment. Moreover, temporal control of metabolic processes is critical to align required metabolic processes to the feeding and fasting cycles. Disruption of this complex spatiotemporal hepatic organization impairs key metabolic processes with both local and systemic consequences. Many metabolic diseases, such as nonalcoholic steatohepatitis and diabetes, are associated with impaired metabolic liver zonation. Recent technological advances shed new light on the spatiotemporal gene expression networks controlling liver function and how their deregulation may be involved in a large variety of diseases. We summarize the current knowledge about spatiotemporal metabolic liver zonation and consequences on liver pathobiology.

Characterization of the humanized FRG mouse model and development of an AAV-LK03 variant with improved liver lobular biodistribution

Recent clinical successes have intensified interest in using adeno-associated virus (AAV) vectors for therapeutic gene delivery. The liver is a key clinical target, given its critical physiological functions and involvement in a wide range of genetic diseases. In the present study, we first investigated the validity of a liver xenograft mouse model repopulated with primary hepatocytes using single-nucleus RNA sequencing (sn-RNA-seq) by studying the transcriptomic profile of human hepatocytes pre- and post-engraftment. Complementary immunofluorescence analyses performed in highly engrafted animals confirmed that the human hepatocytes organize and present appropriate patterns of zone-dependent enzyme expression in this model. Next, we tested a set of rationally designed HSPG de-targeted AAV-LK03 variants for relative transduction performance in human hepatocytes. We used immunofluorescence, next-generation sequencing, and single-nucleus transcriptomics data from highly engrafted FRG mice to demonstrate that the optimally HSPG de-targeted AAV-LK03 displayed a significantly improved lobular transduction profile in this model.

Spatial transcriptome profiling of normal human liver

The comprehensive study of the spatial-cellular anatomy of the human liver is critical to addressing the cellular origins of liver disease. Here we conducted spatial transcriptomics on normal human liver tissue sections, providing detailed information of liver zonation at the transcriptional level. We present 6581 high-quality spots from normal livers of two human donors. In this dataset, cells were mainly hepatocytes, and we classified them into four sub-groups. Collectively, these data provide a reliable reference for studies on spatial heterogeneity of liver lobules.

Measurement(s)	mRNA Expression
Technology Type(s)	spatial transcriptomics
Sample Characteristic - Organism	Homo sapiens

Liver zonation, revisited

The concept of hepatocyte functional zonation is well established, with differences in metabolism and xenobiotic processing determined by multiple factors including oxygen and nutrient levels across the hepatic lobule. However, recent advances in single-cell genomics technologies, including single-cell and nuclei RNA sequencing, and the rapidly evolving fields of spatial transcriptomic and proteomic profiling have greatly increased our understanding of liver zonation. Here we discuss how these transformative experimental strategies are being leveraged to dissect liver zonation at unprecedented resolution and how this new information should facilitate the emergence of novel precision medicine-based therapies for patients with liver disease.

The RSPO-LGR4/5-ZNRF3/RNF43 module in liver homeostasis, regeneration, and disease

WNT/β-catenin signaling plays pivotal roles during liver development, homeostasis, and regeneration. Likewise, its deregulation disturbs metabolic liver zonation and is responsible for the development of a large number of hepatic tumors. Liver fibrosis, which has become a major health burden for society and a hallmark of NASH, can also be promoted by WNT/β-catenin signaling. Upstream regulatory mechanisms controlling hepatic WNT/β-catenin activity may constitute targets for the development of novel therapies addressing these life-threatening conditions. The R-spondin (RSPO)-leucine-rich repeat-containing G protein-coupled receptor (LGR) 4/5-zinc and ring finger (ZNRF) 3/ring finger 43 (RNF43) module is fine-tuning WNT/β-catenin signaling in several tissues and is essential for hepatic WNT/β-catenin activity. In this review article, we recapitulate the role of the RSPO-LGR4/5-ZNRF3/RNF43 module during liver development, homeostasis, metabolic zonation, regeneration, and disease. We further discuss the controversy around LGR5 as a liver stem cell marker.

Pregnane X receptor promotes liver enlargement in mice through the spatial induction of hepatocyte hypertrophy and proliferation

Nuclear receptor pregnane X receptor (PXR) can induce significant liver enlargement through hepatocyte hypertrophy and proliferation. A previous report showed that during the process of PXR-induced liver enlargement, hepatocyte hypertrophy occurs around the central vein (CV) area while hepatocyte proliferation occurs around the portal vein (PV) area. However, the features of this spatial change remain unclear. Therefore, this study aims to explore the features of the spatial changes in hepatocytes in PXR-induced liver enlargement. PXR-induced spatial changes in hepatocyte hypertrophy and proliferation were confirmed in C57BL/6 mice. The liver was perfused with digitonin to destroy the hepatocytes around the CV or PV areas, and then the regional expression of proteins related to hepatocyte hypertrophy and proliferation was further measured. The results showed that the expression of PXR downstream proteins, such as cytochrome P450 (CYP) 3A11, CYP2B10, P-glycoprotein (P-gp) and organ anion transporting polypeptide 2 (OATP2) was upregulated around the CV area, while the expression of proliferation-related proteins such as cyclin B1 (CCNB1), cyclin D1 (CCND1) and serine/threonine NIMA-related kinase 2 (NEK2) was upregulated around the PV area. At the same time, the expression of cyclin-dependent kinase inhibitors such as retinoblastoma-like protein 2 (RBL2), cyclin-dependent kinase inhibitor 1B (CDKN1B) and CDKN1A was downregulated around the PV area. This study demonstrated that the spatial change in PXR-induced hepatocyte hypertrophy and proliferation is associated with the regional expression of PXR downstream targets and proliferation-related proteins and the regional distribution of triglycerides (TGs). These findings provide new insight into the understanding of PXR-induced hepatomegaly.

Transcriptional switch of hepatocytes initiates macrophage recruitment and T-cell suppression in endotoxemia

BACKGROUND & AIMS

The liver plays crucial roles in the regulation of immune defense during acute systemic infections. However, the roles of liver cellular clusters and intercellular communication in the progression of endotoxemia have not been well-characterized.

METHODS

Single-cell RNA sequencing analysis was performed, and the transcriptomes of 19,795 single liver cells from healthy and endotoxic mice were profiled. The spatial and temporal changes in hepatocytes and non-parenchymal cell types were validated by multiplex immunofluorescence staining, bulk transcriptomic sequencing, or flow cytometry. Furthermore, we used an adeno-associated virus delivery system to confirm the major mechanisms mediating myeloid cell infiltration and T-cell suppression in septic murine liver.

RESULTS

We identified a proinflammatory hepatocyte (PIH) subpopulation that developed primarily from periportal hepatocytes and to a lesser extent from pericentral hepatocytes and played key immunoregulatory roles in endotoxemia. Multicellular cluster modeling of ligand-receptor interactions revealed that PIHs play a crucial role in the recruitment of macrophages via the CCL2-CCR2 interaction. Recruited macrophages (RMs) released cytokines (e.g., IL6, TNFα, and IL17) to induce the expression of inhibitory ligands, such as PD-L1, on hepatocytes. Subsequently, RM-stimulated hepatocytes led to the suppression of CD4⁺ and memory T-cell subsets partly via the PD-1/PD-L1 interaction in endotoxemia. Furthermore, sinusoidal endothelial cells expressed the highest levels of proapoptotic and inflammatory genes around the periportal zone. This pattern of gene expression facilitated increases in the number of fenestrations and infiltration of immune cells in the periportal zone.

CONCLUSIONS

Our study elucidates unanticipated aspects of the cellular and molecular effects of endotoxemia on liver cells at the single-cell level and provides a conceptual framework for the development of novel therapeutic approaches for acute infection.

LAY SUMMARY

The liver plays a crucial role in the regulation of immune defense during acute systemic infections. We identified a proinflammatory hepatocyte subpopulation and demonstrated that the interactions of this subpopulation with recruited macrophages are pivotal in the immune response during endotoxemia. These novel findings provide a conceptual framework for the discovery of rational therapeutic targets in acute infection.

Liver regeneration: Cellular origin and molecular mechanisms

The liver is known as an organ with high proliferation potential. Clarifying the cellular origin and deepening the understanding of liver regeneration mechanisms will help provide new directions for the treatment of liver disease. With the development and application of lineage tracing technology, the specific distribution and dynamic changes of hepatocyte subpopulations in homeostasis and liver injury have been illustrated. Self-replication of hepatocytes is responsible for the maintenance of liver function and mass under homeostasis. The compensatory proliferation of remaining hepatocytes is the main mechanism of liver regeneration following acute and chronic liver injury. Transdifferentiation between hepatocytes and cholangiocytes has been recognized upon severe chronic liver injury. Wnt/β-catenin, Hippo/YAP and Notch signalling play essential roles in the maintenance of homeostatic liver and hepatocyte-to-cholangiocyte conversion under liver injury. In this review, we summarized the recent studies on cell origin of newly generated hepatocytes and the underlying mechanisms of liver regeneration in homeostasis and liver injury.

From Seeing to Simulating: A Survey of Imaging Techniques and Spatially-Resolved Data for Developing Multiscale Computational Models of Liver Regeneration

Liver regeneration, which leads to the re-establishment of organ mass, follows a specifically organized set of biological processes acting on various time and length scales. Computational models of liver regeneration largely focused on incorporating molecular and signaling detail have been developed by multiple research groups in the recent years. These modeling efforts have supported a synthesis of disparate experimental results at the molecular scale. Incorporation of tissue and organ scale data using noninvasive imaging methods can extend these computational models towards a comprehensive accounting of multiscale dynamics of liver regeneration. For instance, microscopy-based imaging methods provide detailed histological information at the tissue and cellular scales. Noninvasive imaging methods such as ultrasound, computed tomography and magnetic resonance imaging provide morphological and physiological features including volumetric measures over time. In this review, we discuss multiple imaging modalities capable of informing computational models of liver regeneration at the organ-, tissue- and cellular level. Additionally, we discuss available software and algorithms, which aid in the analysis and integration of imaging data into computational models. Such models can be generated or tuned for an individual patient with liver disease. Progress towards integrated multiscale models of liver regeneration can aid in prognostic tool development for treating liver disease.

Amino acid metabolism, transport and signalling in the liver revisited

The liver controls the systemic exposure of amino acids entering via the gastro-intestinal tract. For most amino acids except branched chain amino acids, hepatic uptake is very efficient. This implies that the liver orchestrates amino acid metabolism and also controls systemic amino acid exposure. Although many amino acid transporters have been identified, cloned and investigated with respect to substrate specificity, transport mechanism, and zonal distribution, which of these players are involved in hepatocellular amino acid transport remains unclear. Here, we aim to provide a review of current insight into the molecular machinery of hepatic amino acid transport. Furthermore, we place this information in a comprehensive overview of amino acid transport, signalling and metabolism.

Temporal analyses of postnatal liver development and maturation by single-cell transcriptomics

The postnatal development and maturation of the liver, the major metabolic organ, are inadequately understood. We have analyzed 52,834 single-cell transcriptomes and identified 31 cell types or states in mouse livers at postnatal days 1, 3, 7, 21, and 56. We observe unexpectedly high levels of hepatocyte heterogeneity in the developing liver and the progressive construction of the zonated metabolic functions from pericentral to periportal hepatocytes, which is orchestrated with the development of sinusoid endothelial, stellate, and Kupffer cells. Trajectory and gene regulatory analyses capture 36 transcription factors, including a circadian regulator, Bhlhe40, in programming liver development. Remarkably, we identified a special group of macrophages enriched at day 7 with a hybrid phenotype of macrophages and endothelial cells, which may regulate sinusoidal construction and Treg-cell function. This study provides a comprehensive atlas that covers all hepatic cell types and is instrumental for further dissection of liver development, metabolism, and disease.

Hepatocyte glutathione S-transferase mu 2 prevents non-alcoholic steatohepatitis by suppressing ASK1 signaling

BACKGROUND & AIMS

Non-alcoholic fatty liver disease (NAFLD) has become the most common chronic liver disease worldwide. The advanced stage of NAFLD, non-alcoholic steatohepatitis (NASH), has been recognized as a leading cause of end-stage liver injury for which there are no FDA-approved therapeutic options. Glutathione S-transferase Mu 2 (GSTM2) is a phase II detoxification enzyme. However, the roles of GSTM2 in NASH have not been elucidated.

METHODS

Multiple RNA-seq analyses were used to identify hepatic GSTM2 expression in NASH. In vitro and in vivo gain- or loss-of-function approaches were used to investigate the role and molecular mechanism of GSTM2 in NASH.

RESULTS

We identified GSTM2 as a sensitive responder and effective suppressor of NASH progression. GSTM2 was significantly downregulated during NASH progression. Hepatocyte GSTM2 deficiency markedly aggravated insulin resistance, hepatic steatosis, inflammation and fibrosis induced by a high-fat diet and a high-fat/high-cholesterol diet. Mechanistically, GSTM2 sustained MAPK pathway signaling by directly interacting with apoptosis signal-regulating kinase 1 (ASK1). GSTM2 directly bound to the N-terminal region of ASK1 and inhibited ASK1 N-terminal dimerization to subsequently repress ASK1 phosphorylation and the activation of its downstream JNK/p38 signaling pathway under conditions of metabolic dysfunction.

CONCLUSIONS

These data demonstrated that hepatocyte GSTM2 is an endogenous suppressor that protects against NASH progression by blocking ASK1 N-terminal dimerization and phosphorylation. Activating GSTM2 holds promise as a therapeutic strategy for NASH.

CLINICAL TRIAL NUMBER

IIT-2021-277.

LAY SUMMARY

New therapeutic strategies for non-alcoholic steatohepatitis are urgently needed. We identified that the protein GSTM2 exerts a protective effect in response to metabolic stress. Therapies that aim to increase the activity of GSTM2 could hold promise for the treatment of non-alcoholic steatohepatitis.

Role and Regulation of Wnt/β-Catenin in Hepatic Perivenous Zonation and Physiological Homeostasis

Metabolic heterogeneity or functional zonation is a key characteristic of the liver that allows different metabolic pathways to be spatially regulated within the hepatic system and together contribute to whole body homeostasis. These metabolic pathways are segregated along the portocentral axis of the liver lobule into three hepatic zones: periportal, intermediate or midzonal, and perivenous. The liver performs complementary or opposing metabolic functions within different hepatic zones while synergistic functions are regulated by overlapping zones, thereby maintaining the overall physiological stability. The Wnt/β-catenin signaling pathway is well known for its role in liver growth, development, and regeneration. In addition, the Wnt/β-catenin pathway plays a fundamental and dominant role in hepatic zonation and signals to orchestrate various functions of liver metabolism and pathophysiology. The β-catenin protein is the central player in the Wnt/β-catenin signaling cascade, and its activation is crucial for metabolic patterning of the liver. However, dysregulation of Wnt/β-catenin signaling is also implicated in different liver pathologies, including those associated with metabolic syndrome. β-Catenin is preferentially localized in the central region of the hepatic lobule surrounding the central vein and regulates multiple functions of this region. This review outlines the role of Wnt/β-catenin signaling pathway in controlling the different metabolic processes surrounding the central vein and its relation to liver homeostasis and dysfunction.

Spatial Transcriptomics to define transcriptional patterns of zonation and structural components in the mouse liver

Reconstruction of heterogeneity through single cell transcriptional profiling has greatly advanced our understanding of the spatial liver transcriptome in recent years. However, global transcriptional differences across lobular units remain elusive in physical space. Here, we apply Spatial Transcriptomics to perform transcriptomic analysis across sectioned liver tissue. We confirm that the heterogeneity in this complex tissue is predominantly determined by lobular zonation. By introducing novel computational approaches, we enable transcriptional gradient measurements between tissue structures, including several lobules in a variety of orientations. Further, our data suggests the presence of previously transcriptionally uncharacterized structures within liver tissue, contributing to the overall spatial heterogeneity of the organ. This study demonstrates how comprehensive spatial transcriptomic technologies can be used to delineate extensive spatial gene expression patterns in the liver, indicating its future impact for studies of liver function, development and regeneration as well as its potential in pre-clinical and clinical pathology.

Characterization of hepatic zonation in mice by mass-spectrometric and antibody-based proteomics approaches

Periportal and perivenous hepatocytes show zonal heterogeneity in metabolism and signaling. Here, hepatic zonation in mouse liver was analyzed by non-targeted mass spectrometry (MS) and by the antibody-based DigiWest technique, yielding a comprehensive overview of protein expression in periportal and perivenous hepatocytes. Targeted immunoaffinity-based proteomics were used to substantiate findings related to drug metabolism. 165 (MS) and 82 (DigiWest) zonated proteins were identified based on the selected criteria for statistical significance, including 7 (MS) and 43 (DigiWest) proteins not identified as zonated before. New zonated proteins especially comprised kinases and phosphatases related to growth factor-dependent signaling, with mainly periportal localization. Moreover, the mainly perivenous zonation of a large panel of cytochrome P450 enzymes was characterized. DigiWest data were shown to complement the MS results, substantially improving possibilities to bioinformatically identify zonated biological processes. Data mining revealed key regulators and pathways preferentially active in either periportal or perivenous hepatocytes, with β-catenin signaling and nuclear xeno-sensing receptors as the most prominent perivenous regulators, and several kinase- and G-protein-dependent signaling cascades active mainly in periportal hepatocytes. In summary, the present data substantially broaden our knowledge of hepatic zonation in mouse liver at the protein level.

Survival-Assured Liver Injury Preconditioning (SALIC) Enables Robust Expansion of Human Hepatocytes in Fah-/- Rag2-/- IL2rg-/- Rats

Although liver-humanized animals are desirable tools for drug development and expansion of human hepatocytes in large quantities, their development is restricted to mice. In animals larger than mice, a precondition for efficient liver humanization remains preliminary because of different xeno-repopulation kinetics in livers of larger sizes. Since rats are ten times larger than mice and widely used in pharmacological studies, liver-humanized rats are more preferable. Here, Fah-/- Rag2-/- IL2rg-/- (FRG) rats are generated by CRISPR/Cas9, showing accelerated liver failure and lagged liver xeno-repopulation compared to FRG mice. A survival-assured liver injury preconditioning (SALIC) protocol, which consists of retrorsine pretreatment and cycling 2-(2-nitro-4-trifluoromethylbenzoyl)-1,3-cyclohexanedione (NTBC) administration by defined concentrations and time intervals, is developed to reduce the mortality of FRG rats and induce a regenerative microenvironment for xeno-repopulation. Human hepatocyte repopulation is boosted to 31 ± 4% in rat livers at 7 months after transplantation, equivalent to approximately a 1200-fold expansion. Human liver features of transcriptome and zonation are reproduced in humanized rats. Remarkably, they provide sufficient samples for the pharmacokinetic profiling of human-specific metabolites. This model is thus preferred for pharmacological studies and human hepatocyte production. SALIC may also be informative to hepatocyte transplantation in other large-sized species.

Liver Zonation - Revisiting Old Questions With New Technologies

Despite the ever-increasing prevalence of non-alcoholic fatty liver disease (NAFLD), the etiology and pathogenesis remain poorly understood. This is due, in part, to the liver's complex physiology and architecture. The liver maintains glucose and lipid homeostasis by coordinating numerous metabolic processes with great efficiency. This is made possible by the spatial compartmentalization of metabolic pathways a phenomenon known as liver zonation. Despite the importance of zonation to normal liver function, it is unresolved if and how perturbations to liver zonation can drive hepatic pathophysiology and NAFLD development. While hepatocyte heterogeneity has been identified over a century ago, its examination had been severely hindered due to technological limitations. Recent advances in single cell analysis and imaging technologies now permit further characterization of cells across the liver lobule. This review summarizes the advances in examining liver zonation and elucidating its regulatory role in liver physiology and pathology. Understanding the spatial organization of metabolism is vital to further our knowledge of liver disease and to provide targeted therapeutic avenues.

Metabolic Contributions of Wnt Signaling: More Than Controlling Flight

The canonical Wnt signaling pathway is ubiquitous throughout the body and influences a diverse array of physiological processes. Following the initial discovery of the Wnt signaling pathway during wing development in Drosophila melanogaster, it is now widely appreciated that active Wnt signaling in mammals is necessary for the development and growth of various tissues involved in whole-body metabolism, such as brain, liver, pancreas, muscle, and adipose. Moreover, elegant gain- and loss-of-function studies have dissected the tissue-specific roles of various downstream effector molecules in the regulation of energy homeostasis. This review attempts to highlight and summarize the contributions of the Wnt signaling pathway and its downstream effectors on whole-body metabolism and their influence on the development of metabolic diseases, such as diabetes and obesity. A better understanding of the Wnt signaling pathway in these tissues may aid in guiding the development of future therapeutics to treat metabolic diseases.

Metabolic nuclear receptors coordinate energy metabolism to regulate Sox9+ hepatocyte fate

Recent research has indicated the adult liver Sox9⁺ cells located in the portal triads contribute to the physiological maintenance of liver mass and injury repair. However, the physiology and pathology regulation mechanisms of adult liver Sox9⁺ cells remain unknown. Here, PPARα and FXR bound to the shared site in Sox9 promoter with opposite transcriptional outputs. PPARα activation enhanced the fatty acid β-oxidation, oxidative phosphorylation (OXPHOS), and adenosine triphosphate (ATP) production, thus promoting proliferation and differentiation of Sox9⁺ hepatocytes along periportal (PP)-perivenous (PV) axis. However, FXR activation increased glycolysis but decreased OXPHOS and ATP production, therefore preventing proliferation of Sox9⁺ hepatocytes along PP-PV axis by promoting Sox9⁺ hepatocyte self-renewal. Our research indicates that metabolic nuclear receptors play critical roles in liver progenitor Sox9⁺ hepatocyte homeostasis to initiate or terminate liver injury-induced cell proliferation and differentiation, suggesting that PPARα and FXR are potential therapeutic targets for modulating liver regeneration.

•

PPARα promotes Sox9 expression and FXR inhibits Sox9 expression

•

PPARα promotes proliferation and differentiation of Sox9⁺ hepatocytes

•

FXR promotes Sox9⁺ hepatocyte self-renewal

•

PPARα and FXR coordinate energy metabolism to regulate Sox9⁺ hepatocyte fate

Spatial Reconstruction of the Early Hepatic Transcriptomic Landscape After an Acetaminophen Overdose Using Single-Cell RNA-Sequencing

An acetaminophen (APAP) overdose is the most common cause of acute liver failure in the United States. A hallmark characteristic of APAP hepatotoxicity is centrilobular necrosis. General, innate mechanisms such as lower amounts of GSH and higher cytochrome P450 2e1 expression in pericentral (PC) hepatocytes are known to contribute to the differences in susceptibility to cell injury between periportal (PP) hepatocytes and PC hepatocytes. Although a sequence of molecular events involving formation of the reactive metabolite N-acetyl-p-benzoquinone imine, GSH depletion, oxidative stress, and c-Jun N-terminal kinase activation define the early cell stress trajectory following APAP exposure, their activation in PC versus PP hepatocytes is not well characterized. By using single-cell RNA-sequencing, we provide the first reconstruction of the early transcriptomic APAP liver lobule after validation of our methodology using human liver single-cell RNA-sequencing data. Two hours after APAP treatment, we find that PP hepatocytes progress along the APAP stress axis to oxidative stress, before resolving injury due to innate and adaptive mechanisms. However, PC hepatocytes continue along this stress axis as indicated by activation of mitogen-activated protein kinase genes, which is absent in PP hepatocytes. We also identify a population of glutamine synthetase enriched PC hepatocytes in close proximity to the central vein, where a stepwise induction of a stress program culminated in cell death. Collectively, these findings elucidate a molecular sequence of events distinguishing the differential response to APAP exposure between PP and PC hepatocytes and identify a subset of uniquely susceptible PC hepatocytes.

Hormonal regulation of glycine decarboxylase and its relationship to oxidative stress

In both humans and rodent models, circulating glycine levels are significantly reduced in obesity, glucose intolerance, type II diabetes, and non-alcoholic fatty liver disease. The glycine cleavage system and its rate-limiting enzyme, glycine decarboxylase (GLDC), is a major determinant of plasma glycine levels. The goals of this study were to determine if the increased expression of GLDC contributes to the reduced plasma glycine levels seen in disease states, to characterize the hormonal regulation of GLDC gene expression, and to determine if altered GLDC expression has physiological effects that might affect the development of diabetes. The findings presented here show that hepatic GLDC gene expression is elevated in mouse models of obesity and diabetes, as well as by fasting. We demonstrated that GLDC gene expression is strongly regulated by the metabolic hormones glucagon and insulin, and we identified the signaling pathways involved in this regulation. Finally, we found that GLDC expression is linked to glutathione levels, with increased expression associated with elevated levels of glutathione and reduced expression associated with a suppression of glutathione and increased cellular ROS levels. These findings suggest that the hormonal regulation of GLDC contributes not only to the changes in circulating glycine levels seen in metabolic disease, but also affects glutathione production, possibly as a defense against metabolic disease-associated oxidative stress.

Single-cell RNA sequencing of human liver reveals hepatic stellate cell heterogeneity

Background & Aims

The multiple vital functions of the human liver are performed by highly specialised parenchymal and non-parenchymal cells organised in complex collaborative sinusoidal units. Although crucial for homeostasis, the cellular make-up of the human liver remains to be fully elucidated. Here, single-cell RNA-sequencing was used to unravel the heterogeneity of human liver cells, in particular of hepatocytes (HEPs) and hepatic stellate cells (HSCs).

Method

The transcriptome of ~25,000 freshly isolated human liver cells was profiled using droplet-based RNA-sequencing. Recently published data sets and RNA in situ hybridisation were integrated to validate and locate newly identified cell populations.

Results

In total, 22 cell populations were annotated that reflected the heterogeneity of human parenchymal and non-parenchymal liver cells. More than 20,000 HEPs were ordered along the portocentral axis to confirm known, and reveal previously undescribed, zonated liver functions. The existence of 2 subpopulations of human HSCs with unique gene expression signatures and distinct intralobular localisation was revealed (i.e. portal and central vein-concentrated GPC3⁺ HSCs and perisinusoidally located DBH⁺ HSCs). In particular, these data suggest that, although both subpopulations collaborate in the production and organisation of extracellular matrix, GPC3⁺ HSCs specifically express genes involved in the metabolism of glycosaminoglycans, whereas DBH⁺ HSCs display a gene signature that is reminiscent of antigen-presenting cells.

Conclusions

This study highlights metabolic zonation as a key determinant of HEP transcriptomic heterogeneity and, for the first time, outlines the existence of heterogeneous HSC subpopulations in the human liver. These findings call for further research on the functional implications of liver cell heterogeneity in health and disease.

Lay summary

This study resolves the cellular landscape of the human liver in an unbiased manner and at high resolution to provide new insights into human liver cell biology. The results highlight the physiological heterogeneity of human hepatic stellate cells.

•

A cell atlas from the near-native transcriptome of >25,000 human liver cells is presented.

•

Hepatocytes were ordered along the portocentral axis to reveal previously undescribed gene expression patterns and zonated liver functions.

•

Two subpopulations of human hepatic stellate cells (HSCs) are reported, characterised by different spatial distribution in the native tissue.

•

Characteristic gene signatures of HSC subpopulations are suggestive of far-reaching functional differences.

Functional Consequences of Metabolic Zonation in Murine Livers: Insights for an Old Story

BACKGROUND AND AIMS

Zone-dependent differences in expression of metabolic enzymes along the portocentral axis of the acinus are a long-known feature of liver metabolism. A prominent example is the preferential localization of the enzyme, glutamine synthetase, in pericentral hepatocytes, where it converts potentially toxic ammonia to the valuable amino acid, glutamine. However, with the exception of a few key regulatory enzymes, a comprehensive and quantitative assessment of zonal differences in the abundance of metabolic enzymes and, much more important, an estimation of the associated functional differences between portal and central hepatocytes is missing thus far.

APPROACH AND RESULTS

We addressed this problem by establishing a method for the separation of periportal and pericentral hepatocytes that yields sufficiently pure fractions of both cell populations. Quantitative shotgun proteomics identified hundreds of differentially expressed enzymes in the two cell populations. We used zone-specific proteomics data for scaling of the maximal activities to generate portal and central instantiations of a comprehensive kinetic model of central hepatic metabolism (Hepatokin1).

CONCLUSIONS

The model simulations revealed significant portal-to-central differences in almost all metabolic pathways involving carbohydrates, fatty acids, amino acids, and detoxification.

Commensal-driven immune zonation of the liver promotes host defence

The liver connects the intestinal portal vasculature with the general circulation, using a diverse array of immune cells to protect from pathogens that translocate from the gut¹. In liver lobules, blood flows from portal triads that are situated in periportal lobular regions to the central vein via a polarized sinusoidal network. Despite this asymmetry, resident immune cells in the liver are considered to be broadly dispersed across the lobule. This differs from lymphoid organs, in which immune cells adopt spatially biased positions to promote effective host defence^2,3. Here we used quantitative multiplex imaging, genetic perturbations, transcriptomics, infection-based assays and mathematical modelling to reassess the relationship between the localization of immune cells in the liver and host protection. We found that myeloid and lymphoid resident immune cells concentrate around periportal regions. This asymmetric localization was not developmentally controlled, but resulted from sustained MYD88-dependent signalling induced by commensal bacteria in liver sinusoidal endothelial cells, which in turn regulated the composition of the pericellular matrix involved in the formation of chemokine gradients. In vivo experiments and modelling showed that this immune spatial polarization was more efficient than a uniform distribution in protecting against systemic bacterial dissemination. Together, these data reveal that liver sinusoidal endothelial cells sense the microbiome, actively orchestrating the localization of immune cells, to optimize host defence.

Regulation of CAR and PXR Expression in Health and Disease

Pregnane X receptor (PXR, NR1I2) and constitutive androstane receptor (CAR, NR1I3) are members of the nuclear receptor superfamily that mainly act as ligand-activated transcription factors. Their functions have long been associated with the regulation of drug metabolism and disposition, and it is now well established that they are implicated in physiological and pathological conditions. Considerable efforts have been made to understand the regulation of their activity by their cognate ligand; however, additional regulatory mechanisms, among which the regulation of their expression, modulate their pleiotropic effects. This review summarizes the current knowledge on CAR and PXR expression during development and adult life; tissue distribution; spatial, temporal, and metabolic regulations; as well as in pathological situations, including chronic diseases and cancers. The expression of CAR and PXR is modulated by complex regulatory mechanisms that involve the interplay of transcription factors and also post-transcriptional and epigenetic modifications. Moreover, many environmental stimuli affect CAR and PXR expression through mechanisms that have not been elucidated.

β-catenin signaling, the constitutive androstane receptor and their mutual interactions

Aberrant signaling through β-catenin is an important determinant of tumorigenesis in rodents as well as in humans. In mice, xenobiotic activators of the constitutive androstane receptor (CAR), a chemo-sensing nuclear receptor, promote liver tumor growth by means of a non-genotoxic mechanism and, under certain conditions, select for hepatocellular tumors which contain activated β-catenin. In normal hepatocytes, interactions of β-catenin and CAR have been demonstrated with respect to the induction of proliferation and drug metabolism-related gene expression. The molecular details of these interactions are still not well understood. Recently it has been hypothesized that CAR might activate β-catenin signaling, thus providing a possible explanation for some of the observed phenomena. Nonetheless, many aspects of the molecular interplay of the two regulators have still not been elucidated. This review briefly summarizes our current knowledge about the interplay of CAR and β-catenin. By taking into account data and observations obtained with different mouse models and employing different experimental approaches, it is shown that published data also contain substantial evidence that xenobiotic activators of CAR do not activate, or do even inhibit signaling through the β-catenin pathway. The review highlights new aspects of possible ways of interaction between the two signaling cascades and will help to stimulate scientific discussion about the crosstalk of β-catenin signaling and the nuclear receptor CAR.

Generation of metabolically functional hepatocyte-like cells from dedifferentiated fat cells by Foxa2, Hnf4a and Sall1 transduction

Mature adipocyte-derived dedifferentiated fat (DFAT) cells have been identified to possess similar multipotency to mesenchymal stem cells, but a method for converting DFAT cells into hepatocytes was previously unknown. Here, using comprehensive analysis of gene expression profiles, we have extracted three transcription factors, namely Foxa2, Hnf4a and Sall1 (FHS), that can convert DFAT cells into hepatocytes. Hepatogenic induction has converted FHS-infected DFAT cells into an epithelial-like morphological state and promoted the expression of hepatocyte-specific features. Furthermore, the DFAT-derived hepatocyte-like (D-Hep) cells catalyzed the detoxification of several compounds. These results indicate that the transduction of DFAT cells with three genes, which were extracted by comprehensive gene expression analysis, efficiently generated D-Hep cells with detoxification abilities similar to those of primary hepatocytes. Thus, D-Hep cells may be useful as a new cell source for surrogate hepatocytes and may be applied to drug discovery studies, such as hepatotoxicity screening and drug metabolism tests.

Reproducibility across single-cell RNA-seq protocols for spatial ordering analysis

As newer single-cell protocols generate increasingly more cells at reduced sequencing depths, the value of a higher read depth may be overlooked. Using data from three different single-cell RNA-seq protocols that lend themselves to having either higher read depth (Smart-seq) or many cells (MARS-seq and 10X), we evaluate their ability to recapitulate biological signals in the context of spatial reconstruction. Overall, we find gene expression profiles after spatial reconstruction analysis are highly reproducible between datasets despite being generated by different protocols and using different computational algorithms. While UMI-based protocols such as 10X and MARS-seq allow for capturing more cells, Smart-seq's higher sensitivity and read-depth allow for analysis of lower expressed genes and isoforms. Additionally, we evaluate trade-offs for each protocol by performing subsampling analyses and find that optimizing the balance between sequencing depth and number of cells within a protocol is necessary for efficient use of resources. Our analysis emphasizes the importance of selecting a protocol based on the biological questions and features of interest.