Search

Skates are cartilaginous fish whose body plan features enlarged wing-like pectoral fins, enabling them to thrive in benthic environments1,2. However, the molecular underpinnings of this unique trait remain unclear. Here we investigate the origin of this phenotypic innovation by developing the little skate Leucoraja erinacea as a genomically enabled model. Analysis of a high-quality chromosome-scale genome sequence for the little skate shows that it preserves many ancestral jawed vertebrate features compared with other sequenced genomes, including numerous ancient microchromosomes. Combining genome comparisons with extensive regulatory datasets in developing fins—including gene expression, chromatin occupancy and three-dimensional conformation—we find skate-specific genomic rearrangements that alter the three-dimensional regulatory landscape of genes that are involved in the planar cell polarity pathway. Functional inhibition of planar cell polarity signalling resulted in a reduction in anterior fin size, confirming that this pathway is a major contributor to batoid fin morphology. We also identified a fin-specific enhancer that interacts with several hoxa genes, consistent with the redeployment of hox gene expression in anterior pectoral fins, and confirmed its potential to activate transcription in the anterior fin using zebrafish reporter assays. Our findings underscore the central role of genome reorganization and regulatory variation in the evolution of phenotypes, shedding light on the molecular origin of an enigmatic trait.

7 pages, 5 figures.-- This article is licensed under a Creative Commons Attribution 4.0 International License, Metamorphosis is a widely studied post-embryonic process in which many tissues undergo dramatic modifications to adapt to the new adult lifestyle. Flatfishes represent a good example of metamorphosis in teleost fishes. During metamorphosis of flatfish, organ regression and neoformation occur, with one of the most notable changes being the migration of one of the eyes to the other side of the body. In order to create a useful and reliable tool to advance the molecular study of metamorphosis in flatfish, we generated a chromatin accessible atlas as well as gene expression profile during four developmental stages ranging from a phylotypic to a post-metamorphic stage. We identified 29,019 differentially accessible chromatin regions and 3,253 differentially expressed genes. We found stage-specific regulatory regions and gene expression profiles, supporting the quality of the results. Our work provides strongly reproducible data for further studies to elucidate the regulatory elements that ensure successful metamorphosis in flatfish species

Coordinated chromatin interactions between enhancers and promoters are critical for gene regulation. The architectural protein CTCF mediates chromatin looping and is enriched at the boundaries of topologically associating domains (TADs), which are sub-megabase chromatin structures. In vitro CTCF depletion leads to a loss of TADs but has only limited effects over gene expression, challenging the concept that CTCF-mediated chromatin structures are a fundamental requirement for gene regulation. However, how CTCF and a perturbed chromatin structure impacts gene expression during development remains poorly understood. Here we link the loss of CTCF and gene regulation during patterning and organogenesis in a ctcf knockout zebrafish model. CTCF absence leads to loss of chromatin structure and affects the expression of thousands of genes, including many developmental regulators. Our results demonstrate the essential role of CTCF in providing the structural context for enhancer-promoter interactions, thus regulating developmental genes., J.L.G.-S. received funding from the ERC (Grant Agreement No. 740041), the Spanish Ministerio de Economía y Competitividad (Grant No. BFU2016-74961-P) and the institutional grant Unidad de Excelencia María de Maeztu (MDM-2016-0687). J.T. was funded by a 2019 Leonardo Grant for Researchers and Cultural Creators, BBVA Foundation. I.I.-A. acknowledges support from the Federation of European Biochemical Societies (FEBS Long-Term Fellowship). M.F. was funded by the European Union’s Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement [#800396] and a Juan de la Cierva-Formación fellow from the Spanish Ministry of Science and Innovation (FJC2018-038233-I). JMS-P was funded by a postdoctoral fellowship from Junta de Andalucía (DOC_00512).

Evolutionary history of deuterostomes has remained unsolved and is intimately related to how chordates originated. Within the clade of deuterostomes, hemichordates and echinoderms (together called Ambulacraria) are sister groups of chordates; comparative studies among the three are insightful for understanding deuterostome evolution. Indirect-developing hemichordates, such as Ptychodera flava, develop into planktonic larvae, which are morphologically similar to echinoderm larvae, before metamorphosing into the adult body plan that retains an anteroposterior polarity homologous to that of chordates. Therefore, deciphering developmental programmes of indirect-developing hemichordates may provide a key for understanding the evolution of deuterostomes and chordate origins. Here, we characterise transcriptomes and chromatin accessibility across multiple developmental stages of P. flava and uncover biphasic developmental programmes controlled by different sets of transcription factors and their corresponding cis-regulatory elements. Transcriptome age and network analyses reveal that the gastrula transcriptome is relatively ancient and highly connected. By comparing developmental transcriptomes of hemichordate, sea urchin and amphioxus, a high conservation of gene expression during gastrulation is identified and extended to neurula stages of amphioxus, in addition to the highly similar larval transcriptomes among the three species. Moreover, P. flava possesses conserved interactions of transcription factors required for the development of echinoderm endomesoderm and chordate axial mesoderm. These results suggest a deuterostome phylotypic stage during gastrulation (corresponding to amphioxus gastrulation and neurulation) controlled by gene regulatory networks with conserved cis-regulatory interactions, and support the hypothesis that the deuterostome ancestor is an indirect-developer.

Skates are cartilaginous fish whose novel body plan features remarkably enlarged wing-like pectoral fins that allow them to thrive in benthic environments. The molecular underpinnings of this unique trait, however, remain elusive. Here we investigate the origin of this phenotypic innovation by developing the little skate Leucoraja erinacea as a genomically enabled model. Analysis of a high-quality chromosome-scale genome sequence for the little skate shows that it preserves many ancestral jawed vertebrate features compared with other sequenced genomes, including numerous ancient microchromosomes. Combining genome comparisons with extensive regulatory datasets in developing fins (gene expression, chromatin occupancy and three-dimensional (3D) conformation) we find skate-specific genomic rearrangements that alter the 3D regulatory landscape of genes involved in the planar cell polarity (PCP) pathway. Functional inhibition of PCP signaling resulted in marked reduction of anterior fin size, confirming this pathway as a major contributor of batoid fin morphology. We also identified a fin-specific enhancer that interacts with 3' HOX genes, consistent with the redeployment of Hox gene expression in anterior pectoral fins, and confirmed the potential of this element to activate transcription in the anterior fin using zebrafish reporter assays. Our findings underscore the central role of genome reorganizations and regulatory variation in the evolution of phenotypes, shedding light on the molecular origin of an enigmatic trait.

Signaling pathways control a large number of gene regulatory networks (GRNs) during animal development, acting as major tools for body plan formation [A. Pires-daSilva, R. J. Sommer, Nat. Rev. Genet. 4, 39-49 (2003)], although only a few of these pathways operate during this period [J. J. Sanz-Ezquerro, A. E. Munsterberg, € S. Stricker, Front. Cell Dev. Biol. 5, 76 (2017)]. Moreover, most of them have been largely conserved during metazoan evolution [L. S. Babonis, M. Q. Martindale, Philos. Trans. R. Soc. Lond. B Biol. Sci. 372, 20150477 (2017)]. How evolution has generated a vast diversity of animal morphologies with such a limited number of tools is still largely unknown. Here, we show that gain of interconnectivity between signaling pathways and the GRNs they control may have critically contributed to the origin of vertebrates. We perturbed the retinoic acid, Wnt, FGF, and Nodal signaling pathways during gastrulation in the invertebrate chordate amphioxus and zebrafish and compared the effects on gene expression and cis-regulatory elements (CREs). We found that multiple developmental genes gain response to these pathways through vertebrate-specific CREs. Moreover, in contrast to amphioxus, many of these CREs responded to multiple pathways in zebrafish, which reflects their high interconnectivity. Furthermore, we found that vertebrate-specific cell types are more enriched in highly interconnected genes than in tissues with more ancient origin. Thus, the increase of CREs in vertebrates integrating inputs from different signaling pathways probably contributed to gene expression complexity and to the formation of new cell types and morphological novelties in this lineage., This project has received funding from the European Research Council under the European Union’s Horizon 2020 research and innovation program (Grant agreement no. 740041) and the Spanish Ministerio de Economía y Competitividad (Grants BFU2016-74961-P and PID2019-103921GB-I00 to J.L.G.-S. and J.J.T.). This work was also supported by the institutional grant Unidad de Excelencia María de Maeztu (Grant MDM-2016-0687 to the Department of Gene regulation and morphogenesis of Centro Andaluz de Biología del Desarrollo). M.F. was funded by the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement [#800396]. M.S., H.E., and S.B. were supported by the CNRS, and H.E. and S.B. additionally by Agence Nationale de la Recherche (ANR) CHORELAND (Grant ANR-16-CE12-0008-01) and the Institut Universitaire de France. Y.-H.S. and J.-K.Y. are supported by intramural funds from Academia Sinica and grants from Ministry of Science and Technology, Taiwan (Grants 110-2326-B-001-006 to Y.-H.S. and 110-2621-B-001-001-MY3 to J.-K.Y.).

North Carolina macular dystrophy (NCMD) is a rare autosomal dominant disease affecting macular development. The disease is caused by non-coding single nucleotide variants (SNVs) in two hotspot regions near PRDM13 and by duplications in two distinct chromosomal loci, overlapping DNase I hypersensitive sites near either PRDM13 or IRX1.To unravel the mechanisms by which these variants cause disease, we first established a genome-wide multi-omics retinal database, RegRet. Integration of UMI-4C profiles we generated on adult human retina then allowed fine-mapping of the interactions of the PRDM13 and IRX1 gene promoters, and the identification of eighteen candidate cis-regulatory elements (cCREs), the activity of which was investigated by luciferase and Xenopus enhancer assays.Next, luciferase assays showed that the non-coding SNVs located in the two hotspot regions of PRDM13 affect cCRE activity, including two novel NCMD-associated non-coding SNVs that we identified. Interestingly, the cCRE containing one of these SNVs was shown to interact with the PRDM13 promoter, demonstrated in vivo activity in Xenopus, and is active at the developmental stage when progenitor cells of the central retina exit mitosis, putting forward this region as a PRDM13 enhancer.Finally, mining of single-cell transcriptional data of embryonic and adult retina revealed the highest expression of PRDM13 and IRX1 when amacrine cells start to synapse with retinal ganglion cells, supporting the hypothesis that altered PRDM13 or IRX1 expression impairs interactions between these cells during retinogenesis.Overall, this study gained insight into the cis-regulatory mechanisms of NCMD and supports that this condition is a retinal enhanceropathy.Graphical abstract

X-linked acrogigantism (X-LAG) is the most severe form of pituitary gigantism and is characterized by aggressive growth hormone (GH)-secreting pituitary tumors that occur in early childhood. X-LAG is associated with chromosome Xq26.3 duplications (the X-LAG locus typically includes VGLL1, CD40LG, ARHGEF6, RBMX, and GPR101) that lead to massive pituitary tumoral expression of GPR101, a novel regulator of GH secretion. The mechanism by which the duplications lead to marked pituitary misexpression of GPR101 alone was previously unclear. Using Hi-C and 4C-seq, we characterized the normal chromatin structure at the X-LAG locus. We showed that GPR101 is located within a topologically associating domain (TAD) delineated by a tissue-invariant border that separates it from centromeric genes and regulatory sequences. Next, using 4C-seq with GPR101, RBMX, and VGLL1 viewpoints, we showed that the duplications in multiple X-LAG-affected individuals led to ectopic interactions that crossed the invariant TAD border, indicating the existence of a similar and consistent mechanism of neo-TAD formation in X-LAG. We then identified several pituitary active cis-regulatory elements (CREs) within the neo-TAD and demonstrated in vitro that one of them significantly enhanced reporter gene expression. At the same time, we showed that the GPR101 promoter permits the incorporation of new regulatory information. Our results indicate that X-LAG is a TADopathy of the endocrine system in which Xq26.3 duplications disrupt the local chromatin architecture forming a neo-TAD. Rewiring GPR101-enhancer interaction within the new regulatory unit is likely to cause the high levels of aberrant expression of GPR101 in pituitary tumors caused by X-LAG., The work was supported by the following funding sources: Fondazione Telethon, Italy grant no. GGP20130 (to G.T.); Society for Endocrinology equipment grant (to G.T.); Intramural Research Program, Eunice Kennedy Shriver National Institute of Child Health & Human Development (NICHD), National Institutes of Health (NIH) Research project Z01-HD008920 (to C.A.S., supporting G.T., F.R.F.); Fonds d’Investissement pour la Recherche Scientifique (FIRS) of the Centre Hospitalier Universitaire de Liège (to A.F.D. and A.B.); the JABBS Foundation, UK (to A.B.); and Novo Nordisk Belgium Educational Grant, Belgium (to A.F.D. and A.B.). M.F. was funded by the European Union’s Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement (#800396) and a Juan de la Cierva-Formación fellowship from the Spanish Ministry of Science and Innovation (FJC2018-038233-I). G.T. was funded by the European Union’s Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement (#843843). A.F.D. and D.A. were supported by Action de Recherche Concertée (ARC) Grant 17/21-01 from Liège University. D.A. was supported by grants from Télévie (7461117 F, 7454719 F) and the Léon Fredericq Foundation, Belgium.

Genome-wide association studies (GWAS) have implicated the IL33 locus in asthma, but the underlying mechanisms remain unclear. Here, we identify a 5 kb region within the GWAS-defined segment that acts as an enhancer-blocking element in vivo and in vitro. Chromatin conformation capture showed that this 5 kb region loops to the IL33 promoter, potentially regulating its expression. We show that the asthma-associated single nucleotide polymorphism (SNP) rs1888909, located within the 5 kb region, is associated with IL33 gene expression in human airway epithelial cells and IL-33 protein expression in human plasma, potentially through differential binding of OCT-1 (POU2F1) to the asthma-risk allele. Our data demonstrate that asthma-associated variants at the IL33 locus mediate allele-specific regulatory activity and IL33 expression, providing a mechanism through which a regulatory SNP contributes to genetic risk of asthma., This work was supported by NIH grants R01 HL118758, R01 HL128075, R01 HL119577, R01 HL085197, U19 AI095230, UG3 OD023282 and UM1 AI114271.

15 Páginas.-- 5 Figuras.-- 3 Tablas, Histone modifications are of paramount importance during plant development. Investigating chromatin remodeling in developing oilseeds sheds light on the molecular mechanisms controlling fatty acid metabolism and facilitates the identification of new functional regions in oil crop genomes. The present study characterizes the epigenetic modifications H3K4me3 in relationship with the expression of fatty acid-related genes and transcription factors in developing sunflower seeds. Two master transcriptional regulators identified in this analysis, VIV1 (homologous to Arabidopsis ABI3) and FUS3, cooperate in the regulation of WRINKLED 1, a transcriptional factor regulating glycolysis, and fatty acid synthesis in developing oilseeds., This research was funded by the Spanish AEI/FEDER (UE), Project AGL2017-83449-R.

A key advance in our understanding of gene regulation came with the finding that the genome undergoes three-dimensional nuclear folding in a genetically determined process. This 3D conformation directly influences the association between enhancers and their target promoters. This complex interplay has been proven to be essential for gene regulation, and genetic variants affecting this process have been associated to human diseases. The development of new technologies that quantify these DNA interactions represented a revolution in the field. High throughput techniques like HiC provide a general picture of chromatin topology. However, they often lack resolution to evidence subtle effects that single nucleotide polymorphisms exert over the contacts between cis-regulatory regions and target promoters. Here we propose a cost-efficient approach to perform allele-specific chromatin conformation analysis. As a proof of concept, we analyzed the impact of a common deletion mapping between SIRPB1 promoter and one of its downstream enhancers., JLG-S acknowledges the Spanish Ministerio de Economía y Competitividad (Grant BFU2016-74961-P), the European Research Council (ERC, grant agreement No 740041) and the institutional grant Unidad de Excelencia María de Maeztu (MDM-2016-0687 to the Department of Gene regulation and morphogenesis of Centro Andaluz de Biología del Desarrollo. Project supported by a 2019 Leonardo Grant for Researchers and Cultural Creators, BBVA Foundation, granted to JJT (The BBVA Foundation accepts no responsibility for the opinions, statements and contents included in the project and/or the results thereof, which are entirely the responsibility of the authors). We wish to dedicate this work to the memory of Jose Luis Gomez-Skarmeta (1966-2020); mentor, collegue and friend.

Exposure to genotoxic stress promotes cell cycle arrest and DNA repair or apoptosis. These ¿life¿ or ¿death¿ cell fate decisions often rely on the activity of the tumor suppressor gene p53. Therefore, the precise regulation of p53 is essential to maintain tissue homeostasis and to prevent cancer development. However, how cell cycle progression has an impact on p53 cell fate decision-making is mostly unknown. In this work, we demonstrate that Drosophila p53 proapoptotic activity can be impacted by the G2/M kinase Cdk1. We find that cell cycle arrested or endocycle-induced cells are refractory to ionizing radiation-induced apoptosis. We show that p53 binding to the regulatory elements of the proapoptotic genes and its ability to activate their expression is compromised in experimentally arrested cells. Our results indicate that p53 genetically and physically interacts with Cdk1 and that p53 proapoptotic role is regulated by the cell cycle status of the cell. We propose a model in which cell cycle progression and p53 proapoptotic activity are molecularly connected to coordinate the appropriate response after DNA damage., FEDER/Ministerio de Ciencia e Innovación-Agencia Estatal de Investigación [No. PGC2018-095144-B-I00 to CE and BFU2014-54153-P to AB), the Fundación Ramón Areces, and institutional grant from Banco de Santander to the CBMSO

CTCF is an 11-zinc-finger DNA-binding protein that acts as a transcriptional repressor and insulator as well as an architectural protein required for 3D genome folding1–5. CTCF mediates long-range chromatin looping and is enriched at the boundaries of topologically associating domains, which are sub-megabase chromatin structures that are believed to facilitate enhancer-promoter interactions within regulatory landscapes 6–12. Although CTCF is essential for cycling cells and developing embryos13,14, its in vitro removal has only modest effects over gene expression5,15, challenging the concept that CTCF-mediated chromatin interactions and topologically associated domains are a fundamental requirement for gene regulation16–18. Here we link the loss of chromatin structure and gene regulation in an in vivo model and during animal development. We generated a ctcf knockout mutant in zebrafish that allows us to monitor the effect of CTCF loss of function during embryo patterning and organogenesis. CTCF absence leads to loss of chromatin structure in zebrafish embryos and affects the expression of thousands of genes, including many developmental genes. In addition, chromatin accessibility, both at CTCF binding sites and cis-regulatory elements, is severely compromised in ctcf mutants. Probing chromatin interactions from developmental genes at high resolution, we further demonstrate that promoters fail to fully establish long-range contacts with their associated regulatory landscapes, leading to altered gene expression patterns and disruption of developmental programs. Our results demonstrate that CTCF and topologically associating domains are essential to regulate gene expression during embryonic development, providing the structural basis for the establishment of developmental gene regulatory landscapes.

Genome-wide association studies (GWAS) have implicated the IL33 locus in asthma, but the underlying mechanisms remain unclear. Here, we identify a 5 kb region within the GWAS-defined segment that acts as a strong regulatory element in vivo and in vitro. Chromatin conformation capture showed that this 5 kb region loops to the IL33 promoter, potentially regulating its expression. We show that genotype at the asthma-associated SNP rs1888909, located within the 5 kb region, is associated with IL33 gene expression in human airway epithelial cells and IL-33 protein expression in human plasma, potentially through differential binding of OCT-1 (POU2F1) to the asthma-risk allele. Our data demonstrate that asthma-associated variants at the IL33 locus mediate allele-specific regulatory activity and IL33 expression, providing a novel mechanism through which a regulatory SNP contributes to genetic risk of asthma.

Genome-wide association studies (GWAS) have implicated the IL33 locus in asthma, but the underlying mechanisms remain unclear. Here, we identify a 5 kb region within the GWAS-defined segment that acts as a strong regulatory element in vivo and in vitro. Chromatin conformation capture showed that this 5 kb region loops to the IL33 promoter, potentially regulating its expression. Supporting this notion, we show that genotype at an asthma-associated SNP, rs1888909, located within the 5 kb region, is associated with IL33 gene expression in human airway epithelial cells and IL-33 protein expression in human plasma, potentially through differential binding of OCT-1 (POU2F1) to the asthma-risk allele. Our data demonstrate that asthma-associated variants at the IL33 locus mediate allele-specific regulatory activity and IL33 expression, providing a novel mechanism through which a regulatory SNP contributes to genetic risk of asthma.

The pancreas is a central organ for human diseases. Most alleles uncovered by genome-wide association studies of pancreatic dysfunction traits overlap with non-coding sequences of DNA. Many contain epigenetic marks of cis-regulatory elements active in pancreatic cells, suggesting that alterations in these sequences contribute to pancreatic diseases. Animal models greatly help to understand the role of non-coding alterations in disease. However, interspecies identification of equivalent cis-regulatory elements faces fundamental challenges, including lack of sequence conservation. Here we combine epigenetic assays with reporter assays in zebrafish and human pancreatic cells to identify interspecies functionally equivalent cis-regulatory elements, regardless of sequence conservation. Among other potential disease-relevant enhancers, we identify a zebrafish ptf1a distal-enhancer whose deletion causes pancreatic agenesis, a phenotype previously found to be induced by mutations in a distal-enhancer of PTF1A in humans, further supporting the causality of this condition in vivo. This approach helps to uncover interspecies functionally equivalent cis-regulatory elements and their potential role in human disease., This study was supported by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (ERC-2015-StG-680156-ZPR and ERC-2016-AdG-740041-EvoLand to J.L.G.-S.). J.B. is supported by an FCT CEEC grant (CEECIND/03482/2018). J.L.G.-S. is supported by the Spanish Ministerio de Economía y Competitividad (BFU2016-74961-P), the Marató TV3 Fundacion (Grant 201611) and the institutional grant Unidad de Excelencia María de Maeztu (MDM-2016-0687). R.B.C. was funded by FCT (ON2201403-CTO-BPD), IBMC (BIM/04293-UID991520-BPD) and EMBO (Short-Term Fellowship). J.Tx. (SFRH/BD/126467/2016), M.D. (SFRH/BD/135957/2018), A.E. (SFRH/BD/147762/2019), and F.J.F. (PD/BD/105745/2014) are PhD fellows from FCT. M.G. was supported by the EnvMetaGen project via the European Union’s Horizon 2020 research and innovation programme (grant 668981). This work was funded by National Funds through FCT—Fundação para a Ciência e a Tecnologia, I.P., under the project UIDB/04293/2020”.

Introductory paragraphThe pancreas is a central organ for human diseases that have a dramatic societal burden, such as pancreatic cancer and diabetes1,2. Non-coding cis-regulatory elements (CREs) of DNA control gene expression3,4, being required for proper pancreas function. Most disease-associated alleles5,6 are non-coding, often overlapping with CREs5, suggesting that alterations in these regulatory sequences contribute to human pancreatic diseases by impairing gene expression. However, functional testing of CREs in vivo is not fully explored. Here we analysed histone modifications, transcription, chromatin accessibility and interactions, to identify zebrafish pancreas CREs and their human functional equivalents, uncovering disease-associated sequences across species. We found a human pancreatic enhancer whose deletion impairs the tumour suppressor gene ARID1A expression, conferring a potential tumour suppressor role to this non-coding sequence. Additionally, we identified a zebrafish ptf1a distal enhancer which deletion generates pancreatic agenesis, demonstrating the causality of this condition in humans7 and the interspecies functional equivalency of enhancers.

Signaling pathways control a large number of gene regulatory networks (GRNs) during animal development, acting as major tools for body plan formation1. Remarkably, in contrast to the large number of transcription factors present in animal genomes, only a few of these pathways operate during development2. Moreover, most of them are largely conserved along metazoan evolution3. How evolution has generated a vast diversity of animal morphologies with such a limited number of tools is still largely unknown. Here we show that gain of interconnectivity between signaling pathways, and the GRNs they control, may have played a critical contribution to the origin of vertebrates. We perturbed the retinoic acid, Wnt, FGF and Nodal signaling pathways during gastrulation in amphioxus and zebrafish and comparatively examined its effects in gene expression and cis-regulatory elements (CREs). We found that multiple developmental genes gain response to these pathways through novel CREs in the vertebrate lineage. Moreover, in contrast to amphioxus, many of these CREs are highly interconnected and respond to multiple pathways in zebrafish. Furthermore, we found that vertebrate-specific cell types are more enriched in highly interconnected genes than those tissues with more ancestral origin. Thus, the increase of CREs in vertebrates integrating inputs from different signaling pathways probably contributed to gene expression complexity and the formation of new cell types and morphological novelties in this lineage.

The transcription factor p63 is an essential regulator of vertebrate ectoderm development, including epidermis, limbs, and craniofacial tissues. Here, we have investigated the evolutionary conservation of p63 binding sites (BSs) between zebrafish and human. First, we have analyzed sequence conservation of p63 BSs by comparing ChIP-seq data from human keratinocytes and zebrafish embryos, observing a very poor conservation. Next, we compared the gene regulatory network orchestrated by p63 in both species and found a high overlap between them, suggesting a high degree of functional conservation during evolution despite sequence divergence and the large evolutionary distance. Finally, we used transgenic reporter assays in zebrafish embryos to functionally validate a set of equivalent p63 BSs from zebrafish and human located close to genes involved in epidermal development. Reporter expression was driven by human and zebrafish BSs to many common tissues related to p63 expression domains. Therefore, we conclude that the gene regulatory network controlled by p63 is highly conserved across vertebrates despite the fact that p63-bound regulatory elements show high divergence., This work was supported by a grant from the Spanish Ministry of Economy and Competitiveness to JT (BFU2014-58449-JIN) and an ERC-Advanced grant from the European Research Council (ERC) under the European Union’s Horizon 2020 Research and Innovation Program (Grant Agreement No. 740041). JS-P was funded by a Juan de la Cierva-Incorporación fellow from the Spanish Ministry of Economy and Competitiveness (IJCI-2016-29884) and LG-F by a grant from La Marató TV3 (Grant No. 201611).

The Author(s) 2020., We investigated how the two rounds of whole-genome duplication that occurred at the base of the vertebrate lineage have impacted ancient microsyntenic associations involving developmental regulators (known as genomic regulatory blocks, GRBs). We showed that the majority of GRBs identified in the last common ancestor of chordates have been maintained as a single copy in humans. We found evidence that dismantling of the duplicated GRB copies occurred early in vertebrate evolution often through the differential retention of the regulatory gene but loss of the bystander gene’s exonic sequences. Despite the large evolutionary scale, the presence of duplicated highly conserved noncoding regions provided unambiguous proof for this scenario for multiple ancient GRBs. Remarkably, the dismantling of ancient GRB duplicates has contributed to the creation of large gene deserts associated with regulatory genes in vertebrates, providing a potentially widespread mechanism for the origin of these enigmatic genomic traits., The research has been funded by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (ERC-StG-LS2-637591 to M.I. and ERC-AdG-LS8-740041 to J.L.G.-S.), the Spanish Ministry of Science and Innovation (BFU2017-89201-P to M.I., RYC-2016-20089 and PGC2018-099392-A-I00 to I.M., BFU2016-74961-P to J.L.G.-S., and BFU2016-81887-REDT/AEI to J.L.G.-S. and M.I.), the People Programme (Marie Curie Actions) of the European Union’s Seventh Framework Programme (FP7/2007-2013) under REA grant agreement 608959, the “Centro de Excelencia Severo Ochoa 2013-2017” (SEV-2012-0208), and the “Unidad de Excelencia María de Maetzu 2017-2021” (MDM-2016-0687). We acknowledge the support of the CERCA Programme/Generalitat de Catalunya and of the Spanish Ministry of Economy, Industry and Competitiveness (MEIC) to the EMBL partnership.

© 2021 by the authors., Metamorphosis is a captivating process of change during which the morphology of the larva is completely reshaped to face the new challenges of adult life. In the case of fish, this process initiated in the brain has traditionally been considered to be a critical rearing point and despite the pioneering molecular work carried out in other flatfishes, the underlying molecular basis is still relatively poorly characterized. Turbot brain transcriptome of three developmental stages (pre-metamorphic, climax of metamorphosis and post-metamorphic) were analyzed to study the gene expression dynamics throughout the metamorphic process. A total of 1570 genes were differentially expressed in the three developmental stages and we found a specific pattern of gene expression at each stage. Unexpectedly, at the climax stage of metamorphosis, we found highly expressed genes related to the immune response, while the biological pathway enrichment analysis in pre-metamorphic and post-metamorphic were related to cell differentiation and oxygen carrier activity, respectively. In addition, our results confirm the importance of thyroid stimulating hormone, increasing its expression during metamorphosis. Based on our findings, we assume that immune system activation during the climax of metamorphosis stage could be related to processes of larval tissue inflammation, resorption and replacement, as occurs in other vertebrates., This research was funded by the MCIN/AEI/10.13039/501100011033 grant number AGL2017-89648P to JR and by “ERDF A way of making Europe”. L. Guerrero-Peña was supported by a pre-doctoral fellowship of the Spanish FPI Research Training Program funded by Spanish Economy and Competitiveness Ministry, grant number (PRE2018-085475). L. Méndez-Martínez was supported by a pre-doctoral fellowship of the Xunta de Galicia, grant number (IN606A-2020/006).

[Background] : The organisation of vertebrate genomes into topologically associating domains (TADs) is believed to facilitate the regulation of the genes located within them. A remaining question is whether TAD organisation is achieved through the interactions of the regulatory elements within them or if these interactions are favoured by the pre-existence of TADs. If the latter is true, the fusion of two independent TADs should result in the rewiring of the transcriptional landscape and the generation of ectopic contacts., [Results]: We show that interactions within the PAX3 and FOXO1 domains are restricted to their respective TADs in normal conditions, while in a patient-derived alveolar rhabdomyosarcoma cell line, harbouring the diagnostic t(2;13)(q35;q14) translocation that brings together the PAX3 and FOXO1 genes, the PAX3 promoter interacts ectopically with FOXO1 sequences. Using a combination of 4C-seq datasets, we have modelled the three-dimensional organisation of the fused landscape in alveolar rhabdomyosarcoma., [Conclusions]: The chromosomal translocation that leads to alveolar rhabdomyosarcoma development generates a novel TAD that is likely to favour ectopic PAX3:FOXO1 oncogene activation in non-PAX3 territories. Rhabdomyosarcomas may therefore arise from cells which do not normally express PAX3. The borders of this novel TAD correspond to the original 5'- and 3'- borders of the PAX3 and FOXO1 TADs, respectively, suggesting that TAD organisation precedes the formation of regulatory long-range interactions. Our results demonstrate that, upon translocation, novel regulatory landscapes are formed allowing new intra-TAD interactions between the original loci involved., This work was supported by The Institute of Cancer Research and Cancer Research UK (Grant C1178/A4520). JJC was funded by grants from Spanish Ministerio de Ciencia e Innovación (BFU2011-22928) and the European Commission (PCIG10-GA-2011-303904). JLGS was funded by grants from Spanish Ministerio de Economía y Competitividad (BFU2013-41322-P) and the Andalusian Government (BIO-396). DPD and IIA were funded by grants from the Spanish Ministerio de Economia y Competitividad () and from the Junta de Andalucia (C2A). BV-B held a studentship from The Institute of Cancer Research., We acknowledge support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).

We investigated how the two rounds of whole genome duplication that occurred at the base of the vertebrate lineage have impacted ancient microsyntenic associations involving developmental regulators (known as genomic regulatory blocks, GRBs). We showed that the majority of GRBs present in the last common ancestor of chordates have been maintained as a single copy in humans. We found evidence that dismantling of the additional GRB copies occurred early in vertebrate evolution often through the differential retention of the regulatory gene but loss of the bystander gene’s exonic sequences. Despite the large evolutionary scale, the presence of duplicated highly conserved non-coding regions provided unambiguous proof for this scenario for dozens of ancient GRBs. Remarkably, the dismantling of ancient GRB duplicates has contributed to the creation of large gene deserts associated with regulatory genes in vertebrates, providing a widespread mechanism for the origin of these enigmatic genomic traits.

The transcription factor p63 is a master regulator of ectoderm development. Although previous studies show that p63 triggers epidermal differentiation in vitro, the roles of p63 in developing embryos remain poorly understood. Here, we use zebrafish embryos to analyze in vivo how p63 regulates gene expression during development. We generate tp63-knock-out mutants that recapitulate human phenotypes and show down-regulated epidermal gene expression. Following p63-binding dynamics, we find two distinct functions clearly separated in space and time. During early development, p63 binds enhancers associated to neural genes, limiting Sox3 binding and reducing neural gene expression. Indeed, we show that p63 and Sox3 are co-expressed in the neural plate border. On the other hand, p63 acts as a pioneer factor by binding non-accessible chromatin at epidermal enhancers, promoting their opening and epidermal gene expression in later developmental stages. Therefore, our results suggest that p63 regulates cell fate decisions during vertebrate ectoderm specification., The transcription factor p63 is a master regulator of ectoderm development. Here the authors show that zebrafish p63 binds enhancers associated with neural genes to limit Sox3 binding and gene expression while also acting as a pioneer factor by promoting chromatin opening at epidermal gene enhancers.

Vertebrates have greatly elaborated the basic chordate body plan and evolved highly distinctive genomes that have been sculpted by two whole-genome duplications. Here we sequence the genome of the Mediterranean amphioxus (Branchiostoma lanceolatum) and characterize DNA methylation, chromatin accessibility, histone modifications and transcriptomes across multiple developmental stages and adult tissues to investigate the evolution of the regulation of the chordate genome. Comparisons with vertebrates identify an intermediate stage in the evolution of differentially methylated enhancers, and a high conservation of gene expression and its cis-regulatory logic between amphioxus and vertebrates that occurs maximally at an earlier mid-embryonic phylotypic period. We analyse regulatory evolution after whole-genome duplications, and find that—in vertebrates—over 80% of broadly expressed gene families with multiple paralogues derived from whole-genome duplications have members that restricted their ancestral expression, and underwent specialization rather than subfunctionalization. Counter-intuitively, paralogues that restricted their expression increased the complexity of their regulatory landscapes. These data pave the way for a better understanding of the regulatory principles that underlie key vertebrate innovations., This research was funded primarily by the European Research Council (ERC) under the European Union’s Horizon 2020 and Seventh Framework Program FP7 research and innovation programs (ERC-AdG-LS8-740041 to J.L.G.-S., ERC-StG-LS2-637591 to M.I., a Marie Sklodowska-Curie Grant (658521) to I.M. and a FP7/2007-2013-ERC-268513 to P.W.H.H.), the Spanish Ministerio de Economía y Competitividad (BFU2016-74961-P to J.L.G.-S., RYC-2016-20089 to I.M., BFU2014-55076-P and BFU2017-89201-P to M.I. and BFU2014-55738-REDT to J.L.G.-S, M.I. and J.R.M.-M), the ‘Centro de Excelencia Severo Ochoa 2013-2017’(SEV-2012-0208), the ‘Unidad de Excelencia María de Maetzu 2017-2021’(MDM-2016-0687), the People Program (Marie Curie Actions) of the European Union’s Seventh Framework Program FP7 under REA grant agreement number 607142 (DevCom) to J.L.G.-S., and the CNRS and the ANR (ANR16-CE12-0008-01) to H.E. O.B. was supported by an Australian Research Council Discovery Early Career Researcher Award (DECRA; DE140101962).

Cis-regulatory changes are arguably the primary evolutionary source of animal morphological diversity. With the recent explosion of genome-wide comparisons of the cis-regulatory content in different animal species is now possible to infer general principles underlying enhancer evolution. However, these studies have also revealed numerous discrepancies and paradoxes, suggesting that the mechanistic causes and modes of cis-regulatory evolution are still not well understood and are probably much more complex than generally appreciated. Here, we argue that the mutational mechanisms and genomic regions generating new regulatory activities must comply with the constraints imposed by the molecular properties of cis-regulatory elements (CREs) and the organizational features of long-range chromatin interactions. Accordingly, we propose a new integrative evolutionary framework for cis-regulatory evolution based on two major premises for the origin of novel enhancer activity: (i) an accessible chromatin environment and (ii) compatibility with the 3D structure and interactions of pre-existing CREs. Mechanisms and DNA sequences not fulfilling these premises, will be less likely to have a measurable impact on gene expression and as such, will have a minor contribution to the evolution of gene regulation. Finally, we discuss current comparative cis-regulatory data under the light of this new evolutionary model, and propose that the two most prominent mechanisms for the evolution of cis-regulatory changes are the overprinting of ancestral CREs and the exaptation of transposable elements., J.J.T. is funded by a contract associated to Spanish government project (BFU2013-41322-P granted to J.L. Gómez-Skarmeta) and I.M. holds a Juan de la Cierva postdoctoral contract.

Evolution of organismal complexity and species diversity depends on the emergence of novel gene functions. Nevertheless, evolution rarely produces novelties from scratch but works on the weak promiscuous preexisting activities or appears by genomic tinkering. We provide evidence of how rearrangement of conserved regulatory blocks can act as a driving force for gene cooption and evolution of novel developmental mechanisms at the base of important ecological adaptations. We gain insight into a crucial system for segregation of neuronal progenitors within the hindbrain: the evolutionary origin of the actomyosin-dependent cell-sorting mechanism, with rac3b as a main effector. We unveil that the rac3b/rfng/sgca regulatory cluster—specifically expressed at boundaries—emerged by establishment of novel long-range cis-regulatory interactions, allowing the evolution of a backup regulatory mechanism for cell segregation.

© 2018 Irastorza-Azcarate et al., The use of 3C-based methods has revealed the importance of the 3D organization of the chromatin for key aspects of genome biology. However, the different caveats of the variants of 3C techniques have limited their scope and the range of scientific fields that could benefit from these approaches. To address these limitations, we present 4Cin, a method to generate 3D models and derive virtual Hi-C (vHi-C) heat maps of genomic loci based on 4C-seq or any kind of 4C-seq-like data, such as those derived from NG Capture-C. 3D genome organization is determined by integrative consideration of the spatial distances derived from as few as four 4C-seq experiments. The 3D models obtained from 4C-seq data, together with their associated vHi-C maps, allow the inference of all chromosomal contacts within a given genomic region, facilitating the identification of Topological Associating Domains (TAD) boundaries. Thus, 4Cin offers a much cheaper, accessible and versatile alternative to other available techniques while providing a comprehensive 3D topological profiling. By studying TAD modifications in genomic structural variants associated to disease phenotypes and performing cross-species evolutionary comparisons of 3D chromatin structures in a quantitative manner, we demonstrate the broad potential and novel range of applications of our method., This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No 740041 to JLGS), the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement #658521 to IM (https://ec.europa.eu/programmes/horizon2020/en/h2020-section/marie-sklodowska-curie-actions), the Spanish Ministerio de Economía y Competitividad (http://www.mineco.gob.es/) (BFU2013-41322-P and BFU2016-74961-P to JLGS; BFU2014-58449-JIN to JJT; BFU2013-40866-P to DPD and IIA) and the Andalusian government (http://www.juntadeandalucia.es/index.html/) (BIO-396 to JLGS; C2A (EE: 2013/2506) to DPD and IIA). RDA holds a FPU fellowship from the Spanish Ministerio de Educación, Cultura y Deporte (MECD) (http://www.mecd.gob.es/portada-mecd/).

There is no obvious morphological counterpart of the autopod (wrist/ankle and digits) in living fishes. Comparative molecular data may provide insight into understanding both the homology of elements and the evolutionary developmental mechanisms behind the fin to limb transition. In mouse limbs the autopod is built by a >late> phase of Hoxd and Hoxa gene expression, orchestrated by a set of enhancers located at the 5? end of each cluster. Despite a detailed mechanistic understanding of mouse limb development, interpretation of Hox expression patterns and their regulation in fish has spawned multiple hypotheses as to the origin and function of >autopod> enhancers throughout evolution. Using phylogenetic footprinting, epigenetic profiling, and transgenic reporters, we have identified and functionally characterized hoxD and hoxA enhancers in the genomes of zebrafish and the spotted gar, Lepisosteus oculatus, a fish lacking the whole genome duplication of teleosts. Gar and zebrafish > autopod> enhancers drive expression in the distal portion of developing zebrafish pectoral fins, and respond to the same functional cues as their murine orthologs. Moreover, gar enhancers drive reporter gene expression in both the wrist and digits of mouse embryos in patterns that are nearly indistinguishable from their murine counterparts. These functional genomic data support the hypothesis that the distal radials of bony fish are homologous to the wrist and/or digits of tetrapods., This study was supported by The Brinson Foundation (N.H.S.), National Science Foundation Doctoral Dissertation Improvement Grant 1311436 (to A.R.G.), Brazilian National Council for Scientific and Technological Development Grants 402754/2012-3 and 477658/2012-1 (to I.S.), National Institutes of Health Grant R01RR020833 (R01OD011116) (to J.H.P), “Initiative Evolutionary Biology” from the Volkswagen Foundation, Germany and Alexander von Humboldt-Foundation (I.B.), and Spanish and Andalusian Governments Grants BFU2010-14839, BFU2013-41322-P, and Proyecto de Excelencia BIO-396 (to J.L.G.-S.).

Atonal homolog 1 (Atoh1) and Neurogenin1 (Neurog1) are basic Helix-Loop-Helix (bHLH) transcription factors crucial for the generation of hair cells (HCs) and neurons in the inner ear. Both genes are induced early in development, but the expression of Atoh1 is counteracted by Neurog1. As a result, HC development is prevented during neurogenesis. This work aimed at understanding the molecular basis of this interaction. Atoh1 regulation depends on a 3’Atoh1-enhancer that is the site for Atoh1 autoregulation. Reporter assays on chick embryos and P19 cells show that Neurog1 hampers the autoactivation of Atoh1, the effect being cell autonomous and independent on Notch activity. Assay for Transposase-Accessible Chromatin with high throughput sequencing (ATAC-Seq) analysis shows that the region B of the 3’Atoh1-enhancer is accessible during development and sufficient for both activation and repression. Neurog1 requires the regions flanking the class A E-box to show its repressor effect, however, it does not require binding to DNA for Atoh1 repression. This depends on the dimerization domains Helix-1 and Helix-2 and the reduction of Atoh1 protein levels. The results point towards the acceleration of Atoh1 mRNA degradation as the potential mechanism for the reduction of Atoh1 levels. Such a mechanism dissociates the prevention of Atoh1 expression in neurosensory progenitors from the unfolding of the neurogenic program., This work was supported by Fundació La Marató-TV3, the Spanish Ministerio de Economía y Competitividad MINECO/FEDER grant BFU2015-67400-P, and Unidad de Excelencia María de Maetzu 2015-19 grant MDM-2014-0370.

et al., Cardiac progenitors are specified early in development and progressively differentiate and mature into fully functional cardiomyocytes. This process is controlled by an extensively studied transcriptional program. However, the regulatory events coordinating the progression of such program from development to maturation are largely unknown. Here, we show that the genome organizer CTCF is essential for cardiogenesis and that it mediates genomic interactions to coordinate cardiomyocyte differentiation and maturation in the developing heart. Inactivation of Ctcf in cardiac progenitor cells and their derivatives in vivo during development caused severe cardiac defects and death at embryonic day 12.5. Genome wide expression analysis in Ctcf mutant hearts revealed that genes controlling mitochondrial function and protein production, required for cardiomyocyte maturation, were upregulated. However, mitochondria from mutant cardiomyocytes do not mature properly. In contrast, multiple development regulatory genes near predicted heart enhancers, including genes in the IrxA cluster, were downregulated in Ctcf mutants, suggesting that CTCF promotes cardiomyocyte differentiation by facilitating enhancer-promoter interactions. Accordingly, loss of CTCF disrupts gene expression and chromatin interactions as shown by chromatin conformation capture followed by deep sequencing. Furthermore, CRISPR-mediated deletion of an intergenic CTCF site within the IrxA cluster alters gene expression in the developing heart. Thus, CTCF mediates local regulatory interactions to coordinate transcriptional programs controlling transitions in morphology and function during heart development., This work was funded by the Spanish Ministerio de Economia y Competitividad (grants BFU2011-23083 and BFU2014-54608-P to MM, BFU2013-41322-P to JLGS and BFU2014-58449-JIN to JJT), the Comunidad Autónoma de Madrid (grant CELLDD-CM to MM), the Andalusian regional government (grant BIO-396 to JLGS), the Universidad Pablo de Olavide (postdoctoral grant to JJT). PDO is funded by the Heart and Stroke Foundation of Canada (G-17-0018613), Operational Funds from the Hospital for Sick Children, the Natural Sciences and Engineering Research Council of Canada (NSERC) (500865 to PDO), and the Canadian Institutes of Health Research (CIHR) (PJT-149046). The CNIC is supported by the Spanish Ministerio de Economia y Competitividad and the Pro-CNIC Foundation, and is a Severo Ochoa Center of Excellence (SEV-2015-0505).

Genome-wide profiling of DNA methylation and histone modifications answered many questions as to how the genes are regulated on a global scale and what their epigenetic makeup is. Yet, little is known about the function of these marks during early vertebrate embryogenesis. Here we provide detailed protocols for ChIP-seq and MethylCap-seq procedures applied to zebrafish (Danio rerio) embryonic material at four developmental stages. As a proof of principle, we have profiled on a global scale a number of post-translational histone modifications including H3K4me1, H3K4me3 and H3K27ac. We demonstrate that these marks are dynamic during early development and that such developmental transitions can be detected by ChIP-seq. In addition, we applied MethylCap-seq to show that developmentally-regulated DNA methylation remodeling can be detected by such a procedure. Our MethylCap-seq data concur with previous DNA methylation studies of early zebrafish development rendering this method highly suitable for the global assessment of DNA methylation in early vertebrate embryos.

Bogdanovic, Ozren et al., The vertebrate body plan and organs are shaped during a conserved embryonic phase called the phylotypic stage. However, the mechanisms that guide the epigenome through this transition and their evolutionary conservation remain elusive. Here we report widespread DNA demethylation of enhancers during the phylotypic period in zebrafish, Xenopus tropicalis and mouse. These enhancers are linked to developmental genes that display coordinated transcriptional and epigenomic changes in the diverse vertebrates during embryogenesis. Binding of Tet proteins to (hydroxy)methylated DNA and enrichment of 5-hydroxymethylcytosine in these regions implicated active DNA demethylation in this process. Furthermore, loss of function of Tet1, Tet2 and Tet3 in zebrafish reduced chromatin accessibility and increased methylation levels specifically at these enhancers, indicative of DNA methylation being an upstream regulator of phylotypic enhancer function. Overall, our study highlights a regulatory module associated with the most conserved phase of vertebrate embryogenesis and suggests an ancient developmental role for Tet dioxygenases., Spanish and Andalusian government grants BFU2013-41322-P and BIO-396 to J.L.G.-S. supported this work. R.L. was supported by an Australian Research Council Future Fellowship (FT120100862) and a Sylvia and Charles Viertel Senior Medical Research Fellowship, and work in the laboratory of R.L. was funded by the Australian Research Council, National Health and Medical Research Council, and the Raine Medical Research Foundation. O.B. is supported by an Australian Research Council Discovery Early Career Researcher Award (DECRA; DE140101962). The laboratory of M.V. is supported by grants from the Netherlands Organisation for Scientific Research (NWO-VIDI; 864.09.003) and Cancer Genomics Netherlands, a European Research Council starting grant (309384) and the European Union Framework Programme 7 Network of Excellence EpiGeneSys. J.R.E. was supported by the Gordon and Betty Moore Foundation (GBMF3034) and is an Investigator of the Howard Hughes Medical Institute. Work in the laboratory of M.M. is funded by grants from the Ministerio de Economia y Competitividad (BFU2011-23083), Comunidad Autónoma de Madrid (CELLDD-CM), and by the Pro-CNIC Foundation. This work has been supported by a grant from the US National Institutes of Health (National Institute of Child Health and Human Development, grant R01HD069344) to G.J.C.V.

Acemel, Rafael D. et al., The HoxA and HoxD gene clusters of jawed vertebrates are organized into bipartite three-dimensional chromatin structures that separate long-range regulatory inputs coming from the anterior and posterior Hox-neighboring regions. This architecture is instrumental in allowing vertebrate Hox genes to pattern disparate parts of the body, including limbs. Almost nothing is known about how these three-dimensional topologies originated. Here we perform extensive 4C-seq profiling of the Hox cluster in embryos of amphioxus, an invertebrate chordate. We find that, in contrast to the architecture in vertebrates, the amphioxus Hox cluster is organized into a single chromatin interaction domain that includes long-range contacts mostly from the anterior side, bringing distant cis-regulatory elements into contact with Hox genes. We infer that the vertebrate Hox bipartite regulatory system is an evolutionary novelty generated by combining ancient long-range regulatory contacts from DNA in the anterior Hox neighborhood with new regulatory inputs from the posterior side., Work was funded by grants from the Ministerio de Economía y Competitividad (BFU2013-41322-P to J.L.G.-S.; Juan de la Cierva postdoctoral contract to I.M.; BFU2014-58449-JIN to J.J.T.); the Andalusian government (BIO-396 to J.L.G.-S.; C2A (EE: 2013/2506) to D.P.D. and I.I.-A.); the European Research Council (ERC; grant 268513) to P.W.H.H. and F.M.; a European Molecular Biology Organization (EMBO) short fellowship to I.M.; the Universidad Pablo de Olavide to J.J.T.; and Conicyt 'Becas Chile' to D.A.

The unicellular ancestor of animals had a complex repertoire of genes linked to multicellular processes. This suggests that changes in the regulatory genome, rather than in gene innovation, were key to the origin of animals. Here, we carry out multiple functional genomic assays in Capsaspora owczarzaki, the unicellular relative of animals with the largest known gene repertoire for transcriptional regulation. We show that changing chromatin states, differential lincRNA expression, and dynamic cis-regulatory sites are associated with life cycle transitions in Capsaspora. Moreover, we demonstrate conservation of animal developmental transcription-factor networks and extensive network interconnection in this premetazoan organism. In contrast, however, Capsaspora lacks animal promoter types, and its regulatory sites are small, proximal, and lack signatures of animal enhancers. Overall, our results indicate that the emergence of animal multicellularity was linked to a major shift in genome cis-regulatory complexity, most notably the appearance of distal enhancer regulation., This work was supported by an Institució Catalana de Recerca i Estudis Avançats contract, a European Research Council Consolidator Grant (ERC-2012-Co-616960), and a grant from Ministerio de Economía y Competitividad (MINECO) (BFU-2011-23434) (to I.R.-T.). We also acknowledge financial support from Secretaria d’Universitats i Recerca del Departament d’Economia i Coneixement de la Generalitat de Catalunya (project 2014 SGR 619). The work in L.D.’s laboratory was supported by grants from the Spanish “Ministerio de Educación y Ciencia” (SAF2013-48926-P), AGAUR, and the European Commission’s 7th Framework Program 4DCellFate (277899). A.S.-P. is supported by an EMBO Long-Term Fellowship (ALTF 841-2014). J.L.G.-S. was funded by grants from Ministerio de Economía y Competitividad (BFU2013-41322-P) and the Andalusian Government (BIO-396). J.J.T. has a postdoctoral grant from the University Pablo de Olavide. The CRG/UPF Proteomics Unit is part of the “Plataforma de Recursos Biomoleculares y Bioinformáticos (ProteoRed)” supported by a grant from Instituto de Salud Carlos III (ISCIII) (PT13/0001).

Accurate transcriptional control of genes is fundamental for the correct functioning of organs and developmental processes. This control depends on the interplay between the promoter of genes and other noncoding sequences, whose interaction is mediated by 3D chromatin arrangements. Thus, the detailed description of transcriptional regulatory landscapes is essential to understand the mechanisms of transcriptional regulation. However, to achieve that, two important challenges have to be faced: (1) the identification of the noncoding sequences that contribute to gene transcription and (2) the association of these sequences to the respective genes they control. In this chapter, we describe two protocols that allow overcoming these important challenges: the assay for transposase-accessible chromatin using sequencing (ATAC-seq) and circularized chromosome conformation capture (4C-seq). ATAC-seq is a very efficient technique that, using a very low number of cells as starting material, allows the identification of active chromatin regions genome wide, whereas 4C-seq detects the subset of sequences that interact specifically with the promoter of a given gene. When combined, both techniques provide a comprehensive snapshot of the regulatory landscapes of developmental genes. The protocols we present here have been optimized for teleost fish samples, zebrafish and medaka, allowing the in-depth study of transcriptional regulation in these two emerging animal models. Given the amenability and easy genetic manipulation of these two experimental systems, we anticipate that they will be important in revealing general principles of the vertebrate regulatory genome., Work was funded by grants from Ministerio de Economía y Competitividad (BFU2013-41322-P and BFU2014-58449-JIN) and the Andalusian Government (BIO-396). JB was supported by an FCT Investigator grant (IF/00654/2013).