Introduction
Liver diseases have become a global health burden in recent decades. For example, cirrhosis causes 1.16 million deaths, and liver cancer causes 0.788 million deaths annually, making them the 11th and 16th most common causes of death.1 More than one-fifth of the Chinese population is affected by liver diseases, particularly hepatitis B viral infections, while cases of nonalcoholic hepatitis and autoimmune hepatitis are increasing.2 Persistent chronic inflammation in the liver is closely associated with cirrhosis and ultimately leads to liver cancer or liver failure, and the main cause of acute liver failure in developing countries is viral hepatitis.3,4 However, effective treatment options for liver diseases are so limited that liver transplantation is regarded as the only viable solution, but its demand fairly exceeds the supply.1 Therefore, better tools are urgently required to elucidate the pathological mechanism and development of liver diseases with a clearer perspective. Their application may contribute to precise, personalized therapy and a decrease in mortality.
Each cell is unique for it occupies an exclusive position in space, carries distinct errors in its copied genome and is subject to programmed and induced changes in gene expression. Yet traditional sequencing is performed on tissue samples or cell populations, in which biological differences between cells can be obscured by averaging or mistaken for technical noise. In contrast, single-cell sequencing technology can isolate several thousand single cells at one time and has the advantage of detecting heterogeneity among individual cells then depicting cellular profiles, which can help us identify the rare populations of cells that have derived significance for disease progression with unprecedented high resolution.5 The workflow of single-cell sequencing experiments and analyses generally includes single-cell suspension preparation, single-cell sorting, amplification and library preparation, high-throughput sequencing, and bioinformatics analysis (Fig. 1). Tang et al.6 proposed the first single-cell transcriptome sequencing (scRNA-seq) technology in 2009. Since then, the single-cell field has advanced rapidly, and several new approaches have been developed, including genome, transcriptome, proteome, spatial transcriptome, epigenome. Single-cell sequencing can not only analyze the heterogeneity of cells with the same phenotype but also obtain the genetic information of hard-to-culture microorganisms and valuable clinical samples.
In this article, we review the current single-cell approaches and discussed their key principles along with their advantages and disadvantages to understand which method is more suitable under what kind of circumstances (Table 1).7–23 We also discuss the novel application of those techniques in studies on liver physiology and pathology and the application prospects and technical improvements in single-cell technology to find effective individualized treatment strategies for patients with liver diseases and reduce the overall morbidity and mortality related to liver diseases.
Table 1Summary of single-cell analysis methods
Single-cell histology | Application | Single-cell sequencing | Advantages | Disadvantages | Ref. |
---|
Transcriptome | Cell type discovery and characterization in diseases; cell development trajectory and plasticity | Smart-seq | Full-length coverage across transcripts | Unstable; transcript bias | 20 |
| | Smart-seq2 | Captures rare genes; commercial off-the-shelf SMARTer kit | Can not analyze nonpolyA RNAs; nonstrand specificity | 8 |
| | Drop-seq | High throughput; low cost | Low sensitivity | 21 |
| | 10× Genomics Chromium | High throughput; high sensitivity; short cycle time; low technical noise | High cost | 22 |
| | snRNA-seq | Allow for frozen samples; relatively simple to operate | Few genes detected; cellular preference | 7 |
Spatial transcriptome | Spatial distribution of heterogeneous cells and spatial patterns of gene expression in different tissues; tumor microenvironment | 10× Visium | Commercial platform is mature and stable; widely used; with supporting verification platform | Low spatial resolution; non-single cell level | 9 |
| | GeoMx DSP | Compatible with FFPE samples and fresh frozen samples; for in-depth analysis of local tissues | Low throughput | 10 |
| | DBiT-seq | High spatial resolution; co-localization of mRNA and protein; simplified workflow | Inability to readily identify cell types | 11 |
Proteome | Disease heterogeneity and treatment resistance; cancer immunotherapy; high throughput drug screening | scWestern | High specificity; | Low protein flux | 23 |
| | SCoPE-MS | Accurately differentiate the proteome of different cells | Low protein flux Complex workflow | 12 |
| | nanoPOTS | Microprotein identification; high quantitative accuracy and protein coverage | | 13 |
Epigenome | Genetic characteristics of disease; deeper understanding of gene regulation from open chromatin; infer development trajectories and identify regulatory elements | scCGI-seq | High accuracy; simplified workflow | Low coverage | 14 |
| | CUT & RUN technique | Low sample requirement; high resolution and accuracy; commercial kits | Highly antibody specificity-dependent | 17 |
| | scATAC-seq | High resolution, sensitivity, and accuracy | Low recovery of DNA fragments | 15 |
| | scChIP-seq | High throughput, high coverage | Highly antibody-dependent | 16 |
Multiomics | Integrated analysis different histological techniques to characterize the state of gene regulation more comprehensively within cells and reveal the mechanisms of disease regulation | scTrio-seq | Simultaneous analysis of genome,transcriptome, and epigenome | Low throughput | 18 |
| | CITE-seq | Simultaneous analysis of transcriptome and protein expression; high throughput | Only surface proteins can be measured | 19 |
An overview of single-cell technologies
scRNA-seq
scRNA-seq measures gene expression at the single-cell level and reveals heterogeneity of gene expression in individual cells or homologous cell types, its core principle refers to tag-based identification of single cells. A unique cell tag attached to a barcode is added to each cell when reverse transcription is performed prior to RNA sequencing. Therefore, even if all cells are mixed and sequenced, RNA fragments carrying the same barcode can be considered as coming from the same cell. By building a library once, we can measure the information of tens of thousands of single cells (Fig. 1). A challenge in conducting research using scRNA-seq is that high quality suspensions based on fresh tissue are not available for all tissues. In contrast, single-nucleus RNA sequencing (snRNA-seq)7 is not limited by dissociation conditions and can be used for analyzing cell and gene expression of primary tissues obtained from frozen samples. Thus, this process complements scRNA-seq, but the nuclear sequencing depth is low. Some studies have compared the differences in sample type, study purpose, and target cell population between scRNA-seq and snRNA-seq for various tissues and diseases.24–26
In addition, immune cell receptor sequencing, which focuses on immune T and B cells, can comprehensively resolve T-cell receptor/B-cell receptor (TCR/BCR) gene rearrangement base sequences and the abundance of each sequence to reveal the mechanism of T/B cell-mediated humoral immune response. Single-cell analysis was first used for examining immune cells around 2013 and early 2014; specifically, researchers used scRNA-seq to generate immune cell profiles from tissues.27,28 Recent advancements in scRNA-seq have enabled cellular transcriptome data and immune repertoire information to be combined for the coverage of full-length mRNA,29 which can be used to determine the functional changes and interactions of subsets of immune cells during pathogen invasion and the infiltration landscape of T cells in the tumor microenvironment.30,31 The most widely used transcriptome sequencing technologies in the market today are droplet-based 10× Genomics32 and Smart-seq28 a single-cell full-length transcriptome sequencing technology (Table 1).
Spatial transcriptome technology (ST)
Compared with bulk RNA-seq and scRNA-seq, ST can be used to analyze the information of tissue sections in space without cell enzymatic hydrolysis and maximum restoration of the spatial specificity of gene expression in individual cells in human organs.33 It has advanced considerably in the past 30 years, initially, transcriptome technology could be used to analyze only a small number of genes, but now it can be used for analyzing the complete transcriptome. From the early technologies, such as single-molecule fluorescence in situ hybridization (smFISH)34 and laser capture microdissection (LCM), to the current technologies, such as 10× Visium,9 GeoMx DSP,10 and DBiT-seq (Table 1).11 More details on ST can be found in a recent review.35 Among all ST, for now, the 10× Visium ST is the most widely used.
Single-cell proteomic analysis
Qualitative and quantitative analysis of the single-cell proteome can avoid processing information on the expression level changes caused by post transcriptional and post translational modifications and efficaciously obtain information on the regulatory nature of cells. No tools are available for analyzing the whole proteome at the single-cell level, but various single-cell proteome analysis methods are available in different subsegments, including methods based on antibody-based strategies, mass spectrometry (SCoPE-MS),12 Proseek multiplex technology based on polymerase chain reaction (commonly known as PCR) immunoassay,36 and others. Recently, nanodroplet processing in one pot for trace samples (nanoPOTS) microsample processing13 combined with LCM and mass spectrometry found that in mouse uterine tissue sections prepared for blastocyst implantation, quantitative cell type-specific images of >2,000 proteins could be generated with a high spatial resolution.37 Protein profiling characterized by captured ion mobility has led to the development of four-dimensional (4D) proteomics, which can perform high-depth quantitative analysis 4D proteomics quantification (4D-LFQ, 4D-DIA), and 4D protein modification omics.38,39
Single-cell epigenomic analysis
Cellular epigenetic heterogeneity at the single-cell level can be used to explain the pathological regulatory mechanisms of diseases in the absence of altered genetic information. The detection of DNA methylation in single cells was first conducted in 2015 when Farlik et al.40 developed genome-wide bisulfite sequencing, which was used to analyze genome-wide methylation sites at the single-cell level, but required many DNA samples. Currently, methylation sequencing targeting genomic regions with high CpG content (scCGI-seq)14 and DNA methylation analysis using indexing methods (sci-MET)41 at the single-cell level are also available. There are significant differences in chromatin accessibility between cells, which often indicate gene regulation functions. For analyzing open chromatin conformations with loose features at the single-cell level, sensitive enzymes, such as the Tn5 transposase-based single-cell assay for transposase-accessible chromatin using sequencing (also referred to as scATAC-seq)15 and the DNase I-based single-cell DNase sequencing (also referred to as scDNase-seq)42 are required. The latter can detect a much wider range of chromatin, up to 350,000 reads/cell. To determine the link between the diversity of histone modification sites and changes in the chromatin state, single-cell histone modification, single-cell chromatin immune-precipitation sequencing (also referred to as scChIP-seq), was used in 2015.16 Subsequently, the cleavage under targets and release using nuclease (CUT & RUN) technique, indexing single-cell immunocleavage (also referred to as iscChICseq), and antibody-guided chromatin tagmentation (also referred to as ACT-seq)17,43,44 were developed to overcome the low throughput characteristics of cells. More information on single-cell epigenetic techniques is summarized in a previous review.45
Single-cell multiomics analysis
In the near future, with the help of single-cell multiomics analysis, the genomic information of cells and the regulation of transcription and translation in a cell can be presented in a comprehensive manner. In 2016, Hou et al.18 developed single-cell triple omics sequencing (also referred to as scTrio-seq), and were the first to simultaneously achieve high-throughput sequencing of genomic copy number variation, DNA methylome, and transcriptome.18 Cellular indexing of transcriptomes and epitopes by sequencing (also referred to as CITE-seq) using oligonucleotide-tagged antibodies can be applied to target cell surface proteins and analyze differences between protein and transcript levels in individual cells.19 Presently, multimodal analysis has greatly improved our ability to determine cellular states. The single-cell approaches described above has been applied in the field of hepatology.
Application of findings from single-cell approaches in liver physiology and pathology
The liver is the largest digestive gland and detoxification organ in the human body and participates in the synthesis and decomposition of various substances.46 The incidence of hepatitis B virus (HBV) related to end-stage liver disease have become the main causes of mortality due to liver diseases in China,47 and early in vitro and in vivo studies have found that hepatitis B virus X protein (HBx) can induce malignant transformation and carcinogenesis in hepatocytes.48 In 2019, integrated proteomic genomics identified the key signaling and metabolic pathways in HBV-associated hepatocellular carcinoma (HCC), where mutant CTNNB1-associated ALDOA phosphorylation was identified to be involved in HCC metabolic reprogramming, and metabolic alterations may be the most important factor for advanced cancer and poor clinical outcomes.49 Studies using single-cell technology have provided valuable information on the characteristics of single cells and their gene-expression profile of normal liver and liver diseases from multiple perspectives (Table 2).50–62 As Wang et al.63 provided a comprehensive single-cell immune profile of HBV-infected patients, revealing that increased Treg and Tex cells are associated with the degree of liver damage, and that immune cell interactions networks may promote persistent HBV infection and liver histopathology. Thus, HBV infection is a dynamic process, and a deeper understanding of cellular interactions and the changes of immune microenvironment will help develop more effective therapeutic targets and biomarkers for immunotherapy in chronic hepatitis B (CHB) patients, thereby slowing down the occurrence of end-stage liver disease (Fig. 2).
Table 2Overview of recent studies on human liver disease at the single-cell level
Liver diseases | Samples | Patients | Methods | Synopsis | Year | Ref. |
---|
Hepatitis | Huh7.5-NTCP cells | none | scRNA-seq(Illumina HiSeq PE150); CRISPR-Cas9 | Rare transcripts from HBV-infected hepatocytes are preferentially enriched; HBV infection did not significantly alter gene expression but treatment with IFN-α altered gene expression in Huh7.5-NTCP cells | 2021 | 50 |
| Peripheral blood | 12 | scATAC-seq | Identified 12 leukocytic clusters corresponding to five cell types; specific cell type of CHB was B-0 and T-3 | 2022 | 51 |
| NASH liver tissue | 10 | scRNA-seq (10×) | Cell type-specific transcriptional features of NASH; two distinct populations of activated hepatic stellate cells associated with fibrosis levels | 2022 | 56 |
Liver fibrosis / Cirrhosis | Healthy, cirrhotic liver tissue; peripheral blood | 10 | scRNA-seq (10×) | Characterizing and comparing the functional diversity of liver biopsy cells from human fibrotic and normal livers; identified scar-associated TREM2+CD9+macrophages, ACKR1+ and PLVAP+ endothelial cells | 2019 | 52 |
| Healthy, HBV-related cirrhosis, alcoholic cirrhosis liver tissue | 15 | scRNA-seq (Illumina HiSeq X) | Different immune profiles between alcoholic and HBV cirrhotic patients; Galectin-9 may be a potential therapeutic target for alcoholic cirrhosis treatment | 2023 | 53 |
ACLF | Healthy, cirrhosis, ACLF liver tissue | 25 | scRNA-seq (Illumina HiSeq X) | Revealed the state of hepatic lymphatic endothelial cells in hepatitis B virus-associated ACLF; monocytes/macrophages infiltrate into ACLF | 2022 | 54 |
HCC | HCC tumor tissue | 8 | scRNA-seq (10×) | Cellular heterogeneity in HCC; intertumoral heterogeneity higher than intratumor heterogeneity; immunocyte characteristics | 2021 | 57 |
| HCC tumor tissue | 1 | scRNA-seq (SMART- seq) | Heterogeneity of CSCs in human HCC; identified CSCs as CD24+/CD133+/EpCAM+/CD45- cells | 2018 | 58 |
| HCC tumor tissue; normal liver tissue; peripheral blood | 8 | scRNA-seq (Illumina HiSeq) | Characteristics of immune cell in HCC; revealed specific cytokines/compartmental factors | 2019 | 59 |
| HCC, iCCA tumor tissue | 7 | scRNA-seq (10×) | Similar developmental trajectories in the same tumor, but heterogeneity in stemness gene expression | 2022 | 60 |
| HCC Tumor tissue; normal liver tissue; peripheral blood | 6 | scRNA-seq (Illumina HiSeq) | Identified 11 populations of T cells in HCC; TME; exhausted CD8+T cells and Tregs preferentially enriched in HCC | 2017 | 55 |
| HCC tumor tissue; normal liver tissue | 18 | scRNA-seq (10×) | Immune characteristics in early-relapse HCC; early-relapse tumors have less regulatory T cells, DCs and increased infiltrating CD8+T cells and higher expressed CD161+CD8+T cell compared to the primary tumors | 2021 | 61 |
| (HCC, iCCA, HH, CHC, ASC, SAR, SLC) Tumor tissue; normal liver tissue; peripheral blood | 124 | scRNA-seq (Illumina HiSeq) | Identified six tumor-associated neutrophils, among them CCL4+ TAN and PD-L1+ TAN could promote tumor growth | 2022 | 62 |
Normal liver
scRNA-seq can be used to construct the complete cellular map of the normal liver based on the information of human liver regarding zonation, development, regeneration, and immune cell heterogeneity. In 2017, Halpern et al.64 performed scRNA-seq to measure the entire transcriptome of thousands of mouse hepatocytes, combining it with marker gene function by smFISH to determine the positional coordinate characteristic of liver lobe positions. The researchers found that about 50% of liver genes had partitioning rules and obtained the spatial zonation profiles of the liver.64 In 2019, a group of researchers performed scRNA-seq to analyze approximately 10,000 cells in normal liver tissue from nine human donors and constructed a detailed cellular landscape of the human liver,47 which covered almost all important liver cell types, including hepatocytes, major liver metabolic cells, vascular endothelial cells, liver macrophages, and other immune cell types. A subsequent study integrated and analyzed scRNA-seq datasets from five independent studies of healthy livers from individuals of different sexes and ages and developed an interactively accessible online browser that could be used to access information on the expression of the gene of interest for each cell (sub) type from 28 healthy livers.65 ScRNA-seq combined with ST has provided a comprehensive understanding of the differentiation of all major cell types in the human fetal/embryonic liver, including the function of fetal liver hematopoietic stem cells.66,67 Such information can be used for developing effective strategies for treating liver diseases and enhancing liver regeneration. Other than that the developmental trajectories of various cell types in the human liver at the single-cell level have been described.68 In 2021, Inverso et al.69 combined transcriptomics and proteomics to analyze and obtain high-resolution data on liver vascular endothelial cells and their hepatic vascular signaling regulation. They found that the receptor tyrosine kinase Tie1 (a protein activation region) regulates vascular activity signals and strongly influences the maintenance of liver regeneration.
Immune cell proteomic analysis associated with liver damage and regeneration were developed later.70,71 In 2018, MacParland et al.72 isolated 8,444 parenchymal cells from five healthy human livers and then performed scRNA-seq analysis to obtain the macrophage population distribution of the human liver and mapped the human liver cell landscape. In addition, researchers also used comprehensive transcriptomes of unbiased scRNA-seq to study liver-resident immune cells (LrICs).73 That study showed the tissue distribution, gene expression, and functional modules of specific immune cell populations (natural killer cells, T cells, B cells, and plasma cells) in human liver and presented a landscape of LrICs to explain how the liver promotes immune defense and tolerance in vivo. Liver sinusoidal endothelial cells can function as host immune defenses by coordinating the proper distribution of intrahepatic resident immune cells.74 The integrated information on the liver and immunity by single-cell analysis can promote the development of new immunotherapies and targeted therapies, such as programmed cell death-1 (PD-1) for antitumor treatment.75 Thus, a complete intrahepatic immune signature is essential for understanding the development of liver disease and conditions under which intrahepatic compartmentalization extends to immune cells (Fig. 2).
Hepatitis B
Hepatitis, for various reasons, can lead to adverse outcomes, including cirrhosis, liver failure, or HCC. Among them, for viral hepatitis B, which is very common in China, which affects 257 million people worldwide, a suitable functional cure is lacking, even though effective viral suppression regimens have been developed.76,77 A study used scRNA-seq to analyze crosstalk between B cell and T cell subpopulations and blood samples from interferon-α (IFN-α)-treated CHB patients.78 The results showed that the immune profile of B cells differed across patients, and the effect of IFN-α treatment also differed. Based on the antigen-presenting and costimulatory functions of B cells in cellular immunity, B cell depletion was found to impair functional T-cell subsets, such as specific CD8+T cells, leading to the ineffective clearance of HBV.78 The ability of different subgroups of CD8+T cells to control the virus in patients with chronic viral hepatitis is closely related to the prognosis of the disease.79,80 The CRISPR-Cas9 technology integrated with scRNA-seq analysis can be combined to interpret host gene expression in HBV-infected hepatocytes, including sequence information of low-abundance genes.50 Interestingly, Mederacke et al.81 analyzed the TCR repertoire of CD4+T, CD8+T, and regulatory T cells from peripheral blood and the transplanted livers of liver transplant recipients and showed that their TCR pools were nearly distinct, and the activated T cells might be released by the graft in an acute cellular response to rejection. The latest study on the application of scATAC-seq in the pathophysiology of hepatitis B infection reported specific immune cell populations and chromatin-accessible specific molecules associated with the HBV infection through the detection of peripheral blood mononuclear cells (PBMCs); for example, the B-0 subgroup was found to be over activated in CHB patients.51
Liver fibrosis and cirrhosis
Liver fibrosis involves complex interactions between multiple mesenchymal cells located in scarred areas. These cells include mesenchymal cells, immune cells, vascular endothelial cells, etc. To determine the specific mechanisms of human liver fibrosis, Ramachandran et al.52 used scRNA-seq to analyze the cellular diversity of more than 105 human cells in normal and fibrotic livers of humans and found that TREM2+/CD9+macrophages and PDGFRα+collagen-producing myofibroblast subpopulations were involved in fibrosis. Additionally, profibrotic atypical chemokine receptor 1 (referred to herein as ACKR1) and plasmalemma vesicle-associated protein (referred to herein as PLVAP) endothelial cells were found to accumulate in the fibrotic ecotone and participate in the migration of immune cells, which present several profibrotic signaling pathways through receptor and ligand interactions. Such interactions, combined with newly identified fibrogenic markers, facilitate the development of rational therapeutic targets.
Cirrhosis is a major factor associated with mortality in nonalcoholic steatohepatitis (NASH). The application of scRNA-seq by Schwabe et al.82 provided new insights into the cellular networks involved in the development of NASH fibrosis. The cellular heterogeneity and intercellular transitions involved in the injurious fibrotic response are becoming clearer at the single-cell level.83 To investigate the immune characteristics of cirrhosis due to different causes, a recent study in China performed scRNA-seq on 15 liver specimens (five each from healthy controls, alcoholic cirrhotic (ALC) patients, and HBV-related cirrhotic patients) and found that the proportion of intrahepatic monocytes and macrophages was higher in the ALC group compared to that in the HBV group, with CD5L+Kupffer cells being the predominant subpopulation in ALC. The proportion of T and B cells was significantly lower, which might be related to the expression of Galectin-9 in Kupffer cells.53 Therefore, the immune profile of ALC patients differs from that of HBV cirrhotic patients, and Galectin-9 might be a therapeutic target for treating ALC.
Acute-on-chronic liver failure (ACLF)
Liver failure occurs due to serious liver damage caused by various factors, especially cirrhosis, which leads to severe impairment or loss of liver function. Thus, a comprehensive and systematic understanding of liver failure and the development of new treatment strategies are extremely important. In 2022, a study on the hepatic inflammatory milieu in ACLF found that hepatic lymphatic vessels (LVs) and lymphatic vessel endothelial cells (LyECs) might be new therapeutic targets for treating ACLF.54 That study used scRNA-seq to analyze nonparenchymal cells from three groups of human liver tissue (healthy participants and cirrhotic and ACLF patients) and found the development of apoptotic and dysfunctional LyECs caused by the release of secretory phosphoprotein 1 (also known as SPP1) from infiltrating monocytes/macrophages in the liver. The study also found that the number of LVs was significantly lower in the ACLF group than that in the cirrhotic group. A study on HBV- associated ACLF combined RNA-seq and differential gene analysis and identified key genes associated with immunometabolic disorders, including FGF19, ADCY8, and KRT17. The expression of these genes was significantly upregulated during the progression of ACLF, and thus, they might be sensitive biomarkers for HBV-associated ACLF.84 A transcriptomic study of PBMCs suggested that significant immune dysregulation leads to the development of ACLF from CHB.85 In a recent study, researchers performed scRNA-seq in 17,310 circulating monocytes from healthy & ACLF patients and found five monocyte subpopulations, including pro-inflammatory monocytes, CD16 monocytes, HLA monocytes, megakaryocyte-like monocytes, and natural killer-like monocytes; thrombospondin-1 (THBS1) from pro-inflammatory monocytes characterized specific types of cells and transcriptional signatures, and its production was speculated to reflect the progression of ACLF and a poor prognosis. A comparison between the ACLF survival group and the death group showed that inflammation-related cytokines were significantly higher in the ACLF death group. The results of that study indicated that a detailed analysis of the characteristic monocyte subpopulations and relevant proinflammatory molecular markers might provide therapeutic strategies to arrest the progression of ACLF.86 Shen et al.87 investigated the dynamics of T cells in ACLF longitudinally by collecting PBMCs from five HBV-ACLF patients and sorting total CD4+T and CD8+T cells. They performed unbiased high-throughput sequencing to analyze the TCRβ complementarity-determining regions 3 (also referred to as CDR3) and found that the diversity of the T cell pool decreased significantly as HBV-ACLF progressed. There was also a significant clonal expansion of CD8+T cells but not CD4+T cells. Thus, the diversity of the CD8+T-TCRβ pool might be a predictor of ACLF prognosis. To summarize, single-cell studies assisted by functional experiments might help in understanding the immunological aspect of liver failure and reveal new methods for treating this disease.
HBV-associated liver cancer
Tumorigenesis is a complex pathological process characterized by intratumoral and intertumor heterogeneity and changes in the immune microenvironment.88 HCC is the most prevalent liver tumor due to HBV infection. Here, we summarized the studies based on scRNA-seq related to human HCC in the last five years (Table 2). Tumor heterogeneity and the complexity of the immune microenvironment associated with HCC are the main reasons for its drug resistance and low clinical cure rate.89 To elucidate the immune landscape of HCC, Zheng et al.55 defined 11 functionally characterized T cell subsets and their characteristic genes by combining TCR-seq and scRNA-seq. The gene known as Layilin is functionally associated with tumor-characteristic CD8+T cells and might serve as a new target for immunotherapy in HCC. Wu et al.90 constructed a complete spatial ecological landscape of tumor nodules in a comprehensive genome-wide transcriptome analysis of primary HCC cells, including normal liver tissue, tumor margins, and tumor regions. They found that the tumor envelope might influence the spatial distribution and growth of tumor cells, and HCC heterogeneity can even exist in a single tumor nodule, which might contribute to the ineffectiveness of drug therapy for HCC. Tumor marginal cells are most responsive to their local microenvironment and are highly aggressive and metastatic.91 In 2022, the relationship between the immunophenotype of human liver-resident lymphocytes and HCC was determined by scRNA-seq.92 The researchers found that the changes in the HCC cell number and subpopulation were associated with the dysregulation of innate lymphoid cells (ILCs); seven characteristic ILCs were identified, including inducible T cell costimulator (ICOS)+ ILC2a, which were enriched in HCC and suggested a poor prognosis.
The accurate identification of specific cell subtypes and their targeted molecules is conducive to the diagnosis and treatment of liver diseases.93,94 For example, in the study of intrahepatic cholangiocarcinoma (iCCA), the tagging function of proteomics to complement other histological examinations was used to successfully identify three distinct functional molecular subtypes, of which POE+/C1QB+tumor-associated macrophages (commonly referred to as TAMs) represented a chronic inflammatory marker that might be developed as an immunotherapeutic target for iCCA.95 In a study on the association between HBV and iCCA, a proteomic analysis of the link between cells from HBV-iCCA patients found that a decrease in intercellular adhesion was more likely to lead to an increase in the epithelial-mesenchymal transition signaling.96 Additionally, abnormal methylation was found to be a marker for the progression of cirrhosis to HCC, and altered epigenetic modifications were found to affect chromatin structure and function in HCC patients. For example, the viral protein HBx was found to be an epigenetic regulator, playing an important role in HBV-associated HCC.97–99
Regarding metastatic liver cancer, an adequate understanding of the immune status of the tumor is needed to detect immune checkpoint suppression. By combining transcriptomic, epigenomic, and proteomic data from five HCC cell lines, Wang et al.100 showed that the metastatic potential of HCC is related to the mesenchymal state of tissues and the proliferation ability of cells. They also found that the hypoxic subtype increases the metastatic ability of cancer cells. The combination of multiomics can help identify the regulatory factors of tumor growth and metastasis and reveal the multifaceted nature of tumor cells and their interaction with the immune microenvironment, which might be used to find therapeutic targets.101,102 Yang et al.103 combined transcriptome and whole exome sequencing in primary and metastatic HCC and identified an immunosuppressive target known as SLC2A1, which acts by regulating SPP1+macrophages and T cells. To summarize, single-cell technologies provide valuable information for the diagnosis and treatment of liver cancers, and some of these strategies might be implemented soon in the clinic.
Challenges and future perspectives
The rapid advancement in single-cell sequencing has promoted the identification of the gene structure and gene expression state of single cells at extremely high resolution. With scRNA-seq having a key role, newly developed single-cell approaches have combined genome, epigenome, proteome, spatial transcriptome, and even multiomics to demonstrate the heterogeneity between cells and the changes in the microenvironment of cell populations. Such information might be used for the individualized diagnosis and treatment of diseases and to develop better precision medicine.104 By focusing on the progression of CHB to end-stage liver disease, our single-cell perspective on the mechanism of disease development can help identify relevant sensitive biomarkers and, thus, arrest disease progression on time.
Advances in technology and analytical techniques have encouraged further development in the field of single-cell technology, but research in this field is still at an early and has not been clinically accessible because of many challenges. First, single-cell studies require fresh liver tissue samples, and the availability and dissociation of samples greatly limit the ability to expand the spectrum of liver diseases.105 Although snRNA-seq can reduce the effects of enzymatic hydrolysis and mechanical pressure and is suitable for sequencing frozen samples and hard-to-dissociate samples,24,26,106 its sequencing depth and cell diversity are limited. Thus, snRNA-seq might not be suitable for studying immune cells. In 2023, Yu et al.107 developed a new technology known as “scONE-seq”, which can simultaneously sequence the DNA and RNA of single cells in frozen and fresh tissues. This technology was used for identifying one kind of cancer cell in the human brain with genomic mutations and near-normal transcriptomes. Another challenge involves the complexity of data analysis; specifically, defining quality control criteria, removing technical artifacts, and interpreting results are performed by a few technical staff members. Additionally, open source software tools are unsuitable for most clinical analyses. Therefore, in the future, data visualization and computational methods might be performed by artificial intelligence, which can be designed to analyze complex datasets quickly, thus allowing the networked dissemination of the results of single-cell analysis to a wider field. In addition, the current single-cell sequencing data visualization mainly applies principal component analysis, t-distributed stochastic neighbor embedding, and uniform manifold approximation and projection dimensionality reduction methods to project the cell population into two-dimensional space,108 but in fact, it is difficult for each dimension to correspond to the last biological structure after dimensionality reduction, which will lead to the loss of spatial information in a certain direction of the cell population, so mapping cells to three-dimensional space one by one is the future improvement of single-cell sequencing. Recently, in order to visualize different dataset types in a unified spatial framework, Voytek109 used a pure data approach of neurological mapping to study neuroscience, which opens new doors for other fields to study the spatial three-dimensional mapping of cells. More computational tools and related methods for ST data analysis and visualization are summarized in this review.110 It must be mentioned that the cost of sequencing is an important factor. Therefore, addressing common issues related to single-cell sequencing technologies, such as reducing sequencing costs, simplifying the sequencing process, enhancing sensitivity, and improving resolution, are areas that need to be focused on to further elucidate the molecular mechanisms of the development and progression of different liver diseases.
Novel single-cell techniques combined with multimodal, integrated analytical approaches have expanded the perspective of analysis at different levels within a cell. For example, the newly described technique for studying histone modifications in single cells, known as single-cell chromatin immune cleavage and unmixing sequencing (or scChIX-seq),111 overcomes the limitation that only one histone modification can be studied in a cell. This technology combines and compares the modification patterns of different histones in a cell to generate corresponding maps. The spatial information of gene expression with high resolution provides accurate information for analyzing the native environment of the lesion. The newly described Spatial Tumor Cell Estimator (commonly referred to as SpaCET)112 was used in 2023 to identify different types of immune and stromal cells from information on the spatial abundance of tumor cells. Thus, it can be used to study the role of immune cells in the progression of tumors. In addition, full-length transcriptome sequencing technology represented by single-molecule real-time sequencing can be used to obtain the full-length information of transcripts without splicing. Such information can be used to accurately identify gene mutations and find new mutants quickly.112 Surprisingly, the recently developed droplet-based technology detecting single nucleus information (also referred to as snRandom-seq) allows to capture full-length total RNAs with random primers and detect those noncoding RNA without polyA tails and nascent RNAs.113
The application of single-cell approaches for identifying and assessing HCC is the most prominent use of this technology. However, HCC is often the end-stage manifestation of chronic liver disease, and its prognosis is extremely poor; also, only a few patients are suitable for liver transplantation.114 In China, CHB, cirrhosis-HCC are part of a typical axis of chronic liver disease. Advancements in single-cell technologies are needed for detecting abnormal biomolecules in the early stages of chronic liver disease, especially in the early stages of hepatitis, to arrest the progression of the disease. Some of the latest studies show promise. For example, studies have shown that core transcription factors and pre-mRNA splicing regulators of the hepatocyte phenotype can predict the regression of chronic liver disease.115 Single-cell analysis in the liver involving parenchymal cells, immune cells, and stromal cells has elucidated the pathological mechanisms of viral infections, parasites, cancers, and other liver injuries.74 A more in-depth single-cell analysis of macrophages from patients with liver metastasis was performed to identify two markers associated with cancer.116 In the future, single-cell techniques should be combined with general clinical examinations, and such techniques should be applied beyond basic research to analyze the pathology and physiology of liver disease in patients with different characteristics. Such steps can facilitate precise and individualized treatment of patients with liver disease.
Conclusions
With scRNA-seq as a core sequencing technology, single-cell approaches have yielded significant discoveries in the physiology and pathology of the liver at an unprecedented high resolution. We comprehensively discuss the application and findings of the single-cell approach in the progression of CHB to liver cancer or liver failure. These continuously optimized single-cell techniques provide invaluable tools for exploring subpopulations of liver cells and developing new biomarkers. At present, experience in single-cell technology is most extensive in HCC, and we look forward to increased use in other liver diseases, reduced cost, and simplifying the procedures. One day, single-cell technology will become a basic tool for clinical diagnosis, achieving precise treatment of patients with liver disease.
Abbreviations
- 4D:
four-dimensional
- 10×:
10× Genomics
- ACLF:
acute-on-chronic liver failure
- ACKR1:
atypical chemokine receptor 1
- ALC:
alcoholic cirrhotic
- CCL4:
C-C motif chemokine ligand 4
- CHB:
chronic hepatitis B
- CHC:
combined hepatocellular-cholangiocarcinoma
- CITE-seq:
cellular indexing of transcriptomes and epitopes by sequencing
- CSCs:
cancer stem cells
- CUT & RUN:
cleavage under targets and release using nuclease
- DCs:
dendritic cells
- EpCAM:
epithelial cell adhesion molecule
- FFPE:
formalin-fixed and paraffin embedded
- HBV:
hepatitis B virus
- HCC:
hepatocellular carcinoma
- HH:
hepatic hemangioma
- iCCA:
intrahepatic cholangiocarcinoma
- IFN-α:
interferon-α
- ILCs:
innate lymphoid cells
- LVs:
lymphatic vessels
- LyECs:
lymphatic vessel endothelial cells
- nanoPOTS:
nanodroplet processing in one pot for trace samples
- NASH:
nonalcoholic steatohepatitis
- PBMCs:
peripheral blood mononuclear cells
- PD-1:
programmed cell death-1
- PDGFRα:
recombinant human platelet-derived growth factor-α
- PD-L1:
programmed cell death-ligand 1
- PLVAP:
plasmalemma vesicle-associated protein
- SAR:
sarcomatoid carcinoma
- scATAC-seq:
single-cell assay for transposase-accessible chromatin using sequencing
- scCGI-seq:
single-cell CpG island methylation sequencing
- scChIP-seq:
single-cell chromatin immune-precipitation sequencing
- scChIX-seq:
single-cell chromatin immune cleavage and unmixing sequencing
- scCOOL-seq:
single-cell chromatin overall omic-scale landscape sequencing
- scDNase-seq:
single-cell DNase sequencing
- sci-MET:
single-cell combinatorial indexing for methylation analysis
- SCoPE-MS:
single-cell proteomics by mass spectrometry
- scRNA-seq:
single-cell RNA sequencing
- scTrio-seq:
single-cell triple omics sequencing
- scWestern:
single-cell western blotting
- SLC:
secondary liver cancer
- smFISH:
single-molecule fluorescence in situ hybridization
- snRNA-seq:
single-nucleus RNA sequencing
- SPP1:
secretory phosphoprotein 1
- ST:
spatial transcriptome technology
- TAMs:
tumor-associated macrophages
- TAN:
tumors-associated neutrophil
- THBS1:
thrombospondin-1
- Tie 1:
tyrosine kinase with immunoglobulin like and EGF like domains 1
- TREM2:
triggering receptor expressed on myeloid cells-2
Declarations
Funding
This work was supported in part by the National Natural Science Foundation of China (32000892); the Guangdong Basic and Applied Basic Research Foundation (2023A1515010143) and Science and Technology Program of Guangzhou, China (202102010155).
Conflict of interest
The authors have no conflict of interests related to this publication.
Authors’ contributions
Made substantial contributions to conceptualizing, visualizing, drafting, reviewing, and editing the review article (GC, WX), made substantial contributions to reviewing of and funding acquisitions for the article (YJ, BL, WX), contributed to discussion and conceptualization (ZL, YC). All authors approved the final version of this manuscript.