A single-cell transcriptomic analysis reveals precise pathways and regulatory mechanisms underlying hepatoblast differentiation

OrganogenesisDB: A Comprehensive Database Exploring the Cell-Type Identities and Gene Expression Dynamics during Organogenesis

Organogenesis, the phase of embryonic development that starts at the end of gastrulation and continues until birth is the critical process for understanding cellular differentiation and maturation during organ development. The rapid development of single-cell transcriptomics technology has led to many novel discoveries in understanding organogenesis while also accumulating a large quantity of data. To fill this gap, OrganogenesisDB (http://organogenesisdb.com/), which is a comprehensive database dedicated to exploring cell-type identification and gene expression dynamics during organogenesis, is developed. OrganogenesisDB contains single-cell RNA sequencing data for more than 1.4 million cells from 49 published datasets spanning various developmental stages. Additionally, 3324 cell markers are manually curated for 1120 cell types across 9 human organs and 4 mouse organs. OrganogenesisDB leverages various analysis tools to assist users in annotating and understanding cell types at different developmental stages and helps in mining and presenting genes that exhibit specific patterns and play key regulatory roles during cell maturation and differentiation. This work provides a critical resource and useful tool for deciphering cell lineage determination and uncovering the mechanisms underlying organogenesis.

TrajectoryGeometry suggests cell fate decisions can involve branches rather than bifurcations

Differentiation of multipotential progenitor cells is a key process in the development of any multi-cellular organism and often continues throughout its life. It is often assumed that a bi-potential progenitor develops along a (relatively) straight trajectory until it reaches a decision point where the trajectory bifurcates. At this point one of two directions is chosen, each direction representing the unfolding of a new transcriptional programme. However, we have lacked quantitative means for testing this model. Accordingly, we have developed the R package TrajectoryGeometry. Applying this to published data we find several examples where, rather than bifurcate, developmental pathways branch. That is, the bipotential progenitor develops along a relatively straight trajectory leading to one of its potential fates. A second relatively straight trajectory branches off from this towards the other potential fate. In this sense only cells that branch off to follow the second trajectory make a 'decision'. Our methods give precise descriptions of the genes and cellular pathways involved in these trajectories. We speculate that branching may be the more common behaviour and may have advantages from a control-theoretic viewpoint.

An increment of diversity method for cell state trajectory inference of time-series scRNA-seq data

The increasing emergence of the time-series single-cell RNA sequencing (scRNA-seq) data, inferring developmental trajectory by connecting transcriptome similar cell states (i.e., cell types or clusters) has become a major challenge. Most existing computational methods are designed for individual cells and do not take into account the available time series information. We present IDTI based on the Increment of Diversity for Trajectory Inference, which combines time series information and the minimum increment of diversity method to infer cell state trajectory of time-series scRNA-seq data. We apply IDTI to simulated and three real diverse tissue development datasets, and compare it with six other commonly used trajectory inference methods in terms of topology similarity and branching accuracy. The results have shown that the IDTI method accurately constructs the cell state trajectory without the requirement of starting cells. In the performance test, we further demonstrate that IDTI has the advantages of high accuracy and strong robustness.

Van Gogh-like 2 is essential for the architectural patterning of the mammalian biliary tree

BACKGROUND & AIMS

In the developing liver, bipotent epithelial progenitor cells undergo lineage segregation to form hepatocytes, which constitute the bulk of the liver parenchyma, and biliary epithelial cells (cholangiocytes), which comprise the bile duct (a complex tubular network that is critical for normal liver function). Notch and TGFβ signalling promote the formation of a sheet of biliary epithelial cells, the ductal plate, that organises into discontinuous tubular structures. How these structures elongate and connect to form a continuous duct remains undefined. We aimed to define the mechanisms by which the ductal plate transitions from a simple sheet of epithelial cells into a complex and connected bile duct.

METHODS

By combining single-cell RNA sequencing of embryonic mouse livers with genetic tools and organoid models we functionally dissected the role of planar cell polarity in duct patterning.

RESULTS

We show that the planar cell polarity protein VANGL2 is expressed late in intrahepatic bile duct development and patterns the formation of cell-cell contacts between biliary cells. The patterning of these cell contacts regulates the normal polarisation of the actin cytoskeleton within biliary cells and loss of Vangl2 function results in the abnormal distribution of cortical actin remodelling, leading to the failure of bile duct formation.

CONCLUSIONS

Planar cell polarity is a critical step in the post-specification sculpture of the bile duct and is essential for establishing normal tissue architecture.

IMPACT AND IMPLICATIONS

Like other branched tissues, such as the lung and kidney, the bile ducts use planar cell polarity signalling to coordinate cell movements; however, how these biochemical signals are linked to ductular patterning remains unclear. Here we show that the core planar cell polarity protein VANGL2 patterns how cell-cell contacts form in the mammalian bile duct and how ductular cells transmit confluent mechanical changes along the length of a duct. This work sheds light on how biological tubes are patterned across mammalian tissues (including within the liver) and will be important in how we promote ductular growth in patients where the duct is mis-patterned or poorly formed.

Optimal sequencing budget allocation for trajectory reconstruction of single cells

BACKGROUND

Charting cellular trajectories over gene expression is key to understanding dynamic cellular processes and their underlying mechanisms. While advances in single-cell RNA-sequencing technologies and computational methods have pushed forward the recovery of such trajectories, trajectory inference remains a challenge due to the noisy, sparse, and high-dimensional nature of single-cell data. This challenge can be alleviated by increasing either the number of cells sampled along the trajectory (breadth) or the sequencing depth, i.e. the number of reads captured per cell (depth). Generally, these two factors are coupled due to an inherent breadth-depth tradeoff that arises when the sequencing budget is constrained due to financial or technical limitations.

RESULTS

Here we study the optimal allocation of a fixed sequencing budget to optimize the recovery of trajectory attributes. Empirical results reveal that reconstruction accuracy of internal cell structure in expression space scales with the logarithm of either the breadth or depth of sequencing. We additionally observe a power law relationship between the optimal number of sampled cells and the corresponding sequencing budget. For linear trajectories, non-monotonicity in trajectory reconstruction across the breadth-depth tradeoff can impact downstream inference, such as expression pattern analysis along the trajectory. We demonstrate these results for five single-cell RNA-sequencing datasets encompassing differentiation of embryonic stem cells, pancreatic beta cells, hepatoblast and multipotent hematopoietic cells, as well as induced reprogramming of embryonic fibroblasts into neurons. By addressing the challenges of single-cell data, our study offers insights into maximizing the efficiency of cellular trajectory analysis through strategic allocation of sequencing resources.

SlowMoMan: a web app for discovery of important features along user-drawn trajectories in 2D embeddings

Motivation

Nonlinear low-dimensional embeddings allow humans to visualize high-dimensional data, as is often seen in bioinformatics, where datasets may have tens of thousands of dimensions. However, relating the axes of a nonlinear embedding to the original dimensions is a nontrivial problem. In particular, humans may identify patterns or interesting subsections in the embedding, but cannot easily identify what those patterns correspond to in the original data.

Results

Thus, we present SlowMoMan (SLOW Motions on MANifolds), a web application which allows the user to draw a one-dimensional path onto a 2D embedding. Then, by back-projecting the manifold to the original, high-dimensional space, we sort the original features such that those most discriminative along the manifold are ranked highly. We show a number of pertinent use cases for our tool, including trajectory inference, spatial transcriptomics, and automatic cell classification.

Availability and implementation

Software: https://yunwilliamyu.github.io/SlowMoMan/; Code: https://github.com/yunwilliamyu/SlowMoMan.

Generation of in vivo-like multicellular liver organoids by mimicking developmental processes: A review

Liver is involved in metabolic reactions, ammonia detoxification, and immunity. Multicellular liver tissue cultures are more desirable for drug screening, disease modeling, and researching transplantation therapy, than hepatocytes monocultures. Hepatocytes monocultures are not stable for long. Further, hepatocyte-like cells induced from pluripotent stem cells and in vivo hepatocytes are functionally dissimilar. Organoid technology circumvents these issues by generating functional ex vivo liver tissue from intrinsic liver progenitor cells and extrinsic stem cells, including pluripotent stem cells. To function as in vivo liver tissue, the liver organoid cells must be arranged precisely in the 3-dimensional space, closely mimicking in vivo liver tissue. Moreover, for long term functioning, liver organoids must be appropriately vascularized and in contact with neighboring epithelial tissues (e.g., bile canaliculi and intrahepatic bile duct, or intrahepatic and extrahepatic bile ducts). Recent discoveries in liver developmental biology allows one to successfully induce liver component cells and generate organoids. Thus, here, in this review, we summarize the current state of knowledge on liver development with a focus on its application in generating different liver organoids. We also cover the future prospects in creating (functionally and structurally) in vivo-like liver organoids using the current knowledge on liver development.

Generation and characterization of mature hepatocyte organoids for liver metabolic studies

Hepatocyte organoids (HOs) generated in vitro are powerful tools for liver regeneration. However, previously reported HOs have mostly been fetal in nature with low expression levels of metabolic genes characteristic of adult liver functions, hampering their application in studies of metabolic regulation and therapeutic testing for liver disorders. Here, we report development of novel culture conditions that combine optimized levels of triiodothyronine (T3) with the removal of growth factors to enable successful generation of mature hepatocyte organoids (MHOs) of both mouse and human origin with metabolic functions characteristic of adult livers. We show that the MHOs can be used to study various metabolic functions including bile and urea production, zonal metabolic gene expression, and metabolic alterations in both alcoholic liver disease and non-alcoholic fatty liver disease, as well as hepatocyte proliferation, injury and cell fate changes. Notably, MHOs derived from human fetal hepatocytes also show improved hepatitis B virus infection. Therefore, these MHOs provide a powerful in vitro model for studies of human liver physiology and diseases. The human MHOs are potentially also a robust research tool for therapeutic development.

A Gold Standard-Derived Modular Barcoding Approach to Cancer Transcriptomics

A challenge with studying cancer transcriptomes is in distilling the wealth of information down into manageable portions of information. In this resource, we develop an approach that creates and assembles cancer type-specific gene expression modules into flexible barcodes, allowing for adaptation to a wide variety of uses. Specifically, we propose that modules derived organically from high-quality gold standards such as The Cancer Genome Atlas (TCGA) can accurately capture and describe functionally related genes that are relevant to specific cancer types. We show that such modules can: (1) uncover novel gene relationships and nominate new functional memberships, (2) improve and speed up analysis of smaller or lower-resolution datasets, (3) re-create and expand known cancer subtyping schemes, (4) act as a "decoder" to bridge seemingly disparate established gene signatures, and (5) efficiently apply single-cell RNA sequencing information to other datasets. Moreover, such modules can be used in conjunction with native spreadsheet program commands to create a powerful and rapid approach to hypothesis generation and testing that is readily accessible to non-bioinformaticians. Finally, we provide tools for users to create and interpret their own modules. Overall, the flexible modular nature of the proposed barcoding provides a user-friendly approach to rapidly decoding transcriptome-wide data for research or, potentially, clinical uses.

Hepatocyte regeneration is driven by embryo-like DNA methylation reprogramming

As a result of partial hepatectomy, the remaining liver tissue undergoes a process of renewed proliferation that leads to rapid regeneration of the liver. By following the early stages of this process, we observed dramatic programmed changes in the DNA methylation profile, characterized by both de novo and demethylation events, with a subsequent return to the original adult pattern as the liver matures. Strikingly, these transient alterations partially mimic the DNA methylation state of embryonic hepatoblasts (E16.5), indicating that hepatocytes actually undergo epigenetic dedifferentiation. Furthermore, Tet2/Tet3-deletion experiments demonstrated that these changes in methylation are necessary for carrying out basic embryonic functions, such as proliferation, a key step in liver regeneration. This implies that unlike tissue-specific regulatory regions that remain demethylated in the adult, early embryonic genes are programmed to first undergo demethylation, followed by remethylation as development proceeds. The identification of this built-in system may open targeting opportunities for regenerative medicine.

Evaluation of different media compositions promoting hepatocyte differentiation in the canine liver organoid model

Organoids are 3-dimensional (3D) self-assembled structures capable of replicating the microanatomy and physiology of the epithelial components of their organ of origin. Adult stem cell (ASC) derived organoids from the liver have previously been shown to differentiate into primarily mature cholangiocytes, and their partial differentiation into functional hepatocytes can be promoted using specific media compositions. While full morphological differentiation of mature hepatocytes from ASCs has not yet been reported for any species, the functional differentiation can be approximated using various media compositions. Six differentiation media formulations from published studies on hepatic organoids were used for the differentiation protocol. Target species for these protocols were humans, mice, cats, and dogs, and encompassed various combinations and concentrations of four major hepatocyte media components: Bone morphogenetic protein 7 (BMP7), Fibroblast Growth Factor 19 (FGF19), Dexamethasone (Dex), and Gamma-Secretase Inhibitor IX (DAPT). Additionally, removing R-spondin from basic organoid media has previously been shown to drive the differentiation of ASC into mature hepatocytes. Differentiation media (N = 20) were designed to encompass combinations of the four major hepatocyte media components. The preferred differentiation of ASC-derived organoids from liver tissue into mature hepatocytes over cholangiocytes was confirmed by albumin production in the culture supernatant. Out of the twenty media compositions tested, six media resulted in the production of the highest amounts of albumin in the supernatant of the organoids. The cell lines cultured using these six media were further characterized via histological staining, transmission electron microscopy, RNA in situ hybridization, analysis of gene expression patterns, immunofluorescence, and label-free proteomics. The results indicate that preferential hepatocyte maturation from canine ADC-derived organoids from liver tissue is mainly driven by Dexamethasone and DAPT components. FGF19 did not enhance organoid differentiation but improved cell culture survival. Furthermore, we confirm that removing R-spondin from the media is crucial for establishing mature hepatic organoid cultures.

Powerful and accurate detection of temporal gene expression patterns from multi-sample multi-stage single-cell transcriptomics data with TDEseq

We present a non-parametric statistical method called TDEseq that takes full advantage of smoothing splines basis functions to account for the dependence of multiple time points in scRNA-seq studies, and uses hierarchical structure linear additive mixed models to model the correlated cells within an individual. As a result, TDEseq demonstrates powerful performance in identifying four potential temporal expression patterns within a specific cell type. Extensive simulation studies and the analysis of four published scRNA-seq datasets show that TDEseq can produce well-calibrated p-values and up to 20% power gain over the existing methods for detecting temporal gene expression patterns.

Endothelial Slc35a1 Deficiency Causes Loss of LSEC Identity and Exacerbates Neonatal Lipid Deposition in the Liver in Mice

BACKGROUND & AIMS

The functional maturation of the liver largely occurs after birth. In the early stages of life, the liver of a newborn encounters enormous high-fat metabolic stress caused by the consumption of breast milk. It is unclear how the maturing liver adapts to high lipid metabolism. Liver sinusoidal endothelial cells (LSECs) play a fundamental role in establishing liver vasculature and are decorated with many glycoproteins on their surface. The Slc35a1 gene encodes a cytidine-5'-monophosphate (CMP)-sialic acid transporter responsible for transporting CMP-sialic acids between the cytoplasm and the Golgi apparatus for protein sialylation. This study aimed to determine whether endothelial sialylation plays a role in hepatic vasculogenesis and functional maturation.

METHODS

Endothelial-specific Slc35a1 knockout mice were generated. Liver tissues were collected for histologic analysis, lipidomic profiling, RNA sequencing, confocal immunofluorescence, and immunoblot analyses.

RESULTS

Endothelial Slc35a1-deficient mice exhibited excessive neonatal hepatic lipid deposition, severe liver damage, and high mortality. Endothelial deletion of Slc35a1 led to sinusoidal capillarization and disrupted hepatic zonation. Mechanistically, vascular endothelial growth factor receptor 2 (VEGFR2) in LSECs was desialylated and VEGFR2 signaling was enhanced in Slc35a1-deficient mice. Inhibition of VEGFR2 signaling by SU5416 alleviated lipid deposition and restored hepatic vasculature in Slc35a1-deficient mice.

CONCLUSIONS

Our findings suggest that sialylation of LSECs is critical for maintaining hepatic vascular development and lipid homeostasis. Targeting VEGFR2 signaling may be a new strategy to prevent liver disorders associated with abnormal vasculature and lipid deposition.

Interrogations of single-cell RNA splicing landscapes with SCASL define new cell identities with physiological relevance

RNA splicing shapes the gene regulatory programs that underlie various physiological and disease processes. Here, we present the SCASL (single-cell clustering based on alternative splicing landscapes) method for interrogating the heterogeneity of RNA splicing with single-cell RNA-seq data. SCASL resolves the issue of biased and sparse data coverage on single-cell RNA splicing and provides a new scheme for classifications of cell identities. With previously published datasets as examples, SCASL identifies new cell clusters indicating potentially precancerous and early-tumor stages in triple-negative breast cancer, illustrates cell lineages of embryonic liver development, and provides fine clusters of highly heterogeneous tumor-associated CD4 and CD8 T cells with functional and physiological relevance. Most of these findings are not readily available via conventional cell clustering based on single-cell gene expression data. Our study shows the potential of SCASL in revealing the intrinsic RNA splicing heterogeneity and generating biological insights into the dynamic and functional cell landscapes in complex tissues.

Sox9 links biliary maturation to branching morphogenesis

Branching morphogenesis couples cellular differentiation with development of tissue architecture. Intrahepatic bile duct (IHBD) morphogenesis is initiated with biliary epithelial cell (BEC) specification and eventually forms a heterogeneous network of large ducts and small ductules. Here, we show that Sox9 is required for developmental establishment of small ductules. IHBDs emerge as a webbed structure by E15.5 and undergo morphological maturation through 2 weeks of age. Developmental knockout of Sox9 leads to decreased postnatal branching morphogenesis, manifesting as loss of ductules in adult livers. In the absence of Sox9, BECs fail to mature and exhibit elevated TGF-β signaling and Activin A. Activin A induces developmental gene expression and morphological defects in BEC organoids and represses ductule formation in postnatal livers. Our data demonstrate that adult IHBD morphology and BEC maturation is regulated by the Sox9-dependent formation of precursors to ductules during development, mediated in part by downregulation of Activin A.

VEGFA mRNA-LNP promotes biliary epithelial cell-to-hepatocyte conversion in acute and chronic liver diseases and reverses steatosis and fibrosis

The liver is known for its remarkable regenerative ability through proliferation of hepatocytes. Yet, during chronic injury or severe hepatocyte death, proliferation of hepatocytes is exhausted. To overcome this hurdle, we propose vascular-endothelial-growth-factor A (VEGFA) as a therapeutic means to accelerate biliary epithelial-cell (BEC)-to-hepatocyte conversion. Investigation in zebrafish establishes that blocking VEGF receptors abrogates BEC-driven liver repair, while VEGFA overexpression promotes it. Delivery of VEGFA via nonintegrative and safe nucleoside-modified mRNA encapsulated into lipid nanoparticles (mRNA-LNPs) in acutely or chronically injured mouse livers induces robust BEC-to-hepatocyte conversion and elimination of steatosis and fibrosis. In human and murine diseased livers, we further identified VEGFA-receptor KDR-expressing BECs associated with KDR-expressing cell-derived hepatocytes. This work defines KDR-expressing cells, most likely being BECs, as facultative progenitors. This study reveals unexpected therapeutic benefits of VEGFA delivered via nucleoside-modified mRNA-LNP, whose safety is widely validated with COVID-19 vaccines, for harnessing BEC-driven repair to potentially treat liver diseases.

Chronic metabolic stress drives developmental programs and loss of tissue functions in non-transformed liver that mirror tumor states and stratify survival

Under chronic stress, cells must balance competing demands between cellular survival and tissue function. In metabolic dysfunction-associated steatotic liver disease (MASLD, formerly NAFLD/NASH), hepatocytes cooperate with structural and immune cells to perform crucial metabolic, synthetic, and detoxification functions despite nutrient imbalances. While prior work has emphasized stress-induced drivers of cell death, the dynamic adaptations of surviving cells and their functional repercussions remain unclear. Namely, we do not know which pathways and programs define cellular responses, what regulatory factors mediate (mal)adaptations, and how this aberrant activity connects to tissue-scale dysfunction and long-term disease outcomes. Here, by applying longitudinal single-cell multi -omics to a mouse model of chronic metabolic stress and extending to human cohorts, we show that stress drives survival-linked tradeoffs and metabolic rewiring, manifesting as shifts towards development-associated states in non-transformed hepatocytes with accompanying decreases in their professional functionality. Diet-induced adaptations occur significantly prior to tumorigenesis but parallel tumorigenesis-induced phenotypes and predict worsened human cancer survival. Through the development of a multi -omic computational gene regulatory inference framework and human in vitro and mouse in vivo genetic perturbations, we validate transcriptional (RELB, SOX4) and metabolic (HMGCS2) mediators that co-regulate and couple the balance between developmental state and hepatocyte functional identity programming. Our work defines cellular features of liver adaptation to chronic stress as well as their links to long-term disease outcomes and cancer hallmarks, unifying diverse axes of cellular dysfunction around core causal mechanisms.

Single-cell and spatial architecture of primary liver cancer

Primary liver cancer (PLC) poses a leading threat to human health, and its treatment options are limited. Meanwhile, the investigation of homogeneity and heterogeneity among PLCs remains challenging. Here, using single-cell RNA sequencing, spatial transcriptomic and bulk multi-omics, we elaborated a molecular architecture of 3 PLC types, namely hepatocellular carcinoma (HCC), intrahepatic cholangiocarcinoma (ICC) and combined hepatocellular-cholangiocarcinoma (CHC). Taking a high-resolution perspective, our observations revealed that CHC cells exhibit internally discordant phenotypes, whereas ICC and HCC exhibit distinct tumor-specific features. Specifically, ICC was found to be the primary source of cancer-associated fibroblasts, while HCC exhibited disrupted metabolism and greater individual heterogeneity of T cells. We further revealed a diversity of intermediate-state cells residing in the tumor-peritumor junctional zone, including a congregation of CPE+ intermediate-state endothelial cells (ECs), which harbored the molecular characteristics of tumor-associated ECs and normal ECs. This architecture offers insights into molecular characteristics of PLC microenvironment, and hints that the tumor-peritumor junctional zone could serve as a targeted region for precise therapeutical strategies.

In-organoid single-cell CRISPR screening reveals determinants of hepatocyte differentiation and maturation

BACKGROUND

Harnessing hepatocytes for basic research and regenerative medicine demands a complete understanding of the genetic determinants underlying hepatocyte differentiation and maturation. Single-cell CRISPR screens in organoids could link genetic perturbations with parallel transcriptomic readout in single cells, providing a powerful method to delineate roles of cell fate regulators. However, a big challenge for identifying key regulators during data analysis is the low expression levels of transcription factors (TFs), which are difficult to accurately estimate due to noise and dropouts in single-cell sequencing. Also, it is often the changes in TF activities in the transcriptional cascade rather than the expression levels of TFs that are relevant to the cell fate transition.

RESULTS

Here, we develop Organoid-based Single-cell CRISPR screening Analyzed with Regulons (OSCAR), a framework using regulon activities as readouts to dissect gene knockout effects in organoids. In adult-stem-cell-derived liver organoids, we map transcriptomes in 80,576 cells upon 246 perturbations associated with transcriptional regulation of hepatocyte formation. Using OSCAR, we identify known and novel positive and negative regulators, among which Fos and Ubr5 are the top-ranked ones. Further single-gene loss-of-function assays demonstrate that Fos depletion in mouse and human liver organoids promote hepatocyte differentiation by specific upregulation of liver metabolic genes and pathways, and conditional knockout of Ubr5 in mouse liver delays hepatocyte maturation.

CONCLUSIONS

Altogether, we provide a framework to explore lineage specifiers in a rapid and systematic manner, and identify hepatocyte determinators with potential clinical applications.

Determination of key events in mouse hepatocyte maturation at the single-cell level

Hepatocytes, the liver's predominant cells, perform numerous essential biological functions. However, crucial events and regulators during hepatocyte maturation require in-depth investigation. In this study, we performed single-cell RNA sequencing (scRNA-seq) and single-nucleus RNA sequencing (snRNA-seq) to explore the precise hepatocyte development process in mice. We defined three maturation stages of postnatal hepatocytes, each of which establishes specific metabolic functions and exhibits distinct proliferation rates. Hepatic zonation is gradually formed during hepatocyte maturation. Hepatocytes or their nuclei with distinct ploidies exhibit zonation preferences in distribution and asynchrony in maturation. Moreover, by combining gene regulatory network analysis with in vivo genetic manipulation, we identified critical maturation- and zonation-related transcription factors. This study not only delineates the comprehensive transcriptomic profiles of hepatocyte maturation but also presents a paradigm to identify genes that function in the development of hepatocyte maturation and zonation by combining genetic manipulation and measurement of coordinates in a single-cell developmental trajectory.

Lineage tracing identifies heterogeneous hepatoblast contribution to cell lineages and postembryonic organ growth dynamics

To meet the physiological demands of the body, organs need to establish a functional tissue architecture and adequate size as the embryo develops to adulthood. In the liver, uni- and bipotent progenitor differentiation into hepatocytes and biliary epithelial cells (BECs), and their relative proportions, comprise the functional architecture. Yet, the contribution of individual liver progenitors at the organ level to both fates, and their specific proportion, is unresolved. Combining mathematical modelling with organ-wide, multispectral FRaeppli-NLS lineage tracing in zebrafish, we demonstrate that a precise BEC-to-hepatocyte ratio is established (i) fast, (ii) solely by heterogeneous lineage decisions from uni- and bipotent progenitors, and (iii) independent of subsequent cell type-specific proliferation. Extending lineage tracing to adulthood determined that embryonic cells undergo spatially heterogeneous three-dimensional growth associated with distinct environments. Strikingly, giant clusters comprising almost half a ventral lobe suggest lobe-specific dominant-like growth behaviours. We show substantial hepatocyte polyploidy in juveniles representing another hallmark of postembryonic liver growth. Our findings uncover heterogeneous progenitor contributions to tissue architecture-defining cell type proportions and postembryonic organ growth as key mechanisms forming the adult liver.

The default and directed pathways of hepatoblast differentiation involve distinct epigenomic mechanisms

The effectiveness of multiomics analyses in defining cell differentiation pathways during development is ambiguous. During liver development, hepatoblasts follow a default or directed pathway to differentiate into hepatocytes or cholangiocytes, respectively, and this provides a practical model to address this issue. Our study discovered that promoter-associated histone modifications and chromatin accessibility dynamics, rather than enhancer-associated histone modifications, effectively delineated the "default vs. directed" process of hepatoblast differentiation. Histone H3K27me3 on bivalent promoters is associated with this asymmetric differentiation strategy in mice and humans. We demonstrated that Ezh2 and Jmjd3 exert opposing regulatory roles in hepatoblast-cholangiocyte differentiation. Additionally, active enhancers, regulated by P300, correlate with the development of both hepatocytes and cholangiocytes. This research proposes a model highlighting the division of labor between promoters and enhancers, with promoter-associated chromatin modifications governing the "default vs. directed" differentiation mode of hepatoblasts, whereas enhancer-associated modifications primarily dictate the progressive development processes of hepatobiliary lineages.

SGAE: single-cell gene association entropy for revealing critical states of cell transitions during embryonic development

The critical point or pivotal threshold of cell transition occurs in early embryonic development when cell differentiation culminates in its transition to specific cell fates, at which the cell population undergoes an abrupt and qualitative shift. Revealing such critical points of cell transitions can track cellular heterogeneity and shed light on the molecular mechanisms of cell differentiation. However, precise detection of critical state transitions proves challenging when relying on single-cell RNA sequencing data due to their inherent sparsity, noise, and heterogeneity. In this study, diverging from conventional methods like differential gene analysis or static techniques that emphasize classification of cell types, an innovative computational approach, single-cell gene association entropy (SGAE), is designed for the analysis of single-cell RNA-seq data and utilizes gene association information to reveal critical states of cell transitions. More specifically, through the translation of gene expression data into local SGAE scores, the proposed SGAE can serve as an index to quantitatively assess the resilience and critical properties of genetic regulatory networks, consequently detecting the signal of cell transitions. Analyses of five single-cell datasets for embryonic development demonstrate that the SGAE method achieves better performance in facilitating the characterization of a critical phase transition compared with other existing methods. Moreover, the SGAE value can effectively discriminate cellular heterogeneity over time and performs well in the temporal clustering of cells. Besides, biological functional analysis also indicates the effectiveness of the proposed approach.

Fetal liver development and implications for liver disease pathogenesis

The metabolic, digestive and homeostatic roles of the liver are dependent on proper crosstalk and organization of hepatic cell lineages. These hepatic cell lineages are derived from their respective progenitors early in organogenesis in a spatiotemporally controlled manner, contributing to the liver's specialized and diverse microarchitecture. Advances in genomics, lineage tracing and microscopy have led to seminal discoveries in the past decade that have elucidated liver cell lineage hierarchies. In particular, single-cell genomics has enabled researchers to explore diversity within the liver, especially early in development when the application of bulk genomics was previously constrained due to the organ's small scale, resulting in low cell numbers. These discoveries have substantially advanced our understanding of cell differentiation trajectories, cell fate decisions, cell lineage plasticity and the signalling microenvironment underlying the formation of the liver. In addition, they have provided insights into the pathogenesis of liver disease and cancer, in which developmental processes participate in disease emergence and regeneration. Future work will focus on the translation of this knowledge to optimize in vitro models of liver development and fine-tune regenerative medicine strategies to treat liver disease. In this Review, we discuss the emergence of hepatic parenchymal and non-parenchymal cells, advances that have been made in in vitro modelling of liver development and draw parallels between developmental and pathological processes.

Single-Cell Heterogeneity in Snake Venom Expression Is Hardwired by Co-Option of Regulators from Progressively Activated Pathways

The ubiquitous cellular heterogeneity underlying many organism-level phenotypes raises questions about what factors drive this heterogeneity and how these complex heterogeneous systems evolve. Here, we use single-cell expression data from a Prairie rattlesnake (Crotalus viridis) venom gland to evaluate hypotheses for signaling networks underlying snake venom regulation and the degree to which different venom gene families have evolutionarily recruited distinct regulatory architectures. Our findings suggest that snake venom regulatory systems have evolutionarily co-opted trans-regulatory factors from extracellular signal-regulated kinase and unfolded protein response pathways that specifically coordinate expression of distinct venom toxins in a phased sequence across a single population of secretory cells. This pattern of co-option results in extensive cell-to-cell variation in venom gene expression, even between tandemly duplicated paralogs, suggesting this regulatory architecture has evolved to circumvent cellular constraints. While the exact nature of such constraints remains an open question, we propose that such regulatory heterogeneity may circumvent steric constraints on chromatin, cellular physiological constraints (e.g., endoplasmic reticulum stress or negative protein-protein interactions), or a combination of these. Regardless of the precise nature of these constraints, this example suggests that, in some cases, dynamic cellular constraints may impose previously unappreciated secondary constraints on the evolution of gene regulatory networks that favors heterogeneous expression.

A Refined Single Cell Landscape of Haematopoiesis in the Mouse Foetal Liver

During prenatal life, the foetal liver is colonised by several waves of haematopoietic progenitors to act as the main haematopoietic organ. Single cell (sc) RNA-seq has been used to identify foetal liver cell types via their transcriptomic signature and to compare gene expression patterns as haematopoietic development proceeds. To obtain a refined single cell landscape of haematopoiesis in the foetal liver, we have generated a scRNA-seq dataset from a whole mouse E12.5 liver that includes a larger number of cells than prior datasets at this stage and was obtained without cell type preselection to include all liver cell populations. We combined mining of this dataset with that of previously published datasets at other developmental stages to follow transcriptional dynamics as well as the cell cycle state of developing haematopoietic lineages. Our findings corroborate several prior reports on the timing of liver colonisation by haematopoietic progenitors and the emergence of differentiated lineages and provide further molecular characterisation of each cell population. Extending these findings, we demonstrate the existence of a foetal intermediate haemoglobin profile in the mouse, similar to that previously identified in humans, and a previously unidentified population of primitive erythroid cells in the foetal liver.

Multiple Facets of Cellular Homeostasis and Regeneration of the Mammalian Liver

Liver regeneration occurs in response to diverse injuries and is capable of functionally reestablishing the lost parenchyma. This phenomenon has been known since antiquity, encapsulated in the Greek myth where Prometheus was to be punished by Zeus for sharing the gift of fire with humanity by having an eagle eat his liver daily, only to have the liver regrow back, thus ensuring eternal suffering and punishment. Today, this process is actively leveraged clinically during living donor liver transplantation whereby up to a two-thirds hepatectomy (resection or removal of part of the liver) on a donor is used for transplant to a recipient. The donor liver rapidly regenerates to recover the lost parenchymal mass to form a functional tissue. This astonishing regenerative process and unique capacity of the liver are examined in further detail in this review.

Rat liver extracellular matrix and perfusion bioreactor culture promote human amnion epithelial cell differentiation towards hepatocyte-like cells

Congenital and chronic liver diseases have a substantial health burden worldwide. The most effective treatment available for these patients is whole organ transplantation; however, due to the severely limited supply of donor livers and the side effects associated with the immunosuppressive regimen required to accept allograft, the mortality rate in patients with end-stage liver disease is annually rising. Stem cell-based therapy aims to provide alternative treatments by either cell transplantation or bioengineered construct transplantation. Human amnion epithelial cells (AEC) are a widely available, ethically neutral source of cells with the plasticity and potential of multipotent stem cells and immunomodulatory properties of perinatal cells. AEC have been proven to be able to achieve functional improvement towards hepatocyte-like cells, capable of rescuing animals with metabolic disorders; however, they showed limited metabolic activities in vitro. Decellularised extracellular matrix (ECM) scaffolds have gained recognition as adjunct biological support. Decellularised scaffolds maintain native ECM components and the 3D architecture instrumental of the organ, necessary to support cells' maturation and function. We combined ECM-scaffold technology with primary human AEC, which we demonstrated being equipped with essential ECM-adhesion proteins, and evaluated the effects on AEC differentiation into functional hepatocyte-like cells (HLC). This novel approach included the use of a custom 4D bioreactor to provide constant oxygenation and media perfusion to cells in 3D cultures over time. We successfully generated HLC positive for hepatic markers such as ALB, CYP3A4 and CK18. AEC-derived HLC displayed early signs of hepatocyte phenotype, secreted albumin and urea, and expressed Phase-1 and -2 enzymes. The combination of liver-specific ECM and bioreactor provides a system able to aid differentiation into HLC, indicating that the innovative perfusion ECM-scaffold technology may support the functional improvement of multipotent and pluripotent stem cells, with important repercussions in the bioengineering of constructs for transplantation.

Distance covariance entropy reveals primed states and bifurcation dynamics in single-cell RNA-Seq data

Cell-fate transitions are fundamental to development and differentiation. Studying them with single-cell omic data is important to advance our understanding of the cell-fate commitment process, yet this remains challenging. Here we present a computational method called DICE, which analyzes the entropy of expression covariation patterns and which is applicable to static and dynamically changing cell populations. Using only single-cell RNA-Seq data, DICE is able to predict multipotent primed states and their regulatory factors, which we subsequently validate with single-cell epigenomic data. DICE reveals that primed states are often defined by epigenetic regulators or pioneer factors alongside lineage-specific transcription factors. In developmental time course single-cell RNA-Seq datasets, DICE can pinpoint the timing of bifurcations more precisely than lineage-trajectory inference algorithms or competing variance-based methods. In summary, by studying the dynamic changes of expression covariation entropy, DICE can help elucidate primed states and bifurcation dynamics without the need for single-cell epigenomic data.

Hippo signaling instructs ectopic but not normal organ growth

The Hippo signaling pathway is widely considered a master regulator of organ growth because of the prominent overgrowth phenotypes caused by experimental manipulation of its activity. Contrary to this model, we show here that removing Hippo transcriptional output did not impair the ability of the mouse liver and Drosophila eyes to grow to their normal size. Moreover, the transcriptional activity of the Hippo pathway effectors Yap/Taz/Yki did not correlate with cell proliferation, and hyperactivation of these effectors induced gene expression programs that did not recapitulate normal development. Concordantly, a functional screen in Drosophila identified several Hippo pathway target genes that were required for ectopic overgrowth but not normal growth. Thus, Hippo signaling does not instruct normal growth, and the Hippo-induced overgrowth phenotypes are caused by the activation of abnormal genetic programs.

Identification of fetal liver stroma in spectral cytometry using the parameter autofluorescence

The fetal liver (FL) is the main hematopoietic organ during embryonic development. The FL is also the unique anatomical site where hematopoietic stem cells expand before colonizing the bone marrow, where they ensure life-long blood cell production and become mostly resting. The identification of the different cell types that comprise the hematopoietic stroma in the FL is essential to understand the signals required for the expansion and differentiation of the hematopoietic stem cells. We used a panel of monoclonal antibodies to identify FL stromal cells in a 5-laser equipped spectral flow cytometry (FCM) analyzer. The "Autofluorescence Finder" of SONY ID7000 software identified two distinct autofluorescence emission spectra. Using autofluorescence as a fluorescence parameter we could assign the two autofluorescent signals to three distinct cell types and identified surface markers that characterize these populations. We found that one autofluorescent population corresponds to hepatoblast-like cells and cholangiocytes whereas the other expresses mesenchymal transcripts and was identified as stellate cells. Importantly, after birth, autofluorescence becomes the unique identifying property of hepatoblast-like cells because mature cholangiocytes are no longer autofluorescent. These results show that autofluorescence used as a parameter in spectral FCM is a useful tool to identify new cell subsets that are difficult to analyze in conventional FCM.

Comparative Transcriptome Analysis Provides Novel Molecular Events for the Differentiation and Maturation of Hepatocytes during the Liver Development of Zebrafish

The liver plays an essential role in multiple biological functions including metabolism, detoxification, digestion, coagulation, and homeostasis in vertebrates. The specification and differentiation of embryonic hepatoblasts, the proliferation of hepatocytes, and the hepatic tissue architecture are well documented, but molecular events governing the maturation of hepatocytes during liver development remain largely unclear. In this study, we performed a comparative transcriptome analysis of hepatocytes that were sorted by flow cytometry from developing zebrafish embryos at 60, 72, and 96 hpf. We identified 667 up-regulated and 3640 down-regulated genes in hepatocytes between 60 and 72 hpf, 606 up-regulated and 3924 down-regulated genes between 60 and 96 hpf, and 1693 up-regulated genes and 1508 down-regulated genes between 72 and 96 hpf. GO enrichment analysis revealed that key biological processes, cellular components, and molecular functions in hepatocytes between 60 to 72 hpf, such as cell cycle, DNA replication, DNA repair, RNA processing, and transcription regulation, are mainly associated with the proliferation of hepatocytes. In addition to biological processes, cellular components, and molecular functions for cell proliferation, molecular functions for carbohydrate metabolism were enriched in hepatocytes during 72 to 96 hpf. KEGG enrichment analysis identified key signaling pathways, such as cell cycle, RNA degradation, ubiquitin-mediated proteolysis, ErbB and Hedgehog signaling, basal transcription factors, Wnt signaling, and glycan degradation, which are closely associated with cell proliferation or carbohydrate metabolism in hepatocytes between 60 to 72 hpf. Newly enriched signaling pathways in hepatocytes during 72 to 96 hpf include metabolisms of pyrimidine, purine, nicotinate and nicotinamide, caffeine, glycine, serine and threonine, ABC transporters, and p53 signaling that function in metabolisms of lipid, protein and energy, cellular secretion, or detoxification, indicating the functional maturation of hepatocytes between 72 to 96 hpf. These findings provide novel clues for further understanding the functional differentiation and maturation of hepatocytes during liver development.

SCDD: a novel single-cell RNA-seq imputation method with diffusion and denoising

Single-cell sequencing technologies are widely used to discover the evolutionary relationships and the differences in cells. Since dropout events may frustrate the analysis, many imputation approaches for single-cell RNA-seq data have appeared in previous attempts. However, previous imputation attempts usually suffer from the over-smooth problem, which may bring limited improvement or negative effect for the downstream analysis of single-cell RNA-seq data. To solve this difficulty, we propose a novel two-stage diffusion-denoising method called SCDD for large-scale single-cell RNA-seq imputation in this paper. We introduce the diffusion i.e. a direct imputation strategy using the expression of similar cells for potential dropout sites, to perform the initial imputation at first. After the diffusion, a joint model integrated with graph convolutional neural network and contractive autoencoder is developed to generate superposition states of similar cells, from which we restore the original states and remove the noise introduced by the diffusion. The final experimental results indicate that SCDD could effectively suppress the over-smooth problem and remarkably improve the effect of single-cell RNA-seq downstream analysis, including clustering and trajectory analysis.

Cross-Organ Transcriptomic Comparison Reveals Universal Factors During Maturation

Various cell types can be derived from stem cells. However, these cells are immature and do not match their adult counterparts in functional capabilities, limiting their use in disease modeling and cell therapies. Thus, it is crucial to understand the mechanisms of maturation in vivo. However, it is unknown if there are genes and pathways conserved across organs during maturation. To address this, we performed a time-series analysis of the transcriptome of the mouse heart, brain, liver, and kidney and analyzed their trajectories over time. In addition, gene regulatory networks were reconstructed to determine overlapping expression patterns. Based on these, we identified commonly upregulated and downregulated pathways across all four organs. Key upstream regulators were also predicted based on the temporal expression of downstream genes. These findings suggest the presence of universal regulators during organ maturation, which may help us develop a general strategy to mature stem cell-derived cells in vitro.

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

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

•

Identified minimal regulatory network to model liver development and regeneration

•

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

•

Ability to explain physiological spatial patterning of liver cell types

•

Decoded strategies for efficient reprogramming among liver cell phenotypes

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

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

Multiscale Methods for Signal Selection in Single-Cell Data

Analysis of single-cell transcriptomics often relies on clustering cells and then performing differential gene expression (DGE) to identify genes that vary between these clusters. These discrete analyses successfully determine cell types and markers; however, continuous variation within and between cell types may not be detected. We propose three topologically motivated mathematical methods for unsupervised feature selection that consider discrete and continuous transcriptional patterns on an equal footing across multiple scales simultaneously. Eigenscores (eigi) rank signals or genes based on their correspondence to low-frequency intrinsic patterning in the data using the spectral decomposition of the Laplacian graph. The multiscale Laplacian score (MLS) is an unsupervised method for locating relevant scales in data and selecting the genes that are coherently expressed at these respective scales. The persistent Rayleigh quotient (PRQ) takes data equipped with a filtration, allowing the separation of genes with different roles in a bifurcation process (e.g., pseudo-time). We demonstrate the utility of these techniques by applying them to published single-cell transcriptomics data sets. The methods validate previously identified genes and detect additional biologically meaningful genes with coherent expression patterns. By studying the interaction between gene signals and the geometry of the underlying space, the three methods give multidimensional rankings of the genes and visualisation of relationships between them.

Human Cytomegalovirus-IE2 Affects Embryonic Liver Development and Survival in Transgenic Mouse

BACKGROUND & AIMS

Congenital human cytomegalovirus (HCMV) infection is a common cause of liver injury. The major immediate-early protein 2 (IE2) of HCMV is critical for the progression of HCMV infection. As a result of species isolation, there are no animal models suitable for HCMV infection, which aimed to study the long-term effects of IE2 on embryonic liver development in vivo. Hence, this study aimed to investigate the role of IE2 in liver development using a transgenosis mouse model.

METHODS

Rosa26-Loxp-STOP-Loxp (LAS)-IE2^+/-, cre mice that could specifically and stably express IE2 in the liver, were constructed. Phenotypic analysis, immunolocalization studies, messenger RNA analyses, transcriptome sequencing, and flow cytometry analysis were performed on Rosa26-LSL-IE2^+/-, cre mice during hepatogenesis.

RESULTS

Rosa26-LSL-IE2^+/-, cre mice could consistently express IE2 at different embryonic stages in vivo. With the development of Rosa26-LSL-IE2^+/-, cre embryos from embryonic day 17.5 to postnatal day 1, progressive liver hypoplasia and embryonic deaths were observed. Furthermore, molecular evidence that IE2 expression inhibited hepatocyte proliferation, increased cell apoptosis, and impaired hepatocyte maturation was provided.

CONCLUSIONS

Rosa26-LSL-IE2^+/-, cre mice could stably express IE2 in the liver. IE2 expression resulted in embryonic liver hypoplasia by disrupting hepatic morphogenesis and hepatocyte maturation, which may be responsible for embryonic deaths. This study is helpful in understanding the mechanism of liver injuries induced by HCMV infection.

Hippo signalling in the liver: role in development, regeneration and disease

The Hippo signalling pathway has emerged as a major player in many aspects of liver biology, such as development, cell fate determination, homeostatic function and regeneration from injury. The regulation of Hippo signalling is complex, with activation of the pathway by diverse upstream inputs including signals from cellular adhesion, mechanotransduction and crosstalk with other signalling pathways. Pathological activation of the downstream transcriptional co-activators yes-associated protein 1 (YAP) and transcriptional co-activator with PDZ-binding motif (TAZ, encoded by WWTR1), which are negatively regulated by Hippo signalling, has been implicated in multiple aspects of chronic liver disease, such as the development of liver fibrosis and tumorigenesis. Thus, development of pharmacological inhibitors of YAP-TAZ signalling has been an area of great interest. In this Review, we summarize the diverse roles of Hippo signalling in liver biology and highlight areas where outstanding questions remain to be investigated. Greater understanding of the mechanisms of Hippo signalling in liver function should help facilitate the development of novel therapies for the treatment of liver disease.

Prognostic Pathways Guide Drug Indications in Pan-Cancers

Pathway-level analysis is a powerful approach enabling the interpretation of post-genomic data at a higher level than that of individual molecules. Molecular-targeted therapy focusing on cascade signaling pathways has become a new paradigm in anticancer therapy, instead of a single protein. However, the approaches to narrowing down the long list of biological pathways are limited. Here, we proposed a strategy for in silico Drug Prescription on biological pathways across pan-Cancers (CDP), by connecting drugs to candidate pathways. Applying on a list of 120 traditional Chinese medicines (TCM), we especially identified the “TCM–pathways–cancers” triplet and constructed it into a heterogeneous network across pan-cancers. Applying them into TCMs, the computational prescribing methods deepened the understanding of the efficacy of TCM at the molecular level. Further applying them into Western medicines, CDP could promote drug reposition avoiding time-consuming developments of new drugs.

Epigenetics of alcohol-related liver diseases

Alcohol-related liver disease (ARLD) is a primary cause of chronic liver disease in the United States. Despite advances in the diagnosis and management of ARLD, it remains a major public health problem associated with significant morbidity and mortality, emphasising the need to adopt novel approaches to the study of ARLD and its complications. Epigenetic changes are increasingly being recognised as contributing to the pathogenesis of multiple disease states. Harnessing the power of innovative technologies for the study of epigenetics (e.g., next-generation sequencing, DNA methylation assays, histone modification profiling and computational techniques like machine learning) has resulted in a seismic shift in our understanding of the pathophysiology of ARLD. Knowledge of these techniques and advances is of paramount importance for the practicing hepatologist and researchers alike. Accordingly, in this review article we will summarise the current knowledge about alcohol-induced epigenetic alterations in the context of ARLD, including but not limited to, DNA hyper/hypo methylation, histone modifications, changes in non-coding RNA, 3D chromatin architecture and enhancer-promoter interactions. Additionally, we will discuss the state-of-the-art techniques used in the study of ARLD (e.g. single-cell sequencing). We will also highlight the epigenetic regulation of chemokines and their proinflammatory role in the context of ARLD. Lastly, we will examine the clinical applications of epigenetics in the diagnosis and management of ARLD.

TraSig: inferring cell-cell interactions from pseudotime ordering of scRNA-Seq data

A major advantage of single cell RNA-sequencing (scRNA-Seq) data is the ability to reconstruct continuous ordering and trajectories for cells. Here we present TraSig, a computational method for improving the inference of cell-cell interactions in scRNA-Seq studies that utilizes the dynamic information to identify significant ligand-receptor pairs with similar trajectories, which in turn are used to score interacting cell clusters. We applied TraSig to several scRNA-Seq datasets and obtained unique predictions that improve upon those identified by prior methods. Functional experiments validate the ability of TraSig to identify novel signaling interactions that impact vascular development in liver organoids.Software https://github.com/doraadong/TraSig .

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

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

Quantitative modeling identifies critical cell mechanics driving bile duct lumen formation

Biliary ducts collect bile from liver lobules, the smallest functional and anatomical units of liver, and carry it to the gallbladder. Disruptions in this process caused by defective embryonic development, or through ductal reaction in liver disease have a major impact on life quality and survival of patients. A deep understanding of the processes underlying bile duct lumen formation is crucial to identify intervention points to avoid or treat the appearance of defective bile ducts. Several hypotheses have been proposed to characterize the biophysical mechanisms driving initial bile duct lumen formation during embryogenesis. Here, guided by the quantification of morphological features and expression of genes in bile ducts from embryonic mouse liver, we sharpened these hypotheses and collected data to develop a high resolution individual cell-based computational model that enables to test alternative hypotheses in silico. This model permits realistic simulations of tissue and cell mechanics at sub-cellular scale. Our simulations suggest that successful bile duct lumen formation requires a simultaneous contribution of directed cell division of cholangiocytes, local osmotic effects generated by salt excretion in the lumen, and temporally-controlled differentiation of hepatoblasts to cholangiocytes, with apical constriction of cholangiocytes only moderately affecting luminal size.

The initial step in bile duct development is the formation of a biliary lumen, a process which involves several cellular mechanisms, such as cell division and polarization, and secretion of fluid. However, how these mechanisms are orchestrated in time and space is difficult to understand. Here, we built a computational model of biliary lumen formation which represents every cell and its function in detail. With the model we can simulate the effect of biophysical aspects that affect duct formation. We have tested the individual and combined effects of directed cell division, apical constriction, and osmotic effects on lumen expansion by varying the parameters that control their relative strength. Our simulations suggest that successful bile duct lumen formation requires the simultaneous contribution of directed cell division of cholangiocytes, local osmotic effects generated by salt excretion in the lumen, and temporally-controlled differentiation of hepatoblasts to cholangiocytes, with apical constriction of cholangiocytes only moderately affecting luminal size.

Role of YAP1 Signaling in Biliary Development, Repair, and Disease

Yes-associated protein 1 (YAP1) is a transcriptional coactivator that activates transcriptional enhanced associate domain transcription factors upon inactivation of the Hippo signaling pathway, to regulate biological processes like proliferation, survival, and differentiation. YAP1 is most prominently expressed in biliary epithelial cells (BECs) in normal adult livers and during development. In the current review, we will discuss the multiple roles of YAP1 in the development and morphogenesis of bile ducts inside and outside the liver, as well as in orchestrating the cholangiocyte repair response to biliary injury. We will review how biliary repair can occur through the process of hepatocyte-to-BEC transdifferentiation and how YAP1 is pertinent to this process. We will also discuss the liver's capacity for metabolic reprogramming as an adaptive mechanism in extreme cholestasis, such as when intrahepatic bile ducts are absent due to YAP1 loss from hepatic progenitors. Finally, we will discuss the roles of YAP1 in the context of pediatric pathologies afflicting bile ducts, such as Alagille syndrome and biliary atresia. In conclusion, we will comprehensively discuss the spatiotemporal roles of YAP1 in biliary development and repair after biliary injury while describing key interactions with other well-known developmental pathways.

RNF43/ZNRF3 loss predisposes to hepatocellular-carcinoma by impairing liver regeneration and altering the liver lipid metabolic ground-state

RNF43/ZNRF3 negatively regulate WNT signalling. Both genes are mutated in several types of cancers, however, their contribution to liver disease is unknown. Here we describe that hepatocyte-specific loss of Rnf43/Znrf3 results in steatohepatitis and in increase in unsaturated lipids, in the absence of dietary fat supplementation. Upon injury, Rnf43/Znrf3 deletion results in defective hepatocyte regeneration and liver cancer, caused by an imbalance between differentiation/proliferation. Using hepatocyte-, hepatoblast- and ductal cell-derived organoids we demonstrate that the differentiation defects and lipid alterations are, in part, cell-autonomous. Interestingly, ZNRF3 mutant liver cancer patients present poorer prognosis, altered hepatic lipid metabolism and steatohepatitis/NASH signatures. Our results imply that RNF43/ZNRF3 predispose to liver cancer by controlling the proliferative/differentiation and lipid metabolic state of hepatocytes. Both mechanisms combined facilitate the progression towards malignancy. Our findings might aid on the management of those RNF43/ZNRF3 mutated individuals at risk of developing fatty liver and/or liver cancer.

Harnessing orthogonal recombinases to decipher cell fate with enhanced precision

Precisely deciphering the cellular plasticity in vivo is essential in understanding many key biological processes. Site-specific recombinases are genetic tools used for in vivo lineage tracing and gene manipulation. Conventional Cre-loxP, Dre-rox, and Flp-frt technologies form the orthogonal recombination systems that can also be used in combination to increase the precision. As such, more than one marker gene can be targeted for lineage tracing, studying cellular heterogeneity, recording cellular activities, or even genome editing. Their combinatory use has recently resolved some controversies in defining cellular fate plasticity. Focusing on cell fate studies, we introduce the design principles of orthogonal recombinases-based strategies, describe some working examples in resolving cell fate-related controversies, and discuss some of their technical strengths and limits.

Genomic characterization of rare molecular subclasses of hepatocellular carcinoma

Primary liver cancer, consisting of both cholangiocarcinoma (CCA) and hepatocellular carcinoma (HCC), is the second leading cause of cancer deaths worldwide. Our goal is to genomically characterize rare HCC subclasses to provide insight into disease biology. Leveraging The Cancer Genome Atlas (TCGA) to perform a combined analysis of CCA (n = 36) and HCC (n = 275), we integrated multiple genomic platforms, to assess transcriptional profiles, mutational signatures, and copy number patterns to uncover underlying etiology and linage specific patterns. We identified two molecular classes distinct from prototypical HCC tumors. The first, CCA-Like, although histologically indistinguishable from HCC, had enrichment of CCA mutations (IDH1, BAP1), mutational signatures, and transcriptional patterns (SOX9, KRT19). CCA-Like, however, retained a copy number landscape similar to HCC, suggesting a hepatocellular linage. The second, Blast-Like, is enriched in TP53 mutations, HBV infection, exposure related mutational signatures and transcriptionally similar to hepatoblasts. Although these subclasses are molecularly distinct, they both have a worse progression-free survival compared to classical HCC tumors, yet are clinically treated the same. The identification of and characterization of CCA-Like and Blast-Like subclasses advance our knowledge of HCC as well as represents an urgent need for the identification of class specific biomarkers and targeted therapy.

Jeffrey Damrauer, Markia Smith et al. used existing datasets from cholangiocarcinoma (CCA) and hepatocellular carcinoma (HCC) to characterize two subsets of HCC distinct from prototypical HCC tumors, based on comprehensive analysis of molecular data. The two classes differed from HCC by their copy number, gene expression and mutational signature and exhibited worse progression free survival, highlighting the need to identify class-specific biomarkers and develop targeted therapies for these forms of cancer.