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

Cell Type Resolved Expression of Duplicate Genes Retained From Whole Genome Duplication in Atlantic salmon

The functional and evolutionary outcomes of whole genome duplication (WGD) events are driven by global remodeling of gene expression. Most investigations of gene expression changes following WGD have applied bulk transcriptomics using tissue samples, thus failing to resolve affected cell types. Here, we leverage single-cell transcriptomics of liver tissue in Atlantic salmon (Salmo salar L.) to quantify cell-specific expression and transcriptional responses to a bacterial infection with Aeromonas salmonicida for thousands of duplicate gene pairs (ohnologs) retained from WGD ancestral to all salmonids. The major liver cell types showed hundreds of differentially expressed ohnolog pairs, with hepatocytes showing the greatest number and immune cells the least number of uniquely differentially expressed pairs. Many more differentially expressed ohnolog pairs were identified after accounting for cell type heterogeneity within a cell lineage, despite a reduction in statistical power. The degree of conservation in ohnolog expression responses to bacterial infection also varied significantly among cell types, both in terms of the number of differentially expressed pairs and the direction of responses. Overall, this study highlights the importance of resolving cell-specific gene expression to understand the functional and evolutionary outcomes of WGD events.

Transcriptomic characterization of transitioning cell types in the skin of Atlantic salmon

BACKGROUND

The skin maintains the body's integrity and serves as the first line of defence against pathogens, stressors and mechanical injuries. Despite the global significance of salmon in aquaculture, how the transcriptomic profile of cells varies during wound healing remains unexplored. Teleost's skin contains adult pluripotent cells that differentiate into various tissues, including bone, cartilage, tendon, ligament, adipose, dermis, muscle and connective tissue within the skin. These cells are pivotal for preserving the integrity of skin tissue throughout an organism's lifespan and actively participate in the wound healing processes. In this study, we characterize the transcriptomic profiles of putative mesenchymal stromal cells (fibroblast-like adult stem cells) in healthy Atlantic salmon tissue and during the wound healing process.

RESULTS

Single-nucleus sequencing and spatial transcriptomics were used to detect transcriptomic changes occurring during wound healing that are commonly associated with mesenchymal stromal cells. We followed the transcriptomic activity of these cells during an in vivo wound healing time course study showing that these cells become more transcriptionally active during the remodelling stage of wound healing. The changes detected give insights into the potential differentiation pathways leading to osteogenic and fibroblast lineages in the skin of Atlantic salmon.

CONCLUSIONS

We chart the transcriptomic activity of subclusters of putative differentiating stromal cells during the process of wound healing for the first time, revealing different spatial niches of the various putative MSC subclusters, and setting the stage for further investigation of the manipulation of transitioning cell types to improve fish health.

Single-cell transcriptomic dynamics of scallop heart reveals the heterogeneous response to heat stress

BACKGROUND

Animals with open circulatory systems are highly vulnerable to environmental temperature fluctuations, making them particularly threatened by global warming. However, research on the cellular heterogeneity of heart responses to elevated temperatures in animals with open circulatory systems remains limited.

RESULTS

Here, we conducted a comprehensive investigation of the morphology, metabolism and scRNA-seq of the heart in a molluscan model, Argopecten irradians, under heat stress. Our results unraveled that the severity of cardiac structure damage increased progressively with rising temperature, accompanied by widespread mitochondrial dysfunction and neurohumoral response. We identified two subpopulations within cardiomyocytes (CMs), including ventricular myocytes (VMs) and atrial myocytes (AMs), which exhibited specialized functional roles in response to thermal stress. Specifically, AMs enhanced cell-cell communications with the immune-like cells and fibroblasts to contribute to maintaining cardiac homeostasis under heat stress. Whereas, VMs displayed enhanced energy supply and differentiation potential to withstand thermal challenges. Furthermore, RNA interference targeting the most heat-responsive gene, PLRP2-like, resulted in a significant reduction in heat tolerance and triglyceride accumulation in scallops.

CONCLUSIONS

Our study investigated the heterogeneous response of the scallop heart to high temperatures, revealing distinct response patterns between VMs and AMs. We further identified a key gene, AiPLRP2-like, which exhibits unique cellular localization patterns compared to its mammalian counterpart and may play a pivotal role in regulating cardiac thermotolerance in organisms with open circulatory systems. These findings provide novel insights into the theoretical framework and evolutionary adaptations of marine invertebrate hearts in response to environmental temperature fluctuations.

Advancing genetic improvement in the omics era: status and priorities for United States aquaculture

BACKGROUND

The innovations of the "Omics Era" have ushered in significant advancements in genetic improvement of agriculturally important animal species through transforming genetics, genomics and breeding strategies. These advancements were often coordinated, in part, by support provided over 30 years through the 1993-2023 National Research Support Project 8 (NRSP8, National Animal Genome Research Program, NAGRP) and affiliate projects focused on enabling genomic discoveries in livestock, poultry, and aquaculture species. These significant and parallel advances demand strategic planning of future research priorities. This paper, as an output from the May 2023 Aquaculture Genomics, Genetics, and Breeding Workshop, provides an updated status of genomic resources for United States aquaculture species, highlighting major achievements and emerging priorities.

MAIN TEXT

Finfish and shellfish genome and omics resources enhance our understanding of genetic architecture and heritability of performance and production traits. The 2023 Workshop identified present aims for aquaculture genomics/omics research to build on this progress: (1) advancing reference genome assembly quality; (2) integrating multi-omics data to enhance analysis of production and performance traits; (3) developing resources for the collection and integration of phenomics data; (4) creating pathways for applying and integrating genomics information across animal industries; and (5) providing training, extension, and outreach to support the application of genome to phenome. Research focuses should emphasize phenomics data collection, artificial intelligence, identifying causative relationships between genotypes and phenotypes, establishing pathways to apply genomic information and tools across aquaculture industries, and an expansion of training programs for the next-generation workforce to facilitate integration of genomic sciences into aquaculture operations to enhance productivity, competitiveness, and sustainability.

CONCLUSION

This collective vision of applying genomics to aquaculture breeding with focus on the highlighted priorities is intended to facilitate the continued advancement of the United States aquaculture genomics, genetics and breeding research community and industries. Critical challenges ahead include the practical application of genomic tools and analytical frameworks beyond academic and research communities that require collaborative partnerships between academia, government, and industry. The scope of this review encompasses the use of omics tools and applications in the study of aquatic animals cultivated for human consumption in aquaculture settings throughout their life-cycle.

Mechanisms of immune responses in Acanthopagrus latus to Streptococcus agalactiae infection revealed by transcriptomic analysis

Acanthopagrus latus (yellowfin seabream) is an economically important fish in the southeast coastal sea of China. Its slower growth rate makes it more prone to diseases in the cultivation period, leading to substantial economic losses. Epidemiological investigations have indicated that Streptococcus agalactiae is one of the most common Gram-positive pathogens, which have garnered increasing attention due to its high contagion and lethality rates in A. latus. In this work, an infection model of yellowfin seabream was established with an intraperitoneal injection of S. agalactiae. Clinical sign observations and various analyses, including histological examination, serum biochemical index assessment, immune-related enzyme level measurement, and transcriptome analysis of tissues (liver and intestine) with obvious clinical signs, were conducted for revealing the effects of S. agalactiae infection and immune response mechanisms in yellowfin seabream. The results indicate that evident clinical signs and multi-tissue damages with the notable changes in indices and significant increase in immune-related enzyme levels in the serum occurred in infected fish. RNA sequencing analysis identified 1130 differentially expressed genes (DEGs) in the liver and 1218 DEGs in the intestine, which were involved in multiple immune- and metabolism-related pathways via KEGG enrichment analysis. The transcriptomic results were further corroborated by quantitative real-time RT-PCR (qRT-PCR) tests of some specific immune-related genes. These findings provide new insights into the molecular immune mechanisms in yellowfin seabream following S. agalactiae infection and offer valuable reference data for disease prevention and molecular breeding (i.e., selective breeding through developing molecular markers of key genes).

Exploring the Unknown: How Can We Improve Single-cell RNAseq Cell Type Annotations in Non-model Organisms?

Single-cell RNA sequencing (scRNAseq) is a powerful tool to describe cell types in multicellular organisms across the animal kingdom. In standard scRNAseq analysis pipelines, clusters of cells with similar transcriptional signatures are given cell type labels based on marker genes that infer specialized known characteristics. Since these analyses are designed for model organisms, such as humans and mice, problems arise when attempting to label cell types of distantly related, non-model species that have unique or divergent cell types. Consequently, this leads to limited discovery of novel species-specific cell types and potential mis-annotation of cell types in non-model species while using scRNAseq. To address this problem, we discuss recently published approaches that help annotate scRNAseq clusters for any non-model organism. We first suggest that annotating with an evolutionary context of cell lineages will aid in the discovery of novel cell types and provide a marker-free approach to compare cell types across distantly related species. Secondly, machine learning has greatly improved bioinformatic analyses, so we highlight some open-source programs that use reference-free approaches to annotate cell clusters. Lastly, we propose the use of unannotated genes as potential cell markers for non-model organisms, as many do not have fully annotated genomes and these data are often disregarded. Improving single-cell annotations will aid the discovery of novel cell types and enhance our understanding of non-model organisms at a cellular level. By unifying approaches to annotate cell types in non-model organisms, we can increase the confidence of cell annotation label transfer and the flexibility to discover novel cell types.

Methionine restriction reduced growth performance and exacerbated lipid deposition in hybrid grouper (Epinephelus fuscoguttatus ♀ × E. lanceolatus ♂) under high-lipid diets by suppressing the lipid degradation pathways with the single-cell RNA-seq analysis

The experiment mainly focused on the liver of the hybrid grouper (Epinephelus fuscoguttatus ♀ × E. lanceolatus ♂) with low methionine levels by single-cell RNA-seq under high-lipid diets. Both weight gain rate (WGR) and specific growth rate (SGR) in the MR group were obviously lower than the C group, and the intraperitoneal fat (IPF) in the MR group was obviously enhanced than HL and C groups, which led to more visible growth inhibition and lipid deposition. By using the scRNA-seq analysis, nine cell types were classified as: liver parenchymal cells, erythrocytes, hepatic stellate cells, cholangiocytes, macrophages, T cells, epidermal cells, eosinophil and fibroblasts. Lipid and carbohydrate metabolisms were mainly enriched in liver parenchymal cells. The lipid degradation pathways were obviously inhibited in the MR group, such as cholesterol degradation pathway. The lipid degradation genes were obviously decreased in the MR group. In conclusion, methionine restriction might suppress the lipid degradation pathways in liver parenchymal cells, led to obvious lipid deposition.

Catching the Big Fish in Big Data: A Meta-Analysis of Zebrafish Kidney scRNA-Seq Datasets Highlights Conserved Molecular Profiles of Macrophages and Neutrophils in Vertebrates

The innate immune system (IIS) is an ancient and essential defense mechanism that protects animals against a wide range of pathogens and diseases. Although extensively studied in mammals, our understanding of the IIS in other taxa remains limited. The zebrafish (Danio rerio) serves as a promising model organism for investigating IIS-related processes, yet the immunogenetics of fish are not fully elucidated. To address this gap, we conducted a meta-analysis of single-cell RNA sequencing (scRNA-seq) datasets from zebrafish kidney marrow, encompassing approximately 250,000 immune cells. Our analysis confirms the presence of key genetic pathways in zebrafish innate immune cells that are similar to those identified in mammals. Zebrafish macrophages specifically express genes encoding cathepsins, major histocompatibility complex class II proteins, integral membrane proteins, and the V-ATPase complex and demonstrate the enrichment of oxidative phosphorylation ferroptosis processes. Neutrophils are characterized by the significant expression of genes encoding actins, cytoskeleton organizing proteins, the Arp2/3 complex, and glycolysis enzymes and have demonstrated their involvement in GnRH and CLR signaling pathways, adherents, and tight junctions. Both macrophages and neutrophils highly express genes of NOD-like receptors, phagosomes, and lysosome pathways and genes involved in apoptosis. Our findings reinforce the idea about the existence of a wide spectrum of immune cell phenotypes in fish since we found only a small number of cells with clear pro- or anti-inflammatory signatures.

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.

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.

Keratinocytes drive the epithelial hyperplasia key to sea lice resistance in coho salmon

BACKGROUND

Salmonid species have followed markedly divergent evolutionary trajectories in their interactions with sea lice. While sea lice parasitism poses significant economic, environmental, and animal welfare challenges for Atlantic salmon (Salmo salar) aquaculture, coho salmon (Oncorhynchus kisutch) exhibit near-complete resistance to sea lice, achieved through a potent epithelial hyperplasia response leading to rapid louse detachment. The molecular mechanisms underlying these divergent responses to sea lice are unknown.

RESULTS

We characterized the cellular and molecular responses of Atlantic salmon and coho salmon to sea lice using single-nuclei RNA sequencing. Juvenile fish were exposed to copepodid sea lice (Lepeophtheirus salmonis), and lice-attached pelvic fin and skin samples were collected 12 h, 24 h, 36 h, 48 h, and 60 h after exposure, along with control samples. Comparative analysis of control and treatment samples revealed an immune and wound-healing response that was common to both species, but attenuated in Atlantic salmon, potentially reflecting greater sea louse immunomodulation. Our results revealed unique but complementary roles of three layers of keratinocytes in the epithelial hyperplasia response leading to rapid sea lice rejection in coho salmon. Our results suggest that basal keratinocytes direct the expansion and mobility of intermediate and, especially, superficial keratinocytes, which eventually encapsulate the parasite.

CONCLUSIONS

Our results highlight the key role of keratinocytes in coho salmon's sea lice resistance and the diverged biological response of the two salmonid host species when interacting with this parasite. This study has identified key pathways and candidate genes that could be manipulated using various biotechnological solutions to improve Atlantic salmon sea lice resistance.

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.

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.