The Hox code responsible for the patterning of the anterior vertebrae in zebrafish.

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.

Teleost Hox code defines regional identities competent for the formation of dorsal and anal fins.

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.

Antagonistic interactions safeguard mitotic propagation of genetic and epigenetic information in zebrafish.

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.

Maturation of the medaka immune system depends on reciprocal interactions between the microbiota and the intestinal tract.

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.

Tnpo3 controls splicing of the pre-mRNA encoding the canonical TCR α chain of iNKT cells.

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.

Tnpo3 enables EBF1 function in conditions of antagonistic Notch signaling.

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.

Clonal dynamics underlying the skewed CD4/CD8 ratio of mouse thymocytes revealed by TCR-independent barcoding.

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.

Limits to in vivo fate changes of epithelia in thymus and parathyroid by ectopic expression of transcription factors Gcm2 and Foxn1.

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.

Evolutionarily conserved role of hps1 in melanin production and blood coagulation in medaka fish.

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.

Genetic landscape of T cells identifies synthetic lethality for T-ALL.

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.

Antigen receptor repertoires of one of the smallest known vertebrates.

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.

Transgenerational inheritance of impaired larval T cell development in zebrafish.

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.

Epigenetic Protection of Vertebrate Lymphoid Progenitor Cells by Dnmt1.

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.

Lymphocyte-Specific Function of the DNA Polymerase Epsilon Subunit Pole3 Revealed by Neomorphic Alleles.

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.

Cytidine deaminase 2 is required for VLRB antibody gene assembly in lampreys.

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.

Evolutionary transition from degenerate to nonredundant cytokine signaling networks supporting intrathymic T cell development.

In mammals, T cell development critically depends on the IL-7 cytokine signaling pathway. Here we describe the identification of the zebrafish ortholog of mammalian IL-7 based on chromosomal localization, deduced protein sequence, and expression patterns. To examine the biological role of il7 in teleosts, we generated an il7 allele lacking most of its coding exons using CRISPR/Cas9-based mutagenesis. il7-deficient animals are viable and exhibit no obvious signs of immune disorder. With respect to intrathymic T cell development, il7 deficiency is associated with only a mild reduction of thymocyte numbers, contrasting with a more pronounced impairment of T cell development in il7r-deficient fish. Genetic interaction studies between il7 and il7r mutants, and il7 and crlf2(tslpr) mutants suggest the contribution of additional, as-yet unidentified cytokines to intrathymic T cell development. Such activities were also ascertained for other cytokines, such as il2 and il15, collectively indicating that in contrast to the situation in mammals, T cell development in the thymus of teleosts is driven by a degenerate multicomponent network of γc cytokines; this explains why deficiencies of single components have little detrimental effect. In contrast, the dependence on a single cytokine in the mammalian thymus has catastrophic consequences in cases of congenital deficiencies in genes affecting the IL-7 signaling pathway. We speculate that the transition from a degenerate to a nonredundant cytokine network supporting intrathymic T cell development emerged as a consequence of repurposing evolutionarily ancient constitutive cytokine pathways for regulatory functions in the mammalian peripheral immune system.

A small fish model for quantitative analysis of radiation effects using visualized thymus responses in GFP transgenic medaka.

Purpose: The aim of this study was to establish a new method of real-time, in vivo detection of radiation damage and recovery. Methods: The thymus was observed under fluorescent light in a green fluorescent protein transgenic medaka. After irradiation, medaka thymus images were analyzed to quantify the effects of radiation by measuring changes in thymus size. A single acute irradiation of X-rays (0-30 Gy) or heavy Fe ions (0-10 Gy) was delivered to the medaka. Images were captured 0, 1, 2, 3, 5, 7, 11, and 21 d after irradiation. Dose-response assessment was conducted to provide a direct measurement of the effects of the radiation. Conclusion: A biomonitoring system to detect the effects of radiation in real time was established. Using this system, the threshold doses for the induction of thymic atrophy by acute X-rays and Fe ions were 2-5 Gy and 0.5-1 Gy, respectively. The Relative Biological Effectiveness (RBE) of Fe-ion to X-rays was estimated to be around 3. This system may be used to evaluate the risk from concurrent exposure to hazards, such as chemicals and radiation, and for aging research.

Expansions, diversification, and interindividual copy number variations of AID/APOBEC family cytidine deaminase genes in lampreys.

Cytidine deaminases of the AID/APOBEC family catalyze C-to-U nucleotide transitions in mRNA or DNA. Members of the APOBEC3 branch are involved in antiviral defense, whereas AID contributes to diversification of antibody repertoires in jawed vertebrates via somatic hypermutation, gene conversion, and class switch recombination. In the extant jawless vertebrate, the lamprey, two members of the AID/APOBEC family are implicated in the generation of somatic diversity of the variable lymphocyte receptors (VLRs). Expression studies linked CDA1 and CDA2 genes to the assembly of VLRA/C genes in T-like cells and the VLRB genes in B-like cells, respectively. Here, we identify and characterize several CDA1-like genes in the larvae of different lamprey species and demonstrate that these encode active cytidine deaminases. Structural comparisons of the CDA1 variants highlighted substantial differences in surface charge; this observation is supported by our finding that the enzymes require different conditions and substrates for optimal activity in vitro. Strikingly, we also found that the number of CDA-like genes present in individuals of the same species is variable. Nevertheless, irrespective of the number of different CDA1-like genes present, all lamprey larvae have at least one functional CDA1-related gene encoding an enzyme with predicted structural and chemical features generally comparable to jawed vertebrate AID. Our findings suggest that, similar to APOBEC3 branch expansion in jawed vertebrates, the AID/APOBEC family has undergone substantial diversification in lamprey, possibly indicative of multiple distinct biological roles.