T lymphocytes are classified into CD4+ and CD8+ T cells based on the glycoproteins expressed on their cell surface. Teleosts possess two CD4 paralogs, CD4-1 and CD4-2; however, their functional roles in T cell differentiation remain unclear. We generated fish mutant lines lacking the two cd4 paralog genes and genes encoding the two components of the CD8 co-receptor. The resulting phenotypes were analyzed using polyclonal antisera specific to medaka (Oryzias latipes) CD4-1 and CD8α, raised using a DNA immunization strategy. Our results indicate that cd4-2 deficiency does not affect the development of CD4-1+ T cells, providing genetic evidence for their distinct functions. Interestingly, in cd4-1 deficient mutants, the smaller compartment of CD4-1+ T cells was accompanied by a relative increase of CD8α+ T cells. Reciprocally, the lack or reduction of CD8α+ T cells caused by the deletion of cd8 genes was accompanied by a greater preponderance of CD4-1+ T cells. For the first time, our study provides direct genetic evidence for the interconnected differentiation programs of CD4+ and CD8+ T cells in teleost.

Edwardsiella piscicida is a potent intracellular pathogen that infects a wide variety of freshwater and marine fish, causing systemic hemorrhagic sepsis that results in lethality. Although vaccines against E. piscicida have been developed for many fish species, a direct comparison between innate and adaptive immunity in protecting fish from E. piscicida infection has not been performed. We established an il2rb mutant lacking natural killer (NK) cells, innate type lymphocytes, showing that NK cells are required for protection against E. piscicida infection in adult medaka. In addition, we elucidated the role of the TLR signaling pathway transmitted by Myd88 in resistance to E. piscicida infection in larvae and adults. These results indicate the requirement of innate immunity in resistance to this intracellular pathogen, providing insights into vaccination strategies for fish.

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.

The stability of cellular phenotypes in developing organisms depends on error-free transmission of epigenetic and genetic information during mitosis. Methylation of cytosine residues in genomic DNA is a key epigenetic mark that modulates gene expression and prevents genome instability. Here, we report on a genetic test of the relationship between DNA replication and methylation in the context of the developing vertebrate organism instead of cell lines. Our analysis is based on the identification of hypomorphic alleles of dnmt1 , encoding the DNA maintenance methylase Dnmt1, and pole1 , encoding the catalytic subunit of leading-strand DNA polymerase epsilon holoenzyme (Pole). Homozygous dnmt1 mutants exhibit genome-wide DNA hypomethylation, whereas the pole1 mutation is associated with increased DNA methylation levels. In dnmt1 / pole1 double-mutant zebrafish larvae, DNA methylation levels are restored to near normal values, associated with partial rescue of mutant-associated transcriptional changes and phenotypes. Hence, a balancing antagonism between DNA replication and maintenance methylation buffers against replicative errors contributing to the robustness of vertebrate development.

The dorsal and anal fins can vary widely in position and length along the anterior-posterior axis in teleost fishes. However, the molecular mechanisms underlying the diversification of these fins remain unknown. Here, we used genetic approaches in zebrafish and medaka, in which the relative positions of the dorsal and anal fins are opposite, to demonstrate the crucial role of hox genes in the patterning of the teleost posterior body, including the dorsal and anal fins. By the CRISPR-Cas9-induced frameshift mutations and positional cloning of spontaneous dorsalfinless medaka, we show that various hox mutants exhibit the absence of dorsal or anal fins, or a stepwise posterior extension of these fins, with vertebral abnormalities. Our results indicate that multiple hox genes, primarily from hoxc-related clusters, encompass the regions responsible for the dorsal and anal fin formation along the anterior-posterior axis. These results further suggest that shifts in the anterior boundaries of hox expression which vary among fish species, lead to diversification in the position and size of the dorsal and anal fins, similar to how modulations in Hox expression can alter the number of anatomically distinct vertebrae in tetrapods. Furthermore, we show that hox genes responsible for dorsal fin formation are different between zebrafish and medaka. Our results suggest that a novel mechanism has occurred during teleost evolution, in which the gene network responsible for fin formation might have switched to the regulation downstream of other hox genes, leading to the remarkable diversity in the dorsal fin position.

The vertebral column is a characteristic structure of vertebrates. Genetic studies in mice have shown that Hox-mediated patterning plays a key role in specifying discrete anatomical regions of the vertebral column. Expression pattern analyses in several vertebrate embryos have provided correlative evidence that the anterior boundaries of Hox expression coincide with distinct anatomical vertebrae. However, because functional analyses have been limited to mice, it remains unclear which Hox genes actually function in vertebral patterning in other vertebrates. In this study, various zebrafish Hox mutants were generated for loss-of-function phenotypic analysis to functionally decipher the Hox code responsible for the zebrafish anterior vertebrae between the occipital and thoracic vertebrae. We found that Hox genes in HoxB- and HoxC-related clusters participate in regulating the morphology of the zebrafish anterior vertebrae. In addition, medaka hoxc6a was found to be responsible for anterior vertebral identity, as in zebrafish. Based on phenotypic similarities with Hoxc6 knockout mice, our results suggest that the Hox patterning system, including at least Hoxc6, may have been functionally established in the vertebral patterning of the common ancestor of ray-finned and lobe-finned fishes.

The interactions between the host immune system and intestinal microorganisms have been studied in many animals, including fish. However, a detailed analysis has not been performed in medaka, an established fish model for biological studies. Here, we investigated the effect of immunodeficiency on the microbiota composition and the effect of gut bacteria on intestinal epithelial development and immune responses in medaka. Chronological analysis of the intestinal microbiota of interleukin 2 receptor subunit gamma (il2rg) mutant medaka showed a gradual decrease in the evenness of operational taxonomic units, mainly caused by the increased abundance of the Aeromonadaceae family. Exposure of wild-type medaka to high doses of an intestine-derived opportunistic bacterium of the Aeromonadaceae family induced an inflammatory response, suggesting a harmful effect on adult il2rg mutants. In addition, we established germ-free conditions in larval medaka and observed large absorptive vacuoles in intestinal epithelial cells, indicating a block in epithelial maturation. Transcriptome analysis revealed a decrease in the expression of genes involved in the defense response, including the antimicrobial peptide gene hepcidin , whose expression is induced by lipopolysaccharide stimulation in normal larvae. These results show that reciprocal interactions between the microbiome and the intestinal tract are required for the maturation of the medaka immune system.

Unconventional T cells, such as innate natural killer T cells (iNKT) cells, are an important part of vertebrate immune defences. iNKT recognise glycolipids through a T cell receptor (TCR) that is composed of a semi-invariant TCR α chain, paired with a restricted set of TCR β chains. Here, we show that splicing of the cognate Trav11-Traj18-Trac pre-mRNA encoding the characteristic Vα14Jα18 variable region of this semi-invariant TCR depends on the presence of Tnpo3. The Tnpo3 gene encodes a nuclear transporter of the β-karyopherin family whose cargo includes various splice regulators. The block of iNKT cell development in the absence of Tnpo3 can be overcome by transgenic provision of a rearranged Trav11-Traj18-Trac cDNA, indicating that Tnpo3 deficiency does not interfere with the development of iNKT cells per se. Our study thus identifies a role for Tnpo3 in regulating the splicing of the pre-mRNA encoding the cognate TCRα chain of iNKT cells.

The development of the parathyroid and the thymus from the third pharyngeal pouch depends on the activities of the Gcm2 and Foxn1 transcription factors, respectively, whose expression domains sharply demarcate two regions in the developing third pharyngeal pouch. Here, we have generated novel mouse models to examine whether ectopic co-expression of Gcm2 in the thymic epithelium and of Foxn1 in the parathyroid perturbs the establishment of organ fates in vivo. Expression of Gcm2 in the thymic rudiment does not activate a parathyroid-specific expression programme, even in the absence of Foxn1 activity. Co-expression of Foxn1 in the parathyroid fails to impose thymopoietic capacity. We conclude that the actions of Foxn1 and Gcm2 transcription factors are cell context-dependent and that they each require permissive transcription factor landscapes in order to successfully interfere with organ-specific cell fate.

Hermansky-Pudlak syndrome is an autosomal recessive disease characterized by albinism, visual impairment, and blood platelet dysfunction. One of the genes responsible for Hermansky-Pudlak syndrome, hps1, regulates organelle biogenesis and thus plays important roles in melanin production, blood clotting, and the other organelle-related functions in humans and mice. However, the function of hps1 in other species remains poorly understood. In this study, we discovered albino medaka fish during the maintenance of a wild-derived population and identified hps1 as the responsible gene using positional cloning. In addition to the specific absence of melanophore pigmentation, the hps1 mutant showed reduced blood coagulation, suggesting that hps1 is involved in clotting caused by both mammalian platelets and fish thrombocytes. Together, the findings of our study demonstrate that hps1 has an evolutionarily conserved role in melanin production and blood coagulation. In addition, our study presents a useful vertebrate model for understanding the molecular mechanisms of Hermansky-Pudlak syndrome.

T cell differentiation in the thymus generates CD4+ helper and cytotoxic CD8+ cells as the two principal T cell lineages. Curiously, at the end of this complex selection process, CD4+ cells invariably outnumber CD8+ cells. Here, we examine the dynamics of repertoire formation and the emergence of the skewed CD4/CD8 ratio using high-resolution endogenous CRISPR/Cas9 barcoding that indelibly marks immature T cells at the DN2/DN3 pre-TCR stage. In wild-type mice, greater clone size of CD4+ cells and an intrinsically greater probability of Tcr β clonotypes for pMHCII interactions are major contributors to the skewed CD4/CD8 ratio. Clonal perturbations of thymocyte differentiation following the precocious expression of a rearranged iNKT invariant TCR α chain are due to loss of Tcr β clonotypes from the CD4 lineage-committed pre-selection repertoire. The present barcoding scheme offers a novel means to examine the clonal dynamics of lymphocyte differentiation orthogonal to that using TCR clonotypes.

Transcription factor EBF1 (early B cell factor 1) acts as a key regulator of B cell specification. The transcriptional network in which EBF1 operates has been extensively studied; however, the regulation of EBF1 function remains poorly defined. By mass spectrometric analysis of proteins associated with endogenous EBF1 in pro-B cells, we identified the nuclear import receptor Transportin-3 (Tnpo3) and found that it interacts with the immunoglobulin-like fold domain of EBF1. We delineated glutamic acid 271 of EBF1 as a critical residue for the association with Tnpo3. EBF1E271A showed normal nuclear localization; however, it had an impaired B cell programming ability in conditions of Notch signaling, as determined by retroviral transduction of Ebf1 -/- progenitors. By RNA-seq analysis of EBF1E271A-expressing progenitors, we found an up-regulation of T lineage determinants and down-regulation of early B genes, although similar chromatin binding of EBF1E271A and EBF1wt was detected in pro-B cells expressing activated Notch1. B lineage-specific inactivation of Tnpo3 in mice resulted in a block of early B cell differentiation, accompanied by a down-regulation of B lineage genes and up-regulation of T and NK lineage genes. Taken together, our observations suggest that Tnpo3 ensures B cell programming by EBF1 in nonpermissive conditions.

The rules underlying the structure of antigen receptor repertoires are not yet fully defined, despite their enormous importance for the understanding of adaptive immunity. With current technology, the large antigen receptor repertoires of mice and humans cannot be comprehensively studied. To circumvent the problems associated with incomplete sampling, we have studied the immunogenetic features of one of the smallest known vertebrates, the cyprinid fish Paedocypris sp. "Singkep" ("minifish"). Despite its small size, minifish has the key genetic facilities characterizing the principal vertebrate lymphocyte lineages. As described for mammals, the frequency distributions of immunoglobulin and T cell receptor clonotypes exhibit the features of fractal systems, demonstrating that self-similarity is a fundamental property of antigen receptor repertoires of vertebrates, irrespective of body size. Hence, minifish achieve immunocompetence via a few thousand lymphocytes organized in robust scale-free networks, thereby ensuring immune reactivity even when cells are lost or clone sizes fluctuate during immune responses.

To capture the global gene network regulating the differentiation of immature T cells in an unbiased manner, large-scale forward genetic screens in zebrafish were conducted and combined with genetic interaction analysis. After ENU mutagenesis, genetic lesions associated with failure of T cell development were identified by meiotic recombination mapping, positional cloning, and whole genome sequencing. Recessive genetic variants in 33 genes were identified and confirmed as causative by additional experiments. The mutations affected T cell development but did not perturb the development of an unrelated cell type, growth hormone-expressing somatotrophs, providing an important measure of cell-type specificity of the genetic variants. The structure of the genetic network encompassing the identified components was established by a subsequent genetic interaction analysis, which identified many instances of positive (alleviating) and negative (synthetic) genetic interactions. Several examples of synthetic lethality were subsequently phenocopied using combinations of small molecule inhibitors. These drugs not only interfered with normal T cell development, but also elicited remission in a model of T cell acute lymphoblastic leukaemia. Our findings illustrate how genetic interaction data obtained in the context of entire organisms can be exploited for targeted interference with specific cell types and their malignant derivatives.

DNA methylation is a universal epigenetic mechanism involved in regulation of gene expression and genome stability. The DNA maintenance methylase DNMT1 ensures that DNA methylation patterns are faithfully transmitted to daughter cells during cell division. Because loss of DNMT1 is lethal, a pan-organismic analysis of DNMT1 function is lacking. We identified new recessive dnmt1 alleles in medaka and zebrafish and, guided by the structures of mutant proteins, generated a recessive variant of mouse Dnmt1. Each of the three missense mutations studied here distorts the catalytic pocket and reduces enzymatic activity. Because all three DNMT1 mutant animals are viable, it was possible to examine their phenotypes throughout life. The consequences of genome-wide hypomethylation of DNA of somatic tissues in the Dnmt1 mutants are surprisingly mild but consistently affect the development of the lymphoid lineage. Our findings indicate that developing lymphocytes in vertebrates are sensitive to perturbations of DNA maintenance methylation.

Evidence for transgenerational inheritance of epigenetic information in vertebrates is scarce. Aberrant patterns of DNA methylation in gametes may set the stage for transmission into future generations. Here, we describe a viable hypomorphic allele of dnmt1 in zebrafish that causes widespread demethylation of CpG dinucleotides in sperm and somatic tissues. We find that homozygous mutants are essentially normal, with the exception of drastically impaired lymphopoiesis, affecting both larval and adult phases of T cell development. The phenotype of impaired larval (but not adult) T cell development is transmitted to subsequent generations by genotypically wildtype fish. We further find that about 200 differentially methylated regions in sperm DNA of transmitting and non-transmitting males, including hypermethylated sites associated with runx3 and rptor genes, whose reduced activities are associated with impaired larval T cell development. Our results indicate a particular sensitivity of larval T cell development to transgenerationally inherited epimutations.

The antibodies of jawless vertebrates consist of leucine-rich repeat arrays encoded by somatically assembled VLRB genes. It is unknown how the incomplete germline VLRB loci are converted into functional antibody genes during B lymphocyte development in lampreys. In Lampetra planeri larvae lacking the cytidine deaminase CDA2 gene, VLRB assembly fails, whereas the T lineage-associated VLRA and VLRC antigen receptor gene assemblies occur normally. Thus, CDA2 acts in a B cell lineage-specific fashion to support the somatic diversification of VLRB antibody genes. CDA2 is closely related to activation-induced cytidine deaminase (AID), which is essential for the elaboration of immunoglobulin gene repertoires in jawed vertebrates. Our results thus identify a convergent mechanism of antigen receptor gene assembly and diversification that independently evolved in the two sister branches of vertebrates.

Immunodeficiencies are typically caused by loss-of-function mutations in lymphocyte-specific genes. Occasionally, mutations in ubiquitous general-purpose factors, including those affecting essential components of the DNA polymerase epsilon (POLE) holoenzyme, have cell-type-specific consequences. POLE3, one of the four components of the POLE holoenzyme, features a histone fold domain and a unique acidic C terminus, making it a particularly attractive candidate mediating cell type-specific activities of POLE. Mice lacking Pole3 survive up to late embryonic stages, indicating that this subunit is dispensable for DNA replication. The phenotypes of viable hypomorphic and neomorphic alleles are surprisingly tissue restricted and reveal a stage-specific function of the histone fold domain of Pole3 during T and B cell development. Gradual introduction of positively charged residues into the acidic C terminus leads to peripheral lymphopenia of increasing severity. Our findings serve as a paradigm to understand the molecular basis of cell-type-specific non-replicative functions of the ubiquitous POLE complex.