A human liver cell atlas reveals heterogeneity and epithelial progenitors

Achieving a Deeper Understanding of Drug Metabolism and Responses Using Single-Cell Technologies

Recent advancements in single-cell technologies have enabled detection of RNA, proteins, metabolites, and xenobiotics in individual cells, and the application of these technologies has the potential to transform pharmacological research. Single-cell data has already resulted in the development of human and model species cell atlases, identifying different cell types within a tissue, further facilitating the characterization of tumor heterogeneity, and providing insight into treatment resistance. Research discussed in this review demonstrates that distinct cell populations express drug metabolizing enzymes to different extents, indicating there may be variability in drug metabolism not only between organs, but within tissue types. Additionally, we put forth the concept that single-cell analyses can be used to expose underlying variability in cellular response to drugs, providing a unique examination of drug efficacy, toxicity, and metabolism. We will outline several of these techniques: single-cell RNA-sequencing and mass cytometry to characterize and distinguish different cell types, single-cell proteomics to quantify drug metabolizing enzymes and characterize cellular responses to drug, capillary electrophoresis-ultrasensitive laser-induced fluorescence detection and single-probe single-cell mass spectrometry for detection of drugs, and others. Emerging single-cell technologies such as these can comprehensively characterize heterogeneity in both cell-type-specific drug metabolism and response to treatment, enhancing progress toward personalized and precision medicine. SIGNIFICANCE STATEMENT: Recent technological advances have enabled the analysis of gene expression and protein levels in single cells. These types of analyses are important to investigating mechanisms that cannot be elucidated on a bulk level, primarily due to the variability of cell populations within biological systems. Here, we summarize cell-type-specific drug metabolism and how pharmacologists can utilize single-cell approaches to obtain a comprehensive understanding of drug metabolism and cellular heterogeneity in response to drugs.

To be or not to be: The double-edged sword roles of liver progenitor cells

Given the liver's remarkable and unique regenerative capacity, researchers have long focused on liver progenitor cells (LPCs) and liver cancer stem cells (LCSCs). LPCs can differentiate into both hepatocytes and cholangiocytes. However, the mechanism underlying cell conversion and its distinct contribution to liver homeostasis and tumorigenesis remain unclear. In this review, we discuss the complicated conversions involving LPCs and LCSCs. As the critical intermediate state in malignant transformation, LPCs play double-edged sword roles. LPCs are not only involved in hepatic wound-healing responses by supplementing liver cells and bile duct cells in the damaged liver but may transform into LCSCs under dysregulation of key signaling pathways, resulting in refractory malignant liver tumors. Because LPC lineages are temporally and spatially dynamic, we discuss crucial LPC subgroups and summarize regulatory factors correlating with the trajectories of LPCs and LCSCs in the liver tumor microenvironment. This review elaborates on the double-edged sword roles of LPCs to help understand the liver's regenerative potential and tumor heterogeneity. Understanding the sources and transformations of LPCs is essential in determining how to exploit their regenerative capacity in the future.

Single-Cell RNA Transcriptome Profiling of Liver Cells of Short-Term Alcoholic Liver Injury in Mice

Alcoholic liver disease (ALD) is currently considered a global healthcare problem with limited pharmacological treatment options. There are abundant cell types in the liver, such as hepatocytes, endothelial cells, Kupffer cells and so on, but little is known about which kind of liver cells play the most important role in the process of ALD. To obtain a cellular resolution of alcoholic liver injury pathogenesis, 51,619 liver single-cell transcriptomes (scRNA-seq) with different alcohol consumption durations were investigated, 12 liver cell types were identified, and the cellular and molecular mechanisms of the alcoholic liver injury were revealed. We found that more aberrantly differential expressed genes (DEGs) were present in hepatocytes, endothelial cells, and Kupffer cells than in other cell types in alcoholic treatment mice. Alcohol promoted the pathological processes of liver injury; the specific mechanisms involved: lipid metabolism, oxidative stress, hypoxia, complementation and anticoagulation, and hepatocyte energy metabolism on hepatocytes; NO production, immune regulation, epithelial and cell migration on endothelial cells; antigen presentation and energy metabolism on Kupffer cells, based on the GO analysis. In addition, our results showed that some transcription factors (TFs) are activated in alcohol-treated mice. In conclusion, our study improves the understanding of liver cell heterogeneity in alcohol-fed mice at the single-cell level. It has potential value for understanding key molecular mechanisms and improving current prevention and treatment strategies for short-term alcoholic liver injury.

The native liver as inspiration to create superior in vitro hepatic models

Drug induced liver injury (DILI) is one of the major reasons of drug withdrawal during the different phases of drug development. The later in the drug development a drug is discovered to be toxic, the higher the economical as well as the ethical impact will be. In vitro models for early detection of drug liver toxicity are under constant development, however to date a superior model of the liver is still lacking. Ideally, a highly reliable model should be established to maintain the different hepatic cell functionalities to the greatest extent possible, during a period of time long enough to allow for tracking of the toxicity of compounds. In the case of DILI, toxicity can appear even after months of exposure. To reach this goal, an in vitro model should be developed that mimics the in vivo liver environment, function and response to external stimuli. The different approaches for the development of liver models currently used in the field of tissue engineering will be described in this review. Combining different technologies, leading to optimal materials, cells and 3D-constructs will ultimately lead to an ideal superior model that fully recapitulates the 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.

The co-inhibitory receptor TIGIT regulates NK cell function and is upregulated in human intrahepatic CD56bright NK cells

The crosstalk between NK cells and their surrounding environment is enabled through activating and inhibitory receptors, which tightly control NK cell activity. The co-inhibitory receptor TIGIT decreases NK cell cytotoxicity and is involved in NK cell exhaustion, but has also been associated with liver regeneration, highlighting that the contribution of human intrahepatic CD56bright NK cells in regulating tissue homeostasis remains incompletely understood. A targeted single-cell mRNA analysis revealed distinct transcriptional differences between matched human peripheral blood and intrahepatic CD56bright NK cells. Multiparameter flow cytometry identified a cluster of intrahepatic NK cells with overlapping high expression of CD56, CD69, CXCR6, TIGIT and CD96. Intrahepatic CD56bright NK cells also expressed significantly higher protein surface levels of TIGIT, and significantly lower levels of DNAM-1 compared to matched peripheral blood CD56bright NK cells. TIGIT+ CD56bright NK cells showed diminished degranulation and TNF-α production following stimulation. Co-incubation of peripheral blood CD56bright NK cells with human hepatoma cells or primary human hepatocyte organoids resulted in migration of NK cells into hepatocyte organoids and upregulation of TIGIT and downregulation of DNAM-1 expression, in line with the phenotype of intrahepatic CD56bright NK cells. Intrahepatic CD56bright NK cells represent a transcriptionally, phenotypically, and functionally distinct population of NK cells that expresses higher levels of TIGIT and lower levels of DNAM-1 than matched peripheral blood CD56bright NK cells. Increased expression of inhibitory receptors by NK cells within the liver environment can contribute to tissue homeostasis and reduction of liver inflammation.

CD203c is expressed by human fetal hepatoblasts and distinguishes subsets of hepatoblastoma

Background & Aims

Hepatocytic cells found during prenatal development have unique features compared to their adult counterparts, and are believed to be the precursors of pediatric hepatoblastoma. The cell-surface phenotype of hepatoblasts and hepatoblastoma cell lines was evaluated to discover new markers of these cells and gain insight into the development of hepatocytic cells and the phenotypes and origins of hepatoblastoma.

Methods

Human midgestation livers and four pediatric hepatoblastoma cell lines were screened using flow cytometry. Expression of over 300 antigens was evaluated on hepatoblasts defined by their expression of CD326 (EpCAM) and CD14. Also analyzed were hematopoietic cells, expressing CD45, and liver sinusoidal-endothelial cells (LSECs), expressing CD14 but lacking CD45 expression. Select antigens were further examined by fluorescence immunomicroscopy of fetal liver sections. Antigen expression was also confirmed on cultured cells by both methods. Gene expression analysis by liver cells, 6 hepatoblastoma cell lines, and hepatoblastoma cells was performed. Immunohistochemistry was used to evaluate CD203c, CD326, and cytokeratin-19 expression on three hepatoblastoma tumors.

Results

Antibody screening identified many cell surface markers commonly or divergently expressed by hematopoietic cells, LSECs, and hepatoblasts. Thirteen novel markers expressed on fetal hepatoblasts were identified including ectonucleotide pyrophosphatase/phosphodiesterase family member 3 (ENPP-3/CD203c), which was found to be expressed by hepatoblasts with widespread expression in the parenchyma of the fetal liver. In culture CD203c⁺CD326⁺⁺ cells resembled hepatocytic cells with coexpression of albumin and cytokeratin-19 confirming a hepatoblast phenotype. CD203c expression declined rapidly in culture whereas the loss of CD326 was not as pronounced. CD203c and CD326 were co-expressed on a subset of hepatoblastoma cell lines and hepatoblastomas with an embryonal pattern.

Conclusions

CD203c is expressed on hepatoblasts and may play a role in purinergic signaling in the developing liver. Hepatoblastoma cell lines were found to consist of two broad phenotypes consisting of a cholangiocyte-like phenotype that expressed CD203c and CD326 and a hepatocyte-like phenotype with diminished expression of these markers. CD203c was expressed by some hepatoblastoma tumors and may represent a marker of a less differentiated embryonal component.

On triangle inequalities of correlation-based distances for gene expression profiles

BACKGROUND

Distance functions are fundamental for evaluating the differences between gene expression profiles. Such a function would output a low value if the profiles are strongly correlated-either negatively or positively-and vice versa. One popular distance function is the absolute correlation distance, [Formula: see text], where [Formula: see text] is similarity measure, such as Pearson or Spearman correlation. However, the absolute correlation distance fails to fulfill the triangle inequality, which would have guaranteed better performance at vector quantization, allowed fast data localization, as well as accelerated data clustering.

RESULTS

In this work, we propose [Formula: see text] as an alternative. We prove that [Formula: see text] satisfies the triangle inequality when [Formula: see text] represents Pearson correlation, Spearman correlation, or Cosine similarity. We show [Formula: see text] to be better than [Formula: see text], another variant of [Formula: see text] that satisfies the triangle inequality, both analytically as well as experimentally. We empirically compared [Formula: see text] with [Formula: see text] in gene clustering and sample clustering experiment by real-world biological data. The two distances performed similarly in both gene clustering and sample clustering in hierarchical clustering and PAM (partitioning around medoids) clustering. However, [Formula: see text] demonstrated more robust clustering. According to the bootstrap experiment, [Formula: see text] generated more robust sample pair partition more frequently (P-value [Formula: see text]). The statistics on the time a class "dissolved" also support the advantage of [Formula: see text] in robustness.

CONCLUSION

[Formula: see text], as a variant of absolute correlation distance, satisfies the triangle inequality and is capable for more robust clustering.

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

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

ScRNA-seq and ST-seq in liver research

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

Modeling Nonalcoholic Fatty Liver Disease in the Dish Using Human-Specific Platforms: Strategies and Limitations

Nonalcoholic fatty liver disease (NAFLD) is a chronic liver disease affecting multiple cell types of the human liver. The high prevalence of NAFLD and the lack of approved therapies increase the demand for reliable models for the preclinical discovery of drug targets. In the last decade, multiple proof-of-principle studies have demonstrated human-specific NAFLD modeling in the dish. These systems have included technologies based on human induced pluripotent stem cell derivatives, liver tissue section cultures, intrahepatic cholangiocyte organoids, and liver-on-a-chip. These platforms differ in functional maturity, multicellularity, scalability, and spatial organization. Identifying an appropriate model for a specific NAFLD-related research question is challenging. Therefore, we review different platforms for their strengths and limitations in modeling NAFLD. To define the fidelity of the current human in vitro NAFLD models in depth, we define disease hallmarks within the NAFLD spectrum that range from steatosis to severe fibroinflammatory tissue injury. We discuss how the most common methods are efficacious in modeling genetic contributions and aspects of the early NAFLD-related tissue response. We also highlight the shortcoming of current models to recapitulate the complexity of inter-organ crosstalk and the chronic process of liver fibrosis-to-cirrhosis that usually takes decades in patients. Importantly, we provide methodological overviews and discuss implementation hurdles (eg, reproducibility or costs) to help choose the most appropriate NAFLD model for the individual research focus: hepatocyte injury, ductular reaction, cellular crosstalk, or other applications. In sum, we highlight current strategies and deficiencies to model NAFLD in the dish and propose a framework for the next generation of human-specific investigations.

A new glutamine synthetase index to evaluate hepatic lobular restoration in advanced fibrosis during anti-HBV therapy

Hepatic lobular architecture distortion is a deleterious turning point and a crucial histological feature of advanced liver fibrosis in chronic liver diseases. Regression of fibrosis has been documented in chronic hepatitis B (CHB) patients. However, whether lobular architecture could be restored following fibrosis regression after antiviral therapy is still unclear. Glutamine synthetase (GS) is generally expressed by perivenular hepatocytes around hepatic veins (HV). In this study, we defined abnormal lobular architecture (GSPT ) as GS expressing in the vicinity of portal tracts (PT), which denotes parenchymal extinction and lobular collapse. We defined normal lobular architecture (GSHV ) as GS positivity area not approximating PTs. Therefore, we propose a new GS-index, defined as the percentage of GSHV /(GSHV  + GSPT ), to evaluate the extent of architectural disruption and restoration. We evaluated 43 CHB patients with advanced fibrosis (Ishak stage ≥4). Posttreatment liver biopsy was performed after 78 weeks of anti-HBV therapy. The median GS-index improved from 7% (interquartile range [IQR]: 0%-23%) at baseline to 36% (IQR: 20%-57%) at Week 78 (p < 0.001). Totals of 22 patients (51%) had significant GS-index improvement from 0% (IQR: 0%-13%) to 55% (IQR: 44%-81%), while the other half had almost no change between 17% (IQR: 0%-33%) to 20% (IQR: 12%-31%). When GS-index78w  ≥ 50% was used to define hepatic lobular restoration, 37% of patients (16/43) achieved lobular restoration, with much improvement in alanine aminotransferase (ALT) and aspartate aminotransferase (AST) levels (median value of ∆/Baseline in ALT: restored vs. nonrestored was 79.1% vs. 48.8%, p = 0.018; median value of ∆/Baseline in AST: restored vs. nonrestored was 69.1% vs. 32.5%, p = 0.005). More importantly, lobular restoration correlated with fibrosis regression (median value of ∆/Baseline in Ishak stage: restored vs. nonrestored was 25.0% vs. 0%, p = 0.008). Therefore, in the era of antiviral therapy for CHB, restoration of hepatic lobular architecture is achievable in patients with advanced fibrosis. GS-index provides additional insight into fibrosis regression that goes beyond collagen degradation.

Recent Advances in Intrahepatic Biliary Epithelial Heterogeneity

Biliary epithelium (i.e., cholangiocytes) is a heterogeneous population of epithelial cells in the liver, which line small and large bile ducts and have individual responses and functions dependent on size and location in the biliary tract. We discuss the recent findings showing that the intrahepatic biliary tree is heterogeneous regarding (1) morphology and function, (2) hormone expression and signaling (3), response to injury, and (4) roles in liver regeneration. This review overviews the significant characteristics and differences of the small and large cholangiocytes. Briefly, it outlines the in vitro and in vivo models used in the heterogeneity evaluation. In conclusion, future studies addressing biliary heterogeneity's role in the pathogenesis of liver diseases characterized by ductular reaction may reveal novel therapeutic approaches.

HtrA2 Independently Predicts Poor Prognosis and Correlates with Immune Cell Infiltration in Hepatocellular Carcinoma

High-temperature requirement protein A2 (HtrA2), a mitochondrial protein, is related to apoptosis regulation. However, the role of HtrA2 in hepatocellular carcinoma (HCC) remains unclear. In the present study, we explored the prognostic value and expression pattern of HtrA2 in HCC and confirmed its independent value for predicting outcomes via Cox analyses. LinkedOmics and GEPIA2 were used to construct the coexpression and functional networks of HtrA2. Additionally, the data obtained from TCGA was analyzed to investigate the relationship between the infiltration of immune cells and HtrA2 mRNA expression. Finally, the expression pattern of HtrA2 in HCC was confirmed by wet-lab experiments. The results showed high HtrA2 expression (P < 0.001) presented in tumor tissues in TCGA-HCC. Moreover, high HtrA2 expression was confirmed to be associated with poor HCC patient survival (P < 0.05). HtrA2 has also been recognized as an essential risk factor for overall survival (P=0.01, HR = 1.654, 95% CI 1.128-2.425), disease-specific survival (P=0.004, HR = 2.204, 95% CI 1.294-3.753), and progression-free interval (P=0.007, HR = 1.637, 95% CI 1.145-2.341) of HCC. HCC patients with low HtrA2 methylation had worse overall survival than patients with high methylation (P=0.0019). Functional network analysis suggests that HtrA2 regulates mitochondrial homeostasis through pathways involving multiple microRNAs and transcription factors in HCC. In addition, HtrA2 expression correlated with infiltrating levels of multiple immune cell populations. At last, increased expression of HtrA2 in HCC was confirmed using wet-lab experiments. Our study provides evidence that the upregulation of HtrA2 in HCC is an independent predictor of prognosis. Our results provide the foundation for further study on the roles of HtrA2 in HCC tumorigenesis.

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

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

Multicytokine-producing CD4+ T cells characterize the livers of patients with NASH

A role of CD4+ T cells during the progression from nonalcoholic fatty liver disease (NAFLD) to nonalcoholic steatohepatitis (NASH) has been suggested, but which polarization state of these cells characterizes this progression and the development of fibrosis remain unclear. In addition, a gut-liver axis has been suggested to play a role in NASH, but the role of CD4+ T cells in this axis has just begun to be investigated. Combining single-cell RNA sequencing and multiple-parameter flow cytometry, we provide the first cell atlas to our knowledge focused on liver-infiltrating CD4+ T cells in patients with NAFLD and NASH, showing that NASH is characterized by a population of multicytokine-producing CD4+ T cells. Among these cells, only those with a Th17 polarization state were enriched in patients with advanced fibrosis. In parallel, we observed that Bacteroides appeared to be enriched in the intestine of NASH patients and to correlate with the frequency of multicytokine-producing CD4+ T cells. In short, we deliver a CD4+ T cell atlas of NAFLD and NASH, providing the rationale to target CD4+ T cells with a Th17 polarization state to block fibrosis development. Finally, our data offer an early indication to test whether multicytokine-producing CD4+ T cells are part of the gut-liver axis characterizing NASH.

Single-cell RNA sequencing reveals intrahepatic and peripheral immune characteristics related to disease phases in HBV-infected patients

OBJECTIVE

A comprehensive immune landscape for HBV infection is pivotal to achieve HBV cure.

DESIGN

We performed single-cell RNA sequencing of 2 43 000 cells from 46 paired liver and blood samples of 23 individuals, including six immune tolerant, 5 immune active (IA), 3 acute recovery (AR), 3 chronic resolved and 6 HBV-free healthy controls (HCs). Flow cytometry and histological assays were applied in a second HBV cohort for validation.

RESULTS

Both IA and AR were characterised by high levels of intrahepatic exhausted CD8+ T (Tex) cells. In IA, Tex cells were mainly derived from liver-resident GZMK+ effector memory T cells and self-expansion. By contrast, peripheral CX3CR1+ effector T cells and GZMK+ effector memory T cells were the main source of Tex cells in AR. In IA but not AR, significant cell-cell interactions were observed between Tex cells and regulatory CD4+ T cells, as well as between Tex and FCGR3A+ macrophages. Such interactions were potentially mediated through human leukocyte antigen class I molecules together with their receptors CANX and LILRBs, respectively, contributing to the dysfunction of antiviral immune responses. By contrast, CX3CR1+GNLY+ central memory CD8+ T cells were concurrently expanded in both liver and blood of AR, providing a potential surrogate marker for viral resolution. In clinic, intrahepatic Tex cells were positively correlated with serum alanine aminotransferase levels and histological grading scores.

CONCLUSION

Our study dissects the coordinated immune responses for different HBV infection phases and provides a rich resource for fully understanding immunopathogenesis and developing effective therapeutic strategies.

A new and rare type of hepatocellular carcinoma: Survival and gene analysis of portal vein tumour thrombus-type hepatocellular carcinoma

BACKGROUND

Portal vein tumour thrombus (PVTT) in patients with hepatocellular carcinoma (HCC) is known as a major complication associated with poor survival. We clinically defined a new and rare type of HCC, PVTT-type HCC (PVTT-HCC), in a small group of HCC patients with HCC presenting only as PVTT without a demonstrable parenchyma tumour. The clinicopathological and biological features of PVTT-HCC are not clear.

METHODS

The data for patients who had PVTT-HCC with histologically confirmed HCC from January 2004 to December 2012 at the Eastern Hepatobiliary Surgery Hospital were retrospectively analysed. The survival outcomes of patients with PVTT-HCC were compared with those of HCC patients with PVTT (HCC-PVTT). Propensity score matching (PSM) analysis was performed to match patients at a ratio of 1:3. Then, we performed RNA-Seq analysis of liver samples from PVTT-HCC and HCC-PVTT patients to identify and compare differentially expressed genes and biological pathways between the two groups.

RESULTS

We observed and collected 10 rare cases of PVTT-HCC and performed a prospective cohort study to compare overall survival (OS) between PVTT-HCC and HCC-PVTT. PVTT invaded the main portal vein in 10 PVTT-HCC patients. Univariate and multivariate analyses demonstrated that ChildPugh (A/B), different treatments (LR/non-LR), and different groups were independent risk factors for OS. The median OS was 10.3 months (95 % CI = 6.7-13.8) in the HCC-PVTT group and 7.5 months (95 % CI = 2.8-12.1) in the PVTT-HCC group (P = 0.042). From RNA-Seq, 1630 differentially expressed genes were obtained, of which 731 were upregulated and 899 downregulated in PVTT-HCC compared with HCC-PVTT.

CONCLUSIONS

The survival outcomes of patients with PVTT-HCC were worse than those of patients with HCC-PVTT. RNA-Seq demonstrated differential gene expression between PVTT-HCC and HCC-PVTT, indicating that the former may have distinguishing biological characteristics and be a new and rare type of HCC.

Secretin alleviates biliary and liver injury during late-stage primary biliary cholangitis via restoration of secretory processes

BACKGROUND & AIMS

Primary biliary cholangitis (PBC) is characterised by ductopenia, ductular reaction, impairment of anion exchanger 2 (AE2) and the 'bicarbonate umbrella'. Ductulo-canalicular junction (DCJ) derangement is hypothesised to promote PBC progression. The secretin (Sct)/secretin receptor (SR) axis regulates cystic fibrosis transmembrane receptor (CFTR) and AE2, thus promoting choleresis. We evaluated the role of Sct/SR signalling on biliary secretory processes and subsequent injury in a late-stage PBC mouse model and human samples.

METHODS

At 32 weeks of age, female and male wild-type and dominant-negative transforming growth factor beta receptor II (late-stage PBC model) mice were treated with Sct for 1 or 8 weeks. Bulk RNA-sequencing was performed in isolated cholangiocytes from mouse models.

RESULTS

Biliary Sct/SR/CFTR/AE2 expression and bile bicarbonate levels were reduced in late-stage PBC mouse models and human samples. Sct treatment decreased bile duct loss, ductular reaction, inflammation, and fibrosis in late-stage PBC models. Sct reduced hepatic bile acid levels, modified bile acid composition, and restored the DCJ and 'bicarbonate umbrella'. RNA-sequencing identified that Sct promoted mature epithelial marker expression, specifically anterior grade protein 2 (Agr2). Late-stage PBC models and human samples exhibited reduced biliary mucin 1 levels, which were enhanced by Sct treatment.

CONCLUSION

Loss of Sct/SR signalling in late-stage PBC results in a faulty 'bicarbonate umbrella' and reduced Agr2-mediated mucin production. Sct restores cholangiocyte secretory processes and DCJ formation through enhanced mature cholangiocyte phenotypes and bile duct growth. Sct treatment may be beneficial for individuals with late-stage PBC.

IMPACT AND IMPLICATIONS

Secretin (Sct) regulates biliary proliferation and bicarbonate secretion in cholangiocytes via its receptor, SR, and in mouse models and human samples of late-stage primary biliary cholangitis (PBC), the Sct/SR axis is blunted along with loss of the protective 'bicarbonate umbrella'. We found that both short- and long-term Sct treatment ameliorated ductular reaction, immune cell influx, and liver fibrosis in late-stage PBC mouse models. Importantly, Sct treatment promoted bicarbonate and mucin secretion and hepatic bile acid efflux, thus reducing cholestatic and toxic bile acid-associated injury in late-stage PBC mouse models. Our work perpetuates the hypothesis that PBC pathogenesis hinges on secretory defects, and restoration of secretory processes that promote the 'bicarbonate umbrella' may be important for amelioration of PBC-associated damage.

Single-cell RNA sequencing revealed the liver heterogeneity between egg-laying duck and ceased-laying duck

BACKGROUND

In the late phase of production, ducks untimely cease laying, leading to a lower feed conversion. Liver plays a vital role in the synthesis and transport of yolk materials during egg formation in birds. However, the molecular mechanism of liver in ceased-laying duck is far from clear, higher resolution and deeper analysis is needed. Sing-cell RNA-sequencing of 10 × Genomics platform can help to map the liver single cell gene expression atlas of Shaoxing duck and provide new insights into the liver between egg-laying and ceased-laying ducks.

RESULTS

About 20,000 single cells were profiled and 22 clusters were identified. All the clusters were identified as 6 cell types. The dominant cell type is hepatocyte, accounted for about 60% of all the cells. Of note, the heterogeneity of cells between egg-laying duck and ceased-laying duck mainly occurred in hepatocytes. Cells of cluster 3 and 12 were the unique hepatocyte states of egg-laying ducks, while cells of cluster 0 and 15 were the unique hepatocyte states of ceased-laying ducks. The expression mode of yolk precursor transporters, lipid metabolizing enzymes and fibrinogens were different in hepatocytes between egg-laying duck and ceased-laying duck. APOV1, VTG2, VTG1, APOB, RBP, VTDB and SCD might be activated in egg-laying ducks, while APOA1, APOA4, APOC3, FGB and FGG might be activated in ceased-laying ducks.

CONCLUSIONS

Our study further proofs that APOV1 and APOB play key roles in egg production, rather than APOA1 and APOA4. It is also the first to detect a correlation between the higher expression of APOC3, FGB, FGG and ceased-laying in duck.

Biology-inspired data-driven quality control for scientific discovery in single-cell transcriptomics

BACKGROUND

Quality control (QC) of cells, a critical first step in single-cell RNA sequencing data analysis, has largely relied on arbitrarily fixed data-agnostic thresholds applied to QC metrics such as gene complexity and fraction of reads mapping to mitochondrial genes. The few existing data-driven approaches perform QC at the level of samples or studies without accounting for biological variation.

RESULTS

We first demonstrate that QC metrics vary with both tissue and cell types across technologies, study conditions, and species. We then propose data-driven QC (ddqc), an unsupervised adaptive QC framework to perform flexible and data-driven QC at the level of cell types while retaining critical biological insights and improved power for downstream analysis. ddqc applies an adaptive threshold based on the median absolute deviation on four QC metrics (gene and UMI complexity, fraction of reads mapping to mitochondrial and ribosomal genes). ddqc retains over a third more cells when compared to conventional data-agnostic QC filters. Finally, we show that ddqc recovers biologically meaningful trends in gradation of gene complexity among cell types that can help answer questions of biological interest such as which cell types express the least and most number of transcripts overall, and ribosomal transcripts specifically.

CONCLUSIONS

ddqc retains cell types such as metabolically active parenchymal cells and specialized cells such as neutrophils which are often lost by conventional QC. Taken together, our work proposes a revised paradigm to quality filtering best practices-iterative QC, providing a data-driven QC framework compatible with observed biological diversity.

A monoclonal antibody targeting nonjunctional claudin-1 inhibits fibrosis in patient-derived models by modulating cell plasticity

Tissue fibrosis is a key driver of end-stage organ failure and cancer, overall accounting for up to 45% of deaths in developed countries. There is a large unmet medical need for antifibrotic therapies. Claudin-1 (CLDN1) is a member of the tight junction protein family. Although the role of CLDN1 incorporated in tight junctions is well established, the function of nonjunctional CLDN1 (njCLDN1) is largely unknown. Using highly specific monoclonal antibodies targeting a conformation-dependent epitope of exposed njCLDN1, we show in patient-derived liver three-dimensional fibrosis and human liver chimeric mouse models that CLDN1 is a mediator and target for liver fibrosis. Targeting CLDN1 reverted inflammation-induced hepatocyte profibrogenic signaling and cell fate and suppressed the myofibroblast differentiation of hepatic stellate cells. Safety studies of a fully humanized antibody in nonhuman primates did not reveal any serious adverse events even at high steady-state concentrations. Our results provide preclinical proof of concept for CLDN1-specific monoclonal antibodies for the treatment of advanced liver fibrosis and cancer prevention. Antifibrotic effects in lung and kidney fibrosis models further indicate a role of CLDN1 as a therapeutic target for tissue fibrosis across organs. In conclusion, our data pave the way for further therapeutic exploration of CLDN1-targeting therapies for fibrotic diseases in patients.

Obeticholic acid and 5β-cholanic acid 3 exhibit anti-tumor effects on liver cancer through CXCL16/CXCR6 pathway

Hepatocellular carcinoma (HCC) is the most common type of liver malignancy with a high incidence and mortality rate. Previous in vitro and in vivo studies have confirmed that liver sinusoidal endothelial cells (LSEC) secrete CXCL16, which acts as a messenger to increase the hepatic accumulation of CXCR6⁺ natural killer T (NKT) cells and exert potent antitumor effects. However, evidence for this process in humans is lacking and its clinical significance is still unclear. In this study, by dissecting the human HCC single-cell RNA-seq data, we verified this process through cellphoneDB. NKT cells in patients with high expression of CXCL16 exhibited a higher activation state and produced more interferon-γ (IFN-γ) compared with those with low expression. We next investigated the signaling pathways between activated (CD69 high) and unactivated NKT cells (CD69 low) using NKT cell-developmental trajectories and functional enrichment analyses. In vivo experiments, we found that farnesoid X receptor agonist (obeticholic acid) combined with the takeda G protein coupled receptor 5 antagonist (5β-cholanic acid 3) exhibited significant tumor suppressive effects in the orthotopic liver tumor model and this result may be related to the CXCL16/CXCR6 axis. In conclusion, our study provides the basis and potential strategies for HCC immunotherapy based on NKT cells.

Improving the performance of single-cell RNA-seq data mining based on relative expression orderings

The advent of single-cell RNA-sequencing (scRNA-seq) provides an unprecedented opportunity to explore gene expression profiles at the single-cell level. However, gene expression values vary over time and under different conditions even within the same cell. There is an urgent need for more stable and reliable feature variables at the single-cell level to depict cell heterogeneity. Thus, we construct a new feature matrix called the delta rank matrix (DRM) from scRNA-seq data by integrating an a priori gene interaction network, which transforms the unreliable gene expression value into a stable gene interaction/edge value on a single-cell basis. This is the first time that a gene-level feature has been transformed into an interaction/edge-level for scRNA-seq data analysis based on relative expression orderings. Experiments on various scRNA-seq datasets have demonstrated that DRM performs better than the original gene expression matrix in cell clustering, cell identification and pseudo-trajectory reconstruction. More importantly, the DRM really achieves the fusion of gene expressions and gene interactions and provides a method of measuring gene interactions at the single-cell level. Thus, the DRM can be used to find changes in gene interactions among different cell types, which may open up a new way to analyze scRNA-seq data from an interaction perspective. In addition, DRM provides a new method to construct a cell-specific network for each single cell instead of a group of cells as in traditional network construction methods. DRM's exceptional performance is due to its extraction of rich gene-association information on biological systems and stable characterization of cells.

Enzymes of Fibrosis in Chronic Liver Disease

Introduction: Liver fibrosis has been extensively studied at the cellular and molecular level, but very few data exist on the final enzymatic stages of collagen synthesis (prolyl hydroxylase, PH) and degradation (matrix metalloproteinases, MMPs), particularly in primary biliary cholangitis (PBC). Aim: We studied enzyme activities in liver tissue from patients with chronic liver diseases and compared them to normal livers. Patients: Eighteen patients with PBC of early and late stages (Ludwig’s classification) and seven on treatment with ursodeoxycholate (UDCA) were studied and compared to 34 patients with alcoholic liver disease (ALD), 25 patients with chronic viral liver disease and five normal biopsies. Sera were available from a total of 140 patients. Methods: The tritiated water released from the tritiated proline was measured in PH assessment. 14C intact and heat-denatured collagen substrates were used to measure collagenase and gelatinases, respectively. 3H Elastin was the substrate for elastase. In serum, ELISAs were used for MMP-1, TIMP-1, and TIMP-2 measurements while MMP-2 and MMP-9 were estimated by zymography. Results: PH was significantly increased in early and late PBC. Collagenase was reduced only in the late stages (p < 0.01), where the ratio PH/collagenase was increased. UDCA treatment restored values to almost normal. Gelatinases were reduced in late stages (p < 0.05). In contrast to PBC and ALD fibrosis, collagen synthesis is not increased in viral fibrosis. The balance shifted towards collagen deposition due to reduced degradation. Interestingly, gelatinolytic activity is not impaired in ALD. Elastase was similar to controls in all diseases studied. TIMP-1 was reduced in early PBC and viral and alcoholic hepatitis and cirrhosis (p < 0.001). Conclusions: (1) There is evidence that collagen synthesis increases in the early stages of PBC, but the collagenolytic mechanism may compensate for the increased synthesis. (2) In viral disease, fibrosis may be due to decreased degradation rather than increased synthesis. (3) The final biochemical stages of liver fibrosis may be quantitatively different according to underlying etiology.

Comparison of cell type distribution between single-cell and single-nucleus RNA sequencing: enrichment of adherent cell types in single-nucleus RNA sequencing

Single-cell ribonucleic acid (RNA) sequencing (scRNA-seq) is an effective technique for estimating the cellular composition and transcriptional profiles of individual cells from fresh tissue. Single-nucleus RNA sequencing (snRNA-seq) is necessary to perform this type of analysis in frozen or difficult-to-dissociate tissues, which cannot be subjected to scRNA-seq. This difference in the state of tissues leads to variation in cell-type distributions among each platform. To identify the characteristics of these methods and their differences, scRNA-seq and snRNA-seq were performed in parallel for colon and liver tissues. The two platforms revealed similar diversity but different proportions of cell types in matched tissues. The proportions of epithelial cells in the colon and hepatocytes in the liver were relatively high in snRNA-seq and that of immune cells was relatively high in scRNA-seq. This difference could be explained by variations in the expression scores of adhesion genes due to the disruption of the cytoplasmic contents during scRNA-seq. The enrichment of epithelial cells in the colon resulted in a discrepancy in the differentiation of epithelial cells. This enrichment was also well matched with the images of hematoxylin and eosin staining and the estimated distribution of cell types in bulk RNA sequencing. These results showed that snRNA-seq could be used to analyze tissues that cannot be subjected to scRNA-seq and provides more information in specific cell type analysis.

Single-cell transcriptomes from turtle livers reveal sensitivity of hepatic immune cells to bacteria-infection

The liver is important in the synthesis, metabolism and storage of nutrients, detoxification and immune response of the body, and the liver immune response against exogenous pathogens from the intestinal tract plays a key role in the immune activities. However, the cellular composition of the liver immune atlas remains sparsely studied in reptiles. We used single-cell RNA sequencing to identify the cellular profile of the liver of the Chinese soft-shelled turtle (Pelodiscus sinensis). We obtained the transcriptional landscape based on 9938 cells from the fractionation of fresh hepatic tissues from two individuals, uninfected and infected with bacteria (Aeromonas hydrophila). We identified seven hepatic immune cell subsets, including plasma, erythroid, T/NK, B, endothelial, dendritic and Kupffer cells. Bacteria-infection altered the number of liver immune cells, as revealed by the fact that the infected turtle had more plasma, endothelial and Kupffer cells and fewer T/NK, dendritic and erythroid cells than did the uninfected turtle. Our study is the first to provide a comprehensive view of the hepatic immune landscape of P. sinensis at the single-cell resolution that outlines the characteristics of immune cells in the turtle liver and provides a liver transcriptome baseline for turtle immunology.

Modelling metabolic diseases and drug response using stem cells and organoids

Metabolic diseases, including obesity, diabetes mellitus and cardiovascular disease, are a major threat to health in the modern world, but efforts to understand the underlying mechanisms and develop rational treatments are limited by the lack of appropriate human model systems. Notably, advances in stem cell and organoid technology allow the generation of cellular models that replicate the histological, molecular and physiological properties of human organs. Combined with marked improvements in gene editing tools, human stem cells and organoids provide unprecedented systems for studying mechanisms of metabolic diseases. Here, we review progress made over the past decade in the generation and use of stem cell-derived metabolic cell types and organoids in metabolic disease research, especially obesity and liver diseases. In particular, we discuss the limitations of animal models and the advantages of stem cells and organoids, including their application to metabolic diseases. We also discuss mechanisms of drug action, understanding the efficacy and toxicity of existing therapies, screening for new treatments and pursuing personalized therapies. We highlight the potential of combining stem cell-derived organoids with gene editing and functional genomics to revolutionize the approach to finding treatments for metabolic diseases.

Reflected stemness as a potential driver of the tumour microenvironment

A fundamental requirement for cancer initiation is the activation of developmental programmes by mutant cells. Oncogenic signals often confer an undifferentiated, stem cell-like phenotype that supports the long-term proliferative potential of cancer cells. Although cancer is a genetically driven disease, mutations in cancer-driver genes alone are insufficient for tumour formation, and the proliferation of cells harbouring oncogenic mutations depends on their microenvironment. In this Opinion article we discuss how the reprogrammed status of cancer cells not only represents the essence of their tumorigenicity but triggers 'reflected stemness' in their surrounding normal counterparts. We propose that this reciprocal interaction underpins the establishment of the tumour microenvironment (TME).

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

Background

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

Methods

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

Results

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

Conclusions

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

Liver tumour immune microenvironment subtypes and neutrophil heterogeneity

The heterogeneity of the tumour immune microenvironment (TIME), organized by various immune and stromal cells, is a major contributing factor of tumour metastasis, relapse and drug resistance^1-3, but how different TIME subtypes are connected to the clinical relevance in liver cancer remains unclear. Here we performed single-cell RNA-sequencing (scRNA-seq) analysis of 189 samples collected from 124 patients and 8 mice with liver cancer. With more than 1 million cells analysed, we stratified patients into five TIME subtypes, including immune activation, immune suppression mediated by myeloid or stromal cells, immune exclusion and immune residence phenotypes. Different TIME subtypes were spatially organized and associated with chemokine networks and genomic features. Notably, tumour-associated neutrophil (TAN) populations enriched in the myeloid-cell-enriched subtype were associated with an unfavourable prognosis. Through in vitro induction of TANs and ex vivo analyses of patient TANs, we showed that CCL4⁺ TANs can recruit macrophages and that PD-L1⁺ TANs can suppress T cell cytotoxicity. Furthermore, scRNA-seq analysis of mouse neutrophil subsets revealed that they are largely conserved with those of humans. In vivo neutrophil depletion in mouse models attenuated tumour progression, confirming the pro-tumour phenotypes of TANs. With this detailed cellular heterogeneity landscape of liver cancer, our study illustrates diverse TIME subtypes, highlights immunosuppressive functions of TANs and sheds light on potential immunotherapies targeting TANs.

Immunobiology of the biliary tract system

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

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

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

Challenges and considerations for single-cell and spatially resolved transcriptomics sample collection during spaceflight

Single-cell RNA sequencing (scRNA-seq) and spatially resolved transcriptomics (SRT) have experienced rapid development in recent years. The findings of spaceflight-based scRNA-seq and SRT investigations are likely to improve our understanding of life in space and our comprehension of gene expression in various cell systems and tissue dynamics. However, compared to their Earth-based counterparts, gene expression experiments conducted in spaceflight have not experienced the same pace of development. Out of the hundreds of spaceflight gene expression datasets available, only a few used scRNA-seq and SRT. In this perspective piece, we explore the growing importance of scRNA-seq and SRT in space biology and discuss the challenges and considerations relevant to robust experimental design to enable growth of these methods in the field.

In vivo partial reprogramming by bacteria promotes adult liver organ growth without fibrosis and tumorigenesis

Ideal therapies for regenerative medicine or healthy aging require healthy organ growth and rejuvenation, but no organ-level approach is currently available. Using Mycobacterium leprae (ML) with natural partial cellular reprogramming capacity and its animal host nine-banded armadillos, we present an evolutionarily refined model of adult liver growth and regeneration. In infected armadillos, ML reprogram the entire liver and significantly increase total liver/body weight ratio by increasing healthy liver lobules, including hepatocyte proliferation and proportionate expansion of vasculature, and biliary systems. ML-infected livers are microarchitecturally and functionally normal without damage, fibrosis, or tumorigenesis. Bacteria-induced reprogramming reactivates liver progenitor/developmental/fetal genes and upregulates growth-, metabolism-, and anti-aging-associated markers with minimal change in senescence and tumorigenic genes, suggesting bacterial hijacking of homeostatic, regeneration pathways to promote de novo organogenesis. This may facilitate the unraveling of endogenous pathways that effectively and safely re-engage liver organ growth, with broad therapeutic implications including organ regeneration and rejuvenation.

Engineered matrix microenvironments reveal the heterogeneity of liver sinusoidal endothelial cell phenotypic responses

Fibrosis is one of the hallmarks of chronic liver disease and is associated with aberrant wound healing. Changes in the composition of the liver microenvironment during fibrosis result in a complex crosstalk of extracellular cues that promote altered behaviors in the cell types that comprise the liver sinusoid, particularly liver sinusoidal endothelial cells (LSECs). Recently, it has been observed that LSECs may sustain injury before other fibrogenesis-associated cells of the sinusoid, implicating LSECs as key actors in the fibrotic cascade. A high-throughput cellular microarray platform was used to deconstruct the collective influences of defined combinations of extracellular matrix (ECM) proteins, substrate stiffness, and soluble factors on primary human LSEC phenotype in vitro. We observed remarkable heterogeneity in LSEC phenotype as a function of stiffness, ECM, and soluble factor context. LYVE-1 and CD-31 expressions were highest on 1 kPa substrates, and the VE-cadherin junction localization was highest on 25 kPa substrates. Also, LSECs formed distinct spatial patterns of LYVE-1 expression, with LYVE-1+ cells observed in the center of multicellular domains, and pattern size regulated by microenvironmental context. ECM composition also influenced a substantial dynamic range of expression levels for all markers, and the collagen type IV was observed to promote elevated expressions of LYVE-1, VE-cadherin, and CD-31. These studies highlight key microenvironmental regulators of LSEC phenotype and reveal unique spatial patterning of the sinusoidal marker LYVE-1. Furthermore, these data provide insight into understanding more precisely how LSECs respond to fibrotic microenvironments, which will aid drug development and identification of targets to treat liver fibrosis.

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

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

Identification of a rare Gli1+ progenitor cell population contributing to liver regeneration during chronic injury

In adults, hepatocytes are mainly replenished from the existing progenitor pools of hepatocytes and cholangiocytes during chronic liver injury. However, it is unclear whether other cell types in addition to classical hepatocytes and cholangiocytes contribute to hepatocyte regeneration after chronic liver injuries. Here, we identified a new biphenotypic cell population that contributes to hepatocyte regeneration during chronic liver injuries. We found that a cell population expressed Gli1 and EpCAM (EpCAM⁺Gli1⁺), which was further characterized with both epithelial and mesenchymal identities by single-cell RNA sequencing. Genetic lineage tracing using dual recombinases revealed that Gli1⁺ nonhepatocyte cell population could generate hepatocytes after chronic liver injury. EpCAM⁺Gli1⁺ cells exhibited a greater capacity for organoid formation with functional hepatocytes in vitro and liver regeneration upon transplantation in vivo. Collectively, these findings demonstrate that EpCAM⁺Gli1⁺ cells can serve as a new source of liver progenitor cells and contribute to liver repair and regeneration.

De novo analysis of bulk RNA-seq data at spatially resolved single-cell resolution

Uncovering the tissue molecular architecture at single-cell resolution could help better understand organisms' biological and pathological processes. However, bulk RNA-seq can only measure gene expression in cell mixtures, without revealing the transcriptional heterogeneity and spatial patterns of single cells. Herein, we introduce Bulk2Space ( https://github.com/ZJUFanLab/bulk2space ), a deep learning framework-based spatial deconvolution algorithm that can simultaneously disclose the spatial and cellular heterogeneity of bulk RNA-seq data using existing single-cell and spatial transcriptomics references. The use of bulk transcriptomics to validate Bulk2Space unveils, in particular, the spatial variance of immune cells in different tumor regions, the molecular and spatial heterogeneity of tissues during inflammation-induced tumorigenesis, and spatial patterns of novel genes in different cell types. Moreover, Bulk2Space is utilized to perform spatial deconvolution analysis on bulk transcriptome data from two different mouse brain regions derived from our in-house developed sequencing approach termed Spatial-seq. We have not only reconstructed the hierarchical structure of the mouse isocortex but also further annotated cell types that were not identified by original methods in the mouse hypothalamus.

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

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

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

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

Long-term and non-invasive in vivo tracking of DiD dye-labeled human hepatic progenitors in chronic liver disease models

BACKGROUND

Chronic liver diseases (CLD) are the major public health burden due to the continuous increasing rate of global morbidity and mortality. The inherent limitations of organ transplantation have led to the development of stem cell-based therapy as a supportive and promising therapeutic option. However, identifying the fate of transplanted cells in vivo represents a crucial obstacle.

AIM

To evaluate the potential applicability of DiD dye as a cell labeling agent for long-term, and non-invasive in vivo tracking of transplanted cells in the liver.

METHODS

Magnetically sorted, epithelial cell adhesion molecule positive (1 × 10⁶ cells/mL) fetal hepatic progenitor cells were labeled with DiD dye and transplanted into the livers of CLD-severe combined immunodeficiency (SCID) mice. Near-infrared (NIR) imaging was performed for in vivo tracking of the DiD-labeled transplanted cells along with colocalization of hepatic markers for up to 80 d. The existence of human cells within mouse livers was identified using Alu polymerase chain reaction and sequencing.

RESULTS

NIR fluorescence imaging of CLD-SCID mice showed a positive fluorescence signal of DiD at days 7, 15, 30, 45, 60, and 80 post-transplantation. Furthermore, positive staining of cytokeratin, c-Met, and albumin colocalizing with DiD fluorescence clearly demonstrated that the fluorescent signal of hepatic markers emerged from the DiD-labeled transplanted cells. Recovery of liver function was also observed with serum levels of glutamic-oxaloacetic transaminase, glutamate-pyruvate transaminase, and bilirubin. The detection of human-specific Alu sequence from the transplanted mouse livers provided evidence for the survival of transplanted cells at day 80.

CONCLUSION

DiD-labeling is promising for long-term and non-invasive in vivo cell tracking, and understanding the regenerative mechanisms incurred by the transplanted cells.

Multiplexed Digital Spatial Protein Profiling Reveals Distinct Phenotypes of Mononuclear Phagocytes in Livers with Advanced Fibrosis

Background and Aims: Intrahepatic mononuclear phagocytes (MPs) are critical for the initiation and progression of liver fibrosis. In this study, using multiplexed digital spatial protein profiling, we aimed to derive a unique protein signature predicting advanced liver fibrosis. Methods: Snap-frozen liver tissues from various chronic liver diseases were subjected to spatially defined protein-based multiplexed profiling (Nanostring GeoMX^TM). A single-cell RNA sequencing analysis was performed using Gene Expression Omnibus (GEO) datasets from normal and cirrhotic livers. Results: Sixty-four portal regions of interest (ROIs) were selected for the spatial profiling. Using the results from the CD68⁺ area, a highly sensitive and specific immune-related protein signature (CD68, HLA-DR, OX40L, phospho-c-RAF, STING, and TIM3) was developed to predict advanced (F3 and F4) fibrosis. A combined analysis of single-cell RNA sequencing data from GEO datasets (GSE136103) and spatially-defined, protein-based multiplexed profiling revealed that most proteins upregulated in F0-F2 livers in portal CD68⁺ cells were specifically marked in tissue monocytes, whereas proteins upregulated in F3 and F4 livers were marked in scar-associated macrophages (SAMacs) and tissue monocytes. Internal validation using mRNA expression data with the same cohort tissues demonstrated that mRNA levels for TREM2, CD9, and CD68 are significantly higher in livers with advanced fibrosis. Conclusions: In patients with advanced liver fibrosis, portal MPs comprise of heterogeneous populations composed of SAMacs, Kupffer cells, and tissue monocytes. This is the first study that used spatially defined protein-based multiplexed profiling, and we have demonstrated the critical difference in the phenotypes of portal MPs between livers with early- or late-stage fibrosis.

The Intratumoral Bacterial Metataxonomic Signature of Hepatocellular Carcinoma

Microbiota is implicated in hepatocellular carcinoma (HCC). The spectrum of intratumoral microbiota associated with HCC progression remains elusive. Fluorescence in situ hybridization revealed that microbial DNAs were distributed in the cytosol of liver hepatocytes and erythrocytes. Viable anaerobic or aerobic bacteria were recovered in HCC tissues by fresh tissue culture. We performed a comprehensive DNA sequencing of bacterial 16S rRNA genes in 156 samples from 28 normal liver, 64 peritumoral, and 64 HCC tissues, and the DNA sequencing yielded 4.2 million high-quality reads. Both alpha and beta diversity in peritumor and HCC microbiota were increased compared to normal controls. The most predominant phyla in HCC were Patescibacteria, Proteobacteria, Bacteroidota, Firmicutes, and Actinobacteriota. phyla of Proteobacteria, Firmicutes, and Actinobacteriota, and classes of Bacilli and Actinobacteria, were consistently enriched in peritumor and HCC tissues, while Gammaproteobacteria was especially abundant in HCC tissues compared to normal controls. Streptococcaceae and Lactococcus were the marker taxa of HCC cirrhosis. The Staphylococcus branch and Caulobacter branch were selectively enriched in HBV-negative HCCs. The abundance of Proteobacteria, Gammaproteobacteria, Firmicutes, Actinobacteriota, and Saccharimonadia were associated with the clinicopathological features of HCC patients. The inferred functions of different taxa were changed between the microbiota of normal liver and peritumor/HCC. Random forest machine learning achieved great discriminative performance in HCC prediction (area under the curve [AUC] = 1.00 in the training cohort, AUC = 0.950 for top five class signature, and AUC = 0.943 for the top 50 operational taxonomy units [OTUs] in the validation cohort). Our analysis highlights the complexity and diversity of the liver and HCC microbiota and established a specific intratumoral microbial signature for the potential prediction of HCC. IMPORTANCE Gut microbiome is an important regulator of hepatic inflammation, detoxification, and immunity, and contributes to the carcinogenesis of liver cancer. Intratumoral bacteria are supposed to be closer to the tumor cells, forming a microenvironment that may be relevant to the pathological process of hepatocellular carcinoma (HCC). However, the presence of viable intratumoral bacteria remains unclear. It is worth exploring whether the metataxonomic characteristics of intratumoral bacteria can be used as a potential marker for HCC prediction. Here, we present the first evidence of the existence of viable intratumoral bacteria in HCC using the tissue culture method. We revealed that microbial DNAs were distributed in the cytosol of liver hepatocytes and erythrocytes. We analyzed the diversity, structure, and abundance of normal liver and HCC microbiota. We built a machine learning model for HCC prediction using intratumoral bacterial features. We show that specific taxa represent potential targets for both therapeutic and diagnostic interventions.

Spatial transcriptome profiling of normal human liver

The comprehensive study of the spatial-cellular anatomy of the human liver is critical to addressing the cellular origins of liver disease. Here we conducted spatial transcriptomics on normal human liver tissue sections, providing detailed information of liver zonation at the transcriptional level. We present 6581 high-quality spots from normal livers of two human donors. In this dataset, cells were mainly hepatocytes, and we classified them into four sub-groups. Collectively, these data provide a reliable reference for studies on spatial heterogeneity of liver lobules.

Measurement(s)	mRNA Expression
Technology Type(s)	spatial transcriptomics
Sample Characteristic - Organism	Homo sapiens

Discovery of a MUC3B gene reconstructs the membrane mucin gene cluster on human chromosome 7

Human tissue surfaces are coated with mucins, a family of macromolecular sugar-laden proteins serving diverse functions from lubrication to the formation of selective biochemical barriers against harmful microorganisms and molecules. Membrane mucins are a distinct group of mucins that are attached to epithelial cell surfaces where they create a dense glycocalyx facing the extracellular environment. All mucin proteins carry long stretches of tandemly repeated sequences that undergo extensive O-linked glycosylation to form linear mucin domains. However, the repetitive nature of mucin domains makes them prone to recombination and renders their genetic sequences particularly difficult to read with standard sequencing technologies. As a result, human mucin genes suffer from significant sequence gaps that have hampered the investigation of gene function in health and disease. Here we leveraged a recent human genome assembly to characterize a previously unmapped MUC3B gene located at the q22 locus on chromosome 7, within a cluster of four structurally related membrane mucin genes that we name the MUC3 cluster. We found that MUC3B shares high sequence identity with the known MUC3A gene and that the two genes are governed by evolutionarily conserved regulatory elements. Furthermore, we show that MUC3A, MUC3B, MUC12, and MUC17 in the human MUC3 cluster are expressed in intestinal epithelial cells (IECs). Our results complete existing genetic gaps in the MUC3 cluster which is a conserved genetic unit in vertebrates. We anticipate our results to be the starting point for the detection of disease-associated polymorphisms in the human MUC3 cluster. Moreover, our study provides the basis for the exploration of intestinal mucin gene function in widely used experimental models such as human intestinal organoids and genetic mouse models.

Online single-cell data integration through projecting heterogeneous datasets into a common cell-embedding space

Computational tools for integrative analyses of diverse single-cell experiments are facing formidable new challenges including dramatic increases in data scale, sample heterogeneity, and the need to informatively cross-reference new data with foundational datasets. Here, we present SCALEX, a deep-learning method that integrates single-cell data by projecting cells into a batch-invariant, common cell-embedding space in a truly online manner (i.e., without retraining the model). SCALEX substantially outperforms online iNMF and other state-of-the-art non-online integration methods on benchmark single-cell datasets of diverse modalities, (e.g., single-cell RNA sequencing, scRNA-seq, single-cell assay for transposase-accessible chromatin use sequencing, scATAC-seq), especially for datasets with partial overlaps, accurately aligning similar cell populations while retaining true biological differences. We showcase SCALEX's advantages by constructing continuously expandable single-cell atlases for human, mouse, and COVID-19 patients, each assembled from diverse data sources and growing with every new data. The online data integration capacity and superior performance makes SCALEX particularly appropriate for large-scale single-cell applications to build upon previous scientific insights.

A comprehensive transcriptomic comparison of hepatocyte model systems improves selection of models for experimental use

The myriad of available hepatocyte in vitro models provides researchers the possibility to select hepatocyte-like cells (HLCs) for specific research goals. However, direct comparison of hepatocyte models is currently challenging. We systematically searched the literature and compared different HLCs, but reported functions were limited to a small subset of hepatic functions. To enable a more comprehensive comparison, we developed an algorithm to compare transcriptomic data across studies that tested HLCs derived from hepatocytes, biliary cells, fibroblasts, and pluripotent stem cells, alongside primary human hepatocytes (PHHs). This revealed that no HLC covered the complete hepatic transcriptome, highlighting the importance of HLC selection. HLCs derived from hepatocytes had the highest transcriptional resemblance to PHHs regardless of the protocol, whereas the quality of fibroblasts and PSC derived HLCs varied depending on the protocol used. Finally, we developed and validated a web application (HLCompR) enabling comparison for specific pathways and addition of new HLCs. In conclusion, our comprehensive transcriptomic comparison of HLCs allows selection of HLCs for specific research questions and can guide improvements in culturing conditions.