Abstract Zebrafish, a popular model for embryonic development and for modelling human diseases, has so far lacked a systematic functional annotation programme akin to those in other animal models. To address this, we formed the international DANIO-CODE consortium and created the first 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 reanalysed published genomics data, which we used to improve existing annotations and show its utility in experimental design. We identified over 140,000 cis-regulatory elements in development, including novel classes with distinct features dependent on their activity in time and space. We delineated the distinction between regulatory elements active during zygotic genome activation and those active during organogenesis, identifying new aspects of how they relate to each other. Finally, we matched regulatory elements and epigenomic landscapes between zebrafish and mouse and predict functional relationships between them beyond sequence similarity, extending the utility of zebrafish developmental genomics to mammals.

ABSTRACT Mammalian sex is determined by opposing networks of ovarian and testicular genes that are well characterized. However, its epigenetic regulation is still largely unknown, thus limiting our understanding of a fundamental process for species propagation. Here we explore the 3D chromatin landscape of sex determination in vivo , by profiling FACS-sorted embryonic mouse gonadal populations, prior and after sex determination, in both sexes. We integrate Hi-C with ChIP-seq experiments using METALoci , a novel genome spatial auto-correlation analysis that identifies 3D enhancer hubs across the genome. We uncover a prominent rewiring of chromatin interactions during sex determination, affecting the enhancer hubs of hundreds of genes that display temporal- and sex-specific expression. Moreover, the identification of the 3D enhancer hubs allows the reconstruction of regulatory networks, revealing key transcription factors involved in sex determination. By combining predictive approaches and validations in transgenic mice we identify a novel Fgf9 regulatory hub, deletion of which results in male-to-female sex reversal with the upregulation of ovarian-specific markers and the initiation of meiosis. Thus, spatial auto-correlation analysis is an effective strategy to identify regulatory networks associated to biological processes and to further characterize the functional role of the 3D genome.

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 folding 1–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 embryos 13,14 , its in vitro removal has only modest effects over gene expression 5,15 , challenging the concept that CTCF-mediated chromatin interactions and topologically associated domains are a fundamental requirement for gene regulation 16–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.

Abstract Background Amphioxus are non-vertebrate chordates characterized by a slow morphological and molecular evolution. They share the basic chordate body-plan and genome organization with vertebrates but lack their 2 R whole-genome duplications and their developmental complexity. For these reasons, amphioxus are frequently used as an outgroup to study vertebrate genome evolution and Evo-Devo. Aside from whole-genome duplications, genes continuously duplicate on a smaller scale. Small-scale duplicated genes can be found in both amphioxus and vertebrate genomes, while only the vertebrate genomes have duplicated genes product of their 2R whole-genome duplications. Here, we explore the history of small-scale gene duplications in the amphioxus lineage and compare it to small- and large-scale gene duplication history in vertebrates. Results We present a study of the European amphioxus ( Branchiostoma lanceolatum ) gene duplications thanks to a new, high-quality genome reference. We find that, despite its overall slow molecular evolution, the amphioxus lineage has had a history of small-scale duplications similar to the one observed in vertebrates. We find parallel gene duplication profiles between amphioxus and vertebrates, and conserved functional constraints in gene duplication. Moreover, amphioxus gene duplicates show levels of expression and patterns of functional specialization similar to the ones observed in vertebrate duplicated genes. We also find strong conservation of gene synteny between two distant amphioxus species, B. lanceolatum and B. floridae , with two major chromosomal rearrangements. Conclusions In contrast to their slower molecular and morphological evolution, amphioxus’ small-scale gene duplication history resembles that of the vertebrate lineage both in quantitative and in functional terms.

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.

Pluripotent cells are a transient population present in the early mammalian embryo dependent on transcription factors, such as OCT4 and NANOG, which maintain pluripotency while simultaneously suppressing lineage specification. Interestingly, these factors are not exclusive to uncommitted cells, but are also expressed during early phases of differentiation. However, their role in the transition from pluripotency to lineage specification is largely unknown. Using genetic models for controlled Oct4 or Nanog expression during postimplantation stages, we found that pluripotency factors play a dual role in regulating key lineage specifiers, initially repressing their expression and later being required for their proper activation. We show that the HoxB cluster is coordinately regulated in this way by OCT4 binding sites located at the 3' end of the cluster. Our results show that core pluripotency factors are not limited to maintaining the pre-committed epiblast, but are also necessary for the proper deployment of subsequent developmental programs.

Abstract Signaling pathways control a large number of gene regulatory networks (GRNs) during animal development, acting as major tools for body plan formation 1 . Remarkably, in contrast to the large number of transcription factors present in animal genomes, only a few of these pathways operate during development 2 . Moreover, most of them are largely conserved along metazoan evolution 3 . 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 Regulators of Complement Activation (RCA) gene cluster comprises several tandemly arranged genes which share functions in the innate immune system. RCA members, such as complement receptor 2 (CR2), are well-established susceptibility genes in complex autoimmune diseases. Altered expression of RCA genes has been demonstrated at both the functional and genetic level, but the mechanism underlying their regulation are not fully characterised. We aimed to investigate the structural organisation of the RCA gene cluster to identify key regulatory elements that influence the expression of CR2 and other genes in this immunomodulatory region. Using 4C, we captured extensive CTCF-mediated chromatin looping across the RCA gene cluster in B cells and showed these were organised into two topological associated domains (TADs). Interestingly, the inter-TAD boundary was located within the CR1 gene at a well-characterised segmental duplication. Additionally, we mapped numerous gene-gene and gene-enhancer interactions across the region, revealing extensive co-regulation. Importantly, we identified an intergenic enhancer and functionally demonstrated this element upregulates two RCA members (CR2 and CD55) in B cells. We have uncovered novel, long-range mechanisms whereby SLE susceptibility may be influenced by genetic variants, highlighting the important contribution of chromatin topology to gene regulation and complex genetic disease.

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.

ABSTRACT Vertebrate genomes organize into topologically associating domains (TADs), delimited by boundaries that insulate regulatory elements from non-target genes. However, how boundary function is established is not well understood. Here, we combine genome-wide analyses and transgenic mouse assays to dissect the regulatory logic of clustered-CTCF boundaries in vivo , interrogating their function at multiple levels: chromatin interactions, transcription and phenotypes. Individual CTCF binding sites (CBS) deletions revealed that the characteristics of specific sites can outweigh other factors like CBS number and orientation. Combined deletions demonstrated that CBS cooperate redundantly and provide boundary robustness. We show that divergent CBS signatures are not strictly required for effective insulation and that chromatin loops formed by non-convergently oriented sites could be mediated by a loop interference mechanism. Further, we observe that insulation strength constitutes a quantitative modulator of gene expression and phenotypes. Our results highlight the modular nature of boundaries and their control over developmental processes.