Search

Background: The spinal cord is a crucial part of the vertebrate CNS, controlling movements and receiving and processing sensory information from the trunk and limbs. However, there is much we do not know about how this essential organ develops. Here, we describe expression of 21 transcription factors and one transcriptional regulator in zebrafish spinal cord., Results: We analyzed the expression of aurkb, foxb1a, foxb1b, her8a, homeza, ivns1abpb, mybl2b, myt1a, nr2f1b, onecut1, sall1a, sall3a, sall3b, sall4, sox2, sox19b, sp8b, tsc22d1, wdhd1, zfhx3b, znf804a, and znf1032 in wild-type and MIB E3 ubiquitin protein ligase 1 zebrafish embryos. While all of these genes are broadly expressed in spinal cord, they have distinct expression patterns from one another. Some are predominantly expressed in progenitor domains, and others in subsets of post-mitotic cells. Given the conservation of spinal cord development, and the transcription factors and transcriptional regulators that orchestrate it, we expect that these genes will have similar spinal cord expression patterns in other vertebrates, including mammals and humans., Conclusions: Our data identify 22 different transcriptional regulators that are strong candidates for playing different roles in spinal cord development. For several of these genes, this is the first published description of their spinal cord expression., (© 2024 The Author(s). Developmental Dynamics published by Wiley Periodicals LLC on behalf of American Association for Anatomy.)

Introduction: Congenital hypopituitarism (CH) is characterized by the deficiency of pituitary hormones. Among CH patients, 85% lack a molecular diagnosis. Whole Exome Sequencing (WES) identified a homozygous variant (c.865G>A, p.Val289Ile) in the CDH2 gene, responsible for N-Cadherin production, crucial for cell-cell adhesion. Predicted to be likely pathogenic, the variant was found in a patient deficient in GH , TSH , ACTH , and LH/FSH . Its impact on cell adhesion was confirmed in L1 fibroblast cell lines., Objective: Create a cdh2 knockdown in zebrafish for investigating its role in pituitary development through growth hormone and transcription factors expression., Methods: Utilized pET28B-RfxCas13d-His plasmid for Cas13 mRNA production via in vitro transcription, guiding Cas13 to cdh2 with three RNAs. Injected the complex into single-cell embryos for analysis up to 96 hpf. Assessed gene expression of cdh2 , prop1 , pit1 , and gh1 using RT-qPCR. Evaluated cdh2 protein expression through the western blot technique., Results: Knockdown animals displayed developmental delay. The cdh2 expression decreased by 75% within 24 hours, rebounded by 48 hours, and reached wild-type levels by 96 hpf. gh1 expression decreased at 48h but increased by 96 hpf, aligning with WT. No significant differences in prop1 and pit1 expression were observed., Conclusion: Our findings underscore cdh2 's role in pituitary development and hormonal regulation, offering insights for developmental biology research., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Ventura Fernandes, Junqueira, MacRae and Silveira de Carvalho.)

Isl1 has been described as an embryonic master control gene expressed in the pericloacal mesenchyme. Deletion of Isl1 from the genital mesenchyme in mice leads to an ectopic urethral opening and epispadias-like phenotype. Using genome wide association methods, we identified ISL1 as the key susceptibility gene for classic bladder exstrophy (CBE), comprising epispadias and exstrophy of the urinary bladder. The most significant marker (rs6874700) identified in our recent GWAS meta-analysis achieved a p value of 1.48 × 10 - 24 within the ISL1 region. In silico analysis of rs6874700 and all other genome-wide significant markers in Linkage Disequilibrium (LD) with rs6874700 (D' = 1.0; R 2  > 0.90) revealed marker rs2303751 (p value 8.12 × 10 - 20 ) as the marker with the highest regulatory effect predicted. Here, we describe a novel 1.2 kb intragenic promoter residing between 6.2 and 7.4 kb downstream of the ISL1 transcription starting site, which is located in the reverse DNA strand and harbors a binding site for EZH2 at the exact region of marker rs2303751. We show, that EZH2 silencing in HEK cells reduces ISL1 expression. We show that ezh2 -/- knockout (KO) zebrafish larvae display tissues specificity of ISL1 regulation with reduced expression of Isl1 in the pronephric region of zebrafish larvae. In addition, a shorter and malformed nephric duct is observed in ezh2 -/- ko zebrafish Tg(wt1ß:eGFP) reporter lines. Our study shows, that Ezh2 is a key regulator of Isl1 during urinary tract formation and suggests tissue specific ISL1 dysregulation as an underlying mechanism for CBE formation., (© 2024. The Author(s).)

Prenatal lethality associated with mouse knockout of Mettl16, a recently identified RNA N6-methyladenosine (m 6 A) methyltransferase, has hampered characterization of the essential role of METTL16-mediated RNA m 6 A modification in early embryonic development. Here, using cross-species single-cell RNA sequencing analysis, we found that during early embryonic development, METTL16 is more highly expressed in vertebrate hematopoietic stem and progenitor cells (HSPCs) than other methyltransferases. In Mettl16-deficient zebrafish, proliferation capacity of embryonic HSPCs is compromised due to G1/S cell cycle arrest, an effect whose rescue requires Mettl16 with intact methyltransferase activity. We further identify the cell-cycle transcription factor mybl2b as a directly regulated by Mettl16-mediated m 6 A modification. Mettl16 deficiency resulted in the destabilization of mybl2b mRNA, likely due to lost binding by the m 6 A reader Igf2bp1 in vivo. Moreover, we found that the METTL16-m 6 A-MYBL2-IGF2BP1 axis controlling G1/S progression is conserved in humans. Collectively, our findings elucidate the critical function of METTL16-mediated m 6 A modification in HSPC cell cycle progression during early embryonic development., (© 2024. The Author(s).)

Brain pericytes are one of the critical cell types that regulate endothelial barrier function and activity, thus ensuring adequate blood flow to the brain. The genetic pathways guiding undifferentiated cells into mature pericytes are not well understood. We show here that pericyte precursor populations from both neural crest and head mesoderm of zebrafish express the transcription factor nkx3.1 develop into brain pericytes. We identify the gene signature of these precursors and show that an nkx3.1-, foxf2a-, and cxcl12b-expressing pericyte precursor population is present around the basilar artery prior to artery formation and pericyte recruitment. The precursors later spread throughout the brain and differentiate to express canonical pericyte markers. Cxcl12b-Cxcr4 signaling is required for pericyte attachment and differentiation. Further, both nkx3.1 and cxcl12b are necessary and sufficient in regulating pericyte number as loss inhibits and gain increases pericyte number. Through genetic experiments, we have defined a precursor population for brain pericytes and identified genes critical for their differentiation., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Ahuja et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Transcription factors (TFs) are crucial epigenetic regulators, which enable cells to dynamically adjust gene expression in response to environmental signals. Computational procedures like digital genomic footprinting on chromatin accessibility assays such as ATACseq can be used to identify bound TFs in a genome-wide scale. This method utilizes short regions of low accessibility signals due to steric hindrance of DNA bound proteins, called footprints (FPs), which are combined with motif databases for TF identification. However, while over 1600 TFs have been described in the human genome, only ~ 700 of these have a known binding motif. Thus, a substantial number of FPs without overlap to a known DNA motif are normally discarded from FP analysis. In addition, the FP method is restricted to organisms with a substantial number of known TF motifs. Here we present DENIS (DE Novo motIf diScovery), a framework to generate and systematically investigate the potential of de novo TF motif discovery from FPs. DENIS includes functionality (1) to isolate FPs without binding motifs, (2) to perform de novo motif generation and (3) to characterize novel motifs. Here, we show that the framework rediscovers artificially removed TF motifs, quantifies de novo motif usage during an early embryonic development example dataset, and is able to analyze and uncover TF activity in organisms lacking canonical motifs. The latter task is exemplified by an investigation of a scATAC-seq dataset in zebrafish which covers different cell types during hematopoiesis., (© 2024. The Author(s).)

Background: Elongation of the spinal cord is dependent on neural development from neuromesodermal progenitors in the tail bud. We previously showed the involvement of the Oct4-type gene, pou5f3, in this process in zebrafish mainly by dominant-interference gene induction, but, to compensate for the limitation of this transgene approach, mutant analysis was indispensable. pou5f3 involvement in the signaling pathways was another unsolved question., Results: We examined the phenotypes of pou5f3 mutants and the effects of Pou5f3 activation by the tamoxifen-ERT2 system in the posterior neural tube, together confirming the involvement of pou5f3. The reporter assays using P19 cells implicated tail bud-related transcription factors in pou5f3 expression. Regulation of tail bud development by retinoic acid (RA) signaling was confirmed by treatment of embryos with RA and the synthesis inhibitor, and in vitro reporter assays further showed that RA signaling regulated pou5f3 expression. Importantly, the expression of the RA degradation enzyme gene, cyp26a1, was down-regulated in embryos with disrupted pou5f3 activity., Conclusions: The involvement of pou5f3 in spinal cord extension was supported by using mutants and the gain-of-function approach. Our findings further suggest that pou5f3 regulates the RA level, contributing to neurogenesis in the posterior neural tube., (© 2023 The Authors. Developmental Dynamics published by Wiley Periodicals LLC on behalf of American Association for Anatomy.)

The aim of this work is the full characterization of all the nocturnin ( noc ) paralogues expressed in a teleost, the goldfish. An in silico analysis of the evolutive origin of noc in Osteichthyes is performed, including the splicing variants and new paralogues appearing after teleostean 3R genomic duplication and the cyprinine 4Rc. After sequencing the full-length mRNA of goldfish, we obtained two isoforms for noc-a ( noc-aa and noc-ab ) with two splice variants (I and II), and only one for noc-b ( noc-bb ) with two transcripts (II and III). Using the splicing variant II, the prediction of the secondary and tertiary structures renders a well-conserved 3D distribution of four α-helices and nine β-sheets in the three noc isoforms. A synteny analysis based on the localization of noc genes in the patrilineal or matrilineal subgenomes and a phylogenetic tree of protein sequences were accomplished to stablish a classification and a long-lasting nomenclature of noc in goldfish, and valid to be extrapolated to allotetraploid Cyprininae. Finally, both goldfish and zebrafish showed a broad tissue expression of all the noc paralogues. Moreover, the enriched expression of specific paralogues in some tissues argues in favour of neo- or subfunctionalization., Competing Interests: The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Monocyte chemoattractant protein-induced protein 1 (MCPIP1), also called Regnase-1, is an RNase that has been described as a key negative modulator of inflammation. MCPIP1 also controls numerous tumor-related processes, such as proliferation, apoptosis and differentiation. In this study, we utilized a zebrafish model to investigate the role of Mcpip1 during embryogenic development. Our results demonstrated that during embryogenesis, the expression of the zc3h12a gene encoding Mcpip1 undergoes dynamic changes. Its transcript levels gradually increase from the 2-cell stage to the spherical stage and then decrease rapidly. We further found that ectopic overexpression of wild-type Mcpip1 but not the catalytically inactive mutant form resulted in an embryonic lethal phenotype in zebrafish embryos (24 hpf). At the molecular level, transcriptomic profiling revealed extensive changes in the expression of genes encoding proteins important in the endoplasmic reticulum stress response and in protein folding as well as involved in the formation of primary germ layer, mesendoderm and endoderm development, heart morphogenesis and cell migration. Altogether, our results demonstrate that the expression of zc3h12a must be tightly controlled during the first cell divisions of zebrafish embryos and that a rapid decrease in its mRNA expression is an important factor promoting proper embryo development., (© 2023. Springer Nature Limited.)

Neural progenitor cells are self-renewable, proliferative, and multipotent cell populations that generate diverse types of neurons and glia to build the nervous system. Transcription factors play critical roles in regulating various cellular processes; however, the transcription factors that regulate the development of neural progenitors are yet to be identified. In the present study, we demonstrated that zebrafish etv5a is expressed in the neural progenitor cells of the neuroectoderm. Downregulation of endogenous Etv5a function by etv5a morpholino or an etv5a dominant-negative variant increased the proliferation of sox2 -positive neural progenitor cells, accompanied by inhibition of neurogenesis and gliogenesis. These phenotypes in Etv5a-depleted embryos could be rescued by a co-injection with etv5a cRNA. Etv5a overexpression reduced sox2 expression. Direct binding of Etv5a to the regulatory elements of sox2 was affirmed by chromatin immunoprecipitation. These data revealed that Etv5a directly suppressed sox2 expression to reduce the proliferation of neural progenitor cells. In addition, the expression of foxm1 , a putative target gene of Etv5a and a direct upstream transcription factor of sox2 , was upregulated in Etv5a-deficient embryos. Moreover, the suppression of Foxm1 function by the foxm1 dominant-negative construct nullified the phenotype of upregulated sox2 expression caused by Etv5a deficiency. Overall, our results indicated that Etv5a regulates the expression of sox2 via direct binding to the sox2 promoter and indirect regulation by inhibiting foxm1 expression. Hence, we revealed the role of Etv5a in the transcriptional hierarchy that regulates the proliferation of neural progenitor cells.

Background: In vertebrates, the development of the inner ear is a delicate process, whereas its relating molecular pathways are still poorly understood. LMO4 , an LIM domain-only transcriptional regulator, is drawing an increasing amount of interest for its multiple roles regarding human embryonic development and the modulation of ototoxic side effects of cisplatin including cochlear apoptosis and hearing loss. The aim of the present study is to further explore the role of lmo4a in zebrafish inner ear development and thus explore its functional role., Methods: The Spatial Transcript Omics DataBase was referred to in order to evaluate the expression of lmo4a during the first 24 h of zebrafish development. In situ hybridization was applied to validate and extend the expression profile of lmo4a to 3 days post-fertilization. The morpholino (MO) knockdown and CRISPR/Cas9 knockout (KO) of lmo4a was applied. Morphological analyses of otic vesical, hair cells, statoacoustic ganglion and semicircular canals were conducted. The swimming pattern of lmo4a KO and MO zebrafish was tracked. In situ hybridization was further applied to verify the expression of genes of the related pathways. Rescue of the phenotype was attempted by blockage of the bmp pathway via heat shock and injection of Dorsomorphin., Results: lmo4a is constitutively expressed in the otic placode and otic vesicle during the early stages of zebrafish development. Knockdown and knockout of lmo4a both induced smaller otocysts, less hair cells, immature statoacoustic ganglion and malformed semicircular canals. Abnormal swimming patterns could be observed in both lmo4a MO and KO zebrafish. eya1 in preplacodal ectoderm patterning was downregulated. bmp2 and bmp4 expressions were found to be upregulated and extended in lmo4a morphants, and blockage of the Bmp pathway partially rescued the vestibular defects., Conclusions: We concluded that lmo4a holds a regulative effect on the Bmp pathway and is required for the normal development of zebrafish inner ear. Our study pointed out the conservatism of LMO4 in inner ear development between mammals and zebrafish as well as shed more light on the molecular mechanisms behind it. Further research is needed to distinguish the relationships between lmo4 and the Bmp pathway, which may lead to diagnostic and therapeutic approaches towards human inner ear malformation.

The hematopoietic niche is a supportive microenvironment composed of distinct cell types, including specialized vascular endothelial cells that directly interact with hematopoietic stem and progenitor cells (HSPCs). The molecular factors that specify niche endothelial cells and orchestrate HSPC homeostasis remain largely unknown. Using multi-dimensional gene expression and chromatin accessibility analyses in zebrafish, we define a conserved gene expression signature and cis-regulatory landscape that are unique to sinusoidal endothelial cells in the HSPC niche. Using enhancer mutagenesis and transcription factor overexpression, we elucidate a transcriptional code that involves members of the Ets, Sox, and nuclear hormone receptor families and is sufficient to induce ectopic niche endothelial cells that associate with mesenchymal stromal cells and support the recruitment, maintenance, and division of HSPCs in vivo. These studies set forth an approach for generating synthetic HSPC niches, in vitro or in vivo, and for effective therapies to modulate the endogenous niche., Competing Interests: Declaration of interests L.I.Z. is a founder and stockholder of Fate Therapeutics, CAMP4 Therapeutics, and Scholar Rock. He is a consultant for Celularity., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Pancreatic ε-cells producing ghrelin are one type of endocrine cell found in islets, which have been shown to influence other intra-islet cells, especially in regulating the function of β cells. However, the role of such cells during β-cell regeneration is currently unknown. Here, using a zebrafish nitroreductase (NTR)-mediated β-cell ablation model, we reveal that ghrelin-positive ε-cells in the pancreas act as contributors to neogenic β-cells after extreme β-cell loss. Further studies show that the overexpression of ghrelin or the expansion of ε-cells potentiates β-cell regeneration. Lineage tracing confirms that a proportion of embryonic ε-cells can transdifferentiate to β-cells, and that the deletion of Pax4 enhances this transdifferentiation of ε-cells to β-cells. Mechanistically, Pax4 binds to the ghrelin regulatory region and represses its transcription. Thus, deletion of Pax4 derepresses ghrelin expression and causes producing more ghrelin-positive cells, enhancing the transdifferentiation of ε-cells to β-cells and consequently potentiating β-cell regeneration. Our findings reveal a previously unreported role for ε-cells during zebrafish β-cell regeneration, indicating that Pax4 regulates ghrelin transcription and mediates the conversion of embryonic ε-cells to β-cells after extreme β-cell loss., Competing Interests: Competing interests The authors declare no competing or financial interests., (© 2023. Published by The Company of Biologists Ltd.)

Median fins are thought to be ancestors of paired fins which in turn give rise to limbs in tetrapods. However, the developmental mechanisms of median fins remain largely unknown. Nonsense mutation of the T-box transcription factor eomesa in zebrafish results in a phenotype without dorsal fin. Compared to zebrafish, the common carp undergo an additional round of whole genome duplication, acquiring an extra copy of protein-coding genes. To verify the function of eomesa genes in common carp, we established a biallelic gene editing technology in this tetraploidy fish through simultaneous disruption of two homologous genes, eomesa1 and eomesa2. We targeted four sites located upstream or within the sequences encoding the T-box domain. Sanger sequencing data indicated the average knockout efficiency was around 40% at T1-T3 sites and 10% at T4 site in embryos at 24 hours post fertilization. The individual editing efficiency was high to about 80% at T1-T3 sites and low to 13.3% at T4 site in larvae at 7 days post fertilization. Among 145 mosaic F0 examined at four months old, three individuals (Mutant 1-3) showed varying degrees of maldevelopment in the dorsal fin and loss of anal fin. Genotyping showed the genomes of all three mutants were disrupted at T3 sites. The null mutation rates on the eomesa1 and eomesa2 loci were 0% and 60% in Mutant 1, 66.7% and 100% in Mutant 2, and 90% and 77.8% in Mutant 3, respectively. In conclusion, we demonstrated a role of eomesa in the formation and development of median fins in Oujiang color common carp and established an method that simultaneously disrupt two homologous genes with one gRNA, which would be useful in genome editing in other polyploidy fishes., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Song et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Mutations in the mitochondrial or nuclear genomes are associated with a diverse group of human disorders characterized by impaired mitochondrial respiration. Within this group, an increasing number of mutations have been identified in nuclear genes involved in mitochondrial RNA biology. The TEFM gene encodes the mitochondrial transcription elongation factor responsible for enhancing the processivity of mitochondrial RNA polymerase, POLRMT. We report for the first time that TEFM variants are associated with mitochondrial respiratory chain deficiency and a wide range of clinical presentations including mitochondrial myopathy with a treatable neuromuscular transmission defect. Mechanistically, we show muscle and primary fibroblasts from the affected individuals have reduced levels of promoter distal mitochondrial RNA transcripts. Finally, tefm knockdown in zebrafish embryos resulted in neuromuscular junction abnormalities and abnormal mitochondrial function, strengthening the genotype-phenotype correlation. Our study highlights that TEFM regulates mitochondrial transcription elongation and its defect results in variable, tissue-specific neurological and neuromuscular symptoms., (© 2023. The Author(s).)

During development, retinal progenitors navigate a complex landscape of fate decisions to generate the major cell classes necessary for proper vision. Transcriptional regulation is critical to generate diversity within these major cell classes. Here, we aim to provide the resources and techniques required to identify transcription factors necessary to generate and maintain diversity in photoreceptor subtypes, which are critical for vision. First, we generate a key resource: a high-quality and deep transcriptomic profile of each photoreceptor subtype in adult zebrafish. We make this resource openly accessible, easy to explore, and have integrated it with other currently available photoreceptor transcriptomic datasets. Second, using our transcriptomic profiles, we derive an in-depth map of expression of transcription factors in photoreceptors. Third, we use efficient CRISPR-Cas9 based mutagenesis to screen for null phenotypes in F0 larvae (F0 screening) as a fast, efficient, and versatile technique to assess the involvement of candidate transcription factors in the generation of photoreceptor subtypes. We first show that known phenotypes can be easily replicated using this method: loss of S cones in foxq2 mutants and loss of rods in nr2e3 mutants. We then identify novel functions for the transcription factor Tbx2, demonstrating that it plays distinct roles in controlling the generation of all photoreceptor subtypes within the retina. Our study provides a roadmap to discover additional factors involved in this process. Additionally, we explore four transcription factors of unknown function (Skor1a, Sall1a, Lrrfip1a, and Xbp1), and find no evidence for their involvement in the generation of photoreceptor subtypes. This dataset and screening method will be a valuable way to explore the genes involved in many other essential aspects of photoreceptor biology., Competing Interests: JA, VK, LP, HK, KK, WL No competing interests declared

Cell identity is governed by the complex regulation of gene expression, represented as gene-regulatory networks 1 . Here we use gene-regulatory networks inferred from single-cell multi-omics data to perform in silico transcription factor perturbations, simulating the consequent changes in cell identity using only unperturbed wild-type data. We apply this machine-learning-based approach, CellOracle, to well-established paradigms-mouse and human haematopoiesis, and zebrafish embryogenesis-and we correctly model reported changes in phenotype that occur as a result of transcription factor perturbation. Through systematic in silico transcription factor perturbation in the developing zebrafish, we simulate and experimentally validate a previously unreported phenotype that results from the loss of noto, an established notochord regulator. Furthermore, we identify an axial mesoderm regulator, lhx1a. Together, these results show that CellOracle can be used to analyse the regulation of cell identity by transcription factors, and can provide mechanistic insights into development and differentiation., (© 2023. The Author(s), under exclusive licence to Springer Nature Limited.)

The localization of transcriptional activity in specialized transcription bodies is a hallmark of gene expression in eukaryotic cells. 1 How proteins of the transcriptional machinery come together to form such bodies, however, is unclear. Here, we take advantage of two large, isolated, and long-lived transcription bodies that reproducibly form during early zebrafish embryogenesis to characterize the dynamics of transcription body formation. Once formed, these transcription bodies are enriched for initiating and elongating RNA polymerase II, as well as the transcription factors Nanog and Sox19b. Analyzing the events leading up to transcription, we find that Nanog and Sox19b cluster prior to transcription. The clustering of transcription factors is sequential; Nanog clusters first, and this is required for the clustering of Sox19b and the initiation of transcription. Mutant analysis revealed that both the DNA-binding domain as well as one of the two intrinsically disordered regions of Nanog are required to organize the two bodies of transcriptional activity. Taken together, our data suggest that the clustering of transcription factors dictates the formation of transcription bodies.- 3 How proteins of the transcriptional machinery come together to form such bodies, however, is unclear. Here, we take advantage of two large, isolated, and long-lived transcription bodies that reproducibly form during early zebrafish embryogenesis to characterize the dynamics of transcription body formation. Once formed, these transcription bodies are enriched for initiating and elongating RNA polymerase II, as well as the transcription factors Nanog and Sox19b. Analyzing the events leading up to transcription, we find that Nanog and Sox19b cluster prior to transcription. The clustering of transcription factors is sequential; Nanog clusters first, and this is required for the clustering of Sox19b and the initiation of transcription. Mutant analysis revealed that both the DNA-binding domain as well as one of the two intrinsically disordered regions of Nanog are required to organize the two bodies of transcriptional activity. Taken together, our data suggest that the clustering of transcription factors dictates the formation of transcription bodies., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 Elsevier Inc. All rights reserved.)

The zinc finger transcription factor Ikaros1 (Ikzf1) is required for lymphoid development in mammals. Four zinc fingers constitute its DNA binding domain and two zinc fingers are present in the C-terminal protein interaction module. We describe the phenotypes of zebrafish homozygous for two distinct mutant ikzf1 alleles. The IT325 variant lacks the C-terminal two zinc fingers, whereas the fr105 variant retains only the first zinc finger of the DNA binding domain. An intact ikzf1 gene is required for larval T cell development, whereas low levels of adult lymphoid development recover in the mutants. By contrast, the mutants exhibit a signature of increased myelopoiesis at larval and adult stages. Both mutations stimulate erythroid differentiation in larvae, indicating that the C-terminal zinc fingers negatively regulate the extent of red blood cell production. An unexpected differential effect of the two mutants on adult erythropoiesis suggests a direct requirement of an intact DNA binding domain for entry of progenitors into the red blood cell lineage. Collectively, our results reinforce the biological differences between larval and adult haematopoiesis, indicate a stage-specific function of ikzf1 in regulating the hierarchical bifurcations of differentiation, and assign distinct functions to the DNA binding domain and the C-terminal zinc fingers., (© 2022. The Author(s).)

Technologies for gene activation are valuable tools for the study of gene functions and have a wide range of potential applications in bioengineering and medicine. In contrast to existing methods based on recruiting transcriptional modulators via DNA-binding proteins, we developed a strategy termed Narta (nascent RNA-guided transcriptional activation) to achieve gene activation by recruiting artificial transcription factors (aTFs) to transcription sites through nascent RNAs of the target gene. Using Narta, we demonstrate robust activation of a broad range of exogenous and endogenous genes in various cell types, including zebrafish embryos, mouse and human cells. Importantly, the activation is reversible, tunable and specific. Moreover, Narta provides better activation potency of some expressed genes than CRISPRa and, when used in combination with CRISPRa, has an enhancing effect on gene activation. Quantitative imaging illustrated that nascent RNA-directed aTFs could induce the high-density assembly of coactivators at transcription sites, which may explain the larger transcriptional burst size induced by Narta. Overall, our work expands the gene activation toolbox for biomedical research., (© 2022. The Author(s).)

The pMN domain is a restricted domain in the ventral spinal cord, defined by the expression of the olig2 gene. Though it is known that the pMN progenitor cells can sequentially generate motor neurons and oligodendrocytes, the lineages of these progenitors are controversial and how their progeny are generated is not well understood. Using single-cell RNA sequencing, here, we identified a previously unknown heterogeneity among pMN progenitors with distinct fates and molecular signatures in zebrafish. Notably, we characterized two distinct motor neuron lineages using bioinformatic analysis. We then went on to investigate specific molecular programs that regulate neural progenitor fate transition. We validated experimentally that expression of the transcription factor myt1 (myelin transcription factor 1) and inner nuclear membrane integral proteins lbr (lamin B receptor) were critical for the development of motor neurons and neural progenitor maintenance, respectively. We anticipate that the transcriptome features and molecular programs identified in zebrafish pMN progenitors will not only provide an in-depth understanding of previous findings regarding the lineage analysis of oligodendrocyte progenitor cells and motor neurons but will also help in further understanding of the molecular programming involved in neural progenitor fate transition., Competing Interests: Conflict of interest The authors declare no competing interests., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Members of the FOS protein family regulate gene expression responses to a multitude of extracellular signals and are dysregulated in several pathological states. Whilst mouse genetic models have provided key insights into the tissue-specific functions of these proteins in vivo, little is known about their roles during early vertebrate embryonic development. This study examined the potential of using zebrafish as a model for such studies and, more broadly, for investigating the mechanisms regulating the functions of Fos proteins in vivo. Through phylogenetic and sequence analysis, we identified six zebrafish FOS orthologues, fosaa, fosab, fosb, fosl1a, fosl1b, and fosl2 , which show high conservation in key regulatory domains and post-translational modification sites compared to their equivalent human proteins. During embryogenesis, zebrafish fos genes exhibit both overlapping and distinct spatiotemporal patterns of expression in specific cell types and tissues. Most fos genes are also expressed in a variety of adult zebrafish tissues. As in humans, we also found that expression of zebrafish FOS orthologs is induced by oncogenic BRAF-ERK signalling in zebrafish melanomas. These findings suggest that zebrafish represent an alternate model to mice for investigating the regulation and functions of Fos proteins in vertebrate embryonic and adult tissues, and cancer.

Many regulators controlling arterial identity are well described; however, transcription factors that promote vein identity and vascular patterning have remained largely unknown. We previously identified the transcription factors Islet2 (Isl2) and Nr2f1b required for specification of the vein and tip cell identity mediated by notch pathway in zebrafish. However, the interaction between Isl2 and Nr2f1b is not known. In this study, we report that Nr2f2 plays minor roles on vein and intersegmental vessels (ISV) growth and dissect the genetic interactions among the three transcription factors Isl2, Nr2f1b, and Nr2f2 using a combinatorial knockdown strategy. The double knockdown of isl2/nr2f1b, isl2/nr2f2, and nr2f1b/nr2f2 showed the enhanced defects in vasculature including less completed ISV, reduced veins, and ISV cells. We further tested the genetic relationship among these three transcription factors. We found isl2 can regulate the expression of nr2f1b and nr2f2, suggesting a model where Isl2 functions upstream of Nr2f1b and Nr2f2. We hypothsized that Isl2 and Nr2f1b can function together through cis-regulatory binding motifs. In-vitro luciferase assay results, we showed that Isl2 and Nr2f1b can cooperatively enhance gene expression. Moreover, co-immunoprecipitation results indicated that Isl2 and Nr2f1b interact physically. Together, we showed that the interaction of the Nr2f1b and Nr2f2 transcription factors in combination with the Islet2 play coordinated roles in the vascular development of zebrafish., Competing Interests: Declaration of competing interest The authors declare that they have no conflict of interest., (Copyright © 2022 Elsevier Inc. All rights reserved.)

The limited regenerative capacity in mammals has serious implications for cardiac tissue damage. Meanwhile, zebrafish has a high regenerative capacity, but the regulation of the heart healing process has yet to be elucidated. The dynamic nature of cardiac regeneration requires consideration of the inherent temporal dimension of this process. Here, we conducted a systematic review to find genes that define the regenerative cell state of the zebrafish heart. We then performed an in silico temporal gene regulatory network analysis using transcriptomic data from the zebrafish heart regenerative process obtained from databases. In this analysis, the genes found in the systematic review were used to represent the final cell state of the transition process from a non-regenerative cell state to a regenerative state. We found 135 transcription factors driving the cellular state transition process during zebrafish cardiac regeneration, including Hand2, Nkx2.5, Tbx20, Fosl1, Fosb, Junb, Vdr, Wt1, and Tcf21 previously reported for playing a key role in tissue regeneration. Furthermore, we demonstrate that most regulators are activated in the first days post-injury, indicating that the transition from a non-regenerative to a regenerative state occurs promptly., (Copyright © 2022 Elsevier B.V. All rights reserved.)

Gene expression is tightly regulated during hematopoiesis. Recent studies have suggested that RNA polymerase II (Pol II) promoter proximal pausing, a temporary stalling downstream of the promoter region after initiation, plays a critical role in regulating the expression of various genes in metazoans. However, the function of proximal pausing in hematopoietic gene regulation remains largely unknown. The negative elongation factor (NELF) complex is a key factor important for this proximal pausing. Previous studies have suggested that NELF regulates granulocytic differentiation in vitro, but its in vivo function during hematopoiesis remains uncharacterized. Here, we generated the zebrafish mutant for one NELF complex subunit Nelfb using the CRISPR-Cas9 technology. We found that the loss of nelfb selectively induced excessive granulocytic development during primitive and definitive hematopoiesis. The loss of nelfb reduced hematopoietic progenitor cell formation and did not affect erythroid development. Moreover, the accelerated granulocytic differentiation and reduced progenitor cell development could be reversed by inhibiting Pol II elongation. Further experiments demonstrated that the other NELF complex subunits (Nelfa and Nelfe) played similar roles in controlling granulocytic development. Together, our studies suggested that NELF is critical in controlling the proper granulocytic development in vivo, and that promoter proximal pausing might help maintain the undifferentiated state of hematopoietic progenitor cells.

The T-box containing family member, TBX5, has been shown to play important functional roles in the pectoral appendages of a variety of vertebrate species. While a single TBX5 gene exists in all tetrapods studied to date, the zebrafish genome retains two paralogues, designated as tbx5a and tbx5b, resulting from a whole genome duplication in the teleost lineage. Zebrafish deficient in tbx5a lack pectoral fin buds, whereas zebrafish deficient in tbx5b exhibit misshapen pectoral fins, showing that both paralogues function in fin development. The mesenchymal cells of the limb/fin bud are derived from the Lateral Plate Mesoderm (LPM). Previous fate mapping work in zebrafish has shown that wildtype (wt) fin field cells are initially located adjacent to somites (s)1-4. The wt fin field cells migrate in opposing diagonal directions placing the limb bud between s2-3 and lateral to the main body. To better characterize tbx5 paralogue functions in zebrafish, time-lapse analyses of the migrations of fin bud precursors under conditions of tbx5a knock-down, tbx5b knock-down and double-knock-down were performed. Our data suggest that zebrafish tbx5a and tbx5b have functionally separated migration direction vectors, that when combined recapitulate the migration of the wt fin field. We and others have shown that loss of Tbx5a function abolishes an fgf24 signaling cue resulting in fin field cells failing to converge in an Antero-Posterior (AP) direction and migrating only in a mediolateral (ML) direction. We show here that loss of Tbx5b function affects initial ML directed movements so that fin field cells fail to migrate laterally but continue to converge along the AP axis. Furthermore, fin field cells in the double Tbx5a/Tbx5b knock-down zebrafish do not engage in directed migrations along either the ML or AP axis. Therefore, these two paralogues may be acting to instruct separate vectors of fin field migration in order to direct proper fin bud formation., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Palatogenesis is affected by many factors, including gene polymorphisms and exposure to toxic chemicals during sensitive developmental periods. Cleft palate is one of the most common congenital anomalies, and ongoing efforts to elucidate the molecular mechanisms underlying palatogenesis are providing useful insights to reduce the risk of this disorder. To identify novel potential regulators of palatogenesis, we analyzed public transcriptome datasets from a mouse model of cleft palate caused by selective deletion of transforming growth factor-β (TGFβ) receptor type 2 in cranial neural crest cells. We identified the homeobox transcription factor Mohawk (Mkx) as a gene downregulated in the maxilla of TGFβ knockout mice compared with wild-type mice. To examine the role of mkx in palatogenesis, we used CRISPR/Cas9 editing to generate zebrafish with impaired expression of mkxa and mkxb, the zebrafish homologs of Mkx. We found that mkx crispants expressed reduced levels of gli1, a critical transcription factor in the Sonic hedgehog (SHH) signaling pathway that plays an important role in the regulation of palatogenesis. Furthermore, we found that mkxa -/- zebrafish were more susceptible than mkxa +/+ zebrafish to the deleterious effects of cyclopamine, an inhibitor of SHH signaling, on upper jaw development. These results suggest that Mkx may be involved in palatogenesis regulated by TGFβ and SHH signaling, and that impairment in Mkx function may be related to the etiology of cleft palate., (© 2021 Japanese Teratology Society.)

Mitochondrial disorders are collectively common, genetically heterogeneous disorders in both pediatric and adult populations. They are caused by molecular defects in oxidative phosphorylation, failure of essential bioenergetic supply to mitochondria, and apoptosis. Here, we present three affected individuals from a consanguineous family of Pakistani origin with variable seizures and intellectual disability. Both females display primary ovarian insufficiency (POI), while the male shows abnormal sex hormone levels. We performed whole exome sequencing and identified a recessive missense variant c.694C > T, p.Arg232Cys in TFAM that segregates with disease. TFAM (mitochondrial transcription factor A) is a component of the mitochondrial replisome machinery that maintains mtDNA transcription and replication. In primary dermal fibroblasts, we show depletion of mtDNA and significantly altered mitochondrial function and morphology. Moreover, we observed reduced nucleoid numbers with significant changes in nucleoid size or shape in fibroblasts from an affected individual compared to controls. We also investigated the effect of tfam impairment in zebrafish; homozygous tfam mutants carrying an in-frame c.141_149 deletion recapitulate the mtDNA depletion and ovarian dysgenesis phenotypes observed in affected humans. Together, our genetic and functional data confirm that TFAM plays a pivotal role in gonad development and expands the repertoire of mitochondrial disease phenotypes., (© 2021. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

The neural crest (NC) is a transient multipotent cell population that migrates extensively to produce a remarkable array of vertebrate cell types. NC cell specification progresses in an anterior to posterior fashion, resulting in distinct, axial-restricted subpopulations. The anterior-most, cranial, population of NC is specified as gastrulation concludes and neurulation begins, while more posterior populations become specified as the body elongates. The mechanisms that govern development of the more posterior NC cells remain incompletely understood. Here, we report a key role for zebrafish Cdx4, a homeodomain transcription factor, in the development of posterior NC cells. We demonstrate that cdx4 is expressed in trunk NC cell progenitors, directly binds NC cell-specific enhancers in the NC GRN, and regulates expression of the key NC development gene foxd3 in the posterior body. Moreover, cdx4 mutants show disruptions to the segmental pattern of trunk NC cell migration due to loss of normal leader/follower cell dynamics. Finally, using cell transplantation to generate chimeric specimens, we show that Cdx4 does not function in the paraxial mesoderm-the environment adjacent to which crest migrates-to influence migratory behaviors. We conclude that cdx4 plays a critical, and likely tissue autonomous, role in the establishment of trunk NC migratory behaviors. Together, our results indicate that cdx4 functions as an early NC specifier gene in the posterior body of zebrafish embryos., Competing Interests: Declaration of competing interest The authors declare no competing financial interests., (Copyright © 2021 Elsevier Inc. All rights reserved.)

The homeodomain-containing transcription factors dlx1a, dlx2a, dlx5a and dlx6a are expressed in the zebrafish brain in overlapping patterns and are important in vertebrate development. Previous work in mice have suggested the overlapping expression pattern is in part due to cross-regulatory interactions among the aforementioned dlx genes. However, the extent of these interactions and whether they are conserved among vertebrates remains to be determined. Through whole-mount in situ hybridization in zebrafish dlx mutants produced by CRISPR-Cas9 mutagenesis, cross-regulatory interactions between dlx1a, dlx2a, dlx5a and dlx6a were examined from 24 to 72 h post fertilization (hpf). Notably, and different from previous work done in mouse, zebrafish dlx2a -/- mutants continue to express dlx5a until 72hpf, whereas deletion of both enhancers within the dlx5a/dlx6a locus resulted in delayed dlx5a/dlx6a expression and relative increased dlx2a expression. These results suggest alternative regulatory elements and pathways exist to mediate dlx expression in zebrafish and may highlight evolutionary differences in gene interactions between vertebrates., (Copyright © 2021 Elsevier B.V. All rights reserved.)

Zygotic genome activation (ZGA) initiates regionalized transcription underlying distinct cellular identities. ZGA is dependent upon dynamic chromatin architecture sculpted by conserved DNA-binding proteins. However, the direct mechanistic link between the onset of ZGA and the tissue-specific transcription remains unclear. Here, we have addressed the involvement of chromatin organizer Satb2 in orchestrating both processes during zebrafish embryogenesis. Integrative analysis of transcriptome, genome-wide occupancy and chromatin accessibility reveals contrasting molecular activities of maternally deposited and zygotically synthesized Satb2. Maternal Satb2 prevents premature transcription of zygotic genes by influencing the interplay between the pluripotency factors. By contrast, zygotic Satb2 activates transcription of the same group of genes during neural crest development and organogenesis. Thus, our comparative analysis of maternal versus zygotic function of Satb2 underscores how these antithetical activities are temporally coordinated and functionally implemented highlighting the evolutionary implications of the biphasic and bimodal regulation of landmark developmental transitions by a single determinant., (© 2021. The Author(s).)

The transcription factor Nkx3.2 (Bapx1) is an important chondrocyte maturation inhibitor. Previous Nkx3.2 knockdown and overexpression studies in non-mammalian gnathostomes have focused on its role in primary jaw joint development, while the function of this gene in broader skeletal development is not fully described. We generated a mutant allele of nkx3.2 in zebrafish with CRISPR/Cas9 and applied a range of techniques to characterize skeletal phenotypes at developmental stages from larva to adult, revealing loss of the jaw joint, fusions in bones of the occiput, morphological changes in the Weberian apparatus, and the loss or deformation of bony elements derived from basiventral cartilages of the vertebrae. Axial phenotypes are reminiscent of Nkx3.2 knockout in mammals, suggesting that the function of this gene in axial skeletal development is ancestral to osteichthyans. Our results highlight the broad role of nkx3.2 in zebrafish skeletal development and its context-specific functions in different skeletal elements., Competing Interests: The authors have declared that no competing interests exist.

Proper function of the vertebrate skeleton requires the development of distinct articulating embryonic cartilages. Irx transcription factors are arranged in co-regulated clusters that are expressed in the developing skeletons of the face and appendages. IrxB cluster genes are required for the separation of toes in mice and formation of the hyoid joint in zebrafish, yet whether Irx genes have broader roles in skeletal development remains unclear. Here, we perform a comprehensive loss-of-function analysis of all 11 Irx genes in zebrafish. We uncover conserved requirements for IrxB genes in formation of the fish and mouse scapula. In the face, we find a requirement for IrxAb genes and irx7 in formation of anterior neural crest precursors of the jaw, and for IrxBa genes in formation of endodermal pouches and gill cartilages. We also observe extensive joint loss and cartilage fusions in animals with combinatorial losses of Irx clusters, with in vivo imaging revealing that at least some of these fusions arise through inappropriate chondrogenesis. Our analysis reveals diverse roles for Irx genes in the formation and later segmentation of the facial skeleton., Competing Interests: Competing interests The authors declare no competing or financial interests., (© 2021. Published by The Company of Biologists Ltd.)

During development, cells respond rapidly to intra- and intercellular signals, which induce signaling cascades regulating the activity of transcription factors at the transcriptional and/or post-translational level. The transcription factor ISL1 plays a key role in second heart field development and cardiac differentiation, and its mRNA levels are tightly regulated during cardiogenesis. Here, we show that a BMP-p38 MAPK signaling axis controls ISL1 protein function at the post-translational level. BMP-mediated activation of p38 MAPK leads to ISL1 phosphorylation at S269 by p38, which prevents ISL1 degradation and ensures its transcriptional activity during cardiogenesis. Interfering with p38 MAPK signaling leads to the degradation of ISL1 by the proteasome, resulting in defects in cardiomyocyte differentiation and impaired zebrafish and mouse heart morphogenesis and function. Given the critical role of the tight control of ISL1 activity during cardiac lineage diversification, modulation of BMP4-p38 MAPK signaling could direct differentiation into specific cardiac cell subpopulations., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: hes genes are chordate homologs of Drosophila genes, hairy and enhancer of split, which encode a basic helix-loop-helix (bHLH) transcriptional repressor with a WRPW motif. Various developmental functions of hes genes, including early embryogenesis and neurogenesis, have been elucidated in vertebrates. However, their orthologous relationships remain unclear partly because of less conservation of relatively short amino acid sequences, the fact that the genome was not analyzed as it is today, and species-specific genome duplication. This results in complicated gene names in vertebrates, which are not consistent in orthologs. We previously revealed that Xenopus frogs have two clusters of hes5, named "the hes5.1 cluster" and "the hes5.3 cluster", but the origin and the conservation have not yet been revealed., Results: Here, we elucidated the orthologous and paralogous relationships of all hes genes of human, mouse, chicken, gecko, zebrafish, medaka, coelacanth, spotted gar, elephant shark and three species of frogs, Xenopus tropicalis (X. tropicalis), X. laevis, Nanorana parkeri, by phylogenetic and synteny analyses. Any duplicated hes5 were not found in mammals, whereas hes5 clusters in teleost were conserved although not as many genes as the three frog species. In addition, hes5 cluster-like structure was found in the elephant shark genome, but not found in cyclostomata., Conclusion: These data suggest that the hes5 cluster existed in the gnathostome ancestor but became a single gene in mammals. The number of hes5 cluster genes were specifically large in frogs., (© 2021. The Author(s).)

Tissue factor pathway inhibitor (TFPI) is an anticoagulant protein that inhibits factor VIIa and Xa in the coagulation cascade. It has been shown that forkhead box P3 protein is a TFPI transcriptional repressor. However, there are no studies on chromatin remodeling that control TFPI expression. We hypothesized that the genome-wide knockdowns of the chromatin binding and regulatory proteins (CBRPs) in zebrafish could identify novel tfpia gene regulators. As an initial step, we selected 69 CBRP genes from the list of zebrafish thrombocyte-expressed genes. We then performed a 3-gene piggyback knockdown screen of these 69 genes, followed by quantification of tfpia mRNA levels. The results revealed that knockdown of brd7, ing2, ing3, ing4, and suz12b increased tfpia mRNA levels. The simultaneous knockdown of these 5 genes also increased tfpia mRNA levels. We also performed individual gene and simultaneous 5-gene knockdowns on the 5 genes in zebrafish larvae. We found that after laser injury, it took a longer time for the formation of the thrombus to occlude the caudal vessel compared to the control larvae. We then treated the larvae and adults with a chemical UNC6852 known to proteolytically degrade polycomb repressor complex 2, where SUZ12 is a member, and observed prolongation of time to occlude (TTO) the caudal vein after laser injury and increased tfpia mRNA levels in larvae and adults, respectively. In summary, our results have identified novel epigenetic regulators for tfpia and exploited this information to discover a drug that enhances tfpia mRNA levels and prolongation of TTO. This discovery provides the basis for testing whether UNC6852 could be used as an antithrombotic drug. This approach could be used to study the regulation of other plasma proteins, including coagulant and anticoagulant factors., (© 2021. The Author(s).)

Cone dystrophies are a rare subgroup of inherited retinal dystrophies and hallmarked by color vision defects, low or decreasing visual acuity and central vision loss, nystagmus and photophobia. Applying genome-wide linkage analysis and array comparative genome hybridization, we identified a locus for autosomal dominant cone dystrophy on chromosome 16q12 in four independent multigeneration families. The locus is defined by duplications of variable size with a smallest region of overlap of 608 kb affecting the IRXB gene cluster and encompasses the genes IRX5 and IRX6. IRX5 and IRX6 belong to the Iroquois (Iro) protein family of homeodomain-containing transcription factors involved in patterning and regionalization of embryonic tissue in vertebrates, including the eye and the retina. All patients presented with a unique progressive cone dystrophy phenotype hallmarked by early tritanopic color vision defects. We propose that the disease underlies a misregulation of the IRXB gene cluster on chromosome 16q12 and demonstrate that overexpression of Irx5a and Irx6a, the two orthologous genes in zebrafish, results in visual impairment in 5-day-old zebrafish larvae., (© The Author(s) 2021. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Osteoporosis is an increasing burden on public health as the world-wide population ages and effective therapeutics are severely needed. Two pathways with high potential for osteoporosis treatment are the retinoic acid (RA) and endocannabinoid system (ECS) signaling pathways. We sought to elucidate the roles that these pathways play in bone development and maturation. Here, we use chemical treatments to modulate the RA and ECS pathways at distinct early, intermediate, and late times bone development in zebrafish. We further assessed osteoclast activity later in zebrafish and medaka. Finally, by combining sub-optimal doses of AR and ECS modulators, we show that enhancing RA signaling or reducing the ECS promote bone formation and decrease osteoclast abundance and activity. These data demonstrate that RA signaling and the ECS can be combined as sub-optimal doses to influence bone growth and may be key targets for potential therapeutics., (Copyright © 2021 Elsevier B.V. All rights reserved.)

Embryos must communicate instructions to their constituent cells over long distances. These instructions are often encoded in the concentration of signals called morphogens. In the textbook view, morphogen molecules diffuse from a localized source to form a concentration gradient, and target cells adopt fates by measuring the local morphogen concentration. However, natural patterning systems often incorporate numerous co-factors and extensive signaling feedback, suggesting that embryos require additional mechanisms to generate signaling patterns. Here, we examine the mechanisms of signaling pattern formation for the mesendoderm inducer Nodal during zebrafish embryogenesis. We find that Nodal signaling activity spans a normal range in the absence of signaling feedback and relay, suggesting that diffusion is sufficient for Nodal gradient formation. We further show that the range of endogenous Nodal ligands is set by the EGF-CFC co-receptor Oep: in the absence of Oep, Nodal activity spreads to form a nearly uniform distribution throughout the embryo. In turn, increasing Oep levels sensitizes cells to Nodal ligands. We recapitulate these experimental results with a computational model in which Oep regulates the diffusive spread of Nodal ligands by setting the rate of capture by target cells. This model predicts, and we confirm in vivo, the surprising observation that a failure to replenish Oep transforms the Nodal signaling gradient into a travelling wave. These results reveal that patterns of Nodal morphogen signaling are shaped by co-receptor-mediated restriction of ligand spread and sensitization of responding cells., Competing Interests: NL, AC, PA, AS No competing interests declared, (© 2021, Lord et al.)

Transcription and metabolism both influence cell function, but dedicated transcriptional control of metabolic pathways that regulate cell fate has rarely been defined. We discovered, using a chemical suppressor screen, that inhibition of the pyrimidine biosynthesis enzyme dihydroorotate dehydrogenase (DHODH) rescues erythroid differentiation in bloodless zebrafish moonshine (mon) mutant embryos defective for transcriptional intermediary factor 1 gamma ( tif1γ ). This rescue depends on the functional link of DHODH to mitochondrial respiration. The transcription elongation factor TIF1γ directly controls coenzyme Q (CoQ) synthesis gene expression. Upon tif1γ loss, CoQ levels are reduced, and a high succinate/α-ketoglutarate ratio leads to increased histone methylation. A CoQ analog rescues mon 's bloodless phenotype. These results demonstrate that mitochondrial metabolism is a key output of a lineage transcription factor that drives cell fate decisions in the early blood lineage., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Precise control of gene expression is crucial to ensure proper development and biological functioning of an organism. Enhancers are non-coding DNA elements which play an essential role in regulating gene expression. They contain specific sequence motifs serving as binding sites for transcription factors which interact with the basal transcription machinery at their target genes. Heart development is regulated by intricate gene regulatory network ensuring precise spatiotemporal gene expression program. Mutations affecting enhancers have been shown to result in devastating forms of congenital heart defect. Therefore, identifying enhancers implicated in heart biology and understanding their mechanism is key to improve diagnosis and therapeutic options. Despite their crucial role, enhancers are poorly studied, mainly due to a lack of reliable way to identify them and determine their function. Nevertheless, recent technological advances have allowed rapid progress in enhancer discovery. Model organisms such as the zebrafish have contributed significant insights into the genetics of heart development through enabling functional analyses of genes and their regulatory elements in vivo. Here, we summarize the current state of knowledge on heart enhancers gained through studies in model organisms, discuss various approaches to discover and study their function, and finally suggest methods that could further advance research in this field.

Aims/hypothesis: Variants close to the VPS13C/C2CD4A/C2CD4B locus are associated with altered risk of type 2 diabetes in genome-wide association studies. While previous functional work has suggested roles for VPS13C and C2CD4A in disease development, none has explored the role of C2CD4B., Methods: CRISPR/Cas9-induced global C2cd4b-knockout mice and zebrafish larvae with c2cd4a deletion were used to study the role of this gene in glucose homeostasis. C2 calcium dependent domain containing protein (C2CD)4A and C2CD4B constructs tagged with FLAG or green fluorescent protein were generated to investigate subcellular dynamics using confocal or near-field microscopy and to identify interacting partners by mass spectrometry., Results: Systemic inactivation of C2cd4b in mice led to marked, but highly sexually dimorphic changes in body weight and glucose homeostasis. Female C2cd4b mice displayed unchanged body weight compared with control littermates, but abnormal glucose tolerance (AUC, p = 0.01) and defective in vivo, but not in vitro, insulin secretion (p = 0.02). This was associated with a marked decrease in follicle-stimulating hormone levels as compared with wild-type (WT) littermates (p = 0.003). In sharp contrast, male C2cd4b null mice displayed essentially normal glucose tolerance but an increase in body weight (p < 0.001) and fasting blood glucose (p = 0.003) after maintenance on a high-fat and -sucrose diet vs WT littermates. No metabolic disturbances were observed after global inactivation of C2cd4a in mice, or in pancreatic beta cell function at larval stages in C2cd4a null zebrafish. Fasting blood glucose levels were also unaltered in adult C2cd4a-null fish. C2CD4B and C2CD4A were partially localised to the plasma membrane, with the latter under the control of intracellular Ca 2+ . Binding partners for both included secretory-granule-localised PTPRN2/phogrin., Conclusions/interpretation: Our studies suggest that C2cd4b may act centrally in the pituitary to influence sex-dependent circuits that control pancreatic beta cell function and glucose tolerance in rodents. However, the absence of sexual dimorphism in the impact of diabetes risk variants argues for additional roles for C2CD4A or VPS13C in the control of glucose homeostasis in humans., Data Availability: The datasets generated and/or analysed during the current study are available in the Biorxiv repository ( www.biorxiv.org/content/10.1101/2020.05.18.099200v1 ). RNA-Seq (GSE152576) and proteomics (PXD021597) data have been deposited to GEO ( www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE152576 ) and ProteomeXchange ( www.ebi.ac.uk/pride/archive/projects/PXD021597 ) repositories, respectively.

In this study, we elucidate a single cell resolution fate map in the zebrafish in a sub-section of the anterior Lateral Plate Mesoderm (aLPM) at 18 hpf. Our results show that this tissue is not organized into segregated regions but gives rise to intermingled pericardial sac, peritoneum, pharyngeal arch and cardiac precursors. We further report upon asymmetrical contributions of lateral aLPM-derived heart precursors-specifically that twice as many heart precursors arise from the right side versus the left side of the embryo. Cell tracking analyses and large-scale cell labeling of the lateral aLPM corroborate these differences and show that the observed asymmetries are dependent upon Tbx5a expression. Previously, it was shown that cardiac looping was affected in Tbx5a knock-down and knock-out zebrafish (Garrity et al., 2002; Parrie et al., 2013); our present data also implicate tbx5a function in cell specification, establishment and maintenance of cardiac left-right asymmetry., (Copyright © 2021. Published by Elsevier Inc.)

The homeodomain transcription factor SHOX2 is involved in the development and function of the heart's primary pacemaker, the sinoatrial node (SAN), and has been associated with cardiac conduction-related diseases such as atrial fibrillation and sinus node dysfunction. To shed light on Shox2-dependent genetic processes involved in these diseases, we established a murine embryonic stem cell (ESC) cardiac differentiation model to investigate Shox2 pathways in SAN-like cardiomyocytes. Differential RNA-seq-based expression profiling of Shox2 +/+ and Shox2 -/- ESCs revealed 94 dysregulated transcripts in Shox2 -/- ESC-derived SAN-like cells. Of these, 15 putative Shox2 target genes were selected for further validation based on comparative expression analysis with SAN- and right atria-enriched genes. Network-based analyses, integrating data from the Mouse Organogenesis Cell Atlas and the Ingenuity pathways, as well as validation in mouse and zebrafish models confirmed a regulatory role for the novel identified Shox2 target genes including Cav1, Fkbp10, Igfbp5, Mcf2l and Nr2f2. Our results indicate that genetic networks involving SHOX2 may contribute to conduction traits through the regulation of these genes., (Copyright © 2021. Published by Elsevier B.V.)

Cell types are the basic building units of multicellular life, with extensive diversities. The evolution of cell types is a crucial layer of comparative cell biology but is thus far not comprehensively studied. We define a compendium of cell atlases using single-cell RNA-seq (scRNA-seq) data from seven animal species and construct a cross-species cell-type evolutionary hierarchy. We present a roadmap for the origin and diversity of major cell categories and find that muscle and neuron cells are conserved cell types. Furthermore, we identify a cross-species transcription factor (TF) repertoire that specifies major cell categories. Overall, our study reveals conservation and divergence of cell types during animal evolution, which will further expand the landscape of comparative genomics., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Genetic variants associated with type 2 diabetes (T2D) risk affect gene regulation in metabolically relevant tissues, such as pancreatic islets. Here, we investigated contributions of regulatory programs active during pancreatic development to T2D risk. Generation of chromatin maps from developmental precursors throughout pancreatic differentiation of human embryonic stem cells (hESCs) identifies enrichment of T2D variants in pancreatic progenitor-specific stretch enhancers that are not active in islets. Genes associated with progenitor-specific stretch enhancers are predicted to regulate developmental processes, most notably tissue morphogenesis. Through gene editing in hESCs, we demonstrate that progenitor-specific enhancers harboring T2D-associated variants regulate cell polarity genes LAMA1 and CRB2 . Knockdown of lama1 or crb2 in zebrafish embryos causes a defect in pancreas morphogenesis and impairs islet cell development. Together, our findings reveal that a subset of T2D risk variants specifically affects pancreatic developmental programs, suggesting that dysregulation of developmental processes can predispose to T2D., Competing Interests: RG, AW, JC, JL, NW, SK, JW, YQ, JY, AA, BR, PD, MS No competing interests declared, KG This author consults for Genentech., (© 2021, Geusz et al.)

The regulated expansion of chondrocytes within growth plates and joints ensures proper skeletal development through adulthood. Mutations in the transcription factor NKX3.2 underlie spondylo-megaepiphyseal-metaphyseal dysplasia (SMMD), which is characterized by skeletal defects including scoliosis, large epiphyses, wide growth plates and supernumerary distal limb joints. Whereas nkx3.2 knockdown zebrafish and mouse Nkx3.2 mutants display embryonic lethal jaw joint fusions and skeletal reductions, respectively, they lack the skeletal overgrowth seen in SMMD patients. Here, we report adult viable nkx3.2 mutant zebrafish displaying cartilage overgrowth in place of a missing jaw joint, as well as severe dysmorphologies of the facial skeleton, skullcap and spine. In contrast, cartilage overgrowth and scoliosis are absent in rare viable nkx3.2 knockdown animals that lack jaw joints, supporting post-embryonic roles for Nkx3.2. Single-cell RNA-sequencing and in vivo validation reveal increased proliferation and upregulation of stress-induced pathways, including prostaglandin synthases, in mutant chondrocytes. By generating a zebrafish model for the skeletal overgrowth defects of SMMD, we reveal post-embryonic roles for Nkx3.2 in dampening proliferation and buffering the stress response in joint-associated chondrocytes., Competing Interests: Competing interestsThe authors declare no competing or financial interests., (© 2021. Published by The Company of Biologists Ltd.)

LIM homeobox 4 (LHX4) is expressed in the photoreceptors (PRs) of the outer nuclear layer (ONL) and bipolar cells (BCs) of the inner nuclear layer (INL) in mouse and chicken retina. It regulates the subtype-specific development of rod BCs and cone BCs in the mouse retina. However, no report has been published on its expression and function in the zebrafish retina. In this study, we assessed the expression of Lhx4 using in situ hybridization (ISH) technique and explored its role in zebrafish (Danio rerio) retinal development via morpholino (MO) technology. We found that the expression of lhx4 in the zebrafish retina begins 48 h post-fertilization (hpf) and is continuously expressed in the ONL and INL. A zebrafish model constructed with lhx4 knockdown in the eyes through vivo-MO revealed that: lhx4 knockdown inhibits the differentiation of Parvalbumin + amacrine cells (ACs) and Rhodopsin + rod photoreceptors (RPs), enhances the expression of visual system homeobox 2 (vsx2); and damages the responses of zebrafish to light stimulus, without affecting the differentiation of OFF-BCs and rod BCs, and apoptosis in the retina. These findings reveal that lhx4 regulates neural differentiation in the retina and visual function during zebrafish embryonic development.

Rax ( Rx ) genes encode paired-type homeodomain-containing transcription factors present in virtually all metazoan groups. In vertebrates, studies in fish, amphibian, chick and mouse models have revealed that these genes play important roles in the development of structures located at the anterior portion of the central nervous system, in particular the eyes, the hypothalamus and the pituitary gland. In addition, human patients with eye and brain defects carry mutations in the two human Rax paralogues, RAX and RAX2 . Here, we review work done in the last years on Rax genes, focusing especially on the function that mouse Rax and its zebrafish homologue, rx3 , play in hypothalamic and pituitary development. Work on both of these model organisms indicate that Rax genes are necessary for the patterning, growth and differentiation of the hypothalamus, in particular the ventro-tuberal and dorso-anterior hypothalamus, where they effect their action by controlling expression of the secreted signalling protein, Sonic hedgehog (Shh). In addition, Rax/rx3 mutations disturb the development of the pituitary gland, mimicking phenotypes observed in human subjects carrying mutations in the RAX gene. Thus, along with their crucial role in eye morphogenesis, Rax genes play a conserved role in the development of the hypothalamus and adjacent structures in the vertebrate clade.

Gene expression is an inherently stochastic process 1,2 ; however, organismal development and homeostasis require cells to coordinate the spatiotemporal expression of large sets of genes. In metazoans, pairs of co-expressed genes often reside in the same chromosomal neighbourhood, with gene pairs representing 10 to 50% of all genes, depending on the species 3-6 . Because shared upstream regulators can ensure correlated gene expression, the selective advantage of maintaining adjacent gene pairs remains unknown 6 . Here, using two linked zebrafish segmentation clock genes, her1 and her7, and combining single-cell transcript counting, genetic engineering, real-time imaging and computational modelling, we show that gene pairing boosts correlated transcription and provides phenotypic robustness for the formation of developmental patterns. Our results demonstrate that the prevention of gene pairing disrupts oscillations and segmentation, and the linkage of her1 and her7 is essential for the development of the body axis in zebrafish embryos. We predict that gene pairing may be similarly advantageous in other organisms, and our findings could lead to the engineering of precise synthetic clocks in embryos and organoids.