Single-cell RNA-Seq analysis reveals dynamic trajectories during mouse liver development

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.

Unveiling inflammatory and prehypertrophic cell populations as key contributors to knee cartilage degeneration in osteoarthritis using multi-omics data integration

OBJECTIVES

Single-cell and spatial transcriptomics analysis of human knee articular cartilage tissue to present a comprehensive transcriptome landscape and osteoarthritis (OA)-critical cell populations.

METHODS

Single-cell RNA sequencing and spatially resolved transcriptomic technology have been applied to characterise the cellular heterogeneity of human knee articular cartilage which were collected from 8 OA donors, and 3 non-OA control donors, and a total of 19 samples. The novel chondrocyte population and marker genes of interest were validated by immunohistochemistry staining, quantitative real-time PCR, etc. The OA-critical cell populations were validated through integrative analyses of publicly available bulk RNA sequencing data and large-scale genome-wide association studies.

RESULTS

We identified 33 cell population-specific marker genes that define 11 chondrocyte populations, including 9 known populations and 2 new populations, that is, pre-inflammatory chondrocyte population (preInfC) and inflammatory chondrocyte population (InfC). The novel findings that make this an important addition to the literature include: (1) the novel InfC activates the mediator MIF-CD74; (2) the prehypertrophic chondrocyte (preHTC) and hypertrophic chondrocyte (HTC) are potentially OA-critical cell populations; (3) most OA-associated differentially expressed genes reside in the articular surface and superficial zone; (4) the prefibrocartilage chondrocyte (preFC) population is a major contributor to the stratification of patients with OA, resulting in both an inflammatory-related subtype and a non-inflammatory-related subtype.

CONCLUSIONS

Our results highlight InfC, preHTC, preFC and HTC as potential cell populations to target for therapy. Also, we conclude that profiling of those cell populations in patients might be used to stratify patient populations for defining cohorts for clinical trials and precision medicine.

Characterization of dUTPase expression in mouse postnatal development and adult neurogenesis

The enzyme dUTPase has an essential role in maintaining genomic integrity. In mouse, nuclear and mitochondrial isoforms of the enzyme have been described. Here we present the isoform-specific mRNA expression levels in different murine organs during development using RT-qPCR. In this study, we analyzed organs of 14.5-day embryos and of postnatal 2-, 4-, 10-week- and 13-month-old mice. We demonstrate organ-, sex- and developmental stage-specific differences in the mRNA expression levels of both isoforms. We found high mRNA expression level of the nuclear isoform in the embryo brain, and the expression level remained relatively high in the adult brain as well. This was surprising, since dUTPase is known to play an important role in proliferating cells, and mass production of neural cells is completed by adulthood. Thus, we investigated the pattern of the dUTPase protein expression specifically in the adult brain with immunostaining and found that dUTPase is present in the germinative zones, the subventricular and the subgranular zones, where neurogenesis occurs and in the rostral migratory stream where neuroblasts migrate to the olfactory bulb. These novel findings suggest that dUTPase may have a role in cell differentiation and indicate that accurate dTTP biosynthesis can be vital, especially in neurogenesis.

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.

Low-dose biliatresone treatment of pregnant mice causes subclinical biliary disease in their offspring: Evidence for a spectrum of neonatal injury

Biliary atresia is a neonatal disease characterized by damage, inflammation, and fibrosis of the liver and bile ducts and by abnormal bile metabolism. It likely results from a prenatal environmental exposure that spares the mother and affects the fetus. Our aim was to develop a model of fetal injury by exposing pregnant mice to low-dose biliatresone, a plant toxin implicated in biliary atresia in livestock, and then to determine whether there was a hepatobiliary phenotype in their pups. Pregnant mice were treated orally with 15 mg/kg/d biliatresone for 2 days. Histology of the liver and bile ducts, serum bile acids, and liver immune cells of pups from treated mothers were analyzed at P5 and P21. Pups had no evidence of histological liver or bile duct injury or fibrosis at either timepoint. In addition, growth was normal. However, serum levels of glycocholic acid were elevated at P5, suggesting altered bile metabolism, and the serum bile acid profile became increasingly abnormal through P21, with enhanced glycine conjugation of bile acids. There was also immune cell activation observed in the liver at P21. These results suggest that prenatal exposure to low doses of an environmental toxin can cause subclinical disease including liver inflammation and aberrant bile metabolism even in the absence of histological changes. This finding suggests a wide potential spectrum of disease after fetal biliary injury.

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.

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.

GITRL impairs hepatocyte repopulation by liver progenitor cells to aggravate inflammation and fibrosis by GITR+CD8+ T lymphocytes in CDE Mice

As an alternative pathway for liver regeneration, liver progenitor cells and their derived ductular reaction cells increase during the progression of many chronic liver diseases. However, the mechanism underlying their hepatocyte repopulation after liver injury remains unknown. Here, we conducted progenitor cell lineage tracing in mice and found that fewer than 2% of hepatocytes were derived from liver progenitor cells after 9 weeks of injury with a choline-deficient diet supplemented with ethionine (CDE), and this percentage increased approximately three-fold after 3 weeks of recovery. We also found that the proportion of liver progenitor cells double positive for the ligand of glucocorticoid-induced tumour necrosis factor receptor (GITRL, also called Tnfsf18) and SRY-related HMG box transcription 9 (Sox9) among nonparenchymal cells increased time-dependently upon CDE injury and reduced after recovery. When GITRL was conditionally knocked out from hepatic progenitor cells, its expression in nonparenchymal cells was downregulated by approximately fifty percent, and hepatocyte repopulation increased by approximately three folds. Simultaneously, conditional knockout of GITRL reduced the proportion of liver-infiltrating CD8+ T lymphocytes and glucocorticoid-induced tumour necrosis factor receptor (GITR)-positive CD8+ T lymphocytes. Mechanistically, GITRL stimulated cell proliferation but suppressed the differentiation of liver progenitor organoids into hepatocytes, and CD8+ T cells further reduced their hepatocyte differentiation by downregulating the Wnt/β-catenin pathway. Therefore, GITRL expressed by liver progenitor cells impairs hepatocyte differentiation, thus hindering progenitor cell-mediated liver regeneration.

Sox7-positive endothelial progenitors establish coronary arteries and govern ventricular compaction

The cardiac endothelium influences ventricular chamber development by coordinating trabeculation and compaction. However, the endothelial-specific molecular mechanisms mediating this coordination are not fully understood. Here, we identify the Sox7 transcription factor as a critical cue instructing cardiac endothelium identity during ventricular chamber development. Endothelial-specific loss of Sox7 function in mice results in cardiac ventricular defects similar to non-compaction cardiomyopathy, with a change in the proportions of trabecular and compact cardiomyocytes in the mutant hearts. This phenotype is paralleled by abnormal coronary artery formation. Loss of Sox7 function disrupts the transcriptional regulation of the Notch pathway and connexins 37 and 40, which govern coronary arterial specification. Upon Sox7 endothelial-specific deletion, single-nuclei transcriptomics analysis identifies the depletion of a subset of Sox9/Gpc3-positive endocardial progenitor cells and an increase in erythro-myeloid cell lineages. Fate mapping analysis reveals that a subset of Sox7-null endothelial cells transdifferentiate into hematopoietic but not cardiomyocyte lineages. Our findings determine that Sox7 maintains cardiac endothelial cell identity, which is crucial to the cellular cross-talk that drives ventricular compaction and coronary artery development.

scLiverDB: a Database of Human and Mouse Liver Transcriptome Landscapes at Single-Cell Resolution

The liver is critical for the digestive and immune systems. Although the physiology and pathology of liver have been well studied and many scRNA-seq data are generated, a database and landscape for characterizing cell types and gene expression in different liver diseases or developmental stages at single-cell resolution are lacking. Hence, scLiverDB is developed, a specialized database for human and mouse liver transcriptomes to unravel the landscape of liver cell types, cell heterogeneity and gene expression at single-cell resolution across various liver diseases/cell types/developmental stages. To date, 62 datasets including 9,050 samples and 1,741,734 cells is curated. A uniform workflow is used, which included quality control, dimensional reduction, clustering, and cell-type annotation to analyze datasets on the same platform; integrated manual and automatic methods for accurate cell-type identification and provided a user-friendly web interface with multiscale functions. There are two case studies to show the usefulness of scLiverDB, which identified the LTB (lymphotoxin Beta) gene as a potential biomarker of lymphoid cells differentiation and showed the expression changes of Foxa3 (forkhead box A3) in liver chronic progressive diseases. This work provides a crucial resource to resolve molecular and cellular information in normal, diseased, and developing human and mouse livers.

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.

Genome-wide mapping of cancer dependency genes and genetic modifiers of chemotherapy in high-risk hepatoblastoma

A lack of relevant genetic models and cell lines hampers our understanding of hepatoblastoma pathogenesis and the development of new therapies for this neoplasm. Here, we report an improved MYC-driven hepatoblastoma-like murine model that recapitulates the pathological features of embryonal type of hepatoblastoma, with transcriptomics resembling the high-risk gene signatures of the human disease. Single-cell RNA-sequencing and spatial transcriptomics identify distinct subpopulations of hepatoblastoma cells. After deriving cell lines from the mouse model, we map cancer dependency genes using CRISPR-Cas9 screening and identify druggable targets shared with human hepatoblastoma (e.g., CDK7, CDK9, PRMT1, PRMT5). Our screen also reveals oncogenes and tumor suppressor genes in hepatoblastoma that engage multiple, druggable cancer signaling pathways. Chemotherapy is critical for human hepatoblastoma treatment. A genetic mapping of doxorubicin response by CRISPR-Cas9 screening identifies modifiers whose loss-of-function synergizes with (e.g., PRKDC) or antagonizes (e.g., apoptosis genes) the effect of chemotherapy. The combination of PRKDC inhibition and doxorubicin-based chemotherapy greatly enhances therapeutic efficacy. These studies provide a set of resources including disease models suitable for identifying and validating potential therapeutic targets in human high-risk hepatoblastoma.

Monitoring Cultured Rat Hepatocytes Using RNA-Seq In Vitro

Compared to other techniques, RNA sequencing (RNA-Seq) has the advantage of having details of the expression abundance of all transcripts in a single run. In this study, we used RNA-Seq to monitor the maturity and dynamic characteristics of in vitro hepatocyte cultures. Hepatocytes, including mature hepatocytes and small hepatocytes, were analyzed in vitro using RNA-Seq and quantitative polymerase chain reaction (qPCR). The results demonstrated that the gene expression profiles measured by RNA-Seq showed a similar trend to the expression profiles measured by qPCR, and can be used to infer the success of in vitro hepatocyte cultures. The results of the differential analysis, which compared mature hepatocytes against small hepatocytes, revealed 836 downregulated and 137 upregulated genes. In addition, the success of the hepatocyte cultures could be explained by the gene list screened from the adopted gene enrichment test. In summary, we demonstrated that RNA-Seq could become an effective method for monitoring the whole transcriptome of hepatocyte cultures and provide a more comprehensive list of factors related to the differentiation of small hepatocytes into mature hepatocytes. This monitoring system not only shows high potential in medical applications but may also be a novel method for the clinical diagnosis of liver-related diseases.

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.

Inflammatory and immune etiology of type 2 diabetes

Type 2 diabetes (T2D) represents a global threat affecting millions of patients worldwide. However, its causes remain incompletely dissected and we lack the tools to predict which individuals will develop T2D. Although there is a clear proven clinical association of T2D with metabolic disorders such as obesity and nonalcoholic fatty liver disease (NAFLD), the existence of a significant number of nondiabetic obese subjects suggests yet-uncovered features of such relationships. Here, we propose that a significant proportion of individuals may harbor an immune profile that renders them susceptible to developing T2D. We note the heterogeneity of circulating monocytes and tissue macrophages in organs that are key to metabolic disorders such as liver, white adipose tissue (WAT), and endocrine pancreas, as well as their contribution to T2D genesis.

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.

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.

acorde unravels functionally interpretable networks of isoform co-usage from single cell data

Alternative splicing (AS) is a highly-regulated post-transcriptional mechanism known to modulate isoform expression within genes and contribute to cell-type identity. However, the extent to which alternative isoforms establish co-expression networks that may be relevant in cellular function has not been explored yet. Here, we present acorde, a pipeline that successfully leverages bulk long reads and single-cell data to confidently detect alternative isoform co-expression relationships. To achieve this, we develop and validate percentile correlations, an innovative approach that overcomes data sparsity and yields accurate co-expression estimates from single-cell data. Next, acorde uses correlations to cluster co-expressed isoforms into a network, unraveling cell type-specific alternative isoform usage patterns. By selecting same-gene isoforms between these clusters, we subsequently detect and characterize genes with co-differential isoform usage (coDIU) across cell types. Finally, we predict functional elements from long read-defined isoforms and provide insight into biological processes, motifs, and domains potentially controlled by the coordination of post-transcriptional regulation. The code for acorde is available at https://github.com/ConesaLab/acorde .

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.

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.

Application of single-cell RNA sequencing technology in liver diseases: a narrative review

Objective

This review aimed to summarize the application of single-cell transcriptome sequencing technology in liver diseases.

Background

The increasing application of single-cell ribonucleic acid (RNA) sequencing (scRNA-seq) in life science and biomedical research has greatly improved our understanding of cellular heterogeneity in immunology, oncology, and developmental biology. scRNA-seq has proven to be a powerful tool for identifying and classifying cell subsets, characterizing rare or small cell subsets and tracking cell differentiation along the dynamic cell stages. Globally, liver disease has high rates of morbidity and mortality, and its exact pathological mechanism remains unclear, current treatment options are limited to clearance of the underlying cause or liver transplantation, which cannot overwhelm and cure liver diseases. scRNA-seq provides many novel insights for healthy and diseased livers.

Methods

In this review, we searched for related articles in the PubMed database and summarized the advances of scRNA-seq in revealing the molecular mechanisms of liver development, regeneration, and disease. We also discussed the challenges and future application potential of scRNA-seq, which is expected to enhance the ability to explore the field of liver research and accelerate the clinical application of liver precision medicine.

Conclusions

With the continuous improvement of scRNA-seq technology, scRNA-seq is expected to unlock new avenues for liver biology exploration, liver disease diagnosis, and personalized treatment, which will pave the way for breakthrough innovation in personalized medicine.

Heterogeneity of Hepatic Stellate Cells in Fibrogenesis of the Liver: Insights from Single-Cell Transcriptomic Analysis in Liver Injury

Background & Aims: Liver fibrosis is a pathological healing process resulting from hepatic stellate cell (HSC) activation and the generation of myofibroblasts from activated HSCs. The precise underlying mechanisms of liver fibrogenesis are still largely vague due to lack of understanding the functional heterogeneity of activated HSCs during liver injury. Approach and Results: In this study, to define the mechanism of HSC activation, we performed the transcriptomic analysis at single-cell resolution (scRNA-seq) on HSCs in mice treated with carbon tetrachloride (CCl₄). By employing LRAT-Cre:Rosa26^mT/mG mice, we were able to isolate an activated GFP-positive HSC lineage derived cell population by fluorescence-activated cell sorter (FACS). A total of 8 HSC subpopulations were identified based on an unsupervised analysis. Each HSC cluster displayed a unique transcriptomic profile, despite all clusters expressing common mouse HSC marker genes. We demonstrated that one of the HSC subpopulations expressed high levels of mitosis regulatory genes, velocity, and monocle analysis indicated that these HSCs are at transitioning and proliferating phases at the beginning of HSCs activation and will eventually give rise to several other HSC subtypes. We also demonstrated cell clusters representing HSC-derived mature myofibroblast populations that express myofibroblasts hallmark genes with unique contractile properties. Most importantly, we found a novel HSC cluster that is likely to be critical in liver regeneration, immune reaction, and vascular remodeling, in which the unique profiles of genes such as Rgs5, Angptl6, and Meg3 are highly expressed. Lastly, we demonstrated that the heterogeneity of HSCs in the injured mouse livers is closely similar to that of cirrhotic human livers. Conclusions: Collectively, our scRNA-seq data provided insight into the landscape of activated HSC populations and the dynamic transitional pathway from HSC to myofibroblasts in response to liver injury.

Emerging Single-Cell Technological Approaches to Investigate Chromatin Dynamics and Centromere Regulation in Human Health and Disease

Epigenetic regulators play a crucial role in establishing and maintaining gene expression states. To date, the main efforts to study cellular heterogeneity have focused on elucidating the variable nature of the chromatin landscape. Specific chromatin organisation is fundamental for normal organogenesis and developmental homeostasis and can be affected by different environmental factors. The latter can lead to detrimental alterations in gene transcription, as well as pathological conditions such as cancer. Epigenetic marks regulate the transcriptional output of cells. Centromeres are chromosome structures that are epigenetically regulated and are crucial for accurate segregation. The advent of single-cell epigenetic profiling has provided finer analytical resolution, exposing the intrinsic peculiarities of different cells within an apparently homogenous population. In this review, we discuss recent advances in methodologies applied to epigenetics, such as CUT&RUN and CUT&TAG. Then, we compare standard and emerging single-cell techniques and their relevance for investigating human diseases. Finally, we describe emerging methodologies that investigate centromeric chromatin specification and neocentromere formation.

Nuclear Organization during Hepatogenesis in Zebrafish Requires Uhrf1

Acquisition of cellular fate during development is initiated and maintained by well-coordinated patterns of gene expression that are dictated by the epigenetic landscape and genome organization in the nucleus. While the epigenetic marks that mediate developmental gene expression patterns during organogenesis have been well studied, less is known about how epigenetic marks influence nuclear organization during development. This study examines the relationship between nuclear structure, chromatin accessibility, DNA methylation, and gene expression during hepatic outgrowth in zebrafish larvae. We investigate the relationship between these features using mutants that lack DNA methylation. Hepatocyte nuclear morphology was established coincident with hepatocyte differentiation at 80 h post-fertilization (hpf), and nuclear shape and size continued to change until the conclusion of outgrowth and morphogenesis at 120 hpf. Integrating ATAC-Seq analysis with DNA methylation profiling of zebrafish livers at 120 hpf showed that closed and highly methylated chromatin occupies most transposable elements and that open chromatin correlated with gene expression. DNA hypomethylation, due to mutation of genes encoding ubiquitin-like, containing PHD and RING Finger Domains 1 (uhrf1) and DNA methyltransferase (dnmt1), did not block hepatocyte differentiation, but had dramatic effects on nuclear organization. Hepatocytes in uhrf1 mutants have large, deformed nuclei with multiple nucleoli, downregulation of nucleolar genes, and a complete lack of the nuclear lamina. Loss of lamin B2 staining was phenocopied by dnmt1 mutation. Together, these data show that hepatocyte nuclear morphogenesis coincides with organ morphogenesis and outgrowth, and that DNA methylation directs chromatin organization, and, in turn, hepatocyte nuclear shape and size during liver development.

Single-nucleus RNA-seq2 reveals functional crosstalk between liver zonation and ploidy

Single-cell RNA-seq reveals the role of pathogenic cell populations in development and progression of chronic diseases. In order to expand our knowledge on cellular heterogeneity, we have developed a single-nucleus RNA-seq2 method tailored for the comprehensive analysis of the nuclear transcriptome from frozen tissues, allowing the dissection of all cell types present in the liver, regardless of cell size or cellular fragility. We use this approach to characterize the transcriptional profile of individual hepatocytes with different levels of ploidy, and have discovered that ploidy states are associated with different metabolic potential, and gene expression in tetraploid mononucleated hepatocytes is conditioned by their position within the hepatic lobule. Our work reveals a remarkable crosstalk between gene dosage and spatial distribution of hepatocytes.

Compensatory hepatic adaptation accompanies permanent absence of intrahepatic biliary network due to YAP1 loss in liver progenitors

Yes-associated protein 1 (YAP1) regulates cell plasticity during liver injury, regeneration, and cancer, but its role in liver development is unknown. We detect YAP1 activity in biliary cells and in cells at the hepatobiliary bifurcation in single-cell RNA sequencing analysis of developing livers. Deletion of Yap1 in hepatoblasts does not impair Notch-driven SOX9+ ductal plate formation but does prevent the formation of the abutting second layer of SOX9+ ductal cells, blocking the formation of a patent intrahepatic biliary tree. Intriguingly, these mice survive for 8 months with severe cholestatic injury and without hepatocyte-to-biliary transdifferentiation. Ductular reaction in the perihilar region suggests extrahepatic biliary proliferation, likely seeking the missing intrahepatic biliary network. Long-term survival of these mice occurs through hepatocyte adaptation via reduced metabolic and synthetic function, including altered bile acid metabolism and transport. Overall, we show YAP1 as a key regulator of bile duct development while highlighting a profound adaptive capability of hepatocytes.

TNFAIP8 Regulates Intestinal Epithelial Cell Differentiation and May Alter Terminal Differentiation of Secretory Progenitors

The intestine is a highly proliferative dynamic environment that relies on constant self-renewal of the intestinal epithelium to maintain homeostasis. Tumor necrosis factor-alpha-induced protein 8 (TNFAIP8 or TIPE0) is a regulator of PI3K-mediated signaling. By binding to PIP2 and PIP3, TIPE family members locally activate PI3K activity while globally inhibiting PI3K activity through sequestration of membranous PIP2. Single-cell RNA sequencing survey of Tipe0-/- small intestine was used to investigate the role of TIPE0 in intestinal differentiation. Tipe0-/- intestinal cells were shown to shift towards an undifferentiated state, with the notable exception of goblet cells. Additionally, three possible novel regulators of terminal cell fate decisions in the secretory lineage were identified: Nupr1, Kdm4a, and Gatad1. We propose that these novel regulators drive changes involved in goblet cell (Nupr1) or tuft cell (Kdm4a and Gatad1) fate commitment and that TIPE0 may play a role in orchestrating terminal differentiation.

Time-course single-cell RNA sequencing reveals transcriptional dynamics and heterogeneity of limbal stem cells derived from human pluripotent stem cells

BACKGROUND

Human pluripotent stem cell-derived limbal stem cells (hPSC-derived LSCs) provide a promising cell source for corneal transplants and ocular surface reconstruction. Although recent efforts in the identification of LSC markers have increased our understanding of the biology of LSCs, much more remains to be characterized in the developmental origin, cell fate determination, and identity of human LSCs. The lack of knowledge hindered the establishment of efficient differentiation protocols for generating hPSC-derived LSCs and held back their clinical application.

RESULTS

Here, we performed a time-course single-cell RNA-seq to investigate transcriptional heterogeneity and expression changes of LSCs derived from human embryonic stem cells (hESCs). Based on current protocol, expression heterogeneity of reported LSC markers were identified in subpopulations of differentiated cells. EMT has been shown to occur during differentiation process, which could possibly result in generation of untargeted cells. Pseudotime trajectory analysis revealed transcriptional changes and signatures of commitment of hESCs-derived LSCs and their progeny-the transit amplifying cells.

CONCLUSION

Single-cell RNA-seq revealed time-course expression changes and significant transcriptional heterogeneity during hESC-derived LSC differentiation in vitro. Our results demonstrated candidate developmental trajectory and several new candidate markers for LSCs, which could facilitate elucidating the identity and developmental origin of human LSCs in vivo.

Advances in Liver Cancer Stem Cell Isolation and their Characterization

Over the last decade research on cancer stem cells (CSC) significantly contributed to a better understanding of tumor biology. Given their similarity to normal stem cells, i.e. self-renewal and pluripotency the need arises to develop robust protocols for the isolation and characterization of CSCs. As with other malignancies, hepatic tumors are composed of a heterogeneous population of cells including liver cancer stem cells (LCSC). Yet, a precise understanding of why stem cells become cancerous is still lacking. There is unmet need to develop robust protocols for the successful isolation of LCSCs from human tissue resection material as to assist in the development of molecular targeted therapies. Here we review the research progress made in the isolation and characterization of LCSCs by considering a wide range of cell surface markers and sorting methods, as applied to side populations, microsphere cultures and the gradient centrifugation method. We emphasize the different fluorescence activated cell sorting methods and the possibility to enrich LCSCs by immunomagnetic beads. We review the specificity of functional assays by considering ABCG transporter and ALDH1 enzyme activities and evaluate the in vivo tumorigenicity of LCSCs in highly sensitive bioassays. Finally, we evaluate different LCSC markers in association with viral and non-viral liver disease and explore the potential of novel drug delivery systems targeting CD133, EpCAM, CD13 and CD90 for the development of molecular targeted therapies. Graphical Abstract.

Novel approaches to liver disease diagnosis and modeling

Lack of a prompt and accurate diagnosis remains on top of the list of challenges faced by patients with rare liver diseases. Although rare liver diseases affect a significant percentage of the population as a group, when taken singularly they represent unique diseases and the approaches used for diagnosis of common liver diseases are insufficient. However, the development of new methods for the acquisition of molecular and clinical data (i.e., genomic, proteomics, metabolomics) and computational tools for their analysis and integration, together with advances in modeling diseases using stem cell-based technology [i.e., induced pluripotent stem cells (iPSCs) and tissue organoids] represent a promising and powerful tool to improve the clinical management of these patients. This is the goal of precision medicine, a novel approach of modern medicine that aims at delivering a specific treatment based on disease-specific biological insights and individual profile. This review will discuss the application and advances of these technologies and how they represent a new opportunity in hepatology.

Single-cell omics analysis reveals functional diversification of hepatocytes during liver regeneration

Adult liver has enormous regenerative capacity; it can regenerate after losing two-thirds of its mass while sustaining essential metabolic functions. How the liver balances dual demands for increased proliferative activity with maintenance of organ function is unknown but essential to prevent liver failure. Using partial hepatectomy (PHx) in mice to model liver regeneration, we integrated single-cell RNA- and ATAC-Seq to map state transitions in approximately 13,000 hepatocytes at single-cell resolution as livers regenerated, and validated key findings with IHC, to uncover how the organ regenerates hepatocytes while simultaneously fulfilling its vital tissue-specific functions. After PHx, hepatocytes rapidly and transiently diversified into multiple distinct populations with distinct functional bifurcation: some retained the chromatin landscapes and transcriptomes of hepatocytes in undamaged adult livers, whereas others transitioned to acquire chromatin landscapes and transcriptomes of fetal hepatocytes. Injury-related signaling pathways known to be critical for regeneration were activated in transitioning hepatocytes, and the most fetal-like hepatocytes exhibited chromatin landscapes that were enriched with transcription factors regulated by those pathways.

Embryonic liver developmental trajectory revealed by single-cell RNA sequencing in the Foxa2eGFP mouse

The liver and gallbladder are among the most important internal organs derived from the endoderm, yet the development of the liver and gallbladder in the early embryonic stages is not fully understood. Using a transgenic Foxa2^eGFP reporter mouse line, we performed single-cell full-length mRNA sequencing on endodermal and hepatic cells isolated from ten embryonic stages, ranging from E7.5 to E15.5. We identified the embryonic liver developmental trajectory from gut endoderm to hepatoblasts and characterized the transcriptome of the hepatic lineage. More importantly, we identified liver primordium as the nascent hepatic progenitors with both gut and liver features and documented dynamic gene expression during the epithelial-hepatic transition (EHT) at the stage of liver specification during E9.5–11.5. We found six groups of genes switched on or off in the EHT process, including diverse transcripitional regulators that had not been previously known to be expressed during EHT. Moreover, we identified and revealed transcriptional profiling of gallbladder primordium at E9.5. The present data provides a high-resolution resource and critical insights for understanding the liver and gallbladder development.

The authors report a single cell-resolution gene expression atlas for the developing mouse liver and gallbladder using a transgenic Foxa2^eGFP mouse line. By tracing the development of cells from gut endoderm to hepatoblasts they identify key transcriptional changes during liver specification.

Causal network perturbations for instance-specific analysis of single cell and disease samples

MOTIVATION

Complex diseases involve perturbation in multiple pathways and a major challenge in clinical genomics is characterizing pathway perturbations in individual samples. This can lead to patient-specific identification of the underlying mechanism of disease thereby improving diagnosis and personalizing treatment. Existing methods rely on external databases to quantify pathway activity scores. This ignores the data dependencies and that pathways are incomplete or condition-specific.

RESULTS

ssNPA is a new approach for subtyping samples based on deregulation of their gene networks. ssNPA learns a causal graph directly from control data. Sample-specific network neighborhood deregulation is quantified via the error incurred in predicting the expression of each gene from its Markov blanket. We evaluate the performance of ssNPA on liver development single-cell RNA-seq data, where the correct cell timing is recovered; and two TCGA datasets, where ssNPA patient clusters have significant survival differences. In all analyses ssNPA consistently outperforms alternative methods, highlighting the advantage of network-based approaches.

AVAILABILITY AND IMPLEMENTATION

http://www.benoslab.pitt.edu/Software/ssnpa/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Single-Nuclei RNA Sequencing Assessment of the Hepatic Effects of 2,3,7,8-Tetrachlorodibenzo-p-dioxin

BACKGROUND AND AIMS

Characterization of cell specific transcriptional responses to hepatotoxicants is lost in the averages of bulk RNA-sequencing (RNA-seq). Single-cell/nuclei RNA-seq technologies enable the transcriptomes of individual cell (sub)types to be assessed within the context of in vivo models.

METHODS

Single-nuclei RNA-sequencing (snSeq) of frozen liver samples from male C57BL/6 mice gavaged with sesame oil vehicle or 30 μg/kg 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) every 4 days for 28 days was used to demonstrate the application of snSeq for the evaluation of xenobiotics.

RESULTS

A total of 19,907 genes were detected across 16,015 nuclei from control and TCDD-treated livers. Eleven cell (sub)types reflected the expected cell diversity of the liver including distinct pericentral, midzonal, and periportal hepatocyte subpopulations. TCDD altered relative proportions of cell types and elicited cell-specific gene expression profiles. For example, macrophages increased from 0.5% to 24.7%, while neutrophils were only present in treated samples, consistent with histological evaluation. The number of differentially expressed genes (DEGs) in each cell type ranged from 122 (cholangiocytes) to 7625 (midcentral hepatocytes), and loosely correlated with the basal expression level of Ahr, the canonical mediator of TCDD and related compounds. In addition to the expected functions within each cell (sub)types, RAS signaling and related pathways were specifically enriched in nonparenchymal cells while metabolic process enrichment occurred primarily in hepatocytes. snSeq also identified the expansion of a Kupffer cell subtype highly expressing Gpnmb, as reported in a dietary NASH model.

CONCLUSIONS

We show that snSeq of frozen liver samples can be used to assess cell-specific transcriptional changes and population shifts in models of hepatotoxicity when examining freshly isolated cells is not feasible.

Investigating higher-order interactions in single-cell data with scHOT

Single-cell genomics has transformed our ability to examine cell fate choice. Examining cells along a computationally ordered 'pseudotime' offers the potential to unpick subtle changes in variability and covariation among key genes. We describe an approach, scHOT-single-cell higher-order testing-which provides a flexible and statistically robust framework for identifying changes in higher-order interactions among genes. scHOT can be applied for cells along a continuous trajectory or across space and accommodates various higher-order measurements including variability or correlation. We demonstrate the use of scHOT by studying coordinated changes in higher-order interactions during embryonic development of the mouse liver. Additionally, scHOT identifies subtle changes in gene-gene correlations across space using spatially resolved transcriptomics data from the mouse olfactory bulb. scHOT meaningfully adds to first-order differential expression testing and provides a framework for interrogating higher-order interactions using single-cell data.

Comparative analysis of cell lineage differentiation during hepatogenesis in humans and mice at the single-cell transcriptome level

During embryogenesis, the liver is the site of hepatogenesis and hematopoiesis and contains many cell lineages derived from the endoderm and mesoderm. However, the characteristics and developmental programs of many of these cell lineages remain unclear, especially in humans. Here, we performed single-cell RNA sequencing of whole human and mouse fetal livers throughout development. We identified four cell lineage families of endoderm-derived, erythroid, non-erythroid hematopoietic, and mesoderm-derived non-hematopoietic cells, and defined the developmental pathways of the major cell lineage families. In both humans and mice, we identified novel markers of hepatic lineages and an ID3⁺ subpopulation of hepatoblasts as well as verified that hepatoblast differentiation follows the “default-directed” model. Additionally, we found that human but not mouse fetal hepatocytes display heterogeneity associated with expression of metabolism-related genes. We described the developmental process of erythroid progenitor cells during human and mouse hematopoiesis. Moreover, despite the general conservation of cell differentiation programs between species, we observed different cell lineage compositions during hematopoiesis in the human and mouse fetal livers. Taken together, these results reveal the dynamic cell landscape of fetal liver development and illustrate the similarities and differences in liver development between species, providing an extensive resource for inducing various liver cell lineages in vitro.

Erythroid Lineage Cells in the Liver: Novel Immune Regulators and Beyond

The lineage of the erythroid cell has been revisited in recent years. Instead of being classified as simply inert oxygen carriers, emerging evidence has shown that they are a tightly regulated in immune potent population with potential developmental plasticity for lineage crossing. Erythroid cells have been reported to exert immune regulatory function through secreted cytokines, or cell-cell contact, depending on the conditions of the microenvironment and disease models. In this review, we explain the natural history of erythroid cells in the liver through a developmental lens, as it offers perspectives into newly recognized roles of this lineage in liver biology. Here, we review the known immune roles of erythroid cells and discuss the mechanisms in the context of disease models and stages. Then, we explore the capability of erythroid lineage as a cell source for regenerative medicine. We propose that the versatile lineage of erythroid cells provides an underappreciated and potentially promising area for basic and translational research in the field of liver disease.

Advances of single-cell genomics and epigenomics in human disease: where are we now?

Cellular heterogeneity is revolutionizing the way to study, monitor and dissect complex diseases. This has been possible with the technological and computational advances associated to single-cell genomics and epigenomics. Deeper understanding of cell-to-cell variation and its impact on tissue function will open new avenues for early disease detection, accurate diagnosis and personalized treatments, all together leading to the next generation of health care. This review focuses on the recent discoveries that single-cell genomics and epigenomics have facilitated in the context of human health. It highlights the potential of single-cell omics to further advance the development of personalized treatments and precision medicine in cancer, diabetes and chronic age-related diseases. The promise of single-cell technologies to generate new insights about the differences in function between individual cells is just emerging, and it is paving the way for identifying biomarkers and novel therapeutic targets to tackle age, complex diseases and understand the effect of life style interventions and environmental factors.

Single-Cell Transcriptome Analysis of Mouse Liver Cell-Specific Tropism and Transcriptional Dysregulation Following Intravenous Administration of AAVrh.10 Vectors

Capitalizing on liver tropism of adeno-associated viral (AAV) vectors, intravenous vector administration is commonly used to genetically modify hepatocytes, a strategy currently in clinical trials for a number of liver-based hereditary disorders. Although hepatocytes are known to exhibit extensive phenotypic heterogeneity influenced by liver zonation and dietary cycle, there is little data available for the tropism capacity, as well as the potential transcriptional dysregulation, of AAV vectors for specific liver cell types. To assess these issues, we employed single-cell RNA sequencing of the mouse liver after intravenous administration of the liver tropic AAVrh.10 vector to characterize cell-specific AAV-mediated transgene expression and transcriptome dysregulation. Wild-type 8-week-old male C57Bl/6 mice under normal feed cycle were randomly divided into three groups and intravenously administered phosphate-buffered saline (PBS), AAVrh.10Null (no transgene), or AAVrh.10mCherry (marker gene). Overall, a total of 46,500 liver cells were sequenced. The single-cell transcriptomic profiles were grouped into three separate clusters of hepatocytes (Ttr-enriched "Hep1," Tat-enriched "Hep2," and Alb-enriched "Hep3") and multiple other cell types. The hepatocyte diversity was driven by glucose and lipid homeostasis signaling. Assessment of the transgene expression demonstrated that AAVrh.10 is primarily Hep1-tropic, with a 10-gene signature positively correlated with AAVrh.10-mediated transgene expression. The transgene expression was less in Hep2 and Hep3 cells with a high receptor tyrosine kinase phenotype. Importantly, AAVrh.10 vector interactions with the liver markedly altered the transcriptional patterns of all cell types, with modified genes enriched in pathways of complement and coagulation cascade, cytochrome P450, peroxisome, antigen processing and presentation, and endoplasmic reticulum protein processing. These observations provide insights into the liver cell-specific consequences of AAV-mediated liver gene transfer, far beyond the well-known organ-specific expression of the vector-delivered transgene.

Mutational and transcriptomic landscapes of a rare human prostate basal cell carcinoma

BACKGROUND

As a rare subtype of prostate carcinoma, basal cell carcinoma (BCC) has not been studied extensively and thus lacks systematic molecular characterization.

METHODS

Here, we applied single-cell genomic amplification and RNA-Seq to a specimen of human prostate BCC (CK34βE12⁺ /P63⁺ /PAP^- /PSA^- ). The mutational landscape was obtained via whole exome sequencing of the amplification mixture of 49 single cells, and the transcriptomes of 69 single cells were also obtained.

RESULTS

The five putative driver genes mutated in BCC are CASC5, NUTM1, PTPRC, KMT2C, and TBX3, and the top three nucleotide substitutions are C>T, T>C, and C>A, similar to common prostate cancer. The distribution of the variant allele frequency values indicated that these single cells are from the same tumor clone. The 69 single cells were clustered into tumor, stromal, and immune cells based on their global transcriptomic profiles. The tumor cells specifically express basal cell markers like KRT5, KRT14, and KRT23 and epithelial markers EPCAM, CDH1, and CD24. The transcription factor covariance network analysis showed that the BCC tumor cells have distinct regulatory networks. By comparison with current prostate cancer datasets, we found that some of the bulk samples exhibit basal cell signatures. Interestingly, at single-cell resolution the gene expression patterns of prostate BCC tumor cells show uniqueness compared with that of common prostate cancer-derived circulating tumor cells.

CONCLUSIONS

This study, for the first time, discloses the comprehensive mutational and transcriptomic landscapes of prostate BCC, which lays a foundation for the understanding of its tumorigenesis mechanism and provides new insights into prostate cancers in general.

A time-resolved multi-omic atlas of the developing mouse liver

Liver organogenesis and development are composed of a series of complex, well-orchestrated events. Identifying key factors and pathways governing liver development will help elucidate the physiological and pathological processes including those of cancer. We conducted multidimensional omics measurements including protein, mRNA, and transcription factor (TF) DNA-binding activity for mouse liver tissues collected from embryonic day 12.5 (E12.5) to postnatal week 8 (W8), encompassing major developmental stages. These data sets reveal dynamic changes of core liver functions and canonical signaling pathways governing development at both mRNA and protein levels. The TF DNA-binding activity data set highlights the importance of TF activity in early embryonic development. A comparison between mouse liver development and human hepatocellular carcinoma (HCC) proteomic profiles reveal that more aggressive tumors are characterized with the activation of early embryonic development pathways, whereas less aggressive ones maintain liver function-related pathways that are elevated in the mature liver. This work offers a panoramic view of mouse liver development and provides a rich resource to explore in-depth functional characterization.

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

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

scReClassify: post hoc cell type classification of single-cell rNA-seq data

Background

Single-cell RNA-sequencing (scRNA-seq) is a fast emerging technology allowing global transcriptome profiling on the single cell level. Cell type identification from scRNA-seq data is a critical task in a variety of research such as developmental biology, cell reprogramming, and cancers. Typically, cell type identification relies on human inspection using a combination of prior biological knowledge (e.g. marker genes and morphology) and computational techniques (e.g. PCA and clustering). Due to the incompleteness of our current knowledge and the subjectivity involved in this process, a small amount of cells may be subject to mislabelling.

Results

Here, we propose a semi-supervised learning framework, named scReClassify, for ‘post hoc’ cell type identification from scRNA-seq datasets. Starting from an initial cell type annotation with potentially mislabelled cells, scReClassify first performs dimension reduction using PCA and next applies a semi-supervised learning method to learn and subsequently reclassify cells that are likely mislabelled initially to the most probable cell types. By using both simulated and real-world experimental datasets that profiled various tissues and biological systems, we demonstrate that scReClassify is able to accurately identify and reclassify misclassified cells to their correct cell types.

Conclusions

scReClassify can be used for scRNA-seq data as a post hoc cell type classification tool to fine-tune cell type annotations generated by any cell type classification procedure. It is implemented as an R package and is freely available from https://github.com/SydneyBioX/scReClassify

Single-Cell RNA Sequencing-Based Computational Analysis to Describe Disease Heterogeneity

The trillions of cells in the human body can be viewed as elementary but essential biological units that achieve different body states, but the low resolution of previous cell isolation and measurement approaches limits our understanding of the cell-specific molecular profiles. The recent establishment and rapid growth of single-cell sequencing technology has facilitated the identification of molecular profiles of heterogeneous cells, especially on the transcription level of single cells [single-cell RNA sequencing (scRNA-seq)]. As a novel method, the robustness of scRNA-seq under changing conditions will determine its practical potential in major research programs and clinical applications. In this review, we first briefly presented the scRNA-seq-related methods from the point of view of experiments and computation. Then, we compared several state-of-the-art scRNA-seq analysis frameworks mainly by analyzing their performance robustness on independent scRNA-seq datasets for the same complex disease. Finally, we elaborated on our hypothesis on consensus scRNA-seq analysis and summarized the potential indicative and predictive roles of individual cells in understanding disease heterogeneity by single-cell technologies.

Lgr5+ stem and progenitor cells reside at the apex of a heterogeneous embryonic hepatoblast pool

During mouse embryogenesis, progenitors within the liver known as hepatoblasts give rise to adult hepatocytes and cholangiocytes. Hepatoblasts, which are specified at E8.5-E9.0, have been regarded as a homogeneous progenitor population that initiate differentiation from E13.5. Recently, scRNA-seq analysis has identified sub-populations of transcriptionally distinct hepatoblasts at E11.5. Here, we show that hepatoblasts are not only transcriptionally but also functionally heterogeneous, and that a subpopulation of E9.5-E10.0 hepatoblasts exhibit a previously unidentified early commitment to cholangiocyte fate. Importantly, we also identify a subpopulation constituting 2% of E9.5-E10.0 hepatoblasts that express the adult stem cell marker Lgr5, and generate both hepatocyte and cholangiocyte progeny that persist for the lifespan of the mouse. Combining lineage tracing and scRNA-seq, we show that Lgr5 marks E9.5-E10.0 bipotent liver progenitors residing at the apex of a hepatoblast hierarchy. Furthermore, isolated Lgr5⁺ hepatoblasts can be clonally expanded in vitro into embryonic liver organoids, which can commit to either hepatocyte or cholangiocyte fates. Our study demonstrates functional heterogeneity within E9.5 hepatoblasts and identifies Lgr5 as a marker for a subpopulation of bipotent liver progenitors.

scMerge leverages factor analysis, stable expression, and pseudoreplication to merge multiple single-cell RNA-seq datasets

Concerted examination of multiple collections of single-cell RNA sequencing (RNA-seq) data promises further biological insights that cannot be uncovered with individual datasets. Here we present scMerge, an algorithm that integrates multiple single-cell RNA-seq datasets using factor analysis of stably expressed genes and pseudoreplicates across datasets. Using a large collection of public datasets, we benchmark scMerge against published methods and demonstrate that it consistently provides improved cell type separation by removing unwanted factors; scMerge can also enhance biological discovery through robust data integration, which we show through the inference of development trajectory in a liver dataset collection.

Application of qPCR assays based on haloacids transporter gene dehp2 for discrimination of Burkholderia and Paraburkholderia

Background

A major facilitator superfamily transporter Dehp2 was recently shown to be playing an important role in transport and biodegradation of haloacids in Paraburkholderia caribensis MBA4, and Dehp2 is phylogenetically conserved in Burkholderia sensu lato.

Results

We designed both Burkholderia sensu stricto-specific and Paraburkholderia-specific qPCR assays based on dehp2 and 16S rRNA, and validated the qPCR assays in 12 bacterial strains. The qPCR assays could detect single species of Burkholderia sensu stricto or Paraburkholderia with high sensitivity and discriminate them in mixtures with high specificity over a wide dynamic range of relative concentrations. At relatively lower cost compared with sequencing-based approach, the qPCR assays will facilitate discrimination of Burkholderia sensu stricto and Paraburkholderia in a large number of samples.

Conclusions

For the first time, we report the utilization of a haloacids transporter gene for discriminative purpose in Burkholderia sensu lato. This enables not only quick decision on proper handling of putative pathogenic samples in Burkholderia sensu stricto group but also future exploitation of relevant species in Paraburkholderia group for haloacids biodegradation purposes.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1411-0) contains supplementary material, which is available to authorized users.