Single-cell transcriptional analysis of normal, aberrant, and malignant hematopoiesis in zebrafish

Stat3 Regulates Developmental Hematopoiesis and Impacts Myeloid Cell Function via Canonical and Non-Canonical Modalities

INTRODUCTION

Signal transducer and activator of transcription (STAT) 3 is extensively involved in the development, homeostasis, and function of immune cells, with STAT3 disruption associated with human immune-related disorders. The roles ascribed to STAT3 have been assumed to be due to its canonical mode of action as an inducible transcription factor downstream of multiple cytokines, although alternative noncanonical functional modalities have also been identified. The relative involvement of each mode was further explored in relevant zebrafish models.

METHODS

Genome editing with CRISPR/Cas9 was used to generate mutants of the conserved zebrafish Stat3 protein: a loss of function knockout (KO) mutant and a mutant lacking C-terminal sequences including the transactivation domain (ΔTAD). Lines harboring these mutations were analyzed with respect to blood and immune cell development and function in comparison to wild-type zebrafish.

RESULTS

The Stat3 KO mutant showed perturbation of hematopoietic lineages throughout primitive and early definitive hematopoiesis. Neutrophil numbers did not increase in response to lipopolysaccharide (LPS) or granulocyte colony-stimulating factor (G-CSF) and their migration was significantly diminished, the latter correlating with abrogation of the Cxcl8b/Cxcr2 pathway, with macrophage responses perturbed. Intriguingly, many of these phenotypes were not shared by the Stat3 ΔTAD mutant. Indeed, only neutrophil and macrophage development were disrupted in these mutants with responsiveness to LPS and G-CSF maintained, and neutrophil migration actually increased.

CONCLUSION

This study has identified roles for zebrafish Stat3 within hematopoietic stem cells impacting multiple lineages throughout primitive and early definitive hematopoiesis, myeloid cell responses to G-CSF and LPS and neutrophil migration. Many of these roles showed conservation, but notably several involved noncanonical modalities, providing additional insights for relevant diseases.

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.

Janus Kinase 3 (JAK3): A Critical Conserved Node in Immunity Disrupted in Immune Cell Cancer and Immunodeficiency

The Janus kinase (JAK) family is a small group of protein tyrosine kinases that represent a central component of intracellular signaling downstream from a myriad of cytokine receptors. The JAK3 family member performs a particularly important role in facilitating signal transduction for a key set of cytokine receptors that are essential for immune cell development and function. Mutations that impact JAK3 activity have been identified in a number of human diseases, including somatic gain-of-function (GOF) mutations associated with immune cell malignancies and germline loss-of-function (LOF) mutations associated with immunodeficiency. The structure, function and impacts of both GOF and LOF mutations of JAK3 are highly conserved, making animal models highly informative. This review details the biology of JAK3 and the impact of its perturbation in immune cell-related diseases, including relevant animal studies.

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.

Lineage skewing and genome instability underlie marrow failure in a zebrafish model of GATA2 deficiency

Inherited bone marrow failure associated with heterozygous mutations in GATA2 predisposes toward hematological malignancies, but the mechanisms remain poorly understood. Here, we investigate the mechanistic basis of marrow failure in a zebrafish model of GATA2 deficiency. Single-cell transcriptomics and chromatin accessibility assays reveal that loss of gata2a leads to skewing toward the erythroid lineage at the expense of myeloid cells, associated with loss of cebpa expression and decreased PU.1 and CEBPA transcription factor accessibility in hematopoietic stem and progenitor cells (HSPCs). Furthermore, gata2a mutants show impaired expression of npm1a, the zebrafish NPM1 ortholog. Progressive loss of npm1a in HSPCs is associated with elevated levels of DNA damage in gata2a mutants. Thus, Gata2a maintains myeloid lineage priming through cebpa and protects against genome instability and marrow failure by maintaining expression of npm1a. Our results establish a potential mechanism underlying bone marrow failure in GATA2 deficiency.

Single-cell transcriptome analysis reveals a cellular immune response in freshwater dark sleeper (Odontobutis potamophila) after infection with Aeromonas veronii

The bacterium Aeromonas veronii is a co-pathogenic species that can negatively impact the health of both humans and aquatic animals. In this study, we used single-cell transcriptome analysis (scRNA-seq) to investigate the effects of infection with A. veronii on head kidney cells and the regulation of gene expression in the dark sleeper (Odontobutis potamophila). scRNA-seq was used to assess the effects of infection with A. veronii in O. potamophila B cells, endothelial cells, macrophages, and granulocytes, and differential enrichment analysis of gene expression in B cells and granulocytes was performed. The analyses revealed a significant increase in neutrophils and decrease in eosinophils in granulocytes infected with A. veronii. Activation of neutrophils enhanced ribosome biogenesis by up-regulating the expression of RPS12 and RPL12 to fight against invading pathogens. Crucial pro-inflammatory mediators IL1B, IGHV1-4, and the major histocompatibility class II genes MHC2A and MHC2DAB, which are involved in virulence processes, were upregulated, suggesting that A. veronii activates an immune response that presents antigens and activates immunoglobulin receptors in B cells. These cellular immune responses triggered by infection with A. veronii enriched the available scRNA-seq data for teleosts, and these results are important for understanding the evolution of cellular immune defense and functional differentiation of head kidney cells.

A zebrafish model of growth hormone insensitivity syndrome with immune dysregulation 1 (GHISID1)

Signal transducer and activator of transcription (STAT) proteins act downstream of cytokine receptors to facilitate changes in gene expression that impact a range of developmental and homeostatic processes. Patients harbouring loss-of-function (LOF) STAT5B mutations exhibit postnatal growth failure due to lack of responsiveness to growth hormone as well as immune perturbation, a disorder called growth hormone insensitivity syndrome with immune dysregulation 1 (GHISID1). This study aimed to generate a zebrafish model of this disease by targeting the stat5.1 gene using CRISPR/Cas9 and characterising the effects on growth and immunity. The zebrafish Stat5.1 mutants were smaller, but exhibited increased adiposity, with concomitant dysregulation of growth and lipid metabolism genes. The mutants also displayed impaired lymphopoiesis with reduced T cells throughout the lifespan, along with broader disruption of the lymphoid compartment in adulthood, including evidence of T cell activation. Collectively, these findings confirm that zebrafish Stat5.1 mutants mimic the clinical impacts of human STAT5B LOF mutations, establishing them as a model of GHISID1.

Zebrafish: A Relevant Genetic Model for Human Primary Immunodeficiency (PID) Disorders?

Primary immunodeficiency (PID) disorders, also commonly referred to as inborn errors of immunity, are a heterogenous group of human genetic diseases characterized by defects in immune cell development and/or function. Since these disorders are generally uncommon and occur on a variable background profile of potential genetic and environmental modifiers, animal models are critical to provide mechanistic insights as well as to create platforms to underpin therapeutic development. This review aims to review the relevance of zebrafish as an alternative genetic model for PIDs. It provides an overview of the conservation of the zebrafish immune system and details specific examples of zebrafish models for a multitude of specific human PIDs across a range of distinct categories, including severe combined immunodeficiency (SCID), combined immunodeficiency (CID), multi-system immunodeficiency, autoinflammatory disorders, neutropenia and defects in leucocyte mobility and respiratory burst. It also describes some of the diverse applications of these models, particularly in the fields of microbiology, immunology, regenerative biology and oncology.

Transgenic IDH2R172K and IDH2R140Q zebrafish models recapitulated features of human acute myeloid leukemia

Isocitrate dehydrogenase 2 (IDH2) mutations occur in more than 15% of cytogenetically normal acute myeloid leukemia (CN-AML) but comparative studies of their roles in leukemogenesis have been scarce. We generated zebrafish models of IDH2^R172K and IDH2^R140Q AML and reported their pathologic, functional and transcriptomic features and therapeutic responses to target therapies. Transgenic embryos co-expressing FLT3^ITD and IDH2 mutations showed accentuation of myelopoiesis. As these embryos were raised to adulthood, full-blown leukemia ensued with multi-lineage dysplasia, increase in myeloblasts and marrow cellularity and splenomegaly. The leukemia cells were transplantable into primary and secondary recipients and resulted in more aggressive disease. Tg(Runx1:FLT3^ITDIDH2^R172K) but not Tg(Runx1:FLT3^ITDIDH2^R140Q) zebrafish showed an increase in T-cell development at embryonic and adult stages. Single-cell transcriptomic analysis revealed increased myeloid skewing, differentiation blockade and enrichment of leukemia-associated gene signatures in both zebrafish models. Tg(Runx1:FLT3^ITDIDH2^R172K) but not Tg(Runx1:FLT3^ITDIDH2^R140Q) zebrafish showed an increase in interferon signals at the adult stage. Leukemic phenotypes in both zebrafish could be ameliorated by quizartinib and enasidenib. In conclusion, the zebrafish models of IDH2 mutated AML recapitulated the morphologic, clinical, functional and transcriptomic characteristics of human diseases, and provided the prototype for developing zebrafish leukemia models of other genotypes that would become a platform for high throughput drug screening.

Socs3b regulates the development and function of innate immune cells in zebrafish

Introduction

Suppressor of cytokine signaling 3 (SOCS3) is a critical component of the negative feedback regulation that controls signaling by cytokines and other factors thereby ensuring that important processes such as hematopoiesis and inflammation occur at appropriate levels.

Methods

To gain further insights into SOCS3 function, the zebrafish socs3b gene was investigated through analysis of a knockout line generated using CRISPR/Cas9-mediated genome editing.

Results

Zebrafish socs3b knockout embryos displayed elevated numbers of neutrophils during primitive and definitive hematopoiesis but macrophage numbers were not altered. However, the absence of socs3b reduced neutrophil functionality but enhanced macrophage responses. Adult socs3b knockout zebrafish displayed reduced survival that correlated with an eye pathology involving extensive infiltration of neutrophils and macrophages along with immune cell dysregulation in other tissues.

Discussion

These findings identify a conserved role for Socs3b in the regulation of neutrophil production and macrophage activation.

Zebrafish Model of Severe Combined Immunodeficiency (SCID) Due to JAK3 Mutation

JAK3 is principally activated by members of the interleukin-2 receptor family and plays an essential role in lymphoid development, with inactivating JAK3 mutations causing autosomal-recessive severe combined immunodeficiency (SCID). This study aimed to generate an equivalent zebrafish model of SCID and to characterize the model across the life-course. Genome editing of zebrafish jak3 created mutants similar to those observed in human SCID. Homozygous jak3 mutants showed reduced embryonic T lymphopoiesis that continued through the larval stage and into adulthood, with B cell maturation and adult NK cells also reduced and neutrophils impacted. Mutant fish were susceptible to lymphoid leukemia. This model has many of the hallmarks of human SCID resulting from inactivating JAK3 mutations and will be useful for a variety of pre-clinical applications.

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.

Cytokine Receptor-Like Factor 3 (CRLF3) Contributes to Early Zebrafish Hematopoiesis

Cytokine receptor-like factor 3 (CRLF3) is an ancient protein conserved across metazoans that contains an archetypal cytokine receptor homology domain (CHD). This domain is found in cytokine receptors present in bilateria, including higher vertebrates, that play key roles in a variety of developmental and homeostatic processes, particularly relating to blood and immune cells. However, understanding of CRLF3 itself remains very limited. This study aimed to investigate this evolutionarily significant protein by studying its embryonic expression and function in early development, particularly of blood and immune cells, using zebrafish as a model. Expression of crlf3 was identified in mesoderm-derived tissues in early zebrafish embryos, including the somitic mesoderm and both anterior and posterior lateral plate mesoderm. Later expression was observed in the thymus, brain, retina and exocrine pancreas. Zebrafish crlf3 mutants generated by genome editing technology exhibited a significant reduction in primitive hematopoiesis and early definitive hematopoiesis, with decreased early progenitors impacting on multiple lineages. No other obvious phenotypes were observed in the crlf3 mutants.

In vivo impact of JAK3 A573V mutation revealed using zebrafish

Background

Janus kinase 3 (JAK3) acts downstream of the interleukin-2 (IL-2) receptor family to play a pivotal role in the regulation of lymphoid cell development. Activating JAK3 mutations are associated with a number of lymphoid and other malignancies, with mutations within the regulatory pseudokinase domain common.

Methods

The pseudokinase domain mutations A572V and A573V were separately introduced into the highly conserved zebrafish Jak3 and transiently expressed in cell lines and zebrafish embryos to examine their activity and impact on early T cells. Genome editing was subsequently used to introduce the A573V mutation into the zebrafish genome to study the effects of JAK3 activation on lymphoid cells in a physiologically relevant context throughout the life-course.

Results

Zebrafish Jak3 A573V produced the strongest activation of downstream STAT5 in vitro and elicited a significant increase in T cells in zebrafish embryos. Zebrafish carrying just a single copy of the Jak3 A573V allele displayed elevated embryonic T cells, which continued into adulthood. Hematopoietic precursors and NK cells were also increased, but not B cells. The lymphoproliferative effects of Jak3 A573V in embryos was shown to be dependent on zebrafish IL-2Rγc, JAK1 and STAT5B equivalents, and could be suppressed with the JAK3 inhibitor Tofacitinib.

Conclusions

This study demonstrates that a single JAK3 A573V allele expressed from the endogenous locus was able to enhance lymphopoiesis throughout the life-course, which was mediated via an IL-2Rγc/JAK1/JAK3/STAT5 signaling pathway and was sensitive to Tofacitinib. This extends our understanding of oncogenic JAK3 mutations and creates a novel model to underpin further translational investigations.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00018-022-04361-8.

Dietary Inclusion of Seabuckthorn (Hippophae rhamnoides) Mitigates Foodborne Enteritis in Zebrafish Through the Gut-Liver Immune Axis

To help prevent foodborne enteritis in aquaculture, several feed additives, such as herbal medicine, have been added to fish diets. Predictions of effective herb medicines for treating fish foodborne enteritis from key regulated DEGs (differentially expressed genes) in transcriptomic data can aid in the development of feed additives using the Traditional Chinese Medicine Integrated Database. Seabuckthorn has been assessed as a promising candidate for treating grass carp soybean-induced enteritis (SBMIE). In the present study, the SBMIE zebrafish model was used to assess seabuckthorn's therapeutic or preventative effects. The results showed that intestinal and hepatic inflammation was reduced when seabuckthorn was added, either pathologically (improved intestinal villi morphology, less oil-drops) or growth-related (body fat deposition). Moreover, seabuckthorn may block the intestinal p53 signaling pathway, while activating the PPAR signaling pathway and fatty acid metabolism in the liver. 16S rRNA gene sequencing results also indicated a significant increase in OTU numbers and skewed overlapping with the fish meal group following the addition of seabuckthorn. Additionally, there were signs of altered gut microbiota taxa composition, particularly for reduced TM7, Sphingomonas, and Shigella, following the addition of seabuckthorn. Hindgut imaging of fluorescent immune cells in SBMIE larvae revealed the immune regulatory mechanisms at the cellular level. Seabuckthorn may significantly inhibit the inflammatory gathering of neutrophils, macrophages, and mature T cells, as well as cellular protrusions' formation. On the other hand, in larvae, seabuckthorn inhibited the inflammatory aggregation of lck+ T cells but not immature lymphocytes, indicating that it affected intestinal adaptive immunity. Although seabuckthorn did not affect the distribution of intestinal CD4+ cells, the number of hepatic CD4+ cells were reduced in fish from the seabuckthorn supplementation group. Thus, the current data indicate that seabuckthorn may alleviate foodborne gut-liver symptoms by enhancing intestinal mucosal immunity and microbiota while simultaneously inhibiting hepatic adipose disposition, making it a potential additive for preventing fish foodborne gut-liver symptoms.

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.

Generation and Characterization of a Zebrafish IL-2Rγc SCID Model

The IL-2 family of cytokines act via receptor complexes that share the interleukin-2 receptor gamma common (IL-2Rγc) chain to play key roles in lymphopoiesis. Inactivating IL-2Rγc mutations results in severe combined immunodeficiency (SCID) in humans and other species. This study sought to generate an equivalent zebrafish SCID model. The zebrafish il2rga gene was targeted for genome editing using TALENs and presumed loss-of-function alleles analyzed with respect to immune cell development and impacts on intestinal microbiota and tumor immunity. Knockout of zebrafish Il-2rγc.a resulted in a SCID phenotype, including a significant reduction in T cells, with NK cells also impacted. This resulted in dysregulated intestinal microbiota and defective immunity to tumor xenotransplants. Collectively, this establishes a useful zebrafish SCID model.

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.

Shifting the focus of zebrafish toward a model of the tumor microenvironment

Cancer cells exist in a complex ecosystem with numerous other cell types in the tumor microenvironment (TME). The composition of this tumor/TME ecosystem will vary at each anatomic site and affects phenotypes such as initiation, metastasis, and drug resistance. A mechanistic understanding of the large number of cell-cell interactions between tumor and TME requires models that allow us to both characterize as well as genetically perturb this complexity. Zebrafish are a model system optimized for this problem, because of the large number of existing cell-type-specific drivers that can label nearly any cell in the TME. These include stromal cells, immune cells, and tissue resident normal cells. These cell-type-specific promoters/enhancers can be used to drive fluorophores to facilitate imaging and also CRISPR cassettes to facilitate perturbations. A major advantage of the zebrafish is the ease by which large numbers of TME cell types can be studied at once, within the same animal. While these features make the zebrafish well suited to investigate the TME, the model has important limitations, which we also discuss. In this review, we describe the existing toolset for studying the TME using zebrafish models of cancer and highlight unique biological insights that can be gained by leveraging this powerful resource.

Resistance to inflammation underlies enhanced fitness in clonal hematopoiesis

[Figure: see text].

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.

Single-cell ATAC-seq reveals GATA2-dependent priming defect in myeloid and a maturation bottleneck in lymphoid lineages

Germline heterozygous mutations in GATA2 are associated with a syndrome characterized by cytopenias, atypical infections, and increased risk of hematologic malignancies. Here, we generated a zebrafish mutant of gata2b that recapitulated the myelomonocytopenia and B-cell lymphopenia of GATA2 deficiency syndrome. Using single-cell assay for transposase accessible chromatin with sequencing of marrow cells, we showed that loss of gata2b led to contrasting alterations in chromosome accessibility in early myeloid and lymphoid progenitors, associated with defects in gene expression. Within the myeloid lineage in gata2b mutant zebrafish, we identified an attenuated myeloid differentiation with reduced transcriptional priming and skewing away from the monocytic program. In contrast, in early lymphoid progenitors, gata2b loss led to accumulation of B-lymphoid transcription factor accessibility coupled with increased expression of the B-cell lineage-specification program. However, gata2b mutant zebrafish had incomplete B-cell lymphopoiesis with loss of lineage-specific transcription factor accessibility in differentiating B cells, in the context of aberrantly reduced oxidative metabolic pathways. Our results establish that transcriptional events in early progenitors driven by Gata2 are required to complete normal differentiation.

Tipping the Scales With Zebrafish to Understand Adaptive Tumor Immunity

The future of improved immunotherapy against cancer depends on an in-depth understanding of the dynamic interactions between the immune system and tumors. Over the past two decades, the zebrafish has served as a valuable model system to provide fresh insights into both the development of the immune system and the etiologies of many different cancers. This well-established foundation of knowledge combined with the imaging and genetic capacities of the zebrafish provides a new frontier in cancer immunology research. In this review, we provide an overview of the development of the zebrafish immune system along with a side-by-side comparison of its human counterpart. We then introduce components of the adaptive immune system with a focus on their roles in the tumor microenvironment (TME) of teleosts. In addition, we summarize zebrafish models developed for the study of cancer and adaptive immunity along with other available tools and technology afforded by this experimental system. Finally, we discuss some recent research conducted using the zebrafish to investigate adaptive immune cell-tumor interactions. Without a doubt, the zebrafish will arise as one of the driving forces to help expand the knowledge of tumor immunity and facilitate the development of improved anti-cancer immunotherapy in the foreseeable future.

New tools for 'ZEBRA-FISHING'

Zebrafish has been established as a classical model for developmental studies, yet in the past years, with the explosion of novel technological methods, the use of zebrafish as a model has expanded. One of the prominent fields that took advantage of zebrafish as a model organism early on is hematopoiesis, the process of blood cell generation from hematopoietic stem and progenitor cells (HSPCs). In zebrafish, HSPCs are born early during development in the aorta-gonad-mesonephros region and then translocate to the caudal hematopoietic tissue, where they expand and finally take residence in the kidney marrow. This journey is tightly regulated at multiple levels from extracellular signals to chromatin. In order to delineate the mechanistic underpinnings of this process, next-generation sequencing techniques could be an important ally. Here, we describe genome-wide approaches that have been undertaken to delineate zebrafish hematopoiesis.

The chromatin remodeler Brg1 is required for formation and maintenance of hematopoietic stem cells

Hematopoietic stem and progenitor cells (HSPCs) have the ability to self-renew and differentiate into various blood cells, thus playing an important role in maintenance of lifelong hematopoiesis. Brahma-related gene 1 (BRG1), which acts as the ATP subunit of mammalian SWI-SNF-related chromatin remodeling complexes, is involved in human acute myeloid leukemia and highly expresses in short-term HSPCs. But its role and regulatory mechanism for HSPC development have not yet been well established. Here, we generated a brg1 knockout zebrafish model using TALEN technology. We found that in brg1^-/- embryo, the primitive hematopoiesis remained well, while definitive hematopoiesis formation was significantly impaired. The number of hemogenic endothelial cells was decreased, further affecting definitive hematopoiesis with reduced myeloid and lymphoid cells. During embryogenesis, the nitric oxide (NO) microenvironment in brg1^-/- embryo was seriously damaged and the reduction of HSPCs could be partially rescued by a NO donor. Chromatin immunoprecipitation (ChIP) assays showed that BRG1 could bind to the promoter of KLF2 and trigger its transcriptional activity of NO synthase. Our findings show that Brg1 promotes klf2a expression in hemogenic endothelium and highlight a novel mechanism for HSPC formation and maintenance.

Single-cell RNA-seq unveils critical regulators of human FOXP3+ regulatory T cell stability

The heterogeneity and plasticity of T lymphocytes is critical for determining immune response outcomes. Functional regulatory T (Treg) cells are commonly characterized by stable FOXP3 expression and have reported to exhibit heterogeneous phenotypes under inflammatory conditions. However, the interplay between inflammation and Treg cell suppressive activity still remains elusive. Here, we utilized single-cell RNA sequencing to investigate how human Treg cells respond to the pro-inflammatory cytokine interleukin-6 (IL-6). We observed that Treg cells divided into two subpopulations after IL-6 stimulation. TIGIT^- unstable Treg cells lost FOXP3 expression and gained an effector-like T cell phenotype, whereas TIGIT⁺ Treg cells retained robust suppressive function. Single cell transcriptome analysis revealed a spectrum of cellular states of IL-6-stimulated Treg cells and how cytochrome P450 family 1 subfamily A member 1 (CYP1A1) is a crucial regulator of Treg cell suppressive capability and stability. CYP1A1-deficient human Treg cells developed a Th17-like phenotype after IL-6 stimulation. Our findings implicate CYP1A1 as a previously unidentified regulator of Treg cells that may have target potential for clinical application for biotherapies.

Cell Type Purification by Single-Cell Transcriptome-Trained Sorting

Much of current molecular and cell biology research relies on the ability to purify cell types by fluorescence-activated cell sorting (FACS). FACS typically relies on the ability to label cell types of interest with antibodies or fluorescent transgenic constructs. However, antibody availability is often limited, and genetic manipulation is labor intensive or impossible in the case of primary human tissue. To date, no systematic method exists to enrich for cell types without a priori knowledge of cell-type markers. Here, we propose GateID, a computational method that combines single-cell transcriptomics with FACS index sorting to purify cell types of choice using only native cellular properties such as cell size, granularity, and mitochondrial content. We validate GateID by purifying various cell types from zebrafish kidney marrow and the human pancreas to high purity without resorting to specific antibodies or transgenes.

Zebrafish B cell acute lymphoblastic leukemia: new findings in an old model

Acute lymphoblastic leukemia (ALL) is the most common pediatric, and ninth most common adult, cancer. ALL can develop in either B or T lymphocytes, but B-lineage ALL (B-ALL) exceeds T-ALL clinically. As for other cancers, animal models allow study of the molecular mechanisms driving ALL. Several zebrafish (Danio rerio) T-ALL models have been reported, but until recently, robust D. rerio B-ALL models were not described. Then, D. rerio B-ALL was discovered in two related zebrafish transgenic lines; both were already known to develop T-ALL. Here, we report new B-ALL findings in one of these models, fish expressing transgenic human MYC (hMYC). We describe B-ALL incidence in a large cohort of hMYC fish, and show B-ALL in two new lines where T-ALL does not interfere with B-ALL detection. We also demonstrate B-ALL responses to steroid and radiation treatments, which effect ALL remissions, but are usually followed by prompt relapses. Finally, we report gene expression in zebrafish B lymphocytes and B-ALL, in both bulk samples and single B- and T-ALL cells. Using these gene expression profiles, we compare differences between the two new D. rerio B-ALL models, which are both driven by transgenic mammalian MYC oncoproteins. Collectively, these new data expand the utility of this new vertebrate B-ALL model.

Single-cell qPCR demonstrates that Repsox treatment changes cell fate from endoderm to neuroectoderm and disrupts epithelial-mesenchymal transition

A definitive endodermal cell lineage is a prerequisite for the efficient generation of mature endoderm derivatives that give rise to organs, such as the pancreas and liver. We previously reported that the induction of mesenchymal definitive endoderm cells depends on autocrine TGF-β signaling and that pharmacological blockage of TGF-β signaling by Repsox disrupts endoderm specification. The definitive endoderm arises from a primitive streak, which depends largely on TGF-β signaling. If the TGF-β pathway is blocked by Repsox, cell fate after the primitive streak induction is so-far unknown. We report here, that an induced primitive streak cell-population contained many T/SOX2 co-expressing cells, and subsequent inhibition of TGF-β signaling by Repsox promoted neuroectodermal cell fate, which was characterized using single-cell qPCR analysis and immunostaining. The process of epithelial-to-mesenchymal transition, which is inherent to the process of definitive endoderm differentiation, was also disrupted upon Repsox treatment. Our findings may provide a new approach to produce neural progenitors.

Visualizing Engrafted Human Cancer and Therapy Responses in Immunodeficient Zebrafish

Xenograft cell transplantation into immunodeficient mice has become the gold standard for assessing pre-clinical efficacy of cancer drugs, yet direct visualization of single-cell phenotypes is difficult. Here, we report an optically-clear prkdc-/-, il2rga-/- zebrafish that lacks adaptive and natural killer immune cells, can engraft a wide array of human cancers at 37°C, and permits the dynamic visualization of single engrafted cells. For example, photoconversion cell-lineage tracing identified migratory and proliferative cell states in human rhabdomyosarcoma, a pediatric cancer of muscle. Additional experiments identified the preclinical efficacy of combination olaparib PARP inhibitor and temozolomide DNA-damaging agent as an effective therapy for rhabdomyosarcoma and visualized therapeutic responses using a four-color FUCCI cell-cycle fluorescent reporter. These experiments identified that combination treatment arrested rhabdomyosarcoma cells in the G2 cell cycle prior to induction of apoptosis. Finally, patient-derived xenografts could be engrafted into our model, opening new avenues for developing personalized therapeutic approaches in the future.

Single-cell trajectories reconstruction, exploration and mapping of omics data with STREAM

Single-cell transcriptomic assays have enabled the de novo reconstruction of lineage differentiation trajectories, along with the characterization of cellular heterogeneity and state transitions. Several methods have been developed for reconstructing developmental trajectories from single-cell transcriptomic data, but efforts on analyzing single-cell epigenomic data and on trajectory visualization remain limited. Here we present STREAM, an interactive pipeline capable of disentangling and visualizing complex branching trajectories from both single-cell transcriptomic and epigenomic data. We have tested STREAM on several synthetic and real datasets generated with different single-cell technologies. We further demonstrate its utility for understanding myoblast differentiation and disentangling known heterogeneity in hematopoiesis for different organisms. STREAM is an open-source software package.

Purification of high-quality RNA from a small number of fluorescence activated cell sorted zebrafish cells for RNA sequencing purposes

Background

Transgenic zebrafish lines with the expression of a fluorescent reporter under the control of a cell-type specific promoter, enable transcriptome analysis of FACS sorted cell populations. RNA quality and yield are key determinant factors for accurate expression profiling. Limited cell number and FACS induced cellular stress make RNA isolation of sorted zebrafish cells a delicate process. We aimed to optimize a workflow to extract sufficient amounts of high-quality RNA from a limited number of FACS sorted cells from Tg(fli1a:GFP) zebrafish embryos, which can be used for accurate gene expression analysis.

Results

We evaluated two suitable RNA isolation kits (the RNAqueous micro and the RNeasy plus micro kit) and determined that sorting cells directly into lysis buffer is a critical step for success. For low cell numbers, this ensures direct cell lysis, protects RNA from degradation and results in a higher RNA quality and yield. We showed that this works well up to 0.5× dilution of the lysis buffer with sorted cells. In our sort settings, this corresponded to 30,000 and 75,000 cells for the RNAqueous micro kit and RNeasy plus micro kit respectively. Sorting more cells dilutes the lysis buffer too much and requires the use of a collection buffer. We also demonstrated that an additional genomic DNA removal step after RNA isolation is required to completely clear the RNA from any contaminating genomic DNA. For cDNA synthesis and library preparation, we combined SmartSeq v4 full length cDNA library amplification, Nextera XT tagmentation and sample barcoding. Using this workflow, we were able to generate highly reproducible RNA sequencing results.

Conclusions

The presented optimized workflow enables to generate high quality RNA and allows accurate transcriptome profiling of small populations of sorted zebrafish cells.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5608-2) contains supplementary material, which is available to authorized users.

Technical Advances in Single-Cell RNA Sequencing and Applications in Normal and Malignant Hematopoiesis

Single-cell RNA sequencing (scRNA-seq) has been tremendously developed in the past decade owing to overcoming challenges associated with isolation of massive quantities of single cells. Previously, cell heterogeneity and low quantities of available biological material posed significant difficulties to scRNA-seq. Cell-to-cell variation and heterogeneity are fundamental and intrinsic characteristics of normal and malignant hematopoietic cells; this heterogeneity has often been ignored in omics studies. The application of scRNA-seq has profoundly changed our comprehension of many biological phenomena, including organ development and carcinogenesis. Hematopoiesis, is actually a maturation process for more than ten distinct blood and immune cells, and is thought to be critically involved in hematological homeostasis and in sustaining the physiological functions. However, aberrant hematopoiesis directly leads to hematological malignancy, and a deeper understanding of malignant hematopoiesis will provide deeper insights into diagnosis and prognosis for patients with hematological malignancies. Here, we aim to review the recent technical progress and future prospects for scRNA-seq, as applied in physiological and malignant hematopoiesis, in efforts to further understand the hematopoietic hierarchy and to illuminate personalized therapy and precision medicine approaches used in the clinical treatment of hematological malignancies.

Genetic Models of Leukemia in Zebrafish

The zebrafish animal model is gaining increasing popularity as a tool for studying human disease. Over the past 15 years, many models of leukemia and other hematological malignancies have been developed in the zebrafish. These confer some significant advantages over similar models in other animals and systems, representing a powerful resource for investigation of the molecular basis of human leukemia. This review discusses the various zebrafish models of lymphoid and myeloid leukemia available, the major discoveries that have been made possible by them, and opportunities for future exploration.

Simultaneous B and T cell acute lymphoblastic leukemias in zebrafish driven by transgenic MYC: implications for oncogenesis and lymphopoiesis

Precursor-B cell acute lymphoblastic leukemia (pre-B ALL) is the most common pediatric cancer, but there are no useful zebrafish pre-B ALL models. We describe the first highly- penetrant zebrafish pre-B ALL, driven by human MYC. Leukemias express B lymphoblast-specific genes and are distinct from T cell ALL (T-ALL)—which these fish also develop. Zebrafish pre-B ALL shares in vivo features and expression profiles with human pre-B ALL, and these profiles differ from zebrafish T-ALL or normal B and T cells. These animals also exhibit aberrant lymphocyte development. As the only robust zebrafish pre-B ALL model and only example where T-ALL also develops, this model can reveal differences between MYC-driven pre-B vs. T-ALL and be exploited to discover novel pre-B ALL therapies.

Organism-Level Tumor Models in Zebrafish Danio rerio

Development and implementation of adequate organism-level models is one of the key elements in biomedical research that focuses on experimental oncology. Over the last decade, studies using Zebrafish (Danio rerio) have gained in popularity in this area of research. This review describes the various approaches that have been used in developing highly effective models for oncological (clinical term, better cancer or tumor) studies based on D. rerio. Priority is given to transplantation models of cancer and their application to optically transparent D. rerio lines, including clonal ones, and utilization tumors of various origins bearing fluorescent labels. The combination of tumor transplantation at organism-level models in transparent clonal D. rerio lines with fluorescent microscopy, FACS-fractionation of tumor cell subsets, and transcription analysis can result in one of the most promising research approaches in providing new information on tumor formation and growth.

Whole-organism clone tracing using single-cell sequencing

Embryonic development is a crucial period in the life of a multicellular organism, during which limited sets of embryonic progenitors produce all cells in the adult body. Determining which fate these progenitors acquire in adult tissues requires the simultaneous measurement of clonal history and cell identity at single-cell resolution, which has been a major challenge. Clonal history has traditionally been investigated by microscopically tracking cells during development, monitoring the heritable expression of genetically encoded fluorescent proteins and, more recently, using next-generation sequencing technologies that exploit somatic mutations, microsatellite instability, transposon tagging, viral barcoding, CRISPR-Cas9 genome editing and Cre-loxP recombination. Single-cell transcriptomics provides a powerful platform for unbiased cell-type classification. Here we present ScarTrace, a single-cell sequencing strategy that enables the simultaneous quantification of clonal history and cell type for thousands of cells obtained from different organs of the adult zebrafish. Using ScarTrace, we show that a small set of multipotent embryonic progenitors generate all haematopoietic cells in the kidney marrow, and that many progenitors produce specific cell types in the eyes and brain. In addition, we study when embryonic progenitors commit to the left or right eye. ScarTrace reveals that epidermal and mesenchymal cells in the caudal fin arise from the same progenitors, and that osteoblast-restricted precursors can produce mesenchymal cells during regeneration. Furthermore, we identify resident immune cells in the fin with a distinct clonal origin from other blood cell types. We envision that similar approaches will have major applications in other experimental systems, in which the matching of embryonic clonal origin to adult cell type will ultimately allow reconstruction of how the adult body is built from a single cell.

Single-cell qPCR facilitates the optimization of hematopoietic differentiation in hPSCs/OP9 coculture system

Human pluripotent stem cells (hPSCs)/OP9 coculture system is a widely used hematopoietic differentiation approach. The limited understanding of this process leads to its low efficiency. Thus, we used single-cell qPCR to reveal the gene expression profiles of individual CD34⁺ cells from different stages of differentiation. According to the dynamic gene expression of hematopoietic transcription factors, we overexpressed specific hematopoietic transcription factors (Gata2, Lmo2, Etv2, ERG, and SCL) at an early stage of hematopoietic differentiation. After overexpression, we generated more CD34⁺ cells with normal expression level of CD43 and CD31, which are used to define various hematopoietic progenitors. Furthermore, these CD34⁺ cells possessed normal differentiation potency in colony-forming unit assays and normal gene expression profiles. In this study, we demonstrated that single-cell qPCR can provide guidance for optimization of hematopoietic differentiation and transient overexpression of selected hematopoietic transcription factors can enhance hematopoietic differentiation.

Modeling Myeloid Malignancies Using Zebrafish

Human myeloid malignancies represent a substantial disease burden to individuals, with significant morbidity and death. The genetic underpinnings of disease formation and progression remain incompletely understood. Large-scale human population studies have identified a high frequency of potential driver mutations in spliceosomal and epigenetic regulators that contribute to malignancies, such as myelodysplastic syndromes (MDS) and leukemias. The high conservation of cell types and genes between humans and model organisms permits the investigation of the underlying mechanisms of leukemic development and potential therapeutic testing in genetically pliable pre-clinical systems. Due to the many technical advantages, such as large-scale screening, lineage-tracing studies, tumor transplantation, and high-throughput drug screening approaches, zebrafish is emerging as a model system for myeloid malignancies. In this review, we discuss recent advances in MDS and leukemia using the zebrafish model.

Identification and characterization of T reg-like cells in zebrafish

Kasheta et al. report the identification of T reg–like cells in zebrafish, a means to track and live-image these cells, and foxp3a-deficient mutants that display lymphoproliferation, severe inflammation, and other hallmarks of T reg deficiency syndromes.

Regulatory T (T reg) cells are a specialized sublineage of T lymphocytes that suppress autoreactive T cells. Functional studies of T reg cells in vitro have defined multiple suppression mechanisms, and studies of T reg–deficient humans and mice have made clear the important role that these cells play in preventing autoimmunity. However, many questions remain about how T reg cells act in vivo. Specifically, it is not clear which suppression mechanisms are most important, where T reg cells act, and how they get there. To begin to address these issues, we sought to identify T reg cells in zebrafish, a model system that provides unparalleled advantages in live-cell imaging and high-throughput genetic analyses. Using a FOXP3 orthologue as a marker, we identified CD4-enriched, mature T lymphocytes with properties of T reg cells. Zebrafish mutant for foxp3a displayed excess T lymphocytes, splenomegaly, and a profound inflammatory phenotype that was suppressed by genetic ablation of lymphocytes. This study identifies T reg–like cells in zebrafish, providing both a model to study the normal functions of these cells in vivo and mutants to explore the consequences of their loss.

Dissecting hematopoietic and renal cell heterogeneity in adult zebrafish at single-cell resolution using RNA sequencing

The work by Tang et al. provides a comprehensive, single-cell, transcriptomic analysis of kidney and blood cells from the adult zebrafish, identifying novel cell types, including two classes of NK immune cells, classically defined and erythroid-primed hematopoietic stem and progenitor cells, mucin-secreting kidney cells, and kidney stem/progenitor cells.

Recent advances in single-cell, transcriptomic profiling have provided unprecedented access to investigate cell heterogeneity during tissue and organ development. In this study, we used massively parallel, single-cell RNA sequencing to define cell heterogeneity within the zebrafish kidney marrow, constructing a comprehensive molecular atlas of definitive hematopoiesis and functionally distinct renal cells found in adult zebrafish. Because our method analyzed blood and kidney cells in an unbiased manner, our approach was useful in characterizing immune-cell deficiencies within DNA–protein kinase catalytic subunit (prkdc), interleukin-2 receptor γ a (il2rga), and double-homozygous–mutant fish, identifying blood cell losses in T, B, and natural killer cells within specific genetic mutants. Our analysis also uncovered novel cell types, including two classes of natural killer immune cells, classically defined and erythroid-primed hematopoietic stem and progenitor cells, mucin-secreting kidney cells, and kidney stem/progenitor cells. In total, our work provides the first, comprehensive, single-cell, transcriptomic analysis of kidney and marrow cells in the adult zebrafish.

Zebrafish B Cell Development without a Pre-B Cell Stage, Revealed by CD79 Fluorescence Reporter Transgenes

CD79a and CD79b proteins associate with Ig receptors as integral signaling components of the B cell Ag receptor complex. To study B cell development in zebrafish, we isolated orthologs of these genes and performed in situ hybridization, finding that their expression colocalized with IgH-μ in the kidney, which is the site of B cell development. CD79 transgenic lines were made by linking the promoter and upstream regulatory segments of CD79a and CD79b to enhanced GFP to identify B cells, as demonstrated by PCR analysis of IgH-μ expression in sorted cells. We crossed these CD79-GFP lines to a recombination activating gene (Rag)2:mCherry transgenic line to identify B cell development stages in kidney marrow. Initiation of CD79:GFP expression in Rag2:mCherry⁺ cells and the timing of Ig H and L chain expression revealed simultaneous expression of both IgH-μ- and IgL-κ-chains, without progressing through the stage of IgH-μ-chain alone. Rag2:mCherry⁺ cells without CD79:GFP showed the highest Rag1 and Rag2 mRNAs compared with CD79a and CD79b:GFP⁺ B cells, which showed strongly reduced Rag mRNAs. Thus, B cell development in zebrafish does not go through a Rag^hi CD79⁺IgH-μ⁺ pre-B cell stage, different from mammals. After the generation of CD79:GFP⁺ B cells, decreased CD79 expression occurred upon differentiation to Ig secretion, as detected by alteration from membrane to secreted IgH-μ exon usage, similar to in mammals. This confirmed a conserved role for CD79 in B cell development and differentiation, without the requirement of a pre-B cell stage in zebrafish.

Identification of innate lymphoid cells in single-cell RNA-Seq data

Innate lymphoid cells (ILCs) consist of natural killer (NK) cells and non-cytotoxic ILCs that are broadly classified into ILC1, ILC2, and ILC3 subtypes. These cells recently emerged as important early effectors of innate immunity for their roles in tissue homeostasis and inflammation. Over the last few years, ILCs have been extensively studied in mouse and human at the functional and molecular level, including gene expression profiling. However, sorting ILCs with flow cytometry for gene expression analysis is a delicate and time-consuming process. Here we propose and validate a novel framework for studying ILCs at the transcriptomic level using single-cell RNA-Seq data. Our approach combines unsupervised clustering and a new cell type classifier trained on mouse ILC gene expression data. We show that this approach can accurately identify different ILCs, especially ILC2 cells, in human lymphocyte single-cell RNA-Seq data. Our new model relies only on genes conserved across vertebrates, thereby making it in principle applicable in any vertebrate species. Considering the rapid increase in throughput of single-cell RNA-Seq technology, our work provides a computational framework for studying ILC2 cells in single-cell transcriptomic data and may help exploring their conservation in distant vertebrate species.

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

The immune system of vertebrate species consists of many different cell types that have distinct functional roles and are subject to different evolutionary pressures. Here, we first analyzed conservation of genes specific for all major immune cell types in human and mouse. Our results revealed higher gene turnover and faster evolution of trans-membrane proteins in NK cells compared with other immune cell types, and especially T cells, but similar conservation of nuclear and cytoplasmic protein coding genes. To validate these findings in a distant vertebrate species, we used single-cell RNA sequencing of lck:GFP cells in zebrafish and obtained the first transcriptome of specific immune cell types in a nonmammalian species. Unsupervised clustering and single-cell TCR locus reconstruction identified three cell populations, T cells, a novel type of NK-like cells, and a smaller population of myeloid-like cells. Differential expression analysis uncovered new immune-cell-specific genes, including novel immunoglobulin-like receptors, and neofunctionalization of recently duplicated paralogs. Evolutionary analyses confirmed the higher gene turnover of trans-membrane proteins in NK cells compared with T cells in fish species, suggesting that this is a general property of immune cell types across all vertebrates.