[en] 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.
Rhinn, M. & Dollé, P. Retinoic acid signalling during development. Development 139, 843–858 (2012).
Niederreither, K. & Dollé, P. Retinoic acid in development: Towards an integrated view. Nat. Rev. Genet. 9, 541–553 (2008).
Ghyselinck, N. B. & Duester, G. Retinoic acid signaling pathways. Development 146, dev167502 (2019).
Shimozono, S., Iimura, T., Kitaguchi, T., Higashijima, S. & Miyawaki, A. Visualization of an endogenous retinoic acid gradient across embryonic development. Nature 496, 363–366 (2013).
Larsen, H. L. & Grapin-Botton, A. The molecular and morphogenetic basis of pancreas organogenesis. Semin. Cell Dev. Biol. 66, 51–68 (2017).
Zorn, A. M. & Wells, J. M. Vertebrate endoderm development and organ formation. Annu. Rev. Cell Dev. Biol. 25, 221 (2009).
Stafford, D. & Prince, V. E. Retinoic acid signaling is required for a critical early step in zebrafish pancreatic development. Curr. Biol. 12, 1215–1220 (2002).
Stafford, D., Hornbruch, A., Mueller, P. R. & Prince, V. E. A conserved role for retinoid signaling in vertebrate pancreas development. Dev. Genes Evol. 214, 432–441 (2004).
Bayha, E., Jørgensen, M. C., Serup, P. & Grapin-Botton, A. Retinoic acid signaling organizes endodermal organ specification along the entire antero-posterior axis. PLoS One 4, e5845 (2009).
Chen, Y. et al. Retinoic acid signaling is essential for pancreas development and promotes endocrine at the expense of exocrine cell differentiation in Xenopus. Dev. Biol. 271, 144–160 (2004).
Molotkov, A., Molotkova, N. & Duester, G. Retinoic acid generated byRaldh2 in mesoderm is required for mouse dorsal endodermal pancreas development. Dev. Dyn. 232, 950–957 (2005).
Martín, M. et al. Dorsal pancreas agenesis in retinoic acid-deficient Raldh2 mutant mice. Dev. Biol. 284, 399–411 (2005).
Pagliuca, F. W. et al. Generation of functional human pancreatic β cells in vitro. Cell 159, 428–439 (2014).
Stafford, D. et al. Retinoids signal directly to zebrafish endoderm to specify insulin-expressing beta-cells. Development 133, 949–956 (2006).
Song, J., Kim, H. J., Gong, Z., Liu, N.-A. & Lin, S. Vhnf1 acts downstream of Bmp, Fgf, and RA signals to regulate endocrine beta cell development in zebrafish. Dev. Biol. 303, 561–575 (2007).
Hernandez, R. E. et al. vhnf1 integrates global RA patterning and local FGF signals to direct posterior hindbrain development in zebrafish. Development 131, 4511–4520 (2004).
Nolte, C., De Kumar, B. & Krumlauf, R. Hox genes: Downstream “effectors” of retinoic acid signaling in vertebrate embryogenesis. Genesis 57, 1–17 (2019).
Gere-Becker, M. B., Pommerenke, C., Lingner, T. & Pieler, T. Retinoic acid-induced expression of Hnf1b and Fzd4 is required for pancreas development in Xenopus laevis. Development 145, dev161372 (2018).
Chatagnon, A. et al. RAR/RXR binding dynamics distinguish pluripotency from differentiation associated cis-regulatory elements. Nucleic Acids Res. 43, 4833–4854 (2015).
Paschaki, M. et al. Transcriptomic analysis of murine embryos lacking endogenous retinoic acid signaling. PLoS One 8, e62274 (2013).
Arima, K. et al. Global analysis of RAR-responsive genes in the Xenopus neurula using cDNA microarrays. Dev. Dyn. 232, 414–431 (2005).
Samarut, E. et al. Retinoic acid receptor subtype-specific transcriptotypes in the early zebrafish embryo. Mol. Endocrinol. 28, 260–272 (2014).
Feng, L., Hernandez, R. E., Waxman, J. S., Yelon, D. & Moens, C. B. Dhrs3a regulates retinoic acid biosynthesis through a feedback inhibition mechanism. Dev. Biol. 338, 1–14 (2010).
Cunningham, T. J., Colas, A. & Duester, G. Early molecular events during retinoic acid induced differentiation of neuromesodermal progenitors. Biol. Open 5, 1821–1833 (2016).
Berenguer, M., Meyer, K. F., Yin, J. & Duester, G. Discovery of genes required for body axis and limb formation by global identification of retinoic acid-regulated epigenetic marks. PLoS Biol. 18, e3000719 (2020).
Zhang, T., Guo, X. & Chen, Y. Retinoic acid-activated Ndrg1a represses Wnt/β-catenin signaling to allow Xenopus pancreas, oesophagus, stomach, and duodenum specification. PLoS One 8, e65058 (2013).
White, R. J., Nie, Q., Lander, A. D. & Schilling, T. F. Complex regulation of cyp26a1 creates a robust retinoic acid gradient in the zebrafish embryo. PLoS Biol. 5, e304 (2007).
Schilling, T. F., Nie, Q. & Lander, A. D. Dynamics and precision in retinoic acid morphogen gradients. Curr. Opin. Genet. Dev. 22, 562–569 (2012).
Kinkel, M. D. et al. Cyp26 enzymes function in endoderm to regulate pancreatic field size. Proc. Natl. Acad. Sci. U.S.A. 106, 7864–7869 (2009).
Delacroix, L. et al. Cell-specific interaction of retinoic acid receptors with target genes in mouse embryonic fibroblasts and embryonic stem cells. Mol. Cell. Biol. 30, 231–244 (2010).
Mahony, S. et al. Ligand-dependent dynamics of retinoic acid receptor binding during early neurogenesis. Genome Biol. 12, R2 (2011).
Moutier, E. et al. Retinoic acid receptors recognize the mouse genome through binding elements with diverse spacing and topology. J. Biol. Chem. 287, 26328 (2012).
Buenrostro, J. D., Giresi, P. G., Zaba, L. C., Chang, H. Y. & Greenleaf, W. J. Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position. Nat. Methods 10, 1213–1218 (2013).
Hiller, M. et al. Computational methods to detect conserved non-genic elements in phylogenetically isolated genomes: Application to zebrafish. Nucleic Acids Res. 41, e151 (2013).
Quillien, A. et al. Robust identification of developmentally active endothelial enhancers in zebrafish using FANS-assisted ATAC-Seq. Cell Rep. 20, 709–720 (2017).
Paik, E. J. et al. A Cdx4-Sall4 regulatory module controls the transition from mesoderm formation to embryonic hematopoiesis. Stem Cell Rep. 1, 425–436 (2013).
Love, C. E. & Prince, V. E. Expression and retinoic acid regulation of the zebrafish nr2f orphan nuclear receptor genes. Dev. Dyn. 241, 1603–1615 (2012).
Germain, P. et al. Differential action on coregulator interaction defines inverse retinoid agonists and neutral antagonists. Chem. Biol. 16, 479–489 (2009).
Nadendla, E. et al. An unexpected mode of binding defines BMS948 as a full retinoic acid receptor β (RARβ, NR1B2) selective agonist. PLoS One 10, e0123195 (2015).
Yang, H. et al. A map of cis-regulatory elements and 3D genome structures in zebrafish. Nature 588, 337–343 (2020).
Zaret, K. S. Pioneer transcription factors initiating gene network changes. Annu. Rev. Genet. 54, 367–385 (2020).
Sun, Z. & Hopkins, N. vhnf1, the MODY5 and familial GCKD-associated gene, regulates regional specification of the zebrafish gut, pronephros, and hindbrain. Genes Dev. 15, 3217–3229 (2001).
Haumaitre, C. et al. Lack of TCF2/vHNF1 in mice leads to pancreas agenesis. Proc. Natl. Acad. Sci. U.S.A. 102, 1490–1495 (2005).
Carrasco, M., Delgado, I., Soria, B., Martín, F. & Rojas, A. GATA4 and GATA6 control mouse pancreas organogenesis. J. Clin. Investig. 122, 3504–3515 (2012).
Allen, H. L. et al. GATA6 haploinsufficiency causes pancreatic agenesis in humans. Nat. Genet. 44, 20–22 (2011).
Bossard, P. & Zaret, K. S. GATA transcription factors as potentiators of gut endoderm differentiation. Development 125, 4909–4917 (1998).
Mizoguchi, T. et al. Sdf1/Cxcr4 signaling controls the dorsal migration of endodermal cells during zebrafish gastrulation. Development 135, 2521–2529 (2008).
Lavergne, A. et al. Pancreatic and intestinal endocrine cells in zebrafish share common transcriptomic signatures and regulatory programmes. BMC Biol. 18, 109 (2020).
Dobin, A. et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21 (2013).
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).
Morley, R. H. et al. A gene regulatory network directed by zebrafish No tail accounts for its roles in mesoderm formation. Proc. Natl. Acad. Sci. U.S.A. 106, 3829–3834 (2009).
Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359 (2012).
Zhang, Y. et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 9, R137 (2008).
Heinz, S. et al. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol. Cell 38, 576–589 (2010).
Yu, G., Wang, L.-G. & He, Q.-Y. ChIPseeker: an R/Bioconductor package for ChIP peak annotation, comparison and visualization. Bioinformatics 31, 2382–2383 (2015).
Buenrostro, J. D., Wu, B., Chang, H. Y. & Greenleaf, W. J. ATAC-seq: A method for assaying chromatin accessibility genome-wide. Curr. Protoc. Mol. Biol. 109, 21.29.1-21.29.9 (2015).
Ou, J. et al. ATACseqQC: A bioconductor package for post-alignment quality assessment of ATAC-seq data. BMC Genomics 19, 169 (2018).
Zhu, L. J. et al. ChIPpeakAnno: A Bioconductor package to annotate ChIP-seq and ChIP-chip data. BMC Bioinform. 11, 237 (2010).
Ross-Innes, C. S. et al. Differential oestrogen receptor binding is associated with clinical outcome in breast cancer. Nature 481, 389–393 (2012).
Fisher, S. et al. Evaluating the biological relevance of putative enhancers using Tol2 transposon-mediated transgenesis in zebrafish. Nat. Protoc. 1, 1297–1305 (2006).