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

Single-cell atlas of rainbow trout peripheral blood leukocytes and profiling of their early response to infectious pancreatic necrosis virus

The recent development of single cell sequencing technologies has revolutionized the state-of-art of cell biology, allowing the simultaneous measurement of thousands of genes in single cells. This technology has been applied to study the transcriptome of single cells in homeostasis and also in response to pathogenic exposure, greatly increasing our knowledge of the immune response to infectious agents. Yet the number of these studies performed in aquacultured fish species is still very limited. Thus, in the current study, we have used the 10x Genomics single cell RNA sequencing technology to study the response of rainbow trout (Oncorhynchus mykiss) peripheral blood leukocytes (PBLs) to infectious pancreatic necrosis virus (IPNV), an important trout pathogen. The study allowed us to obtain a transcriptomic profile of 12 transcriptionally distinct leukocyte cell subpopulations that included four different subsets of B cells, T cells, monocytes, two populations of dendritic-like cells (DCs), hematopoietic progenitor cells, non-specific cytotoxic cells (NCC), neutrophils and thrombocytes. The transcriptional pattern of these leukocyte subpopulations was compared in PBL cultures that had been exposed in vitro to IPNV for 24 h and mock-infected cultures. Our results revealed that monocytes and neutrophils showed the highest number of upregulated protein-coding genes in response to IPNV. Interestingly, IgM+IgD+ and IgT+ B cells also upregulated an important number of genes to the virus, but a much fainter response was observed in ccl4 + or plasma-like cells (irf4 + cells). A substantial number of protein-coding genes and genes coding for ribosomal proteins were also transcriptionally upregulated in response to IPNV in T cells and thrombocytes. Interestingly, although genes coding for ribosomal proteins were regulated in all affected PBL subpopulations, the number of such genes transcriptionally regulated was higher in IgM+IgD+ and IgT+ B cells. A further analysis dissected which of the regulated genes were common and which were specific to the different cell clusters, identifying eight genes that were transcriptionally upregulated in all the affected groups. The data provided constitutes a comprehensive transcriptional perspective of how the different leukocyte populations present in blood respond to an early viral encounter in fish.

Single-Cell Transcriptomic Analysis Revealed the Cell Population Changes and Cell-Cell Communication in the Liver of a Carnivorous Fish in Response to High-Carbohydrate Diet

BACKGROUND

Carnivorous fish have a low carbohydrate utilization ability, and the physiologic and molecular basis of glucose intolerance has not been fully illustrated.

OBJECTIVES

This study aimed to use largemouth bass as a model to investigate the possible mechanism of glucose intolerance in carnivorous fish with the help of single-nuclei RNA sequencing (snRNA-seq).

METHODS

Two diets were formulated, a low-carbohydrate (LC) diet and a high-carbohydrate (HC) diet. The feeding trial lasted for 6 wk, and then, growth performance, biochemical parameters, liver histology, and snRNA-seq were performed.

RESULTS

Growth performance of fish was not affected by the HC diet, while liver glucolipid metabolism disorder and liver injury were observed. A total of 13,247 and 12,848 cells from the liver derived from 2 groups were isolated and sequenced, and 7 major liver cell types were annotated by the marker genes. Hepatocytes and cholangiocytes were lower and hepatic stellate cells (HSCs) and immune cells were higher in the HC group than those in the LC group. Reclustering analysis identified 7 subtypes of hepatocytes and immune cells, respectively. The HSCs showed more cell communication with other cell types, and periportal hepatocytes showed more cell communication with other hepatocyte subtypes. Cell-cell communication mainly focused on cell junction-related signaling pathways. Uncovered by the pseudotime analysis, midzonal hepatocytes were differentiated into 2 major branches-biliary epithelial hepatocytes and hepatobiliary hybrid progenitor. Cell junction and liver fibrosis-related genes were highly expressed in the HC group. HC diet induced the activation of HSCs and, therefore, led to the liver fibrosis of largemouth bass.

CONCLUSIONS

HC diet induces liver glucolipid metabolism disorder and liver injury of largemouth bass. The increase and activation of HSCs might be the main reason for the liver injury. In adaption to HC diet, midzonal hepatocytes differentiates into 2 major branches-biliary epithelial hepatocytes and hepatobiliary hybrid progenitors.

Single-nucleus transcriptome reveals cell dynamic response of liver during the late chick embryonic development

The late embryonic development of the liver, a major metabolic organ, remains poorly characterized at single cell resolution. Here, we used single-nucleus RNA-sequencing (snRNA-seq) to characterize the chicken liver cells at 2 embryonic development time points (E14 and D1). We uncovered 8 cell types including hepatocytes, endothelial cells, hepatic stellate cells, erythrocytes, cholangiocytes, kupffer cells, mesothelial cells, and lymphocytes. And we discovered significant differences in the abundance of different cell types between E14 and D1. Moreover, we characterized the heterogeneity of hepatocytes, endothelial cells, and mesenchymal cells based on the gene regulatory networks of each clusters. Trajectory analyses revealed 128 genes associated with hepatocyte development and function, including apolipoprotein genes involved hepatic lipid metabolism and NADH dehydrogenase subunits involved hepatic oxidative phosphorylation. Furthermore, we identified the differentially expressed genes (DEGs) between E14 and D1 at the cellular levels, which contribute to changes in liver development and function. These DEGs were significantly enriched in PPAR signaling pathways and lipid metabolism related pathways. Our results presented the single-cell mapping of chick embryonic liver at late stages of development and demonstrated the metabolic changes across the 2 age stages at the cellular level, which can help to further study the molecular development mechanism of embryonic liver.

Season affects the estrogen system and the immune response of common carp

The physiology of ectothermic animals, including fish, is strictly regulated by season-related external factors such as temperature or photoperiod. The immune response and the production of hormones, such as estrogens, are therefore also subject to seasonal changes. This study in common carp aimed to determine how the season affects the estrogen system and the immune response, including the antibacterial response during Aeromonas salmonicida infection. We compared the immune reaction in spring and autumn in the head kidney and liver and found that carp have higher levels of blood 17β-estradiol in autumn, while in the liver of these fish there is a higher constitutive expression of genes encoding vitellogenin, estrogen receptors and Cyp19 aromatase than in spring. Fish sampled in autumn also exhibited higher expression of immune-related genes in the liver. In contrast, in the head kidney from fish sampled in the autumn, the expression of genes encoding estrogen receptors and aromatase was lower than in spring, and a similar profile of expression was also measured in the head kidney for inos, arginases and il-10. In turn, during bacterial infection, we observed higher upregulation of the expression of inos, il-12p35, ifnγ-2, arginase 2 and il-10 in the liver of carp sampled in spring. In the liver of carp infected in spring a higher upregulation of the expression of the genes encoding CRPs was observed compared to fish infected during autumn. The opposite trend occurred in the head kidney, where the upregulation of the expression of the genes involved in the immune response was higher in fish infected in autumn than in those infected in spring. During the infection, also season-dependent changes occurred in the estrogen system. In conclusion, we demonstrated that season differentially affects the estrogenic and immune activity of the head kidney and liver. These results reinforce our previous findings that the endocrine and immune systems cooperate in maintaining homeostasis and fighting infection.

Mapping the cellular landscape of Atlantic salmon head kidney by single cell and single nucleus transcriptomics

Single-cell transcriptomics is the current gold standard for global gene expression profiling, not only in mammals and model species, but also in non-model fish species. This is a rapidly expanding field, creating a deeper understanding of tissue heterogeneity and the distinct functions of individual cells, making it possible to explore the complexities of immunology and gene expression on a highly resolved level. In this study, we compared two single cell transcriptomic approaches to investigate cellular heterogeneity within the head kidney of healthy farmed Atlantic salmon (Salmo salar). We compared 14,149 cell transcriptomes assayed by single cell RNA-seq (scRNA-seq) with 18,067 nuclei transcriptomes captured by single nucleus RNA-Seq (snRNA-seq). Both approaches detected eight major cell populations in common: granulocytes, heamatopoietic stem cells, erythrocytes, mononuclear phagocytes, thrombocytes, B cells, NK-like cells, and T cells. Four additional cell types, endothelial, epithelial, interrenal, and mesenchymal cells, were detected in the snRNA-seq dataset, but appeared to be lost during preparation of the single cell suspension submitted for scRNA-seq library generation. We identified additional heterogeneity and subpopulations within the B cells, T cells, and endothelial cells, and revealed developmental trajectories of heamatopoietic stem cells into differentiated granulocyte and mononuclear phagocyte populations. Gene expression profiles of B cell subtypes revealed distinct IgM and IgT-skewed resting B cell lineages and provided insights into the regulation of B cell lymphopoiesis. The analysis revealed eleven T cell sub-populations, displaying a level of T cell heterogeneity in salmon head kidney comparable to that observed in mammals, including distinct subsets of cd4/cd8-negative T cells, such as tcrγ positive, progenitor-like, and cytotoxic cells. Although snRNA-seq and scRNA-seq were both useful to resolve cell type-specific expression in the Atlantic salmon head kidney, the snRNA-seq pipeline was overall more robust in identifying several cell types and subpopulations. While scRNA-seq displayed higher levels of ribosomal and mitochondrial genes, snRNA-seq captured more transcription factor genes. However, only scRNA-seq-generated data was useful for cell trajectory inference within the myeloid lineage. In conclusion, this study systematically outlines the relative merits of scRNA-seq and snRNA-seq in Atlantic salmon, enhances understanding of teleost immune cell lineages, and provides a comprehensive list of markers for identifying major cell populations in the head kidney with significant immune relevance.

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.

Single cell genomics as a transformative approach for aquaculture research and innovation

Single cell genomics encompasses a suite of rapidly maturing technologies that measure the molecular profiles of individual cells within target samples. These approaches provide a large up-step in biological information compared to long-established 'bulk' methods that profile the average molecular profiles of all cells in a sample, and have led to transformative advances in understanding of cellular biology, particularly in humans and model organisms. The application of single cell genomics is fast expanding to non-model taxa, including aquaculture species, where numerous research applications are underway with many more envisaged. In this review, we highlight the potential transformative applications of single cell genomics in aquaculture research, considering barriers and potential solutions to the broad uptake of these technologies. Focusing on single cell transcriptomics, we outline considerations for experimental design, including the essential requirement to obtain high quality cells/nuclei for sequencing in ectothermic aquatic species. We further outline data analysis and bioinformatics considerations, tailored to studies with the under-characterized genomes of aquaculture species, where our knowledge of cellular heterogeneity and cell marker genes is immature. Overall, this review offers a useful source of knowledge for researchers aiming to apply single cell genomics to address biological challenges faced by the global aquaculture sector though an improved understanding of cell biology.

Novel insights into the cytokine network of rainbow trout Oncorhynchus mykiss using cell lines and primary leukocyte populations

Cytokines are small proteins that regulate innate and adaptive immune responses and are released by both immune and non-immune cell types. In the current study, the constitutive and induced gene expression profiles of a suite of proinflammatory and regulatory cytokines was examined comparatively in eight rainbow trout (Oncorhynchus mykiss) cell lines, in order to establish the cytokine repertoires of these different cell types, especially the understudied non-immune cells. They included three epithelial cell lines (RTgut, RTgill, and RTL), one endothelial cell line (RTH), one fibroblast cell line (RTG-2), two stromal cell lines (TSS and TPS-2) and one monocyte/macrophage-like cell line (RTS-11). Three types of primary leukocytes (derived from blood, spleen and head kidney) of trout were also included in the analysis, to allow comparison to the repertoires expressed in T cells, as a major source of cytokines in immune responses. The major findings are: 1) IL-2A, IL-2B, IL-4/13B1, IL-4/13B2, IL-10b, P40B1, P28B, IL-17A/F1b, TNF-α3, TNF-α4, IFNγ1, CCL20L2b and CCL20L3a are expressed mainly in leukocytes but IL-17 N, IL-17D, IL-20 and CCL20L1b2 are not expressed in these cells. Hence future studies in these cell lines will help establish their function in fish; 2) Some of the cytokines were differentially expressed in the cell lines, revealing the potential role of these cell types in aspects of trout mucosal and inflammatory immune responses, 3) Similar cell types grouped together in the cell cluster analysis, including the leukocyte cluster, stromal cell cluster, and epithelial and endothelial cell cluster. Taken together, this investigation of these trout cell lines forms a good database for studying the function of cytokines not expressed in isolated leukocytes or that are preferentially expressed in the cell lines. Furthermore, the cytokine expression analysis undertaken confirmed the phenotypic relationship of these cell types at the molecular level.