Multiomic atlas with functional stratification and developmental dynamics of zebrafish cis-regulatory elements.

Zebrafish, a popular organism for studying embryonic development and for modeling human diseases, has so far lacked a systematic functional annotation program akin to those in other animal models. To address this, we formed the international DANIO-CODE consortium and created a central repository to store and process zebrafish developmental functional genomic data. Our data coordination center ( https://danio-code.zfin.org ) combines a total of 1,802 sets of unpublished and re-analyzed published genomic data, which we used to improve existing annotations and show its utility in experimental design. We identified over 140,000 cis-regulatory elements throughout development, including classes with distinct features dependent on their activity in time and space. We delineated the distinct distance topology and chromatin features between regulatory elements active during zygotic genome activation and those active during organogenesis. Finally, we matched regulatory elements and epigenomic landscapes between zebrafish and mouse and predicted functional relationships between them beyond sequence similarity, thus extending the utility of zebrafish developmental genomics to mammals.

Identification of an evolutionarily conserved domain in Neurod1 favouring enteroendocrine versus goblet cell fate.

ARP/ASCL transcription factors are key determinants of cell fate specification in a wide variety of tissues, coordinating the acquisition of generic cell fates and of specific subtype identities. How these factors, recognizing highly similar DNA motifs, display specific activities, is not yet fully understood. To address this issue, we overexpressed different ARP/ASCL factors in zebrafish ascl1a-/- mutant embryos to determine which ones are able to rescue the intestinal secretory lineage. We found that Ascl1a/b, Atoh1a/b and Neurod1 factors are all able to trigger the first step of the secretory regulatory cascade but distinct secretory cells are induced by these factors. Indeed, Neurod1 rescues the enteroendocrine lineage while Ascl1a/b and Atoh1a/b rescue the goblet cells. Gain-of-function experiments with Ascl1a/Neurod1 chimeric proteins revealed that the functional divergence is encoded by a 19-aa ultra-conserved element (UCE), present in all Neurod members but absent in the other ARP/ASCL proteins. Importantly, inserting the UCE into the Ascl1a protein reverses the rescuing capacity of this Ascl1a chimeric protein that cannot rescue the goblet cells anymore but can efficiently rescue the enteroendocrine cells. This novel domain acts indeed as a goblet cell fate repressor that inhibits gfi1aa expression, known to be important for goblet cell differentiation. Deleting the UCE domain of the endogenous Neurod1 protein leads to an increase in the number of goblet cells concomitant with a reduction of the enteroendocrine cells, phenotype also observed in the neurod1 null mutant. This highlights the crucial function of the UCE domain for NeuroD1 activity in the intestine. As Gfi1 acts as a binary cell fate switch in several tissues where Neurod1 is also expressed, we can envision a similar role of the UCE in other tissues, allowing Neurod1 to repress Gfi1 to influence the balance between cell fates.

A δ-cell subpopulation with a pro-ß-cell identity contributes to efficient age-independent recovery in a zebrafish model of diabetes.

Restoring damaged ß-cells in diabetic patients by harnessing the plasticity of other pancreatic cells raises the questions of the efficiency of the process and of the functionality of the new Insulin-expressing cells. To overcome the weak regenerative capacity of mammals, we used regeneration-prone zebrafish to study ß-cells arising following destruction. We show that most new insulin cells differ from the original ß-cells as they coexpress Somatostatin and Insulin. These bihormonal cells are abundant, functional and able to normalize glycemia. Their formation in response to ß-cell destruction is fast, efficient, and age-independent. Bihormonal cells are transcriptionally close to a subset of δ-cells that we identified in control islets and that are characterized by the expression of somatostatin 1.1 (sst1.1) and by genes essential for glucose-induced Insulin secretion in ß-cells such as pdx1, slc2a2 and gck. We observed in vivo the conversion of monohormonal sst1.1-expressing cells to sst1.1+ ins + bihormonal cells following ß-cell destruction. Our findings support the conclusion that sst1.1 δ-cells possess a pro-ß identity enabling them to contribute to the neogenesis of Insulin-producing cells during regeneration. This work unveils that abundant and functional bihormonal cells benefit to diabetes recovery in zebrafish.

Identification of downstream effectors of retinoic acid specifying the zebrafish pancreas by integrative genomics.

Retinoic acid (RA) is a key signal for the specification of the pancreas. Still, the gene regulatory cascade triggered by RA in the endoderm remains poorly characterized. In this study, we investigated this regulatory network in zebrafish by combining RNA-seq, RAR ChIP-seq and ATAC-seq assays. By analysing the effect of RA and of the RA receptor (RAR) inverse-agonist BMS493 on the transcriptome and on the chromatin accessibility of endodermal cells, we identified a large set of genes and regulatory regions regulated by RA signalling. RAR ChIP-seq further defined the direct RAR target genes in zebrafish, including hox genes as well as several pancreatic regulators like mnx1, insm1b, hnf1ba and gata6. Comparison of zebrafish and murine RAR ChIP-seq data highlighted the conserved direct target genes and revealed that some RAR sites are under strong evolutionary constraints. Among them, a novel highly conserved RAR-induced enhancer was identified downstream of the HoxB locus and driving expression in the nervous system and in the gut in a RA-dependent manner. Finally, ATAC-seq data unveiled the role of the RAR-direct targets Hnf1ba and Gata6 in opening chromatin at many regulatory loci upon RA treatment.

Pancreatic and intestinal endocrine cells in zebrafish share common transcriptomic signatures and regulatory programmes.

BACKGROUND: Endocrine cells of the zebrafish digestive system play an important role in regulating metabolism and include pancreatic endocrine cells (PECs) clustered in the islets of Langerhans and the enteroendocrine cells (EECs) scattered in the intestinal epithelium. Despite EECs and PECs are being located in distinct organs, their differentiation involves shared molecular mechanisms and transcription factors. However, their degree of relatedness remains unexplored. In this study, we investigated comprehensively the similarity of EECs and PECs by defining their transcriptomic landscape and comparing the regulatory programmes controlled by Pax6b, a key player in both EEC and PEC differentiations. RESULTS: RNA sequencing was performed on EECs and PECs isolated from wild-type and pax6b mutant zebrafish. Data mining of wild-type zebrafish EEC data confirmed the expression of orthologues for most known mammalian EEC hormones, but also revealed the expression of three additional neuropeptide hormones (Proenkephalin-a, Calcitonin-a and Adcyap1a) not previously reported to be expressed by EECs in any species. Comparison of transcriptomes from EECs, PECs and other zebrafish tissues highlights a very close similarity between EECs and PECs, with more than 70% of genes being expressed in both endocrine cell types. Comparison of Pax6b-regulated genes in EECs and PECs revealed a significant overlap. pax6b loss-of-function does not affect the total number of EECs and PECs but instead disrupts the balance between endocrine cell subtypes, leading to an increase of ghrelin- and motilin-like-expressing cells in both the intestine and pancreas at the expense of other endocrine cells such as beta and delta cells in the pancreas and pyyb-expressing cells in the intestine. Finally, we show that the homeodomain of Pax6b is dispensable for its action in both EECs and PECs. CONCLUSION: We have analysed the transcriptomic landscape of wild-type and pax6b mutant zebrafish EECs and PECs. Our study highlights the close relatedness of EECs and PECs at the transcriptomic and regulatory levels, supporting the hypothesis of a common phylogenetic origin and underscoring the potential implication of EECs in metabolic diseases such as type 2 diabetes.

Dual-initiation promoters with intertwined canonical and TCT/TOP transcription start sites diversify transcript processing.

Variations in transcription start site (TSS) selection reflect diversity of preinitiation complexes and can impact on post-transcriptional RNA fates. Most metazoan polymerase II-transcribed genes carry canonical initiation with pyrimidine/purine (YR) dinucleotide, while translation machinery-associated genes carry polypyrimidine initiator (5'-TOP or TCT). By addressing the developmental regulation of TSS selection in zebrafish we uncovered a class of dual-initiation promoters in thousands of genes, including snoRNA host genes. 5'-TOP/TCT initiation is intertwined with canonical initiation and used divergently in hundreds of dual-initiation promoters during maternal to zygotic transition. Dual-initiation in snoRNA host genes selectively generates host and snoRNA with often different spatio-temporal expression. Dual-initiation promoters are pervasive in human and fruit fly, reflecting evolutionary conservation. We propose that dual-initiation on shared promoters represents a composite promoter architecture, which can function both coordinately and divergently to diversify RNAs.

Codon-specific translation reprogramming promotes resistance to targeted therapy.

Reprogramming of mRNA translation has a key role in cancer development and drug resistance 1 . However, the molecular mechanisms that are involved in this process remain poorly understood. Wobble tRNA modifications are required for specific codon decoding during translation2,3. Here we show, in humans, that the enzymes that catalyse modifications of wobble uridine 34 (U34) tRNA (U34 enzymes) are key players of the protein synthesis rewiring that is induced by the transformation driven by the BRAF V600E oncogene and by resistance to targeted therapy in melanoma. We show that BRAF V600E -expressing melanoma cells are dependent on U34 enzymes for survival, and that concurrent inhibition of MAPK signalling and ELP3 or CTU1 and/or CTU2 synergizes to kill melanoma cells. Activation of the PI3K signalling pathway, one of the most common mechanisms of acquired resistance to MAPK therapeutic agents, markedly increases the expression of U34 enzymes. Mechanistically, U34 enzymes promote glycolysis in melanoma cells through the direct, codon-dependent, regulation of the translation of HIF1A mRNA and the maintenance of high levels of HIF1α protein. Therefore, the acquired resistance to anti-BRAF therapy is associated with high levels of U34 enzymes and HIF1α. Together, these results demonstrate that U34 enzymes promote the survival and resistance to therapy of melanoma cells by regulating specific mRNA translation.

Nifurpirinol: A more potent and reliable substrate compared to metronidazole for nitroreductase-mediated cell ablations.

The zebrafish is a popular animal model with well-known regenerative capabilities. To study regeneration in this fish, the nitroreductase/metronidazole-mediated system is widely used for targeted ablation of various cell types. Nevertheless, we highlight here some variability in ablation efficiencies with the metronidazole prodrug that led us to search for a more efficient and reliable compound. Herein, we present nifurpirinol, another nitroaromatic antibiotic, as a more potent prodrug compared to metronidazole to trigger cell-ablation in nitroreductase expressing transgenic models. We show that nifurpirinol induces robust and reliable ablations at concentrations 2,000 fold lower than metronidazole and three times below its own toxic concentration. We confirmed the efficiency of nifurpirinol in triggering massive ablation of three different cell types: the pancreatic beta cells, osteoblasts, and dopaminergic neurons. Our results identify nifurpirinol as a very potent prodrug for the nitroreductase-mediated ablation system and suggest that its use could be extended to many other cell types, especially if difficult to ablate, or when combined pharmacological treatments are desired.

Tcf7l2 plays pleiotropic roles in the control of glucose homeostasis, pancreas morphology, vascularization and regeneration.

Type 2 diabetes (T2D) is a disease characterized by impaired insulin secretion. The Wnt signaling transcription factor Tcf7l2 is to date the T2D-associated gene with the largest effect on disease susceptibility. However, the mechanisms by which TCF7L2 variants affect insulin release from ß-cells are not yet fully understood. By taking advantage of a tcf7l2 zebrafish mutant line, we first show that these animals are characterized by hyperglycemia and impaired islet development. Moreover, we demonstrate that the zebrafish tcf7l2 gene is highly expressed in the exocrine pancreas, suggesting potential bystander effects on ß-cell growth, differentiation and regeneration. Finally, we describe a peculiar vascular phenotype in tcf7l2 mutant larvae, characterized by significant reduction in the average number and diameter of pancreatic islet capillaries. Overall, the zebrafish Tcf7l2 mutant, characterized by hyperglycemia, pancreatic and vascular defects, and reduced regeneration proves to be a suitable model to study the mechanism of action and the pleiotropic effects of Tcf7l2, the most relevant T2D GWAS hit in human populations.

Transcriptome analysis of pancreatic cells across distant species highlights novel important regulator genes.

BACKGROUND: Defining the transcriptome and the genetic pathways of pancreatic cells is of great interest for elucidating the molecular attributes of pancreas disorders such as diabetes and cancer. As the function of the different pancreatic cell types has been maintained during vertebrate evolution, the comparison of their transcriptomes across distant vertebrate species is a means to pinpoint genes under strong evolutionary constraints due to their crucial function, which have therefore preserved their selective expression in these pancreatic cell types. RESULTS: In this study, RNA-sequencing was performed on pancreatic alpha, beta, and delta endocrine cells as well as the acinar and ductal exocrine cells isolated from adult zebrafish transgenic lines. Comparison of these transcriptomes identified many novel markers, including transcription factors and signaling pathway components, specific for each cell type. By performing interspecies comparisons, we identified hundreds of genes with conserved enriched expression in endocrine and exocrine cells among human, mouse, and zebrafish. This list includes many genes known as crucial for pancreatic cell formation or function, but also pinpoints many factors whose pancreatic function is still unknown. A large set of endocrine-enriched genes can already be detected at early developmental stages as revealed by the transcriptomic profiling of embryonic endocrine cells, indicating a potential role in cell differentiation. The actual involvement of conserved endocrine genes in pancreatic cell differentiation was demonstrated in zebrafish for myt1b, whose invalidation leads to a reduction of alpha cells, and for cdx4, selectively expressed in endocrine delta cells and crucial for their specification. Intriguingly, comparison of the endocrine alpha and beta cell subtypes from human, mouse, and zebrafish reveals a much lower conservation of the transcriptomic signatures for these two endocrine cell subtypes compared to the signatures of pan-endocrine and exocrine cells. These data suggest that the identity of the alpha and beta cells relies on a few key factors, corroborating numerous examples of inter-conversion between these two endocrine cell subtypes. CONCLUSION: This study highlights both evolutionary conserved and species-specific features that will help to unveil universal and fundamental regulatory pathways as well as pathways specific to human and laboratory animal models such as mouse and zebrafish.

Habenular Neurogenesis in Zebrafish Is Regulated by a Hedgehog, Pax6 Proneural Gene Cascade.

The habenulae are highly conserved nuclei in the dorsal diencephalon that connect the forebrain to the midbrain and hindbrain. These nuclei have been implicated in a broad variety of behaviours in humans, primates, rodents and zebrafish. Despite this, the molecular mechanisms that control the genesis and differentiation of neural progenitors in the habenulae remain relatively unknown. We have previously shown that, in zebrafish, the timing of habenular neurogenesis is left-right asymmetric and that in the absence of Nodal signalling this asymmetry is lost. Here, we show that habenular neurogenesis requires the homeobox transcription factor Pax6a and the redundant action of two proneural bHLH factors, Neurog1 and Neurod4. We present evidence that Hedgehog signalling is required for the expression of pax6a, which is in turn necessary for the expression of neurog1 and neurod4. Finally, we demonstrate by pharmacological inhibition that Hedgehog signalling is required continuously during habenular neurogenesis and by cell transplantation experiments that pathway activation is required cell autonomously. Our data sheds light on the mechanism underlying habenular development that may provide insights into how Nodal signalling imposes asymmetry on the timing of habenular neurogenesis.

Phenotypic and biomarker evaluation of zebrafish larvae as an alternative model to predict mammalian hepatotoxicity.

Zebrafish phenotypic assays have shown promise to assess human hepatotoxicity, though scoring of liver morphology remains subjective and difficult to standardize. Liver toxicity in zebrafish larvae at 5 days was assessed using gene expression as the biomarker approach, complementary to phenotypic analysis and analytical data on compound uptake. This approach aimed to contribute to improved hepatotoxicity prediction, with the goal of identifying biomarker(s) as a step towards the development of transgenic models for prioritization. Morphological effects of hepatotoxic compounds (acetaminophen, amiodarone, coumarin, methapyrilene and myclobutanil) and saccharin as the negative control were assessed after exposure in zebrafish larvae. The hepatotoxic compounds induced the expected zebrafish liver degeneration or changes in size, whereas saccharin did not have any phenotypic adverse effect. Analytical methods based on liquid chromatography-mass spectrometry were optimized to measure stability of selected compounds in exposure medium and internal concentration in larvae. All compounds were stable, except amiodarone for which precipitation was observed. There was a wide variation between the levels of compound in the zebrafish larvae with a higher uptake of amiodarone, methapyrilene and myclobutanil. Detection of hepatocyte markers (CP, CYP3A65, GC and TF) was accomplished by in situ hybridization of larvae to coumarin and myclobutanil and confirmed by real-time reverse transcription-quantitative polymerase chain reaction. Experiments showed decreased expression of all markers. Next, other liver-specific biomarkers (i.e. FABP10a and NR1H4) and apoptosis (i.e. CASP-3 A and TP53) or cytochrome P450-related (CYP2K19) and oxidoreductase activity-related (ZGC163022) genes, were screened. Links between basic mechanisms of liver injury and results of biomarker responses are described. Copyright © 2016 John Wiley & Sons, Ltd.

ADAMTS3 activity is mandatory for embryonic lymphangiogenesis and regulates placental angiogenesis.

The only documented activity of a subclass of ADAMTS proteases comprising ADAMTS2, 3 and 14 is the cleavage of the aminopropeptide of fibrillar procollagens. A limited number of in vitro studies suggested that ADAMTS3 is mainly responsible for procollagen II processing in cartilage. Here, we created an ADAMTS3 knockout mouse (Adamts3(-/-)) model to determine in vivo the actual functions of ADAMTS3. Heterozygous Adamts3(+/-) mice were viable and fertile, but their intercrosses demonstrated lethality of Adamts3(-/-) embryos after 15 days of gestation. Procollagens I, II and III processing was unaffected in these embryos. However, a massive lymphedema caused by the lack of lymphatics development, an abnormal blood vessel structure in the placenta and a progressive liver destruction were observed. These phenotypes are most probably linked to dysregulation of the VEGF-C pathways. This study is the first demonstration that an aminoprocollagen peptidase is crucial for developmental processes independently of its primary role in collagen biology and has physiological functions potentially involved in several human diseases related to angiogenesis and lymphangiogenesis.

Copy number variants in patients with intellectual disability affect the regulation of ARX transcription factor gene.

Protein-coding mutations in the transcription factor-encoding gene ARX cause various forms of intellectual disability (ID) and epilepsy. In contrast, variations in surrounding non-coding sequences are correlated with milder forms of non-syndromic ID and autism and had suggested the importance of ARX gene regulation in the etiology of these disorders. We compile data on several novel and some already identified patients with or without ID that carry duplications of ARX genomic region and consider likely genetic mechanisms underlying the neurodevelopmental defects. We establish the long-range regulatory domain of ARX and identify its brain region-specific autoregulation. We conclude that neurodevelopmental disturbances in the patients may not simply arise from increased dosage due to ARX duplication. This is further exemplified by a small duplication involving a non-functional ARX copy, but with duplicated enhancers. ARX enhancers are located within a 504-kb region and regulate expression specifically in the forebrain in developing and adult zebrafish. Transgenic enhancer-reporter lines were used as in vivo tools to delineate a brain region-specific negative and positive autoregulation of ARX. We find autorepression of ARX in the telencephalon and autoactivation in the ventral thalamus. Fluorescently labeled brain regions in the transgenic lines facilitated the identification of neuronal outgrowth and pathfinding disturbances in the ventral thalamus and telencephalon that occur when arxa dosage is diminished. In summary, we have established a model for how breakpoints in long-range gene regulation alter the expression levels of a target gene brain region-specifically, and how this can cause subtle neuronal phenotypes relating to the etiology of associated neuropsychiatric disease.

Progenitor potential of nkx6.1-expressing cells throughout zebrafish life and during beta cell regeneration.

BACKGROUND: In contrast to mammals, the zebrafish has the remarkable capacity to regenerate its pancreatic beta cells very efficiently. Understanding the mechanisms of regeneration in the zebrafish and the differences with mammals will be fundamental to discovering molecules able to stimulate the regeneration process in mammals. To identify the pancreatic cells able to give rise to new beta cells in the zebrafish, we generated new transgenic lines allowing the tracing of multipotent pancreatic progenitors and endocrine precursors. RESULTS: Using novel bacterial artificial chromosome transgenic nkx6.1 and ascl1b reporter lines, we established that nkx6.1-positive cells give rise to all the pancreatic cell types and ascl1b-positive cells give rise to all the endocrine cell types in the zebrafish embryo. These two genes are initially co-expressed in the pancreatic primordium and their domains segregate, not as a result of mutual repression, but through the opposite effects of Notch signaling, maintaining nkx6.1 expression while repressing ascl1b in progenitors. In the adult zebrafish, nkx6.1 expression persists exclusively in the ductal tree at the tip of which its expression coincides with Notch active signaling in centroacinar/terminal end duct cells. Tracing these cells reveals that they are able to differentiate into other ductal cells and into insulin-expressing cells in normal (non-diabetic) animals. This capacity of ductal cells to generate endocrine cells is supported by the detection of ascl1b in the nkx6.1:GFP ductal cell transcriptome. This transcriptome also reveals, besides actors of the Notch and Wnt pathways, several novel markers such as id2a. Finally, we show that beta cell ablation in the adult zebrafish triggers proliferation of ductal cells and their differentiation into insulin-expressing cells. CONCLUSIONS: We have shown that, in the zebrafish embryo, nkx6.1+ cells are bona fide multipotent pancreatic progenitors, while ascl1b+ cells represent committed endocrine precursors. In contrast to the mouse, pancreatic progenitor markers nkx6.1 and pdx1 continue to be expressed in adult ductal cells, a subset of which we show are still able to proliferate and undergo ductal and endocrine differentiation, providing robust evidence of the existence of pancreatic progenitor/stem cells in the adult zebrafish. Our findings support the hypothesis that nkx6.1+ pancreatic progenitors contribute to beta cell regeneration. Further characterization of these cells will open up new perspectives for anti-diabetic therapies.

Ascl1b and Neurod1, instead of Neurog3, control pancreatic endocrine cell fate in zebrafish.

BACKGROUND: NEUROG3 is a key regulator of pancreatic endocrine cell differentiation in mouse, essential for the generation of all mature hormone producing cells. It is repressed by Notch signaling that prevents pancreatic cell differentiation by maintaining precursors in an undifferentiated state. RESULTS: We show that, in zebrafish, neurog3 is not expressed in the pancreas and null neurog3 mutant embryos do not display any apparent endocrine defects. The control of endocrine cell fate is instead fulfilled by two basic helix-loop-helix factors, Ascl1b and Neurod1, that are both repressed by Notch signaling. ascl1b is transiently expressed in the mid-trunk endoderm just after gastrulation and is required for the generation of the first pancreatic endocrine precursor cells. Neurod1 is expressed afterwards in the pancreatic anlagen and pursues the endocrine cell differentiation program initiated by Ascl1b. Their complementary role in endocrine differentiation of the dorsal bud is demonstrated by the loss of all hormone-secreting cells following their simultaneous inactivation. This defect is due to a blockage of the initiation of endocrine cell differentiation. CONCLUSIONS: This study demonstrates that NEUROG3 is not the unique pancreatic endocrine cell fate determinant in vertebrates. A general survey of endocrine cell fate determinants in the whole digestive system among vertebrates indicates that they all belong to the ARP/ASCL family but not necessarily to the Neurog3 subfamily. The identity of the ARP/ASCL factor involved depends not only on the organ but also on the species. One could, therefore, consider differentiating stem cells into insulin-producing cells without the involvement of NEUROG3 but via another ARP/ASCL factor.

The bHLH transcription factor Ascl1a is essential for the specification of the intestinal secretory cells and mediates Notch signaling in the zebrafish intestine.

Notch signaling has a fundamental role in stem cell maintenance and in cell fate choice in the intestine of different species. Canonically, Notch signaling represses the expression of transcription factors of the achaete-scute like (ASCL) or atonal related protein (ARP) families. Identifying the ARP/ASCL genes expressed in the gastrointestinal tract is essential to build the regulatory cascade controlling the differentiation of gastrointestinal progenitors into the different intestinal cell types. The expression of the ARP/ASCL factors was analyzed in zebrafish to identify, among all the ARP/ASCL factors found in the zebrafish genome, those expressed in the gastrointestinal tract. ascl1a was found to be the earliest factor detected in the intestine. Loss-of-function analyses using the pia/ascl1a mutant, revealed that ascl1a is crucial for the differentiation of all secretory cells. Furthermore, we identify a battery of transcription factors expressed during secretory cell differentiation and downstream of ascl1a. Finally, we show that the repression of secretory cell fate by Notch signaling is mediated by the inhibition of ascl1a expression. In conclusion, this work identifies Ascl1a as a key regulator of the secretory cell lineage in the zebrafish intestine, playing the same role as Atoh1 in the mouse intestine. This highlights the diversity in the ARP/ASCL family members acting as cell fate determinants downstream from Notch signaling.

Pax4 is not essential for beta-cell differentiation in zebrafish embryos but modulates alpha-cell generation by repressing arx gene expression.

BACKGROUND: Genetic studies in mouse have demonstrated the crucial function of PAX4 in pancreatic cell differentiation. This transcription factor specifies ß- and δ-cell fate at the expense of α-cell identity by repressing Arx gene expression and ectopic expression of PAX4 in α-cells is sufficient to convert them into ß-cells. Surprisingly, no Pax4 orthologous gene can be found in chicken and Xenopus tropicalis raising the question of the function of pax4 gene in lower vertebrates such as in fish. In the present study, we have analyzed the expression and the function of the orthologous pax4 gene in zebrafish. RESULTS: pax4 gene is transiently expressed in the pancreas of zebrafish embryos and is mostly restricted to endocrine precursors as well as to some differentiating δ- and ε-cells but was not detected in differentiating ß-cells. pax4 knock-down in zebrafish embryos caused a significant increase in α-cells number while having no apparent effect on ß- and δ-cell differentiation. This rise of α-cells is due to an up-regulation of the Arx transcription factor. Conversely, knock-down of arx caused to a complete loss of α-cells and a concomitant increase of pax4 expression but had no effect on the number of ß- and δ-cells. In addition to the mutual repression between Arx and Pax4, these two transcription factors negatively regulate the transcription of their own gene. Interestingly, disruption of pax4 RNA splicing or of arx RNA splicing by morpholinos targeting exon-intron junction sites caused a blockage of the altered transcripts in cell nuclei allowing an easy characterization of the arx- and pax4-deficient cells. Such analyses demonstrated that arx knock-down in zebrafish does not lead to a switch of cell fate, as reported in mouse, but rather blocks the cells in their differentiation process towards α-cells. CONCLUSIONS: In zebrafish, pax4 is not required for the generation of the first ß- and δ-cells deriving from the dorsal pancreatic bud, unlike its crucial role in the differentiation of these cell types in mouse. On the other hand, the mutual repression between Arx and Pax4 is observed in both mouse and zebrafish. These data suggests that the main original function of Pax4 during vertebrate evolution was to modulate the number of pancreatic α-cells and its role in ß-cells differentiation appeared later in vertebrate evolution.

Fast homozygosity mapping and identification of a zebrafish ENU-induced mutation by whole-genome sequencing.

Forward genetics using zebrafish is a powerful tool for studying vertebrate development through large-scale mutagenesis. Nonetheless, the identification of the molecular lesion is still laborious and involves time-consuming genetic mapping. Here, we show that high-throughput sequencing of the whole zebrafish genome can directly locate the interval carrying the causative mutation and at the same time pinpoint the molecular lesion. The feasibility of this approach was validated by sequencing the m1045 mutant line that displays a severe hypoplasia of the exocrine pancreas. We generated 13 Gb of sequence, equivalent to an eightfold genomic coverage, from a pool of 50 mutant embryos obtained from a map-cross between the AB mutant carrier and the WIK polymorphic strain. The chromosomal region carrying the causal mutation was localized based on its unique property to display high levels of homozygosity among sequence reads as it derives exclusively from the initial AB mutated allele. We developed an algorithm identifying such a region by calculating a homozygosity score along all chromosomes. This highlighted an 8-Mb window on chromosome 5 with a score close to 1 in the m1045 mutants. The sequence analysis of all genes within this interval revealed a nonsense mutation in the snapc4 gene. Knockdown experiments confirmed the assertion that snapc4 is the gene whose mutation leads to exocrine pancreas hypoplasia. In conclusion, this study constitutes a proof-of-concept that whole-genome sequencing is a fast and effective alternative to the classical positional cloning strategies in zebrafish.