RNA-Seq of Single Fish Cells - Seeking Out the Leukocytes Mediating Immunity in Teleost Fishes

Integrative RNA, miRNA, and 16S rRNA sequencing reveals immune-related regulation network for glycinin-induced enteritis in hybrid yellow catfish, Pelteobagrus fulvidraco ♀ × Pelteobagrus vachelli ♂

Glycinin-induced foodborne enteritis is a significant obstacle that hinders the healthy development of the aquatic industry. Glycinin causes growth retardation and intestinal damage in hybrid yellow catfish (Pelteobagrus fulvidraco ♀ × Pelteobagrus vachelli ♂), but its immune mechanisms are largely unknown. In the current study, five experimental diets containing 0% (CK), 1.74% (G2), 3.57% (G4), 5.45% (G6), and 7.27% (G8) immunological activity of glycinin were fed to juvenile hybrid yellow catfish to reveal the mechanism of the intestinal immune response to glycinin through RNA and microRNA (miRNA) sequencing and to explore the interrelation between immune molecules and intestinal microbiota. The results demonstrated that glycinin content in the posterior intestine increased significantly and linearly with the rise of dietary glycinin levels. More than 5.45% of dietary glycinin significantly reduced the nutritional digestion and absorption function of the posterior intestine. Notably, an obvious alteration in the expression levels of inflammatory genes (tnf-α, il-1β, il-15, and tgf-β1) of the posterior intestine was observed when dietary glycinin exceeded 3.57%. Sequencing results of RNA and miRNA deciphered 4,246 differentially expressed genes (DEGs) and 28 differentially expressed miRNAs (DEmiRNAs) between the CK and G6 groups. Furthermore, enrichment analysis of DEGs and DEmiRNA target genes exhibited significant responses of the MAPK, NF-κB, and WNT pathways following experimental fish exposure to 5.45% dietary glycinin. Additionally, at the level of 3.57% in the diet, glycinin obviously inhibited the increase of microbiota, especially potential probiotics such as Ruminococcus bromii, Bacteroides plebeius, Faecalibacterium prausnitzii, and Clostridium clostridioforme. In sum, 5.45% dietary glycinin through the MAPK/NF-κB/WNT pathway induces enteritis, and inflammatory conditions could disrupt micro-ecological equilibrium through miRNA secreted by the host in hybrid yellow catfish. This study constitutes a comprehensive transcriptional perspective of how intestinal immunity responds to excessive glycinin in fish intestines.

The alleviative effect of Bacillus subtilis-supplemented diet against Vibrio cholerae infection in Nile tilapia (Oreochromis niloticus)

Bacterial illness causes detrimental impacts on fish health and survival and finally economic losses for the aquaculture industry. Antibiotic medication causes microbial resistance, so alternative control strategies should be applied. In this work, we investigated the probiotic-medicated diet as an alternative control approach for antibiotics in treating Vibrio cholerae infection in Nile tilapia (Oreochromis niloticus). One hundred eighty fish (50 ± 2.5 g Mean ± SD) were allocated into six groups in glass aquariums (96 L) in triplicate for 10 days. Groups 1 (G1), G2, and G 3 were intraperitoneally (IP) injected with 0.5 mL sterilized tryptic soy broth and fed on a basal diet, basal diet contained B. subtilis (BS) (1 × 10 5 CFU/ kg-1 diet), and basal diet contained trimethoprim-sulfamethoxazole (TMP-SMX) (1.5 g/kg-1 diet), respectively. Additionally, G4, G5, and G6 were IP challenged with 0.5 mL of V. cholerae (1.5 × 107 CFU) and received the same feeding regime as G 1 to 3, respectively. The results exhibited that the V. cholera-infected fish exhibited skin hemorrhage, fin rot, and the lowest survival (63.33%). Additionally, lowered immune-antioxidant biomarkers (white blood cells count, serum bactericidal activity, phagocytic activity, phagocytic index, and lysozymes) with higher lipid peroxidation marker (malondialdehyde) were consequences of V. cholerae infection. Noteworthy, fish-fed therapeutic diets fortified with BS and TMP-SMX showed a substantial amelioration in the clinical signs and survival. The BS diet significantly improved (P < 0.05) the immune-antioxidant indices of the infected fish compared to the TMP-SMX diet. The current findings supported the use of a BS-enriched diet as an eco-friendly approach for the control of V. cholerae in O. niloticus.

The killifish germline regulates longevity and somatic repair in a sex-specific manner

Classical evolutionary theories propose tradeoffs among reproduction, damage repair and lifespan. However, the specific role of the germline in shaping vertebrate aging remains largely unknown. In this study, we used the turquoise killifish (Nothobranchius furzeri) to genetically arrest germline development at discrete stages and examine how different modes of infertility impact life history. We first constructed a comprehensive single-cell gonadal atlas, providing cell-type-specific markers for downstream phenotypic analysis. We show here that germline depletion-but not arresting germline differentiation-enhances damage repair in female killifish. Conversely, germline-depleted males instead showed an extension in lifespan and rejuvenated metabolic functions. Through further transcriptomic analysis, we highlight enrichment of pro-longevity pathways and genes in germline-depleted male killifish and demonstrate functional conservation of how these factors may regulate longevity in germline-depleted Caenorhabditis elegans. Our results, therefore, demonstrate that different germline manipulation paradigms can yield pronounced sexually dimorphic phenotypes, implying alternative responses to classical evolutionary tradeoffs.

The killifish germline regulates longevity and somatic repair in a sex-specific manner

Classical evolutionary theories propose tradeoffs between reproduction, damage repair, and lifespan. However, the specific role of the germline in shaping vertebrate aging remains largely unknown. Here, we use the turquoise killifish ( N. furzeri ) to genetically arrest germline development at discrete stages, and examine how different modes of infertility impact life-history. We first construct a comprehensive single-cell gonadal atlas, providing cell-type-specific markers for downstream phenotypic analysis. Next, we show that germline depletion - but not arresting germline differentiation - enhances damage repair in female killifish. Conversely, germline-depleted males instead showed an extension in lifespan and rejuvenated metabolic functions. Through further transcriptomic analysis, we highlight enrichment of pro-longevity pathways and genes in germline-depleted male killifish and demonstrate functional conservation of how these factors may regulate longevity in germline-depleted C. elegans . Our results therefore demonstrate that different germline manipulation paradigms can yield pronounced sexually dimorphic phenotypes, implying alternative responses to classical evolutionary tradeoffs.

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.

Cell atlas of the Atlantic salmon spleen reveals immune cell heterogeneity and cell-specific responses to bacterial infection

The spleen is a conserved secondary lymphoid organ that emerged in parallel to adaptive immunity in early jawed vertebrates. Recent studies have applied single cell transcriptomics to reveal the cellular composition of spleen in several species, cataloguing diverse immune cell types and subpopulations. In this study, 51,119 spleen nuclei transcriptomes were comprehensively investigated in the commercially important teleost Atlantic salmon (Salmo salar L.), contrasting control animals with those challenged with the bacterial pathogen Aeromonas salmonicida. We identified clusters of nuclei representing the expected major cell types, namely T cells, B cells, natural killer-like cells, granulocytes, mononuclear phagocytes, endothelial cells, mesenchymal cells, erythrocytes and thrombocytes. We discovered heterogeneity within several immune lineages, providing evidence for resident macrophages and melanomacrophages, infiltrating monocytes, several candidate dendritic cell subpopulations, and B cells at distinct stages of differentiation, including plasma cells and an igt + subset. We provide evidence for twelve candidate T cell subsets, including cd4+ T helper and regulatory T cells, one cd8+ subset, three γδT subsets, and populations double negative for cd4 and cd8. The number of genes showing differential expression during the early stages of Aeromonas infection was highly variable across immune cell types, with the largest changes observed in macrophages and infiltrating monocytes, followed by resting mature B cells. Our analysis provides evidence for a local inflammatory response to infection alongside B cell maturation in the spleen, and upregulation of ccr9 genes in igt + B cells, T helper and cd8+ cells, and monocytes, consistent with the recruitment of immune cell populations to the gut to deal with Aeromonas infection. Overall, this study provides a new cell-resolved perspective of the immune actions of Atlantic salmon spleen, highlighting extensive heterogeneity hidden to bulk transcriptomics. We further provide a large catalogue of cell-specific marker genes that can be leveraged to further explore the function and structural organization of the salmonid immune system.

Single-cell RNA-seq mapping of chicken peripheral blood leukocytes

BACKGROUND

Single-cell transcriptomics provides means to study cell populations at the level of individual cells. In leukocyte biology this approach could potentially aid the identification of subpopulations and functions without the need to develop species-specific reagents. The present study aimed to evaluate single-cell RNA-seq as a tool for identification of chicken peripheral blood leukocytes. For this purpose, purified and thrombocyte depleted leukocytes from 4 clinically healthy hens were subjected to single-cell 3' RNA-seq. Bioinformatic analysis of data comprised unsupervised clustering of the cells, and annotation of clusters based on expression profiles. Immunofluorescence phenotyping of the cell preparations used was also performed.

RESULTS

Computational analysis identified 31 initial cell clusters and based on expression of defined marker genes 28 cluster were identified as comprising mainly B-cells, T-cells, monocytes, thrombocytes and red blood cells. Of the remaining clusters, two were putatively identified as basophils and eosinophils, and one as proliferating cells of mixed origin. In depth analysis on gene expression profiles within and between the initial cell clusters allowed further identification of cell identity and possible functions for some of them. For example, analysis of the group of monocyte clusters revealed subclusters comprising heterophils, as well as putative monocyte subtypes. Also, novel aspects of TCRγ/δ + T-cell subpopulations could be inferred such as evidence of at least two subtypes based on e.g., different expression of transcription factors MAF, SOX13 and GATA3. Moreover, a novel subpopulation of chicken peripheral B-cells with high SOX5 expression was identified. An overall good correlation between mRNA and cell surface phenotypic cell identification was shown.

CONCLUSIONS

Taken together, we were able to identify and infer functional aspects of both previously well known as well as novel chicken leukocyte populations although some cell types. e.g., T-cell subtypes, proved more challenging to decipher. Although this methodology to some extent is limited by incomplete annotation of the chicken genome, it definitively has benefits in chicken immunology by expanding the options to distinguish identity and functions of immune cells also without access to species specific reagents.

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.

NITR12+ NK Cells Release Perforin to Mediate IgMhi B Cell Killing in Turbot (Scophthalmus maximus)

B lymphocytes engaged in humoral immunity play a critical role in combating pathogenic infections; however, the mechanisms of NK cells in regulating the responses of B cells remain largely unknown. In the present study, we established an Edwardsiella piscicida infection model in turbot (Scophthalmus maximus) and found that the production of IgM was decreased. Meanwhile, through establishing the head kidney-derived lymphocyte infection model, we revealed that the impairment of IgMhi B cells was associated with bacterial infection-induced perforin production. Interestingly, we reveal that perforin production in NK cells is tightly regulated by an inhibitory novel immune-type receptor, NITR12. Moreover, we confirm that inhibiting NITR12 can result in elevated perforin production, engaging the impairment of IgMhi B cells. Taken together, these findings demonstrate an innovative strategy of NK cells in mediating B lymphocyte killing in turbot and suggest that relieving NK cells through NITR12 might be the target for the development of efficacious vaccines.

CD4+ T lymphocyte responses to viruses and virus-relevant stimuli in teleost fish

Fish diseases caused by viruses are a major threat to aquaculture. Development of disease protection strategies for sustainable fish aquaculture requires a better understanding of the immune mechanisms involved in antiviral defence. The innate and adaptive arms of the vertebrate immune system collaborate to mount an effective defence against viral pathogens. The T lymphocyte components of the adaptive immune system, comprising two major classes (helper T, Th or CD4+ and cytotoxic T lymphocytes, CTLs or CD8+ T cells), are responsible for cell-mediated immune responses. In particular, CD4+ T cells and their different subsets orchestrate the actions of various other immune cells during immune responses, making CD4+ T cells central drivers of responses to pathogens and vaccines. CD4+ T cells are also present in teleost fish. Here we review the literature that reported the use of antibodies against CD4 in a few teleost fish species and transcription profiling of Th cell-relevant genes in the context of viral infections and virus-relevant immunomodulation. Studies reveal massive CD4+ T cell proliferation and expression of key cytokines, transcription factors, and effector molecules that evoke mammalian Th cell responses. We also discuss gaps in the current understanding and evaluation of teleost CD4+ T cell responses and how development and application of novel tools and approaches to interrogate such responses could bridge these gaps. A greater understanding of fish Th cell responses will further illuminate the evolution of vertebrate adaptive immunity, inform strategies to address viral infections in aquaculture, and could further foster fish as model organisms.

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.

Single-cell transcriptome, phagocytic activity and immunohistochemical analysis of crucian carp (Carassius auratus) in response to Rahnella aquatilis infection

In teleost fish, kidney is an important immune and hematopoietic organ with multiple physiological functions. However, the immune cells and cellular markers of kidney require further elucidation in crucian carp (C. auratus). Here we report on the single-cell transcriptional landscape in posterior kidney, immunohistochemical and phagocytic features of C. auratus with R. aquatilis infection. The results showed that a total of 18 cell populations were identified for the main immune cells such as monocytes/macrophages (Mo/Mφ), dendritic cells (DCs), B cells, T cells, granulocytes and hematopoietic progenitor cells (HPCs). Pseudo-time trajectory analysis was reconstructed for the immune cells using Monocle2 to obtain additional insights into their developmental lineage relationships. In the detected tissues (liver, spleen, kidney, intestine, skin, and gills) of infected fish exhibited positive immunohistochemical staining with prepared for antibody to R. aquatilis. Apoptotic cells were fluorescently demonstrated by TUNEL assay, and bacterial phagocytic activity were observed for neutrophils and Mo/Mφ cells, respectively. Moreover, a similar up-ward/down-ward expression trend of the selected immune and inflammatory genes was found in the kidney against R. aquatilis infection, which were significantly involved in TLR/NLR, ECM adhesion, phago-lysosome, apoptosis, complement and coagulation pathways. To our knowledge, this is the first report on the detailed characterization of immune cells and host-R. aquatilis interaction, which will contribute to understanding on the biology of renal immune cells and repertoire of potential markers in cyprinid fish species.

Screening of potential immune-related genes expressed during sepsis using gene sequencing technology

To screen potential pivotal targets in sepsis through peripheral blood. Septic patients (n = 23) and healthy volunteers (n = 10) were enrolled according to SEPSIS 3.0. Peripheral blood was collected within 24 h of enrollment, RNA-seq was performed on the peripheral blood. The sequencing data was screened for DEGs (p < 0.01; logFC ≥ 2). PPI, WGCNA and survival curve analysis were used to identify potential targets. Then, 5 PBMC samples were conducted by single-cell sequencing for cell lineage location. Finally, mouse sepsis model and clinic samples were performed to verify the targets gene using RNA-seq and RT-PCR, respectively. Compared to the control group, 1007 DEGs were found in septic group. BCL9L, BCL11B, CD247, CD96, MAFG and SAMD3 were in the core of network. These six genes correlated to the survival rate of septic patients and they were mainly expressed in T cells, except that MAFG was located in monocyte cell. The expression levels of six key genes were confirmed by animal and clinical samples. BCL9L, BCL11B, CD247, CD96 and SAMD3 were decreased in sepsis and mainly expressed in the T cell; while MAFG increased in sepsis and localizes to monocytes. These genes may be therapeutic targets for sepsis.

Mass Spectrometric Identification of Licania rigida Benth Leaf Extracts and Evaluation of Their Therapeutic Effects on Lipopolysaccharide-Induced Inflammatory Response

Licania rigida Benth has been evaluated as an alternative drug to treat diseases associated with inflammatory processes. This study evaluated the anti-inflammatory effects of aqueous and hydroalcoholic leaf extracts of L. rigida with inflammation induced by lipopolysaccharides in in vitro and in vivo inflammation models. The phytochemical profile of the extracts, analyzed by ultra-fast liquid chromatography coupled with tandem mass spectrometry, revealed the presence of gallic and ellagic acids in both extracts, whereas isovitexin, ferulate, bulky amino acids (e.g., phenylalanine), pheophorbide, lactic acid, and pyridoxine were detected in the hydroalcoholic extract. The extracts displayed the ability to modulate in vitro and in vivo inflammatory responses, reducing approximately 50% of pro-inflammatory cytokine secretion (TNF-α, IL-1β, and IL-6), and inhibiting both NO production and leukocyte migration by approximately 30 and 40% at 100 and 500 µg/mL, respectively. Overall, the results highlight and identify, for the first time, the ability of L. rigida leaf extract to modulate inflammatory processes. These data suggest that the leaf extracts of this plant have potential in the development of herbal formulations for the treatment of inflammation.

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.