- Research
- Open access
- Published:
Single-nucleus transcriptome sequencing reveals hepatic cell atlas in pigs
BMC Genomics volume 24, Article number: 770 (2023)
Abstract
Background
As the largest substantive organ of animals, the liver plays an essential role in the physiological processes of digestive metabolism and immune defense. However, the cellular composition of the pig liver remains poorly understood. This investigation used single-nucleus RNA sequencing technology to identify cell types from liver tissues of pigs, providing a theoretical basis for further investigating liver cell types in pigs.
Results
The analysis revealed 13 cells clusters which were further identified 7 cell types including endothelial cells, T cells, hepatocytes, Kupffer cells, stellate cells, B cells, and cholangiocytes. The dominant cell types were endothelial cells, T cells and hepatocytes in the liver tissue of Dahe pigs and Dahe black pigs, which accounts for about 85.76% and 82.74%, respectively. The number of endothelial cells was higher in the liver tissue of Dahe pigs compared to Dahe black pigs, while the opposite tendency was observed for T cells. Moreover, functional enrichment analysis demonstrated that the differentially expressed genes in pig hepatic endothelial cells were significantly enriched in the protein processing in endoplasmic reticulum, MAPK signaling pathway, and FoxO signaling pathway. Functional enrichment analysis demonstrated that the differentially expressed genes in pig hepatic T cells were significantly enriched in the thyroid hormone signaling pathway, B cell receptor signaling pathway, and focal adhesion. Functional enrichment analysis demonstrated that the differentially expressed genes in pig hepatic hepatocytes were significantly enriched in the metabolic pathways.
Conclusions
In summary, this study provides a comprehensive cell atlas of porcine hepatic tissue. The number, gene expression level and functional characteristics of each cell type in pig liver tissue varied between breeds.
Background
Pigs are one of the most farmed livestock in the world. The meat of pigs is widely eaten by people across the world [1]. At the same time, the physiological structure of pigs is similar to that of humans, so it is an ideal animal model [2]. The liver is the largest digestive gland in the animal body and plays an important role in the physiological processes of digestion, metabolism, and immune defense [3,4,5]. Furthermore, the liver is a visceral organ that is capable of remarkable natural regeneration after tissue loss and retain basic metabolic functions [6,7,8]. However, the liver cellular landscape has barely been explored at single-cell resolution, which limits our molecular understanding of pig liver structure and function.
Traditional transcriptomics detects the average value of all cells in a sample, and it is difficult to give a clear answer at the cellular scale [9]. Single-nucleus RNA sequencing (snRNA-seq), enables robust and unbiased exploration of individual cell states and types, yielding new insights into tissue biology [10,11,12]. Therefore, snRNA-seq is a promising approach to investigate the transcriptome of individual cells in tissue [13]. Recent studies have shown that snRNA-seq has emerged as a complementary approach to investigate complex tissues at single-cell level in mouse and human samples [14,15,16,17,18,19,20]. However, the application of snRNA-seq in pigs is relatively limited. Recent snRNA-seq analysis were focused on the lung, cerebral cortex, hypothalamus, and peripheral blood in pigs [21,22,23]. There is limited information on the liver atlas of pigs.
Therefore, this study used snRNA-seq to identify cell types, explore the biological function of pig liver deeply, and describe the transcriptome characteristics of cells comprehensively. The findings revealed the transcriptional landscape of the pig liver and the effect of pig breeds on liver cell types and functions, providing a theoretical basis for future in-depth research on pig liver function.
Materials and methods
Animal feed and liver tissue collection
In the Dahe Black Pig Research Institute of Fuyuan County, six male Dahe pigs and six male Dahe black pigs of the same batch with similar body weight were randomly selected. Pigs of the same breed come from the same maternal parent and parity. Each pig was raised in a single pen and all pigs were fed the same basic diet without any antibiotics. When pigs reached 194 days of age, the right lobe of the liver was collected from every animal, a part was fixed with 4% paraformaldehyde fixative solution for histological verification, and the rest was frozen in liquid nitrogen for the preparation of single nuclear suspension and RNA extraction.
Liver morphology and cell number
Liver samples were removed from the 4% paraformaldehyde fixative solution and embedded in paraffin. Each liver was sliced to a thickness of 4 mm. From each block, two consecutive histological sections randomly positioned within the block and mounted on adhesion microscope slides. Sections were stained with hematoxylin and eosin (Beijing Solarbio Science&Technology Co., Ltd., Beijing, China) for analysis. Periodic acid-Schiff staining (PAS) (Shanghai Yuanye Biotechnology Co., Ltd., Shanghai, China) was used to highlight the contours of the individual cell and to count the endothelial cell, hepatocyte, and lymphocyte.
Cell suspension preparation
The liver tissue samples for snRNA-seq were selected from one male Dahe pigs and one male Dahe black pigs with body weights closest to the average body weight. Samples were subjected to nuclear isolation, sequencing, and library preparation following the 10X Genomics protocol. Approximately 500 mg of hepatic tissue was dissociated into a singular nuclear suspension via tissue homogenization in a chilled lysis buffer (0.25 M sucrose, 5 mM CaCl2, 3 mM MgAc2, 10 mM Tris–HCl pH 8.0, 1 mM DTT, 0.1 mM EDTA, 1 × Protease Inhibitor and 1U/µL RiboLock RNase Inhibitor (Thermo Scientifc, cat no. O0381) with pestle strokes. The resulting homogenate was subsequently filtered through a 70 µm cell strainer, yielding a nuclear fraction collected in a 50 ml centrifuge tube, with a 1 ml volume. This nuclear fraction was mixed with an equal volume of 50% iodixanol solution (0.16 M sucrose, 10 mM NaCl2, 3 mM MgCl2, 10 mM Tris–HCl pH 7.4, 1 U/µL RiboLock RNase Inhibitor, 1 mM DTT and 0.1 mM PMSF Protease Inhibito (Thermo Scientifc, cat no. 36978), yielding a final concentration of 25%, and supplemented with 1 mL of 33% iodixanol solution at the tube's base and 30% iodixanol solution at the top. The solution underwent inversion mixing 10 times before centrifugation at 500 × g for 8 min at 4℃, subsequent to myelin layer removal from the gradient’s apex. Nuclei were harvested from the 30% iodixanol interface, resuspended in nuclear wash and resuspension buffers, and centrifuge (0.04% bovine serum albumin, 0.2 U/µL RiboLock RNase inhibitor, 500 mM mannitol and 0.1 mM PMSF protease inhibitor in PBS) at 500 × g and 4℃ for 5 min. Filtration through a 40 µm cell filter eliminated cell debris and sizable aggregates. The nuclei’s total count, concentration, and integrity ratio were ascertained via hemocytometer-assisted microscopic examination of trypan-stained samples. Ultimately, the nuclear concentration was regulated to 700–1,200 nuclei/µL, with nuclei inspected using the 10X Chromium platform.
snRNA-seq library preparation and sequencing
Cellular suspensions were loaded on a 10X Genomics GemCode Single-cell instrument that generates single-cell Gel Bead-In-EMlusion (GEMs). Using Chromium Next GEM Single Cell 3'Reagent Kit v3.1 for Library Generation and Sequencing of cDNA. Upon dissolution of the Gel Bead in a GEM, primers containing (i) an Illumina® R1 sequence (read 1 sequencing primer), (ii) a 16 nt 10 × Barcode, (iii) a 10 nt Unique Molecular Identifier (UMI), and (iv) a poly-dT primer sequence were released and mixed with cell lysate and Master Mix [24].
Silane magnetic beads were used to remove leftover biochemical reagents and primers from the post GEM reaction mixture. Full-length, barcoded cDNAs were then amplified by PCR to generate sufficient mass for library construction. R1 (read 1 primer sequence) were added to the molecules during GEM incubation. P5, P7, a sample index, and R2 (read 2 primer sequence) were added during library construction via End Repair, A-tailing, Adaptor Ligation, and PCR. The final libraries contained the P5 and P7 primers used in Illumina bridge amplification. A Single Cell 5’ Library comprised standard Illumina paired-end constructs which begin and end with P5 and P7. The Single Cell 5′ 16 bp 10X Barcode and 10 bp UMI were encoded in Read 1, while Read 2 was used to sequence the cDNA fragment. Sample index sequences were incorporated as the i7 index read. Read 1 and Read 2 were standard Illumina® sequencing primer sites used in paired-end sequencing [25].
snRNA-seq data processing and quality control
FASTQ files were processed with 10X Genomics Cell Ranger (v3.1.0) using the default parameters and aligned to the Ensembl_release104 Sscrofa 11.1 (NCBI Accession AEMK00000000.2) reference genomes. The gene expression matrix was processed and analyzed by Seurat package (v3.2.1). To filter out low-quality cells, the investigation removed reads from cells in which less than 240 and over 3,600 genes were detected, cells with more than 16,000 UMIs, and cells with more than 10% share of mitochondrial gene expression.
Cell clustering by Seurat
Using the R package Seurat [26], multiple criteria were employed to filter cells, eliminating multi-cellular entities and cells with suboptimal states, while retaining high-quality cells. The gene expression levels were normalized using a log transformation method, with “mitochondrial genes” and “cell cycle scoring” as regression variables. The resulting high-quality cell data was subjected to data integration and batch effect correction using Harmony [27]. Soft k-means clustering was applied to the dimension-reduced data.
To mitigate noise interference from individual gene expression levels in single-cell transcriptomics data, the vst algorithm was applied for high-variable gene selection. Principal Component Analysis (PCA) was performed on the integrated data to reduce the dimensionality to 50 principal components for capturing the main information. The parameter k.param was set to 20, and the Annoy algorithm with Euclidean distance metric was employed to calculate cell-to-cell distances. The Louvain [28] algorithm was used for clustering the dimension-reduced data with a resolution parameter set to 50, partitioning all cells in the single-cell transcriptome into distinct cell subpopulations, facilitating subsequent analyses.
Differentially expressed genes (DEGs)
Expression value of each gene in given cluster were compared against the rest of cells using Wilcoxon rank sum test [29]. Significant up-regulated and down-regulated genes were identified using a number of criteria. First, genes had to be at least 1.28-fold overexpressed in the target cluster. Second, genes had to be expressed in more than 25% of the cells belonging to the target cluster. Third, P-value is less than 0.05. The False Discovery Rate (FDR) correction was used for multiple test corrections.
Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) analyses of cell types
To annotate the function of these DEGs, GO analysis was conducted by using the GOseq software for each of the three main categories: biological process, cellular component and molecular function. Firstly, all peak related genes were mapped to GO terms in the Gene Ontology database (http://www.geneontology.org/), gene numbers were calculated for every term, significantly enriched GO terms in differentially expressed genes comparing to the genome background were defined by hypergeometric test. The calculating formula of P-value is:
Here N is the number of all genes with GO annotation; n is the number of differentially expressed genes in N; M is the number of all genes that are annotated to the certain GO terms; m is the number of differentially expressed genes in M. The calculated P-values were corrected using FDR, taking FDR ≤ 0.05 as a threshold. GO terms meeting this condition were defined as significantly enriched GO terms in differentially expressed genes.
Genes usually interact with each other to play roles in certain biological functions. Pathway- based analysis helps to further understand genes biological functions. KEGG is the major public pathway-related database [30,31,32]. KEGG pathway enrichment analysis identified significantly enriched metabolic pathways or signal transduction pathways in differentially expressed genes comparing with the whole genome background results [33]. The calculating formula is the same as that in GO analysis. Here N is the number of all transcripts that with KEGG annotation, n is the number of differentially expressed genes in N, M is the number of all transcripts annotated to specific pathways, and m is number of differentially expressed genes in M. The calculated P-values were corrected using FDR, taking FDR ≤ 0.05 as a threshold. Pathways meeting this condition were defined as significantly enriched in differentially expressed genes.
Statistical analysis
Statistical analysis was performed using GraphPad Prism (v8.4.0) [34]. The Shapiro–Wilk test, Wilcoxon rank-sum test, Pearson’s chi-square test, Unpaired student t test and Log-rank test were used in this study. Hazard ratio (HR) and confidence interval (CI) were calculated from Cox proportional hazards regression models. All statistical tests were two-sided, and significant differences between each two groups were indicated by *P < 0.05 and **P < 0.01.
Results
Liver tissue structure and cell number
The body weight of the Dahe pigs (DH) were higher than the Dahe black pigs (DHB), while the liver weight and liver index of DH were lower than the DHB (P < 0.05) (Fig. 1a). Conversely, there was no significant difference in crude fat (EE) content between Dahe pig liver (DHL) and Dahe black pig liver (DHBL) (P > 0.05) (Fig. 1b). Liver tissues samples displayed a normal structure without obvious inflammatory cell infiltration and fibrous tissue deposition in the portal area (Fig. 1c). The number of liver tissue cells was verified using PAS, focusing on endothelial cells, lymphocytes, and hepatocytes, which were the highest numbers. Among them, there was no significant difference in the number of endothelial cells and hepatocytes of the DHL and DHBL, while the number of lymphocytes was significantly higher in the DHBL than that in the DHL (P < 0.05) (Fig. 1d).
Landscape of single-nucleus transcriptome in the porcine hepatic tissue
The study used snRNA-seq of liver samples to characterize cellular heterogeneity (Fig. 2a). 1,647 and 1,481 cells were captured from DHL sample and DHBL sample for library construction and paired-end sequencing. Reads from the cells with gene numbers less than 240 and over 3,600 (Fig. 3a), with more than 16,000 UMIs (Fig. 3b), and with more than 10% mitochondrial gene expression (Fig. 3c) were removed. Finally, 1,355 and 1,171 high-quality cells were obtained from the DHL and DHBL (Fig. 3d-f). 407,797,286 and 403,918,586 sequencing reads were obtained from DHL and DHBL, with 96.9% valid barcodes. 91.5% of the sequencing reads could be mapped to the DHL genome, and 92.3% of the sequencing reads could be mapped to the DHBL genome (Table S1). This study identified 23,402 genes in these 1,355 cells of the DHL, with 1,890 median genes per cell, and the median UMI count per cell was 3,767. Moreover, the study identified 23,102 genes in these 1,171 cells of DHBL, with 1,902 median genes per cell, and the median UMI count per cell was 3,934 (Table 1).
This study classified cell types for all samples together and based on T-distributed stochastic neighbor embedding (t-SNE) dimensionality reduction and unsupervised cell clustering. Thirteen cell clusters were identified based on the expressed unique transcriptional profiles (Fig. 2b). The number of expressed genes in each cluster ranged from 9,141 to 20,115 (Table S1). The cell numbers distributed in each cluster ranged from 17 to 822 (Fig. 2c) and exhibited differences between the DHL and DHBL (Fig. 2d, Table S2). Clusters 0 and 1 were the two clusters with higher number of cells in the DHL and DHBL. Cell number of clusters 0, 2, 5, 6, 7, 9, and 11 in the DHL were higher than that in the DHBL, while Cell number of clusters 1, 3, 4, 8, 10, and 12 in the DHL were lower than that in the DHBL. Additionally, the gene expression heat map was generated for the top 5 marker genes in 13 clusters and identified the genes which were significantly enriched in each of the 13 clusters (Fig. 2e and Table S3). The results indicated certain demarcation boundaries between each cluster (Fig. 2e). Among them, cluster 1 and 2 were inferred to be a cell type preliminarily. Cluster 3 and 5 were inferred to be a cell type preliminarily.
Identification of cell types in the porcine liver tissue
Subsequent annotations were made based on the expression of cell type specific marker genes in relevant liver studies [35,36,37,38,39,40]. The thirteen cell clusters were identified as seven cell types (Fig. 4a), including endothelial cells, T cells, hepatocytes, Kupffer cells, stellate cells, B cells, and cholangiocytes (Fig. 4b). Endothelial cells were marked by adhesion G protein-coupled receptor (ADGRF5) and kinase insert domain receptor (KDR), and T cells were marked by cluster of differentiation 8 subunit alpha (CD8A) (Figs. 4c, d). Hepatocytes and Kupffer cells were marked by glucose-6-phosphatase catalytic subunit 1 (G6PC1) and cluster of differentiation 163 (CD163), respectively. Stellate cells and B cells were marked by collagen type I alpha 1 chain (COL1A1) and BCL11 transcription factor A (BCL11A), respectively. Cholangiocytes were marked by PKHD1 ciliary IPT domain containing fibrocystin (PKHD1) (Fig. 4c, d). The expression levels of each marker gene in each cell type are shown in Table 2.
According to the circle plots, clusters 0, 6, 7, and 11 were classified as endothelial cells. Clusters 1, 2, and 10 were classified as T cells. Clusters 3 and 5 were classified as hepatocytes. Cluster 4 and 8 were classified as Kupffer cells and stellate cells, respectively. Clusters 9 and 12 were classified as B cells and cholangiocytes, respectively (Fig. 4e). A gene expression heat map of marker genes across the seven cell types was generated (Fig. 4f, Table S4). In addition, potential new marker genes for each cell type were found basing on their the highest gene expression level, such as LIM domain binding 2 (LDB2) and RNA binding motif single stranded interacting protein 3 (RBMS3) in endothelial cells, dedicator of cytokinesis 10 (DOCK10) and src kinase associated phosphoprotein 1 (SKAP1) in T cells, carbamoyl-phosphate synthase 1 (CPS1) and alpha-1-antichymotrypsin 2 (SERPINA3-2) in hepatocytes, pleckstrin and Sec7 domain containing 3 (PSD3) in Kupffer cells, matrix Gla protein (MGP) and four and a half LIM domains 2 (FHL2) in stellate cells, EBF1 and histone deacetylase 9 (HDAC9) in B cells, and dachshund family transcription factor 1 (DACH1) in cholangiocytes (Figs. 5, 6, Table S5).
Number and functional enrichment of cells in the porcine liver tissue
The study identified 699 endothelial cells, 244 T cells, 219 hepatocytes, 110 Kupffer cells, 35 stellate cells, 41 B cells, and 7 cholangiocytes in the DHL. In contrast, the study indentified 326 endothelial cells, 421 T cells, 222 hepatocytes, 114 Kupffer cells, 46 stellate cells, 32 B cells, and 10 cholangiocytes in the DHBL (Fig. 7a, b, Table S6). There were differences in the number of each cell type (Fig. 7c, Table S6). Moreover, there were DEGs between each cell type, and the number of up-regulated and down-regulated genes were shown (Fig. 7d). DEGs were obtained by comparing the gene expression between DHBL and DHL, with FC > 1.28 and p value < 0.05. In this study, endothelial cells, T cells and hepatocytes were the dominant cell types in the liver tissue of two pig breeds. These three cell types accounted for 85.76% and 82.74% of the DHL and DHBL, respectively. The up-regulated gene numbers in endothelial cells, T cells and hepatocytes were 985, 556, and 744 in DHBL compared to DHL, respectively. The down-regulated gene numbers in endothelial cells, T cells and hepatocytes were 959, 372, and 815 in DHBL compared to DHL, respectively (Table S7). To gain further understanding of the difference in DHL and DHBL, the study focused on the top 20 pathways.
In endothelial cells, 1,944 DEGs were identified and a volcano plot of DEGs was generated (Fig. 8a). A total of 694 significantly enriched GO terms were identified, including 48 cellular component terms, 527 biological process terms, and 64 molecular function terms (Table S8). Among the top 20 terms ranked in significance, 13 terms belong to cellular components, 4 terms belong to biological process, and 3 terms belong to molecular function (Fig. 8a). A total of 20 significantly enriched KEGG pathway were identified, 10 pathways were significantly down-regulated in endothelial cell of DHBL compared to DHL, including focal adhesion, endocytosis, phosphatidylinositol signaling system, MAPK signaling pathway, FoxO signaling pathway, sphingolipid signaling pathway, phospholipase D signaling pathway, protein processing in endoplasmic reticulum, platelet activation, Fc gamma R-mediated phagocytosis, aldosterone synthesis and secretion (Table S9).
In T cells, 928 DEGs were identified and a volcano plot of DEGs was generated (Fig. 8b). A total of 498 significantly enriched GO terms were identified, including 51 cellular component terms, 399 biological process terms, and 48 molecular function terms (Table S10). Among the top 20 terms ranked in significance, 10 terms belong to cellular components, 8 terms belong to biological process, and 2 terms belong to molecular function (Fig. 8b). A total of 34 significantly enriched KEGG pathway were identified, 9 pathways were significantly up-regulated in T cells of DHBL compared to DHL. Only MAPK signaling pathway was significantly down-regulated in the T cell of the DHBL compared to the DHL (Table S11).
In hepatocytes, 1,559 DEGs were identified and a volcano plot of DEGs was generated (Fig. 8c). A total of 484 significant enriched GO terms were identified, including 77 cellular component terms, 328 biological process terms, and 79 molecular function terms (Table S12). Among the top 20 terms ranked in significance, 12 terms belong to cellular components, 7 terms belong to biological process, and 1 term belong to molecular function (Fig. 8c). A total of 30 significantly enriched KEGG pathway were identified, 27 pathways were significantly up-regulated in hepatocytes of DHBL compared to DHL, including peroxisome, ferroptosis, metabolic pathways. Protein processing in the endoplasmic reticulum and protein export were down-regulated in DHBL hepatocytes compared to DHL, respectively (Table S13).
Discussion
The liver serves as an important metabolic organ and affects a variety of physiological functions in pigs. The molecular characteristics and functional properties of organs are determined by their constituent cell types [20]. Nevertheless, the cellular composition of pig liver remains inadequately understood. The recent advent of sensitive snRNA-seq methods has enabled the research of cell types in animal tissues [41]. To date, single-nucleus studies of human and mouse liver have been reported [14, 15, 42,43,44], but research of single-nucleus in porcine hepatic tissue have not been documented yet. In future research efforts, the importance of pigs is rapidly growing, because of their high homology with humans. Pigs are not only economically essential livestock, but also in anatomy, physiology, biochemistry, and drug metabolism [45,46,47,48]. Consequently, in this study, the snRNA-seq system was used to analyze cell types and functions in porcine liver tissue.
At single-cell resolution, pig liver tissue was found to contain endothelial cells, T cells, hepatocytes, Kupffer cells, stellate cells, B cells and cholangiocytes. These seven cell types are also present in human and mouse livers, suggesting conservation of major liver tissue components across species [49]. Endothelial cells form the inner wall of hepatic sinusoid, facilitating the passage of blood components and promoting cellular uptake of essential substances and secretion discharge [50]. T cells and B cells perform immune functions, while hepatocytes participate in various biological processes, from protein synthesis and lipid metabolism to detoxification [43]. Stellate cells have many functions, such as storing vitamin A, regulating hemodynamics regulation, supporting liver regeneration, and regulating immunity [51]. Kupffer cells are specialized hepatic macrophages [52], whereas cholangiocytes are bile duct-specific epithelial cells involved in bile secretion. The study revealed specific marker genes for each cell type. For example, DOCK10 was abundantly expressed in normal T cells, and SKAP1 acted as an immune cell adaptor, connecting the T cell receptor to LFA-1-facilitated “inside-out” signaling involved in T-cell adhesion. As a result, these genes are considered as marker genes for T cells [53, 54]. This facilitates more complex cell classification and provides new perspectives for future research on hepatic cells.
Endothelial cells, T cells and hepatocytes represent the abundant cell populations in this study. A report on mouse livers demonstrated that the highest number of hepatocytes [20]. However, in this study, endothelial cells were the most abundant, suggesting that their numbers may vary, although hepatic cell types exhibit a degree of conservation across species. Research has indicated that endothelial cells participate in angiogenesis, contraction, and vasodilation process [55]. ADGRF5 and KDR genes are essential for vascular development and maintenance. In this study, ADGRF5 and KDR genes were highly expressed in endothelial cells, leading us to speculate that porcine hepatic tissue may be more involved in vascular-related physiological processes. There are studies on the classification of cells in mouse liver, where endothelial cells use the CD31 gene as a surface marker [56]. MacParland observed three endothelial cell populations, among which the most abundant endothelial cell cluster displayed enriched expression of F8, PECAM1 [57]. The marker genes were used to identify liver cell types, which came from the reports in the literature and first identified in this study.
In this study, Dahe pigs were used, which are representative of typical local breeds in Yunnan, China. Dahe black pigs are a crossbreed using the Duroc × Dahe breeding scheme through five generations of selection [58]. Both varieties have the advantages of tender meat and delicious taste [59]. However, Dahe pigs grow slower and have a lower lean meat percentage, compared to Dahe black pigs [60, 61]. This study revealed differences in the cell number between the liver tissue of Dahe pig and Dahe black pig. The number of endothelial cells in liver tissue of Dahe pig was higher than that of Dahe black pig, and the opposite was true for T cells. Liver tissues of two pig breeds were analyzed by PAS staining. The numbers of endothelial cells were similar because the study observed the interior of hepatic lobules. Endothelial cells were primarily concentrated in hepatic sinusoids and interlobular connective tissue. In this area, the study observed significant differences between Dahe pigs and Dahe black pigs. The results for hepatocytes and lymphocytes were consistent with the sequencing findings, indicating that the number of hepatocytes in the liver was not significantly different, while the number of lymphocytes in Dahe black pigs was significantly higher than that in Dahe pigs.
The liver is a major metabolic organ, and its endothelial cells play a role in clearance. Sinusoidal endothelial cells, one of the most endocytotic cells in humans, display multiple scavenger receptors on their cell surface [62]. Therefore, the clearance of extracellular material depends on the health of endothelial cells and the degree of endocytic function [63]. Hepatic endothelial cells also effectively regulate the exchange of substances between hepatic sinusoidal blood flow and surrounding tissues [4]. Recent studies have demonstrated that endothelial cells and the immune environment was critical for hepatic homeostasis [64, 65]. This study provides a holistic new perspective on the processes that endothelial cell function. Endothelial cells adherently senses exogenous substances, which are then endocytosed and absorbed into the cell. Subsequently, autophagy occurs through signal transduction pathways, facilitating growth, metabolism, and immunity functions. This study revealed that the number of endothelial cells in Dahe pig liver tissue was higher than in Dahe black pig. Among the 20 pathways enriched with differentially expressed genes in hepatic endothelial cells of Dahe pig liver and Dahe black pig liver, 3 pathways were up-regulated, and 10 pathways were down-regulated of Dahe black pig liver tissue compared to Dahe pig. This observation may indicate a positive correlation between cell number and functional diversity. This study indicates that changes in the number of hepatic cells were observed after pig hybridization.
The animal liver, an organ with multiple immune functions, serves as a sentinel for the human immune system [66]. In recent years, the crosstalk between the liver and the immune system has been uncovered through the study of hepatic snRNA-seq [67, 68]. T cells play a central role in adaptive immunity, and their activation involves spatially and temporally coordinated signaling processes across multiple time and length scales [69]. This study identified T cells in the livers of two pig breeds, with T cells constituting the second largest cell population. Notably, the number of T cells in Dahe black pig liver tissue was higher than in Dahe pig the immune and endocrine functions of the liver of Dahe black pigs are potentially stronger. In the future, we need to do more research to prove this conclusion.
Hepatocytes are the major parenchymal cells of the liver. Hepatocytes play critical roles in liver homeostasis and disease development [70]. These cells are responsible for the majority of hepatic metabolic, biosynthetic, biodegradable, and secretory functions [71]. And hepatocytes are involved in biological processes ranging from protein synthesis and lipid metabolism to exogenous and endogenous detoxification [43]. Many functions of hepatocytes require close cooperation between cell adhesion molecules, cell junctions, cytoskeleton, extracellular matrix, and intracellular trafficking machinery [72]. The study revealed only small differences in the number of hepatocytes in the liver tissue of Dahe pig and Dahe black pig. Among the 30 pathways enriched with significantly different genes in hepatocytes of Dahe pig liver and Dahe black pig liver, 27 pathways were up-regulated, and 2 pathways were down-regulated of Dahe black pig liver tissue compared to Dahe pig. These results indicated that hepatocytes of Dahe black pigs exhibited more abundant metabolic and biosynthetic functions, especially in protein metabolism, lipid metabolism, amino acid metabolism, and carbohydrate metabolism. This is supported by a recent report demonstrating specificity in liver lipid metabolism between different pig breeds [73]. The findings highlight breed differences in liver tissue between Dahe pigs (a Chinese indigenous breed) and Dahe black pigs (a crossbreed using the Duroc × Dahe breeding scheme through five generations of selection). Variations in the number of cell types were detected between the two pig breeds. Considering the functional enrichment of differentially expressed genes, it can be inferred that these cell types may exhibit mutual compensatory effects.
A limitation of this study is that the small number of liver samples sequenced. So, interpretations of some functions of the liver may not generalize. In essence, leveraging single-nucleus resolution and cost-effective UMI-based approaches to enhance throughput, the study findings also furnish foundational data for future investigations in this field. In the future, functional validation of target genes will be carried out, with the initial goal of conducting validation on cells.
Conclusions
In summary, this study establishes a comprehensive single-nucleus atlas landscape and identifies potential new marker genes for each cell type. The number, gene expression level and functional characteristics of each cell type in pig liver tissue varied between breeds. These high-quality snRNA-seq data serve as a valuable resource for future studies on porcine hepatic function and may provide informative support for human hepatic health and immunity.
Availability of data and materials
Data is contained within the article or supplementary material. Additional data that support the findings of this study are available from the corresponding author upon reasonable request. The accession number for the liver single nucleus RNA sequencing data reported in this paper is PRJNA956984.
Abbreviations
- ADGRF5 :
-
Adhesion G protein-coupled receptor
- BCL11A :
-
BCL11 transcription factor A
- CD8A :
-
Cluster of Differentiation 8 subunit alpha
- CD163 :
-
Cluster of Differentiation 163
- CPS1 :
-
Carbamoyl-phosphate synthase 1
- COL1A1 :
-
Collagen type I alpha 1 chain
- DH:
-
Dahe pig
- DHB:
-
Dahe black pig
- DHL:
-
Dahe pig liver
- DHBL:
-
Dahe black pig liver
- DEGs:
-
Differentially expressed genes
- DOCK10 :
-
Dedicator of cytokinesis 10
- DACH1 :
-
Dachshund family transcription factor 1
- EBF1 :
-
EBF transcription factor 1
- FDR:
-
False Discovery Rate
- FHL2 :
-
Four and a half LIM domains 2
- GO:
-
Gene ontology
- GEMs:
-
Gel Bead-In-EMlusion
- G6PC1 :
-
Glucose-6-phosphatase catalytic subunit 1
- HDAC9 :
-
Histone deacetylase 9
- KEGG:
-
Kyoto Encyclopedia of Genes and Genomes
- KDR :
-
Kinase insert domain receptor
- LDB2 :
-
LIM domain binding 2
- MGP :
-
Matrix Gla protein
- PAS:
-
Periodic acid-Schiff staining
- PCA:
-
Principal component analysis
- PSD3 :
-
Pleckstrin and Sec7 domain containing 3
- PKHD1 :
-
PKHD1 ciliary IPT domain containing fibrocystin
- RBMS3 :
-
RNA binding motif single stranded interacting protein 3
- SKAP1 :
-
Src kinase associated phosphoprotein 1
- snRNA-seq:
-
Single-nucleus RNA sequencing
- SNN:
-
Shared-nearest neighbour
- SERPINA3-2 :
-
Alpha-1-antichymotrypsin 2
- tSNE:
-
T-distributed stochastic neighbor embedding
- UMI:
-
Unique Molecular Identifier
References
Livestock of the World. Pigs. 2023. http://www.livestockoftheworld.com/Pigs/. Accessed 29 June 2023.
Pabst R. The pig as a model for immunology research. Cell Tissue Res. 2020;380(2):287–304.
Krenkel O, Tacke F. Liver macrophages in tissue homeostasis and disease. Nat Rev Immunol. 2017;17(5):306–21.
Poisson J, Lemoinne S, Boulanger C, Durand F, Moreau R, Valla D, et al. Liver sinusoidal endothelial cells: Physiology and role in liver diseases. J Hepatol. 2017;66(1):212–27.
Monga SP, Ko SJ, Hu SK, Singh SC, Poddar MS. Single-cell spatial transcriptomics analysis of a regenerating liver. FASEB J. 2022;36(S1):R3631.
Shi JH, Line PD. Hallmarks of postoperative liver regeneration: An updated insight on the regulatory mechanisms. J Gastroenterol Hepatol. 2020;35(6):960–6.
Racanelli V, Rehermann B. The liver as an immunological organ. Hepatology. 2006;43(2):S54-62.
Michalopoulos GK. Hepatostat: Liver regeneration and normal liver tissue maintenance. Hepatology. 2017;65(4):1384–92.
Chambers DC, Carew AM, Lukowski SW, Powell JE. Transcriptomics and single-cell RNA-sequencing. Respirology. 2019;24(1):29–36.
Kulkarni A, Anderson AG, Merullo DP, Genevieve K. Beyond bulk: a review of single cell transcriptomics methodologies and applications. Curr Opin Biotechnol. 2019;58:129–36.
Moncada R, Barkley D, Wagner F, Chiodin M, Devlin JC, Baron M, et al. Integrating microarray-based spatial transcriptomics and single-cell RNA-seq reveals tissue architecture in pancreatic ductal adenocarcinomas. Nat Biotechnol. 2020;38:333–42.
MacParland SA, Liu JC, Ma XZ, Innes BT, Bartczak AM, Gage BK, et al. Single cell RNA sequencing of human liver reveals distinct intrahepatic macrophage populations. Naturre Community. 2018;9:4383.
Krishnaswami SR, Grindberg RV, Novotny M, Venepally P, Lacar B, Lacar K, et al. Using single nuclei for RNA-seq to capture the transcriptome of postmortem neurons. Nat Protoc. 2016;11:499–524.
Chen P, Wang Y, Li J, Bo X, Wang J, Nan L, et al. Diversity and intratumoral heterogeneity in human gallbladder cancer progression revealed by single-cell RNA sequencing. Clin Transl Med. 2021;11(6): e462.
Aizarani N, Saviano A, Mailly L, Durand S, Herman JS, Pessaux P, et al. A human liver cell atlas reveals heterogeneity and epithelial progenitors. Nature. 2019;572(7768):199–204.
Cang Z, Nie Q. Inferring spatial and signaling relationships between cells from single cell transcriptomic data. Nat Commun. 2020;11(1):2084.
Koenitzer JR, Wu H, Atkinson JJ, Brody SL, Humphreys BD. Single-nucleus RNA-sequencing profiling of mouse lung. reduced dissociation bias and improved rare cell-type detection compared with single-cell RNA sequencing. Am J Respir Cell Mol Biol. 2020;63:739–47.
Denisenko E, Guo BB, Jones M, Hou R, Kock LD, Lassmann T, Poppe D, et al. Systematic assessment of tissue dissociation and storage biases in single-cell and single-nucleus RNA-seq workflows. Genome Biol. 2020;21:130.
Lake BB, Chen S, Hoshi M, Plongthongkum N, Salamon D, Knoten A, et al. A single-nucleus RNA-sequencing pipeline to decipher the molecular anatomy and pathophysiology of human kidneys. Nat Commun. 2019;10:2832.
Chen T, Oh S, Gregory S, Shen X, Diehl AM. Single-cell omics analysis reveals functional diversification of hepatocytes during liver regeneration. JCI Insight. 2020;5(22): e141024.
Zhang L, Zhu J, Wang H, Xia J, Liu P, Chen F, et al. A high-resolution cell atlas of the domestic pig lung and an online platform for exploring lung single-cell data. J Genet Genomics. 2021;48(5):411–25.
Zhu J, Chen F, Luo L, Wu W, Dai J, Zhong J, et al. Single-cell atlas of domestic pig cerebral cortex and hypothalamus. Sci Bull (Beijing). 2021;66(14):1448–61.
Herrera UJ, Wiarda JE, Sivasankaran SK, Daharsh L, Liu H, Byrne KA, et al. Reference Transcriptomes of Porcine Peripheral Immune Cells Created Through Bulk and Single-Cell RNA Sequencing. Front Genet. 2021;12: 689406.
Zheng G, Terry J, Belgrader P, Ryvkin P, Bent Z, Wilson R, et al. Massively parallel digital transcriptional profiling of single cells. Nat Commun. 2017;8:14049.
Slyper M, Porter CBM, Ashenberg O, Waldman J, Drokhlyansky E, Wakiro I, Smillie C, et al. A single-cell and single-nucleus RNA-Seq toolbox for fresh and frozen human tumors. Nat Med. 2020;26:792–802.
Butler A, Hoffman P, Smibert P, Papalexi E, Satija R. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat Biotechnol. 2018;36(5):411–20.
Korsunsky I, Millard N, Fan J, Slowikowski K, Zhang F, Wei K, et al. Fast, sensitive and accurate integration of single-cell data with Harmony. Nat Methods. 2019;16(12):1289–96.
Rotta R, Noack A. Multilevel Local Search Algorithms for Modularity Clustering. J. Exp. Alg orithmics. 2011;16, 2.3.
Camp JG, Sekine K, Gerber T, Loeffler WH, Binder H, Gac M, et al. Multilineage communication regulates human liver bud development from pluripotency. Nature. 2017;546(7659):533–8.
Kanehisa M, Goto S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000;28(1):27–30.
Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.
Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51(D1):D587–92.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Gen. 2000;25(1):25–9.
GraphPad Software, Inc. Prism (version 8.4.0). 2018. San Diego, CA: GraphPad Software, Inc.
Diamanti K, Inda Díaz JS, Raine A, Pan G, Wadelius C, Cavalli M. Single nucleus transcriptomics data integration recapitulates the major cell types in human liver. Hepatol Res. 2021;51(2):233–8.
Xiong X, Kuang H, Ansari S, Liu T, Gong J, Wang S, et al. Landscape of Intercellular Crosstalk in Healthy and NASH Liver Revealed by Single-Cell Secretome Gene Analysis. Mol Cell. 2019;75(3):644-660.e5.
Castillo CP, Wanek AG, Kumar P, An X, Elsegeiny W, Horne W, et al. Intestinal IL-17R Signaling Constrains IL-18-Driven Liver Inflammation by the Regulation of Microbiome-Derived Products. Cell Rep. 2019;29(8):2270-2283.e7.
Nault R, Fader KA, Bhattacharya S, Zacharewski TR. Single-Nuclei RNA Sequencing Assessment of the Hepatic Effects of 2,3,7,8-Tetrachlorodibenzo-p-dioxin. Cell Mol Gastroenterol Hepatol. 2021;11(1):147–59.
Hou X, Yang Y, Li P, Zeng Z, Hu W, Zhe R, et al. Integrating Spatial Transcriptomics and Single-Cell RNA-seq Reveals the Gene Expression Profling of the Human Embryonic Liver. Front Cell Dev Biol. 2021;9: 652408.
Richter ML, Deligiannis IK, Yin K, Danese A, Lleshi E, Coupland P, Vallejos CA. Single-nucleus RNA-seq2 reveals functional crosstalk between liver zonation and ploidy. Nat Commun. 2021;12(1):4264.
Grün D, van Oudenaarden A. Design and analysis of single-cell sequencing experiments. Cell. 2015;163(4):799–810.
Ramachandran P, Dobie R, Wilson JR, Dora EF, Henderson BEP, Luu NT, et al. Resolving the fibrotic niche of human liver cirrhosis at single-cell level. Nature. 2019;575(7783):512–8.
Cavalli M, Diamanti K, Pan G, Spalinskas R, Kumar C, Deshmukh AS, et al. A Multi-Omics Approach to Liver Diseases: Integration of Single Nuclei Transcriptomics with Proteomics and HiCap Bulk Data in Human Liver. OMICS. 2020;24(4):180–94.
Zheng H, Pomyen Y, Hernandez MO, Li CY, Livak F, Tang W, et al. Single-cell analysis reveals cancer stem cell heterogeneity in hepatocellular carcinoma. Hepatology. 2018;68(1):127–40.
Tandukar S, Hariharan S. Xenotransplantation. Organogenesis. 2018;14(7163):159–62.
Miller ER, Ullrey DE. The pig as a model for human nutrition. Annu Rev Nutr. 1987;7:361–82.
Yang YG, Sykes M. Xenotransplantation: current status and a perspective on the future. Nat Rev Immunol. 2007;7(7):519–31.
Vodicka P, Smetana KJ, Dvoránková B, Emerick T, Xu YZ, Ourednik J, Ourednik V, Motlík J. The miniature pig as an animal model in biomedical research. Ann N Y Acad Sci. 2005;1049:161–71.
Liang Y, Kaneko K, Xin B, Lee J, Sun X, Zhang K, et al. Temporal analyses of postnatal liver development and maturation by single-cell transcriptomics. Dev Cell. 2022;57(3):398-414.e5.
Du YQ, Sun C, Huang SM, Yang DL, Wu J. Liver sinusoidal endothelial cells regulate adaptive immune tolerance in the liver. Zhonghua Gan Zang Bing Za Zhi. 2017;25(4):309–12.
Hellerbrand C. Hepatic stellate cells–the pericytes in the liver. Pflugers Arch. 2013;465(6):775–8.
Gilboa D, HO Y, DU N, BC N, HB S, Reuveni D, et al. Erythropoietin enhances Kupffer cell number and activity in the challenged liver. Sci Rep. 2017;7(1):10379.
Alcaraz MJ, Ruiz LN, Sebastián RS, Majado MJ, González GC, Bernardo MV, Alvarez MR, Parrado A. Human and mouse DOCK10 splicing isoforms with alternative first coding exon usage are differentially expressed in T and B lymphocytes. Hum Immunol. 2011;72(7):531–7.
Smith X, Taylor A, Rudd CE. T-cell immune adaptor SKAP1 regulates the induction of collagen-induced arthritis in mice. Immunol Lett. 2016;176:122–7.
Goldman O, Han S, Hamou W, JV V, Uzan G, Lickert H, et al. Endoderm generates endothelial cells during liver development. Stem Cell Reports. 2014;3(4):556–65.
Halpern KB, Shenhav R, Massalha H, Toth B, Egozi A, Massasa EE, et al. Paired-cell sequencing enables spatial gene expression mapping of liver endothelial cells. Nat Biotechnol. 2018;36(10):962–70.
MacParland SA, Liu JC, Ma XZ, Innes BT, Bartczak AM, Gage BK, et al. Single cell RNA sequencing of human liver reveals distinct intrahepatic macrophage populations. Nat Commun. 2018;9(1):4383.
Jiang Y, Zhu L, Li X, Shi T. Evaluation of the Chinese indigenous pig breed Dahe and crossbred Dawu for growth and carcass characteristics, organ weight, meat quality and intramuscular fatty acid and amino acid composition. Animal. 2021;5(9):1485–92.
Shi Y, Li X, Huang A. A metabolomics-based approach investigates volatile flavor formation and characteristic compounds of the Dahe black pig dry-cured ham. Meat Sci. 2019;158: 107904.
Zhang X, Zhao SM, Ge CR, Lai H, Li CQ, Gao SZ. Effects of dietary energy levels on growth performance and carcass composition of Wujin pigs. CJAN. 2008;20(5):489–500.
Liu LY, Wang J, Huang Y, Pan HB, Zhang X, Huang ZX, et al. The effect of dietary protein levels on the expression of genes coding for four selected protein translation initiation factors in muscle tissue of Wujin pig. J Anim Physiol Anim Nutr (Berl). 2014;98(2):310–7.
Smedsrod B, Pertoft H, Gustafson S, Laurentet TC. Scavenger functions of the liver endothelial cell. Biochem J. 1990;266(2):313–27.
Harris EN, Kyosseva SV, Weigel JA, Weigel PH. Expression, processing, and glycosaminoglycan binding activity of the recombinant human 315-kDa hyaluronic acid receptor for endocytosis (HARE). J Biol Chem. 2007;282(5):2785–97.
Wang J, Hu W, Shen Z, Liu T, Dai W, Shen B, et al. Dissecting the single-cell transcriptomeunderlying chronic liver injury. Mol Ther Nucleic Acids. 2021;26:1364–73.
Shetty S, Lalor PF, Adams DH. Liver sinusoidal endothelial cells - gatekeepers of hepatic immunity. Nat Rev Gastroenterol Hepatol. 2018;15(9):555–67.
Chaudhry S, Emond J, Griesemer A. Immune Cell Trafficking to the Liver. Transplantation. 2019;103(7):1323–37.
Jiang YJ, Gao JF, Lin LH, Li H, Meng QG, Qu YF, et al. Single-cell transcriptomes from turtle livers reveal sensitivity of hepatic immune cells to bacteria-infection. Fish Shellfish Immunol. 2022;131:847–54.
Willem VDB, Simone R, Jerome S, Curtis C, Auda AE, Fabio L. Single-Cell Transcriptome Analysis of T Cells. Methods Mol Biol. 2019;2048:155–205.
Yu Y, Smoligovets AA, Groves JT. Modulation of T cell signaling by the actin cytoskeleton. J Cell Sci. 2013;126(Pt 5):1049–58.
Rychkov GY, Barritt GJ. Expression and function of TRP channels in liver cells. Adv Exp Med Biol. 2011;704:667–86.
Chang N, Tian L, Ji XF, Zhou X, Hou L, Zhao XH, et al. Single-Cell Transcriptomes Reveal Characteristic Features of Mouse Hepatocytes with Liver Cholestatic Injury. Cells. 2019;8(9):1069.
Gissen P, Arias IM. Structural and functional hepatocyte polarity and liver disease. J Hepatol. 2015;63(4):1023–37.
Yu L, Tai L, Gao J, Sun M, Liu SQ, Huang TD, et al. A New lncRNA, lnc-LLMA, Regulates Lipid Metabolism in Pig Hepatocytes. DNA Cell Biol. 2022;41(2):202–14.
Acknowledgements
Thanks for the sequencing platform and bioinformation analysis of Gene Denovo Biotechnology Co., Ltd (Guangzhou, China).
Funding
This work was supported by the Major Science and Technology Project of Yunnan Province (202202AE090032 and 202102AA310054), National Natural Science Foundation of China (32360808, 31760645, 31260592, 31060331), Technological Innovation Talent Program (2020FA011), State School Cooperation (2020ZXND02), Scientific Research Fund of Yunnan Provincial Department of Education (2023Y0975).
Author information
Authors and Affiliations
Contributions
Conceptualization, S.M.Z., J.H.Z., H.X. and W.Z.L.; methodology, S.M.Z., W.Z.L., J.H.Z., and X.C.G.; software, J.H.Z.; validation, S.M.Z., L.L.Y., and Q.Y.L.; formal analysis, J.H.Z.; investigation, X.C.G.; resources, S.M.Z., H.Y.Z., H.J.W., and H.X.; data curation, H.Y.Z., H.J.W., W.J.C, and H.X.; writing—original draft preparation, J.H.Z.; writing—review and editing, J.H.Z., L.L.Y., and S.M.Z.; visualization, S.M.Z.; supervision, H.Y.Z., H.J.W., W.Z.L., and Y.X.X.; project administration, S.M.Z.; funding acquisition, H.J.W., and S.M.Z.. All authors have read and agreed to the published version of the manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
The experiments were performed according to the ARRIVE guidelines (https://arriveguidelines.org) and approved by the Ethics Committee of Experimental Animal of Yunnan Agricultural University (Approval Code: 202105017, Approval Date: 15 May, 2021).
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1. Table S1.
Gene number in each cluster.
Additional file 2. Table S2.
Cell number in each cluster.
Additional file 3. Table S3.
The top five expression genes from each cluster.
Additional file 4. Table S4.
The marker genes of each cell type.
Additional file 5. Table S5.
The new marker genes of each cell type.
Additional file 6. Table S6.
Cell number in each cell type.
Additional file 7. Table S7.
Number of differential genes per cell type.
Additional file 8. Table S8.
Endothelial cells differential gene GO enrichment.
Additional file 9. Table S9.
Endothelial cells differential gene KEGG significant enrichment(DHL-VS-DHBL).
Additional file 10. Table S10.
T cells differential gene GO enrichment.
Additional file 11. Table S11.
T cells differential gene KEGG significant enrichment(DHL-VS-DHBL).
Additional file 12. Table S12.
Hepatocytes differential gene GO enrichment.
Additional file 13. Table S13.
Hepatocytes differential gene KEGG significant enrichment(DHL-VS-DHBL).
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Zhu, JH., Guan, XC., Yi, LL. et al. Single-nucleus transcriptome sequencing reveals hepatic cell atlas in pigs. BMC Genomics 24, 770 (2023). https://doi.org/10.1186/s12864-023-09765-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12864-023-09765-9