Single-cell spatial reconstruction reveals global division of labour in the mammalian liver. Nature. 2017;542:352-356. [PMID: 28166538 PMCID: PMC5321580 DOI: 10.1038/nature21065]

Isolation of Nuclei from Human Snap-frozen Liver Tissue for Single-nucleus RNA Sequencing

Single-nucleus RNA sequencing (snRNA-seq) provides a powerful tool for studying cell type composition in heterogenous tissues. The liver is a vital organ composed of a diverse set of cell types; thus, single-cell technologies could greatly facilitate the deconvolution of liver tissue composition and various downstream omics analyses at the cell-type level. Applying single-cell technologies to fresh liver biopsies can, however, be very challenging, and snRNA-seq of snap-frozen liver biopsies requires some optimization given the high nucleic acid content of the solid liver tissue. Therefore, an optimized protocol for snRNA-seq specifically targeted for the use of frozen liver samples is needed to improve our understanding of human liver gene expression at the cell-type resolution. We present a protocol for performing nuclei isolation from snap-frozen liver tissues, as well as guidance on the application of snRNA-seq. We also provide guidance on optimizing the protocol to different tissue and sample types.

ScRNA-seq and ST-seq in liver research

Spatial transcriptomics, which combine gene expression data with spatial information, has quickly expanded in recent years. With application of this method in liver research, our knowledge about liver development, regeneration, and diseases have been greatly improved. While this field is moving forward, a variety of problems still need to be addressed, including sensitivity, limited capacity to obtain exact single-cell information, data processing methods, as well as others. Methods like single-cell RNA sequencing (scRNA-seq) are usually used together with spatial transcriptome sequencing (ST-seq) to clarify cell-specific gene expression. In this review, we explore how advances of scRNA-seq and ST-seq, especially ST-seq, will pave the way to new opportunities to investigate fundamental questions in liver research. Finally, we will discuss the strengths, limitations, and future perspectives of ST-seq in liver research.

Loss of liver function in chronic liver disease: An identity crisis

Adult hepatocyte identity is constructed throughout embryonic development and fine-tuned after birth. A multinodular network of transcription factors, along with pre-mRNA splicing regulators, define the transcriptome, which encodes the proteins needed to perform the complex metabolic and secretory functions of the mature liver. Transient hepatocellular dedifferentiation can occur as part of the regenerative mechanisms triggered in response to acute liver injury. However, persistent downregulation of key identity genes is now accepted as a strong determinant of organ dysfunction in chronic liver disease, a major global health burden. Therefore, the identification of core transcription factors and splicing regulators that preserve hepatocellular phenotype, and a thorough understanding of how these networks become disrupted in diseased hepatocytes, is of high clinical relevance. In this context, we review the key players in liver differentiation and discuss in detail critical factors, such as HNF4α, whose impairment mediates the breakdown of liver function. Moreover, we present compelling experimental evidence demonstrating that restoration of core transcription factor expression in a chronically injured liver can reset hepatocellular identity, improve function and ameliorate structural abnormalities. The possibility of correcting the phenotype of severely damaged and malfunctional livers may reveal new therapeutic opportunities for individuals with cirrhosis and advanced liver disease.

Liver Organoids, Novel and Promising Modalities for Exploring and Repairing Liver Injury

The past decades have witnessed great advances in organoid technology. Liver is the biggest solid organ, performing multifaceted physiological functions. Nowadays, liver organoids have been applied in many fields including pharmaceutical research, precision medicine and disease models. Compared to traditional 2-dimensional cell line cultures and animal models, liver organoids showed the unique advantages. More importantly, liver organoids can well model the features of the liver and tend to be novel and promising modalities for exploring liver injury, thus finding potential treatment targets and repairing liver injury. In this review, we reviewed the history of the development of liver organoids and summarized the application of liver organoids and recent studies using organoids to explore and further repair the liver injury. These novel modalities could provide new insights about the process of liver injury.

Delineating Plasmodium liver infection across space and time

The liver is a major entry point and gatekeeper for invasive pathogens. However, high-resolution, spatiotemporal transcriptomic analysis of host-pathogen interactions has remained challenging. Afriat et al. have deconvoluted Plasmodium berghei liver-stage maturation at an unprecedented scale and discovered molecular signatures of heterogeneity during pre-erythrocytic development of malarial parasites.

Intermittent fasting induces rapid hepatocyte proliferation to restore the hepatostat in the mouse liver

Nutrient availability fluctuates in most natural populations, forcing organisms to undergo periods of fasting and re-feeding. It is unknown how dietary changes influence liver homeostasis. Here, we show that a switch from ad libitum feeding to intermittent fasting (IF) promotes rapid hepatocyte proliferation. Mechanistically, IF-induced hepatocyte proliferation is driven by the combined action of systemic FGF15 and localized WNT signaling. Hepatocyte proliferation during periods of fasting and re-feeding re-establishes a constant liver-to-body mass ratio, thus maintaining the hepatostat. This study provides the first example of dietary influence on adult hepatocyte proliferation and challenges the widely held view that liver tissue is mostly quiescent unless chemically or mechanically injured.

Single-cell transcriptomics shows dose-dependent disruption of hepatic zonation by TCDD in mice

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) dose-dependently induces the development of hepatic fat accumulation and inflammation with fibrosis in mice initially in the portal region. Conversely, differential gene and protein expression is first detected in the central region. To further investigate cell-specific and spatially resolved dose-dependent changes in gene expression elicited by TCDD, single-nuclei RNA sequencing and spatial transcriptomics were used for livers of male mice gavaged with TCDD every 4 days for 28 days. The proportion of 11 cell (sub)types across 131 613 nuclei dose-dependently changed with 68% of all portal and central hepatocyte nuclei in control mice being overtaken by macrophages following TCDD treatment. We identified 368 (portal fibroblasts) to 1339 (macrophages) differentially expressed genes. Spatial analyses revealed initial loss of portal identity that eventually spanned the entire liver lobule with increasing dose. Induction of R-spondin 3 (Rspo3) and pericentral Apc, suggested dysregulation of the Wnt/β-catenin signaling cascade in zonally resolved steatosis. Collectively, the integrated results suggest disruption of zonation contributes to the pattern of TCDD-elicited NAFLD pathologies.

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.

Murine breast cancers disorganize the liver transcriptome in a zonated manner

The spatially organized gene expression program within the liver specifies hepatocyte functions according to their relative distances to the bloodstream (i.e., zonation), contributing to liver homeostasis. Despite the knowledge that solid cancers remotely disrupt liver homeostasis, it remains unexplored whether solid cancers affect liver zonation. Here, using spatial transcriptomics, we thoroughly investigate the abundance and zonation of hepatic genes in cancer-bearing mice. We find that breast cancers affect liver zonation in various distinct manners depending on biological pathways. Aspartate metabolism and triglyceride catabolic processes retain relatively intact zonation patterns, but the zonation of xenobiotic catabolic process genes exhibits a strong disruption. The acute phase response is induced in zonated manners. Furthermore, we demonstrate that breast cancers activate innate immune cells in particular neutrophils in distinct zonated manners, rather than in a uniform fashion within the liver. Collectively, breast cancers disorganize hepatic transcriptomes in zonated manners, thereby disrupting zonated functions of the liver.

Mapping transcriptional heterogeneity and metabolic networks in fatty livers at single-cell resolution

Non-alcoholic fatty liver disease is a heterogeneous disease with unclear underlying molecular mechanisms. Here, we perform single-cell RNA sequencing of hepatocytes and hepatic non-parenchymal cells to map the lipid signatures in mice with non-alcoholic fatty liver disease (NAFLD). We uncover previously unidentified clusters of hepatocytes characterized by either high or low srebp1 expression. Surprisingly, the canonical lipid synthesis driver Srebp1 is not predictive of hepatic lipid accumulation, suggestive of other drivers of lipid metabolism. By combining transcriptional data at single-cell resolution with computational network analyses, we find that NAFLD is associated with high constitutive androstane receptor (CAR) expression. Mechanistically, CAR interacts with four functional modules: cholesterol homeostasis, bile acid metabolism, fatty acid metabolism, and estrogen response. Nuclear expression of CAR positively correlates with steatohepatitis in human livers. These findings demonstrate significant cellular differences in lipid signatures and identify functional networks linked to hepatic steatosis in mice and humans.

HepaRG cells undergo increased levels of post-differentiation patterning in physiologic conditions when maintained as 3D cultures in paper-based scaffolds

Monolayer cultures of hepatocytes lack many aspects of the liver sinusoid, including a tissue-level organization that results from extracellular matrix interactions and gradients of soluble molecules that span from the portal triad to the central vein. We measured the activity and transcript levels of drug-metabolizing enzymes in HepaRG cells maintained in three different culture configurations: as monolayers, seeded onto paper scaffolds that were pre-loaded with a collagen matrix, and when seeded directly into the paper scaffolds as a cell-laden gel. Drug metabolism was significantly decreased in the presence of the paper scaffolds compared to monolayer configurations when cells were exposed to standard culture conditions. Despite this decreased function, transcript levels suggest the cells undergo increased polarization and adopt a biliary-like character in the paper scaffolds, including the increased expression of transporter proteins (e.g., ABCB11 and SLOC1B1) and the KRT19 cholangiocyte marker. When exposed to representative periportal or perivenous culture conditions, we observed in vivo zonal-like patterns, including increased cytochrome P450 (CYP) activity and transcript levels in the perivenous condition. This increased CYP activity is more pronounced in the laden configuration, supporting the need to include multiple aspects of the liver microenvironment to observe the post-differentiation processing of hepatocytes.

Liver sinusoidal endothelial cells induce BMP6 expression in response to non-transferrin-bound iron

Homeostatic adaptation to systemic iron overload involves transcriptional induction of bone morphogenetic protein 6 (BMP6) in liver sinusoidal endothelial cells (LSECs). BMP6 is then secreted to activate signaling of the iron hormone hepcidin (HAMP) in neighboring hepatocytes. To explore the mechanism of iron sensing by LSECs, we generated TfrcTek-Cre mice with endothelial cell-specific ablation of transferrin receptor 1 (Tfr1). We also used control Tfrcfl/fl mice to characterize the LSEC-specific molecular responses to iron using single-cell transcriptomics. TfrcTek-Cre animals tended to have modestly increased liver iron content (LIC) compared with Tfrcfl/fl controls but expressed physiological Bmp6 and Hamp messenger RNA (mRNA). Despite a transient inability to upregulate Bmp6, they eventually respond to iron challenges with Bmp6 and Hamp induction, yet occasionally to levels slightly lower relative to LIC. High dietary iron intake triggered the accumulation of serum nontransferrin bound iron (NTBI), which significantly correlated with liver Bmp6 and Hamp mRNA levels and elicited more profound alterations in the LSEC transcriptome than holo-transferrin injection. This culminated in the robust induction of Bmp6 and other nuclear factor erythroid 2-related factor 2 (Nrf2) target genes, as well as Myc target genes involved in ribosomal biogenesis and protein synthesis. LSECs and midzonal hepatocytes were the most responsive liver cells to iron challenges and exhibited the highest expression of Bmp6 and Hamp mRNAs, respectively. Our data suggest that during systemic iron overload, LSECs internalize NTBI, which promotes oxidative stress and thereby transcriptionally induces Bmp6 via Nrf2. Tfr1 appears to contribute to iron sensing by LSECs, mostly under low iron conditions.

HepaRG cells undergo increased levels of post-differentiation patterning in physiologic conditions when maintained as 3D cultures in paper-based scaffolds

Monolayer cultures of hepatocytes lack many aspects of the liver sinusoid, including a tissue-level organization that results from extracellular matrix interactions and gradients of soluble molecules that span from the portal triad to the central vein. We measured the activity and transcript levels of drug-metabolizing enzymes in HepaRG cells maintained in three different culture configurations: as monolayers, seeded onto paper scaffolds that were pre-loaded with a collagen matrix, and when seeded directly into the paper scaffolds as a cell-laden gel. Drug metabolism was significantly decreased in the presence of the paper scaffolds compared to monolayer configurations when cells were exposed to standard culture conditions. Despite this decreased function, transcript levels suggest the cells undergo increased polarization and adopt a biliary-like character in the paper scaffolds, including the increased expression of transporter proteins (e.g., ABCB11 and SLOC1B1) and the KRT19 cholangiocyte marker. When exposed to representative periportal or perivenous culture conditions, we observed in vivo zonal-like patterns, including increased cytochrome P450 (CYP) activity and transcript levels in the perivenous condition. This increased CYP activity is more pronounced in the laden configuration, supporting the need to include multiple aspects of the liver microenvironment to observe the post-differentiation processing of hepatocytes.

Catecholamine induces Kupffer cell apoptosis via growth differentiation factor 15 in alcohol-associated liver disease

Chronic alcohol consumption often induces hepatic steatosis but rarely causes severe inflammation in Kupffer cells (KCs) despite the increased hepatic influx of lipopolysaccharide (LPS), suggesting the presence of a veiled tolerance mechanism. In addition to LPS, the liver is affected by several gut-derived neurotransmitters through the portal blood, but the effects of catecholamines on KCs have not been clearly explored in alcohol-associated liver disease (ALD). Hence, we investigated the regulatory roles of catecholamine on inflammatory KCs under chronic alcohol exposure. We discovered that catecholamine levels were significantly elevated in the cecum, portal blood, and liver tissues of chronic ethanol-fed mice. Increased catecholamines induced mitochondrial translocation of cytochrome P450 2E1 in perivenous hepatocytes expressing the β2-adrenergic receptor (ADRB2), leading to the enhanced production of growth differentiation factor 15 (GDF15). Subsequently, GDF15 profoundly increased ADRB2 expression in adjacent inflammatory KCs to facilitate catecholamine/ADRB2-mediated apoptosis. Single-cell RNA sequencing of KCs confirmed the elevated expression of Adrb2 and apoptotic genes after chronic ethanol intake. Genetic ablation of Adrb2 or hepatic Gdf15 robustly decreased the number of apoptotic KCs near perivenous areas, exacerbating alcohol-associated inflammation. Consistently, we found that blood and stool catecholamine levels and perivenous GDF15 expression were increased in patients with early-stage ALD along with an increase in apoptotic KCs. Our findings reveal a novel protective mechanism against ALD, in which the catecholamine/GDF15 axis plays a critical role in KC apoptosis, and identify a unique neuro-metabo-immune axis between the gut and liver that elicits hepatoprotection against alcohol-mediated pathogenic challenges.

Biliary epithelial cells are facultative liver stem cells during liver regeneration in adult zebrafish

The liver is a highly regenerative organ, yet the presence of a dedicated stem cell population remains controversial. Here, we interrogate a severe hepatocyte injury model in adult zebrafish to define that regeneration involves a stem cell population. After near-total hepatocyte ablation, single-cell transcriptomic and high-resolution imaging analyses throughout the entire regenerative timeline reveal that biliary epithelial cells undergo transcriptional and morphological changes to become hepatocytes. As a population, biliary epithelial cells give rise to both hepatocytes and biliary epithelial cells. Biliary epithelial cells proliferate and dedifferentiate to express hepatoblast transcription factors prior to hepatocyte differentiation. This process is characterized by increased MAPK, PI3K, and mTOR signaling, and chemical inhibition of these pathways impairs biliary epithelial cell proliferation and fate conversion. We conclude that, upon severe hepatocyte ablation in the adult liver, biliary epithelial cells act as facultative liver stem cells in an EGFR-PI3K-mTOR-dependent manner.

TCDD dysregulation of lncRNA expression, liver zonation and intercellular communication across the liver lobule

The persistent environmental aryl hydrocarbon receptor agonist and hepatotoxin TCDD (2,3,7,8-tetrachlorodibenzo- p -dioxin) induces hepatic lipid accumulation (steatosis), inflammation (steatohepatitis) and fibrosis. Thousands of liver-expressed, nuclear-localized lncRNAs with regulatory potential have been identified; however, their roles in TCDD-induced hepatoxicity and liver disease are unknown. We analyzed single nucleus (sn)RNA-seq data from control and chronic TCDD-exposed mouse liver to determine liver cell-type specificity, zonation and differential expression profiles for thousands of IncRNAs. TCDD dysregulated >4,000 of these lncRNAs in one or more liver cell types, including 684 lncRNAs specifically dysregulated in liver non-parenchymal cells. Trajectory inference analysis revealed major disruption by TCDD of hepatocyte zonation, affecting >800 genes, including 121 IncRNAs, with strong enrichment for lipid metabolism genes. TCDD also dysregulated expression of >200 transcription factors, including 19 Nuclear Receptors, most notably in hepatocytes and Kupffer cells. TCDD-induced changes in cellâ€"cell communication patterns included marked decreases in EGF signaling from hepatocytes to non-parenchymal cells and increases in extracellular matrix-receptor interactions central to liver fibrosis. Gene regulatory networks constructed from the snRNA-seq data identified TCDD-exposed liver network-essential lncRNA regulators linked to functions such as fatty acid metabolic process, peroxisome and xenobiotic metabolic. Networks were validated by the striking enrichments that predicted regulatory IncRNAs showed for specific biological pathways. These findings highlight the power of snRNA-seq to discover functional roles for many xenobiotic-responsive lncRNAs in both hepatocytes and liver non-parenchymal cells and to elucidate novel aspects of foreign chemical-induced hepatotoxicity and liver disease, including dysregulation of intercellular communication within the liver lobule.

Hepatocyte Smoothened Activity Controls Susceptibility to Insulin Resistance and Nonalcoholic Fatty Liver Disease

BACKGROUND & AIMS

Nonalcoholic steatohepatitis (NASH), a leading cause of cirrhosis, strongly associates with the metabolic syndrome, an insulin-resistant proinflammatory state that disrupts energy balance and promotes progressive liver degeneration. We aimed to define the role of Smoothened (Smo), an obligatory component of the Hedgehog signaling pathway, in controlling hepatocyte metabolic homeostasis and, thereby, susceptibility to NASH.

METHODS

We conditionally deleted Smo in hepatocytes of healthy chow-fed mice and performed metabolic phenotyping, coupled with single-cell RNA sequencing (RNA-seq), to characterize the role of hepatocyte Smo in regulating basal hepatic and systemic metabolic homeostasis. Liver RNA-seq datasets from 2 large human cohorts were also analyzed to define the relationship between Smo and NASH susceptibility in people.

RESULTS

Hepatocyte Smo deletion inhibited the Hedgehog pathway and promoted fatty liver, hyperinsulinemia, and insulin resistance. We identified a plausible mechanism whereby inactivation of Smo stimulated the mTORC1-SREBP1c signaling axis, which promoted lipogenesis while inhibiting the hepatic insulin cascade. Transcriptomics of bulk and single Smo-deficient hepatocytes supported suppression of insulin signaling and also revealed molecular abnormalities associated with oxidative stress and mitochondrial dysfunction. Analysis of human bulk RNA-seq data revealed that Smo expression was (1) highest in healthy livers, (2) lower in livers with NASH than in those with simple steatosis, (3) negatively correlated with markers of insulin resistance and liver injury, and (4) declined progressively as fibrosis severity worsened.

CONCLUSIONS

The Hedgehog pathway controls insulin sensitivity and energy homeostasis in adult livers. Loss of hepatocyte Hedgehog activity induces hepatic and systemic metabolic stress and enhances susceptibility to NASH by promoting hepatic lipoxicity and insulin resistance.

Decoding liver fibrogenesis with single-cell technologies

Liver fibrogenesis is a highly dynamic and complex process that drives the progression of chronic liver disease toward liver failure and end-stage liver diseases. Despite decades of intense studies, the cellular and molecular mechanisms underlying liver fibrogenesis remain elusive, and no approved therapies to treat liver fibrosis are currently available. The rapid development of single-cell RNA sequencing (scRNA-seq) technologies allows the characterization of cellular alterations under healthy and diseased conditions at an unprecedented resolution. In this Review, we discuss how the scRNA-seq studies are transforming our understanding of the regulatory mechanisms of liver fibrosis. We specifically emphasize discoveries on disease-relevant cell subpopulations, molecular events, and cell interactions on cell types including hepatocytes, liver sinusoidal endothelial cells, myofibroblasts, and macrophages. These discoveries have uncovered critical pathophysiological changes during liver fibrogenesis. Further efforts are urged to fully understand the functional contributions of these changes to liver fibrogenesis, and to translate the new knowledge into effective therapeutic approaches.

Dissection of cellular and molecular mechanisms of aristolochic acid-induced hepatotoxicity via single-cell transcriptomics

Background

Aristolochic acids (AAs), a class of carcinogenic and mutagenic natural products from Aristolochia and Asarum plants, are well-known to be responsible for inducing nephrotoxicity and urothelial carcinoma. Recently, accumulating evidence suggests that exposure to AAs could also induce hepatotoxicity and even hepatocellular carcinoma, though the mechanisms are poorly defined.

Methods

Here, we aimed to dissect the underlying cellular and molecular mechanisms of aristolochic acid I (AAI)-induced hepatotoxicity by using advanced single-cell RNA sequencing (scRNA-seq) and proteomics techniques. We established the first single-cell atlas of mouse livers in response to AAI.

Results

In hepatocytes, our results indicated that AAI activated NF-κB and STAT3 signaling pathways, which may contribute to the inflammatory response and apoptosis. In liver sinusoidal endothelial cells (LSECs), AAI activated multiple oxidative stress and inflammatory associated signaling pathways and induced apoptosis. Importantly, AAI induced infiltration of cytotoxic T cells and activation of proinflammatory macrophage and neutrophil cells in the liver to produce inflammatory cytokines to aggravate inflammation.

Conclusions

Collectively, our study provides novel knowledge of AAs-induced molecular characteristics of hepatotoxicity at a single-cell level and suggests future treatment options for AAs associated hepatotoxicity.

Determination of the nanoparticle- and cell-specific toxicological mechanisms in 3D liver spheroids using scRNAseq analysis

Engineered nanomaterials (ENMs) are commonly used in consumer products, allowing exposure to target organs such as the lung, liver, and skin that could lead to adverse health effects in humans. To better reflect on toxicological effects in liver cells, it is important to consider the contribution of hepatocyte morphology, function, and intercellular interactions in a dynamic 3D microenvironment. Herein, we used a 3D liver spheroid model containing hepatocyte and Kupffer cells (KCs) to study the effects of three different material compositions, namely vanadium pentoxide (V2O5), titanium dioxide (TiO2), or graphene oxide (GO). Additionally, we used single-cell RNA sequencing (scRNAseq) to determine the nanoparticle (NP) and cell-specific toxicological responses. A general finding was that hepatocytes exhibit more variation in gene expression and adaptation of signaling pathways than KCs. TNF-α production tied to the NF-κB pathway was a commonly affected pathway by all NPs while impacts on the metabolic function of hepatocytes were unique to V2O5. V2O5 NPs also showed the largest number of differentially expressed genes in both cell types, many of which are related to pro-inflammatory and apoptotic response pathways. There was also evidence of mitochondrial ROS generation and caspase-1 activation after GO and V2O5 treatment, in association with cytokine production. All considered, this study provides insight into the impact of nanoparticles on gene responses in key liver cell types, providing us with a scRNAseq platform that can be used for high-content screening of nanomaterial impact on the liver, for use in biosafety and biomedical applications.

Immunobiology of the biliary tract system

The biliary tract is a complex tubular organ system spanning from the liver to the duodenum. It is the site of numerous acute and chronic disorders, many of unknown origin, that are often associated with cancer development and for which there are limited treatment options. Cholangiocytes with proinflammatory capacities line the lumen and specialised types of immune cells reside in close proximity. Recent technological breakthroughs now permit spatiotemporal assessments of immune cells within distinct niches and have increased our understanding of immune cell tissue residency. In this review, a comprehensive overview of emerging knowledge on the immunobiology of the biliary tract system is provided, with a particular emphasis on the role of distinct immune cells in biliary disorders.

SOTIP is a versatile method for microenvironment modeling with spatial omics data

The rapidly developing spatial omics generated datasets with diverse scales and modalities. However, most existing methods focus on modeling dynamics of single cells while ignore microenvironments (MEs). Here we present SOTIP (Spatial Omics mulTIPle-task analysis), a versatile method incorporating MEs and their interrelationships into a unified graph. Based on this graph, spatial heterogeneity quantification, spatial domain identification, differential microenvironment analysis, and other downstream tasks can be performed. We validate each module's accuracy, robustness, scalability and interpretability on various spatial omics datasets. In two independent mouse cerebral cortex spatial transcriptomics datasets, we reveal a gradient spatial heterogeneity pattern strongly correlated with the cortical depth. In human triple-negative breast cancer spatial proteomics datasets, we identify molecular polarizations and MEs associated with different patient survivals. Overall, by modeling biologically explainable MEs, SOTIP outperforms state-of-art methods and provides some perspectives for spatial omics data exploration and interpretation.

The crosstalk between parenchymal cells and macrophages: A keeper of tissue homeostasis

Non-parenchymal cells (NPCs) and parenchymal cells (PCs) collectively perform tissue-specific functions. PCs play significant roles and continuously adjust the intrinsic functions and metabolism of organs. Tissue-resident macrophages (TRMs) are crucial members of native NPCs in tissues and are essential for immune defense, tissue repair and development, and homeostasis maintenance. As a plastic-phenotypic and prevalent cluster of NPCs, TRMs dynamically assist PCs in functioning by producing cytokines, inflammatory and anti-inflammatory signals, growth factors, and proteolytic enzymes. Furthermore, the PCs of tissues modulate the functional activity and polarization of TRMs. Dysregulation of the PC-TRM crosstalk axis profoundly impacts many essential physiological functions, including synaptogenesis, gastrointestinal motility and secretion, cardiac pulsation, gas exchange, blood filtration, and metabolic homeostasis. This review focuses on the PC-TRM crosstalk in mammalian vital tissues, along with their interactions with tissue homeostasis maintenance and disorders. Thus, this review highlights the fundamental biological significance of the regulatory network of PC-TRM in tissue homeostasis.

scSemiGAN: a single-cell semi-supervised annotation and dimensionality reduction framework based on generative adversarial network

MOTIVATION

Cell-type annotation plays a crucial role in single-cell RNA-seq (scRNA-seq) data analysis. As more and more well-annotated scRNA-seq reference data are publicly available, automatical label transference algorithms are gaining popularity over manual marker gene-based annotation methods. However, most existing methods fail to unify cell-type annotation with dimensionality reduction and are unable to generate deep latent representation from the perspective of data generation.

RESULTS

In this article, we propose scSemiGAN, a single-cell semi-supervised cell-type annotation and dimensionality reduction framework based on a generative adversarial network, to overcome these challenges, modeling scRNA-seq data from the aspect of data generation. Our proposed scSemiGAN is capable of performing deep latent representation learning and cell-type label prediction simultaneously. Through extensive comparison with four state-of-the-art annotation methods on diverse simulated and real scRNA-seq datasets, scSemiGAN achieves competitive or superior performance in multiple downstream tasks including cell-type annotation, latent representation visualization, confounding factor removal and enrichment analysis.

AVAILABILITY AND IMPLEMENTATION

The code and data of scSemiGAN are available on GitHub: https://github.com/rafa-nadal/scSemiGAN.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Single-cell and spatial transcriptomics of the infarcted heart define the dynamic onset of the border zone in response to mechanical destabilization

The border zone (BZ) of the infarcted heart is a geographically complex and biologically enigmatic interface separating poorly perfused infarct zones (IZs) from remote zones (RZs). The cellular and molecular mechanisms of myocardial BZs are not well understood because microdissection inevitably combines them with uncontrolled amounts of RZs and IZs. Here, we use single-cell/nucleus RNA sequencing, spatial transcriptomics and multiplexed RNA fluorescence in situ hybridization to redefine the BZ based on cardiomyocyte transcriptomes. BZ1 (Nppa + Xirp2 -) forms a hundreds-of-micrometer-thick layer of morphologically intact cells adjacent to RZs that are detectable within an hour of injury. Meanwhile, BZ2 (Nppa + Xirp2 +) forms a near-single-cell-thick layer of morphologically distorted cardiomyocytes at the IZ edge that colocalize with matricellular protein-expressing myofibroblasts and express predominantly mechanotransduction genes. Surprisingly, mechanical injury alone is sufficient to induce BZ genes. We propose a 'loss of neighbor' hypothesis to explain how ischemic cell death mechanically destabilizes the BZ to induce its transcriptional response.

A spatiotemporally resolved single-cell atlas of the Plasmodium liver stage

Malaria infection involves an obligatory, yet clinically silent liver stage^1,2. Hepatocytes operate in repeating units termed lobules, exhibiting heterogeneous gene expression patterns along the lobule axis³, but the effects of hepatocyte zonation on parasite development at the molecular level remain unknown. Here we combine single-cell RNA sequencing⁴ and single-molecule transcript imaging⁵ to characterize the host and parasite temporal expression programmes in a zonally controlled manner for the rodent malaria parasite Plasmodium berghei ANKA. We identify differences in parasite gene expression in distinct zones, including potentially co-adaptive programmes related to iron and fatty acid metabolism. We find that parasites develop more rapidly in the pericentral lobule zones and identify a subpopulation of periportally biased hepatocytes that harbour abortive infections, reduced levels of Plasmodium transcripts and parasitophorous vacuole breakdown. These 'abortive hepatocytes', which appear predominantly with high parasite inoculum, upregulate immune recruitment and key signalling programmes. Our study provides a resource for understanding the liver stage of Plasmodium infection at high spatial resolution and highlights the heterogeneous behaviour of both the parasite and the host hepatocyte.

A single-nucleus and spatial transcriptomic atlas of the COVID-19 liver reveals topological, functional, and regenerative organ disruption in patients

The molecular underpinnings of organ dysfunction in acute COVID-19 and its potential long-term sequelae are under intense investigation. To shed light on these in the context of liver function, we performed single-nucleus RNA-seq and spatial transcriptomic profiling of livers from 17 COVID-19 decedents. We identified hepatocytes positive for SARS-CoV-2 RNA with an expression phenotype resembling infected lung epithelial cells. Integrated analysis and comparisons with healthy controls revealed extensive changes in the cellular composition and expression states in COVID-19 liver, reflecting hepatocellular injury, ductular reaction, pathologic vascular expansion, and fibrogenesis. We also observed Kupffer cell proliferation and erythrocyte progenitors for the first time in a human liver single-cell atlas, resembling similar responses in liver injury in mice and in sepsis, respectively. Despite the absence of a clinical acute liver injury phenotype, endothelial cell composition was dramatically impacted in COVID-19, concomitantly with extensive alterations and profibrogenic activation of reactive cholangiocytes and mesenchymal cells. Our atlas provides novel insights into liver physiology and pathology in COVID-19 and forms a foundational resource for its investigation and understanding.

A single-nucleus and spatial transcriptomic atlas of the COVID-19 liver reveals topological, functional, and regenerative organ disruption in patients

The molecular underpinnings of organ dysfunction in acute COVID-19 and its potential long-term sequelae are under intense investigation. To shed light on these in the context of liver function, we performed single-nucleus RNA-seq and spatial transcriptomic profiling of livers from 17 COVID-19 decedents. We identified hepatocytes positive for SARS-CoV-2 RNA with an expression phenotype resembling infected lung epithelial cells. Integrated analysis and comparisons with healthy controls revealed extensive changes in the cellular composition and expression states in COVID-19 liver, reflecting hepatocellular injury, ductular reaction, pathologic vascular expansion, and fibrogenesis. We also observed Kupffer cell proliferation and erythrocyte progenitors for the first time in a human liver single-cell atlas, resembling similar responses in liver injury in mice and in sepsis, respectively. Despite the absence of a clinical acute liver injury phenotype, endothelial cell composition was dramatically impacted in COVID-19, concomitantly with extensive alterations and profibrogenic activation of reactive cholangiocytes and mesenchymal cells. Our atlas provides novel insights into liver physiology and pathology in COVID-19 and forms a foundational resource for its investigation and understanding.

Advances in single-cell sequencing technology in the field of hepatocellular carcinoma

Tumors are a class of diseases characterized by altered genetic information and uncontrolled growth. Sequencing technology provide researchers with a better way to explore specific tumor pathogenesis. In recent years, single-cell sequencing technology has shone in tumor research, especially in the study of liver cancer, revealing phenomena that were unexplored by previous studies. Single-cell sequencing (SCS) is a technique for sequencing the cellular genome, transcriptome, epigenome, proteomics, or metabolomics after dissociation of tissues into single cells. Compared with traditional bulk sequencing, single-cell sequencing can dissect human tumors at single-cell resolution, finely delineate different cell types, and reveal the heterogeneity of tumor cells. In view of the diverse pathological types and complex pathogenesis of hepatocellular carcinoma (HCC), the study of the heterogeneity among tumor cells can help improve its clinical diagnosis, treatment and prognostic judgment. On this basis, SCS has revolutionized our understanding of tumor heterogeneity, tumor immune microenvironment, and clonal evolution of tumor cells. This review summarizes the basic process and development of single-cell sequencing technology and its increasing role in the field of hepatocellular carcinoma.

Transcriptomic Analysis of Cirrhosis-Like Hepatocellular Carcinoma Reveals Distinct Molecular Characteristics and Pathologic Staging Implications

OBJECTIVES

Cirrhosis-like hepatocellular carcinoma (CL-HCC) is a rare hepatocellular malignancy characterized by multiple tumor nodules that clinically, radiologically, macroscopically, and microscopically mimic cirrhosis. We aimed to elucidate the molecular biology of CL-HCC and determine tumor nodule clonality.

METHODS

We performed RNA sequencing on formalin-fixed, paraffin-embedded tissue from confirmed CL-HCC cases (n = 6), along with corresponding nonneoplastic hepatic tissue (n = 4) when available. Transcriptomes from our previous work on steatohepatitic hepatocellular carcinoma and The Cancer Genome Atlas (TCGA) were used for comparison purposes.

RESULTS

Histologically, CL-HCC displayed innumerable nodules and extensive vascular invasion. Intratumoral nodule comparison indicated that the multiple nodules were all clonally related, not independent primaries. The unique histomorphologic appearance corresponded with a distinct transcriptome compared with other HCCs, including fibrolamellar HCC (n = 6), steatohepatitic HCC (n = 8), and conventional HCC in TCGA (n = 424). Tumor-normal gene expression analysis revealed consistent overexpression of several genes involved in degradation of tissue matrix. No recurrent translocations or point mutations were identified. CL-HCC showed a gene expression profile indicative of zone 2 hepatocytes.

CONCLUSIONS

CL-HCC's distinctive clinicopathologic features correspond to a unique gene expression profile, increased expression of invasive markers, features of zone 2 hepatocytes, and features suggestive of intratumoral nodule monoclonality.

Integrative systems analysis identifies genetic and dietary modulators of bile acid homeostasis

Bile acids (BAs) are complex and incompletely understood enterohepatic-derived hormones that control whole-body metabolism. Here, we profiled postprandial BAs in the liver, feces, and plasma of 360 chow- or high-fat-diet-fed BXD male mice and demonstrated that both genetics and diet strongly influence BA abundance, composition, and correlation with metabolic traits. Through an integrated systems approach, we mapped hundreds of quantitative trait loci that modulate BAs and identified both known and unknown regulators of BA homeostasis. In particular, we discovered carboxylesterase 1c (Ces1c) as a genetic determinant of plasma tauroursodeoxycholic acid (TUDCA), a BA species with established disease-preventing actions. The association between Ces1c and plasma TUDCA was validated using data from independent mouse cohorts and a Ces1c knockout mouse model. Collectively, our data are a unique resource to dissect the physiological importance of BAs as determinants of metabolic traits, as underscored by the identification of CES1C as a master regulator of plasma TUDCA levels.

Variant calling enhances the identification of cancer cells in single-cell RNA sequencing data

Single-cell RNA-sequencing is an invaluable research tool that allows for the investigation of gene expression in heterogeneous cancer cell populations in ways that bulk RNA-seq cannot. However, normal (i.e., non tumor) cells in cancer samples have the potential to confound the downstream analysis of single-cell RNA-seq data. Existing methods for identifying cancer and normal cells include copy number variation inference, marker-gene expression analysis, and expression-based clustering. This work aims to extend the existing approaches for identifying cancer cells in single-cell RNA-seq samples by incorporating variant calling and the identification of putative driver alterations. We found that putative driver alterations can be detected in single-cell RNA-seq data obtained with full-length transcript technologies and noticed that a subset of cells in tumor samples are enriched for putative driver alterations as compared to normal cells. Furthermore, we show that the number of putative driver alterations and inferred copy number variation are not correlated in all samples. Taken together, our findings suggest that augmenting existing cancer-cell filtering methods with variant calling and analysis can increase the number of tumor cells that can be confidently included in downstream analyses of single-cell full-length transcript RNA-seq datasets.

The Slc25a47 locus is a novel determinant of hepatic mitochondrial function implicated in liver fibrosis

BACKGROUND & AIMS

Transporters of the SLC25 mitochondrial carrier superfamily bridge cytoplasmic and mitochondrial metabolism by channeling metabolites across mitochondrial membranes and are pivotal for metabolic homeostasis. Despite their physiological relevance as gatekeepers of cellular metabolism, most of the SLC25 family members remain uncharacterized. We undertook a comprehensive tissue distribution analysis of all Slc25 family members across metabolic organs and identified SLC25A47 as a liver-specific mitochondrial carrier.

METHODS

We used a murine loss-of-function model to unravel the role of this transporter in mitochondrial and hepatic homeostasis. We performed extensive metabolic phenotyping and molecular characterization of newly generated Slc25a47^hep-/- and Slc25a47-Fgf21^hep-/- mice.

RESULTS

Slc25a47^hep-/- mice displayed a wide variety of metabolic abnormalities, as a result of sustained energy deficiency in the liver originating from impaired mitochondrial respiration. This mitochondrial phenotype was associated with an activation of the mitochondrial stress response (MSR) in the liver, and the development of fibrosis, which was exacerbated upon feeding a high-fat high-sucrose diet. The MSR induced the secretion of several mitokines, amongst which FGF21 played a preponderant role on systemic physiology. To dissect the FGF21-dependent and -independent physiological changes induced in Slc25a47^hep-/- mice, we generated a double Slc25a47-Fgf21^hep-/- mouse model and demonstrated that several aspects of the hypermetabolic state were driven by hepatic secretion of FGF21. On the other hand, the metabolic fuel inflexibility observed in Slc25a47^hep-/- mice could not be rescued with the genetic removal of Fgf21.

CONCLUSION

Collectively, our data place the Slc25a47 locus at the center of mitochondrial homeostasis, which upon dysfunction triggers robust liver-specific and systemic adaptive stress responses. The prominent role of the Slc25a47 locus in hepatic fibrosis identifies this carrier, or its transported metabolite, as a potential target for therapeutic intervention.

LAY SUMMARY

Herein, we report the importance of a locus containing a liver-specific gene coding for a mitochondrial transport protein called SLC25A47. Mitochondria are the powerhouses of cells. They are crucial for metabolism and energy generation. We show that mice with genetic disruption of the Slc25a47 locus cannot maintain mitochondrial homeostasis (balance), leading to wide-ranging problems in the liver that have far-reaching physiological consequences.

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.

Polygenic enrichment distinguishes disease associations of individual cells in single-cell RNA-seq data

Single-cell RNA sequencing (scRNA-seq) provides unique insights into the pathology and cellular origin of disease. We introduce single-cell disease relevance score (scDRS), an approach that links scRNA-seq with polygenic disease risk at single-cell resolution, independent of annotated cell types. scDRS identifies cells exhibiting excess expression across disease-associated genes implicated by genome-wide association studies (GWASs). We applied scDRS to 74 diseases/traits and 1.3 million single-cell gene-expression profiles across 31 tissues/organs. Cell-type-level results broadly recapitulated known cell-type-disease associations. Individual-cell-level results identified subpopulations of disease-associated cells not captured by existing cell-type labels, including T cell subpopulations associated with inflammatory bowel disease, partially characterized by their effector-like states; neuron subpopulations associated with schizophrenia, partially characterized by their spatial locations; and hepatocyte subpopulations associated with triglyceride levels, partially characterized by their higher ploidy levels. Genes whose expression was correlated with the scDRS score across cells (reflecting coexpression with GWAS disease-associated genes) were strongly enriched for gold-standard drug target and Mendelian disease genes.

Androgen receptor functions in pericentral hepatocytes to decrease gluconeogenesis and avoid hyperglycemia and obesity in male mice

BACKGROUND

Although the impact of hepatic androgen receptor (AR) pathway on liver pathogenesis was documented, its physiological function in normal liver is remained unclear. This study aims to investigate if hepatic AR acts on metabolism, the major liver function, using a hepatic-specific AR-transgenic (H-ARTG) mouse model.

METHODS

We established the albumin promoter driven H-ARTG mice and included wild type (WT) and H-ARKO mice for study. The body weight, specific metabolic parameters and results from various tolerance tests were compared in different groups of mice fed a chow diet, from 2 to 18 months of age. Glucose feeding and insulin treatment were used to study the expression and zonal distribution pattern of AR and related genes in liver at different prandial stages.

RESULTS

The body weight of H-ARTG mice fed a chow diet was 15 % lower than that of wild-type mice, preceded by lower blood glucose and liver triglyceride levels caused by AR reduced hepatic gluconeogenesis. The opposite phenotypes identified in H-ARKO and castrated H-ARTG mice support the critical role of activated AR in decreasing gluconeogenesis and triglyceride levels in liver. Hepatic AR acting by enhancing the expression of cytosolic glycerol-3-phosphate dehydrogenase (cGPDH), a key of glycerophosphate shuttle, was identified as one mechanism to decrease gluconeogenesis from glycerol. We further found AR normally expressed in zone 3 of hepatic lobules. Its level fluctuates dependent on the demand of glucose, decreased by fasting but increased by glucose uptake or insulin stimulation.

CONCLUSION

AR is a newly identified zone 3 hepatic gene with function in reducing blood glucose and body weight in mice. It suggests that stabilization of hepatic AR is a new direction to prevent hyperglycemia, obesity and nonalcoholic fatty liver disease (NAFLD) in males.

Simulation of the crosstalk between glucose and acetaminophen metabolism in a liver zonation model

The liver metabolizes a variety of substances that sometimes interact and regulate each other. The modeling of a single cell or a single metabolic pathway does not represent the complexity of the organ, including metabolic zonation (heterogeneity of functions) along with liver sinusoids. Here, we integrated multiple metabolic pathways into a single numerical liver zonation model, including drug and glucose metabolism. The model simulated the time-course of metabolite concentrations by the combination of dynamic simulation and metabolic flux analysis and successfully reproduced metabolic zonation and localized hepatotoxicity induced by acetaminophen (APAP). Drug metabolism was affected by nutritional status as the glucuronidation reaction rate changed. Moreover, sensitivity analysis suggested that the reported metabolic characteristics of obese adults and healthy infants in glucose metabolism could be associated with the metabolic features of those in drug metabolism. High activities of phosphoenolpyruvate carboxykinase (PEPCK) and glucose-6-phosphate phosphatase in obese adults led to increased APAP oxidation by cytochrome P450 2E1. In contrast, the high activity of glycogen synthase and low activities of PEPCK and glycogen phosphorylase in healthy infants led to low glucuronidation and high sulfation rates of APAP. In summary, this model showed the effects of glucose metabolism on drug metabolism by integrating multiple pathways into a single liver metabolic zonation model.

Liver group 2 innate lymphoid cells regulate blood glucose levels through IL-13 signaling and suppression of gluconeogenesis

The liver stores glycogen and releases glucose into the blood upon increased energy demand. Group 2 innate lymphoid cells (ILC2) in adipose and pancreatic tissues are known for their involvement in glucose homeostasis, but the metabolic contribution of liver ILC2s has not been studied in detail. Here we show that liver ILC2s are directly involved in the regulation of blood glucose levels. Mechanistically, interleukin (IL)-33 treatment induces IL-13 production in liver ILC2s, while directly suppressing gluconeogenesis in a specific Hnf4a/G6pc-high primary hepatocyte cluster via Stat3. These hepatocytes significantly interact with liver ILC2s via IL-13/IL-13 receptor signaling. The results of transcriptional complex analysis and GATA3-ChIP-seq, ATAC-seq, and scRNA-seq trajectory analyses establish a positive regulatory role for the transcription factor GATA3 in IL-13 production by liver ILC2s, while AP-1 family members are shown to suppress IL-13 release. Thus, we identify a regulatory role and molecular mechanism by which liver ILC2s contribute to glucose homeostasis.

Besides hepatocytes, resident immune cells of the liver are also contributing to the body’s energy homeostasis. Here authors show that group 2 innate lymphoid cells interact with a specific set of hepatocytes in suppressing gluconeogenesis and regulate blood glucose levels via Interleukin-13 signalling.

Liver macrophages in health and disease

Single-cell and spatial transcriptomic technologies have revealed an underappreciated heterogeneity of liver macrophages. This has led us to rethink the involvement of macrophages in liver homeostasis and disease. Identification of conserved gene signatures within these cells across species and diseases is enabling the correct identification of specific macrophage subsets and the generation of more specific tools to track and study the functions of these cells. Here, we discuss what is currently known about the definitions of these different macrophage populations, the markers that can be used to identify them, how they are wired within the liver, and their functional specializations in health and disease.

Spatially resolved characterization of tissue metabolic compartments in fasted and high-fat diet livers

Cells adapt their metabolism to physiological stimuli, and metabolic heterogeneity exists between cell types, within tissues, and subcellular compartments. The liver plays an essential role in maintaining whole-body metabolic homeostasis and is structurally defined by metabolic zones. These zones are well-understood on the transcriptomic level, but have not been comprehensively characterized on the metabolomic level. Mass spectrometry imaging (MSI) can be used to map hundreds of metabolites directly from a tissue section, offering an important advance to investigate metabolic heterogeneity in tissues compared to extraction-based metabolomics methods that analyze tissue metabolite profiles in bulk. We established a workflow for the preparation of tissue specimens for matrix-assisted laser desorption/ionization (MALDI) MSI that can be implemented to achieve broad coverage of central carbon, nucleotide, and lipid metabolism pathways. Herein, we used this approach to visualize the effect of nutrient stress and excess on liver metabolism. Our data revealed a highly organized metabolic tissue compartmentalization in livers, which becomes disrupted under high fat diet. Fasting caused changes in the abundance of several metabolites, including increased levels of fatty acids and TCA intermediates while fatty livers had higher levels of purine and pentose phosphate-related metabolites, which generate reducing equivalents to counteract oxidative stress. This spatially conserved approach allowed the visualization of liver metabolic compartmentalization at 30 μm pixel resolution and can be applied more broadly to yield new insights into metabolic heterogeneity in vivo.

Neuregulin 4 suppresses NASH-HCC development by restraining tumor-prone liver microenvironment

The mammalian liver comprises heterogeneous cell types within its tissue microenvironment that undergo pathophysiological reprogramming in disease states, such as non-alcoholic steatohepatitis (NASH). Patients with NASH are at an increased risk for the development of hepatocellular carcinoma (HCC). However, the molecular and cellular nature of liver microenvironment remodeling that links NASH to liver carcinogenesis remains obscure. Here, we show that diet-induced NASH is characterized by the induction of tumor-associated macrophage (TAM)-like macrophages and exhaustion of cytotoxic CD8+ T cells in the liver. The adipocyte-derived endocrine factor Neuregulin 4 (NRG4) serves as a hormonal checkpoint that restrains this pathological reprogramming during NASH. NRG4 deficiency exacerbated the induction of tumor-prone liver immune microenvironment and NASH-related HCC, whereas transgenic NRG4 overexpression elicited protective effects in mice. In a therapeutic setting, recombinant NRG4-Fc fusion protein exhibited remarkable potency in suppressing HCC and prolonged survival in the treated mice. These findings pave the way for therapeutic intervention of liver cancer by targeting the NRG4 hormonal checkpoint.

NOD-like receptor protein 3 activation causes spontaneous inflammation and fibrosis that mimics human NASH

BACKGROUND AND AIMS

The NOD-like receptor protein 3 (NLRP3) inflammasome is a central contributor to human acute and chronic liver disease, yet the molecular and cellular mechanisms by which its activation precipitates injury remain incompletely understood. Here, we present single cell transcriptomic profiling of livers from a global transgenic tamoxifen-inducible constitutively activated Nlrp3A350V mutant mouse, and we investigate the changes in parenchymal and nonparenchymal liver cell gene expression that accompany inflammation and fibrosis.

APPROACH AND RESULTS

Our results demonstrate that NLRP3 activation causes chronic extramedullary myelopoiesis marked by myeloid progenitors that differentiate into proinflammatory neutrophils, monocytes, and monocyte-derived macrophages. We observed prominent neutrophil infiltrates with increased Ly6gHI and Ly6gINT cells exhibiting transcriptomic signatures of granulopoiesis typically found in the bone marrow. This was accompanied by a marked increase in Ly6cHI monocytes differentiating into monocyte-derived macrophages that express transcriptional programs similar to macrophages of NASH models. NLRP3 activation also down-regulated metabolic pathways in hepatocytes and shifted hepatic stellate cells toward an activated profibrotic state based on expression of collagen and extracellular matrix regulatory genes.

CONCLUSIONS

These results define the single cell transcriptomes underlying hepatic inflammation and fibrosis precipitated by NLRP3 activation. Clinically, our data support the notion that NLRP3-induced mechanisms should be explored as therapeutic target in NASH-like inflammation.

Logic of the Temporal Compartmentalization of the Hepatic Metabolic Cycle

The mammalian liver must cope with various metabolic and physiological changes that normally recur every day and result primarily from rest-activity and fasting-feeding cycles. In this article, I present evidence supporting a temporal compartmentalization of rhythmic hepatic metabolic processes into four main clusters: regulation of energy homeostasis, maintenance of information integrity, immune response, and genetic information flow. I further review literatures and discuss how both the circadian and the newly discovered 12-h ultradian clock work together to regulate these four temporally separated processes in mouse liver, which, interestingly, is largely uncoupled from the liver zonation regulation.

Dynamics of hepatocyte-cholangiocyte cell-fate decisions during liver development and regeneration

The immense regenerative potential of the liver is attributed to the ability of its two key cell types – hepatocytes and cholangiocytes – to trans-differentiate to one another either directly or through intermediate progenitor states. However, the dynamic features of decision-making between these cell-fates during liver development and regeneration remains elusive. Here, we identify a core gene regulatory network comprising c/EBPα, TGFBR2, and SOX9 which is multistable in nature, enabling three distinct cell states – hepatocytes, cholangiocytes, and liver progenitor cells (hepatoblasts/oval cells) – and stochastic switching among them. Predicted expression signature for these three states are validated through multiple bulk and single-cell transcriptomic datasets collected across developmental stages and injury-induced liver repair. This network can also explain the experimentally observed spatial organization of phenotypes in liver parenchyma and predict strategies for efficient cellular reprogramming. Our analysis elucidates how the emergent dynamics of underlying regulatory networks drive diverse cell-fate decisions in liver development and regeneration.

•

Identified minimal regulatory network to model liver development and regeneration

•

Changes in phenotypic landscapes by in-silico perturbations of regulatory networks

•

Ability to explain physiological spatial patterning of liver cell types

•

Decoded strategies for efficient reprogramming among liver cell phenotypes

Young mice administered adult doses of AAV5-hFVIII-SQ achieve therapeutic factor VIII expression into adulthood

Valoctocogene roxaparvovec (AAV5-hFVIII-SQ) gene transfer provided reduced bleeding for adult clinical trial participants with severe hemophilia A. However, pediatric outcomes are unknown. Using a mouse model of hemophilia A, we investigated the effect of vector dose and age at treatment on transgene production and persistence. We dosed AAV5-hFVIII-SQ to neonatal and adult mice based on body weight or at a fixed dose and assessed human factor VIII-SQ variant (hFVIII-SQ) expression through 16 weeks. AAV5-hFVIII-SQ dosed per body weight in neonatal mice did not result in meaningful plasma hFVIII-SQ protein levels in adulthood. When treated with the same total vector genomes per mouse as adult mice, neonates maintained hFVIII-SQ expression into adulthood, although plasma levels were 3- to 4-fold lower versus mice dosed as adults. Mice <1 week old initially exhibited high hFVIII-SQ plasma levels and maintained meaningful levels into adulthood, despite a partial decline potentially due to age-related body mass and blood volume increases. Spatial transduction patterns differed between mice dosed as neonates versus adults. No features of hepatotoxicity or endoplasmic reticulum stress were observed with dosing at any age. These data suggest that young mice require the same total vector genomes as adult mice to sustain hFVIII-SQ plasma levels.

Distinct phenotypic spectra of hepatocellular carcinoma in liver-specific tumor suppressor-deficient hepatitis B virus transgenic mice

The hepatitis B virus (HBV) transgenic mouse model was used to interrogate the origins of HCC heterogeneity. HBV biosynthesis was used as a marker of liver tumor heterogeneity. Principal component and correlation analysis of HBV and cellular transcript levels demonstrated major differences within and between the gene expression profiles of Apc-deficient, Apc-deficient Pten-deficient, and Pten-deficient HCC. Hence, both oncogenic stimuli and zonal hepatocyte properties determine heterogeneous HCC phenotypes. Additionally, Apc-deficient HCC display decreased expression of Apob, Otc and Tet2 relative to Pten-deficient HCC and control liver tissue suggesting their gene products may represent markers of Apc-deficient HCC. A subset of human HCC with mutations in the β-catenin gene (CTNNB1) displayed a gene expression profile similar to that observed in the mouse Apc-deficient HCC indicating this model of liver cancer may be useful for interrogating the molecular properties of these tumors and their potential therapeutic vulnerabilities.

A strategy of vascular-targeted therapy for liver fibrosis

BACKGROUND AND AIMS

No effective treatments are available for liver fibrosis. Angiogenesis is deeply involved in liver fibrogenesis. However, current controversial results suggest it is difficult to treat liver fibrosis through vascular targeting. There are three different microvessels in liver: portal vessels, liver sinusoids, and central vessels. The changes and roles for each of the three different vessels during liver fibrogenesis are unclear. We propose that they play different roles during liver fibrogenesis, and a single vascular endothelial cell (EC) regulator is not enough to fully regulate these three vessels to treat liver fibrosis. Therefore, a combined regulation of multiple different EC regulatory signaling pathway may provide new strategies for the liver fibrosis therapy. Herein, we present a proof-of-concept strategy by combining the regulation of leukocyte cell-derived chemotaxin 2 (LECT2)/tyrosine kinase with immunoglobulin-like and epidermal growth factor-like domains 1 signaling with that of vascular endothelial growth factor (VEGF)/recombinant VEGF (rVEGF) signaling.

APPROACH AND RESULTS

The CCl₄ -induced mouse liver fibrosis model and NASH model were both used. During fibrogenesis, vascular changes occurred at very early stage, and different liver vessels showed different changes and played different roles: decreased portal vessels, increased sinusoid capillarization and the increased central vessels the increase of portal vessels alleviates liver fibrosis, the increase of central vessels aggravates liver fibrosis, and the increase of sinusoid capillarization aggravates liver fibrosis. The combinational treatment of adeno-associated viral vector serotype 9 (AAV9)-LECT2-short hairpin RNA (shRNA) and rVEGF showed improved therapeutic effects, but it led to serious side effects. The combination of AAV9-LECT2-shRNA and bevacizumab showed both improved therapeutic effects and decreased side effects.

CONCLUSIONS

Liver vascular changes occurred at very early stage of fibrogenesis. Different vessels play different roles in liver fibrosis. The combinational treatment of AAV9-LECT2-shRNA and bevacizumab could significantly improve the therapeutic effects on liver fibrosis.

Accurate inference of genome-wide spatial expression with iSpatial

Spatially resolved transcriptomic analyses can reveal molecular insights underlying tissue structure and context-dependent cell-cell or cell-environment interaction. Because of the current technical limitation, obtaining genome-wide spatial transcriptome at single-cell resolution is challenging. Here, we developed a new algorithm named iSpatial to derive the spatial pattern of the entire transcriptome by integrating spatial transcriptomic and single-cell RNA-seq datasets. Compared to other existing methods, iSpatial has higher accuracy in predicting gene expression and spatial distribution. Furthermore, it reduces false-positive and false-negative signals in the original datasets. By testing iSpatial with multiple spatial transcriptomic datasets, we demonstrate its wide applicability to datasets from different tissues and by different techniques. Thus, we provide a computational approach to reveal spatial organization of the entire transcriptome at single-cell resolution. With numerous high-quality datasets available in the public domain, iSpatial provides a unique way to understand the structure and function of complex tissues and disease processes.