Epigenomic map of human liver reveals principles of zonated morphogenic and metabolic control

Cell atlas of the regenerating human liver after portal vein embolization

The liver has the remarkable capacity to regenerate. In the clinic, regeneration is induced by portal vein embolization, which redirects portal blood flow, resulting in liver hypertrophy in locations with increased blood supply, and atrophy of embolized segments. Here, we apply single-cell and single-nucleus transcriptomics on healthy, hypertrophied, and atrophied patient-derived liver samples to explore cell states in the regenerating liver. Our data unveils pervasive upregulation of genes associated with developmental processes, cellular adhesion, and inflammation in post-portal vein embolization liver, disrupted portal-central hepatocyte zonation, and altered cell subtype composition of endothelial and immune cells. Interlineage crosstalk analysis reveals mesenchymal cells as an interaction hub between immune and endothelial cells, and highlights the importance of extracellular matrix proteins in liver regeneration. Moreover, we establish tissue-scale iterative indirect immunofluorescence imaging for high-dimensional spatial analysis of perivascular microenvironments, uncovering changes to tissue architecture in regenerating liver lobules. Altogether, our data is a rich resource revealing cellular and histological changes in human liver regeneration.

A spatiotemporal atlas of mouse liver homeostasis and regeneration

The mechanism by which mammalian liver cell responses are coordinated during tissue homeostasis and perturbation is poorly understood, representing a major obstacle in our understanding of many diseases. This knowledge gap is caused by the difficulty involved with studying multiple cell types in different states and locations, particularly when these are transient. We have combined Stereo-seq (spatiotemporal enhanced resolution omics-sequencing) with single-cell transcriptomic profiling of 473,290 cells to generate a high-definition spatiotemporal atlas of mouse liver homeostasis and regeneration at the whole-lobe scale. Our integrative study dissects in detail the molecular gradients controlling liver cell function, systematically defining how gene networks are dynamically modulated through intercellular communication to promote regeneration. Among other important regulators, we identified the transcriptional cofactor TBL1XR1 as a rheostat linking inflammation to Wnt/β-catenin signaling for facilitating hepatocyte proliferation. Our data and analytical pipelines lay the foundation for future high-definition tissue-scale atlases of organ physiology and malfunction.

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.

Mitochondrial heterogeneity and adaptations to cellular needs

Although it is well described that mitochondria are at the epicentre of the energy demands of a cell, it is becoming important to consider how each cell tailors its mitochondrial composition and functions to suit its particular needs beyond ATP production. Here we provide insight into mitochondrial heterogeneity throughout development as well as in tissues with specific energy demands and discuss how mitochondrial malleability contributes to cell fate determination and tissue remodelling.

Emerging Roles of Spatial Transcriptomics in Liver Research

Spatial transcriptomics, leveraging sequencing- and imaging-based techniques, has emerged as a groundbreaking technology for mapping gene expression within the complex architectures of tissues. This approach provides an in-depth understanding of cellular and molecular dynamics across various states of healthy and diseased livers. Through the integration of sophisticated bioinformatics strategies, it enables detailed exploration of cellular heterogeneity, transitions in cell states, and intricate cell-cell interactions with remarkable precision. In liver research, spatial transcriptomics has been particularly revelatory, identifying distinct zonated functions of hepatocytes that are crucial for understanding the metabolic and detoxification processes of the liver. Moreover, this technology has unveiled new insights into the pathogenesis of liver diseases, such as the role of lipid-associated macrophages in steatosis and endothelial cell signals in liver regeneration and repair. In the domain of liver cancer, spatial transcriptomics has proven instrumental in delineating intratumor heterogeneity, identifying supportive microenvironmental niches and revealing the complex interplay between tumor cells and the immune system as well as susceptibility to immune checkpoint inhibitors. In conclusion, spatial transcriptomics represents a significant advance in hepatology, promising to enhance our understanding and treatment of liver diseases.

Histone proteoform analysis reveals epigenetic changes in adult mouse brown adipose tissue in response to cold stress

BACKGROUND

Regulation of the thermogenic response by brown adipose tissue (BAT) is an important component of energy homeostasis with implications for the treatment of obesity and diabetes. Our preliminary analyses of RNA-Seq data uncovered many nodes representing epigenetic modifiers that are altered in BAT in response to chronic thermogenic activation. Thus, we hypothesized that chronic thermogenic activation broadly alters epigenetic modifications of DNA and histones in BAT.

RESULTS

Motivated to understand how BAT function is regulated epigenetically, we developed a novel method for the first-ever unbiased top-down proteomic quantitation of histone modifications in BAT and validated our results with a multi-omic approach. To test our hypothesis, wildtype male C57BL/6J mice were housed under chronic conditions of thermoneutral temperature (TN, 28°C), mild cold/room temperature (RT, 22°C), or severe cold (SC, 8°C) and BAT was analyzed for DNA methylation and histone modifications. Methylation of promoters and intragenic regions in genomic DNA decrease in response to chronic cold exposure. Integration of DNA methylation and RNA expression datasets suggest a role for epigenetic modification of DNA in regulation of gene expression in response to cold. In response to cold housing, we observe increased bulk acetylation of histones H3.2 and H4, increased histone H3.2 proteoforms with di- and trimethylation of lysine 9 (K9me2 and K9me3), and increased histone H4 proteoforms with acetylation of lysine 16 (K16ac) in BAT.

CONCLUSIONS

Our results reveal global epigenetically-regulated transcriptional "on" and "off" signals in murine BAT in response to varying degrees of chronic cold stimuli and establish a novel methodology to quantitatively study histones in BAT, allowing for direct comparisons to decipher mechanistic changes during the thermogenic response. Additionally, we make histone PTM and proteoform quantitation, RNA splicing, RRBS, and transcriptional footprint datasets available as a resource for future research.

Dynamic, IPSC-derived hepatic tissue tri-culture system for the evaluation of liver physiologyin vitro

Availability of hepatic tissue for the investigation of metabolic processes is severely limited. While primary hepatocytes or animal models are widely used in pharmacological applications, a change in methodology towards more sustainable and ethical assays is highly desirable. Stem cell derived hepatic cells are generally regarded as a viable alternative for the above model systems, if current limitations in functionality and maturation can be overcome. By combining microfluidic organ-on-a-chip technology with individually differentiated, multicellular hepatic tissue fractions, we aim to improve overall functionality of hepatocyte-like cells, as well as evaluate cellular composition and interactions with non-parenchymal cell populations towards the formation of mature liver tissue. Utilizing a multi-omic approach, we show the improved maturation profiles of hepatocyte-like cells maintained in a dynamic microenvironment compared to standard tissue culture setups without continuous perfusion. In order to evaluate the resulting tissue, we employ single cell sequencing to distinguish formed subpopulations and spatial localization. While cellular input was strictly defined based on established differentiation protocols of parenchyma, endothelial and stellate cell fractions, resulting hepatic tissue was shown to comprise a complex mixture of epithelial and non-parenchymal fractions with specific local enrichment of phenotypes along the microchannel. Following this approach, we show the importance of passive, paracrine developmental processes in tissue formation. Using such complex tissue models is a crucial first step to develop stem cell-derivedin vitrosystems that can compare functionally with currently used pharmacological and toxicological applications.

Alterations in zonal distribution and plasma membrane localization of hepatocyte bile acid transporters in patients with NAFLD

BACKGROUND

NAFLD is highly prevalent with limited treatment options. Bile acids (BAs) increase in the systemic circulation and liver during NAFLD progression. Changes in plasma membrane localization and zonal distribution of BA transporters can influence transport function and BA homeostasis. However, a thorough characterization of how NAFLD influences these factors is currently lacking. This study aimed to evaluate the impact of NAFLD and the accompanying histologic features on the functional capacity of key hepatocyte BA transporters across zonal regions in human liver biopsies.

METHODS

A novel machine learning image classification approach was used to quantify relative zonal abundance and plasma membrane localization of BA transporters (bile salt export pump [BSEP], sodium-taurocholate cotransporting polypeptide, organic anion transporting polypeptide [OATP] 1B1 and OATP1B3) in non-diseased (n = 10), NAFL (n = 9), and NASH (n = 11) liver biopsies. Based on these data, membrane-localized zonal abundance (MZA) measures were developed to estimate transporter functional capacity.

RESULTS

NAFLD diagnosis and histologic scoring were associated with changes in transporter membrane localization and zonation. Increased periportal BSEPMZA (mean proportional difference compared to non-diseased liver of 0.090) and decreased pericentral BSEPMZA (-0.065) were observed with NASH and also in biopsies with higher histologic scores. Compared to Non-diseased Liver, periportal OATP1B3MZA was increased in NAFL (0.041) and NASH (0.047). Grade 2 steatosis (mean proportional difference of 0.043 when compared to grade 0) and grade 1 lobular inflammation (0.043) were associated with increased periportal OATP1B3MZA.

CONCLUSIONS

These findings provide novel mechanistic insight into specific transporter alterations that impact BA homeostasis in NAFLD. Changes in BSEPMZA likely contribute to altered BA disposition and pericentral microcholestasis previously reported in some patients with NAFLD. BSEPMZA assessment could inform future development and optimization of NASH-related pharmacotherapies.

A spatial map of hepatic mitochondria uncovers functional heterogeneity shaped by nutrient-sensing signaling

In the liver, mitochondria are exposed to different concentrations of nutrients due to their spatial positioning across the periportal and pericentral axis. How the mitochondria sense and integrate these signals to respond and maintain homeostasis is not known. Here, we combine intravital microscopy, spatial proteomics, and functional assessment to investigate mitochondrial heterogeneity in the context of liver zonation. We find that periportal and pericentral mitochondria are morphologically and functionally distinct; beta-oxidation is elevated in periportal regions, while lipid synthesis is predominant in the pericentral mitochondria. In addition, comparative phosphoproteomics reveals spatially distinct patterns of mitochondrial composition and potential regulation via phosphorylation. Acute pharmacological modulation of nutrient sensing through AMPK and mTOR shifts mitochondrial phenotypes in the periportal and pericentral regions, linking nutrient gradients across the lobule and mitochondrial heterogeneity. This study highlights the role of protein phosphorylation in mitochondrial structure, function, and overall homeostasis in hepatic metabolic zonation. These findings have important implications for liver physiology and disease.

Single-cell spatial multi-omics and deep learning dissect enhancer-driven gene regulatory networks in liver zonation

In the mammalian liver, hepatocytes exhibit diverse metabolic and functional profiles based on their location within the liver lobule. However, it is unclear whether this spatial variation, called zonation, is governed by a well-defined gene regulatory code. Here, using a combination of single-cell multiomics, spatial omics, massively parallel reporter assays and deep learning, we mapped enhancer-gene regulatory networks across mouse liver cell types. We found that zonation affects gene expression and chromatin accessibility in hepatocytes, among other cell types. These states are driven by the repressors TCF7L1 and TBX3, alongside other core hepatocyte transcription factors, such as HNF4A, CEBPA, FOXA1 and ONECUT1. To examine the architecture of the enhancers driving these cell states, we trained a hierarchical deep learning model called DeepLiver. Our study provides a multimodal understanding of the regulatory code underlying hepatocyte identity and their zonation state that can be used to engineer enhancers with specific activity levels and zonation patterns.

Epigenetic-based antioxidant effect of an ethanolic extract of Corylus avellana L. on THLE-2 human primary hepatocytes

The ethanolic extract of Corylus avellana L hazelnut, prepared in our laboratories, has been previously characterized by liquid chromatography coupled to high resolution mass spectrometry. We here aimed at testing the antioxidant effect of such extract in H2O2-challenged THLE-2 human primary hepatocytes and verified whether it might be based on microRNA-34b/c expression changes. We here demonstrate that miR-34b/miR-34c undergo significant stimulation (≥2-fold change, p < 0.05) in THLE-2 when treated for 72h with not-toxic hazelnut concentrations (0.04-0.4 mg/ml), when compared with 0.06% ethanol control. When administered with H2O2 (1000-2000 µM, 24h), THLE-2 are significantly protected from oxidative stress if pre-treated with hazelnut, the H2O2-driven cytotoxicity and reactive oxygen species generation being recovered by hazelnut extract, through miR-34b/c stimulation. Although preliminary, our findings pave the way for further preclinical studies aimed at validating the possible health-related application of hazelnut matrix, and/or its metabolites, as powerful epigenetic-based drugs, food supplements or nutraceuticals.

A spatial map of hepatic mitochondria uncovers functional heterogeneity shaped by nutrient-sensing signaling

In the liver, mitochondria are exposed to different concentrations of nutrients due to their spatial positioning across the periportal (PP) and pericentral (PC) axis. How these mitochondria sense and integrate these signals to respond and maintain homeostasis is not known. Here, we combined intravital microscopy, spatial proteomics, and functional assessment to investigate mitochondrial heterogeneity in the context of liver zonation. We found that PP and PC mitochondria are morphologically and functionally distinct; beta-oxidation was elevated in PP regions, while lipid synthesis was predominant in the PC mitochondria. In addition, comparative phosphoproteomics revealed spatially distinct patterns of mitochondrial composition and potential regulation via phosphorylation. Acute pharmacological modulation of nutrient sensing through AMPK and mTOR shifted mitochondrial phenotypes in the PP and PC regions, linking nutrient gradients across the lobule and mitochondrial heterogeneity. This study highlights the role of protein phosphorylation in mitochondrial structure, function, and overall homeostasis in hepatic metabolic zonation. These findings have important implications for liver physiology and disease.

Scalable production of tissue-like vascularized liver organoids from human PSCs

The lack of physiological parity between 2D cell culture and in vivo culture has led to the development of more organotypic models, such as organoids. Organoid models have been developed for a number of tissues, including the liver. Current organoid protocols are characterized by a reliance on extracellular matrices (ECMs), patterning in 2D culture, costly growth factors and a lack of cellular diversity, structure, and organization. Current hepatic organoid models are generally simplistic and composed of hepatocytes or cholangiocytes, rendering them less physiologically relevant compared to native tissue. We have developed an approach that does not require 2D patterning, is ECM independent, and employs small molecules to mimic embryonic liver development that produces large quantities of liver-like organoids. Using single-cell RNA sequencing and immunofluorescence, we demonstrate a liver-like cellular repertoire, a higher order cellular complexity, presenting with vascular luminal structures, and a population of resident macrophages: Kupffer cells. The organoids exhibit key liver functions, including drug metabolism, serum protein production, urea synthesis and coagulation factor production, with preserved post-translational modifications such as N-glycosylation and functionality. The organoids can be transplanted and maintained long term in mice producing human albumin. The organoids exhibit a complex cellular repertoire reflective of the organ and have de novo vascularization and liver-like function. These characteristics are a prerequisite for many applications from cellular therapy, tissue engineering, drug toxicity assessment, and disease modeling to basic developmental biology.

Human liver single nuclear RNA sequencing implicates BMPR2, GDF15, arginine, and estrogen in portopulmonary hypertension

Portopulmonary hypertension (PoPH) is a type of pulmonary vascular disease due to portal hypertension that exhibits high morbidity and mortality. The mechanisms driving disease are unknown, and transcriptional characteristics unique to the PoPH liver remain unexplored. Here, we apply single nuclear RNA sequencing to compare cirrhotic livers from patients with and without PoPH. We identify characteristics unique to PoPH in cells surrounding the central hepatic vein, including increased growth differentiation factor signaling, enrichment of the arginine biosynthesis pathway, and differential expression of the bone morphogenic protein type II receptor and estrogen receptor type I genes. These results provide insight into the transcriptomic characteristics of the PoPH liver and mechanisms by which PoPH cellular dysfunction might contribute to pulmonary vascular remodeling.

Multivariate Evaluation of DNA Quality Differences in Different Preanalytical Procedures in Mouse Livers

Successful histogenetic research relies on proper handling of tissue samples to maximize DNA quality. As the largest gland in the body, the liver is particularly sensitive to sample mishandling owing to its enzymatic and transcriptional activity. However, the impact of preanalytical procedures on the quality of extracted liver DNA remains poorly understood. In this study, we assessed the impact of extraction methods, duration of ex vivo liver ischemia, liver storage time, and temperature on extracted DNA quality. Comprehensive parameters such as DNA yields, purity, DNA integrity number, the percentage of double-stranded DNA (%dsDNA), and PCR amplification of the GAPDH gene fragment were assessed to identify the quality of extracted DNA. Our results revealed that these preanalytical processes had little effect on DIN values and PCR efficiency of GAPDH gene fragments for each sample, whereas the DNA yields, purity, and %dsDNAs varied widely across different processes. For liver DNA extraction, RNase is necessary to isolate "pure" DNA, and the presence of RNase could significantly increase the %dsDNA. In addition, significant increases in the yields, purity, and %dsDNA of extracted DNA were observed in the TissueLyser-processed livers compared with the mortar and pestle or shear cell disruption methods. Further investigation revealed that livers experiencing longer periods of ex vivo ischemia resulted in significantly compromised DNA yields, and to obtain sufficient DNA, the ex vivo liver ischemia should be limited to within 30 minutes. Moreover, compared with storage of livers at -80°C, storage of livers in the vapor phase of liquid nitrogen yielded a higher quality of the extracted DNA. Our findings exhibited significant implications for liver-derived DNA quality assessment and management.

Distinct hepatocyte identities in liver homeostasis and regeneration

The process of metabolic liver zonation is spontaneously established by assigning distributed tasks to hepatocytes along the porto-central blood flow. Hepatocytes fulfil critical metabolic functions, while also maintaining hepatocyte mass by replication when needed. Recent technological advances have enabled us to fine-tune our understanding of hepatocyte identity during homeostasis and regeneration. Subsets of hepatocytes have been identified to be more regenerative and some have even been proposed to function like stem cells, challenging the long-standing view that all hepatocytes are similarly capable of regeneration. The latest data show that hepatocyte renewal during homeostasis and regeneration after liver injury is not limited to rare hepatocytes; however, hepatocytes are not exactly the same. Herein, we review the known differences that give individual hepatocytes distinct identities, recent findings demonstrating how these distinct identities correspond to differences in hepatocyte regenerative capacity, and how the plasticity of hepatocyte identity allows for division of labour among hepatocytes. We further discuss how these distinct hepatocyte identities may play a role during liver disease.

Human induced pluripotent stem cell-derived liver-on-a-chip for studying drug metabolism: the challenge of the cytochrome P450 family

The liver is the primary organ responsible for the detoxification and metabolism of drugs. To date, a lack of preclinical models that accurately emulate drug metabolism by the human liver presents a significant challenge in the drug development pipeline, particularly for predicting drug efficacy and toxicity. In recent years, emerging microfluidic-based organ-on-a-chip (OoC) technologies, combined with human induced pluripotent stem cell (hiPSC) technology, present a promising avenue for the complete recapitulation of human organ biology in a patient-specific manner. However, hiPSC-derived organoids and liver-on-a-chip models have so far failed to sufficiently express cytochrome P450 monooxygenase (CYP450) enzymes, the key enzymes involved in first-pass metabolism, which limits the effectiveness and translatability of these models in drug metabolism studies. This review explores the potential of innovative organoid and OoC technologies for studying drug metabolism and discusses their existing drawbacks, such as low expression of CYP450 genes. Finally, we postulate potential approaches for enhancing CYP450 expression in the hope of paving the way toward developing novel, fully representative liver drug-metabolism models.

Sample size calculation for a NanoString GeoMx spatial transcriptomics experiment to study predictors of fibrosis progression in non-alcoholic fatty liver disease

Sample size calculation for spatial transcriptomics is a novel and understudied research topic. Prior publications focused on powering spatial transcriptomics studies to detect specific cell populations or spatially variable expression patterns on tissue slides. However, power calculations for translational or clinical studies often relate to the difference between patient groups, and this is poorly described in the literature. Here, we present a stepwise process for sample size calculation to identify predictors of fibrosis progression in non-alcoholic fatty liver disease as a case study. We illustrate how to infer study hypothesis from prior bulk RNA-sequencing data, gather input requirements and perform a simulation study to estimate required sample size to evaluate gene expression differences between patients with stable fibrosis and fibrosis progressors with NanoString GeoMx Whole Transcriptome Atlas assay.

Intra-pituitary follicle-stimulating hormone signaling regulates hepatic lipid metabolism in mice

Inter-organ communication is a major hallmark of health and is often orchestrated by hormones released by the anterior pituitary gland. Pituitary gonadotropes secrete follicle-stimulating hormone (FSH) and luteinizing hormone (LH) to regulate gonadal function and control fertility. Whether FSH and LH also act on organs other than the gonads is debated. Here, we find that gonadotrope depletion in adult female mice triggers profound hypogonadism, obesity, glucose intolerance, fatty liver, and bone loss. The absence of sex steroids precipitates these phenotypes, with the notable exception of fatty liver, which results from ovary-independent actions of FSH. We uncover paracrine FSH action on pituitary corticotropes as a mechanism to restrain the production of corticosterone and prevent hepatic steatosis. Our data demonstrate that functional communication of two distinct hormone-secreting cell populations in the pituitary regulates hepatic lipid metabolism.

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.

Impact of Liver Inflammation on Bile Acid Side Chain Shortening and Amidation

Bile acid (BA) synthesis from cholesterol by hepatocytes is inhibited by inflammatory cytokines. Whether liver inflammation also affects BA side chain shortening and conjugation was investigated. In human liver cell lines (IHH, HepG2, and HepaRG), agonists of nuclear receptors including the farnesoid X receptor (FXR), liver X receptor (LXR), and peroxisome proliferator-activated receptors (PPARs) did not affect the expression of BA-related peroxisomal enzymes. In contrast, hepatocyte nuclear factor 4α (HNF4α) inhibition down-regulated acyl-CoA oxidase 2 (ACOX2). ACOX2 was repressed by fibroblast growth factor 19 (FGF19), which was prevented by extracellular signal-regulated kinase (ERK) pathway inhibition. These changes were paralleled by altered BA synthesis (HPLC-MS/MS). Cytokines able to down-regulate cholesterol-7α-hydroxylase (CYP7A1) had little effect on peroxisomal enzymes involved in BA synthesis except for ACOX2 and bile acid-CoA:amino acid N-acyltransferase (BAAT), which were down-regulated, mainly by oncostatin M (OSM). This effect was prevented by Janus kinase (JAK) inhibition, which restored BA side chain shortening and conjugation. The binding of OSM to the extracellular matrix accounted for a persistent effect after culture medium replacement. In silico analysis of four databases (n = 201) and a validation cohort (n = 90) revealed an inverse relationship between liver inflammation and ACOX2/BAAT expression which was associated with changes in HNF4α levels. In conclusion, BA side chain shortening and conjugation are inhibited by inflammatory effectors. However, other mechanisms involved in BA homeostasis counterbalance any significant impact on the serum BA profile.

Computational solutions for spatial transcriptomics

Transcriptome level expression data connected to the spatial organization of the cells and molecules would allow a comprehensive understanding of how gene expression is connected to the structure and function in the biological systems. The spatial transcriptomics platforms may soon provide such information. However, the current platforms still lack spatial resolution, capture only a fraction of the transcriptome heterogeneity, or lack the throughput for large scale studies. The strengths and weaknesses in current ST platforms and computational solutions need to be taken into account when planning spatial transcriptomics studies. The basis of the computational ST analysis is the solutions developed for single-cell RNA-sequencing data, with advancements taking into account the spatial connectedness of the transcriptomes. The scRNA-seq tools are modified for spatial transcriptomics or new solutions like deep learning-based joint analysis of expression, spatial, and image data are developed to extract biological information in the spatially resolved transcriptomes. The computational ST analysis can reveal remarkable biological insights into spatial patterns of gene expression, cell signaling, and cell type variations in connection with cell type-specific signaling and organization in complex tissues. This review covers the topics that help choosing the platform and computational solutions for spatial transcriptomics research. We focus on the currently available ST methods and platforms and their strengths and limitations. Of the computational solutions, we provide an overview of the analysis steps and tools used in the ST data analysis. The compatibility with the data types and the tools provided by the current ST analysis frameworks are summarized.

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.

Single-cell atlases: shared and tissue-specific cell types across human organs

The development of single-cell and spatial transcriptomics methods was instrumental in the conception of the Human Cell Atlas initiative, which aims to generate an integrated map of all cells across the human body. These technology advances are bringing increasing depth and resolution to maps of human organs and tissues, as well as our understanding of individual human cell types. Commonalities as well as tissue-specific features of primary and supportive cell types across human organs are beginning to emerge from these human tissue maps. In this Review, we highlight key biological insights obtained from cross-tissue studies into epithelial, fibroblast, vascular and immune cells based on single-cell gene expression data in humans and contrast it with mechanisms reported in mice.

Significance of Liver Zonation in Hepatocellular Carcinoma

Liver zonation is fundamental to normal liver function, and numerous studies have investigated the microstructure of normal liver lobules. However, only a few studies have explored the zonation signature in hepatocellular carcinoma (HCC). In this study, we investigated the significance of liver zonation in HCC with the help of single-cell RNA sequencing (scRNA-seq) and multicolor immunofluorescence staining. Liver zonation-related genes were extracted from the literature, and a three-gene model was established for HCC prognosis. The model reliability was validated using bulk RNA and single-cell RNA-level data, and the underlying biological mechanism was revealed by a functional enrichment analysis. The results showed that the signaling pathways of high-risk groups were similar to those of perivenous zones in the normal liver, indicating the possible regulating role of hypoxia in HCC zonation. Furthermore, the co-staining results showed that the low-grade tumors lost their zonation features whereas the high-grade tumors lost the expression of zonation-related genes, which supported the results obtained from the sequencing data.

Single cell- and spatial 'Omics revolutionize physiology

Single cell multi- 'Omics and Spatial Transcriptomics are prominent technological highlights of recent years, and both fields still witness a ceaseless firework of novel approaches for high resolution profiling of additional omics layers. As all life processes in organs and organisms are based on the functions of their fundamental building blocks, the individual cells and their interactions, these methods are of utmost worth for the study of physiology in health and disease. Recent discoveries on embryonic development, tumor immunology, the detailed cellular composition and function of complex tissues like for example the kidney or the brain, different roles of the same cell type in different organs, the oncogenic program of individual tumor entities, or the architecture of immunopathology in infected tissue are based on single cell and spatial transcriptomics experiments. In this review, we will give a broad overview of technological concepts for single cell and spatial analysis, showing both advantages and limitations, and illustrate their impact with some particularly impressive case studies.

Cell-Type Resolved Insights into the Cis-Regulatory Genome of NAFLD

The prevalence of non-alcoholic fatty liver disease (NAFLD) is increasing rapidly, and unmet treatment can result in the development of hepatitis, fibrosis, and liver failure. There are difficulties involved in diagnosing NAFLD early and for this reason there are challenges involved in its treatment. Furthermore, no drugs are currently approved to alleviate complications, a fact which highlights the need for further insight into disease mechanisms. NAFLD pathogenesis is associated with complex cellular changes, including hepatocyte steatosis, immune cell infiltration, endothelial dysfunction, hepatic stellate cell activation, and epithelial ductular reaction. Many of these cellular changes are controlled by dramatic changes in gene expression orchestrated by the cis-regulatory genome and associated transcription factors. Thus, to understand disease mechanisms, we need extensive insights into the gene regulatory mechanisms associated with tissue remodeling. Mapping cis-regulatory regions genome-wide is a step towards this objective and several current and emerging technologies allow detection of accessible chromatin and specific histone modifications in enriched cell populations of the liver, as well as in single cells. Here, we discuss recent insights into the cis-regulatory genome in NAFLD both at the organ-level and in specific cell populations of the liver. Moreover, we highlight emerging technologies that enable single-cell resolved analysis of the cis-regulatory genome of the liver.

Selective PPARα Modulator Pemafibrate and Sodium-Glucose Cotransporter 2 Inhibitor Tofogliflozin Combination Treatment Improved Histopathology in Experimental Mice Model of Non-Alcoholic Steatohepatitis

Ballooning degeneration of hepatocytes is a major distinguishing histological feature of non-alcoholic steatosis (NASH) progression that can lead to cirrhosis and hepatocellular carcinoma (HCC). In this study, we evaluated the effect of the selective PPARα modulator (SPPARMα) pemafibrate (Pema) and sodium-glucose cotransporter 2 (SGLT2) inhibitor tofogliflozin (Tofo) combination treatment on pathological progression in the liver of a mouse model of NASH (STAM) at two time points (onset of NASH progression and HCC survival). At both time points, the Pema and Tofo combination treatment significantly alleviated hyperglycemia and hypertriglyceridemia. The combination treatment significantly reduced ballooning degeneration of hepatocytes. RNA-seq analysis suggested that Pema and Tofo combination treatment resulted in an increase in glyceroneogenesis, triglyceride (TG) uptake, lipolysis and liberated fatty acids re-esterification into TG, lipid droplet (LD) formation, and Cidea/Cidec ratio along with an increased number and reduced size and area of LDs. In addition, combination treatment reduced expression levels of endoplasmic reticulum stress-related genes (Ire1a, Grp78, Xbp1, and Phlda3). Pema and Tofo treatment significantly improved survival rates and reduced the number of tumors in the liver compared to the NASH control group. These results suggest that SPPARMα and SGLT2 inhibitor combination therapy has therapeutic potential to prevent NASH-HCC progression.

Zonation in NASH - A key paradigm for understanding pathophysiology and clinical outcomes

Non-alcoholic fatty liver disease (NAFLD) exists as a spectrum ranging from simple steatosis to histologically defined hepatocyte injury and inflammatory changes that define steatohepatitis (NASH), and increase risk for fibrosis. Although zonal differences in NASH have not been systematically studied, periportal involvement has been associated with worse metabolic outcomes and more hepatic fibrosis as compared to pericentral disease. These data suggest that hepatic zonation of disease may influence the diversity of clinical presentations. Similarly, several randomized clinical trials suggest a differential response based on zonation of disease, with preferential effects on periportal (cysteamine) or pericentral disease (obeticholic acid, pioglitazone). Intriguingly, morphogenic pathways known to affect zonal development and maintenance - WNT/β-Catenin, Hedgehog, HIPPO/Yap/TAZ and Notch - have been implicated in NASH pathogenesis, and nuclear hormone receptors downstream of potential NASH therapeutics show zonal preferences. In this review, we summarize these data and propose that patient-specific activation of these pathways may explain the variability in clinical presentation, and the zone-specific response observed in clinical trials.

Integrating single-cell and spatial transcriptomics to elucidate intercellular tissue dynamics

Single-cell RNA sequencing (scRNA-seq) identifies cell subpopulations within tissue but does not capture their spatial distribution nor reveal local networks of intercellular communication acting in situ. A suite of recently developed techniques that localize RNA within tissue, including multiplexed in situ hybridization and in situ sequencing (here defined as high-plex RNA imaging) and spatial barcoding, can help address this issue. However, no method currently provides as complete a scope of the transcriptome as does scRNA-seq, underscoring the need for approaches to integrate single-cell and spatial data. Here, we review efforts to integrate scRNA-seq with spatial transcriptomics, including emerging integrative computational methods, and propose ways to effectively combine current methodologies.

SEAM is a spatial single nuclear metabolomics method for dissecting tissue microenvironment

Spatial metabolomics can reveal intercellular heterogeneity and tissue organization. Here we report on the spatial single nuclear metabolomics (SEAM) method, a flexible platform combining high-spatial-resolution imaging mass spectrometry and a set of computational algorithms that can display multiscale and multicolor tissue tomography together with identification and clustering of single nuclei by their in situ metabolic fingerprints. We first applied SEAM to a range of wild-type mouse tissues, then delineated a consistent pattern of metabolic zonation in mouse liver. We further studied the spatial metabolic profile in the human fibrotic liver. We discovered subpopulations of hepatocytes with special metabolic features associated with their proximity to the fibrotic niche, and validated this finding by spatial transcriptomics with Geo-seq. These demonstrations highlighted SEAM's ability to explore the spatial metabolic profile and tissue histology at the single-cell level, leading to a deeper understanding of tissue metabolic organization.

Signalling pathways and transcriptional regulators orchestrating liver development and cancer

Liver development is controlled by key signals and transcription factors that drive cell proliferation, migration, differentiation and functional maturation. In the adult liver, cell maturity can be perturbed by genetic and environmental factors that disrupt hepatic identity and function. Developmental signals and fetal genetic programmes are often dysregulated or reactivated, leading to dedifferentiation and disease. Here, we highlight signalling pathways and transcriptional regulators that drive liver cell development and primary liver cancers. We also discuss emerging models derived from pluripotent stem cells, 3D organoids and bioengineering for improved studies of signalling pathways in liver cancer and regenerative medicine.

Epigenomic and transcriptional profiling identifies impaired glyoxylate detoxification in NAFLD as a risk factor for hyperoxaluria

Epigenetic modifications (e.g. DNA methylation) in NAFLD and their contribution to disease progression and extrahepatic complications are poorly explored. Here, we use an integrated epigenome and transcriptome analysis of mouse NAFLD hepatocytes and identify alterations in glyoxylate metabolism, a pathway relevant in kidney damage via oxalate release-a harmful waste product and kidney stone-promoting factor. Downregulation and hypermethylation of alanine-glyoxylate aminotransferase (Agxt), which detoxifies glyoxylate, preventing excessive oxalate accumulation, is accompanied by increased oxalate formation after metabolism of the precursor hydroxyproline. Viral-mediated Agxt transfer or inhibiting hydroxyproline catabolism rescues excessive oxalate release. In human steatotic hepatocytes, AGXT is also downregulated and hypermethylated, and in NAFLD adolescents, steatosis severity correlates with urinary oxalate excretion. Thus, this work identifies a reduced capacity of the steatotic liver to detoxify glyoxylate, triggering elevated oxalate, and provides a mechanistic explanation for the increased risk of kidney stones and chronic kidney disease in NAFLD patients.

Nonalcoholic fatty liver disease stratification by liver lipidomics

Nonalcoholic fatty liver disease (NAFLD) is a common metabolic dysfunction leading to hepatic steatosis. However, NAFLD's global impact on the liver lipidome is poorly understood. Using high-resolution shotgun mass spectrometry, we quantified the molar abundance of 316 species from 22 major lipid classes in liver biopsies of 365 patients, including nonsteatotic patients with normal or excessive weight, patients diagnosed with NAFL (nonalcoholic fatty liver) or NASH (nonalcoholic steatohepatitis), and patients bearing common mutations of NAFLD-related protein factors. We confirmed the progressive accumulation of di- and triacylglycerols and cholesteryl esters in the liver of NAFL and NASH patients, while the bulk composition of glycerophospho- and sphingolipids remained unchanged. Further stratification by biclustering analysis identified sphingomyelin species comprising n24:2 fatty acid moieties as membrane lipid markers of NAFLD. Normalized relative abundance of sphingomyelins SM 43:3;2 and SM 43:1;2 containing n24:2 and n24:0 fatty acid moieties, respectively, showed opposite trends during NAFLD progression and distinguished NAFL and NASH lipidomes from the lipidome of nonsteatotic livers. Together with several glycerophospholipids containing a C22:6 fatty acid moiety, these lipids serve as markers of early and advanced stages of NAFL.

Loss of hepatic Mboat7 leads to liver fibrosis

OBJECTIVE

The rs641738C>T variant located near the membrane-bound O-acyltransferase domain containing 7 (MBOAT7) locus is associated with fibrosis in liver diseases, including non-alcoholic fatty liver disease (NAFLD), alcohol-related liver disease, hepatitis B and C. We aim to understand the mechanism by which the rs641738C>T variant contributes to pathogenesis of NAFLD.

DESIGN

Mice with hepatocyte-specific deletion of MBOAT7 (Mboat7Δhep) were generated and livers were characterised by histology, flow cytometry, qPCR, RNA sequencing and lipidomics. We analysed the association of rs641738C>T genotype with liver inflammation and fibrosis in 846 NAFLD patients and obtained genotype-specific liver lipidomes from 280 human biopsies.

RESULTS

Allelic imbalance analysis of heterozygous human liver samples pointed to lower expression of the MBOAT7 transcript on the rs641738C>T haplotype. Mboat7Δhep mice showed spontaneous steatosis characterised by increased hepatic cholesterol ester content after 10 weeks. After 6 weeks on a high fat, methionine-low, choline-deficient diet, mice developed increased hepatic fibrosis as measured by picrosirius staining (p<0.05), hydroxyproline content (p<0.05) and transcriptomics, while the inflammatory cell populations and inflammatory mediators were minimally affected. In a human biopsied NAFLD cohort, MBOAT7 rs641738C>T was associated with fibrosis (p=0.004) independent of the presence of histological inflammation. Liver lipidomes of Mboat7Δhep mice and human rs641738TT carriers with fibrosis showed increased total lysophosphatidylinositol levels. The altered lysophosphatidylinositol and phosphatidylinositol subspecies in MBOAT7Δhep livers and human rs641738TT carriers were similar.

CONCLUSION

Mboat7 deficiency in mice and human points to an inflammation-independent pathway of liver fibrosis that may be mediated by lipid signalling and a potentially targetable treatment option in NAFLD.

PPARs in liver physiology

Peroxisome proliferator-activated receptors (PPARs) are nuclear receptors and transcriptional modulators with crucial functions in hepatic and whole-body energy homeostasis. Besides their well-documented roles in lipid and glucose metabolism, emerging evidence also implicate PPARs in the control of other processes such as inflammatory responses. Recent technological advances, such as single-cell RNA sequencing, have allowed to unravel an unexpected complexity in the regulation of PPAR expression, activity and downstream signaling. Here we provide an overview of the latest advances in the study of PPARs in liver physiology, with a specific focus on formerly neglected aspects of PPAR regulation, such as tissular zonation, cellular heterogeneity, circadian rhythms, sexual dimorphism and species-specific features.

Hepatocyte size fractionation allows dissection of human liver zonation

Human hepatocytes show marked differences in cell size, gene expression, and function throughout the liver lobules, an arrangement termed liver zonation. However, it is not clear if these zonal size differences, and the associated phenotypic differences, are retained in isolated human hepatocytes, the "gold standard" for in vitro studies of human liver function. Here, we therefore explored size differences among isolated human hepatocytes and investigated whether separation by size can be used to study liver zonation in vitro. We used counterflow centrifugal elutriation to separate cells into different size fractions and analyzed them with label-free quantitative proteomics, which revealed an enrichment of 151 and 758 proteins (out of 5163) in small and large hepatocytes, respectively. Further analysis showed that protein abundances in different hepatocyte size fractions recapitulated the in vivo expression patterns of previously described zonal markers and biological processes. We also found that the expression of zone-specific cytochrome P450 enzymes correlated with their metabolic activity in the different fractions. In summary, our results show that differences in hepatocyte size matches zonal expression patterns, and that our size fractionation approach can be used to study zone-specific liver functions in vitro.

Epigenetics of Hepatic Insulin Resistance

Insulin resistance (IR) is largely recognized as a unifying feature that underlies metabolic dysfunction. Both lifestyle and genetic factors contribute to IR. Work from recent years has demonstrated that the epigenome may constitute an interface where different signals may converge to promote IR gene expression programs. Here, we review the current knowledge of the role of epigenetics in hepatic IR, focusing on the roles of DNA methylation and histone post-translational modifications. We discuss the broad epigenetic changes observed in the insulin resistant liver and its associated pathophysiological states and leverage on the wealth of 'omics' studies performed to discuss efforts in pinpointing specific loci that are disrupted by these changes. We envision that future studies, with increased genomic resolution and larger cohorts, will further the identification of biomarkers of early onset hepatic IR and assist the development of targeted interventions. Furthermore, there is growing evidence to suggest that persistent epigenetic marks may be acquired over prolonged exposure to disease or deleterious exposures, highlighting the need for preventative medicine and long-term lifestyle adjustments to avoid irreversible or long-term alterations in gene expression.

Mutual Zonated Interactions of Wnt and Hh Signaling Are Orchestrating the Metabolism of the Adult Liver in Mice and Human

The Hedgehog (Hh) and Wnt/β-Catenin (Wnt) cascades are morphogen pathways whose pronounced influence on adult liver metabolism has been identified in recent years. How both pathways communicate and control liver metabolic functions are largely unknown. Detecting core components of Wnt and Hh signaling and mathematical modeling showed that both pathways in healthy liver act largely complementary to each other in the pericentral (Wnt) and the periportal zone (Hh) and communicate mainly by mutual repression. The Wnt/Hh module inversely controls the spatiotemporal operation of various liver metabolic pathways, as revealed by transcriptome, proteome, and metabolome analyses. Shifting the balance to Wnt (activation) or Hh (inhibition) causes pericentralization and periportalization of liver functions, respectively. Thus, homeostasis of the Wnt/Hh module is essential for maintaining proper liver metabolism and to avoid the development of certain metabolic diseases. With caution due to minor species-specific differences, these conclusions may hold for human liver as well.

Targeting Ca2+ Signaling in the Initiation, Promotion and Progression of Hepatocellular Carcinoma

Simple Summary

Liver cancer (hepatocellular carcinoma) is a significant health burden worldwide. It is often not detected until at an advanced stage when there are few treatment options available. Changes in calcium concentrations within liver cancer cells are essential for regulating their growth, death, and migration (metastasis). Our aim was to review published papers which have identified proteins involved in calcium signaling as potential drug targets for the treatment of liver cancer. About twenty calcium signaling proteins were identified, including those involved in regulating calcium concentrations in the cytoplasm, endoplasmic reticulum and mitochondria. A few of these have turned out to be sites of action of natural products previously known to inhibit liver cancer. More systematic studies are now needed to determine which calcium signaling proteins might be used clinically for treatment of liver cancer, especially advanced stage cancers and those resistant to inhibition by current drugs.

Abstract

Hepatocellular carcinoma (HCC) is a considerable health burden worldwide and a major contributor to cancer-related deaths. HCC is often not noticed until at an advanced stage where treatment options are limited and current systemic drugs can usually only prolong survival for a short time. Understanding the biology and pathology of HCC is a challenge, due to the cellular and anatomic complexities of the liver. While not yet fully understood, liver cancer stem cells play a central role in the initiation and progression of HCC and in resistance to drugs. There are approximately twenty Ca²⁺-signaling proteins identified as potential targets for therapeutic treatment at different stages of HCC. These potential targets include inhibition of the self-renewal properties of liver cancer stem cells; HCC initiation and promotion by hepatitis B and C and non-alcoholic fatty liver disease (principally involving reduction of reactive oxygen species); and cell proliferation, tumor growth, migration and metastasis. A few of these Ca²⁺-signaling pathways have been identified as targets for natural products previously known to reduce HCC. Promising Ca²⁺-signaling targets include voltage-operated Ca²⁺ channel proteins (liver cancer stem cells), inositol trisphosphate receptors, store-operated Ca²⁺ entry, TRP channels, sarco/endoplasmic reticulum (Ca²⁺+Mg²⁺) ATP-ase and Ca²⁺/calmodulin-dependent protein kinases. However, none of these Ca²⁺-signaling targets has been seriously studied any further than laboratory research experiments. The future application of more systematic studies, including genomics, gene expression (RNA-seq), and improved knowledge of the fundamental biology and pathology of HCC will likely reveal new Ca²⁺-signaling protein targets and consolidate priorities for those already identified.

Review article: the signalling and functional role of the extracellular matrix in the development of liver fibrosis

BACKGROUND

Patients with liver fibrosis show a large heterogeneity, and for that reason effective treatments are still lacking. Emerging data suggest that there is more to fibrosis than previously understood. Opposed to earlier belief of being a passive scaffold for cells to reside in, the extracellular matrix (ECM) is now known to hold both signalling and functional properties important for the development of fibrosis. The interaction between the ECM and the collagen-producing cells determines the course of the disease but is still poorly understood. Exploring the dynamics of this interplay will aid in the development of effective treatments.

AIM

To summarise and discuss the latest advances in the pathogenesis of liver fibrosis as well as key mediators of early disease progression.

METHODS

Through literature search using databases including PubMed and Google Scholar, manuscripts published between 1961 and 2019 were included to assess both well-established and recent theories of fibrosis development. Both pre-clinical and clinical studies were included.

RESULTS

Fibrosis alters the structure of the ECM releasing signalling fragments with the potential to escalate disease severity. In a diseased liver, hepatic stellate cells and other fibroblasts, together with hepatocytes and sinusoidal cells, produce an excessive amount of collagens. The cell-to-collagen interactions are unique in the different liver aetiologies, generating ECM profiles with considerable patient-monitoring potential.

CONCLUSIONS

The local milieu in the injured area affects the course of fibrosis development in a site-specific manner. Future research should focus on the dissimilarities in the ECM profile between different aetiologies of liver fibrosis.

Cell-type-resolved proteomic analysis of the human liver

BACKGROUND & AIMS

The human liver functions through a complex interplay between parenchymal and non-parenchymal cells. Mass spectrometry-based proteomic analysis of intact tissue has provided an in-depth view of the human liver proteome. However, the predominance of parenchymal cells (hepatocytes) means that the total tissue proteome mainly reflects hepatocyte expression. Here we therefore set out to analyse the proteomes of the major parenchymal and non-parenchymal cell types in the human liver.

METHODS

We applied quantitative label-free proteomic analysis on the major cell types of the human liver: hepatocytes, liver endothelial cells, Kupffer cells and hepatic stellate cells.

RESULTS

We identified 9791 proteins, revealing distinct protein expression profiles across cell types, whose in vivo relevance was shown by the presence of cell-type-specific proteins. Analysis of proteins related to the immune system indicated that mechanisms of immune-mediated liver injury include the involvement of several cell types. Furthermore, in-depth investigation of proteins related to the absorption, distribution, metabolism, excretion and toxicity (ADMET) of xenobiotics showed that ADMET-related tasks are not exclusively confined to hepatocytes, and that non-parenchymal cells may contribute to drug transport and metabolism.

CONCLUSIONS

Overall, the data we provide constitute a unique resource for exploring the proteomes of the major types of human liver cells, which will facilitate an improved understanding of the human liver in health and disease.

From genome sequencing to the discovery of potential biomarkers in liver disease

Chronic liver disease progresses through several stages, fatty liver, steatohepatitis, cirrhosis, and eventually, it leads to hepatocellular carcinoma (HCC) over a long period of time. Since a large proportion of patients with HCC are accompanied by cirrhosis, it is considered to be an important factor in the diagnosis of liver cancer. This is because cirrhosis leads to an irreversible harmful effect, but the early stages of chronic liver disease could be reversed to a healthy state. Therefore, the discovery of biomarkers that could identify the early stages of chronic liver disease is important to prevent serious liver damage. Biomarker discovery at liver cancer and cirrhosis has enhanced the development of sequencing technology. Next generation sequencing (NGS) is one of the representative technical innovations in the biological field in the recent decades and it is the most important thing to design for research on what type of sequencing methods are suitable and how to handle the analysis steps for data integration. In this review, we comprehensively summarized NGS techniques for identifying genome, transcriptome, DNA methylome and 3D/4D chromatin structure, and introduced framework of processing data set and integrating multi-omics data for uncovering biomarkers.

An in Situ Atlas of Mitochondrial DNA in Mammalian Tissues Reveals High Content in Stem and Proliferative Compartments

Mitochondria regulate ATP production, metabolism, and cell death. Alterations in mitochondrial DNA (mtDNA) sequence and copy number are implicated in aging and organ dysfunction in diverse inherited and sporadic diseases. Because most measurements of mtDNA use homogenates of complex tissues, little is known about cell-type-specific mtDNA copy number heterogeneity in normal physiology, aging, and disease. Thus, the precise cell types whose loss of mitochondrial activity and altered mtDNA copy number that result in organ dysfunction in aging and disease have often not been clarified. Here, an in situ hybridization approach to generate a single-cell-resolution atlas of mtDNA content in mammalian tissues was validated. In hierarchically organized self-renewing tissues, higher levels of mtDNA were observed in stem/proliferative compartments compared with differentiated compartments. Striking zonal patterns of mtDNA levels in the liver reflected the known oxygen tension gradient. In the kidney, proximal and distal tubules had markedly higher mtDNA levels compared with cells within glomeruli and collecting duct epithelial cells. In mice, decreased mtDNA levels were visualized in renal tubules as a function of aging, which was prevented by calorie restriction. This study provides a novel approach for quantifying species- and cell-type-specific mtDNA copy number and dynamics in any normal or diseased tissue that can be used for monitoring the effects of interventions in animal and human studies.

Glyphboard: Visual Exploration of High-Dimensional Data Combining Glyphs with Dimensionality Reduction

Rigorous data science is interdisciplinary at its core. In order to make sense of high-dimensional data, data scientists need to enter into a dialogue with domain experts. We present Glyphboard, a visualization tool that aims to support this dialogue. Glyphboard is a zoomable user interface that combines well-known methods such as dimensionality reduction and glyph-based visualizations in a novel, seamless, and integrated tool. While the dimensionality reduction affords a quick overview over the data, glyph-based visualizations are able to show the most relevant dimensions in the data set at one glance. We contribute an open-source prototype of Glyphboard, a general exchange format for high-dimensional data, and a case study with nine data scientists and domain experts from four exemplary domains in order to evaluate how the different visualization and interaction features of Glyphboard are used.

The Epigenetic Drug Discovery Landscape for Metabolic-associated Fatty Liver Disease

Despite decades of research, effective therapies for metabolic (dysfunction)-associated fatty liver disease (MAFLD) are lacking. An increasing body of evidence suggests that epigenetic dysregulation is frequent in MAFLD, and orchestrates many aspects of its development and progression. Furthermore, the high plasticity of epigenetic modifications in response to environmental cues renders epigenetics a novel area for therapeutic drug discovery. Over recent years, several epigenetics-based drugs and diagnostic biomarkers have entered clinical development and/or obtained regulatory approval. Here, we review recent advances in our understanding of epigenetic regulation and programming during MAFLD, including DNA methylation, histone modifications, chromatin remodelling, transcriptional control, and noncoding (nc)RNAs. We also discuss the potential translational implications and challenges of epigenetics in the context of MAFLD.

Liver gene regulatory networks: Contributing factors to nonalcoholic fatty liver disease

Metabolic diseases such as nonalcoholic fatty liver disease (NAFLD) result from complex interactions between intrinsic and extrinsic factors, including genetics and exposure to obesogenic environments. These risk factors converge in aberrant gene expression patterns in the liver, which are underlined by altered cis-regulatory networks. In homeostasis and in disease states, liver cis-regulatory networks are established by coordinated action of liver-enriched transcription factors (TFs), which define enhancer landscapes, activating broad gene programs with spatiotemporal resolution. Recent advances in DNA sequencing have dramatically expanded our ability to map active transcripts, enhancers and TF cistromes, and to define the 3D chromatin topology that contains these elements. Deployment of these technologies has allowed investigation of the molecular processes that regulate liver development and metabolic homeostasis. Moreover, genomic studies of NAFLD patients and NAFLD models have demonstrated that the liver undergoes pervasive regulatory rewiring in NAFLD, which is reflected by aberrant gene expression profiles. We have therefore achieved an unprecedented level of detail in the understanding of liver cis-regulatory networks, particularly in physiological conditions. Future studies should aim to map active regulatory elements with added levels of resolution, addressing how the chromatin landscapes of different cell lineages contribute to and are altered in NAFLD and NAFLD-associated metabolic states. Such efforts would provide additional clues into the molecular factors that trigger this disease. This article is categorized under: Biological Mechanisms > Metabolism Biological Mechanisms > Regulatory Biology Laboratory Methods and Technologies > Genetic/Genomic Methods.

Graphene Oxide Causes Disordered Zonation Due to Differential Intralobular Localization in the Liver

The liver is the primary organ to sequester nanodrugs, representing a substantial hurdle for drug delivery and raising toxicity concerns. However, the mechanistic details underlying the liver sequestration and effects on the liver are still elusive. The difficulty in studying the liver lies in its complexity, which is structured with stringently organized anatomical units called lobules. Graphene oxide (GO) has attracted attention for its applications in biomedicine, especially as a nanocarrier; however, its sequestration and effects in the liver, the major enrichment and metabolic organ, are less understood. Herein, we unveiled the differential distribution of GO in lobules in the liver, with a higher amount surrounding portal triad zones than the central vein zones. Strikingly, liver zonation patterns also changed, as reflected by changes in vital zonated genes involved in hepatocyte integrity and metabolism, leading to compromised hepatic functions. RNA-Seq and DNA methylation sequencing analyses unraveled that GO-induced changes in liver functional zonation could be ascribed to dysregulation of key signaling pathways governing liver zonation at not only mRNA transcriptions but also DNA methylation imprinting patterns, partially through TET-dependent signaling. Together, this study reveals the differential GO distribution pattern in liver lobules and pinpoints the genetic and epigenetic mechanisms in GO-induced liver zonation alterations.

The complement system in liver diseases: Evidence-based approach and therapeutic options

Complement is usually seen to largely originate from the liver to accomplish its tasks systemically - its return to the production site has long been underestimated. Recent progress in genomics, therapeutic effects on complement, standardised possibilities in medical laboratory tests and involvement of complosome brings the complement system with its three major functions of opsonization, cytolysis and phagocytosis back to liver biology and pathology. The LOINC™ system features 20 entries for the C3 component of complement to anticipate the application of artificial intelligence data banks algorythms of which are fed with patient-specific data connected to standard lab assays for liver function. These advancements now lead to increased vigilance by clinicians. This reassessment article will further elucidate the distribution of synthesis sites to the three germ layer-derived cell systems and the role complement now known to play in embryogenesis, senescence, allotransplantation and autoimmune disease. This establishes the liver as part of the gastro-intestinal system in connection with nosological entities never thought of, such as the microbiota-liver-brain axis. In neurological disease etiology infectious and autoimmune hepatitis play an important role in the context of causative viz reactive complement activation. The mosaic of autoimmunity, i.e. multiple combinations of the many factors producing varying clinical pictures, leads to the manifold facets of liver autoimmunity.