Single-cell transcriptomics reveals conserved cell identities and fibrogenic phenotypes in zebrafish and human liver

Leveraging zebrafish to investigate pancreatic development, regeneration, and diabetes

The zebrafish has become an outstanding model for studying organ development and tissue regeneration, which is prominently leveraged for studies of pancreatic development, insulin-producing β-cells, and diabetes. Although studied for more than two decades, many aspects remain elusive and it has only recently been possible to investigate these due to technical advances in transcriptomics, chemical-genetics, genome editing, drug screening, and in vivo imaging. Here, we review recent findings on zebrafish pancreas development, β-cell regeneration, and how zebrafish can be used to provide novel insights into gene functions, disease mechanisms, and therapeutic targets in diabetes, inspiring further use of zebrafish for the development of novel therapies for diabetes.

Single-Cell RNA Sequencing in Organ and Cell Transplantation

Single-cell RNA sequencing is a high-throughput novel method that provides transcriptional profiling of individual cells within biological samples. This method typically uses microfluidics systems to uncover the complex intercellular communication networks and biological pathways buried within highly heterogeneous cell populations in tissues. One important application of this technology sits in the fields of organ and stem cell transplantation, where complications such as graft rejection and other post-transplantation life-threatening issues may occur. In this review, we first focus on research in which single-cell RNA sequencing is used to study the transcriptional profile of transplanted tissues. This technology enables the analysis of the donor and recipient cells and identifies cell types and states associated with transplant complications and pathologies. We also review the use of single-cell RNA sequencing in stem cell implantation. This method enables studying the heterogeneity of normal and pathological stem cells and the heterogeneity in cell populations. With their remarkably rapid pace, the single-cell RNA sequencing methodologies will potentially result in breakthroughs in clinical transplantation in the coming years.

Comprehensive analysis of the miRNA-mRNA regulatory network involved in spontaneous recovery of an H2O2-induced zebrafish cataract model

OBJECTIVE

To identify the hub miRNAs and mRNAs contributing to the spontaneous recovery of an H2O2-induced zebrafish cataract model.

METHODS

Zebrafishes were divided into three groups, i.e., Group A, which included normal control fish (day 0), and Groups B and C, where fish were injected with 2.5% hydrogen peroxide into the anterior chamber and reared for 14 and 30 days, respectively. Fish eyes were examined by stereomicroscope photography and optical coherence tomography (OCT). RNA profiles of fish lenses were detected by RNA sequencing. Differentially expressed genes (DEGs) and differentially expressed miRNAs (DEmiRs) were identified among three groups. The DEGs and DEmiRs, which changed in opposite positions between "B vs. A" and "C vs. B" were defined as ODGs (opposite positions changed DEGs) and ODmiRs (opposite positions changed DEmiRs). Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) analysis were carried out by R language. The protein-protein interaction network (PPI) was constructed using STRING. Potential targets of miRNAs were obtained using miRanda. miRNA-mRNA networks were constructed by Cytoscape.

RESULTS

The fish lens opacity formed on day 14 and recovered to transparent on day 30 after injection. Compared to group B, 1366 DEGs and 54 DEmiRs were identified in group C. "C vs. B" DEGs were enriched in gene clusters related to development and oxidative phosphorylation. Target genes of DEmiRs were enriched in clusters such as development and cysteine metabolism. Among three groups, 786 ODGs and 27 ODmiRs were identified, and 480 ODGs were predicted as targets of ODmiRs. Target ODGs were enriched in pathways related to methionine metabolism, ubiquitin, sensory system development, and structural constituents of the eye lens. In addition, we established an ODmiRs-ODGs regulation network.

CONCLUSION

We identified several hub mRNAs and altered miRNAs in the formation and reversal of zebrafish cataracts. These hub miRNAs/mRNAs could be potential targets for the non-surgical treatment of ARC.

Comparative Analysis of Acute Kidney Injury Models and Related Fibrogenic Responses: Convergence on Methylation Patterns Regulated by Cold Shock Protein

BACKGROUND

Fibrosis is characterized by excessive extracellular matrix formation in solid organs, disrupting tissue architecture and function. The Y-box binding protein-1 (YB-1) regulates fibrosis-related genes (e.g., Col1a1, Mmp2, and Tgfβ1) and contributes significantly to disease progression. This study aims to identify fibrogenic signatures and the underlying signaling pathways modulated by YB-1.

METHODS

Transcriptomic changes associated with matrix gene patterns in human chronic kidney diseases and murine acute injury models were analyzed with a focus on known YB-1 targets. Ybx1-knockout mouse strains (Ybx1ΔRosaERT+TX and Ybx1ΔLysM) were subjected to various kidney injury models. Fibrosis patterns were characterized by histopathological staining, transcriptome analysis, qRT-PCR, methylation analysis, zymography, and Western blotting.

RESULTS

Integrative transcriptomic analyses revealed that YB-1 is involved in several fibrogenic signatures related to the matrisome, the WNT, YAP/TAZ, and TGFß pathways, and regulates Klotho expression. Changes in the methylation status of the Klotho promoter by specific methyltransferases (DNMT) are linked to YB-1 expression, extending to other fibrogenic genes. Notably, kidney-resident cells play a significant role in YB-1-modulated fibrogenic signaling, whereas infiltrating myeloid immune cells have a minimal impact.

CONCLUSIONS

YB-1 emerges as a master regulator of fibrogenesis, guiding DNMT1 to fibrosis-related genes. This highlights YB-1 as a potential target for epigenetic therapies interfering in this process.

Enantioselective Toxicity and Potential Endocrine-Disruptive Effects of the Insecticides Flufiprole and Ethiprole on Danio rerio

Phenylpyrazole insecticides are widely used as chiral pesticides. However, the enantioselective toxicity and potential endocrine-disrupting effects of these insecticides on aquatic organisms remain unclear. Herein, the enantioselective toxicity and potential endocrine-disrupting effects of flufiprole and ethiprole were investigated by using zebrafish embryos/larvae as a model. The acute toxicity of R-flufiprole and R-ethiprole toward zebrafish embryos and larvae was 1.8-3.1-fold higher than that of the S-configuration. Additionally, R-flufiprole and R-ethiprole had a greater effect on the expression of genes related to the hypothalamus-pituitary-gonad axis in zebrafish compared with the S-configuration. Nevertheless, both S-flufiprole and S-ethiprole exhibited a greater interference effect on the expression of genes related to the hypothalamus-pituitary-thyroid axis and a greater teratogenic effect on zebrafish than the R-configuration. Thus, this study demonstrates that both flufiprole and ethiprole exhibit enantioselective acute toxicity and developmental toxicity toward zebrafish. Furthermore, those pesticides potentially possess enantioselective endocrine-disrupting effects.

Fishing Innate Immune System Properties through the Transcriptomic Single-Cell Data of Teleostei

The innate immune system is the first line of defense in multicellular organisms. Danio rerio is widely considered a promising model for IIS-related research, with the most amount of scRNAseq data available among Teleostei. We summarized the scRNAseq and spatial transcriptomics experiments related to the IIS for zebrafish and other Teleostei from the GEO NCBI and the Single-Cell Expression Atlas. We found a considerable number of scRNAseq experiments at different stages of zebrafish development in organs such as the kidney, liver, stomach, heart, and brain. These datasets could be further used to conduct large-scale meta-analyses and to compare the IIS of zebrafish with the mammalian one. However, only a small number of scRNAseq datasets are available for other fish (turbot, salmon, cavefish, and dark sleeper). Since fish biology is very diverse, it would be a major mistake to use zebrafish alone in fish immunology studies. In particular, there is a special need for new scRNAseq experiments involving nonmodel Teleostei, e.g., long-lived species, cancer-resistant fish, and various fish ecotypes.

The knockout of cytoglobin 1 in zebrafish (Danio rerio) alters lipid metabolism, iron homeostasis and oxidative stress response

Cytoglobin (Cygb) is an evolutionary ancient heme protein with yet unclear physiological function(s). Mammalian Cygb is ubiquitously expressed in all tissues and is proposed to be involved in reactive oxygen species (ROS) detoxification, nitric oxide (NO) metabolism and lipid-based signaling processes. Loss-of-function studies in mouse associate Cygb with apoptosis, inflammation, fibrosis, cardiovascular dysfunction or oncogenesis. In zebrafish (Danio rerio), two cygb genes exist, cytoglobin 1 (cygb1) and cytoglobin 2 (cygb2). Both have different coordination states and distinct expression sites within zebrafish tissues. The biological roles of the cygb paralogs are largely uncharacterized. We used a CRISPR/Cas9 genome editing approach and generated a knockout of the penta-coordinated cygb1 for in vivo analysis. Adult male cygb1 knockouts develop phenotypic abnormalities, including weight loss. To identify the molecular mechanisms underlying the occurrence of these phenotypes and differentiate between function and effect of the knockout we compared the transcriptomes of cygb1 knockout at different ages to age-matched wild-type zebrafish. We found that immune regulatory and cell cycle regulatory transcripts (e.g. tp53) were up-regulated in the cygb1 knockout liver. Additionally, the expression of transcripts involved in lipid metabolism and transport, the antioxidative defense and iron homeostasis was affected in the cygb1 knockout. Cygb1 may function as an anti-inflammatory and cytoprotective factor in zebrafish liver, and may be involved in lipid-, iron-, and ROS-dependent signaling.

The C-type lectin COLEC10 is predominantly produced by hepatic stellate cells and involved in the pathogenesis of liver fibrosis

Hepatic stellate cell is one of the major nonparenchymal cell types in liver. It has been proved the hepatic stellate cells are activated upon liver injury and produce excessive extracellular matrix to induce liver fibrosis. Single-cell RNA sequencing has been introduced to identify the subpopulations and function of hepatic stellate cells for its remarkable resolution of representation of single-cell transcriptome. According to the re-analysis of single-cell RNA sequencing data and pseudotime trajectory inference, we have found the C-type lectins including Colec10 and Colec11 are not produced by hepatocytes but predominantly produced by hepatic stellate cells, especially quiescent ones in the mice livers. In addition, the expression of Colec10 is decreased in the fibrotic livers of CCl4-challenged mice. COLEC10 is also mainly expressed in the hepatic stellate cells of human livers and the expression of COLEC10 is decreased with the progression of liver fibrosis. The bulk RNA sequencing data of the lentivirus transfected LX-2 cells indicates the function of COLEC10 is associated with inflammation, angiogenesis and extracellular matrix alteration. Surprisingly, the in vitro overexpression of COLEC10 in LX-2 cells promotes the mRNA expression of extracellular matrix components including COL1A1, COL1A2 and COL3A1 and the extracellular matrix degradation enzyme MMP2. To further investigate the role of COLEC10 in the pathogenesis of liver fibrosis, the serum concentration of COLEC10 in patients with chronic liver disease and healthy donors is measured. The serum concentration of COLEC10 is elevated in the patients with chronic liver disease compared to the healthy donors and positively correlated with serum concentration of the D-dimer but not the most of liver function markers. Altogether, we conclude that the C-type lectin COLEC10 is predominantly produced by the hepatic stellate cells and involved in the pathogenesis of liver fibrosis.

Functional genomics in Spiralia

Our understanding of the mechanisms that modulate gene expression in animals is strongly biased by studying a handful of model species that mainly belong to three groups: Insecta, Nematoda and Vertebrata. However, over half of the animal phyla belong to Spiralia, a morphologically and ecologically diverse animal clade with many species of economic and biomedical importance. Therefore, investigating genome regulation in this group is central to uncovering ancestral and derived features in genome functioning in animals, which can also be of significant societal impact. Here, we focus on five aspects of gene expression regulation to review our current knowledge of functional genomics in Spiralia. Although some fields, such as single-cell transcriptomics, are becoming more common, the study of chromatin accessibility, DNA methylation, histone post-translational modifications and genome architecture are still in their infancy. Recent efforts to generate chromosome-scale reference genome assemblies for greater species diversity and optimise state-of-the-art approaches for emerging spiralian research systems will address the existing knowledge gaps in functional genomics in this animal group.

Widespread sex dimorphism across single-cell transcriptomes of adult African turquoise killifish tissues

The African turquoise killifish (Nothobranchius furzeri), the shortest-lived vertebrate that can be bred in captivity, is an emerging model organism for aging research. Here, we describe a multitissue, single-cell gene expression atlas of female and male blood, kidney, liver, and spleen. We annotate 22 cell types, define marker genes, and infer differentiation trajectories. We find pervasive sex-dimorphic gene expression across cell types. Sex-dimorphic genes tend to be linked to lipid metabolism, consistent with clear differences in lipid storage in female vs. male turquoise killifish livers. We use machine learning to predict sex using single-cell gene expression and identify potential markers for molecular sex identity. As a proof of principle, we show that our atlas can be used to deconvolute existing bulk RNA sequencing (RNA-seq) data to obtain accurate estimates of cell type proportions. This atlas can be a resource to the community that could be leveraged to develop cell-type-specific expression in transgenic animals.

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

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

Zebrafish as a Useful Model System for Human Liver Disease

Liver diseases represent a significant global health challenge, thereby necessitating extensive research to understand their intricate complexities and to develop effective treatments. In this context, zebrafish (Danio rerio) have emerged as a valuable model organism for studying various aspects of liver disease. The zebrafish liver has striking similarities to the human liver in terms of structure, function, and regenerative capacity. Researchers have successfully induced liver damage in zebrafish using chemical toxins, genetic manipulation, and other methods, thereby allowing the study of disease mechanisms and the progression of liver disease. Zebrafish embryos or larvae, with their transparency and rapid development, provide a unique opportunity for high-throughput drug screening and the identification of potential therapeutics. This review highlights how research on zebrafish has provided valuable insights into the pathological mechanisms of human liver disease.

Ticlopidine induces embryonic development toxicity and hepatotoxicity in zebrafish by upregulating the oxidative stress signaling pathway

Ticlopidine exerts its anti-platelet effects mainly by antagonizing platelet p2y12 receptors. Previously, a few studies have shown that ticlopidine can induce liver injury, but the exact mechanism of hepatotoxicity remains unclear. Oxidative stress, metabolic disorders, hepatocyte apoptosis, lipid peroxidation, and inflammatory responses can all lead to hepatic liver damage, which can cause hepatotoxicity. In this study, in order to deeply explore the potential molecular mechanisms of ticlopidine -induced hepatotoxicity, we used zebrafish as a model organism to comprehensively evaluate the hepatotoxicity of ticlopidine and its associated mechanism. Three days post-fertilization, zebrafish larvae were exposed to varying concentrations (1.5, 1.75 and 2 μg/mL) of ticlopidine for 72 h, in contrast, adult zebrafish were exposed exposure to 4 μg/mL of ticlopidine for 28 days. Ticlopidine-exposed zebrafish larvae showed changes in liver morphology, shortened body length, and delayed development of the swim bladder development. Liver tissues of ticlopidine-exposed zebrafish larvae and adults stained with Hematoxylin & Eosin revealed vacuolization and increased cellular interstitial spaces in liver tissues. Furthermore, using Oil Red O and periodic acid-Schiff staining methods and evaluating different metabolic enzymes of ticlopidine-exposed zebrafish larvae and adults suggested abnormal liver metabolism and liver injury in both ticlopidine-exposed zebrafish larvae and adults. Ticlopidine also significantly elevated inflammation and oxidative stress and reduced hepatocyte proliferation. During the rescue intervention using N-acetylcysteine, we observed significant improvement in ticlopidine-induced morphological changes in the liver, shortened body length, delayed swim bladder development, and proliferation of liver tissues showed significant improvement. In conclusion, ticlopidine might inhibit normal development and liver proliferation in zebrafish by upregulation of oxidative stress levels, thus leading to embryonic developmental toxicity and hepatotoxicity. In this study, we used zebrafish as a model organism to elucidate the developmental toxicity and hepatotoxicity induced by ticlopidine upregulation of oxidative stress signaling pathway in zebrafish, providing a theoretical basis for clinical application.

Liver in infections: a single-cell and spatial transcriptomics perspective

The liver is an immune organ that plays a vital role in the detection, capture, and clearance of pathogens and foreign antigens that invade the human body. During acute and chronic infections, the liver transforms from a tolerant to an active immune state. The defence mechanism of the liver mainly depends on a complicated network of intrahepatic and translocated immune cells and non-immune cells. Therefore, a comprehensive liver cell atlas in both healthy and diseased states is needed for new therapeutic target development and disease intervention improvement. With the development of high-throughput single-cell technology, we can now decipher heterogeneity, differentiation, and intercellular communication at the single-cell level in sophisticated organs and complicated diseases. In this concise review, we aimed to summarise the advancement of emerging high-throughput single-cell technologies and re-define our understanding of liver function towards infections, including hepatitis B virus, hepatitis C virus, Plasmodium, schistosomiasis, endotoxemia, and corona virus disease 2019 (COVID-19). We also unravel previously unknown pathogenic pathways and disease mechanisms for the development of new therapeutic targets. As high-throughput single-cell technologies mature, their integration into spatial transcriptomics, multiomics, and clinical data analysis will aid in patient stratification and in developing effective treatment plans for patients with or without liver injury due to infectious diseases.

Transgenerational exposure to ocean acidification impacts the hepatic transcriptome of European sea bass (Dicentrarchus labrax)

Physiological effects of ocean acidification associated with elevated CO2 concentrations in seawater is the subject of numerous studies in teleost fish. While the short time within-generation impact of ocean acidification (OA) on acid-base exchange and energy metabolism is relatively well described, the effects associated with transgenerational exposure to OA are much less known. Yet, the impacts of OA can vary in time with the potential for acclimation or adaptation of a species. Previous studies in our lab demonstrated that transgenerational exposure to OA had extensive effects on the transcriptome of the olfactory epithelium of European sea bass (Dicentrarchus labrax), especially on genes related to ion balance, energy metabolism, immune system, synaptic plasticity, neuron excitability and wiring. In the present study, we complete the previous work by investigating the effect of transgenerational exposure to OA on the hepatic transcriptome of European sea bass. Differential gene expression analysis was performed by RNAseq technology on RNA extracted from the liver of two groups of 18 months F2 juveniles that had been exposed since spawning to the same AO conditions as their parents (F1) to either actual pH or end-of-century predicted pH levels (IPCC RCP8.5), respectively. Here we show that transgenerational exposure to OA significantly impacts the expression of 236 hepatic transcripts including genes mainly involved in inflammatory/immune responses but also in carbohydrate metabolism and cellular homeostasis. Even if this transcriptomic impact is relatively limited compared to what was shown in the olfactory system, this work confirmed that fish transgenerationally exposed to OA exhibit molecular regulation of processes related to metabolism and inflammation. Also, our data expand the up-regulation of a key gene involved in different physiological pathways including calcium homeostasis (i.e. pthr1), which we already observed in the olfactory epithelium, to the liver. Even if our experimental design does not allow to discriminate direct within F2 generation effects from transgenerational plasticity, these results offer the perspective of more functional analyses to determine the potential physiological impact of OA exposure on fish physiology with ecological relevance.

Resveratrol attenuated fatty acid synthesis through MAPK-PPAR pathway in red tilapia

High-fat (HF) diets have been shown to cause hepatic impairment in fish species, but the mode of action, especially the pathways involved, has not yet been determined. In this study, the effects of resveratrol (RES) supplementation on the hepatic structure and fat metabolism of red tilapia (Oreochromis niloticus) were determined. Based on transcriptome and proteomics results, RES was found to promote fatty acid β-oxidation in the blood, liver, and liver cells associated with apoptosis and the MAPK/PPAR signaling pathway. RES supplementation was found to alter the expression of genes related to apoptosis and fatty acid pathways like blood itga6a and armc5 which were upregulated and downregulated respectively by high-fat feeding while ggh and ensonig00000008711 increased and decreased, respectively, with RES addition. Relative to the PPAR signaling pathway, fabp10a and acbd7 showed a reverse U-shaped tendency, both in different treatments and at different times. Proteomics results demonstrated that MAPK/PPAR, carbon/glyoxylate, dicarboxylate/glycine serine, and threonine/drug-other enzymes/beta-alanine metabolism pathways in the RES group were significantly affected, and Fasn and Acox1 decreased and increased, respectively, with RES addition. Seven subgroups were obtained using scRNA-seq, and enrichment analysis showed that the PPAR signaling pathway was upregulated with RES supplementation. RES significantly increased the expression of the marked genes (pck1) ensonig00000037711, fbp10a, granulin, hbe1, and zgc:136461, which are liver cell-specific genes. In conclusion, RES resulted in significantly enriched DGEs associated with fat metabolism and synthesis via the MAPK-PPAR signaling pathway.

Widespread sex-dimorphism across single-cell transcriptomes of adult African turquoise killifish tissues

The African turquoise killifish (Nothobranchius furzeri), the shortest-lived vertebrate that can be bred in captivity, is an emerging model organism to study vertebrate aging. Here we describe the first multi-tissue, single-cell gene expression atlas of female and male turquoise killifish tissues comprising immune and metabolic cells from the blood, kidney, liver, and spleen. We were able to annotate 22 distinct cell types, define associated marker genes, and infer differentiation trajectories. Using this dataset, we found pervasive sex-dimorphic gene expression across cell types, especially in the liver. Sex-dimorphic genes tended to be involved in processes related to lipid metabolism, and indeed, we observed clear differences in lipid storage in female vs. male turquoise killifish livers. Importantly, we use machine-learning to predict sex using single-cell gene expression in our atlas and identify potential transcriptional markers for molecular sex identity in this species. As proof-of-principle, we show that our atlas can be used to deconvolute existing liver bulk RNA-seq data in this species to obtain accurate estimates of cell type proportions across biological conditions. We believe that this single-cell atlas can be a resource to the community that could notably be leveraged to identify cell type-specific genes for cell type-specific expression in transgenic animals.

Single-cell analysis of shared signatures and transcriptional diversity during zebrafish development

During development, animals generate distinct cell populations with specific identities, functions, and morphologies. We mapped transcriptionally distinct populations across 489,686 cells from 62 stages during wild-type zebrafish embryogenesis and early larval development (3-120 hours post-fertilization). Using these data, we identified the limited catalog of gene expression programs reused across multiple tissues and their cell-type-specific adaptations. We also determined the duration each transcriptional state is present during development and suggest new long-term cycling populations. Focused analyses of non-skeletal muscle and the endoderm identified transcriptional profiles of understudied cell types and subpopulations, including the pneumatic duct, individual intestinal smooth muscle layers, spatially distinct pericyte subpopulations, and homologs of recently discovered human best4+ enterocytes. The transcriptional regulators of these populations remain unknown, so we reconstructed gene expression trajectories to suggest candidates. To enable additional discoveries, we make this comprehensive transcriptional atlas of early zebrafish development available through our website, Daniocell.

In vivo imaging of calcium dynamics in zebrafish hepatocytes

BACKGROUND AND AIMS

Hepatocytes were the first cell type for which oscillations of cytoplasmic calcium levels in response to hormones were described. Since then, investigation of calcium dynamics in liver explants and culture has greatly increased our understanding of calcium signaling. A bottleneck, however, exists in observing calcium dynamics in a noninvasive manner because of the optical inaccessibility of the mammalian liver. Here, we aimed to take advantage of the transparency of the zebrafish larvae to image hepatocyte calcium dynamics in vivo at cellular resolution.

APPROACH AND RESULTS

We developed a transgenic model expressing a calcium sensor, GCaMP6s, specifically in zebrafish hepatocytes. Using this, we provide a quantitative assessment of intracellular calcium dynamics during multiple contexts, including growth, feeding, ethanol-induced stress, and cell ablation. Specifically, we show that synchronized calcium oscillations are present in vivo , which are lost upon starvation. Starvation induces lipid accumulation in the liver. Feeding recommences calcium waves in the liver, but in a spatially restricted manner, as well as resolves starvation-induced hepatic steatosis. By using a genetically encoded scavenger for calcium, we show that dampening of calcium signaling accelerates the accumulation of starvation-related lipid droplets in the liver. Furthermore, ethanol treatment, as well as cell ablation, induces calcium flux, but with different dynamics. The former causes asynchronous calcium oscillations, whereas the latter leads to a single calcium spike.

CONCLUSIONS

We demonstrate the presence of oscillations, waves, and spikes in vivo . Calcium waves are present in response to nutrition and negatively regulate starvation-induced accumulation of lipid droplets.

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

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

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

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

Single cell transcriptomics of Atlantic salmon (Salmo salar L.) liver reveals cellular heterogeneity and immunological responses to challenge by Aeromonas salmonicida

The liver is a multitasking organ with essential functions for vertebrate health spanning metabolism and immunity. In contrast to mammals, our understanding of liver cellular heterogeneity and its role in regulating immunological status remains poorly defined in fishes. Addressing this knowledge gap, we generated a transcriptomic atlas of 47,432 nuclei isolated from the liver of Atlantic salmon (Salmo salar L.) contrasting control fish with those challenged with a pathogenic strain of Aeromonas salmonicida, a problematic bacterial pathogen in global aquaculture. We identified the major liver cell types and their sub-populations, revealing poor conservation of many hepatic cell marker genes utilized in mammals, while identifying novel heterogeneity within the hepatocyte, lymphoid, and myeloid lineages. This included polyploid hepatocytes, multiple T cell populations including γδ T cells, and candidate populations of monocytes/macrophages and dendritic cells. A dominant hepatocyte population radically remodeled its transcriptome following infection to activate the acute phase response and other defense functions, while repressing routine functions such as metabolism. These defense-specialized hepatocytes showed strong activation of genes controlling protein synthesis and secretion, presumably to support the release of acute phase proteins into circulation. The infection response further involved up-regulation of numerous genes in an immune-cell specific manner, reflecting functions in pathogen recognition and killing, antigen presentation, phagocytosis, regulation of inflammation, B cell differentiation and T cell activation. Overall, this study greatly enhances our understanding of the multifaceted role played by liver immune and non-immune cells in host defense and metabolic remodeling following infection and provides many novel cell-specific marker genes to empower future studies of this organ in fishes.

Cellular and transcriptional heterogeneity in the intrahepatic biliary epithelium

Epithelial tissues comprise heterogeneous cellular subpopulations, which often compartmentalize specialized functions like absorption and secretion to distinct cell types. In the liver, hepatocytes and biliary epithelial cells (BECs; also called cholangiocytes) are the two major epithelial lineages and play distinct roles in (1) metabolism, protein synthesis, detoxification, and (2) bile transport and modification, respectively. Recent technological advances, including single cell transcriptomic assays, have shed new light on well-established heterogeneity among hepatocytes, endothelial cells, and immune cells in the liver. However, a "ground truth" understanding of molecular heterogeneity in BECs has remained elusive, and the field currently lacks a set of consensus biomarkers for identifying BEC subpopulations. Here, we review long-standing definitions of BEC heterogeneity as well as emerging studies that aim to characterize BEC subpopulations using next generation single cell assays. Understanding cellular heterogeneity in the intrahepatic bile ducts holds promise for expanding our foundational mechanistic knowledge of BECs during homeostasis and disease.