The little skate genome and the evolutionary emergence of wing-like fins.

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.

Evolution of 3D chromatin organization at different scales.

Most animal genomes fold in 3D chromatin domains called topologically associated domains (TADs) that facilitate interactions between cis-regulatory elements (CREs) and promoters. Owing to their critical role in the control of developmental gene expression, we explore how TADs have shaped animal evolution. In the light of recent studies that profile TADs in disparate animal lineages, we discuss their phylogenetic distribution and the mechanisms that underlie their formation. We present evidence indicating that TADs are plastic entities composed of genomic strata of different ages: ancient cores are combined with newer regions and brought into extant TADs through genomic rearrangements. We highlight that newly incorporated TAD strata enable the establishment of new CRE-promoter interactions and in turn new expression patterns that can drive phenotypical innovation. We further highlight how subtle changes in chromatin folding may fine-tune the expression levels of developmental genes and hold a potential for evolutionary significance.

Parallel evolution of amphioxus and vertebrate small-scale gene duplications.

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.

In vivo dissection of a clustered-CTCF domain boundary reveals developmental principles of regulatory insulation.

Vertebrate genomes organize into topologically associating domains, delimited by boundaries that insulate regulatory elements from nontarget 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-CCCTC-binding factor (CTCF) boundaries in vivo, interrogating their function at multiple levels: chromatin interactions, transcription and phenotypes. Individual CTCF binding site (CBS) deletions revealed that the characteristics of specific sites can outweigh other factors such as CBS number and orientation. Combined deletions demonstrated that CBSs 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 nonconvergently 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.

Pluripotency factors regulate the onset of Hox cluster activation in the early embryo.

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.

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

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

Regulatory Architecture of the RCA Gene Cluster Captures an Intragenic TAD Boundary, CTCF-Mediated Chromatin Looping and a Long-Range Intergenic Enhancer.

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.

MacroH2As regulate enhancer-promoter contacts affecting enhancer activity and sensitivity to inflammatory cytokines.

MacroH2A histone variants have a function in gene regulation that is poorly understood at the molecular level. We report that macroH2A1.2 and macroH2A2 modulate the transcriptional ground state of cancer cells and how they respond to inflammatory cytokines. Removal of macroH2A1.2 and macroH2A2 in hepatoblastoma cells affects the contact frequency of promoters and distal enhancers coinciding with changes in enhancer activity or preceding them in response to the cytokine tumor necrosis factor alpha. Although macroH2As regulate genes in both directions, they globally facilitate the nuclear factor κB (NF-κB)-mediated response. In contrast, macroH2As suppress the response to the pro-inflammatory cytokine interferon gamma. MacroH2A2 has a stronger contribution to gene repression than macroH2A1.2. Taken together, our results suggest that macroH2As have a role in regulating the response of cancer cells to inflammatory signals on the level of chromatin structure. This is likely relevant for the interaction of cancer cells with immune cells of their microenvironment.

Multidimensional chromatin profiling of zebrafish pancreas to uncover and investigate disease-relevant enhancers.

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.

Gain of gene regulatory network interconnectivity at the origin of vertebrates.

SignificanceIn this manuscript, we address an essential question in developmental and evolutionary biology: How have changes in gene regulatory networks contributed to the invertebrate-to-vertebrate transition? To address this issue, we perturbed four signaling pathways critical for body plan formation in the cephalochordate amphioxus and in zebrafish and compared the effects of such perturbations on gene expression and gene regulation in both species. Our data reveal that many developmental genes have gained response to these signaling pathways in the vertebrate lineage. Moreover, we show that the interconnectivity between these pathways is much higher in zebrafish than in amphioxus. We conclude that this increased signaling pathway complexity likely contributed to vertebrate morphological novelties during evolution.

CTCF knockout in zebrafish induces alterations in regulatory landscapes and developmental gene expression.

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.

Ancient Genomic Regulatory Blocks Are a Source for Regulatory Gene Deserts in Vertebrates after Whole-Genome Duplications.

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.

Pioneer and repressive functions of p63 during zebrafish embryonic ectoderm specification.

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.

Reprogramming Nuclear Architecture: Just a TAD.

In this issue of Cell Stem Cell, Laugsch et al. (2019) use direct reprogramming, epigenetics, and chromatin architecture studies to demonstrate that developmental defects observed in a BOFS patient are caused by reduced expression of TFAP2A in neural crest cells due to the spatial separation of the promoter from its neural crest enhancers.

Amphioxus functional genomics and the origins of vertebrate gene regulation.

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.

Computer simulation study of early bacterial biofilm development.

Most bacteria form organized sessile communities, known as biofilms. Their ubiquity and relevance have stimulated the development of efficient mathematical models able to predict biofilm evolution and characteristics at different conditions. Here we present a study of the early stages of bacterial biofilm formation modeled by means of individual cell-based computer simulation. Simulation showed that clusters with different degrees of internal and orientational order were formed as a function of the aspect ratio of the individual particles and the relation between the diffusion and growth rates. Analysis of microscope images of early biofilm formation by the Gram-negative bacterium Pseudomonas putida at varying diffusion rates revealed a good qualitative agreement with the simulation results. Our model is a good predictor of microcolony morphology during early biofilm development, showing that the competition between diffusion and growth rates is a key aspect in the formation of stable biofilm microcolonies.

A common copy-number variant within SIRPB1 correlates with human Out-of-Africa migration after genetic drift correction.

Previous reports have proposed that personality may have played a role on human Out-Of-Africa migration, pinpointing some genetic variants that were positively selected in the migrating populations. In this work, we discuss the role of a common copy-number variant within the SIRPB1 gene, recently associated with impulsive behavior, in the human Out-Of-Africa migration. With the analysis of the variant distribution across forty-two different populations, we found that the SIRPB1 haplotype containing duplicated allele significantly correlated with human migratory distance, being one of the few examples of positively selected loci found across the human world colonization. Circular Chromosome Conformation Capture (4C-seq) experiments from the SIRPB1 promoter revealed important 3D modifications in the locus depending on the presence or absence of the duplication variant. In addition, a 3' enhancer showed neural activity in transgenic models, suggesting that the presence of the CNV may compromise the expression of SIRPB1 in the central nervous system, paving the way to construct a molecular explanation of the SIRPB1 variants role in human migration.

4Cin: A computational pipeline for 3D genome modeling and virtual Hi-C analyses from 4C data.

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.

Regulatory landscape fusion in rhabdomyosarcoma through interactions between the PAX3 promoter and FOXO1 regulatory elements.

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.

Topologically associated domains: a successful scaffold for the evolution of gene regulation in animals.

The evolution of gene regulation is considered one of the main drivers causing the astonishing morphological diversity in the animal kingdom. Gene regulation in animals heavily depends upon cis-regulatory elements, discrete pieces of DNA that interact with target promoters to regulate gene expression. In the last years, Chromosome Conformation Capture experiments (4C-seq, 5C, and HiC) in several organisms have shown that the genomes of many bilaterian animals are organized in the 3D chromatin space in compartments called topologically associated domains (TADs). The appearance of the architectural protein CTCF in the bilaterian ancestor likely facilitated the origin of this chromatin 3D organization. TADs play a critical role favoring the contact of cis-regulatory elements with their proper target promoters (that often lay within the same TAD) and preventing undesired regulatory interactions with promoters located in neighboring TADs. We propose that TAD may have had a major influence in the history of the evolution of gene regulation. They have contributed to the increment of regulatory complexity in bilaterians by allowing newly evolved cis-regulatory elements to find target promoters in a range of hundreds of kilobases. In addition, they have conditioned the mechanisms of evolution of gene regulation. These mechanisms include the appearance, removal, or relocation of TAD borders. Such architectural changes have been able to wire or unwire promoters with different sets of cis-regulatory elements in a single mutational event. We discuss the contribution of these architectural changes to the generation of critical genomic 3D structures required for new regulatory mechanisms associated to morphological novelties. WIREs Dev Biol 2017, 6:e265. doi: 10.1002/wdev.265 For further resources related to this article, please visit the WIREs website.