A human liver cell atlas reveals heterogeneity and epithelial progenitors

Single-cell multi-omics in the study of digestive system cancers

Digestive system cancers are prevalent diseases with a high mortality rate, posing a significant threat to public health and economic burden. The diagnosis and treatment of digestive system cancer confront conventional cancer problems, such as tumor heterogeneity and drug resistance. Single-cell sequencing (SCS) emerged at times required and has developed from single-cell RNA-seq (scRNA-seq) to the single-cell multi-omics era represented by single-cell spatial transcriptomics (ST). This article comprehensively reviews the advances of single-cell omics technology in the study of digestive system tumors. While analyzing and summarizing the research cases, vital details on the sequencing platform, sample information, sampling method, and key findings are provided. Meanwhile, we summarize the commonly used SCS platforms and their features, as well as the advantages of multi-omics technologies in combination. Finally, the development trends and prospects of the application of single-cell multi-omics technology in digestive system cancer research are prospected.

Farnesoid X receptor: From Structure to Function and Its Pharmacology in Liver Fibrosis

The farnesoid X receptor (FXR), a ligand-activated transcription factor, plays a crucial role in regulating bile acid metabolism within the enterohepatic circulation. Beyond its involvement in metabolic disorders and immune imbalances affecting various tissues, FXR is implicated in microbiota modulation, gut-to-brain communication, and liver disease. The liver, as a pivotal metabolic and detoxification organ, is susceptible to damage from factors such as alcohol, viruses, drugs, and high-fat diets. Chronic or recurrent liver injury can culminate in liver fibrosis, which, if left untreated, may progress to cirrhosis and even liver cancer, posing significant health risks. However, therapeutic options for liver fibrosis remain limited in terms of FDA-approved drugs. Recent insights into the structure of FXR, coupled with animal and clinical investigations, have shed light on its potential pharmacological role in hepatic fibrosis. Progress has been achieved in both fundamental research and clinical applications. This review critically examines recent advancements in FXR research, highlighting challenges and potential mechanisms underlying its role in liver fibrosis treatment.

The Pathophysiology of Portal Hypertension

This article reviews the pathophysiology of portal hypertension that includes multiple mechanisms internal and external to the liver. This article starts with a review of literature describing the cellular and molecular mechanisms of portal hypertension, microvascular thrombosis, sinusoidal venous congestion, portal angiogenesis, vascular hypocontractility, and hyperdynamic circulation. Mechanotransduction and the gut-liver axis, which are newer areas of research, are reviewed. Dysfunction of this axis contributes to chronic liver injury, inflammation, fibrosis, and portal hypertension. Sequelae of portal hypertension are discussed in subsequent studies.

Establishing an oxidative stress mitochondria-related prognostic model in hepatocellular carcinoma based on multi-omics characteristics and machine learning computational framework

Hepatocellular carcinoma (HCC) has high incidence and mortality rates worldwide. Damaged mitochondria are characterized by the overproduction of reactive oxygen species (ROS), which can promote cancer development. The prognostic value of the interplay between mitochondrial function and oxidative stress in HCC requires further investigation. Gene expression data of HCC samples were collected from The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO) and International Cancer Genome Consortium (ICGC). We screened prognostic oxidative stress mitochondria-related (OSMT) genes at the bulk transcriptome level. Based on multiple machine learning algorithms, we constructed a consensus oxidative stress mitochondria-related signature (OSMTS), which contained 26 genes. In addition, we identified six of these genes as having a suitable prognostic value for OSMTS to reduce the difficulty of clinical application. Univariate and multivariate analyses verified the OSMTS as an independent prognostic factor for overall survival (OS) in HCC patients. The OSMTS-related nomogram demonstrated to be a powerful tool for the clinical diagnosis of HCC. We observed differences in biological function and immune cell infiltration in the tumor microenvironment between the high- and low-risk groups. The highest expression of the OSMTS was detected in hepatocytes at the single-cell transcriptome level. Hepatocytes in the high- and low-risk groups differed significantly in terms of biological function and intercellular communication. Moreover, at the spatial transcriptome level, high expression of OSMTS was mainly in regions enriched in hepatocytes and B cells. Potential drugs targeting specific risk subgroups were identified. Our study revealed that the OSMTS can serve as a promising tool for prognosis prediction and precise intervention in HCC patients.

Research trends and hotspots in the immune microenvironment related to hepatocellular carcinoma: A bibliometric and visualization study

BACKGROUND

The immune microenvironment (IME) in hepatocellular carcinoma (HCC) plays a pivotal role in determining patient outcomes and responses to treatment. This area is witnessing rapid growth in research interest. However, there is a lack of comprehensive bibliometric analyses that dissect trends and potential focal points in this field.

AIM

To explore the evolution of research on the IME in HCC from January 1, 2004, to December 31, 2023, using bibliometric methodologies.

METHODS

English articles and reviews concerning the IME of HCC were retrieved from the Web of Science Core Collection with a search date of December 31, 2023. The R package Bibliometrix was employed to compute basic bibliometric characteristics, illustrate collaborations among countries and authors, and create a three-field diagram illustrating the connections between authors, affiliations, and keywords. Analyses of country and institutional co-authorship, as well as keyword co-occurrence, were conducted using VOSviewer. Additionally, CiteSpace was utilized for the cite burst analysis of keywords and cited literature.

RESULTS

The study encompassed 3125 documents in the research areas related to HCC of IME, revealing a substantial and continuous increase in the annual publication trend over time. China and Fudan University emerged as leading contributors, with 2103 and 165 publications, respectively. Frontiers in immunology was the most prolific journal in this domain. Among the top ten researchers in the field, eight are based in China. Key research terms identified include tumour microenvironment, expression, immunotherapy, and prognosis.

CONCLUSION

The relationship between HCC and IME is receiving increasing attention, and related research is in a highly developed stage. Key focus areas, including IME and immune checkpoint inhibitors, immunotherapy are poised to be central to future research endeavors, offering promising pathways for further exploration.

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.

Synchronized lineage tracing of cell membranes and nuclei by dual recombinases and dual fluorescent

Genetic lineage tracing has been widely employed to investigate cell lineages and fate. However, conventional reporting systems often label the entire cytoplasm, making it challenging to discern cell boundaries. Additionally, single Cre-loxP recombination systems have limitations in tracing specific cell populations. This study proposes three reporting systems that utilize Cre, Dre and Dre + Cre mediated recombination. These systems incorporate tdTomato expression on the cell membrane and PhiYFP expression within the cell nucleus, allowing for clear observation of the cell nucleus and membrane. The efficacy of these systems is successfully demonstrated by labeling cardiomyocytes and hepatocytes. The potential for dynamic visualization of the cell membrane is showcased using intravital imaging microscopy or three-dimensional imaging. Furthermore, by combining this dual recombinase system with the ProTracer system, hepatocyte proliferation is traced with enhanced precision. This reporting system holds significant importance in advancing the understanding of cell fate studies in development, homeostasis, and diseases.

Enhanced in situ spatial proteomics by effective combination of MALDI imaging and LC-MS/MS

Highly specialized cells are fundamental for proper functioning of complex organs. Variations in cell-type specific gene expression and protein composition have been linked to a variety of diseases. Investigation of the distinctive molecular makeup of these cells within tissues is therefore critical in biomedical research. Although several technologies have emerged as valuable tools to address this cellular heterogeneity, most workflows lack sufficient in situ resolution and are associated with high cost and extremely long analysis times. Here, we present a combination of experimental and computational approaches that allows a more comprehensive investigation of molecular heterogeneity within tissues than by either shotgun LC-MS/MS or MALDI imaging alone. We applied our pipeline on mouse brain, which contains a wide variety of cell types that not only perform unique functions but also exhibit varying sensitivities to insults. We explored the distinct neuronal populations within the hippocampus, a brain region crucial for learning and memory that is involved in various neurological disorders. As an example, we identified the groups of proteins distinguishing the neuronal populations of dentate gyrus (DG) and the cornu ammonis (CA) in the same brain section. Most of the annotated proteins matched the regional enrichment of their transcripts, thereby validating the method. As the method is highly reproducible, the identification of individual masses through the combination of MALDI-IMS and LC-MS/MS methods can be used for the much faster and more precise interpretation of MALDI-IMS measurements only. This greatly speeds up spatial proteomic analyses and allows the detection of local protein variations within the same population of cells. The method's general applicability has the potential to be used to investigate different biological conditions and tissues and a much higher throughput than other techniques making it a promising approach for clinical routine applications.

Effect of Cellular Senescence in Disease Progression and Transplantation: Immune Cells and Solid Organs

Aging of the world population significantly impacts healthcare globally and specifically, the field of transplantation. Together with end-organ dysfunction and prolonged immunosuppression, age increases the frequency of comorbid chronic diseases in transplant candidates and recipients, contributing to inferior outcomes. Although the frequency of death increases with age, limited use of organs from older deceased donors reflects the concerns about organ durability and inadequate function. Cellular senescence (CS) is a hallmark of aging, which occurs in response to a myriad of cellular stressors, leading to activation of signaling cascades that stably arrest cell cycle progression to prevent tumorigenesis. In aging and chronic conditions, senescent cells accumulate as the immune system's ability to clear them wanes, which is causally implicated in the progression of chronic diseases, immune dysfunction, organ damage, decreased regenerative capacity, and aging itself. The intimate interplay between senescent cells, their proinflammatory secretome, and immune cells results in a positive feedback loop, propagating chronic sterile inflammation and the spread of CS. Hence, senescent cells in organs from older donors trigger the recipient's alloimmune response, resulting in the increased risk of graft loss. Eliminating senescent cells or attenuating their inflammatory phenotype is a novel, potential therapeutic target to improve transplant outcomes and expand utilization of organs from older donors. This review focuses on the current knowledge about the impact of CS on circulating immune cells in the context of organ damage and disease progression, discusses the impact of CS on abdominal solid organs that are commonly transplanted, and reviews emerging therapies that target CS.

Aberrant differentiation and proliferation of hepatocytes in chronic liver injury and liver tumors

Chronic liver injury induces liver cirrhosis and facilitates hepatocarcinogenesis. However, the effects of this condition on hepatocyte proliferation and differentiation are unclear. We showed that rodent hepatocytes display a ductular phenotype when they are cultured within a collagenous matrix. This process involves transdifferentiation without the emergence of hepatoblastic features and is at least partially reversible. During the ductular reaction in chronic liver diseases with progressive fibrosis, some hepatocytes, especially those adjacent to ectopic ductules, demonstrate ductular transdifferentiation, but the majority of increased ductules originate from the existing bile ductular system that undergoes extensive remodeling. In chronic injury, hepatocyte proliferation is weak but sustained, and most regenerative nodules in liver cirrhosis are composed of clonally proliferating hepatocytes, suggesting that a small fraction of hepatocytes maintain their proliferative capacity in chronic injury. In mouse hepatocarcinogenesis models, hepatocytes activate the expression of various fetal/neonatal genes, indicating that these cells undergo dedifferentiation. Hepatocyte-specific somatic integration of various oncogenes in mice demonstrated that hepatocytes may be the cells of origin for a broad spectrum of liver tumors through transdifferentiation and dedifferentiation. In conclusion, the phenotypic plasticity and heterogeneity of mature hepatocytes are important for understanding the pathogenesis of chronic liver diseases and liver tumors.

ScHGSC-IGDC: Identifying genes with differential correlations of high-grade serous ovarian cancer based on single-cell RNA sequencing analysis

Due to the high heterogeneity of ovarian cancer (OC), it occupies the main cause of cancer-related death among women. As the most aggressive and frequent subtype of OC, high-grade serous cancer (HGSC) represents around 70 % of all patients. With the booming progress of single-cell RNA sequencing (scRNA-seq), unique and subtle changes among different cell states have been identified including novel risk genes and pathways. Here, our present study aims to identify differentially correlated core genes between normal and tumor status through HGSC scRNA-seq data analysis. R package high-dimension Weighted Gene Co-expression Network Analysis (hdWGCNA) was implemented for building gene interaction networks based on HGSC scRNA-seq data. DiffCorr was integrated for identifying differentially correlated genes between tumor and their adjacent normal counterparts. Software Cytoscape was implemented for constructing and visualizing biological networks. Real-time qPCR (RT-qPCR) was utilized to confirm expression pattern of new genes. We introduced ScHGSC-IGDC (Identifying Genes with Differential Correlations of HGSC based on scRNA-seq analysis), an in silico framework for identifying core genes in the development of HGSC. We detected thirty-four modules in the network. Scores of new genes with opposite correlations with others such as NDUFS5, TMSB4X, SERPINE2 and ITPR2 were identified. Further survival and literature validation emphasized their great values in the HGSC management. Meanwhile, RT-qPCR verified expression pattern of NDUFS5, TMSB4X, SERPINE2 and ITPR2 in human OC cell lines and tissues. Our research offered novel perspectives on the gene modulatory mechanisms from single cell resolution, guiding network based algorithms in cancer etiology field.

Single-Cell Transcriptomic Analysis Revealed the Cell Population Changes and Cell-Cell Communication in the Liver of a Carnivorous Fish in Response to High-Carbohydrate Diet

BACKGROUND

Carnivorous fish have a low carbohydrate utilization ability, and the physiologic and molecular basis of glucose intolerance has not been fully illustrated.

OBJECTIVES

This study aimed to use largemouth bass as a model to investigate the possible mechanism of glucose intolerance in carnivorous fish with the help of single-nuclei RNA sequencing (snRNA-seq).

METHODS

Two diets were formulated, a low-carbohydrate (LC) diet and a high-carbohydrate (HC) diet. The feeding trial lasted for 6 wk, and then, growth performance, biochemical parameters, liver histology, and snRNA-seq were performed.

RESULTS

Growth performance of fish was not affected by the HC diet, while liver glucolipid metabolism disorder and liver injury were observed. A total of 13,247 and 12,848 cells from the liver derived from 2 groups were isolated and sequenced, and 7 major liver cell types were annotated by the marker genes. Hepatocytes and cholangiocytes were lower and hepatic stellate cells (HSCs) and immune cells were higher in the HC group than those in the LC group. Reclustering analysis identified 7 subtypes of hepatocytes and immune cells, respectively. The HSCs showed more cell communication with other cell types, and periportal hepatocytes showed more cell communication with other hepatocyte subtypes. Cell-cell communication mainly focused on cell junction-related signaling pathways. Uncovered by the pseudotime analysis, midzonal hepatocytes were differentiated into 2 major branches-biliary epithelial hepatocytes and hepatobiliary hybrid progenitor. Cell junction and liver fibrosis-related genes were highly expressed in the HC group. HC diet induced the activation of HSCs and, therefore, led to the liver fibrosis of largemouth bass.

CONCLUSIONS

HC diet induces liver glucolipid metabolism disorder and liver injury of largemouth bass. The increase and activation of HSCs might be the main reason for the liver injury. In adaption to HC diet, midzonal hepatocytes differentiates into 2 major branches-biliary epithelial hepatocytes and hepatobiliary hybrid progenitors.

Role of liver sinusoidal endothelial cell in metabolic dysfunction-associated fatty liver disease

Liver sinusoidal endothelial cells (LSECs) are highly specialized endothelial cells that represent the interface between blood cells on one side and hepatocytes on the other side. LSECs not only form a barrier within the hepatic sinus, but also play important physiological functions such as regulating hepatic vascular pressure, anti-inflammatory and anti-fibrotic. Pathologically, pathogenic factors can induce LSECs capillarization, that is, loss of fenestra and dysfunction, which are conducive to early steatosis, lay the foundation for the progression of metabolic dysfunction-associated fatty liver disease (MAFLD), and accelerate metabolic dysfunction-associated steatohepatitis (MASH) and liver fibrosis. The unique localization, phenotype, and function of LSECs make them potential candidates for reducing liver injury, inflammation, and preventing or reversing fibrosis in the future.

Single-Cell RNA Sequencing Reveals the Spatial Heterogeneity and Functional Alteration of Endothelial Cells in Chronic Hepatitis B Infection

Chronic Hepatitis B virus (CHB) infection is a global health challenge, causing damage ranging from hepatitis to cirrhosis and hepatocellular carcinoma. In our study, single-cell RNA sequencing (scRNA-seq) analysis was performed in livers from mice models with chronic inflammation induced by CHB infection and we found that endothelial cells (ECs) exhibited the largest number of differentially expressed genes (DEGs) among all ten cell types. NF-κB signaling was activated in ECs to induce cell dysfunction and subsequent hepatic inflammation, which might be mediated by the interaction of macrophage-derived and cholangiocyte-derived VISFATIN/Nampt signaling. Moreover, we divided ECs into three subclusters, including periportal ECs (EC_Z1), midzonal ECs (EC_Z2), and pericentral ECs (EC_Z3) according to hepatic zonation. Functional analysis suggested that pericentral ECs and midzonal ECs, instead of periportal ECs, were more vulnerable to HBV infection, as the VISFATIN/Nampt- NF-κB axis was mainly altered in these two subpopulations. Interestingly, pericentral ECs showed increasing communication with macrophages and cholangiocytes via the Nampt-Insr and Nampt-Itga5/Itgb1 axis upon CHB infection, which contribute to angiogenesis and vascular capillarization. Additionally, ECs, especially pericentral ECs, showed a close connection with nature killer (NK) cells and T cells via the Cxcl6-Cxcr6 axis, which is involved in shaping the microenvironment in CHB mice livers. Thus, our study described the heterogeneity and functional alterations of three subclusters in ECs. We revealed the potential role of VISFATIN/Nampt signaling in modulating ECs characteristics and related hepatic inflammation, and EC-derived chemokine Cxcl16 in shaping NK and T cell recruitment, providing key insights into the multifunctionality of ECs in CHB-associated pathologies.

The Advancement and Application of the Single-Cell Transcriptome in Biological and Medical Research

Single-cell RNA sequencing technology (scRNA-seq) has been steadily developing since its inception in 2009. Unlike bulk RNA-seq, scRNA-seq identifies the heterogeneity of tissue cells and reveals gene expression changes in individual cells at the microscopic level. Here, we review the development of scRNA-seq, which has gone through iterations of reverse transcription, in vitro transcription, smart-seq, drop-seq, 10 × Genomics, and spatial single-cell transcriptome technologies. The technology of 10 × Genomics has been widely applied in medicine and biology, producing rich research results. Furthermore, this review presents a summary of the analytical process for single-cell transcriptome data and its integration with other omics analyses, including genomes, epigenomes, proteomes, and metabolomics. The single-cell transcriptome has a wide range of applications in biology and medicine. This review analyzes the applications of scRNA-seq in cancer, stem cell research, developmental biology, microbiology, and other fields. In essence, scRNA-seq provides a means of elucidating gene expression patterns in single cells, thereby offering a valuable tool for scientific research. Nevertheless, the current single-cell transcriptome technology is still imperfect, and this review identifies its shortcomings and anticipates future developments. The objective of this review is to facilitate a deeper comprehension of scRNA-seq technology and its applications in biological and medical research, as well as to identify avenues for its future development in alignment with practical needs.

VEGF-FGF Signaling Activates Quiescent CD63+ Liver Stem Cells to Proliferate and Differentiate

Understanding the liver stem cells (LSCs) holds great promise for new insights into liver diseases and liver regeneration. However, the heterogenicity and plasticity of liver cells have made it controversial. Here, by employing single-cell RNA-sequencing technology, transcriptome features of Krt19+ bile duct lineage cells isolated from Krt19CreERT; Rosa26R-GFP reporter mouse livers are examined. Distinct biliary epithelial cells which include adult LSCs, as well as their downstream hepatocytes and cholangiocytes are identified. Importantly, a novel cell surface LSCs marker, CD63, as well as CD56, which distinguished active and quiescent LSCs are discovered. Cell expansion and bi-potential differentiation in culture demonstrate the stemness ability of CD63+ cells in vitro. Transplantation and lineage tracing of CD63+ cells confirm their contribution to liver cell mass in vivo upon injury. Moreover, CD63+CD56+ cells are proved to be activated LSCs with vigorous proliferation ability. Further studies confirm that CD63+CD56- quiescent LSCs express VEGFR2 and FGFR1, and they can be activated to proliferation and differentiation through combination of growth factors: VEGF-A and bFGF. These findings define an authentic adult liver stem cells compartment, make a further understanding of fate regulation on LSCs, and highlight its contribution to liver during pathophysiologic processes.

A novel efficient strategy to generate liver sinusoidal endothelial cells from human pluripotent stem cells

Liver sinusoidal endothelial cells (LSECs) are highly specialized endothelial cells (ECs) that play an important role in liver development and regeneration. Additionally, it is involved in various pathological processes, including steatosis, inflammation, fibrosis and hepatocellular carcinoma. However, the rapid dedifferentiation of LSECs after culture greatly limits their use in vitro modeling for biomedical applications. In this study, we developed a highly efficient protocol to induce LSEC-like cells from human induced pluripotent stem cells (hiPSCs) in only 8 days. Using single-cell transcriptomic analysis, we identified several novel LSEC-specific markers, such as EPAS1, LIFR, and NID1, as well as several previously revealed markers, such as CLEC4M, CLEC1B, CRHBP and FCN3. These LSEC markers are specifically expressed in our LSEC-like cells. Furthermore, hiPSC-derived cells expressed LSEC-specific proteins and exhibited LSEC-related functions, such as the uptake of acetylated low density lipoprotein (ac-LDL) and immune complex endocytosis. Overall, this study confirmed that our novel protocol allowed hiPSCs to rapidly acquire an LSEC-like phenotype and function in vitro. The ability to generate LSECs efficiently and rapidly may help to more precisely mimic liver development and disease progression in a liver-specific multicellular microenvironment, offering new insights into the development of novel therapeutic strategies.

Single-nucleus transcriptome reveals cell dynamic response of liver during the late chick embryonic development

The late embryonic development of the liver, a major metabolic organ, remains poorly characterized at single cell resolution. Here, we used single-nucleus RNA-sequencing (snRNA-seq) to characterize the chicken liver cells at 2 embryonic development time points (E14 and D1). We uncovered 8 cell types including hepatocytes, endothelial cells, hepatic stellate cells, erythrocytes, cholangiocytes, kupffer cells, mesothelial cells, and lymphocytes. And we discovered significant differences in the abundance of different cell types between E14 and D1. Moreover, we characterized the heterogeneity of hepatocytes, endothelial cells, and mesenchymal cells based on the gene regulatory networks of each clusters. Trajectory analyses revealed 128 genes associated with hepatocyte development and function, including apolipoprotein genes involved hepatic lipid metabolism and NADH dehydrogenase subunits involved hepatic oxidative phosphorylation. Furthermore, we identified the differentially expressed genes (DEGs) between E14 and D1 at the cellular levels, which contribute to changes in liver development and function. These DEGs were significantly enriched in PPAR signaling pathways and lipid metabolism related pathways. Our results presented the single-cell mapping of chick embryonic liver at late stages of development and demonstrated the metabolic changes across the 2 age stages at the cellular level, which can help to further study the molecular development mechanism of embryonic liver.

Helper Lipid-Enhanced mRNA Delivery for Treating Metabolic Dysfunction-Associated Fatty Liver Disease

Lipid nanoparticles (LNPs) represent the forefront of mRNA delivery platforms, yet achieving precise delivery to specific cells remains a challenge. The current targeting strategies complicate the formulation and impede the regulatory approval process. Here, through a straightforward regulation of helper lipids within LNPs, we introduce an engineered LNP designed for targeted delivery of mRNA into hepatocytes for metabolic dysfunction-associated fatty liver disease (MAFLD) treatment. The optimized LNP, supplied with POPC as the helper lipid, exhibits a 2.49-fold increase in mRNA transfection efficiency in hepatocytes compared to that of FDA-approved LNPs. CTP:phosphocholine cytidylyltransferase α mRNA is selected for delivery to hepatocytes through the optimized LNP system for self-calibration of phosphatidylcholine levels to prevent lipid droplet expansion in MAFLD. This strategy effectively regulates lipid homeostasis, while demonstrating proven biosafety. Our results present a mRNA therapy for MAFLD and open a new avenue for discovering potent lipids enabling mRNA delivery to specific cells.

Hepatocytes and the art of killing Plasmodium softly

The Plasmodium parasites that cause malaria undergo asymptomatic development in the parenchymal cells of the liver, the hepatocytes, prior to infecting erythrocytes and causing clinical disease. Traditionally, hepatocytes have been perceived as passive bystanders that allow hepatotropic pathogens such as Plasmodium to develop relatively unchallenged. However, now there is emerging evidence suggesting that hepatocytes can mount robust cell-autonomous immune responses that target Plasmodium, limiting its progression to the blood and reducing the incidence and severity of clinical malaria. Here we discuss our current understanding of hepatocyte cell-intrinsic immune responses that target Plasmodium and how these pathways impact malaria.

Multimodal decoding of human liver regeneration

The liver has a unique ability to regenerate1,2; however, in the setting of acute liver failure (ALF), this regenerative capacity is often overwhelmed, leaving emergency liver transplantation as the only curative option3-5. Here, to advance understanding of human liver regeneration, we use paired single-nucleus RNA sequencing combined with spatial profiling of healthy and ALF explant human livers to generate a single-cell, pan-lineage atlas of human liver regeneration. We uncover a novel ANXA2+ migratory hepatocyte subpopulation, which emerges during human liver regeneration, and a corollary subpopulation in a mouse model of acetaminophen (APAP)-induced liver regeneration. Interrogation of necrotic wound closure and hepatocyte proliferation across multiple timepoints following APAP-induced liver injury in mice demonstrates that wound closure precedes hepatocyte proliferation. Four-dimensional intravital imaging of APAP-induced mouse liver injury identifies motile hepatocytes at the edge of the necrotic area, enabling collective migration of the hepatocyte sheet to effect wound closure. Depletion of hepatocyte ANXA2 reduces hepatocyte growth factor-induced human and mouse hepatocyte migration in vitro, and abrogates necrotic wound closure following APAP-induced mouse liver injury. Together, our work dissects unanticipated aspects of liver regeneration, demonstrating an uncoupling of wound closure and hepatocyte proliferation and uncovering a novel migratory hepatocyte subpopulation that mediates wound closure following liver injury. Therapies designed to promote rapid reconstitution of normal hepatic microarchitecture and reparation of the gut-liver barrier may advance new areas of therapeutic discovery in regenerative medicine.

Angiocrine signaling in sinusoidal homeostasis and liver diseases

The hepatic sinusoids are composed of liver sinusoidal endothelial cells (LSECs), which are surrounded by hepatic stellate cells (HSCs) and contain liver-resident macrophages called Kupffer cells, and other patrolling immune cells. All these cells communicate with each other and with hepatocytes to maintain sinusoidal homeostasis and a spectrum of hepatic functions under healthy conditions. Sinusoidal homeostasis is disrupted by metabolites, toxins, viruses, and other pathological factors, leading to liver injury, chronic liver diseases, and cirrhosis. Alterations in hepatic sinusoids are linked to fibrosis progression and portal hypertension. LSECs are crucial regulators of cellular crosstalk within their microenvironment via angiocrine signaling. This review discusses the mechanisms by which angiocrine signaling orchestrates sinusoidal homeostasis, as well as the development of liver diseases. Here, we summarise the crosstalk between LSECs, HSCs, hepatocytes, cholangiocytes, and immune cells in health and disease and comment on potential novel therapeutic methods for treating liver diseases.

Mapping of mitogen and metabolic sensitivity in organoids defines requirements for human hepatocyte growth

Mechanisms underlying human hepatocyte growth in development and regeneration are incompletely understood. In vitro, human fetal hepatocytes (FH) can be robustly grown as organoids, while adult primary human hepatocyte (PHH) organoids remain difficult to expand, suggesting different growth requirements between fetal and adult hepatocytes. Here, we characterize hepatocyte organoid outgrowth using temporal transcriptomic and phenotypic approaches. FHs initiate reciprocal transcriptional programs involving increased proliferation and repressed lipid metabolism upon initiation of organoid growth. We exploit these insights to design maturation conditions for FH organoids, resulting in acquisition of mature hepatocyte morphological traits and increased expression of functional markers. During PHH organoid outgrowth in the same culture condition as for FHs, the adult transcriptomes initially mimic the fetal transcriptomic signatures, but PHHs rapidly acquire disbalanced proliferation-lipid metabolism dynamics, resulting in steatosis and halted organoid growth. IL6 supplementation, as emerged from the fetal dataset, and simultaneous activation of the metabolic regulator FXR, prevents steatosis and promotes PHH proliferation, resulting in improved expansion of the derived organoids. Single-cell RNA sequencing analyses reveal preservation of their fetal and adult hepatocyte identities in the respective organoid cultures. Our findings uncover mitogen requirements and metabolic differences determining proliferation of hepatocytes changing from development to adulthood.

Periportal macrophages protect against commensal-driven liver inflammation

The liver is the main gateway from the gut, and the unidirectional sinusoidal flow from portal to central veins constitutes heterogenous zones, including the periportal vein (PV) and the pericentral vein zones1-5. However, functional differences in the immune system in each zone remain poorly understood. Here intravital imaging revealed that inflammatory responses are suppressed in PV zones. Zone-specific single-cell transcriptomics detected a subset of immunosuppressive macrophages enriched in PV zones that express high levels of interleukin-10 and Marco, a scavenger receptor that sequesters pro-inflammatory pathogen-associated molecular patterns and damage-associated molecular patterns, and consequently suppress immune responses. Induction of Marco+ immunosuppressive macrophages depended on gut microbiota. In particular, a specific bacterial family, Odoribacteraceae, was identified to induce this macrophage subset through its postbiotic isoallolithocholic acid. Intestinal barrier leakage resulted in inflammation in PV zones, which was markedly augmented in Marco-deficient conditions. Chronic liver inflammatory diseases such as primary sclerosing cholangitis (PSC) and non-alcoholic steatohepatitis (NASH) showed decreased numbers of Marco+ macrophages. Functional ablation of Marco+ macrophages led to PSC-like inflammatory phenotypes related to colitis and exacerbated steatosis in NASH in animal experimental models. Collectively, commensal bacteria induce Marco+ immunosuppressive macrophages, which consequently limit excessive inflammation at the gateway of the liver. Failure of this self-limiting system promotes hepatic inflammatory disorders such as PSC and NASH.

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.

Liver specification of human iPSC-derived endothelial cells transplanted into mouse liver

Background & Aims

Liver sinusoidal endothelial cells (LSECs) are important in liver development, regeneration, and pathophysiology, but the differentiation process underlying their tissue-specific phenotype is poorly understood and difficult to study because primary human cells are scarce. The aim of this study was to use human induced pluripotent stem cell (hiPSC)-derived LSEC-like cells to investigate the differentiation process of LSECs.

Methods

hiPSC-derived endothelial cells (iECs) were transplanted into the livers of Fah-/-/Rag2-/-/Il2rg-/- mice and assessed over a 12-week period. Lineage tracing, immunofluorescence, flow cytometry, plasma human factor VIII measurement, and bulk and single cell transcriptomic analysis were used to assess the molecular and functional changes that occurred following transplantation.

Results

Progressive and long-term repopulation of the liver vasculature occurred as iECs expanded along the sinusoids between hepatocytes and increasingly produced human factor VIII, indicating differentiation into LSEC-like cells. To chart the developmental profile associated with LSEC specification, the bulk transcriptomes of transplanted cells between 1 and 12 weeks after transplantation were compared against primary human adult LSECs. This demonstrated a chronological increase in LSEC markers, LSEC differentiation pathways, and zonation. Bulk transcriptome analysis suggested that the transcription factors NOTCH1, GATA4, and FOS have a central role in LSEC specification, interacting with a network of 27 transcription factors. Novel markers associated with this process included EMCN and CLEC14A. Additionally, single cell transcriptomic analysis demonstrated that transplanted iECs at 4 weeks contained zonal subpopulations with a region-specific phenotype.

Conclusions

Collectively, this study confirms that hiPSCs can adopt LSEC-like features and provides insight into LSEC specification. This humanised xenograft system can be applied to further interrogate LSEC developmental biology and pathophysiology, bypassing current logistical obstacles associated with primary human LSECs.

Impact and implications

Liver sinusoidal endothelial cells (LSECs) are important cells for liver biology, but better model systems are required to study them. We present a pluripotent stem cell xenografting model that produces human LSEC-like cells. A detailed and longitudinal transcriptomic analysis of the development of LSEC-like cells is included, which will guide future studies to interrogate LSEC biology and produce LSEC-like cells that could be used for regenerative medicine.

A spatiotemporal atlas of cholestatic injury and repair in mice

Cholestatic liver injuries, characterized by regional damage around the bile ductular region, lack curative therapies and cause considerable mortality. Here we generated a high-definition spatiotemporal atlas of gene expression during cholestatic injury and repair in mice by integrating spatial enhanced resolution omics sequencing and single-cell transcriptomics. Spatiotemporal analyses revealed a key role of cholangiocyte-driven signaling correlating with the periportal damage-repair response. Cholangiocytes express genes related to recruitment and differentiation of lipid-associated macrophages, which generate feedback signals enhancing ductular reaction. Moreover, cholangiocytes express high TGFβ in association with the conversion of liver progenitor-like cells into cholangiocytes during injury and the dampened proliferation of periportal hepatocytes during recovery. Notably, Atoh8 restricts hepatocyte proliferation during 3,5-diethoxycarbonyl-1,4-dihydro-collidin damage and is quickly downregulated after injury withdrawal, allowing hepatocytes to respond to growth signals. Our findings lay a keystone for in-depth studies of cellular dynamics and molecular mechanisms of cholestatic injuries, which may further develop into therapies for cholangiopathies.

Liver-specific mitochondrial amidoxime-reducing component 1 (Mtarc1) knockdown protects the liver from diet-induced MASH in multiple mouse models

BACKGROUND

Human genetic studies have identified several mitochondrial amidoxime-reducing component 1 (MTARC1) variants as protective against metabolic dysfunction-associated steatotic liver disease. The MTARC1 variants are associated with decreased plasma lipids and liver enzymes and reduced liver-related mortality. However, the role of mARC1 in fatty liver disease is still unclear.

METHODS

Given that mARC1 is mainly expressed in hepatocytes, we developed an N-acetylgalactosamine-conjugated mouse Mtarc1 siRNA, applying it in multiple in vivo models to investigate the role of mARC1 using multiomic techniques.

RESULTS

In ob/ob mice, knockdown of Mtarc1 in mouse hepatocytes resulted in decreased serum liver enzymes, LDL-cholesterol, and liver triglycerides. Reduction of mARC1 also reduced liver weight, improved lipid profiles, and attenuated liver pathological changes in 2 diet-induced metabolic dysfunction-associated steatohepatitis mouse models. A comprehensive analysis of mARC1-deficient liver from a metabolic dysfunction-associated steatohepatitis mouse model by metabolomics, proteomics, and lipidomics showed that Mtarc1 knockdown partially restored metabolites and lipids altered by diet.

CONCLUSIONS

Taken together, reducing mARC1 expression in hepatocytes protects against metabolic dysfunction-associated steatohepatitis in multiple murine models, suggesting a potential therapeutic approach for this chronic liver disease.

Role of calcium integrin-binding protein 1 in the mechanobiology of the liver endothelium

Liver sinusoidal endothelial cells (LSECs) dysfunction is a key process in the development of chronic liver disease (CLD). Progressive scarring increases liver stiffness in a winch-like loop stimulating a dysfunctional liver cell phenotype. Cellular stretching is supported by biomechanically modulated molecular factors (BMMFs) that can translocate into the cytoplasm to support mechanotransduction through cytoskeleton remodeling and gene transcription. Currently, the molecular mechanisms of stiffness-induced LSECs dysfunction remain largely unclear. Here we propose calcium- and integrin-binding protein 1 (CIB1) as BMMF with crucial role in LSECs mechanobiology in CLD. CIB1 expression and translocation was characterized in healthy and cirrhotic human livers and in LSECs cultured on polyacrylamide gels with healthy and cirrhotic-like stiffnesses. Following the modulation of CIB1 with siRNA, the transcriptome was scrutinized to understand downstream effects of CIB1 downregulation. CIB1 expression is increased in LSECs in human cirrhosis. In vitro, CIB1 emerges as an endothelial BMMF. In human umbilical vein endothelial cells and LSECs, CIB1 expression and localization are modulated by stiffness-induced trafficking across the nuclear membrane. LSECs from cirrhotic liver tissue both in animal model and human disease exhibit an increased amount of CIB1 in cytoplasm. Knockdown of CIB1 in LSECs exposed to high stiffness improves LSECs phenotype by regulating the intracellular tension as well as the inflammatory response. Our results demonstrate that CIB1 is a key factor in sustaining cellular tension and stretching in response to high stiffness. CIB1 downregulation ameliorates LSECs dysfunction, enhancing their redifferentiation, and reducing the inflammatory response.

Gut-liver axis: Recent concepts in pathophysiology in alcohol-associated liver disease

The growing recognition of the role of the gut microbiome's impact on alcohol-associated diseases, especially in alcohol-associated liver disease, emphasizes the need to understand molecular mechanisms involved in governing organ-organ communication to identify novel avenues to combat alcohol-associated diseases. The gut-liver axis refers to the bidirectional communication and interaction between the gut and the liver. Intestinal microbiota plays a pivotal role in maintaining homeostasis within the gut-liver axis, and this axis plays a significant role in alcohol-associated liver disease. The intricate communication between intestine and liver involves communication between multiple cellular components in each organ that enable them to carry out their physiological functions. In this review, we focus on novel approaches to understanding how chronic alcohol exposure impacts the microbiome and individual cells within the liver and intestine, as well as the impact of ethanol on the molecular machinery required for intraorgan and interorgan communication.

Liver cancer development driven by the AP-1/c-Jun~Fra-2 dimer through c-Myc

Hepatocellular carcinoma (HCC) is a leading cause of cancer-related death. HCC incidence is on the rise, while treatment options remain limited. Thus, a better understanding of the molecular pathways involved in HCC development has become a priority to guide future therapies. While previous studies implicated the Activator Protein-1 (AP-1) (Fos/Jun) transcription factor family members c-Fos and c-Jun in HCC formation, the contribution of Fos-related antigens (Fra-) 1 and 2 is unknown. Here, we show that hepatocyte-restricted expression of a single chain c-Jun~Fra-2 protein, which functionally mimics the c-Jun/Fra-2 AP-1 dimer, results in spontaneous HCC formation in c-Jun~Fra-2hep mice. Several hallmarks of human HCC, such as cell cycle dysregulation and the expression of HCC markers are observed in liver tumors arising in c-Jun~Fra-2hep mice. Tumorigenesis occurs in the context of mild inflammation, low-grade fibrosis, and Pparγ-driven dyslipidemia. Subsequent analyses revealed increased expression of c-Myc, evidently under direct regulation by AP-1 through a conserved distal 3' enhancer. Importantly, c-Jun~Fra-2-induced tumors revert upon switching off transgene expression, suggesting oncogene addiction to the c-Jun~Fra-2 transgene. Tumors escaping reversion maintained c-Myc and c-Myc target gene expression, likely due to increased c-Fos. Interfering with c-Myc in established tumors using the Bromodomain and Extra-Terminal motif inhibitor JQ-1 diminished liver tumor growth in c-Jun~Fra-2 mutant mice. Thus, our data establish c-Jun~Fra-2hep mice as a model to study liver tumorigenesis and identify the c-Jun/Fra-2-Myc interaction as a potential target to improve HCC patient stratification and/or therapy.

The transcription factor ZNF469 regulates collagen production in liver fibrosis

Non-alcoholic fatty liver disease (NAFLD) - characterized by excess accumulation of fat in the liver - now affects one third of the world's population. As NAFLD progresses, extracellular matrix components including collagen accumulate in the liver causing tissue fibrosis, a major determinant of disease severity and mortality. To identify transcriptional regulators of fibrosis, we computationally inferred the activity of transcription factors (TFs) relevant to fibrosis by profiling the matched transcriptomes and epigenomes of 108 human liver biopsies from a deeply-characterized cohort of patients spanning the full histopathologic spectrum of NAFLD. CRISPR-based genetic knockout of the top 100 TFs identified ZNF469 as a regulator of collagen expression in primary human hepatic stellate cells (HSCs). Gain- and loss-of-function studies established that ZNF469 regulates collagen genes and genes involved in matrix homeostasis through direct binding to gene bodies and regulatory elements. By integrating multiomic large-scale profiling of human biopsies with extensive experimental validation we demonstrate that ZNF469 is a transcriptional regulator of collagen in HSCs. Overall, these data nominate ZNF469 as a previously unrecognized determinant of NAFLD-associated liver fibrosis.

Spatial genomics: mapping human steatotic liver disease

Metabolic dysfunction-associated steatotic liver disease (MASLD, formerly known as non-alcoholic fatty liver disease) is a leading cause of chronic liver disease worldwide. MASLD can progress to metabolic dysfunction-associated steatohepatitis (MASH, formerly known as non-alcoholic steatohepatitis) with subsequent liver cirrhosis and hepatocellular carcinoma formation. The advent of current technologies such as single-cell and single-nuclei RNA sequencing have transformed our understanding of the liver in homeostasis and disease. The next frontier is contextualizing this single-cell information in its native spatial orientation. This understanding will markedly accelerate discovery science in hepatology, resulting in a further step-change in our knowledge of liver biology and pathobiology. In this Review, we discuss up-to-date knowledge of MASLD development and progression and how the burgeoning field of spatial genomics is driving exciting new developments in our understanding of human liver disease pathogenesis and therapeutic target identification.

Bulk and single-cell RNA sequencing analysis with 101 machine learning combinations reveal neutrophil extracellular trap involvement in hepatic ischemia-reperfusion injury and early allograft dysfunction

BACKGROUND

Hepatic ischaemia-reperfusion injury (HIRI) is a major clinical concern during the perioperative period and is closely associated with early allograft dysfunction (EAD), acute rejection (AR) and long-term graft survival. Neutrophil extracellular traps (NETs) are extracellular structures formed by the release of decondensed chromatin and granular proteins following neutrophil stimulation. There is growing evidence that NETs are involved in the progression of various liver transplantation complications, including ischaemia-reperfusion injury (IRI). This study aimed to comprehensively analyse the expression patterns of NET-related genes (NRGs) in HIRI, identify HIRI subtypes with distinct characteristics, and develop a reliable EAD prediction model.

METHODS

Microarray, bulk RNA-seq, and single-cell sequencing datasets were obtained from the GEO database. Initially, differentially expressed NRGs (DE-NRGs) were identified using differential gene expression analyses. We then utilised a non-negative matrix factorisation (NMF) algorithm to classify HIRI samples. Subsequently, we employed machine learning algorithms to screen the hub NRGs related to EAD and developed an EAD prediction model based on these hub NRGs. Concurrently, we assessed the expression patterns of hub NRGs at the single-cell level using the HIRI. Additionally, we validated C5AR1 expression and its effect on HIRI and NETs formation in a rat orthotopic liver transplantation (OLT) model.

RESULTS

In this study, we identified 11 DE-NRGs in the HIRI context. Based on these 11 DE-NRGs, HIRI samples were classified into two distinct clusters. Cluster1 exhibited a low expression of DE-NRGs, minimal neutrophil infiltration, mild inflammation, and a low incidence of EAD. Conversely, Cluster2 displayed the opposite phenotype, with an activated inflammatory subtype and a higher incidence of EAD. Furthermore, an EAD prediction model was developed using the four hub NRGs associated with EAD. Based on risk scores, HIRI samples were classified into high- and low-risk groups. The OLT model confirmed substantial upregulation of C5AR1 expression in the liver tissue, accompanied by increased formation of NETs. Treatment with a C5AR1 antagonist improved liver function, reduced tissue inflammation, and decreased NETs formation.

CONCLUSIONS

This study distinguished two apparent HIRI subtypes, established a predictive model for EAD, and validated the effect of C5AR1 on HIRI. These findings provide novel perspectives for the development of advanced clinical strategies to enhance the outcomes of liver transplant recipients.

ER stress signaling at the interphase between MASH and HCC

HCC is the most frequent primary liver cancer with an extremely poor prognosis and often develops on preset of chronic liver diseases. Major risk factors for HCC include metabolic dysfunction-associated steatohepatitis, a complex multifactorial condition associated with abnormal endoplasmic reticulum (ER) proteostasis. To cope with ER stress, the unfolded protein response engages adaptive reactions to restore the secretory capacity of the cell. Recent advances revealed that ER stress signaling plays a critical role in HCC progression. Here, we propose that chronic ER stress is a common transversal factor contributing to the transition from liver disease (risk factor) to HCC. Interventional strategies to target the unfolded protein response in HCC, such as cancer therapy, are also discussed.

A Universal Strategy to Construct High-Performance Homo- and Heterogeneous Microgel Assembly Bioinks

Three dimensional (3D) extrusion bioprinting aims to replicate the complex architectures and functions of natural tissues and organs. However, the conventional hydrogel and new-emerging microgel bioinks are both difficult in achieving simultaneously high shape-fidelity and good maintenance of cell viability/function, leading to limited amount of qualified hydrogel/microgel bioinks. Herein, a universal strategy is reported to construct high-performance microgel assembly (MA) bioinks by using epigallocatechin gallate-modified hyaluronic acid (HA-EGCG) as coating agent and phenylboronic acid grafted hyaluronic acid (HA-PBA) as assembling agent. HA-EGCG can spontaneously form uniform coating on the microgel surface via mussel-inspired chemistry, while HA-PBA quickly forms dynamic phenylborate bonds with HA-EGCG, conferring the as-prepared MA bioinks with excellent rheological properties, self-healing, and tissue-adhesion. More importantly, this strategy is applicable to various microgel materials, enabling the preparation of homo- and heterogeneous MA (homo-MA and hetero-MA) bioinks and the hierarchical printing of complicated structures with high fidelity by integration of different microgels containing multiple materials/cells in spatial and compositional levels. It further demonstrates the printing of breast cancer organoid in vitro using homo-MA and hetero-MA bioinks and its preliminary application for drug testing. This universal strategy offers a new solution to construct high-performance bioinks for extrusion bioprinting.

Single-Cell RNA Sequencing in Organ and Cell Transplantation

Single-cell RNA sequencing is a high-throughput novel method that provides transcriptional profiling of individual cells within biological samples. This method typically uses microfluidics systems to uncover the complex intercellular communication networks and biological pathways buried within highly heterogeneous cell populations in tissues. One important application of this technology sits in the fields of organ and stem cell transplantation, where complications such as graft rejection and other post-transplantation life-threatening issues may occur. In this review, we first focus on research in which single-cell RNA sequencing is used to study the transcriptional profile of transplanted tissues. This technology enables the analysis of the donor and recipient cells and identifies cell types and states associated with transplant complications and pathologies. We also review the use of single-cell RNA sequencing in stem cell implantation. This method enables studying the heterogeneity of normal and pathological stem cells and the heterogeneity in cell populations. With their remarkably rapid pace, the single-cell RNA sequencing methodologies will potentially result in breakthroughs in clinical transplantation in the coming years.

Single-cell RNA-sequencing of virus-specific cellular immune responses in chronic hepatitis B patients

Chronic hepatitis B (CHB) is a major global health challenge. CHB can be controlled by antivirals but a therapeutic cure is lacking. CHB is characterized by limited HBV-specific T cell reactivity and functionality and expression of inhibitory receptors. The mechanisms driving these T cell phenotypes are only partially understood. Here, we created a single-cell RNA-sequencing dataset of HBV immune responses in patients to contribute to a better understanding of the dysregulated immunity. Blood samples of a well-defined cohort of 21 CHB and 10 healthy controls, including a subset of 5 matched liver biopsies, were collected. scRNA-seq data of total immune cells (55,825) plus sorted HBV-specific (1,963), non-naive (32,773) and PD1+ T cells (96,631) was generated using the 10X Genomics platform (186,123 cells) or the full-length Smart-seq2 protocol (1,069 cells). The shared transcript count matrices of single-cells serve as a valuable resource describing transcriptional changes underlying dysfunctional HBV-related T cell responses in blood and liver tissue and offers the opportunity to identify targets or biomarkers for HBV-related immune exhaustion.

Multimodal mass spectrometry imaging identifies cell-type-specific metabolic and lipidomic variation in the mammalian liver

Spatial single-cell omics provides a readout of biochemical processes. It is challenging to capture the transient lipidome/metabolome from cells in a native tissue environment. We employed water gas cluster ion beam secondary ion mass spectrometry imaging ([H2O]n>28K-GCIB-SIMS) at ≤3 μm resolution using a cryogenic imaging workflow. This allowed multiple biomolecular imaging modes on the near-native-state liver at single-cell resolution. Our workflow utilizes desorption electrospray ionization (DESI) to build a reference map of metabolic heterogeneity and zonation across liver functional units at tissue level. Cryogenic dual-SIMS integrated metabolomics, lipidomics, and proteomics in the same liver lobules at single-cell level, characterizing the cellular landscape and metabolic states in different cell types. Lipids and metabolites classified liver metabolic zones, cell types and subtypes, highlighting the power of spatial multi-omics at high spatial resolution for understanding celluar and biomolecular organizations in the mammalian liver.

A novel hypergraph model for identifying and prioritizing personalized drivers in cancer

Cancer development is driven by an accumulation of a small number of driver genetic mutations that confer the selective growth advantage to the cell, while most passenger mutations do not contribute to tumor progression. The identification of these driver genes responsible for tumorigenesis is a crucial step in designing effective cancer treatments. Although many computational methods have been developed with this purpose, the majority of existing methods solely provided a single driver gene list for the entire cohort of patients, ignoring the high heterogeneity of driver events across patients. It remains challenging to identify the personalized driver genes. Here, we propose a novel method (PDRWH), which aims to prioritize the mutated genes of a single patient based on their impact on the abnormal expression of downstream genes across a group of patients who share the co-mutation genes and similar gene expression profiles. The wide experimental results on 16 cancer datasets from TCGA showed that PDRWH excels in identifying known general driver genes and tumor-specific drivers. In the comparative testing across five cancer types, PDRWH outperformed existing individual-level methods as well as cohort-level methods. Our results also demonstrated that PDRWH could identify both common and rare drivers. The personalized driver profiles could improve tumor stratification, providing new insights into understanding tumor heterogeneity and taking a further step toward personalized treatment. We also validated one of our predicted novel personalized driver genes on tumor cell proliferation by vitro cell-based assays, the promoting effect of the high expression of Low-density lipoprotein receptor-related protein 1 (LRP1) on tumor cell proliferation.

Genetics of liver disease in adults

Chronic liver disease stands as a significant global health problem with an estimated 2 million annual deaths across the globe. Combining the use of next-generation sequencing technologies with evolving knowledge in the interpretation of genetic variation across the human genome is propelling our understanding, diagnosis, and management of both rare and common liver diseases. Here, we review the contribution of risk and protective alleles to common forms of liver disease, the rising number of monogenic diseases affecting the liver, and the role of somatic genetic variants in the onset and progression of oncological and non-oncological liver diseases. The incorporation of genomic information in the diagnosis and management of patients with liver disease is driving the beginning of a new era of genomics-informed clinical hepatology practice, facilitating personalized medicine, and improving patient care.

CSN6-SPOP-HMGCS1 Axis Promotes Hepatocellular Carcinoma Progression via YAP1 Activation

Cholesterol metabolism has important roles in maintaining membrane integrity and countering the development of diseases such as obesity and cancers. Cancer cells sustain cholesterol biogenesis for their proliferation and microenvironment reprograming even when sterols are abundant. However, efficacy of targeting cholesterol metabolism for cancer treatment is always compromised. Here it is shown that CSN6 is elevated in HCC and is a positive regulator of hydroxymethylglutaryl-CoA synthase 1 (HMGCS1) of mevalonate (MVA) pathway to promote tumorigenesis. Mechanistically, CSN6 antagonizes speckle-type POZ protein (SPOP) ubiquitin ligase to stabilize HMGCS1, which in turn activates YAP1 to promote tumor growth. In orthotopic liver cancer models, targeting CSN6 and HMGCS1 hinders tumor growth in both normal and high fat diet. Significantly, HMGCS1 depletion improves YAP inhibitor efficacy in patient derived xenograft models. The results identify a CSN6-HMGCS1-YAP1 axis mediating tumor outgrowth in HCC and propose a therapeutic strategy of targeting non-alcoholic fatty liver diseases- associated HCC.

Digital pathology and spatial omics in steatohepatitis: Clinical applications and discovery potentials

Steatohepatitis with diverse etiologies is the most common histological manifestation in patients with liver disease. However, there are currently no specific histopathological features pathognomonic for metabolic dysfunction-associated steatotic liver disease, alcohol-associated liver disease, or metabolic dysfunction-associated steatotic liver disease with increased alcohol intake. Digitizing traditional pathology slides has created an emerging field of digital pathology, allowing for easier access, storage, sharing, and analysis of whole-slide images. Artificial intelligence (AI) algorithms have been developed for whole-slide images to enhance the accuracy and speed of the histological interpretation of steatohepatitis and are currently employed in biomarker development. Spatial biology is a novel field that enables investigators to map gene and protein expression within a specific region of interest on liver histological sections, examine disease heterogeneity within tissues, and understand the relationship between molecular changes and distinct tissue morphology. Here, we review the utility of digital pathology (using linear and nonlinear microscopy) augmented with AI analysis to improve the accuracy of histological interpretation. We will also discuss the spatial omics landscape with special emphasis on the strengths and limitations of established spatial transcriptomics and proteomics technologies and their application in steatohepatitis. We then highlight the power of multimodal integration of digital pathology augmented by machine learning (ML)algorithms with spatial biology. The review concludes with a discussion of the current gaps in knowledge, the limitations and premises of these tools and technologies, and the areas of future research.

Integrating single-cell RNA sequencing data to genome-wide association analysis data identifies significant cell types in influenza A virus infection and COVID-19

With the global pandemic of COVID-19, the research on influenza virus has entered a new stage, but it is difficult to elucidate the pathogenesis of influenza disease. Genome-wide association studies (GWASs) have greatly shed light on the role of host genetic background in influenza pathogenesis and prognosis, whereas single-cell RNA sequencing (scRNA-seq) has enabled unprecedented resolution of cellular diversity and in vivo following influenza disease. Here, we performed a comprehensive analysis of influenza GWAS and scRNA-seq data to reveal cell types associated with influenza disease and provide clues to understanding pathogenesis. We downloaded two GWAS summary data, two scRNA-seq data on influenza disease. After defining cell types for each scRNA-seq data, we used RolyPoly and LDSC-cts to integrate GWAS and scRNA-seq. Furthermore, we analyzed scRNA-seq data from the peripheral blood mononuclear cells (PBMCs) of a healthy population to validate and compare our results. After processing the scRNA-seq data, we obtained approximately 70 000 cells and identified up to 13 cell types. For the European population analysis, we determined an association between neutrophils and influenza disease. For the East Asian population analysis, we identified an association between monocytes and influenza disease. In addition, we also identified monocytes as a significantly related cell type in a dataset of healthy human PBMCs. In this comprehensive analysis, we identified neutrophils and monocytes as influenza disease-associated cell types. More attention and validation should be given in future studies.

Dynamic YAP expression in the non-parenchymal liver cell compartment controls heterologous cell communication

INTRODUCTION

The Hippo pathway and its transcriptional effectors yes-associated protein (YAP) and transcriptional coactivator with PDZ-binding motif (TAZ) are targets for cancer therapy. It is important to determine if the activation of one factor compensates for the inhibition of the other. Moreover, it is unknown if YAP/TAZ-directed perturbation affects cell-cell communication of non-malignant liver cells.

MATERIALS AND METHODS

To investigate liver-specific phenotypes caused by YAP and TAZ inactivation, we generated mice with hepatocyte (HC) and biliary epithelial cell (BEC)-specific deletions for both factors (YAPKO, TAZKO and double knock-out (DKO)). Immunohistochemistry, single-cell sequencing, and proteomics were used to analyze liver tissues and serum.

RESULTS

The loss of BECs, liver fibrosis, and necrosis characterized livers from YAPKO and DKO mice. This phenotype was weakened in DKO tissues compared to specimens from YAPKO animals. After depletion of YAP in HCs and BECs, YAP expression was induced in non-parenchymal cells (NPCs) in a cholestasis-independent manner. YAP positivity was detected in subgroups of Kupffer cells (KCs) and endothelial cells (ECs). The secretion of pro-inflammatory chemokines and cytokines such as C-X-C motif chemokine ligand 11 (CXCL11), fms-related receptor tyrosine kinase 3 ligand (FLT3L), and soluble intercellular adhesion molecule-1 (ICAM1) was increased in the serum of YAPKO animals. YAP activation in NPCs could contribute to inflammation via TEA domain transcription factor (TEAD)-dependent transcriptional regulation of secreted factors.

CONCLUSION

YAP inactivation in HCs and BECs causes liver damage, and concomitant TAZ deletion does not enhance but reduces this phenotype. Additionally, we present a new mechanism by which YAP contributes to cell-cell communication originating from NPCs.

Single-cell transcriptome landscape of zebrafish liver reveals hepatocytes and immune cell interactions in understanding nonalcoholic fatty liver disease

Nonalcoholic fatty liver disease (NAFLD) is becoming the most common chronic liver disease in the world. Immunity is the major contributing factor in NAFLD; however, the interaction of immune cells and hepatocytes in disease progression has not been fully elucidated. As a popular species for studying NAFLD, zebrafish, whose liver is a complex immune system mediated by immune cells and non-immune cells in maintaining immune tolerance and homeostasis. Understanding the cellular composition and immune environment of zebrafish liver is of great significance for its application in NAFLD. Here, we established a liver atlas that consists of 10 cell types using single-cell RNA sequencing (scRNA-seq). By examining the heterogeneity of hepatocytes and analyzing the expression of NAFLD-associated genes in the specific cluster, we provide a potential target cell model to study NAFLD. Additionally, our analysis identified two subtypes of distinct resident macrophages with inflammatory and non-inflammatory functions and characterized the successive stepwise development of T cell subclusters in the liver. Importantly, we uncovered the possible regulation of macrophages and T cells on target cells of fatty liver by analyzing the cellular interaction between hepatocytes and immune cells. Our data provide valuable information for an in-depth study of immune cells targeting hepatocytes to regulate the immune balance in NAFLD.

Oxygen-dependent histone lysine demethylase 4 restricts hepatitis B virus replication

Mammalian cells have evolved strategies to regulate gene expression when oxygen is limited. Hypoxia-inducible factors (HIF) are the major transcriptional regulators of host gene expression. We previously reported that HIFs bind and activate hepatitis B virus (HBV) DNA transcription under low oxygen conditions; however, the global cellular response to low oxygen is mediated by a family of oxygenases that work in concert with HIFs. Recent studies have identified a role for chromatin modifiers in sensing cellular oxygen and orchestrating transcriptional responses, but their role in the HBV life cycle is as yet undefined. We demonstrated that histone lysine demethylase 4 (KDM4) can restrict HBV, and pharmacological or oxygen-mediated inhibition of the demethylase increases viral RNAs derived from both episomal and integrated copies of the viral genome. Sequencing studies demonstrated that KDM4 is a major regulator of the hepatic transcriptome, which defines hepatocellular permissivity to HBV infection. We propose a model where HBV exploits cellular oxygen sensors to replicate and persist in the liver. Understanding oxygen-dependent pathways that regulate HBV infection will facilitate the development of physiologically relevant cell-based models that support efficient HBV replication.

Biliary stem cells in health and cholangiopathies and cholangiocarcinoma

PURPOSE OF REVIEW

This review discusses evidence regarding progenitor populations of the biliary tree in the tissue regeneration and homeostasis, and the pathobiology of cholangiopathies and malignancies.

RECENT FINDINGS

In embryogenesis biliary multipotent progenitor subpopulation contributes cells not only to the pancreas and gall bladder but also to the liver. Cells equipped with a constellation of markers suggestive of the primitive endodermal phenotype exist in the peribiliary glands, the bile duct glands, of the intra- and extrahepatic bile ducts. These cells are able to be isolated and cultured easily, which demonstrates the persistence of a stable phenotype during in vitro expansion, the ability to self-renew in vitro, and the ability to differentiate between hepatocyte and biliary and pancreatic islet fates.

SUMMARY

In normal human livers, stem/progenitors cells are mostly restricted in two distinct niches, which are the bile ductules/canals of Hering and the peribiliary glands (PBGs) present inside the wall of large intrahepatic bile ducts. The existence of a network of stem/progenitor cell niches within the liver and along the entire biliary tree inform a patho-biological-based translational approach to biliary diseases and cholangiocarcinoma since it poses the basis to understand biliary regeneration after extensive or chronic injuries and progression to fibrosis and cancer.

Defining the T cell transcriptional landscape in pediatric liver transplant rejection at single cell resolution

Acute cellular rejection (ACR) affects >80% of pediatric liver transplant recipients within 5 years, and late ACR is associated with graft failure. Traditional anti-rejection therapy for late ACR is ineffective and has remained unchanged for six decades. Although CD8+ T cells promote late ACR, little has been done to define their specificity and gene expression. Here, we used single-cell sequencing and immune repertoire profiling (10X Genomics) on 30 cryopreserved 16G liver biopsies from 14 patients (5 pre-transplant or with no ACR, 9 with ACR). We identified expanded intragraft CD8+ T cell clonotypes (CD8EXP) and their gene expression profiles in response to anti-rejection treatment. Notably, we found that expanded CD8+ clonotypes (CD8EXP) bore markers of effector and CD56hiCD161- 'NK-like' T cells, retaining their clonotype identity and phenotype in subsequent biopsies from the same patients despite histologic ACR resolution. CD8EXP clonotypes localized to portal infiltrates during active ACR, and persisted in the lobule after histologic ACR resolution. CellPhoneDB analysis revealed differential crosstalk between KC and CD8EXP during late ACR, with activation of the LTB-LTBR pathway and downregulation of TGFß signaling. Therefore, persistently-detected intragraft CD8EXP clones remain active despite ACR treatment and may contribute to long-term allograft fibrosis and failure of operational tolerance.