Search

Analysis of the transcriptomes and chromatin accessibility of multiple developmental stages of the indirect-developing hemichordate Ptychodera flava. After the individual analysis of each sample, integration of ATACseq and RNAseq data was performed in order to reconstruct genetic and regulatory networks. These data were also compared to other species, what could shed light on deuterostome evolution.

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.

[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 2R 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.

DNA methylation [5-methylcytosine (5mC)] is a repressive gene-regulatory mark required for vertebrate embryogenesis. Genomic 5mC is tightly regulated through the action of DNA methyltransferases, which deposit 5mC, and ten-eleven translocation (TET) enzymes, which participate in its active removal through the formation of 5-hydroxymethylcytosine (5hmC). TET enzymes are essential for mammalian gastrulation and activation of vertebrate developmental enhancers; however, to date, a clear picture of 5hmC function, abundance, and genomic distribution in nonvertebrate lineages is lacking. By using base-resolution 5mC and 5hmC quantification during sea urchin and lancelet embryogenesis, we shed light on the roles of nonvertebrate 5hmC and TET enzymes. We find that these invertebrate deuterostomes use TET enzymes for targeted demethylation of regulatory regions associated with developmental genes and show that the complement of identified 5hmC-regulated genes is conserved to vertebrates. This work demonstrates that active 5mC removal from regulatory regions is a common feature of deuterostome embryogenesis suggestive of an unexpected deep conservation of a major gene-regulatory module.

Neurons are a highly specialized cell type only found in metazoans. They can be scattered throughout the body or grouped together, forming ganglia or nerve cords. During embryogenesis, centralized nervous systems develop from the ectoderm, which also forms the epidermis. How pluripotent ectodermal cells are directed toward neural or epidermal fates, and to which extent this process is shared among different animal lineages, are still open questions. Here, by using micromere explants, we were able to define in silico the putative gene regulatory networks (GRNs) underlying the first steps of the epidermis and the central nervous system formation in the cephalochordate amphioxus. We propose that although the signal triggering neural induction in amphioxus (i.e., Nodal) is different from vertebrates, the main transcription factors implicated in this process are conserved. Moreover, our data reveal that transcription factors of the neural program seem to not only activate neural genes but also to potentially have direct inputs into the epidermal GRN, suggesting that the Nodal signal might also contribute to neural fate commitment by repressing the epidermal program. Our functional data on whole embryos support this result and highlight the complex interactions among the transcription factors activated by the signaling pathways that drive ectodermal cell fate choice in chordates.

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 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 NCMD-associated non-coding SNVs that we identified herein. 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, suggesting that this region is 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 provides insight into the cis-regulatory mechanisms of NCMD and supports that this condition is a retinal enhanceropathy.

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.

Pluripotent cells are a transient population of the mammalian embryo dependent on transcription factors, such as OCT4 and NANOG, which maintain pluripotency while suppressing lineage specification. However, these factors are also expressed during early phases of differentiation, and their role in the transition from pluripotency to lineage specification is largely unknown. 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 Oct4 is necessary for activation of HoxB genes during differentiation of embryonic stem cells and in the embryo. In addition, we show that the HoxB cluster is coordinately regulated 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 precommitted epiblast but are also necessary for the proper deployment of subsequent developmental programs.

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 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.

The Regulators of Complement Activation (RCA) gene cluster comprises several tandemly arranged genes with shared functions within the 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 mechanisms 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 topologically associated domains (TADs). Interestingly, an 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 autoimmune disease susceptibility may be influenced by genetic variants, thus highlighting the important contribution of chromatin topology to gene regulation and complex genetic disease.

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

Cis-regulatory elements (CREs) and transcription factors (TFs) associated with them determine temporal and spatial domains of gene expression. Therefore, identification of these CREs and TFs is crucial to elucidating transcriptional programs across taxa. With chromatin accessibility facilitating transcription factor access to DNA, the identification of regions of open chromatin sheds light both on the function of the regulatory elements and their evolution, thus allowing the recognition of potential CREs. Buenrostro and colleagues have developed a novel method for exploring chromatin accessibility: assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq), which can be used for the purpose of identifying putative CREs. This method was shown to have considerable advantages when compared to traditional methods such as sequence conservation analyses or functional assays. Here we present the adaptation of the ATAC-seq method to echinoderm species and discuss how it can be used for CRE discovery.

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.

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. 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. Although CTCF is essential for cycling cells and developing embryos, its in vitro removal has limited effects over gene expression, challenging the concept that CTCF-mediated chromatin interactions are a fundamental requirement for gene regulation. 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 is essential to regulate gene expression during embryonic development and for the establishment of developmental gene regulatory landscapes.

The causes and consequences of genome reduction in animals are unclear because our understanding of this process mostly relies on lineages with often exceptionally high rates of evolution. Here, we decode the compact 73.8-megabase genome of Dimorphilus gyrociliatus, a meiobenthic segmented worm. The D. gyrociliatus genome retains traits classically associated with larger and slower-evolving genomes, such as an ordered, intact Hox cluster, a generally conserved developmental toolkit and traces of ancestral bilaterian linkage. Unlike some other animals with small genomes, the analysis of the D. gyrociliatus epigenome revealed canonical features of genome regulation, excluding the presence of operons and trans-splicing. Instead, the gene-dense D. gyrociliatus genome presents a divergent Myc pathway, a key physiological regulator of growth, proliferation and genome stability in animals. Altogether, our results uncover a conservative route to genome compaction in annelids, reminiscent of that observed in the vertebrate Takifugu rubripes.