Single-cell transcriptome analysis of fish immune cells provides insight into the evolution of vertebrate immune cell types

The inflammasome adaptor pycard is essential for immunity against Mycobacterium marinum infection in adult zebrafish

Inflammasomes regulate the host response to intracellular pathogens including mycobacteria. We have previously shown that the course of Mycobacterium marinum infection in adult zebrafish (Danio rerio) mimics the course of tuberculosis in human. To investigate the role of the inflammasome adaptor pycard in zebrafish M. marinum infection, we produced two zebrafish knockout mutant lines for the pycard gene with CRISPR/Cas9 mutagenesis. Although the zebrafish larvae lacking pycard developed normally and had unaltered resistance against M. marinum, the loss of pycard led to impaired survival and increased bacterial burden in the adult zebrafish. Based on histology, immune cell aggregates, granulomas, were larger in pycard-deficient fish than in wild-type controls. Transcriptome analysis with RNA sequencing of a zebrafish haematopoietic tissue, kidney, suggested a role for pycard in neutrophil-mediated defence, haematopoiesis and myelopoiesis during infection. Transcriptome analysis of fluorescently labelled, pycard-deficient kidney neutrophils identified genes that are associated with compromised resistance, supporting the importance of pycard for neutrophil-mediated immunity against M. marinum. Our results indicate that pycard is essential for resistance against mycobacteria in adult zebrafish.

How fish intestinal cells responded to dietary methylmercury exposure? A single-cell transcriptomic analysis

Fish intestine is not only an important digestive and immune organ, but also serves as the first barrier to defend against methylmercury (MeHg) toxicity. Numerous studies have examined the responses of intestine to MeHg, whereas the heterogeneous responses of intestinal cells have not been addressed. In this study, the gilthead seabream were exposed to dietary MeHg, and the gene expression profiles of different intestinal cell populations were examined using scRNA-seq technique. We demonstrated that among the 14 cell types identified, enterocytes, macrophages, T cells and goblet cells were the primary target cell populations exhibiting specific responses to MeHg. Enterocytes appeared to play the most important role in the MeHg transport across the intestinal epithelium as well as intracellular storage. The immune pathways of macrophages and T cells were suppressed by MeHg, which also interfered with the mucus production and secretion in the goblet cells. Furthermore, MeHg not only affected the cell-cell adhesion of the target cells, but also resulted in disorder of lipid metabolism and immune function, thereby leading to increased susceptibility to pathogenic infections. This study provides an important understanding of the specific responses of intestinal cells to MeHg exposure at the cellular level.

scVAEDer: integrating deep diffusion models and variational autoencoders for single-cell transcriptomics analysis

Discovering a lower-dimensional embedding of single-cell data can improve downstream analysis. The embedding should encapsulate both the high-level features and low-level variations. While existing generative models attempt to learn such low-dimensional representations, they have limitations. Here, we introduce scVAEDer, a scalable deep-learning model that combines the power of variational autoencoders and deep diffusion models to learn a meaningful representation that retains both global structure and local variations. Using the learned embeddings, scVAEDer can generate novel scRNA-seq data, predict perturbation response on various cell types, identify changes in gene expression during dedifferentiation, and detect master regulators in biological processes.

Single-cell transcriptome analysis of medaka lymphocytes reveals absence of fully mature T cells in the thymus and the T-lineage commitment in the kidney

The cellular and molecular mechanisms underlying lymphocyte development are diverse among teleost species. Although recent scRNA-seq analyses of zebrafish hematopoietic cells have advanced our understanding of teleost hematopoiesis, comparative studies using another genetic model, medaka, which is evolutionarily distant among teleosts, is useful for understanding commonality and species-specificity in teleosts. In order to gain insight into how different molecular and cellular mechanisms of lymphocyte development in medaka and zebrafish, we established a recombination activating gene 1 (rag1) mutant medaka, which exhibited defects in V(D)J rearrangement of lymphocyte antigen receptor genes, accordingly lacking mature B and T cells. scRNA-seq analysis of wild type and rag1 mutant lymphocytes in the thymus and kidney characterized the developing stages of T and B cells, and found that most developed cd4+cd8- and cd4-cd8+ single-positive (SP) T-cell populations are absent in the thymus, and identified lymphoid progenitor cells already committed to the T lineage in kidney, implying unique features of medaka lymphocyte development.

Diverse Epithelial Lymphocytes in Zebrafish Revealed Using a Novel Scale Biopsy Method

Zebrafish (Danio rerio) are a compelling model for studying lymphocytes because zebrafish and humans have similar adaptive immune systems, including their lymphocytes. Antibodies that recognize zebrafish proteins are sparse, so many investigators use transgenic, lymphocyte-specific fluorophore-labeled lines. Human and zebrafish lymphocyte types are conserved, but many aspects of zebrafish lymphocyte biology remain uninvestigated, including lymphocytes in peripheral tissues, like epidermis. This study is, to our knowledge, the first study to focus on zebrafish epidermal lymphocytes, using scales. Obtaining zebrafish blood via nonlethal methods is difficult; scales represent a source to longitudinally sample live fish. We developed a novel biopsy technique, collecting scales to analyze epithelial lymphocytes from several transgenic lines. We imaged scales via confocal microscopy and demonstrated multiple lymphocyte types in scales/epidermis, quantifying them flow cytometrically. We profiled gene expression of scale, thymic, and kidney-marrow (analogous to mammalian bone marrow) lymphocytes from the same animals, revealing B- and T-lineage signatures. Single-cell quantitative real-time PCR and RNA sequencing show not only canonical B and T cells but also novel lymphocyte populations not described previously. To validate longitudinal scale biopsies, we serially sampled scales from fish treated with dexamethasone, demonstrating epidermal lymphocyte responses. To analyze cells functionally, we employed a bead-ingestion assay, showing that thymic, marrow, and epidermal lymphocytes have phagocytic activity. In summary, we establish a novel, nonlethal technique to obtain zebrafish lymphocytes, providing the first quantification, expression profiling, and functional data from zebrafish epidermal lymphocytes.

Cross-species single-cell analysis reveals divergence and conservation of peripheral blood mononuclear cells

BACKGROUND

Single-cell transcriptome sequencing (scRNA-seq) has revolutionized the study of immune cells by overcoming the limitations of traditional antibody-based identification and isolation methods. This advancement allows us to obtain comprehensive gene expression profiles from a diverse array of vertebrate species, facilitating the identification of various cell types. Comparative immunology across vertebrates presents a promising approach to understanding the evolution of immune cell types. In this study, we conducted a comparative transcriptome analysis of peripheral blood mononuclear cells (PBMCs) at the single-cell level across 12 species.

RESULTS

Our findings shed light on the cellular compositional features of PBMCs, spanning from fish to mammals. Notably, we identified genes that exhibit vertebrate universality in characterizing immune cells. Moreover, our investigation revealed that monocytes have maintained a conserved transcriptional regulatory program throughout evolution, emphasizing their pivotal role in orchestrating immune cells to execute immune programs.

CONCLUSIONS

This comprehensive analysis provides valuable insights into the evolution of immune cells across vertebrates.

Single-cell transcriptome landscape of the kidney reveals potential innate immune regulation mechanisms in hybrid yellow catfish after Aeromonas hydrophila infection

Aeromonas hydrophila, the pathogen that is the causative agent of motile Aeromonas septicemia (MAS) disease, commonly attacks freshwater fishes, including yellow catfish (Pelteobagrus fulvidraco). Although the kidney is one of the most important organs involved in immunity in fish, its role in disease progression has not been fully elucidated. Understanding the cellular composition and innate immune regulation mechanisms of the kidney of yellow catfish is important for the treatment of MAS. In this study, single-cell RNA sequencing (scRNA-seq) was performed on the kidney of hybrid yellow catfish (Pelteobagrus fulvidraco ♀ × Pelteobagrus vachelli ♂) after A. hydrophila infection. Nine types of kidney cells were identified using marker genes, and a transcription module of marker genes in the main immune cells of hybrid yellow catfish kidney tissue was constructed using in-situ hybridization. In addition, the single-cell transcriptome data showed that the differentially expressed genes of macrophages were primarily enriched in the Toll-like receptor and Nod-like receptor signaling pathways. The expression levels of genes involved in these pathways were upregulated in macrophages following A. hydrophila infection. Transmission electron microscopy and TUNEL analysis revealed the cellular characteristics of macrophages before and after A. hydrophila infection. These data provide empirical support for in-depth research on the role of the kidney in the innate immune response of hybrid yellow catfish.

GATA2 heterozygosity causes an epigenetic feedback mechanism resulting in myeloid and erythroid dysplasia

The transcription factor GATA2 has a pivotal role in haematopoiesis. Heterozygous germline GATA2 mutations result in a syndrome characterized by immunodeficiency, bone marrow failure and predispositions to myelodysplastic syndrome (MDS) and acute myeloid leukaemia. Clinical symptoms in these patients are diverse and mechanisms driving GATA2-related phenotypes are largely unknown. To explore the impact of GATA2 haploinsufficiency on haematopoiesis, we generated a zebrafish model carrying a heterozygous mutation of gata2b (gata2b+/-), an orthologue of GATA2. Morphological analysis revealed myeloid and erythroid dysplasia in gata2b+/- kidney marrow. Because Gata2b could affect both transcription and chromatin accessibility during lineage differentiation, this was assessed by single-cell (sc) RNA-seq and single-nucleus (sn) ATAC-seq. Sn-ATAC-seq showed that the co-accessibility between the transcription start site (TSS) and a -3.5-4.1 kb putative enhancer was more robust in gata2b+/- zebrafish HSPCs compared to wild type, increasing gata2b expression and resulting in higher genome-wide Gata2b motif use in HSPCs. As a result of increased accessibility of the gata2b locus, gata2b+/- chromatin was also more accessible during lineage differentiation. scRNA-seq data revealed myeloid differentiation defects, that is, impaired cell cycle progression, reduced expression of cebpa and cebpb and increased signatures of ribosome biogenesis. These data also revealed a differentiation delay in erythroid progenitors, aberrant proliferative signatures and down-regulation of Gata1a, a master regulator of erythropoiesis, which worsened with age. These findings suggest that cell-intrinsic compensatory mechanisms, needed to obtain normal levels of Gata2b in heterozygous HSPCs to maintain their integrity, result in aberrant lineage differentiation, thereby representing a critical step in the predisposition to MDS.

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.

Direct male development in chromosomally ZZ zebrafish

The genetics of sex determination varies across taxa, sometimes even within a species. Major domesticated strains of zebrafish (Danio rerio), including AB and TU, lack a strong genetic sex determining locus, but strains more recently derived from nature, like Nadia (NA), possess a ZZ male/ZW female chromosomal sex-determination system. AB fish pass through a juvenile ovary stage, forming oocytes that survive in fish that become females but die in fish that become males. To understand mechanisms of gonad development in NA zebrafish, we studied histology and single cell transcriptomics in developing ZZ and ZW fish. ZW fish developed oocytes by 22 days post-fertilization (dpf) but ZZ fish directly formed testes, avoiding a juvenile ovary phase. Gonads of some ZW and WW fish, however, developed oocytes that died as the gonad became a testis, mimicking AB fish, suggesting that the gynogenetically derived AB strain is chromosomally WW. Single-cell RNA-seq of 19dpf gonads showed similar cell types in ZZ and ZW fish, including germ cells, precursors of gonadal support cells, steroidogenic cells, interstitial/stromal cells, and immune cells, consistent with a bipotential juvenile gonad. In contrast, scRNA-seq of 30dpf gonads revealed that cells in ZZ gonads had transcriptomes characteristic of testicular Sertoli, Leydig, and germ cells while ZW gonads had granulosa cells, theca cells, and developing oocytes. Hematopoietic and vascular cells were similar in both sex genotypes. These results show that juvenile NA zebrafish initially develop a bipotential gonad; that a factor on the NA W chromosome, or fewer than two Z chromosomes, is essential to initiate oocyte development; and without the W factor, or with two Z doses, NA gonads develop directly into testes without passing through the juvenile ovary stage. Sex determination in AB and TU strains mimics NA ZW and WW zebrafish, suggesting loss of the Z chromosome during domestication. Genetic analysis of the NA strain will facilitate our understanding of the evolution of sex determination mechanisms.

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.

Zebrafish Thrombocyte Transcriptome Analysis and Functional Genomics

Our laboratory is interested in investigating the maturation process of zebrafish thrombocytes, which are functional equivalents to human platelets. We have adopted the zebrafish model to gain insights into mammalian platelet production, or thrombopoiesis. Notably, zebrafish exhibit two distinct populations of thrombocytes in their circulating blood: young and mature thrombocytes. This observation is intriguing because maturation appears to occur in circulation, yet the precise mechanisms governing this maturation remain elusive. Our goal is to understand the mechanisms underlying thrombocyte maturation by conducting single-cell RNA sequencing (scRNA-Seq) on young and mature thrombocytes, analyzing these transcriptomes to identify genes specific to each thrombocyte population, and elucidating the role of these genes in the maturation process, by quantifying thrombocyte numbers after the piggyback knockdown of each of these genes. In this chapter, we present a comprehensive, step-by-step protocol detailing the multifaceted methodology involved in understanding thrombocyte maturation, which encompasses the collection of zebrafish blood, the separation of young and mature thrombocytes using flow cytometry, scRNA-Seq analysis of these distinct thrombocyte populations, identification of genes specific to young and mature thrombocytes, and subsequent validation through gene knockdown techniques.

Direct Male Development in Chromosomally ZZ Zebrafish

The genetics of sex determination varies across taxa, sometimes even within a species. Major domesticated strains of zebrafish ( Danio rerio ), including AB and TU, lack a strong genetic sex determining locus, but strains more recently derived from nature, like Nadia (NA), possess a ZZ male/ZW female chromosomal sex-determination system. AB strain fish pass through a juvenile ovary stage, forming oocytes that survive in fish that become females but die in fish that become males. To understand mechanisms of gonad development in NA zebrafish, we studied histology and single cell transcriptomics in developing ZZ and ZW fish. ZW fish developed oocytes by 22 days post-fertilization (dpf) but ZZ fish directly formed testes, avoiding a juvenile ovary phase. Gonads of some ZW and WW fish, however, developed oocytes that died as the gonad became a testis, mimicking AB fish, suggesting that the gynogenetically derived AB strain is chromosomally WW. Single-cell RNA-seq of 19dpf gonads showed similar cell types in ZZ and ZW fish, including germ cells, precursors of gonadal support cells, steroidogenic cells, interstitial/stromal cells, and immune cells, consistent with a bipotential juvenile gonad. In contrast, scRNA-seq of 30dpf gonads revealed that cells in ZZ gonads had transcriptomes characteristic of testicular Sertoli, Leydig, and germ cells while ZW gonads had granulosa cells, theca cells, and developing oocytes. Hematopoietic and vascular cells were similar in both sex genotypes. These results show that juvenile NA zebrafish initially develop a bipotential gonad; that a factor on the NA W chromosome or fewer than two Z chromosomes is essential to initiate oocyte development; and without the W factor or with two Z doses, NA gonads develop directly into testes without passing through the juvenile ovary stage. Sex determination in AB and TU strains mimics NA ZW and WW zebrafish, suggesting loss of the Z chromosome during domestication. Genetic analysis of the NA strain will facilitate our understanding of the evolution of sex determination mechanisms.

The Immune Response to the Myxozoan Parasite Myxobolus cerebralis in Salmonids: A Review on Whirling Disease

Salmonids are affected by the economically significant whirling disease (WD) caused by the myxozoan parasite Myxobolus cerebralis. In the past, it was endemic to Eurasia, but it has now spread to different regions of North America, Europe, New Zealand, and South Africa. Among salmonids, rainbow trout is considered the most highly susceptible host. Upon entering to the host's body, the parasite invades the spine and cranium, resulting in whirling behaviour, a blackened tail, and destruction of cartilage. The disease is characterized by the infiltration of numerous inflammatory cells, primarily lymphocytes and macrophages, with the onset of fibrous tissue infiltration. Several efforts have been undertaken to investigate the role of various immune modulatory molecules and immune regulatory genes using advanced molecular methods including flow cytometry and transcriptional techniques. Investigation of the molecular and cellular responses, the role of STAT3 in Th17 cell differentiation, and the inhibitory actions of suppressors of cytokine signaling (SOCS) on interferons and interleukins, as well as the role of natural resistance-associated macrophage proteins (Nramp) in WD have significantly contributed to our understanding of the immune regulation mechanism in salmonids against M. cerebralis. This review thoroughly highlights previous research and discusses potential future directions for understanding the molecular immune response of salmonids and the possible development of prophylactic approaches against WD.

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.

Identification of a pharyngeal mucosal lymphoid organ in zebrafish and other teleosts: Tonsils in fish?

The constant exposure of the fish branchial cavity to aquatic pathogens causes local mucosal immune responses to be extremely important for their survival. Here, we used a marker for T lymphocytes/natural killer (NK) cells (ZAP70) and advanced imaging techniques to investigate the lymphoid architecture of the zebrafish branchial cavity. We identified a sub-pharyngeal lymphoid organ, which we tentatively named "Nemausean lymphoid organ" (NELO). NELO is enriched in T/NK cells, plasma/B cells, and antigen-presenting cells embedded in a network of reticulated epithelial cells. The presence of activated T cells and lymphocyte proliferation, but not V(D)J recombination or hematopoiesis, suggests that NELO is a secondary lymphoid organ. In response to infection, NELO displays structural changes including the formation of T/NK cell clusters. NELO and gill lymphoid tissues form a cohesive unit within a large mucosal lymphoid network. Collectively, we reveal an unreported mucosal lymphoid organ reminiscent of mammalian tonsils that evolved in multiple teleost fish families.

Spatio-temporal patterns of ovarian development and VgR gene silencing reduced fecundity in parthenogenetic Artemia

The halophilic zooplankton brine shrimp Artemia has been used as an experimental animal in multidisciplinary studies. However, the reproductive patterns and its regulatory mechanisms in Artemia remain unclear. In this study, the ovarian development process of parthenogenetic Artemia (A. parthenogenetica) was divided into five stages, and oogenesis or egg formation was identified in six phases. The oogenesis mode was assumed to be polytrophic. We also traced the dynamic translocation of candidate germline stem cells (cGSCs) using EdU labelling and elucidated several key cytological events in oogenesis through haematoxylin and eosin staining and fluorescence imaging. Distinguished from the ovary structure of insects and crustaceans, Artemia germarium originated from ovariole buds and are located at the base of the ovarioles. RNA-seq based on five stages of ovarian development identified 2657 upregulated genes related to reproduction by pair-to-pair comparison. Gbb, Dpp, piwi, vasa, nanos, VgA and VgR genes associated with cGSCs recognition and reproductive development were screened and verified using qPCR. Silencing of the VgR gene in A. parthenogenetica (Ap-VgR) at ovarian development Stage II led to a low level of gene expression (less than 10%) within 5 days, which resulted in variations in oogenesis-related gene expression and significantly inhibited vitellogenesis, impeded oocyte maturation, and eventually decreased the number of offspring. In conclusion, we have illustrated the patterns of ovarian development, outlined the key spatio-temporal features of oogenesis and identified the negative impacts of VgR gene knockdown on oogenesis using A. parthenogenetica as an experimental animal. The findings of this study also lay a foundation for the further study of reproductive biology of invertebrates.

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.

Ligament injury in adult zebrafish triggers ECM remodeling and cell dedifferentiation for scar-free regeneration

After traumatic injury, healing of mammalian ligaments is typically associated with fibrotic scarring as opposed to scar-free regeneration. In contrast, here we show that the ligament supporting the jaw joint of adult zebrafish is capable of rapid and complete scar-free healing. Following surgical transection of the jaw joint ligament, we observe breakdown of ligament tissue adjacent to the cut sites, expansion of mesenchymal tissue within the wound site, and then remodeling of extracellular matrix (ECM) to a normal ligament morphology. Lineage tracing of mature ligamentocytes following transection shows that they dedifferentiate, undergo cell cycle re-entry, and contribute to the regenerated ligament. Single-cell RNA sequencing of the regenerating ligament reveals dynamic expression of ECM genes in neural-crest-derived mesenchymal cells, as well as diverse immune cells expressing the endopeptidase-encoding gene legumain. Analysis of legumain mutant zebrafish shows a requirement for early ECM remodeling and efficient ligament regeneration. Our study establishes a new model of adult scar-free ligament regeneration and highlights roles of immune-mesenchyme cross-talk in ECM remodeling that initiates regeneration.

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.

Single cell transcriptome analyses of the developing zebrafish eye- perspectives and applications

Within a relatively short period of time, single cell transcriptome analyses (SCT) have become increasingly ubiquitous with transcriptomic research, uncovering plentiful details that boost our molecular understanding of various biological processes. Stemming from SCT analyses, the ever-growing number of newly assigned genetic markers increases our understanding of general function and development, while providing opportunities for identifying genes associated with disease. SCT analyses have been carried out using tissue from numerous organisms. However, despite the great potential of zebrafish as a model organism, other models are still preferably used. In this mini review, we focus on eye research as an example of the advantages in using zebrafish, particularly its usefulness for single cell transcriptome analyses of developmental processes. As studies have already shown, the unique opportunities offered by zebrafish, including similarities to the human eye, in combination with the possibility to analyze and extract specific cells at distinct developmental time points makes the model a uniquely powerful one. Particularly the practicality of collecting large numbers of embryos and therefore isolation of sufficient numbers of developing cells is a distinct advantage compared to other model organisms. Lastly, the advent of highly efficient genetic knockouts methods offers opportunities to characterize target gene function in a more cost-efficient way. In conclusion, we argue that the use of zebrafish for SCT approaches has great potential to further deepen our molecular understanding of not only eye development, but also many other organ systems.

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.

A single-cell transcriptome atlas of the maturing zebrafish telencephalon

The zebrafish telencephalon is composed of highly specialized subregions that regulate complex behaviors such as learning, memory, and social interactions. The transcriptional signatures of the neuronal cell types in the telencephalon and the timeline of their emergence from larva to adult remain largely undescribed. Using an integrated analysis of single-cell transcriptomes of approximately 64,000 cells obtained from 6-day-postfertilization (dpf), 15-dpf, and adult telencephalon, we delineated nine main neuronal cell types in the pallium and eight in the subpallium and nominated novel marker genes. Comparing zebrafish and mouse neuronal cell types revealed both conserved and absent types and marker genes. Mapping of cell types onto a spatial larval reference atlas created a resource for anatomical and functional studies. Using this multiage approach, we discovered that although most neuronal subtypes are established early in the 6-dpf fish, some emerge or expand in number later in development. Analyzing the samples from each age separately revealed further complexity in the data, including several cell types that expand substantially in the adult forebrain and do not form clusters at the larval stages. Together, our work provides a comprehensive transcriptional analysis of the cell types in the zebrafish telencephalon and a resource for dissecting its development and function.

Olfactory detection of viruses shapes brain immunity and behavior in zebrafish

Olfactory sensory neurons (OSNs) are constantly exposed to pathogens, including viruses. However, serious brain infection via the olfactory route rarely occurs. When OSNs detect a virus, they coordinate local antiviral immune responses to stop virus progression to the brain. Despite effective immune control in the olfactory periphery, pathogen-triggered neuronal signals reach the CNS via the olfactory bulb (OB). We hypothesized that neuronal detection of a virus by OSNs initiates neuroimmune responses in the OB that prevent pathogen invasion. Using zebrafish ( Danio rerio ) as a model, we demonstrate viral-specific neuronal activation of OSNs projecting into the OB, indicating that OSNs are electrically activated by viruses. Further, behavioral changes are seen in both adult and larval zebrafish after viral exposure. By profiling the transcription of single cells in the OB after OSNs are exposed to virus, we found that both microglia and neurons enter a protective state. Microglia and macrophage populations in the OB respond within minutes of nasal viral delivery followed decreased expression of neuronal differentiation factors and enrichment of genes in the neuropeptide signaling pathway in neuronal clusters. Pituitary adenylate-cyclase-activating polypeptide ( pacap ), a known antimicrobial, was especially enriched in a neuronal cluster. We confirm that PACAP is antiviral in vitro and that PACAP expression increases in the OB 1 day post-viral treatment. Our work reveals how encounters with viruses in the olfactory periphery shape the vertebrate brain by inducing antimicrobial programs in neurons and by altering host behavior.

Ligament injury in adult zebrafish triggers ECM remodeling and cell dedifferentiation for scar-free regeneration

After traumatic injury, healing of mammalian ligaments is typically associated with fibrotic scarring as opposed to scar-free regeneration. In contrast, here we show that the ligament supporting the jaw joint of adult zebrafish is capable of rapid and complete scar-free healing. Following surgical transection of the jaw joint ligament, we observe breakdown of ligament tissue adjacent to the cut sites, expansion of mesenchymal tissue within the wound site, and then remodeling of extracellular matrix (ECM) to a normal ligament morphology. Lineage tracing of mature ligamentocytes following transection shows that they dedifferentiate, undergo cell cycle re-entry, and contribute to the regenerated ligament. Single-cell RNA sequencing of the regenerating ligament reveals dynamic expression of ECM genes in neural-crest-derived mesenchymal cells, as well as diverse immune cells expressing the endopeptidase-encoding gene legumain . Analysis of legumain mutant zebrafish shows a requirement for early ECM remodeling and efficient ligament regeneration. Our study establishes a new model of adult scar-free ligament regeneration and highlights roles of immune-mesenchyme cross-talk in ECM remodeling that initiates regeneration.

Highlights

Rapid regeneration of the jaw joint ligament in adult zebrafishDedifferentiation of mature ligamentocytes contributes to regenerationscRNAseq reveals dynamic ECM remodeling and immune activation during regenerationRequirement of Legumain for ECM remodeling and ligament healing.

Loss of autism-candidate CHD8 perturbs neural crest development and intestinal homeostatic balance

Individuals with mutations in CHD8 present with gastrointestinal complaints, yet the underlying mechanisms are understudied. Here, using a stable constitutive chd8 mutant zebrafish model, we found that the loss of chd8 leads to a reduced number of vagal neural crest cells (NCCs), enteric neural and glial progenitors, emigrating from the neural tube, and that their early migration capability was altered. At later stages, although the intestinal colonization by NCCs was complete, we found the decreased numbers of both serotonin-producing enterochromaffin cells and NCC-derived serotonergic neurons, suggesting an intestinal hyposerotonemia in the absence of chd8 Furthermore, transcriptomic analyses revealed an altered expression of key receptors and enzymes in serotonin and acetylcholine signaling pathways. The tissue examination of chd8 mutants revealed a thinner intestinal epithelium accompanied by an accumulation of neutrophils and the decreased numbers of goblet cells and eosinophils. Last, single-cell sequencing of whole intestines showed a global disruption of the immune balance with a perturbed expression of inflammatory interleukins and changes in immune cell clusters. Our findings propose a causal developmental link between chd8, NCC development, intestinal homeostasis, and autism-associated gastrointestinal complaints.

Hemocytes of bivalve mollusks as cellular models in toxicological studies of metals and metal-based nanomaterials

Understanding the impacts of environmental pollutants on immune systems is indispensable in ecological and health risk assessments due to the significance of normal immunological functions in all living organisms. Bivalves as sentinel organisms with vital ecological importance are widely distributed in aquatic environments and their innate immune systems are the sensitive targets of environmental pollutants. As the central component of innate immunity, bivalve hemocytes are endowed with specialized endolysosomal systems for particle internalization and metal detoxification. These intrinsic biological features make them a unique cellular model for metal- and nano-immunotoxicology research. In this review, we firstly provided a general overview of bivalve's innate immunity and the classification and immune functions of hemocytes. We then summarized the recent progress on the interactions of metals and nanoparticles with bivalve hemocytes, with emphasis on the involvement of hemocytes in metal regulation and detoxification, the interactions of hemocytes and nanoparticles at eco/bio-nano interface and hemocyte-mediated immune responses to the exposure of metals and nanoparticles. Finally, we proposed the key knowledge gaps and future research priorities in deciphering the fundamental biological processes of the interactions of environmental pollutants with the innate immune system of bivalves as well as in developing bivalve hemocytes into a promising cellular model for nano-immuno-safety assessment.

Single-cell analyses reveal early thymic progenitors and pre-B cells in zebrafish

The zebrafish has proven to be a valuable model organism for studying hematopoiesis, but relatively little is known about zebrafish immune cell development and functional diversity. Elucidating key aspects of zebrafish lymphocyte development and exploring the breadth of effector functions would provide valuable insight into the evolution of adaptive immunity. We performed single-cell RNA sequencing on ∼70,000 cells from the zebrafish marrow and thymus to establish a gene expression map of zebrafish immune cell development. We uncovered rich cellular diversity in the juvenile and adult zebrafish thymus, elucidated B- and T-cell developmental trajectories, and transcriptionally characterized subsets of hematopoietic stem and progenitor cells and early thymic progenitors. Our analysis permitted the identification of two dendritic-like cell populations and provided evidence in support of the existence of a pre-B cell state. Our results provide critical insights into the landscape of zebrafish immunology and offer a foundation for cellular and genetic studies.

Multi-tissue scRNA-seq reveals immune cell landscape of turbot (Scophthalmus maximus)

In vertebrates, bony fishes possess not only innate immune cells but also T and B cells that are equivalent to those in mammals. However, the precise sub-cluster of immune cells in teleost fish remains largely unknown. Herein, we developed a dynamic bacterial infection model in turbot (Scophthalmus maximus) and created a fish immune cell landscape (FICL) for a primary lymphoid organ (head kidney), a secondary lymphoid organ (spleen), and barrier tissues (gills and posterior intestine). Moreover, through comprehensive characterization of the expression profiles of 16 clusters, including dendritic cells-like (DCs-like), macrophages (MΦs), neutrophils, NK cells, as well as 12 sub-clusters of T and B cells, we found that CD8+ CTLs, CD4-CD8- T, Th17 and ILC3-2 like cells possess a bifunctional role associated with cytotoxicity and immunoregulation during bacterial infection. To our knowledge, these results could provide a useful resource for a better understanding of immune cells in teleost fish and could act as a comprehensive knowledge base for assessing the evolutionary mechanism of adaptive immunity in vertebrates.

Hematopoietic Stem and Progenitor Cells (HSPCs) and Hematopoietic Microenvironment: Molecular and Bioinformatic Studies of the Zebrafish Models

Hematopoietic stem cells (HSCs) reside in a specialized microenvironment in a peculiar anatomic location which regulates the maintenance of stem cells and controls its functions. Recent scientific progress in experimental technologies have enabled the specific detection of epigenetic factors responsible for the maintenance and quiescence of the hematopoietic niche, which has improved our knowledge of regulatory mechanisms. The aberrant role of RNA-binding proteins and their impact on the disruption of stem cell biology have been reported by a number of recent studies. Despite recent modernization in hematopoietic microenvironment research avenues, our comprehension of the signaling mechanisms and interactive pathways responsible for integration of the hematopoietic niche is still limited. In the past few decades, zebrafish usage with regards to exploratory studies of the hematopoietic niche has expanded our knowledge for deeper understanding of novel cellular interactions. This review provides an update on the functional roles of different genetic and epigenetic factors and molecular signaling events at different sections of the hematopoietic microenvironment. The explorations of different molecular approaches and interventions of latest web-based tools being used are also outlined. This will help us to get more mechanistic insights and develop therapeutic options for the malignancies.

Single-cell transcriptome reveals insights into the development and function of the zebrafish ovary

Zebrafish are an established research organism that has made many contributions to our understanding of vertebrate tissue and organ development, yet there are still significant gaps in our understanding of the genes that regulate gonad development, sex, and reproduction. Unlike the development of many organs, such as the brain and heart that form during the first few days of development, zebrafish gonads do not begin to form until the larval stage (≥5 days post-fertilization). Thus, forward genetic screens have identified very few genes required for gonad development. In addition, bulk RNA-sequencing studies that identify genes expressed in the gonads do not have the resolution necessary to define minor cell populations that may play significant roles in the development and function of these organs. To overcome these limitations, we have used single-cell RNA sequencing to determine the transcriptomes of cells isolated from juvenile zebrafish ovaries. This resulted in the profiles of 10,658 germ cells and 14,431 somatic cells. Our germ cell data represents all developmental stages from germline stem cells to early meiotic oocytes. Our somatic cell data represents all known somatic cell types, including follicle cells, theca cells, and ovarian stromal cells. Further analysis revealed an unexpected number of cell subpopulations within these broadly defined cell types. To further define their functional significance, we determined the location of these cell subpopulations within the ovary. Finally, we used gene knockout experiments to determine the roles of foxl2l and wnt9b for oocyte development and sex determination and/or differentiation, respectively. Our results reveal novel insights into zebrafish ovarian development and function, and the transcriptome profiles will provide a valuable resource for future studies.

Comprehensive transcriptome profiling and functional analysis of the meagre (Argyrosomus regius) immune system

Meagre (Argyrosomus regius) belongs to the family Sciaenidae and is a promising candidate for Mediterranean aquaculture diversification. As a relatively recent species in aquaculture, the physiological consequences of the immune system activation in meagre are understudied. Spleen, as a primary lymphoid organ has an essential role in meagre immune and inflammatory responses. In this study, we have evaluated the in vivo effects of lipopolysaccharide (LPS) on the spleen transcriptome of meagre by RNA-seq analysis at 4 and 24 h after injection.

Characterization of Pipefish Immune Cell Populations Through Single-Cell Transcriptomics

Teleost adaptive immune systems have evolved with more flexibility than previously assumed. A particularly enigmatic system to address immune system modifications in the evolutionary past is represented by the Syngnathids, the family of pipefishes, seahorses and seadragons. These small fishes with their unique male pregnancy have lost the spleen as an important immune organ as well as a functional major histocompatibility class II (MHC II) pathway. How these evolutionary changes have impacted immune cell population dynamics have up to this point remained unexplored. Here, we present the first immune cell repertoire characterization of a syngnathid fish (Syngnathus typhle) using single-cell transcriptomics. Gene expression profiles of individual cells extracted from blood and head-kidney clustered in twelve putative cell populations with eight belonging to those with immune function. Upregulated cell marker genes identified in humans and teleosts were used to define cell clusters. While the suggested loss of CD4+ T-cells accompanied the loss of the MHC II pathway was supported, the upregulation of specific subtype markers within the T-cell cluster indicates subpopulations of regulatory T-cells (il2rb) and cytotoxic T-cells (gzma). Utilizing single-cell RNA sequencing this report is the first to characterize immune cell populations in syngnathids and provides a valuable foundation for future cellular classification and experimental work within the lineage.

Real-time imaging of inflammation and its resolution: It's apparent because it's transparent

The ability to directly observe leukocyte behavior in vivo has dramatically expanded our understanding of the immune system. Zebrafish are particularly amenable to the high-resolution imaging of leukocytes during both homeostasis and inflammation. Due to its natural transparency, intravital imaging in zebrafish does not require any surgical manipulation. As a result, zebrafish are particularly well-suited for the long-term imaging required to observe the temporal and spatial events during the onset and resolution of inflammation. Here, we review major insights about neutrophil and macrophage function gained from real-time imaging of zebrafish. We discuss neutrophil reverse migration, the process whereby neutrophils leave sites of tissue damage and resolve local inflammation. Further, we discuss the current tools available for investigating immune function in zebrafish and how future studies that simultaneously image multiple leukocyte subsets can be used to further dissect mechanisms that regulate both the onset and resolution of inflammation.

Transcriptomic and Metabolomic Analyses Provide Insights into the Growth and Development Advantages of Triploid Apostichopus japonicus

Polyploid breeding is widely used in aquaculture as an important area of new research. We have previously grown Apostichopus japonicus triploids with a growth advantage. The body length, body weight, and aestivation time of triploid and diploid A. japonicus were measured in this study, and the transcriptome and metabolome were used to examine the growth advantage of triploids A. japonicus. The results showed that the proportion of triploid A. japonicus with a body length of 6-12 cm and 12-18 cm was significantly higher than that of diploid A. japonicus, and triploid A. japonicus had a shorter aestivation time (39 days) than diploid (63 days). We discovered 3296 differentially expressed genes (DEGs); 13 DEGs (for example, cyclin-dependent kinase 2) related to growth advantage, immune regulation, and energy storage were screened as potential candidates. According to Gene Ontology (GO) enrichment analysis, DEGs were significantly enriched in the cytoplasm (cellular component), ATP binding process (molecular function), oxidation-reduction process (biological process), and other pathways. According to the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment data, DEGs were significantly enriched in ribosome production and other areas. We discovered 414 significant differential metabolites (SDMs), with 11 important SDMs (for example, nocodazole) linked to a growth advantage. SDMs are significantly enriched in metabolic pathways, as well as other pathways, according to the KEGG enrichment results. According to a combined transcriptome and metabolome analysis, 6 DEGs have regulatory relationships with 11 SDMs, which act on 11 metabolic pathways together. Our results further enrich the biological data of triploid A. japonicus and provide useful resources for genetic improvement of this species.

Leukocyte-Mediated Cardiac Repair after Myocardial Infarction in Non-Regenerative vs. Regenerative Systems

Innate and adaptive leukocytes rapidly mobilize to ischemic tissues after myocardial infarction in response to damage signals released from necrotic cells. Leukocytes play important roles in cardiac repair and regeneration such as inflammation initiation and resolution; the removal of dead cells and debris; the deposition of the extracellular matrix and granulation tissue; supporting angiogenesis and cardiomyocyte proliferation; and fibrotic scar generation and resolution. By organizing and comparing the present knowledge of leukocyte recruitment and function after cardiac injury in non-regenerative to regenerative systems, we propose that the leukocyte response to cardiac injury differs in non-regenerative adult mammals such as humans and mice in comparison to cardiac regenerative models such as neonatal mice and adult zebrafish. Specifically, extensive neutrophil, macrophage, and T-cell persistence contributes to a lengthy inflammatory period in non-regenerative systems for adverse cardiac remodeling and heart failure development, whereas their quick removal supports inflammation resolution in regenerative systems for new contractile tissue formation and coronary revascularization. Surprisingly, other leukocytes have not been examined in regenerative model systems. With this review, we aim to encourage the development of improved immune cell markers and tools in cardiac regenerative models for the identification of new immune targets in non-regenerative systems to develop new therapies.

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

The immune system is a complex and sophisticated biological system, spanning multiple levels of complexity, from the molecular level to that of tissue. Our current understanding of its function and complexity, of the heterogeneity of leukocytes, is a result of decades of concentrated efforts to delineate cellular markers using conventional methods of antibody screening and antigen identification. In mammalian models, this led to in-depth understanding of individual leukocyte subsets, their phenotypes, and their roles in health and disease. The field was further propelled forward by the development of single-cell (sc) RNA-seq technologies, offering an even broader and more integrated view of how cells work together to generate a particular response. Consequently, the adoption of scRNA-seq revealed the unexpected plasticity and heterogeneity of leukocyte populations and shifted several long-standing paradigms of immunology. This review article highlights the unprecedented opportunities offered by scRNA-seq technology to unveil the individual contributions of leukocyte subsets and their crosstalk in generating the overall immune responses in bony fishes. Single-cell transcriptomics allow identifying unseen relationships, and formulating novel hypotheses tailored for teleost species, without the need to rely on the limited number of fish-specific antibodies and pre-selected markers. Several recent studies on single-cell transcriptomes of fish have already identified previously unnoticed expression signatures and provided astonishing insights into the diversity of teleost leukocytes and the evolution of vertebrate immunity. Without a doubt, scRNA-seq in tandem with bioinformatics tools and state-of-the-art methods, will facilitate studying the teleost immune system by not only defining key markers, but also teaching us about lymphoid tissue organization, development/differentiation, cell-cell interactions, antigen receptor repertoires, states of health and disease, all across time and space in fishes. These advances will invite more researchers to develop the tools necessary to explore the immunology of fishes, which remain non-conventional animal models from which we have much to learn.

A Single-Cell Transcriptome Profiling of Anterior Kidney Leukocytes From Nile Tilapia (Oreochromis niloticus)

Teleost fish anterior kidney (AK) is an important hematopoietic organ with multifarious immune cells, which have immune functions comparable to mammalian bone marrow. Myeloid and lymphoid cells locate in the AK, but the lack of useful specific gene markers and antibody-based reagents for the cell subsets makes the identification of the different cell types difficult. Single-cell transcriptome sequencing enables single-cell capture and individual library construction, making the study on the immune cell heterogeneity of teleost fish AK possible. In this study, we examined the transcriptional patterns of 11,388 AK leukocytes using 10× Genomics single-cell RNA sequencing (scRNA-seq). A total of 22 clusters corresponding to five distinct immune cell subsets were identified, which included B cells, T cells, granulocytes, macrophages, and dendritic cells (DCs). However, the subsets of myeloid cells (granulocytes, macrophages, and DCs) were not identified in more detail according to the known specific markers, even though significant differences existed among the clusters. Thereafter, we highlighted the B-cell subsets and identified them as pro/pre B cells, immature/mature B cells, activated B/plasmablasts, or plasma cells based on the different expressions of the transcription factors (TFs) and cytokines. Clustering of the differentially modulated genes by pseudo-temporal trajectory analysis of the B-cell subsets showed the distinct kinetics of the responses of TFs to cell conversion. Moreover, we classified the T cells and discovered that CD3+CD4-CD8-, CD3+CD4+CD8+, CD4+CD8-, and CD4-CD8+ T cells existed in AK, but neither CD4+CD8- nor CD4-CD8+ T cells can be further classified into subsets based on the known TFs and cytokines. Pseudotemporal analysis demonstrated that CD4+CD8- and CD4-CD8+ T cells belonged to different states with various TFs that might control their differentiation. The data obtained above provide a valuable and detailed resource for uncovering the leukocyte subsets in Nile tilapia AK, as well as more potential markers for identifying the myeloid and lymphoid cell types.

Inflammation and convergent placenta gene co-option contributed to a novel reproductive tissue

The evolution of pregnancy exposes parental tissues to new, potentially stressful conditions, which can trigger inflammation.¹ Inflammation is costly^2,3 and can induce embryo rejection, which constrains the evolution of pregnancy.¹ In contrast, inflammation can also promote morphological innovation at the maternal-embryonic interface as exemplified by co-option of pro-inflammatory signaling for eutherian embryo implantation.^1,4,5 Given its dual function, inflammation could be a key process explaining how innovations such as pregnancy and placentation evolved many times convergently. Pelvic brooding ricefishes evolved a novel “plug” tissue,^6,7 which forms inside the female gonoduct after spawning, anchors egg-attaching filaments, and enables pelvic brooders to carry eggs externally until hatching.^6,8 Compared to pregnancy, i.e., internal bearing of embryos, external bearing should alleviate constraints on inflammation in the reproductive tract. We thus hypothesized that an ancestral inflammation triggered by the retention of attaching filaments gave rise to pathways orchestrating plug formation. In line with our hypothesis, histological sections of the developing plug revealed signs of gonoduct injuries by egg-attaching filaments in the pelvic brooding ricefish Oryzias eversi. Tissue-specific transcriptomes showed that inflammatory signaling dominates the plug transcriptome and inflammation-induced genes controlling vital processes for plug development such as tissue growth and angiogenesis were overexpressed in the plug. Finally, mammalian placenta genes were enriched in the plug transcriptome, indicating convergent gene co-option for building, attaching, and sustaining a transient tissue in the female reproductive tract. This study highlights the role of gene co-option and suggests that recruiting inflammatory signaling into physiological processes provides a fast-track to evolutionary innovation.

•

Pelvic brooding induces tissue-specific changes in gene expression

•

Inflammatory signaling characterizes transcriptome of the egg-anchoring plug

•

Similar to embryo implantation, the plug likely evolved from an inflammatory response

•

Mammalian placenta genes were independently co-opted into the plug

High-Resolution, 3D Imaging of the Zebrafish Gill-Associated Lymphoid Tissue (GIALT) Reveals a Novel Lymphoid Structure, the Amphibranchial Lymphoid Tissue

The zebrafish is extensively used as an animal model for human and fish diseases. However, our understanding of the structural organization of its immune system remains incomplete, especially the mucosa-associated lymphoid tissues (MALTs). Teleost MALTs are commonly perceived as diffuse and scattered populations of immune cells throughout the mucosa. Yet, structured MALTs have been recently discovered in Atlantic salmon (Salmo salar L.), including the interbranchial lymphoid tissue (ILT) in the gills. The existence of the ILT was only recently identified in zebrafish and other fish species, highlighting the need for in-depth characterizations of the gill-associated lymphoid tissue (GIALT) in teleosts. Here, using 3-D high-resolution microscopy, we analyze the GIALT of adult zebrafish with an immuno-histology approach that reveals the organization of lymphoid tissues via the labeling of T/NK cells with an antibody directed to a highly conserved epitope on the kinase ZAP70. We show that the GIALT in zebrafish is distributed over at least five distinct sub-regions, an organization found in all pairs of gill arches. The GIALT is diffuse in the pharyngeal part of the gill arch, the interbranchial septum and the filaments/lamellae, and structured in two sub-regions: the ILT, and a newly discovered lymphoid structure located along each side of the gill arch, which we named the Amphibranchial Lymphoid Tissue (ALT). Based on RAG2 expression, neither the ILT nor the ALT constitute additional thymi. The ALT shares several features with the ILT such as presence of abundant lymphoid cells and myeloid cells embedded in a network of reticulated epithelial cells. Further, the ILT and the ALT are also a site for T/NK cell proliferation. Both ILT and ALT show structural changes after infection with Spring Viraemia of Carp Virus (SVCV). Together, these data suggest that ALT and ILT play an active role in immune responses. Comparative studies show that whereas the ILT seems absent in most neoteleosts ("Percomorphs"), the ALT is widely present in cyprinids, salmonids and neoteleosts, suggesting that it constitutes a conserved tissue involved in the protection of teleosts via the gills.

Functional coordination of non-myocytes plays a key role in adult zebrafish heart regeneration

Cardiac regeneration occurs primarily through proliferation of existing cardiomyocytes, but also involves complex interactions between distinct cardiac cell types including non-cardiomyocytes (non-CMs). However, the subpopulations, distinguishing molecular features, cellular functions, and intercellular interactions of non-CMs in heart regeneration remain largely unexplored. Using the LIGER algorithm, we assemble an atlas of cell states from 61,977 individual non-CM scRNA-seq profiles isolated at multiple time points during regeneration. This analysis reveals extensive non-CM cell diversity, including multiple macrophage (MC), fibroblast (FB), and endothelial cell (EC) subpopulations with unique spatiotemporal distributions, and suggests an important role for MC in inducing the activated FB and EC subpopulations. Indeed, pharmacological perturbation of MC function compromises the induction of the unique FB and EC subpopulations. Furthermore, we developed computational algorithm Topologizer to map the topological relationships and dynamic transitions between functional states. We uncover dynamic transitions between MC functional states and identify factors involved in mRNA processing and transcriptional regulation associated with the transition. Together, our single-cell transcriptomic analysis of non-CMs during cardiac regeneration provides a blueprint for interrogating the molecular and cellular basis of this process.

Full-Length Transcriptome: A Reliable Alternative for Single-Cell RNA-Seq Analysis in the Spleen of Teleost Without Reference Genome

Fish is considered as a supreme model for clarifying the evolution and regulatory mechanism of vertebrate immunity. However, the knowledge of distinct immune cell populations in fish is still limited, and further development of techniques advancing the identification of fish immune cell populations and their functions are required. Single cell RNA-seq (scRNA-seq) has provided a new approach for effective in-depth identification and characterization of cell subpopulations. Current approaches for scRNA-seq data analysis usually rely on comparison with a reference genome and hence are not suited for samples without any reference genome, which is currently very common in fish research. Here, we present an alternative, i.e. scRNA-seq data analysis with a full-length transcriptome as a reference, and evaluate this approach on samples from Epinephelus coioides-a teleost without any published genome. We show that it reconstructs well most of the present transcripts in the scRNA-seq data achieving a sensitivity equivalent to approaches relying on genome alignments of related species. Based on cell heterogeneity and known markers, we characterized four cell types: T cells, B cells, monocytes/macrophages (Mo/MΦ) and NCC (non-specific cytotoxic cells). Further analysis indicated the presence of two subsets of Mo/MΦ including M1 and M2 type, as well as four subsets in B cells, i.e. mature B cells, immature B cells, pre B cells and early-pre B cells. Our research will provide new clues for understanding biological characteristics, development and function of immune cell populations of teleost. Furthermore, our approach provides a reliable alternative for scRNA-seq data analysis in teleost for which no reference genome is currently available.

To React or Not to React: The Dilemma of Fish Immune Systems Facing Myxozoan Infections

Myxozoans are microscopic, metazoan, obligate parasites, belonging to the phylum Cnidaria. In contrast to the free-living lifestyle of most members of this taxon, myxozoans have complex life cycles alternating between vertebrate and invertebrate hosts. Vertebrate hosts are primarily fish, although they are also reported from amphibians, reptiles, trematodes, mollusks, birds and mammals. Invertebrate hosts include annelids and bryozoans. Most myxozoans are not overtly pathogenic to fish hosts, but some are responsible for severe economic losses in fisheries and aquaculture. In both scenarios, the interaction between the parasite and the host immune system is key to explain such different outcomes of this relationship. Innate immune responses contribute to the resistance of certain fish strains and species, and the absence or low levels of some innate and regulatory factors explain the high pathogenicity of some infections. In many cases, immune evasion explains the absence of a host response and allows the parasite to proliferate covertly during the first stages of the infection. In some infections, the lack of an appropriate regulatory response results in an excessive inflammatory response, causing immunopathological consequences that are worse than inflicted by the parasite itself. This review will update the available information about the immune responses against Myxozoa, with special focus on T and B lymphocyte and immunoglobulin responses, how these immune effectors are modulated by different biotic and abiotic factors, and on the mechanisms of immune evasion targeting specific immune effectors. The current and future design of control strategies for myxozoan diseases is based on understanding this myxozoan-fish interaction, and immune-based strategies such as improvement of innate and specific factors through diets and additives, host genetic selection, passive immunization and vaccination, are starting to be considered.

Holosteans contextualize the role of the teleost genome duplication in promoting the rise of evolutionary novelties in the ray-finned fish innate immune system

Over 99% of ray-finned fishes (Actinopterygii) are teleosts, a clade that comprises half of all living vertebrate species that have diversified across virtually all fresh and saltwater ecosystems. This ecological breadth raises the question of how the immunogenetic diversity required to persist under heterogeneous pathogen pressures evolved. The teleost genome duplication (TGD) has been hypothesized as the evolutionary event that provided the substrate for rapid genomic evolution and innovation. However, studies of putative teleost-specific innate immune receptors have been largely limited to comparisons either among teleosts or between teleosts and distantly related vertebrate clades such as tetrapods. Here we describe and characterize the receptor diversity of two clustered innate immune gene families in the teleost sister lineage: Holostei (bowfin and gars). Using genomic and transcriptomic data, we provide a detailed investigation of the phylogenetic history and conserved synteny of gene clusters encoding diverse immunoglobulin domain-containing proteins (DICPs) and novel immune-type receptors (NITRs). These data demonstrate an ancient linkage of DICPs to the major histocompatibility complex (MHC) and reveal an evolutionary origin of NITR variable-joining (VJ) exons that predate the TGD by at least 50 million years. Further characterizing the receptor diversity of Holostean DICPs and NITRs illuminates a sequence diversity that rivals the diversity of these innate immune receptor families in many teleosts. Taken together, our findings provide important historical context for the evolution of these gene families that challenge prevailing expectations concerning the consequences of the TGD during actinopterygiian evolution.

Zebrafish larvae as experimental model to expedite the search for new biomarkers and treatments for neonatal sepsis

Neonatal sepsis is a major cause of death and disability in newborns. Commonly used biomarkers for diagnosis and evaluation of treatment response lack sufficient sensitivity or specificity. Additionally, new targets to treat the dysregulated immune response are needed, as are methods to effectively screen drugs for these targets. Available research methods have hitherto not yielded the breakthroughs required to significantly improve disease outcomes, we therefore describe the potential of zebrafish (Danio rerio) larvae as preclinical model for neonatal sepsis. In biomedical research, zebrafish larvae combine the complexity of a whole organism with the convenience and high-throughput potential of in vitro methods. This paper illustrates that zebrafish exhibit an immune system that is remarkably similar to humans, both in terms of types of immune cells and signaling pathways. Moreover, the developmental state of the larval immune system is highly similar to human neonates. We provide examples of zebrafish larvae being used to study infections with pathogens commonly causing neonatal sepsis and discuss known limitations. We believe this species could expedite research into immune regulation during neonatal sepsis and may hold keys for the discovery of new biomarkers and novel treatment targets as well as for screening of targeted drug therapies.

Immunological memory in teleost fish

Immunological memory can be regarded as the key aspect of adaptive immunity, i.e. a specific response to first contact with an antigen, which in mammals is determined by the properties of T, B and NK cells. Re-exposure to the same antigen results in a more rapid response of the activated specific cells, which have a unique property that is the immunological memory acquired upon first contact with the antigen. Such a state of immune activity is also to be understood as related to "altered behavior of the immune system" due to genetic alterations, presumably maintained independently of the antigen. It also indicates a possible alternative mechanism of maintaining the immune state at a low level of the immune response, "directed" by an antigen or dependent on an antigen, associated with repeated exposure to the same antigen from time to time, as well as the concept of innate immune memory, associated with epigenetic reprogramming of myeloid cells, i.e. macrophages and NK cells. Studies on Teleostei have provided evidence for the presence of immunological memory determined by T and B cells and a secondary response stronger than the primary response. Research has also demonstrated that in these animals macrophages and NK-like cells (similar to mammalian NK cells) are able to respond when re-exposed to the same antigen. Regardless of previous reports on immunological memory in teleost fish, many reactions and mechanisms related to this ability require further investigation. The very nature of immunological memory and the activity of cells involved in this process, in particular macrophages and NK-like cells, need to be explained. This paper presents problems associated with adaptive and innate immune memory in teleost fish and characteristics of cells associated with this ability.

Essential role for Gata2 in modulating lineage output from hematopoietic stem cells in zebrafish

The differentiation of hematopoietic stem cells (HSCs) is tightly controlled to ensure a proper balance between myeloid and lymphoid cell output. GATA2 is a pivotal hematopoietic transcription factor required for generation and maintenance of HSCs. GATA2 is expressed throughout development, but because of early embryonic lethality in mice, its role during adult hematopoiesis is incompletely understood. Zebrafish contains 2 orthologs of GATA2: Gata2a and Gata2b, which are expressed in different cell types. We show that the mammalian functions of GATA2 are split between these orthologs. Gata2b-deficient zebrafish have a reduction in embryonic definitive hematopoietic stem and progenitor cell (HSPC) numbers, but are viable. This allows us to uniquely study the role of GATA2 in adult hematopoiesis. gata2b mutants have impaired myeloid lineage differentiation. Interestingly, this defect arises not in granulocyte-monocyte progenitors, but in HSPCs. Gata2b-deficient HSPCs showed impaired progression of the myeloid transcriptional program, concomitant with increased coexpression of lymphoid genes. This resulted in a decrease in myeloid-programmed progenitors and a relative increase in lymphoid-programmed progenitors. This shift in the lineage output could function as an escape mechanism to avoid a block in lineage differentiation. Our study helps to deconstruct the functions of GATA2 during hematopoiesis and shows that lineage differentiation flows toward a lymphoid lineage in the absence of Gata2b.

Single-cell RNA sequencing in cancer: Applications, advances, and emerging challenges

Cancer has become one of the greatest threats to human health, and new technologies are urgently needed to further clarify the mechanisms of cancer so that better detection and treatment strategies can be developed. At present, extensive genomic analysis and testing of clinical specimens shape the insights into carcinoma. Nevertheless, carcinoma of humans is a complex ecosystem of cells, including carcinoma cells and immunity-related and stroma-related subsets, with accurate characteristics obscured by extensive genome-related approaches. A growing body of research shows that sequencing of single-cell RNA (scRNA-seq) is emerging to be an effective way for dissecting human tumor tissue at single-cell resolution, presenting one prominent way for explaining carcinoma biology. This review summarizes the research progress of scRNA-seq in the field of tumors, focusing on the application of scRNA-seq in tumor circulating cells, tumor stem cells, tumor drug resistance, the tumor microenvironment, and so on, which provides a new perspective for tumor research.

Single-cell RNA-seq analysis reveals penaeid shrimp hemocyte subpopulations and cell differentiation process

Crustacean aquaculture is expected to be a major source of fishery commodities in the near future. Hemocytes are key players of the immune system in shrimps; however, their classification, maturation, and differentiation are still under debate. To date, only discrete and inconsistent information on the classification of shrimp hemocytes has been reported, showing that the morphological characteristics are not sufficient to resolve their actual roles. Our present study using single-cell RNA sequencing revealed six types of hemocytes of Marsupenaeus japonicus based on their transcriptional profiles. We identified markers of each subpopulation and predicted the differentiation pathways involved in their maturation. We also predicted cell growth factors that might play crucial roles in hemocyte differentiation. Different immune roles among these subpopulations were suggested from the analysis of differentially expressed immune-related genes. These results provide a unified classification of shrimp hemocytes, which improves the understanding of its immune system.